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.      

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 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.     

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].     

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 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 thepostmating reactions of the female, is a single-copy gene. However, inDrosophila subobscura, the gene is duplicated and both copies aretranscribed. To study the molecular evolution of the duplication, a 2.1- kbfragment encompassing both copies of the duplication was sequenced for 10lines of D. subobscura and one line of Drosophila madeirensis. Estimates ofthe divergence between the two copies of the duplicated region and betweenthe two species studied, D. subobscura and D. madeirensis, revealed thatboth copies of the Acp70a gene had evolved independently since theirduplication. The ratio of nonsynonymous to silent divergence between copieswas generally higher than one. The McDonald and Kreitman test revealed anexcess of nonsynonymous changes fixed since the duplication and before thesplit of the D. subobscura and D. madeirensis lineages. These results pointto natural selection driving protein evolution after the duplication.Specifically, adaptive evolution appears to have caused the initialdifferentiation between copies of the N-terminal parts of the proteins,while purifying selection could be responsible for the high conservation ofthe C- terminal parts.     

Positive selection of Iris, a retroviral envelope-derived host gene in Drosophila melanogaster

Eukaryotic genomes can usurp enzymatic functions encoded by mobile elements for their own use. A particularly interesting kind of acquisition involves the domestication of retroviral envelope genes, which confer infectious membrane-fusion ability to retroviruses. So far, these examples have been limited to vertebrate genomes, including primates where the domesticated envelope is under purifying selection to assist placental function. Here, we show that in Drosophila genomes, a previously unannotated gene (CG4715, renamed Iris) was domesticated from a novel, active Kanga lineage of insect retroviruses at least 25 million years ago, and has since been maintained as a host gene that is expressed in all adult tissues. Iris and the envelope genes from Kanga retroviruses are homologous to those found in insect baculoviruses and gypsy and roo insect retroviruses. Two separate envelope domestications from the Kanga and roo retroviruses have taken place, in fruit fly and mosquito genomes, respectively. Whereas retroviral envelopes are proteolytically cleaved into the ligand-interaction and membrane-fusion domains, Iris appears to lack this cleavage site. In the takahashii/suzukii species groups of Drosophila, we find that Iris has tandemly duplicated to give rise to two genes (Iris-A and Iris-B). Iris-B has significantly diverged from the Iris-A lineage, primarily because of the "invention" of an intron de novo in what was previously exonic sequence. Unlike domesticated retroviral envelope genes in mammals, we find that Iris has been subject to strong positive selection between Drosophila species. The rapid, adaptive evolution of Iris is sufficient to unambiguously distinguish the phylogenies of three closely related sibling species of Drosophila (D. simulans, D. sechellia, and D. mauritiana), a discriminative power previously described only for a putative "speciation gene." Iris represents the first instance of a retroviral envelope-derived host gene outside vertebrates. It is also the first example of a retroviral envelope gene that has been found to be subject to positive selection following its domestication. The unusual selective pressures acting on Iris suggest that it is an active participant in an ongoing genetic conflict. We propose a model in which Iris has "switched sides," having been recruited by host genomes to combat baculoviruses and retroviruses, which employ homologous envelope genes to mediate infection.     

Natural selection and molecular evolution in PTC, a bitter-taste receptor gene

The ability to taste phenylthiocarbamide (PTC) is a classic phenotype that has long been known to vary in human populations. This phenotype is of genetic, epidemiologic, and evolutionary interest because the ability to taste PTC is correlated with the ability to taste other bitter substances, many of which are toxic. Thus, variation in PTC perception may reflect variation in dietary preferences throughout human history and could correlate with susceptibility to diet-related diseases in modern populations. To test R. A. Fisher's long-standing hypothesis that variability in PTC perception has been maintained by balancing natural selection, we examined patterns of DNA sequence variation in the recently identified PTC gene, which accounts for up to 85% of phenotypic variance in the trait. We analyzed the entire coding region of PTC (1,002 bp) in a sample of 330 chromosomes collected from African (n=62), Asian (n=138), European (n=110), and North American (n=20) populations by use of new statistical tests for natural selection that take into account the potentially confounding effects of human population growth. Two intermediate-frequency haplotypes corresponding to "taster" and "nontaster" phenotypes were found. These haplotypes had similar frequencies across Africa, Asia, and Europe. Genetic differentiation between the continental population samples was low (FST=0.056) in comparison with estimates based on other genes. In addition, Tajima's D and Fu and Li's D and F statistics demonstrated a significant deviation from neutrality because of an excess of intermediate-frequency variants when human population growth was taken into account (P<.01). These results combine to suggest that balancing natural selection has acted to maintain "taster" and "nontaster" alleles at the PTC locus in humans.     

