Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

Molecular population genetics of accessory gland protein genes and testis-expressed genes in Drosophila mojavensis and D. arizonae

Molecular population genetic investigation of Drosophila male reproductive genes has focused primarily on melanogaster subgroup accessory gland protein genes (Acp's). Consistent with observations from male reproductive genes of numerous taxa, Acp's evolve more rapidly than nonreproductive genes. However, within the Drosophila genus, large data sets from additional types of male reproductive genes and from different species groups are lacking. Here we report findings from a molecular population genetics analysis of male reproductive genes of the repleta group species, Drosophila arizonae and D. mojavensis. We find that Acp's have dramatically higher average pairwise Ka/Ks (0.93) than testis-enriched genes (0.19) and previously reported melanogaster subgroup Acp's (0.42). Overall, 10 of 19 Acp's have Ka/Ks > 1 either in nonpolarized analyses or in at least one lineage of polarized analyses. Of the nine Acp's for which outgroup data were available, average Ka/Ks was considerably higher in D. mojavensis (2.08) than in D. arizonae (0.87). Contrasts of polymorphism and divergence suggest that adaptive protein evolution at Acp's is more common in D. mojavensis than in D. arizonae.     

Molecular evolution of sex-biased genes in Drosophila

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.     

Inference of positive and negative selection on the 5' regulatory regions of Drosophila genes

Both positive selection and negative selection have been shown to drive the evolution of coding regions. It is of interest to know if the corresponding 5' regions of genes may be subjected to selection of comparable intensities. For such a comparison, we chose the Accessory gland protein (Acp) genes as our test case. About 700 bp and 600 bp for the 5' and coding regions, respectively, of eight previously unstudied genes were sequenced from 21 isogenic lines of D. melanogaster and one line from D. simulans. The ratio of divergence at the amino-acid replacement sites (A) over that at the synonymous sites (S) was twice the ratio for common polymorphism. Interestingly, the 5' region shows the same trend, with the 5'/S divergence ratio being 1.8 times higher than the 5'/S ratio for common polymorphism. There are several possible explanations for the 5'/S ratios, including demography, negative selection, and positive selection. Under normal conditions, positive selection is the most likely explanation. If that is true, about 45 to 50 percent of all fixed differences at both the replacement and 5' sites were adaptive, even though the substitution rate in the former is only half that of the latter (K(A)/K(S) approximately 0.3 vs. K(5')/K(S) approximately 0.6). As previous analyses have indicated, the inclusion of slightly deleterious polymorphism confounds the inference of positive selection. The analysis of published polymorphism data covering 97 verified 5' regions of Drosophila suggests more pronounced selective constraint on the 5' untranslated region and the core promoter (together corresponding to approximately 200 bp in this data set) when compared to the more distal portion of the 5' region of genes.     

Molecular population genetics of male accessory gland proteins in Drosophila

