Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex.

The Acp26Aa and Acp26Ab genes that code for male accessory gland proteins are tandemly arranged in the species of the Drosophila melanogaster complex. An approximately 1.6-kb region encompassing both genes has been sequenced in 10, 24, and 18 lines from Spain, Ivory Coast, and Malawi, respectively; the previously studied 10 lines from North Carolina have also been included in the analyses. A total of 110 nucleotide and 4 length polymorphisms were detected. Silent variation for the whole Acp26A region was slightly higher in African than in non-African populations, while for both genes nonsynonymous variation was similar in all populations studied. Based on Fst estimates no major genetic differentiation was detected between East and West Africa, while in general non-African populations were strongly differentiated from both African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites revealed that directional selection acting on amino acid replacement changes has driven the evolution of the Acp26Aa protein in the last 2.5 myr.     

Sex in Drosophila mauritiana: a very high level of amino acid polymorphism in a male reproductive protein gene, Acp26Aa

Many genes pertaining to male reproductive functions have been shown to evolve rapidly between species, and evidence increasingly suggest the influence of positive Darwinian selection. The accessory gland protein gene (Acp26Aa) of Drosophila is one such example. In order to understand the mechanism of selection, it is often helpful to examine the pattern of polymorphism. We report here that the level of amino acid polymorphism in the N-terminal quarter of Acp26Aa is high in Drosophila melanogaster and is unprecedented in its sibling species Drosophila mauritiana. We postulate that (1) this N-terminal segment may play a role in sperm competition, and (2) D. mauritiana may have been under much more intense sexual selection than other species. Both postulates have important ramifications and deserve to be tested rigorously.     

Positive selection and the molecular evolution of a gene of male reproduction, Acp26Aa of Drosophila

The gene for a male ejaculatory protein, Acp26Aa, in four sibling species of the Drosophila melanogaster subgroup has previously been shown to have a nonsynonymous rate (Ka) of nucleotide substitution that is indistinguishable from the synonymous rate (Ks). By examining this gene in two other species of this subgroup, we found that Ka is generally large and can sometimes be more than twice as large as Ks. This suggests that positive selection may be operating at this locus of male reproduction.      

Evidence of a high rate of selective sweeps in African Drosophila melanogaster

Assessing the rate of evolution depends on our ability to detect selection at several genes simultaneously. We summarize DNA sequence variation data in three new and six previously published data sets from the left arm of the second chromosome of Drosophila melanogaster in a population from West Africa, the presumed area of origin of this species. Four loci [Acp26Aa, Fbp2, Vha68-1, and Su(H)] were previously found to deviate from a neutral mutation-drift equilibrium as a consequence of one or several selective sweeps. Polymorphism data from five loci from intervening regions (dpp, Acp26Ab, Acp29AB, GH10711, and Sos) did not show the characteristic deviation from neutrality caused by local selective sweeps. This genomic region is polymorphic for the In(2L)t inversion. Four loci located near inversion breakpoints [dpp, sos, GH10711, and Su(H)] showed significant structuring between the two arrangements or significant deviation from neutrality in the inverted class, probably as a result of a recent shift in inversion frequency. Overall, these patterns of variation suggest that the four selective events were independent. Six loci were observed with no a priori knowledge of selection, and independent selective sweeps were detected in three of them. This suggests that a large part of the D. melanogaster genome has experienced the effect of positive selection in its ancestral African range.     

Positive selection drives the evolution of the Acp29AB accessory gland protein in Drosophila.

Nucleotide sequence variation at the Acp29AB gene region has been surveyed in Drosophila melanogaster from Spain (12 lines), Ivory Coast (14 lines), and Malawi (13 lines) and in one line of D. simulans. The approximately 1.7-kb region studied encompasses the Acp29AB gene that codes for a male accessory gland protein and its flanking regions. Seventy-seven nucleotide and 8 length polymorphisms were detected. Nonsynonymous polymorphism was an order of magnitude lower than synonymous polymorphism, but still high relative to other non-sex-related genes. In D. melanogaster variation at this region revealed no major genetic differentiation between East and West African populations, while differentiation was highly significant between the European and the two African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites showed an excess of fixed nonsynonymous changes, which indicates that the evolution of the Acp29AB protein has been driven by directional selection at least after the split of the D. melanogaster and D. simulans lineages. The pattern of variation in extant populations of D. melanogaster favors a scenario where the fixation of advantageous replacement substitutions occurred in the early stages of speciation and balancing selection is maintaining variation in this species.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Evolutionary History of the Sex-Peptide (Acp70a) Gene Region in Drosophila Melanogaster