Drosophila chemoreceptor gene evolution: selection, specialization and genome size

Chemoperception plays a key role in adaptation and speciation in animals, and the senses of olfaction and gustation are mediated by gene families which show large variation in repertoire size among species. In Drosophila, there are around 60 loci of each type and it is thought that ecological specialization influences repertoire size, with increased pseudogenization of loci. Here, we analyse the size of the gustatory and olfactory repertoires among the genomes of 12 species of Drosophila . We find that repertoire size varies substantially and the loci are evolving by duplication and pseudogenization, with striking examples of lineage-specific duplication. Selection analyses imply that the majority of loci are subject to purifying selection, but this is less strong in gustatory loci and in loci prone to duplication. In contrast to some other studies, we find that few loci show statistically significant evidence of positive selection. Overall genome size is strongly correlated with the proportion of duplicated chemoreceptor loci, but genome size, specialization and endemism may be interrelated in their influence on repertoire size.     

Positive selection for the male functionality of a co-retroposed gene in the hominoids

Background  

New genes generated by retroposition are widespread in humans and other mammalian species. Usually, this process copies a single parental gene and inserts it into a distant genomic location. However, retroposition of two adjacent parental genes, i.e. co-retroposition, had not been reported until the hominoid chimeric gene, PIPSL, was identified recently. It was shown how two genes linked in tandem (phosphatidylinositol-4-phosphate 5-kinase, type I, alpha, PIP5K1A and proteasome 26S subunit, non-ATPase, 4, PSMD4) could be co-retroposed from a single RNA molecule to form this novel chimeric gene. However, understanding of the origination and biological function of PIPSL requires determination of the coding potential of this gene as well as the evolutionary forces acting on its hominoid copies.     

Positive selection drives the evolution of rhino, a member of the heterochromatin protein 1 family in Drosophila

Heterochromatin comprises a significant component of many eukaryotic genomes. In comparison to euchromatin, heterochromatin is gene poor, transposon rich, and late replicating. It serves many important biological roles, from gene silencing to accurate chromosome segregation, yet little is known about the evolutionary constraints that shape heterochromatin. A complementary approach to the traditional one of directly studying heterochromatic DNA sequence is to study the evolution of proteins that bind and define heterochromatin. One of the best markers for heterochromatin is the heterochromatin protein 1 (HP1), which is an essential, nonhistone chromosomal protein. Here we investigate the molecular evolution of five HP1 paralogs present in Drosophila melanogaster. Three of these paralogs have ubiquitous expression patterns in adult Drosophila tissues, whereas HP1D/rhino and HP1E are expressed predominantly in ovaries and testes respectively. The HP1 paralogs also have distinct localization preferences in Drosophila cells. Thus, Rhino localizes to the heterochromatic compartment in Drosophila tissue culture cells, but in a pattern distinct from HP1A and lysine-9 dimethylated H3. Using molecular evolution and population genetic analyses, we find that rhino has been subject to positive selection in all three domains of the protein: the N-terminal chromo domain, the C-terminal chromo-shadow domain, and the hinge region that connects these two modules. Maximum likelihood analysis of rhino sequences from 20 species of Drosophila reveals that a small number of residues of the chromo and shadow domains have been subject to repeated positive selection. The rapid and positive selection of rhino is highly unusual for a gene encoding a chromosomal protein and suggests that rhino is involved in a genetic conflict that affects the germline, belying the notion that heterochromatin is simply a passive recipient of "junk DNA" in eukaryotic genomes.     

Phenotypic plasticity and selection in Drosophila life-history evolution. I. Nutrition and the cost of reproduction

