High nucleotide sequence variation in a region of low recombination in Drosophila simulans is consistent with the background selection model

We surveyed nucleotide sequence variation at glucose dehydrogenase (Gld), in a region of low recombination on chromosome 3R, from a population sample of Drosophila simulans. The levels of nucleotide variation were surprisingly high. There was no departure from the expectation of a neutral model for the level of polymorphism, indicating no evidence of a selective sweep in this region. There was a significant deficiency of singleton polymorphisms according to the Fu and Li test, although Tajima and Hudson, Kreitman, and Aguade (HKA) tests do not provide evidence of a significant elevation of variation due to balancing selection. Genetic map data for the D. simulans third chromosome were used to calculate expected values of pi for Gld under a current model of background selection, varying the values for the parameter sh (selection coefficient against deleterious mutations). We show that the recombinational landscape of D. simulans is sufficiently different from that of D. melanogaster that we expect higher variation under the background selection model, even when effective population sizes are assumed to be equal. The data for Gld were tested against the predictions using computer simulations of the distribution of the number of segregating sites conditioned on pi. Background selection alone can explain our observations as long as sh is larger than 0.005 and species-level effective population size is assumed to be several- fold larger than in D. melanogaster. Alternatively, the deleterious mutation rate may be smaller in D. simulans, or balancing selection may be acting nearby, thereby reducing the effect of background selection.      

Nucleotide sequence divergence in the A+T-rich region of mitochondrial DNA in Drosophila simulans and Drosophila mauritiana

We have determined the nucleotide sequences of two regions within the A+T-rich region of mitochondrial DNA (mtDNA) in the siIII type of Drosophila simulans and the maI type of D. mauritiana. The sequences of the two regions in siIII and maI are almost identical. The sequences include elements corresponding to the type I and type II repeats elements and the T-stretches as reported in D. melanogaster; an approximately 340-bp region (A region) adjacent to the tRNA(Ile) gene includes a part of the type II repeat element, and an approximately 440- bp region (B region) includes a central portion of the A+T-rich region between the type I and type II repeat arrays. Each sequence of the two species was compared with those of D. melanogaster and D. yakuba. The sequences of the A region are relatively well conserved among the four species. The alignment of the two sequences of the B region with those of D. melanogaster and D. yakuba requires numerous insertions/deletions. For both regions, nucleotide differences between D. simulans or D. mauritiana and D. melanogaster are similar to those between the two and D. yakuba. The tendency is obvious in a subregion within the type II repeat element in the A region. These findings suggest that the rate of nucleotide substitution in the subregion is accelerated in the lineage leading to D. melanogaster. Loss of functional constraint in the stem-loop-forming sequence is proposed for this acceleration.      

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.     

Linkage disequilibrium and recent selection at three immunity receptor loci in Drosophila simulans

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations.     

Genomic Patterns of Geographic Differentiation in Drosophila simulans

Geographic patterns of genetic differentiation have long been used to understand population history and to learn about the biological mechanisms of adaptation. Here we present an examination of genomic patterns of differentiation between northern and southern populations of Australian and North American Drosophila simulans, with an emphasis on characterizing signals of parallel differentiation. We report on the genomic scale of differentiation and functional enrichment of outlier SNPs. While, overall, signals of shared differentiation are modest, we find the strongest support for parallel differentiation in genomic regions that are associated with regulation. Comparisons to Drosophila melanogaster yield potential candidate genes involved in local adaptation in both species, providing insight into common selective pressures and responses. In contrast to D. melanogaster, in D. simulans we observe patterns of variation that are inconsistent with a model of temperate adaptation out of a tropical ancestral range, highlighting potential differences in demographic and colonization histories of this cosmopolitan species pair.     

DNA variation at the rp49 gene region of Drosophila simulans: evolutionary inferences from an unusual haplotype structure.

An approximately 1.3-kb region including the rp49 gene plus its 5' and 3' flanking regions was sequenced in 24 lines of Drosophila simulans (10 from Spain and 14 from Mozambique). Fifty-four nucleotide and 8 length polymorphisms were detected. All nucleotide polymorphisms were silent: 52 in noncoding regions and 2 at synonymous sites in the coding region. Estimated silent nucleotide diversity was similar in both populations (pi = 0.016, for the total sample). Nucleotide variation revealed an unusual haplotype structure showing a subset of 11 sequences with a single polymorphism. This haplotype was present at intermediate frequencies in both the European and the African samples. The presence of such a major haplotype in a highly recombining region is incompatible with the neutral equilibrium model. This haplotype structure in both a derived and a putatively ancestral population can be most parsimoniously explained by positive selection. As the rate of recombination in the rp49 region is high, the target of selection should be close to or within the region studied.     

Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans

Background

Understanding how DNA sequence polymorphism relates to variation in gene expression is essential to connecting genotypic differences with phenotypic differences among individuals. Addressing this question requires linking population genomic data with gene expression variation.

Results

