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.     

Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila: II. Divergence versus polymorphism

The evolution of the gene for a male ejaculatory protein, Acp26Aa, has been shown to be driven by positive selection when nonsibling species in the Drosophila melanogaster subgroup are compared. To know if selection has been operating in the recent past and to understand the details of its dynamics, we obtained DNA sequences of Acp26Aa and the nearby Acp26Ab gene from 39 D. melanogaster chromosomes. Together with the 10 published sequences, we analyzed 49 sequences from five populations in four continents. The southern African population is somewhat differentiated from all other populations, but its nucleotide diversity is lower at these two loci. We find the following results for Acp26Aa: (1) The R: S (replacement : silent changes) ratio is significantly higher in the between-species comparisons than in the within-species data by the McDonald and Kreitman test. Positive selection is probably responsible for the excess of amino acid replacements between species. (2) However, within-species nucleotide diversity is high. Neither the Tajima test nor the Fu and Li test indicates a reduction in nucleotide diversity due to positive selection in the recent past. (3) The newly derived nucleotides in D. melanogaster are at high frequency significantly more often than predicted by the neutral equilibrium. Since the nearby Acp26Ab gene does not show these patterns, these observations cannot be attributed to the characteristics of this chromosomal region. We suggest that positive selection is active, but may be weak, for each amino acid change in the Acp26Aa gene.      

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.     

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.      

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in the postmating reactions of the female, is a single-copy gene. However, in Drosophila subobscura, the gene is duplicated and both copies are transcribed. To study the molecular evolution of the duplication, a 2.1- kb fragment encompassing both copies of the duplication was sequenced for 10 lines of D. subobscura and one line of Drosophila madeirensis. Estimates of the divergence between the two copies of the duplicated region and between the two species studied, D. subobscura and D. madeirensis, revealed that both copies of the Acp70a gene had evolved independently since their duplication. The ratio of nonsynonymous to silent divergence between copies was generally higher than one. The McDonald and Kreitman test revealed an excess of nonsynonymous changes fixed since the duplication and before the split of the D. subobscura and D. madeirensis lineages. These results point to natural selection driving protein evolution after the duplication. Specifically, adaptive evolution appears to have caused the initial differentiation between copies of the N-terminal parts of the proteins, while purifying selection could be responsible for the high conservation of the C- terminal parts.      

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.     

The Drosophila seminal fluid protein Acp26Aa stimulates release of oocytes by the ovary

Mating stimulates the rate of egg-laying by female insects. In Drosophila melanogaster this stimulation is initially caused by seminal fluid molecules transferred from the male (Acps or accessory gland proteins; reviewed in [1] [2] [3]). Egg-laying is a multi-step process. It begins with oocyte release by the ovaries, followed by egg movement down the oviducts and the deposition of eggs onto the substratum. Although two Acps are known to stimulate egg-laying [4] [5], they were detected by assays that do not discriminate between the steps of this process or allow examination of its earliest changes [4] [5] [6] [7]. To determine how egg-laying is regulated, we developed a generally applicable assay to separate the process into quantifiable steps, allowing us to assess the ovulation pattern and rate of egg movement. As the steps are interdependent yet potentially subject to independent controls, we determined the contribution of each step and effector independent of the others. We used a statistical method [8] [9] that separately considers and quantifies each 'path' to a common end. We found that the prohormone-like molecule Acp26Aa [5] [10] stimulates the first step in egg-laying - release of oocytes by the ovary. During mating, Acp26Aa begins to accumulate at the base of the ovaries, a position consistent with action on the ovarian musculature to mediate oocyte release. Understanding how individual Acps regulate egg-laying in fruitflies will help provide a full molecular picture of insects' prodigious fertility, of reproductive hormones, and of the roles of these rapidly evolving proteins [11] [12].     

Nucleotide and copy-number polymorphism at the odorant receptor genes Or22a and Or22b in Drosophila melanogaster

In Drosophila, odorant receptors are encoded by an old and moderatelysized multigene family. Or22a and Or22b are two tandemly arrangedgenes of this family that have proved to be the result of arather young duplication. Nucleotide variation in the regionspanning both duplicates was surveyed in four natural populations(two African and two non-African) of Drosophila melanogasterand also analyzed in species of the melanogaster subgroup. Theintraspecific survey revealed a particular copy-number polymorphismin some of the studied populations, with the two genes (Or22aand Or22b) present in the long variant and a single chimericgene (Or22ab) present in the short variant. Estimated nucleotidediversity was higher in the short than in the long variant,despite the ancestral character of the latter variant in D.melanogaster. The general skew toward low-frequency variantsdetected in the non-African long variant and its reduced levelof silent polymorphism relative to divergence is consistentwith the recent fixation of an advantageous mutation at, ornearby, the Or22 long variant region. The nonnegligible frequencyof the short variant and the presence of a highly divergenthaplotype in the East African sample would point to direct orindirect selection for its maintenance in the species. Therewas evidence for a generally more rapid evolution of the Or22bcopy at both synonymous and nonsynonymous sites. However, anexcess of nonsynonymous substitutions was only detected in theearly history of this copy.     

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.     

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.     

Polymorphism and Locus-Specific Effects on Polymorphism at Microsatellite Loci in Natural Drosophila Melanogaster Populations

