The Rosy Region of Drosophila Melanogaster and Drosophila Simulans. I. Contrasting Levels of Naturally Occurring DNA Restriction Map Variation and Divergence

A 40-kb region around the rosy and snake loci was analyzed for restriction map variation among 60 lines of Drosophila melanogaster and 30 lines of Drosophila simulans collected together at a single locality in Raleigh, North Carolina. DNA sequence variation in D. simulans was estimated to be 6.3 times greater than in D. melanogaster (heterozygosities per nucleotide of 1.9% vs. 0.3%). This result stands in marked contrast to results of studies of phenotypic variation including proteins (allozymes), morphology and chromosome arrangements which are generally less variable and less geographically differentiated in D. simulans. Intraspecific polymorphism is not distributed uniformly over the 40-kb region. The level of heterozygosity per nucleotide varies more than 12-fold across the region in D. simulans, being highest over the hsc2 gene. Similar, though less extreme, variation in heterozygosity is also observed in D. melanogaster. Average interspecific divergence (corrected for intraspecific polymorphism) averaged 3.8%. The pattern of interspecific divergence over the 40-kb region shows some disparities with the spatial distribution of intraspecific variation, but is generally consistent with selective neutrality predictions: the most polymorphic regions within species are generally the most divergent between species. Sequence-length polymorphism is observed for D. melanogaster to be at levels comparable to other gene regions in this species. In contrast, no sequence length variation was observed among D. simulans chromosomes (limit of resolution approximately 100 bp). These data indicate that transposable elements play at best a minor role in the generation of naturally occurring genetic variation in D. simulans compared to D. melanogaster. We hypothesize that differences in species effective population size are the major determinant of the contrasting levels and patterns of DNA sequence and insertion/deletion variation that we report here and the patterns of allozyme and morphological variation and differentiation reported by other workers for these two species.     

Molecular and Phenotypic Variation of the Zw Locus Region in Drosophila Melanogaster

Restriction map polymorphism in a 13-kb region of the Zw locus in Drosophila melanogaster was investigated for 64 X chromosome lines with seven 6-cutter and ten 4-cutter restriction enzymes. A total of 203 restriction sites were scored, of which 20 were found to be polymorphic. The estimated nucleotide variation for this region for overall data (pi = 0.003 and 0.001, and theta = 0.003 and 0.002, for 4-cutter and 6-cutter studies, respectively) was smaller than that reported for most regions studied in D. melanogaster. It was found that the Slow allozyme has a larger nucleotide variation and haplotype diversity than the Fast allozyme. Results suggest the relatively recent divergence of the Fast allozyme from the Slow allozyme. Glucose 6-phosphate dehydrogenase (G6PD) activity was measured as a phenotype of the Zw locus. A significant difference in G6PD activity between allozymes was detected. The between-line effect was highly significant within the Slow allozyme, but was not significant within the Fast allozyme. Although a direct causative link could not be established, these results suggest an association between the amounts of quantitative and molecular genetic variation at the Zw locus region.     

Genetic Variation in Natural Populations of Five Drosophila Species and the Hypothesis of the Selective Neutrality of Protein Polymorphisms

We have studied genetic variation at 30-32 loci coding for enzymes in natural populations of five species of Drosophila. The average proportion of heterozygous loci per individual is 17.7 +/- 0.4%. The average proportion of polymorphic loci per population is 69.2 +/- 2.6% or 49.8 +/- 2.2%, depending on what criterion of polymorphism is used. The following generalizations are advanced: (1) The amount of genetic polymorphism varies considerably from locus to locus. (2) At a given locus, populations of the same species are very similar in the amount and pattern of genetic variation. (3) However, at some loci large differences sometimes occur between local populations of the same species. (4) The amount of variation at a given locus is approximately the same in all five species. (5) When different species are compared, the pattern of the variation is either essentially identical or totally different at a majority of loci. We have tested the hypothesis that protein polymorphisms are selectively neutral by examining four predictions derived from the hypothesis. Our results are at variance with every one of the predictions. We have measured the amount of genetic differentiation, D, between taxa of various degrees of evolutionary divergence. The average value of D is 0.033 for local populations, 0.228 for subspecies, 0.226 for semispecies, 0.538 for sibling species, and 1.214 for morphologically distinguishable species. Our results indicate that a substantial degree of genetic differentiation (22.8 allelic substitutions for every 100 loci) occurs between allopatric populations that have diverged to the point where they might become different species if they were to become sympatric. However, very little additional genetic change is required for the development of complete reproductive isolation. After the speciation process is completed, species continue to diverge genetically from each other.     