Drosophila seminal proteins have an unusually high rate of molecular sequence evolution, suggesting either a high rate of neutral substitution or rapid adaptive evolution. To further quantify patterns of polymorphism and divergence in genes encoding seminal proteins, also called accessory gland proteins (Acp's), we conducted a sequencing survey of 10 Acp genes in samples of Drosophila melanogaster and D. simulans (Acp29AB, Acp32CD, Acp33A, Acp36DE, Acp53Ea, Acp62F, Acp63F, Acp76A, Acp95EF, and Acp98AB). Mean heterozygosity at replacement sites in D. simulans was 0.0074 for Acp genes and 0.0013 for a set of 19 non-Acp genes, and mean melanogaster-simulans divergence at replacement sites was 0.0497 for Acp genes and 0.0107 at non-Acp genes. The elevated divergence of Acp genes is thus accompanied by elevated within-species polymorphism. In addition to the already-reported departures of Acp26A, Acp29AB, and Acp70A from neutrality, our data reject neutrality at Acp29AB and Acp36DE in the direction of excess replacements in interspecific comparisons.     

Molecular population genetics and evolution of Drosophila meiosis genes

While many functional elements of the meiotic process are well characterized in model organisms, the genetic basis of most of the natural phenotypic variation observed in meiotic pathways has not been determined. To begin to address this issue, we characterized patterns of polymorphism and divergence in the protein-coding regions of 33 genes across 31 lines of Drosophila melanogaster and 6 lines of Drosophila simulans. We sequenced genes known to be involved in chromosome segregation, recombination, DNA repair, and related heterochromatin binding. As expected, we found several of the genes to be highly conserved, consistent with purifying selection. However, a subset of genes showed patterns of polymorphism and divergence typical of other types of natural selection. Moreover, several intriguing differences between the two Drosophila lineages were evident: along the D. simulans lineage we consistently found evidence of adaptive protein evolution, whereas along the D. melanogaster lineage several loci exhibited patterns consistent with the maintenance of protein variation.     

Adaptive evolution of newly emerged micro-RNA genes in Drosophila

How often micro-RNA (miRNA) genes emerged and how fast theyevolved soon after their emergence are some of the central questionsin the evolution of miRNAs. Because most known miRNA genes areancient and highly conserved, these questions can be best answeredby identifying newly emerged miRNA genes. Among the 78 miRNAgenes in Drosophila reported before 2007, only 5 are confirmedto be newly emerged in the genus (although many more can befound in the newly reported data set; e.g., Ruby et al. 2007;Stark et al. 2007; Lu et al. 2008). These new miRNA genes haveundergone numerous changes, even in the normally invariant maturesequences. Four of them (the miR-310/311/312/313 cluster, denotedmiR-310s) were duplicated from other conserved miRNA genes.The fifth one (miR-303) appears to be a very young gene, originatingde novo from a non-miRNA sequence recently. We sequenced these5 miRNA genes and their neighboring regions from a worldwidecollection of Drosophila melanogaster lines. The levels of divergenceand polymorphism in these miRNA genes, vis-à-vis thoseof the neighboring DNA sequences, suggest that these 5 genesare evolving adaptively. Furthermore, the polymorphism patternof miR-310s in D. melanogaster is indicative of hitchhikingunder positive selection. Thus, a large number of adaptive changesover a long period of time may be essential for the evolutionof newly emerged miRNA genes.     

Molecular evolution of duplicate copies of genes encoding cytosolic glutamine synthetase in Pisum sativum

Here, we describe two nearly identical expressed genes for cytosolic glutamine synthetase (GS3A and GS3B) in Pisum sativum L. RFLP mapping data indicates that the GS3A and GS3B genes are separate loci located on different chromosomes. DNA sequencing of the GS3A and GS3B genes revealed that the coding regions are 99% identical with only simple nucleotide substitutions resulting in three amino acid differences. Surprisingly, the non-coding regions (5 non-coding leader, the 11 introns, and 3 non-coding tail) all showed a high degree of identity (96%). In these non-coding regions, 25% of the observed differences between the GS3A and GS3B genes were deletions or duplications. The single difference in the 3 non-coding regions of the GS3A and GS3B genes was a 25 bp duplication of an AU-rich element in the GS3B gene. As the GS3B mRNA accumulates to lower levels than the GS3A gene, we tested whether this sequence which resembles an mRNA instability determinant functioned as such in the context of the GS mRNA. Using the GS3B 3 tail as part of a chimeric gene in transgenic plants, we showed that this AU-rich sequence has little effect on transgene mRNA levels. To determine whether the GS3A/GS3B genes represent a recent duplication, we examined GS3-like genes in genomic DNA of ancient relatives of P. sativum. We observed that several members of the Viceae each contain two genomic DNA fragments homologous to the GS3B gene, suggesting that this is an ancient duplication event. Gene conversion has been invoked as a possible mechanism for maintaining the high level of nucleotide similarity found between the GS3A and GS3B genes. Possible evolutionary reasons for the maintenance of these twin GS genes in pea, and the general duplication of genes for cytosolic GS in all plant species are discussed.     

Evidence for male-driven evolution in Drosophila

In several vertebrate taxa studied to date, mutation rates arehigher in males than females (male-driven evolution). The male-to-femalemutation rate () can be estimated by contrasting DNA divergencedata at X-linked, Y-linked, and autosomal loci. Previous studiesin Drosophila, comparing X-linked and autosomal divergence,have found no evidence for male-driven evolution in this genus.Here, I compare levels of nucleotide divergence between homologousX- and Y-linked loci in Drosophila miranda. Using divergenceat both synonymous sites and at short introns, I estimate tobe approximately 2. This study thus provides the first evidencefor male-biased mutation rates outside vertebrates, supportingthe view that DNA sequence evolution is male driven in a widevariety of taxa.     

Ultrastructure and nature of secretory proteins in the male accessory gland of Drosophila funebris

The ultrastructural differentiation of the secretory cells and the nature of secretory proteins in the male accessory gland of Drosophila funebris have been studied by electron-microscopic and immunological methods. (1) In the pupae at $\text{1}\text{1}\text{2}$ days before eclosion, secretory products can be detected in the lumen, even though most glandular cells are at the initial phase of differentiation. At the time of eclosion both main and secondary cells are fully differentiated, but the whole set of five immunologically active proteins are detectable only on the second to third day of adult life. (2) The secondary cells contain giant protein granules, the so-called filamentous bodies, which become partially fused and the filaments assume a twisted form. Randomly dispersed filaments and closely packed filament bundles are also visible in the gland lumen. Antigenic labelling of ultrathin sections and immunoreplica electrophoresis yielded no evidence for the microtubular nature of these filaments. The secretion stored in the lumen contains in addition a large quantity of flocculent proteins which have their origin in the main cells. (3) During the period of high secretory activity in the 7-day-old male flies no vacuolization and disintegration of either the main or secondary cells have been observed. We conclude that both types of cells have the merocrine secretory mechanism. (4) Ultrastructural alterations in the glandular cells confirmed our previous observation that copulation stimulates RNA and protein synthesis.     

Molecular population genetics and the search for adaptive evolution in plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.     

Faster evolution of Z-linked duplicate genes in chicken

It has been shown that duplicate genes on the X chromosome evolve much faster than duplicate genes on autosomes in Drosophila melanogaster.However,whether this phenomenon is general and can be applied to other species is not known.Here we examined this issue in chicken that have heterogametic females(females have ZW sex chromosome).We compared sequence divergence of duplicate genes on the Z chromosome with those on autosomes.We found that duplications on the Z chromosome indeed evolved faster than those on autosomes and show distinct patterns of molecular evolution from autosomal duplications.Examination of the expression of duplicate genes revealed an enrichment of duplications on the Z chromosome having male-biased expression and an enrichment of duplications on the autosomes having female-biased expression.These results suggest an evolutionary trend of the recruitment of duplicate genes towards reproduction-specific function.The faster evolution of duplications on Z than on the autosomes is most likely contributed by the selective forces driving the fixation of adaptive mutations on Z.Therefore,the common phenomena observed in both flies and chicken suggest that duplicate genes on sex chromosomes have distinct dynamics and are more influenced by natural selection than antosomal duplications,regardless of the kind of sex determination systems.     

Multiple origins and rapid evolution of duplicated mitochondrial genes in parthenogenetic geckos (Heteronotia binoei; Squamata, Gekkonidae)

Accumulating evidence for alternative gene orders demonstrates that vertebrate mitochondrial genomes are more evolutionarily dynamic than previously thought. Several lineages of parthenogenetic lizards contain large, tandem duplications that include rRNA, tRNA, and protein-coding genes, as well as the control region. Such duplications are hypothesized as intermediate stages in gene rearrangement, but the early stages of their evolution have not been previously studied. To better understand the evolutionary dynamics of duplicated segments of mitochondrial DNA, we sequenced 10 mitochondrial genomes from recently formed ( approximately 300,000 years ago) hybrid parthenogenetic geckos of the Heteronotia binoei complex and 1 from a sexual form. These genomes included some with an arrangement typical of vertebrates and others with tandem duplications varying in size from 5.7 to 9.4 kb, each with different gene contents and duplication endpoints. These results, together with phylogenetic analyses, indicate independent and frequent origins of the duplications. Small, direct repeats at the duplication endpoints imply slipped-strand error as a mechanism generating the duplications as opposed to a false initiation/termination of DNA replication mechanism that has been invoked to explain duplications in other lizard mitochondrial systems. Despite their recent origin, there is evidence for nonfunctionalization of genes due primarily to deletions, and the observed pattern of gene disruption supports the duplication-deletion model for rearrangement of mtDNA gene order. Conversely, the accumulation of mutations between these recent duplicates provides no evidence for gene conversion, as has been reported in some other systems. These results demonstrate that, despite their long-term stasis in gene content and arrangement in some lineages, vertebrate mitochondrial genomes can be evolutionary dynamic even at short timescales.     

Recent gene conversions between duplicated glutamate decarboxylase genes (gadA and gadB) in pathogenic Escherichia coli

Escherichia coli have evolved adaptive systems to resist strongly acidic habitats in part through the production of 2 biochemically identical isoforms of glutamate decarboxylase (GAD), encoded by the gadA and gadB genes. These genes occur in E. coli and other members of the genospecies (e.g., Shigella spp.) and originated as part of a genomic fitness island acquired early in Escherichia evolution. The present duplicated gad loci are widely spaced on the E. coli chromosome, and the 2 genes are 97% similar in sequence. Comparison of the nucleotide sequences of the gadA and gadB in 16 strains of pathogenic E. coli revealed 3.8% and 5.0% polymorphism in the 2 genes, respectively. Alignment of the homologous genes identified a total of 120 variable sites, including 21 fixed nucleotide differences between the loci within the first 82 codons of the genes. Twenty-three phylogenetically informative sites were polymorphic for the same nucleotides in both genes suggesting recent gene conversions or intergenic recombination. Phylogenetic analysis based on the synonymous substitutions per synonymous site indicated 2 cases in which specific gadA and gadB alleles were more closely related to one another than to other alleles at the corresponding locus. The results indicate that at least 3 gene conversion events have occurred after the gad gene duplication in the evolution of E. coli. Despite multiple gene conversion events, the upstream regulatory regions and the 5' end of each gene remains distinct, suggesting that maintaining functionally different gad genes is important in this acid-resistance mechanism in pathogenic E. coli.