DNA variability and divergence at the notch locus in Drosophila melanogaster and D. simulans: a case of accelerated synonymous site divergence

DNA diversity in two segments of the Notch locus was surveyed in four populations of Drosophila melanogaster and two of D. simulans. In both species we observed evidence of non-steady-state evolution. In D. simulans we observed a significant excess of intermediate frequency variants in a non-African population. In D. melanogaster we observed a disparity between levels of sequence polymorphism and divergence between one of the Notch regions sequenced and other neutral X chromosome loci. The striking feature of the data is the high level of synonymous site divergence at Notch, which is the highest reported to date. To more thoroughly investigate the pattern of synonymous site evolution between these species, we developed a method for calibrating preferred, unpreferred, and equal synonymous substitutions by the effective (potential) number of such changes. In D. simulans, we find that preferred changes per "site" are evolving significantly faster than unpreferred changes at Notch. In contrast we observe a significantly faster per site substitution rate of unpreferred changes in D. melanogaster at this locus. These results suggest that positive selection, and not simply relaxation of constraint on codon bias, has contributed to the higher levels of unpreferred divergence along the D. melanogaster lineage at Notch.     

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.     

Recurrent positive selection at bgcn, a key determinant of germ line differentiation, does not appear to be driven by simple coevolution with its partner protein bam

Surveys of nucleotide sequence polymorphism in Drosophila melanogaster and Drosophila simulans were performed at 2 interacting loci crucial for gametogenesis: bag-of-marbles (bam) and benign gonial cell neoplasm (bgcn). At the polymorphism level, both loci appear to be evolving under the expectations of the neutral theory. However, ratios of polymorphism and divergence for synonymous and nonsynonymous mutations depart significantly from neutral expectations for both loci consistent with a previous observation of positive selection at bam. The deviations suggest either an excess of synonymous polymorphisms or an excess of nonsynonymous fixations at both loci. Synonymous evolution appears to conform to neutrality at bam. At bgcn, there is evidence of positive selection affecting preferred synonymous mutations along the D. simulans lineage. However, there is also a significantly higher rate of nonsynonymous fixations at bgcn within D. simulans. Thus, the deviation from neutrality detected by the McDonald-Kreitman test at these 2 loci is likely due to the selective acceleration of nonsynonymous fixations. Differences in the pattern of amino acid fixations between these 2 interacting proteins suggest that the detected positive selection is not due to a simple model of coevolution.     

Genomewide spatial correspondence between nonsynonymous divergence and neutral polymorphism reveals extensive adaptation in Drosophila

The effect of recurrent selective sweeps is a spatially heterogeneous reduction in neutral polymorphism throughout the genome. The pattern of reduction depends on the selective advantage and recurrence rate of the sweeps. Because many adaptive substitutions responsible for these sweeps also contribute to nonsynonymous divergence, the spatial distribution of nonsynonymous divergence also reflects the distribution of adaptive substitutions. Thus, the spatial correspondence between neutral polymorphism and nonsynonymous divergence may be especially informative about the process of adaptation. Here we study this correspondence using genomewide polymorphism data from Drosophila simulans and the divergence between D. simulans and D. melanogaster. Focusing on highly recombining portions of the autosomes, at a spatial scale appropriate to the study of selective sweeps, we find that neutral polymorphism is both lower and, as measured by a new statistic Q(S), less homogeneous where nonsynonymous divergence is higher and that the spatial structure of this correlation is best explained by the action of strong recurrent selective sweeps. We introduce a method to infer, from the spatial correspondence between polymorphism and divergence, the rate and selective strength of adaptation. Our results independently confirm a high rate of adaptive substitution (approximately 1/3000 generations) and newly suggest that many adaptations are of surprisingly great selective effect (approximately 1%), reducing the effective population size by approximately 15% even in highly recombining regions of the genome.     

Positive and negative selection on noncoding DNA in Drosophila simulans