Nucleotide Variation and Divergence in the Histone Multigene Family in Drosophila Melanogaster

Nucleotide differences in the histone H3 gene family in Drosophila melanogaster were studied on three levels: (1) within a chromosome, (2) within a population and (3) between species (D. melanogaster and Drosophila simulans). The average difference within the H3 gene within a chromosome was 0.0040 per nucleotide site, about 52% of that within a population (0.0077). The proportion of divergent sites between the two species was 0.0575, which is about 8.5 times the difference within a species. The distribution of divergence between species was similar to that of variation within a species. Divergence and variation were noted to be greatest in the 3' noncoding region and least in the coding region. Values intermediate between these were found for the 5' noncoding region. Divergence and variation in silent sites exceeded those in the total coding region, thus indicating possible purifying selection for amino-acid-altering change. Phylogenetic relations among H3 genes and genetic differences on these three levels are evidence for the concerted evolution of the histone gene family. The molecular mechanism by which variation is produced and maintained is discussed.     

Molecular Population Genetics of Sex Determination Genes: The Transformer Gene of Drosophila Melanogaster

The transformer locus (tra) produces an RNA processing protein that alternatively splices the doublesex pre-mRNA in the sex determination hierarchy of Drosophila melanogaster. Comparisons of the tra coding region among Drosophila species have revealed an unusually high degree of divergence in synonymous and nonsynonymous sites. In this study, we tested the hypothesis that the tra gene will be polymorphic in synonymous and nonsynonymous sites within species by investigating nucleotide sequence variation in eleven tra alleles within D. melanogaster. Of the 1063 nucleotides examined, two synonymous sites were polymorphic and no amino acid variation was detected. Three statistical tests were used to detect departures from an equilibrium neutral model. Two tests failed to reject a neutral model of molecular evolution because of low statisitical power associated with low levels of genetic variation (Tajima/Fu and Li). The Hudson, Kreitman, and Aguade test rejected a neutral model when the tra region was compared to the 5'-flanking region of alcohol dehydrogenase (Adh). The lack of variability in the tra gene is consistent with a recent selective sweep of a beneficial allele in or near the tra locus.     

A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster. III. Variations in Genetic Structure and Their Causes between Drosophila melanogaster and Its Sibling Species Drosophila simulans

The natural populations of Drosophila melanogaster and Drosophila simulans were compared for their genetic structure. A total of 114 gene-protein loci were studied in four mainland (from Europe and Africa) and an island (Seychelle) populations of D. simulans and the results were compared with those obtained on the same set of homologous loci in fifteen worldwide populations of D. melanogaster. The main results are as follows: (1) D. melanogaster shows a significantly higher proportion of loci polymorphic than D. simulans (52% vs. 39%, P<0.05), (2) both species have similar mean heterozygosity and mean number of alleles per locus, (3) the two species share some highly polymorphic loci but they do not share loci that show high geographic differentiation, and (4) D. simulans shows significantly less geographic differentiation than D. melanogaster. The differences in genetic differentiation between the two species are limited to loci located on the X and second chromosomes only; loci on the third chromosome show similar level of geographic differentiation in both species. These two species have previously been shown to differ in their pattern of variation for chromosomal polymorphisms, quantitative and physiological characters, two-dimensional electrophoretic (2DE) proteins, middle repetitive DNA and mitochondrial DNA. Variation in niche-widths and/or genetic "strategies" of adaptation appear to be the main causes of differences in the genetic structure of these two species.     

Nucleotide Variation at the Gpdh Locus in the Genus Drosophila