In Drosophila the products of the seminal fluid stimulate oviposition and suppress remating in the female. Of all the accessory gland peptides (Acp's) involved in these two responses, the sex-peptide (coded by the Acp70A gene) is among the best characterized at the functional level. A 1.2-kb fragment encompassing the Acp70A gene of nine lines from a natural population of D. melanogaster and one allele of D. sechellia was sequenced to study the forces shaping nucleotide variation within and between species. The coding region of D. simulans and D. mauritiana was also sequenced. A Ser to Ala replacement polymorphism at the last position of the signal peptide was detected in D. melanogaster. The Ser and Ala alleles are at intermediate frequencies. The level of nucleotide variation is lower for the derived Ala allele, which is compatible with a recent origin and an increase in frequency due to positive selection. Variation at the 5' flanking region is structured in two major highly differentiated haplotypes, whose distribution does not conform to neutral expectations. Selective and/or historical factors could contribute to the observed overall patterning of nucleotide variation at the Acp70A region.     

The Acp26Aa seminal fluid protein is a modulator of early egg hatchability in Drosophila melanogaster

Drosophila melanogaster male accessory gland proteins (Acps) that are transferred in the ejaculate with sperm mediate post-mating competition for fertilizations between males. The actions of Acps include effects on oviposition and ovulation, receptivity and sperm storage. Two Acps that modulate egg production are Acp26Aa (ovulin) and Acp70A (the sex peptide). Acp26Aa acts specifically on the process of ovulation (the release of mature eggs from the ovaries), which is initiated 1.5 h after mating. In contrast, sperm storage can take as long as 6-9 h to complete. Initial ovulations after matings by virgin females will therefore occur before all sperm are fully stored and the extra eggs initially laid as a result of Acp26Aa transfer are expected to be inefficiently fertilized. Acp26Aa-mediated release of existing eggs should not cause a significant energetic cost or lead to a decrease in female lifespan assuming, as seems likely, that the energetic cost of egg laying comes from de novo egg synthesis (oogenesis) rather than from ovulation. We tested these predictions using Acp26Aa(1) mutant males that lack Acp26Aa but are normal for other Acps and Acp26Aa(2) males that transfer a truncated but fully functional Acp26Aa protein. Females mating with Acp26Aa(2) (truncation) males that received functional Acp26Aa produced significantly more eggs following their first matings than did mates of Acp26Aa(1) (null) males. However, as predicted above, these extra eggs, which were laid as a result of Acp26Aa transfer to virgin females, showed significantly lower egg hatchability. Control experiments indicated that this lower hatchability was due to lower rates of fertilization at early post-mating times. There was no drop in egg hatchability in subsequent non-virgin matings. In addition, as predicted above, females that did or did not receive Acp26Aa did not differ in survival, lifetime fecundity or lifetime progeny, indicating that Acp26Aa transfer does not represent a significant energetic cost for females and does not contribute to the survival cost of mating. Acp26Aa appears to remove a block to oogenesis by causing the clearing out of existing mature eggs and, thus, indirectly allowing oogenesis to be initiated immediately after mating. The results show that subtle processes coordinate the stimulation of egg production and sperm storage in mating pairs.     

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.     

The molecular evolution of terminal ear1, a regulatory gene in the genus Zea.

Nucleotide diversity in the terminal ear1 (te1) gene, a regulatory locus hypothesized to be involved in the morphological evolution of maize (Zea mays ssp. mays), was investigated for evidence of past selection. Nucleotide polymorphism in a 1.4-kb region of te1 was analyzed for a sample of 26 sequences isolated from 12 maize lines, five populations of the maize progenitor, Z. mays ssp. parviglumis, six other Zea populations, and two Tripsacum species. Although nucleotide diversity in te1 in maize is reduced relative to ssp. parviglumis, phylogenetic and statistical analyses of the pattern of polymorphism among these sequences provided no evidence of past selection, indicating that the region of the gene studied was probably not involved in maize evolution. The level of reduction in genetic diversity in te1 in maize relative to its progenitor is comparable to that found in previous reports for isozymes and other neutrally evolving maize genes and is consistent with a genome-wide reduction of genetic diversity resulting from a domestication bottleneck. An estimate of the age (1.2-1.4 million yr) of the maize gene pool based on te1 is roughly consistent with previous estimates based on other neutral genes, but may be biased by the apparently slow synonymous substitution rate at te1.     

Detecting the footprint of positive selection in a european population of Drosophila melanogaster: multilocus pattern of variation and distance to coding regions

The effects on nucleotide variation of adaptations to temperate habitats and of the possible bottleneck associated with the origin of European populations of Drosophila melanogaster should be detectable in DNA sequences given the short time elapsed relative to the species population size. We surveyed nucleotide variation in 109 fragments distributed across the X chromosome in a European population of D. melanogaster to detect the footprint of positive selection. Fragments were located primarily in large noncoding regions. Multilocus tests based on Tajima's D statistic revealed a significant departure from neutral expectations in a stationary panmictic population, with an important contribution from both positive and negative D values. A positive relationship between Tajima's D values and distance to coding region was detected, with a comparative excess of significantly negative D values in the subset of fragments closer to coding regions. Also, there was a significant heterogeneity in the polymorphism to divergence ratio, with 12 fragments contributing 42% to the test statistic. Moreover, these fragments were comparatively closer to coding regions. These findings would imply positive selection events, and thus selective sweeps, during the species expansion to Europe.     

Distribution of genetic variation in the growth hormone 1 gene in Atlantic salmon (Salmo salar) populations from Europe and North America

The level and hierarchical distribution of genetic variation in complete sequences of the Atlantic salmon (Salmo salar) growth hormone (GH1) gene were investigated in populations from Europe and North America with a view to inferring the major evolutionary forces affecting genetic variation at this locus. Seventeen polymorphic sites were identified in complete sequences from nine populations, with levels of noncoding (intron and untranslated region sequences) nucleotide diversity being similar to those observed in other species. No variation, however, was observed in exonic sequences, indicating that nucleotide diversity in the Atlantic salmon GH1 gene is three and 25 times less than that estimated for human and Drosophila coding sequences, respectively. This suggests that purifying selection is the predominant contemporary force controlling the molecular evolution of GH1 coding sequences. Comparison of haplotype relationships within and between populations indicated that differentiation between populations from Europe and North America was greater than within-continent comparisons. However, several haplotypes observed in the northernmost European populations were more similar to those observed in North American than to any other haplotypes observed in Europe. This is most likely to be a result of historical, rather than contemporary, gene flow. Neutrality test statistics, such as Tajima's D, were significantly positive in the European populations in which North American-like haplotypes were observed. Although a positive Tajima's D is commonly interpreted as the signal of balancing selection, a more likely explanation in this case is that either historical migration or ascertainment bias, rather than within population local adaptation, has given rise to an excess of intermediate frequency alleles.     

Polymorphism and Divergence at a Drosophila Pseudogene Locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcpψ) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcpψ fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcpψ revealed lower levels of polymorphism ( & = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcpψ within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcpψ lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

S-element insertions are associated with the evolution of the Hsp70 genes in Drosophila melanogaster

The "selfish DNA" theory postulates that transposable elements (TEs) are intragenomic parasites, and that natural selection against deleterious effects associated with their presence is the main force preventing their genomic spread in natural populations. In agreement with this model, TEs in Drosophila melanogaster populations are usually found at low frequencies in most genomic locations. Only a few cases of fixation of TE insertions have been reported, usually in regions of low recombination, where selection is expected to be less effective. Here, we report a population genetics study on the apparent fixation of an S-element in a highly recombining region in two natural populations of D. melanogaster. Three similar fragments of an S-element are inserted into the 5' regions of three members of a heat shock gene family, Hsp70 (Hsp70Aa and Hsp70Ab in polytene chromosome band 87A, and Hsp70Bb in 87C). A PCR-based analysis suggests that the insertions are fixed or at high frequencies in the entire species. A population survey of the levels of nucleotide sequence variation at the insertion site in 87C in two natural populations of D. melanogaster provided evidence for reduced levels of variation in the region, normal levels of recombination, and selection, reflected in a significant departure from neutrality of the variant frequency spectrum. This was particularly strong for the S-element inverted repeats (IRs) and suggests that these are of functional significance for the host.     

Variation in Sperm Displacement and Its Association with Accessory Gland Protein Loci in Drosophila Melanogaster

Genes that influence mating and/or fertilization success may be targets for strong natural selection. If females remate frequently relative to the duration of sperm storage and rate of sperm use, sperm displacement may be an important component of male reproductive success. Although it has long been known that mutant laboratory stocks of Drosophila differ in sperm displacement, the magnitude of the naturally occurring genetic variation in this character has not been systematically quantified. Here we report the results of a screen for variation in sperm displacement among 152 lines of Drosophilia melanogaster that were made homozygous for second and/or third chromosomes recovered from natural populations. Sperm displacement was assayed by scoring the progeny of cn;bw females that had been mated sequentially to cn;bw and tested males in either order. Highly significant differences were seen in both the ability to displace sperm that is resident in the female's reproductive tract and in the ability to resist displacement by subsequent sperm. Most lines exhibited nearly complete displacement, having nearly all progeny sired by the second male, but several lines had as few as half the progeny fathered by the second male. Lines that were identified in the screen for naturally occurring variation in sperm displacement were also characterized for single-strand conformation polymorphisms (SSCP) at seven accessory gland protein (Acp) genes, Glucose dehydrogenase (Gld), and Esterase-6 (Est-6). Acp genes encode proteins that are in some cases known to be transmitted to the female in the seminal fluid and are likely candidates for genes that might mediate the phenomenon of sperm displacement. Significant associations were found between particular Acp alleles at four different loci (Acp26Aa/Ab, Acp29B, Acp36DE and Acp53E) and the ability of males to resist displacement by subsequent sperm. There was no correlation between the ability to displace resident sperm and the ability to resist being displaced by subsequent sperm. This lack of correlation, and the association of Acp alleles with resisting subsequent sperm only, suggests that different mechanisms mediate the two components of sperm displacement.     

Nucleotide variation of the Est-6 gene region in natural populations of Drosophila melanogaster

We have investigated nucleotide polymorphism in the Est-6 gene region in four samples of Drosophila melanogaster derived from natural populations of East Africa (Zimbabwe), Europe (Spain), North America (California), and South America (Venezuela). There are two divergent sequence types in the North and South American samples, which are not perfectly (North America) or not at all (South America) associated with the Est-6 allozyme variation. Less pronounced or no sequence dimorphism occurs in the European and African samples, respectively. The level of nucleotide diversity is highest in the African sample, lower (and similar to each other) in the samples from Europe and North America, and lowest in the sample from South America. The extent of linkage disequilibrium is low in Africa (1.23% significant associations), but much higher in non-African populations (22.59, 21.45, and 37.68% in Europe, North America, and South America, respectively). Tests of neutrality with recombination are significant in non-African samples but not significant in the African sample. We propose that demographic history (bottleneck and admixture of genetically different populations) is the major factor shaping the nucleotide patterns in the Est-6 gene region. However, positive selection modifies the pattern: balanced selection creates elevated levels of nucleotide variation around functionally important (target) polymorphic sites (RsaI-/RsaI+ in the promoter region and F/S in the coding region) in both African and non-African samples; and directional selection, acting during the geographic expansion phase of D. melanogaster, creates an excess of very similar sequences (RsaI- and S allelic lineages, in the promoter and coding regions, respectively) in the non-African samples.     

Adaptive microclimatic evolution of the dehydrin 6 gene in wild barley at ��Evolution Canyon��, Israel

Dehydrins are one of the major stress-induced gene families, and the expression of dehydrin 6 (Dhn6) is strictly related to drought in barley. In order to investigate how the evolution of the Dhn6 gene is associated with adaptation to environmental changes, we examined 48 genotypes of wild barley, Hordeum spontaneum, from "Evolution Canyon" at Mount Carmel, Israel. The Dhn6 sequences of the 48 genotypes were identified, and a recent insertion of 342?bp at 5'UTR was found in the sequences of 11 genotypes. Both nucleotide and haplotype diversity of single nucleotide polymorphism in Dhn6 coding regions were higher on the AS ("African" slope or dry slope) than on the ES ("European" slope or humid slope), and the applied Tajima D and Fu-Li test rejected neutrality of SNP diversity. Expression analysis indicated that the 342?bp insertion at 5'UTR was associated with the earlier up-regulation of Dhn6 after dehydration. The genetic divergence of amino acids sequences indicated significant positive selection of Dhn6 among the wild barley populations. The diversity of Dhn6 in microclimatic divergence slopes suggested that Dhn6 has been subjected to natural selection and adaptively associated with drought resistance of wild barley at "Evolution Canyon".     

Positive Selection on Nucleotide Substitutions and Indels in Accessory Gland Proteins of the Drosophila pseudoobscura Subgroup

Genes encoding reproductive proteins often diverge rapidly due to positive selection on nucleotide substitutions. While this general pattern is well established, the extent to which specific reproductive genes experience similar selection in different clades has been little explored, nor have possible targets of positive selection other than nucleotide substitutions, such as indels, received much attention. Here, we inspect for the signature of positive selection in the genes encoding five accessory gland proteins (Acps) (Acp26Aa, Acp32CD, Acp53Ea, Acp62F, and Acp70A) originally described from Drosophila melanogaster but with recognizable orthologues in the D. pseudoobscura subgroup. We compare patterns of selection within the D. psuedoobscura subgroup to those in the D. melanogaster subgroup. Similar patterns of positive selection were found in Acp26Aa and Acp62F in the two subgroups, while Acp53Ea and Acp70A experienced purifying selection in both subgroups. These proteins have thus remained targets for similar types of selection over long (>21-MY) periods of time. We also found several indel substitutions and polymorphisms in Acp26Aa and Acp32CD. These indels occur in the same regions as positively selected nucleotide substitutions for Acp26Aa in the D. pseudoobscura subgroup but not in the D. melanogaster subgroup. Rates of indel substitution within Acp26Aa in the D. pseudoobscura subgroup were up to several times those in noncoding regions of the Drosophila genome. This suggests that indel substitutions may be under positive selection and may play a key role in the divergence of some Acps. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Willis Swanson]