We have studied the natural variation at microsatellite loci in two African and five non-African populations of Drosophila melanogaster. Ten dinucleotide simple sequence loci were cloned from chromosomally mapped P1 clones and typed for single individuals from isofemale lines of the respective populations. We find that the African populations harbor the largest degree of diversity, while the non-African populations show a lower diversity. This supports previous results that D. melanogaster originated in Africa and spread across the rest of the world in historic times. Using genetic distance measures, we find also a distinct population subdivision between the non-African populations. Most interestingly, we find for some loci in some populations a strongly reduced variability, which cannot be explained by bottleneck effects. Employing a conservative test based on the variance in repeat number, we find that at least one locus in one population deviates significantly from the expectations of mutation-drift equilibrium. We suggest that this may be due to a recent selective sweep in this chromosomal region that may have been caused by a linked locus that was involved in local adaptation of the population.     

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.     

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.     

DNA sequence variation and latitudinal associations in hsp23, hsp26 and hsp27 from natural populations of Drosophila melanogaster

Heat shock genes are considered to be likely candidate genes for environmental stress resistance. Nucleotide variation in the coding sequence of the small heat shock genes (hsps) hsp26 and hsp27 from Drosophila melanogaster was studied in flies originating from the Netherlands and eastern Australia. The hsp26 gene was polymorphic for an insertion/deletion of three extra amino acids and two nonsynonymous changes in all populations. The hsp27 gene exhibited two nonsynonymous changes and three synonymous mutations. The hsp26 polymorphism showed a latitudinal cline along the east coast of Australia. This pattern was not confounded by the fact that the shsps are located in the inversion In(3 L)P which also shows a latitudinal cline in eastern Australia. A similar latitudinal cline was found for the previously described variation in hsp23, while frequencies of hsp27 alleles did not change with latitude. These findings suggest that variation at two of the shsps or closely linked loci are under selection in natural populations of D. melanogaster.     

Analysis of mtDNA variation in African populations reveals the most ancient of all human continent-specific haplogroups.

mtDNA sequence variation was examined in 140 Africans, including Pygmies from Zaire and Central African Republic (C.A.R.) and Mandenkalu, Wolof, and Pular from Senegal. More than 76% of the African mtDNAs (100% of the Pygmies and 67.3% of the Senegalese) formed one major mtDNA cluster (haplogroup L) defined by an African-specific HpaI site gain at nucleotide pair (np) 3592. Additional mutations subdivided haplogroup L into two subhaplogroups, each encompassing both Pygmy and Senegalese mtDNAs. A novel 12-bp homoplasmic insertion in the intergenic region between tRNA(Tyr) and cytochrome oxidase I (COI) genes was also observed in 17.6% of the Pygmies from C.A.R. This insertion is one of the largest observed in human mtDNAs. Another 25% of the Pygmy mtDNAs harbored a 9-bp deletion between the cytochrome oxidase II (COII) and tRNA(Lys) genes, a length polymorphism previously reported in non-African populations. In addition to haplogroup L, other haplogroups were observed in the Senegalese. These haplogroups were more similar to those observed in Europeans and Asians than to haplogroup L mtDNAs, suggesting that the African mtDNAs without the HpaI np 3592 site could be the ancestral types from which European and Asian mtDNAs were derived. Comparison of the intrapopulation sequence divergence in African and non-African populations confirms that African populations exhibit the largest extent of mtDNA variation, a result that further supports the hypothesis that Africans represent the most ancient human group and that all modern humans have a common and recent African origin. The age of the total African variation was estimated to be 101,000-133,000 years before present (YBP), while the age of haplogroup L was estimated at 98,000-130,000 YBP. These values substantially exceed the ages of all Asian- and European-specific mtDNA haplogroups.     

Contrasting patterns of natural variation in global Drosophila melanogaster populations

Despite the popularity of Drosophila melanogaster in functional and evolutionary genetics, the global pattern of natural variation has not yet been comprehensively described in this species. For the first time, we report a combined survey using neutral microsatellites and mitochondrial sequence variation jointly. Thirty-five populations originating from five continents were compared. In agreement with previous microsatellite studies, sub-Saharan African populations were the most variable ones. Consistent with previous reports of a single 'out of Africa' habitat expansion, we found that non-African populations contained a subset of the African alleles. The pattern of variation detected for the mitochondrial sequences differed substantially. The most divergent haplotypes were detected in the Mediterranean region while Africa harboured most haplotypes, which were all closely related. In the light of the well-established African origin of D. melanogaster, our results cast severe doubts about the suitability of mtDNA for biogeographic inference in this model organism.     

Drosophila seminal fluid proteins enter the circulatory system of the mated female fly by crossing the posterior vaginal wall

Seminal fluid proteins from males of many insect species affect the behavior and physiology of their mates. In some cases, these effects result from entry of the proteins into the female's circulatory system. In the fruit fly Drosophila melanogaster, some seminal fluid proteins enter the female's circulatory system after transfer from the male while others remain confined within the reproductive tract. To address where and how seminal fluid proteins enter the hemolymph of the mated female, we compared the kinetics of transfer and localization in mated females of two seminal fluid proteins that enter the hemolymph (Acp26Aa and Acp62F) and one that does not (Acp36DE). We also generated transgenic flies that produce Acp26Aa tagged with Aequorea victoria green fluorescent protein (GFP) to monitor its transfer in vivo. We report that Acps enter the female circulatory system from the posterior vagina immediately after insemination. The ability of Acps to enter the female hemolymph correlates with their ability to cross the intima that lines the posterior vagina. The ventral posterior vagina is structurally unlike other parts of the female reproductive tract in that it lacks muscles. We hypothesize that it has higher permeability thus affording access to the female's circulatory system.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.