Earlier experiments have shown that the evolution of postponed senescent populations can be achieved by selection on either demographic or stress resistance characters. Both types of selection have produced results in which survival characters (stress resistance and longevity) have apparently traded-off against early-life fecundity. Here we present the results of a series of experiments in which an environmental variable — the level of live yeast inoculate applied to the substrate — produces a qualitatively similar phenotypic response: longevity and starvation resistance are enhanced by lower yeast levels, at the expense of fecundity. For the starvation resistance versus fecundity experiments we show a negative and linear relationship between the norms of reaction for each character across a gradient of yeast levels. This phenotypic trade-off is stable across the 20 populations and 4 selection treatments reported on here, and its general agreement with earlier selection results suggests that the evolutionary response and the phenotypically plastic response may share a common physiological basis. However, an important discrepancy in the lifetime fecundity data between the selection response and the dietary manipulations preclude strict analogy. The results broadly conform to a simple “Y-model” of allocation, in which a limited resource is divided between survival and reproduction; here the characters are starvation resistance and longevity versus fecundity.     

Translational selection and molecular evolution

An interplay among experimental studies of protein synthesis, evolutionary theory, and comparisons of DNA sequence data has shed light on the roles of natural selection and genetic drift in ‘silent’ DNA evolution.     

Frequency-dependent selection, metrical characters and molecular evolution

Computer models of selection acting on a quantitative character show that a combination of frequency-dependent and stabilizing selection can maintain many polymorphisms among the genes that determine the character. The models also show that the random order of mutations can give rise to selectively driven stochastic effects that are sometimes more important than random genetic drift. They suggest simple explanations for patterns of divergence between populations and species, and for apparent discrepancies between the rates of morphological and molecular evolution. They point towards a selective theory of 'molecular clocks'.     

Dntf-2r,a young Drosophila retroposed gene with specific male expression under positive Darwinian selection

A direct approach to investigating new gene origination is to examine recently evolved genes. We report a new gene in the Drosophila melanogaster subgroup, Drosophila nuclear transport factor-2-related (Dntf-2r). Its sequence features and phylogenetic distribution indicate that Dntf-2r is a retroposed functional gene and originated in the common ancestor of D. melanogaster, D. simulans, D. sechellia, and D. mauritiana, within the past 3-12 million years (MY). Dntf-2r evolved more rapidly than the parental gene, under positive Darwinian selection as revealed by the McDonald-Kreitman test and other evolutionary analyses. Comparative expression analysis shows that Dntf-2r is male specific whereas the parental gene, Dntf-2, is widely expressed in D. melanogaster. In agreement with its new expression pattern, the Dntf-2r putative promoter sequence is similar to the late testis promoter of beta2-tubulin. We discuss the possibility that the action of positive selection in Dntf-2r is related to its putative male-specific functions.     

Positive selection along the evolution of primate mitogenomes

The mitochondrial genomes of four neotropical primates, Aotus infulatus, Chiropotes israelita, Callimico goeldii and Callicebus lugens were sequenced and annotated. Phylogenetic reconstructions with mitochondrial genes of other 66 primates showed a similar arrangement to a topology based on nuclear genes. Screening for positive selection identified 15 codons in 7 genes along 9 independent lineages, three with two or more genes and five in internal nodes, ruling out false positive estimates. Mitochondrial genes of the electron transport chain (ETC.) complexes evolved with high substitution rates. A study of nuclear ETC. genes might elucidate whether they co-evolved with their mitochondrial counterparts.     

Positive selection, relaxation, and acceleration in the evolution of the human and chimp genome

For years evolutionary biologists have been interested in searching for the genetic bases underlying humanness. Recent efforts at a large or a complete genomic scale have been conducted to search for positively selected genes in human and in chimp. However, recently developed methods allowing for a more sensitive and controlled approach in the detection of positive selection can be employed. Here, using 13,198 genes, we have deduced the sets of genes involved in rate acceleration, positive selection, and relaxation of selective constraints in human, in chimp, and in their ancestral lineage since the divergence from murids. Significant deviations from the strict molecular clock were observed in 469 human and in 651 chimp genes. The more stringent branch-site test of positive selection detected 108 human and 577 chimp positively selected genes. An important proportion of the positively selected genes did not show a significant acceleration in rates, and similarly, many of the accelerated genes did not show significant signals of positive selection. Functional differentiation of genes under rate acceleration, positive selection, and relaxation was not statistically significant between human and chimp with the exception of terms related to G-protein coupled receptors and sensory perception. Both of these were over-represented under relaxation in human in relation to chimp. Comparing differences between derived and ancestral lineages, a more conspicuous change in trends seems to have favored positive selection in the human lineage. Since most of the positively selected genes are different under the same functional categories between these species, we suggest that the individual roles of the alternative positively selected genes may be an important factor underlying biological differences between these species.