There is now a wealth of evidence that some of the most important regions of the genome are found outside those that encode proteins, and noncoding regions of the genome have been shown to be subject to substantial levels of selective constraint, particularly in Drosophila. Recent work has suggested that these regions may also have been subject to the action of positive selection, with large fractions of noncoding divergence having been driven to fixation by adaptive evolution. However, this work has focused on Drosophila melanogaster, which is thought to have experienced a reduction in effective population size (N(e)), and thus a reduction in the efficacy of selection, compared with its closest relative Drosophila simulans. Here, we examine patterns of evolution at several classes of noncoding DNA in D. simulans and find that all noncoding DNA is subject to the action of negative selection, indicated by reduced levels of polymorphism and divergence and a skew in the frequency spectrum toward rare variants. We find that the signature of negative selection on noncoding DNA and nonsynonymous sites is obscured to some extent by purifying selection acting on preferred to unpreferred synonymous codon mutations. We investigate the extent to which divergence in noncoding DNA is inferred to be the product of positive selection and to what extent these inferences depend on selection on synonymous sites and demography. Based on patterns of polymorphism and divergence for different classes of synonymous substitution, we find the divergence excess inferred in noncoding DNA and nonsynonymous sites in the D. simulans lineage difficult to reconcile with demographic explanations.     

Inferring Weak Selection from Patterns of Polymorphism and Divergence at ``silent' Sites in Drosophila DNA

Patterns of codon usage and ``silent'''' DNA divergence suggest that natural selection discriminates among synonymous codons in Drosophila. ``Preferred'''' codons are consistently found in higher frequencies within their synonymous families in Drosophila melanogaster genes. This suggests a simple model of silent DNA evolution where natural selection favors mutations from unpreferred to preferred codons (preferred changes). Changes in the opposite direction, from preferred to unpreferred synonymous codons (unpreferred changes), are selected against. Here, selection on synonymous DNA mutations is investigated by comparing the evolutionary dynamics of these two categories of silent DNA changes. Sequences from outgroups are used to determine the direction of synonymous DNA changes within and between D. melanogaster and Drosophila simulans for five genes. Population genetics theory shows that differences in the fitness effect of mutations can be inferred from the comparison of ratios of polymorphism to divergence. Unpreferred changes show a significantly higher ratio of polymorphism to divergence than preferred changes in the D. simulans lineage, confirming the action of selection at silent sites. An excess of unpreferred fixations in 28 genes suggests a relaxation of selection on synonymous mutations in D. melanogaster. Estimates of selection coefficients for synonymous mutations (3.6 <|N(e)s| < 1.3) in D. simulans are consistent with the reduced efficacy of natural selection (|N(e)s| < 1) in the three- to sixfold smaller effective population size of D. melanogaster. Synonymous DNA changes appear to be a prevalent class of weakly selected mutations in Drosophila.     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

The evolution of antifungal peptides in Drosophila

An essential component of the immune system of animals is the production of antimicrobial peptides (AMPs). In vertebrates and termites the protein sequence of some AMPs evolves rapidly under positive selection, suggesting that they may be coevolving with pathogens. However, antibacterial peptides in Drosophila tend to be highly conserved. We have inferred the selection pressures acting on Drosophila antifungal peptides (drosomycins) from both the divergence of drosomycin genes within and between five species of Drosophila and polymorphism data from Drosophila simulans and D. melanogaster. In common with Drosophila antibacterial peptides, there is no evidence of adaptive protein evolution in any of the drosomycin genes, suggesting that they do not coevolve with pathogens. It is possible that this reflects a lack of specific fungal and bacterial parasites in Drosophila populations. The polymorphism data from both species differed from neutrality at one locus, but this was not associated with changes in the protein sequence. The synonymous site diversity was greater in D. simulans than in D. melanogaster, but the diversity both upstream of the genes and at nonsynonymous sites was similar. This can be explained if both upstream and nonsynonymous mutations are slightly deleterious and are removed more effectively from D. simulans due to its larger effective population size.     

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.     

Lack of Polymorphism on the Drosophila Fourth Chromosome Resulting from Selection