Using whole genome expression data and recent light shotgun genome sequencing of six Drosophila simulans genotypes, we assessed the relationship between expression variation in males and females and nucleotide polymorphism across thousands of loci. By examining sequence polymorphism in gene features, such as untranslated regions and introns, we find that genes showing greater variation in gene expression between genotypes also have higher levels of sequence polymorphism in many gene features. Accordingly, X-linked genes, which have lower sequence polymorphism levels than autosomal genes, also show less expression variation than autosomal genes. We also find that sex-specifically expressed genes show higher local levels of polymorphism and divergence than both sex-biased and unbiased genes, and that they appear to have simpler regulatory regions.

Conclusion

The gene-feature-based analyses and the X-to-autosome comparisons suggest that sequence polymorphism in cis-acting elements is an important determinant of expression variation. However, this relationship varies among the different categories of sex-biased expression, and trans factors might contribute more to male-specific gene expression than cis effects. Our analysis of sex-specific gene expression also shows that female-specific genes have been overlooked in analyses that only point to male-biased genes as having unusual patterns of evolution and that studies of sexually dimorphic traits need to recognize that the relationship between genetic and expression variation at these traits is different from the genome as a whole.     

Patterns of sequence variation in the mitochondrial D-loop region of shrews

Direct sequencing of the mitochondrial displacement loop (D-loop) of shrews (genus Sorex) for the region between the tRNA(Pro) and the conserved sequence block-F revealed variable numbers of 79-bp tandem repeats. These repeats were found in all 19 individuals sequenced, representing three subspecies and one closely related species of the masked shrew group (Sorex cinereus cinereus, S. c. miscix, S. c. acadicus, and S. haydeni) and an outgroup, the pygmy shrew (S. hoyi). Each specimen also possessed an adjacent 76-bp imperfect copy of the tandem repeats. One individual was heteroplasmic for length variants consisting of five and seven copies of the 79-bp tandem repeat. The sequence of the repeats is conducive to the formation of secondary structure. A termination-associated sequence is present in each of the repeats and in a unique sequence region 5' to the tandem array as well. Mean genetic distance between the masked shrew taxa and the pygmy shrew was calculated separately for the unique sequence region, one of the tandem repeats, the imperfect repeat, and these three regions combined. The unique sequence region evolved more rapidly than the tandem repeats or the imperfect repeat. The small genetic distance between pairs of tandem repeats within an individual is consistent with a model of concerted evolution. Repeats are apparently duplicated and lost at a high rate, which tends to homogenize the tandem array. The rate of D- loop sequence divergence between the masked and pygmy shrews is estimated to be 15%-20%/Myr, the highest rate observed in D-loops of mammals. Rapid sequence evolution in shrews may be due either to their high metabolic rate and short generation time or to the presence of variable numbers of tandem repeats.      

Organization of the sex-ratio meiotic drive region in Drosophila simulans

Sex-ratio meiotic drive is the preferential transmission of the X chromosome by XY males, which occurs in several Drosophila species and results in female-biased progeny. Although the trait has long been known to exist, its molecular basis remains completely unknown. Here we report a fine-mapping experiment designed to characterize the major drive locus on a sex-ratio X chromosome of Drosophila simulans originating from the Seychelles (XSR6). This primary locus was found to contain two interacting elements at least, both of which are required for drive expression. One of them was genetically tracked to a tandem duplication containing six annotated genes (Trf2, CG32712, CG12125, CG1440, CG12123, org-1), and the other to a candidate region located approximately 110 kb away and spanning seven annotated genes. RT-PCR showed that all but two of these genes were expressed in the testis of both sex-ratio and standard males. In situ hybridization to polytene chromosomes revealed a complete association of the duplication with the sex-ratio trait in random samples of X chromosomes from Madagascar and Reunion.     

Testing models of selection and demography in Drosophila simulans

We analyze patterns of nucleotide variability at 15 X-linked loci and 14 autosomal loci from a North American population of Drosophila simulans. We show that there is significantly more linkage disequilibrium on the X chromosome than on chromosome arm 3R and much more linkage disequilibrium on both chromosomes than expected from estimates of recombination rates, mutation rates, and levels of diversity. To explore what types of evolutionary models might explain this observation, we examine a model of recurrent, nonoverlapping selective sweeps and a model of a recent drastic bottleneck (e.g., founder event) in the demographic history of North American populations of D. simulans. The simple sweep model is not consistent with the observed patterns of linkage disequilibrium nor with the observed frequencies of segregating mutations. Under a restricted range of parameter values, a simple bottleneck model is consistent with multiple facets of the data. While our results do not exclude some influence of selection on X vs. autosome variability levels, they suggest that demography alone may account for patterns of linkage disequilibrium and the frequency spectrum of segregating mutations in this population of D. simulans.     

Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the approximately 5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was incompatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of D. simulans.     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans.

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Patterns of sequence variability and divergence at the diminutive gene region of Drosophila melanogaster: complex patterns suggest an ancestral selective sweep