The Gpdh locus was sequenced in a broad range of Drosophila species. In contrast to the extreme evolutionary constraint seen at the amino acid level, the synonymous sites evolve at rates comparable to those of other genes. Gpdh nucleotide sequences were used to infer a phylogenetic tree, and the relationships among the species of the obscura group were examined in detail. A survey of nucleotide polymorphism within D. pseudoobscura revealed no amino acid variation in this species. Applying a modified McDonald-Kreitman test, the amino acid divergence between species in the obscura group does not appear to be excessive, implying that drift is adequate to explain the patterns of amino acid change at this locus. In addition, the level of polymorphism at the Gpdh locus in D. pseudoobscura is comparable to that found at other loci, as determined by a Hudson-Kreitman-Aguade test. Thus, the pattern of nucleotide variation within and between species at the Gpdh locus is consistent with a neutral model.     

Unusual pattern of single nucleotide polymorphism at the exuperantia2 locus of Drosophila pseudoobscura

We have investigated the pattern of DNA sequence variation at the exuperantia2 locus in Drosophila pseudoobscura. This adds to the increasing dataset of genetic variation in D. pseudoobscura, a useful model species for evolutionary genetic studies. The level of silent site nucleotide diversity and the divergence from an outgroup Drosophila miranda are comparable with those for other X-linked loci. One peculiar pattern at the exu2 locus of D. pseudoobscura is a complete linkage disequilibrium between two SNPs, one of which is a replacement site. As a result, there are two distinct haplotype groups in our dataset. Based upon the comparisons with the outgroup sequences from D. miranda and Drosophila persimilis, we show that the newly derived haplotype group has lower diversity than the ancestral haplotype group. The pattern of protein evolution at exu2 shows some deviation from the neutral model. Together, these and other characteristics of the exu2 locus suggest the action of selection on the pattern of SNP variation, consistent with a partial selective sweep associated with the newly derived haplotype.     

Contrasting modes of natural selection acting on pigmentation genes in the Drosophila dunni subgroup

Genes that encode for divergent adaptive traits may have genealogies that contrast with those from loci that are not functionally involved in differentiation. Here, we examine DNA sequence variation among the species of the eastern Caribbean Drosophila dunni subgroup at two loci, yellow and dopa decaboxylase (Ddc), which both play integral roles in pigmentation patterning of adult Drosophila. Phylogenetic analyses of these loci produce gene genealogies with topologies that mirror those described for other nuclear genes: the six morphologically distinct species within the subgroup are divided into only three lineages, with one lineage containing four species that share extensive ancestral polymorphism. At the Ddc locus these major lineages are delineated only by silent site variation. We observe a significantly higher rate of synonymous site divergence than non-synonymous divergence, consistent with strong purifying selection acting on the locus. In contrast, the yellow locus exhibits patterns of amino acid divergence and nucleotide diversity that are consistent with recent diversifying selection acting in two different lineages. This selection appears to be targeting amino acid variants in the signal sequence of the Yellow protein, a region which is tightly constrained among members of the larger D. cardini radiation. This result highlights not only the potential importance of yellow in the evolution of divergent pigmentation patterns among members of the D. dunni subgroup, but also hints that variation in signal peptide sequences may play a role in phenotypic diversification.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.     

Structure of Genetic Variation within and between Populations of Mycophagous Drosophila

Patterns of genetic variation within and between populations of five species of mycophagous Drosophila were examined by gel electrophoresis of several polymorphic loci. Populations of the five species could not be shown to be subdivided into sympatric host-adapted races. Statistically significant, but small, between-host differences in gene frequencies were observed at three of 15 loci. Mean gene frequencies at all loci were similar in New York and Tennessee, and, with one exception, relatively little genetic differentiation was observed among study sites within those two regions. Gene frequencies generally were stable over several years of collecting as well. The unpredictable nature of the fungal hosts may preclude the site fidelity and continuity of diversifying selection necessary for adaptive divergence of populations.     

Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

Molecular Evolution between Drosophila Melanogaster and D. Simulans: Reduced Codon Bias, Faster Rates of Amino Acid Substitution, and Larger Proteins in D. Melanogaster