Evolutionary processes can be inferred from comparisons of intraspecific polymorphism and interspecific divergence. We sequenced a 1.1-kb fragment of the cubitus interruptus Dominant (ciD) locus located on the nonrecombining fourth chromosome for ten natural lines of Drosophila melanogaster and nine of Drosophila simulans. We found no polymorphism within D. melanogaster and a single polymorphism within D. simulans; divergence between the species was about 5%. Comparison with the alcohol dehydrogenase gene and its two flanking regions in D. melanogaster, for which comparable data are available, revealed a statistically significant departure from neutrality in all three tests. This lack of polymorphism in the ciD locus may reflect recent positive selective sweeps on the fourth chromosome with extreme hitchhiking generated by the lack of recombination. By simulation, we estimate there to be a 50% chance that the selective sweeps occurred within the past 30,000 years in D. melanogaster and 75,000 in D. simulans.     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Molecular evolution of X-linked accessory gland proteins in Drosophila pseudoobscura

In Drosophila melanogaster and Drosophila simulans, positive Darwinian selection drives high rates of evolution of male reproductive genes, and accessory gland proteins (Acps) in particular. Here, we tested whether 13 X-linked male-specific genes, 4 Acps and 9 non-Acps, are under selective forces in the Drosophila pseudoobscura species group, much as those in the D. melanogaster group. We observed a statistically significant correlation in relative rates of nonsynonymous evolution between the two species groups tested. One Acp examined had a higher rate of nonsynonymous substitution than predicted by a neutral model in both species groups, suggesting its divergence was driven by positive Darwinian selection. To further test for the signature of selection, we examined polymorphism of three Acps within D. pseudoobscura. From this test, no Acp individually bore the signature of positive selection, but the 3 Acps together possessed an excess of nonsynonymous differences between species, relative to polymorphism within species. We conclude that faster evolution of Acps in the D. pseudoobscura group appears to be driven by positive selection, as previously suggested in the D. melanogaster group.     

Molecular population genetics of ref(2)P, a locus which confers viral resistance in Drosophila

The ref(2)P locus (2-54.2) is polymorphic for two allelic forms in natural populations of Drosophila melanogaster, ref(2)Po and ref(2)Pp. The latter allele confers resistance to the rhabdovirus sigma infecting wild populations. Previous work, based on a small sample of prescreened restrictive (resistant) and permissive (susceptible) alleles, identified a large number of amino acid replacement changes (7) relative to synonymous changes (1). Such protein variability could be the result of variation-enhancing selection. To further test the selection hypothesis, we have examined the DNA sequences of ten randomly chosen lines of D. melanogaster and one line of D. simulans. Nine of the ten lines are permissive; D. simulans does not harbor the virus. The melanogaster alleles contain 4 synonymous changes, 19 noncoding changes, and 13 amino acid replacement changes, indicating a relatively high level of polymorphism. Three sequenced restrictive alleles have nearly identical sequences, indicating that they are relatively young. Compared to the permissive alleles, they share only a complex deletion at codon 34, CAG-AAT to GGA, which our analysis indicates to be the site conferring the restrictive phenotype. Patterns of polymorphism and divergence differ from neutral predictions by several criteria for the amino terminal region, which contains the complex deletion (codons 1-91), but not the remainder of the protein (codons 92-599). We find a higher rate of evolution on the D. melanogaster lineage than on the D. simulans lineage. The relatively large amount of both replacement and silent polymorphism in the permissive alleles and the lack of divergence between permissive and restrictive alleles suggests that the sigma virus and ref(2)P may be engaged in an evolutionary race in which new restrictive alleles are continually arising but are relatively short-lived.      

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.     

Synonymous rates at the RpII215 gene of Drosophila: variation among species and across the coding region

The region encompassing the RpII215 gene that encodes the largest component of the RNA polymerase II complex (1889 amino acids) has been sequenced in Drosophila subobscura, D. madeirensis, D. guanche, and D. pseudoobscura. Nonsynonymous divergence estimates (Ka) indicate that this gene has a very low rate of amino acid replacements. Given its low Ka and constitutive expression, synonymous substitution rates are, however, unexpectedly high. Sequence comparisons have allowed the molecular clock hypothesis to be tested. D. guanche is an insular species and it is therefore expected to have a reduced effective size relative to D. subobscura. The significantly higher rate of synonymous substitutions detected in the D. guanche lineage could be explained if synonymous mutations behave as nearly neutral. Significant departure from the molecular clock hypothesis for synonymous and nonsynonymous substitutions was detected when comparing the D. subobscura, D. pseudoobscura, and D. melanogaster lineages. Codon bias and synonymous divergence between D. subobscura and D. melanogaster were negatively correlated across the RpII215 coding region, which indicates that selection coefficients for synonymous mutations vary across the gene. The C-terminal domain (CTD) of the RpII215 protein is structurally and functionally differentiated from the rest of the protein. Synonymous substitution rates were significantly different in both regions, which strongly indicates that synonymous mutations in the CTD and in the non-CTD regions are under detectably different selection coefficients.     