To identify putatively swept regions of the Drosophila melanogaster genome, we performed a microsatellite screen spanning a 260-kb region of the X chromosome in populations from Zimbabwe, Ecuador, the United States, and China. Among the regions identified by this screen as showing a complex pattern of reduced heterozygosity and a skewed frequency spectrum was the gene diminutive (dm). To investigate the microsatellite findings, nucleotide sequence polymorphism data were generated in populations from both China and Zimbabwe spanning a 25-kb region and encompassing dm. Analysis of the sequence data reveals strongly reduced nucleotide variation across the entire gene region in both the non-African and the African populations, an extended haplotype pattern, and structured linkage disequilibrium, as well as a rejection of neutrality in favor of selection using a composite likelihood-ratio test. Additionally, unusual patterns of synonymous site evolution were observed at the second exon of this locus. On the basis of simulation studies as well as recently proposed methods for distinguishing between selection and nonequilibrium demography, we find that this "footprint" is best explained by a selective sweep in the ancestral population, the signal of which has been somewhat blurred via founder effects in the non-African samples.     

Rates of DNA sequence evolution are not sex-biased in Drosophila melanogaster and D. simulans

To determine whether male- or female-biased mutation rates have affected the molecular evolution of Drosophila melanogaster and D. simulans, we calculated the male-to-female ratio of germline cell divisions ([symbol: see text]) from germline generation data and the male-to-female ratio of mutation rate ([symbol: see text]) by comparing chromosomal levels of nucleotide divergence. We found that the ratio of germline cell divisions changes from indicating a weak female bias to indicating a weak male bias as the age of reproduction increases. The range of [symbol: see text] values that we observed, however, does not lead us to expect much, if any, difference in mutation rate between the sexes. Silent and intron nucleotide divergence were compared between nine loci on the X chromosome and nine loci on the second and third chromosomes. The average levels of nucleotide divergence were not significantly different across the chromosomes, although both silent and intron sites show a trend toward slightly more divergence on the X. These results indicate a lack of sex- or chromosome-biased molecular evolution in D. melanogaster and D. simulans.      

The frequency distribution of nucleotide variation in Drosophila simulans.

Patterns of codon bias in Drosophila suggest that silent mutations can be classified into two types: unpreferred (slightly deleterious) and preferred (slightly beneficial). Results of previous analyses of polymorphism and divergence in Drosophila simulans were interpreted as supporting a mutation-selection-drift model in which slightly deleterious, silent mutants make significantly greater contributions to polymorphism than to divergence. Frequencies of unpreferred polymorphisms were inferred to be lower than frequencies of other silent polymorphisms. Here, I analyzed additional D. simulans data to reevaluate the support for these ideas. I found that D. simulans has fixed more unpreferred than preferred mutations, suggesting that this lineage has not been at mutation-selection-drift equilibrium at silent sites. Frequencies of polarized unpreferred polymorphisms are not skewed toward rare alleles. However, frequencies of unpolarized unpreferred codons are lower in high-bias genes than in low-bias genes. This supports the idea that unpreferred codons are borderline deleterious mutations. Purifying selection on silent sites appears to be stronger at twofold-degenerate codons than at fourfold-degenerate codons. Finally, I found that X-linked polymorphisms occur at a higher average frequency than polymorphisms on chromosome arm 3R, even though an average X-linked site is significantly less likely to be polymorphic than an average site on 3R. This result supports a previous analysis of D. simulans indicating different population genetics of X-linked versus autosomal mutations.     

Identification of a locus under complex positive selection in Drosophila simulans by haplotype mapping and composite-likelihood estimation

The recent action of positive selection is expected to influence patterns of intraspecific DNA sequence variation in chromosomal regions linked to the selected locus. These effects include decreased polymorphism, increased linkage disequilibrium, and an increased frequency of derived variants. These effects are all expected to dissipate with distance from the selected locus due to recombination. Therefore, in regions of high recombination, it should be possible to localize a target of selection to a relatively small interval. Previously described patterns of intraspecific variation in three tandemly arranged, testes-expressed genes (janusA, janusB, and ocnus) in Drosophila simulans included all three of these features. Here we expand the original sample and also survey nucleotide polymorphism at three neighboring loci. On the basis of recombination events between derived and ancestral alleles, we localize the target of selection to a 1.5-kb region surrounding janusB. A composite-likelihood-ratio test based on the spatial distribution and frequency of derived polymorphic variants corroborates this result and provides an estimate of the strength of selection. However, the data are difficult to reconcile with the simplest model of positive selection, whereas a new composite-likelihood method suggests that the data are better described by a model in which the selected allele has not yet gone to fixation.     

Genetic mapping of adaptive wing size variation in Drosophila simulans

Many ecologically important traits exhibit latitudinal variation. Body size clines have been described repeatedly in insects across multiple continents, suggesting that similar selective forces are shaping these geographical gradients. It is unknown whether these parallel clinal patterns are controlled by the same or different genetic mechanism(s). We present here, quantitative trait loci (QTL) analysis of wing size variation in Drosophila simulans. Our results show that much of the wing size variation is controlled by a QTL on Chr 3L with relatively minor contribution from other chromosome arms. Comparative analysis of the genomic positions of the QTL indicates that the major QTL on Chr 3 are distinct in D. simulans and D. melanogaster, whereas the QTL on Chr 2R might overlap between species. Our results suggest that parallel evolution of wing size clines could be driven by non-identical genetic mechanisms but in both cases involve a major QTL as well as smaller effects of other genomic regions.