Both natural selection and mutational biases contribute to variation in codon usage bias within Drosophila species. This study addresses the cause of codon bias differences between the sibling species, Drosophila melanogaster and D. simulans. Under a model of mutation-selection-drift, variation in mutational processes between species predicts greater base composition differences in neutrally evolving regions than in highly biased genes. Variation in selection intensity, however, predicts larger base composition differences in highly biased loci. Greater differences in the G+C content of 34 coding regions than 46 intron sequences between D. melanogaster and D. simulans suggest that D. melanogaster has undergone a reduction in selection intensity for codon bias. Computer simulations suggest at least a fivefold reduction in N(e)s at silent sites in this lineage. Other classes of molecular change show lineage effects between these species. Rates of amino acid substitution are higher in the D. melanogaster lineage than in D. simulans in 14 genes for which outgroup sequences are available. Surprisingly, protein sizes are larger in D. melanogaster than in D. simulans in the 34 genes compared between the two species. A substantial fraction of silent, replacement, and insertion/deletion mutations in coding regions may be weakly selected in Drosophila.     

Lack of Genic Similarity between Two Sibling Species of Drosophila as Revealed by Varied Techniques

Acrylamide gel electrophoresis was performed on the enzyme xanthine dehydrogenase in sixty isochromosomal lines of Drosophila persimilis from three geographic populations. Sequential electrophoretic analysis using varied gel concentrations and buffers revealed twenty-three alleles in this species where only five had been described previously. These new electrophoretic techniques also detected a profound increase in divergence of gene frequencies at this locus between D. persimilis and its sibling species D. pseudoobscura. The implications of these results for questions of speciation and the maintenance of genetic variability are discussed.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

Electrophoretic variability in natural populations of Drosophila melanogaster and Drosophila simulans

Genetic variation at 59 gene loci coding for enzymes (50) and larval proteins (9) has been studied in sympatric populations of Drosophila melanogaster and D. simulans from insular and continental origin. The average number of alleles per locus, the mean proportion of polymorphic loci and the mean heterozygosity are similar both within and between species. There are however some significant differences between D. simulans populations in the genotypic frequencies for four polymorphic loci.     

Interaction of the Stubble-Stubbloid Locus and the Broad-Complex of Drosophila Melanogaster

The 2B5 region on the X chromosome of Drosophila melanogaster forms an early ecdysone puff at the end of the third instar. The region is coextensive with a complex genetic locus, the Broad-Complex (BR-C). The BR-C is a regulatory gene that contains two major functional domains, the br domain and the l(1)2Bc domain. BR-C mutants prevent metamorphosis, including morphogenesis of imaginal discs; br mutants prevent elongation and eversion of appendages and l(1)2Bc mutants prevent fusion of the discs. The Stubble-stubbloid (Sb-sbd) locus at 89B9-10 is best known for the effects of its mutants on bristle structure. Mutants of the BR-C and the Sb-sbd locus interact to produce severe malformation of appendages. Viable heteroallelic and homoallelic combinations of Sb-sbd mutants, including loss-of-function mutants, affect the elongation of imaginal disc appendages. Thus, the Sb-sbd(+) product is essential for normal appendage elongation. Sb-sbd mutants, however, do not affect eversion or fusion of discs. Correspondingly, only BR-C mutants deficient in br function interact with Sb-sbd mutants. The interaction occurs in deficiency heterozygotes using single, wild-type doses of the BR-C, of the Sb-sbd locus, or of both loci. These last results are formally consistent with the possibility that the BR-C acts as a positive regulator of the Sb-sbd locus. The data do not exclude other possible nonregulatory interactions between the two loci, e.g., interactions between the products of both genes.     

Clines and adaptive evolution in the methuselah gene region in Drosophila melanogaster

In an effort to characterize further the patterns of selection and adaptive evolution at the methuselah locus in Drosophila species, we extended an analysis of geographical variation to include single nucleotide polymorphisms (SNPs) in adjacent genes on either side of the mth locus, and examined the molecular variation in a neighbouring methuselah paralogue (mth2). An analysis of 13 SNPs spanning a region of nearly 19 kilobases surrounding the mth locus demonstrated that a clinal pattern associated with the most common mth haplotype does not extend to adjacent gene loci, providing compelling evidence that the clinal pattern results from selection on as yet unidentified sites associated with the functional mth locus. mth2 exhibited a significant pattern of adaptive divergence among D. melanogaster, D. simulans and D. yakuba similar to that seen at mth. However, Ka : Ks ratios indicate a difference in levels of functional constraint at the two methuselah, loci with mth2 exhibiting a five- to six-fold reduction in levels of amino acid divergence relative to mth.