High divergence among Drosophila simulans mitochondrial haplogroups arose in midst of long term purifying selection

We characterize the type of selection acting within and among mitochondrial lineages in five closely related Drosophila species. We focus on D. simulans, where three genetically distinct mitochondrial haplogroups show high interhaplogroup divergence and low intrahaplogroup polymorphism. Using maximum likelihood models we find that the branches leading to these three distinct mitochondrial groups show a significantly reduced rate of nonsynonymous relative to synonymous substitution. This interhaplogroup rate is significantly reduced compared to the intrahaplogroup rate, and closely resembles the rate observed between distinct species. The data suggest that slightly deleterious mutations segregating within D. simulans haplogroups are removed by selection prior to their fixation among haplogroups. We explore several hypotheses to explain how lineages within a single species can be compatible with this model of slightly deleterious mutation. The most likely hypothesis is that D. simulans haplogroups have persisted in isolation, perhaps due to association with the bacterial symbiont Wolbachia and/or demographic history, introducing a bias against the fixation of slightly deleterious mutations.     

Rapid evolution and gene-specific patterns of selection for three genes of spermatogenesis in Drosophila

Hybrid males resulting from crosses between closely related species of Drosophila are sterile. The F1 hybrid sterility phenotype is mainly due to defects occurring during late stages of development that relate to sperm individualization, and so genes controlling sperm development may have been subjected to selective diversification between species. It is also possible that genes of spermatogenesis experience selective constraints given their role in a developmental pathway. We analyzed the molecular evolution of three genes playing a role during the sperm developmental pathway in Drosophila at an early (bam), a mid (aly), and a late (dj) stage. The complete coding region of these genes was sequenced in different strains of Drosophila melanogaster and Drosophila simulans. All three genes showed rapid divergence between species, with larger numbers of nonsynonymous to synonymous differences between species than polymorphisms. Although this could be interpreted as evidence for positive selection at all three genes, formal tests of selection do not support such a conclusion. Departures from neutrality were detected only for dj and bam but not aly. The role played by selection is unique and determined by gene-specific characteristics rather than site of expression. In dj, the departure was due to a high proportion of neutral synonymous polymorphisms in D. simulans, and there was evidence of purifying selection maintaining a high lysine amino acid protein content that is characteristic of other DNA-binding proteins. The earliest spermatogenesis gene surveyed, which plays a role in both male and female gametogenesis, was bam, and its significant departure from neutrality was due to an excess of nonsynonymous substitutions between species. Bam is degraded at the end of mitosis, and rapid evolutionary changes among species might be a characteristic shared with other degradable transient proteins. However, the large number of nonsynonymous changes between D. melanogaster and D. simulans and a phylogenetic comparative analysis among species confirms evidence of positive selection driving the evolution of Bam and suggests an yet unknown germ cell line developmental adaptive change between these two species.     

The impact of shared ancestral variation on hybrid male lethality--a 16 codon indel in the Drosophila simulans Lhr gene

The Lethal hybrid rescue (Lhr) gene causes hybrid male lethality in crosses between Drosophila simulans and D. melanogaster. Lhr(2) is a D. simulans allele, which rescues hybrid males. It has been recently proposed that a 16 codon insertion, which distinguishes the D. melanogaster and the canonical D. simulans allele, and is lacking in Lhr(2), may be responsible for the functional divergence of D. melanogaster and D. simulans Lhr alleles. Here, we show that the Lhr(2) allele lacking the insertion represents an ancestral polymorphism segregating at a moderate frequency in D. simulans. Crosses of D. melanogaster females to males from two D. simulans strains carrying this deletion showed a severe deficiency of viable hybrid males. Our results suggest that the absence of this insertion alone is not sufficient to explain functional differences between D. melanogaster and D. simulans Lhr alleles.