Synonymous Substitutions in the Xdh Gene of Drosophila: Heterogeneous Distribution along the Coding Region

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

Differentiation of Muller''s Chromosomal Elements D and E in the Obscura Group of Drosophila

Twenty-two markers located on Muller's elements D or E have been mapped by in situ hybridization in six species of the obscura group of Drosophila and in D. melanogaster. The obscura species can be grouped into a Palearctic cluster (D. subobscura, D. madeirensis and D. guanche) and a Nearctic one (D. pseudoobscura, D. persimilis and D. miranda). Eleven of the probes contain known genes: E74, Acp70A, Est5, hsp28/23, hsp83, emc, hsp70, Xdh, Acph-1, Cec and rp49. The remaining probes are recombinant phages isolated from a D. subobscura genomic library. All these markers hybridize to the putative homologous chromosome or chromosomal arm of elements D and E. Thus, these elements have conserved their genic content during species divergence. Chromosomal homologies proposed previously for each element among the species of the same cluster have been compared with the present results. The distribution of markers within each element has changed considerably as inferred from pairwise comparisons of obscura species included in the two different clusters. Only chromosomal segments defined by closely linked markers have been conserved: one such segment has been detected in element D and three in element E between D. subobscura and D. pseudoobscura.     

Estimates of Linkage Disequilibrium and the Recombination Parameter Determined from Segregating Nucleotide Sites in the Alcohol Dehydrogenase Region of Drosophila Pseudoobscura

The alcohol dehydrogenase (Adh) region of Drosophila pseudoobscura, which includes the two genes Adh and Adh-Dup, was used to examine the pattern and organization of linkage disequilibrium among pairs of segregating nucleotide sites. A collection of 99 strains from the geographic range of D. pseudoobscura were nucleotide-sequenced with polymerase chain reaction-mediated techniques. All pairs of the 359 polymorphic sites in the 3.5-kb Adh region were tested for significant linkage disequilibrium with Fisher's exact test. Of the 74,278 pairwise comparisons of segregating sites, 127 were in significant linkage disequilibrium at the 5% level. The distribution of five linkage disequilibrium estimators D(ij), D(2), r(ij), r(2) and Dij were compared to theoretical distributions. The observed distributions of D(ij), D(2), r(ij) and r(2) were consistent with the theoretical distribution given an infinite sites model. The observed distribution of Dij differed from the theoretical distribution because of an excess of values at -1 and 1. No spatial pattern was observed in the linkage disequilibrium pattern in the Adh region except for two clusters of sites nonrandomly associated in the adult intron and intron 2 of Adh. The magnitude of linkage disequilibrium decreases significantly as nucleotide distance increases, or a distance effect. Adh-Dup had a larger estimate of the recombination parameter, 4Nc, than Adh, where N is the effective population size and c is the recombination rate. A comparison of the mutation and recombination parameters shows that 7-17 recombination events occur for each mutation event. The heterogeneous estimates of the recombination parameter and the inverse relationship between linkage disequilibrium and nucleotide distance are no longer significant when the two clusters of Adh intron sites are excluded from analyses. The most likely explanation for the two clusters of linkage disequilibria is epistatic selection between sites in the cluster to maintain pre-mRNA secondary structure.     

Switch in codon bias and increased rates of amino acid substitution in the Drosophila saltans species group.

We investigated the nucleotide composition of five genes, Xdh, Adh, Sod, Per, and 28SrRNA, in nine species of Drosophila (subgenus Sophophora) and one of Scaptodrosophila. The six species of the Drosophila saltans group markedly differ from the others in GC content and codon use bias. The GC content in the third codon position, and to a lesser extent in the first position and the introns, is higher in the D. melanogaster and D. obscura groups than in the D. saltans group (in Scaptodrosophila it is intermediate but closer to the melanogaster and obscura species). Differences are greater for Xdh than for Adh, Sod, Per, and 28SrRNA, which are functionally more constrained. We infer that rapid evolution of GC content in the saltans lineage is largely due to a shift in mutation pressure, which may have been associated with diminished natural selection due to smaller effective population numbers rather than reduced recombination rates. The rate of GC content evolution impacts the rate of protein evolution and may distort phylogenetic inferences. Previous observations suggesting that GC content evolution is very limited in Drosophila may have been distorted due to the restricted number of genes and species (mostly D. melanogaster) investigated.     

Evolutionary rearrangement of the amylase genomic regions between Drosophila melanogaster and Drosophila pseudoobscura

Two Drosophila pseudoobscura genomic clones have sequence similarity to the Drosophila melanogaster amylase region that maps to the 53CD region on the D. melanogaster cytogenetic map. The two clones with similarity to amylase map to sections 73A and 78C of the D. pseudoobscura third chromosome cytogenetic map. The complete sequences of both the 73A and 78C regions were compared to the D. melanogaster genome to determine if the coding region for amylase is present in both regions and to determine the evolutionary mechanism responsible for the observed distribution of the amylase gene or genes. The D. pseudoobscura 73A and 78C linkage groups are conserved with the D. melanogaster 41E and 53CD regions, respectively. The amylase gene, however, has not maintained its conserved linkage between the two species. These data indicate that amylase has moved via a transposition event in the D. melanogaster or D. pseudoobscura lineage. The predicted genes within the 73A and 78C regions show patterns of molecular evolution in synonymous and nonsynonymous sites that are consistent with previous studies of these two species.     

Thermal adaptive significance of ADH and EST-6 allozymes in Indian geographical populations of Drosophila melanogaster

Indian geographical populations of Drosophila melanogaster exhibit significant correlation (r 0.95) of allelic frequencies at Est -6 and Adh loci with latitude as well as altitude. Est -6^S and AdhF allozymes are well adapted to colder environments while Est -6^F and AdhS are warm adapted. The data on allozymic clines match climatic conditions on the Indian subcontinent. On the basis of multiple regression analysis, one major conclusion is that coefficient of variation of temperature ( T _CV) along latitude/altitude accounts for alterations in allelic frequency at the Adh locus while T _max and T _max explain changes at the Est -6 locus. Thus, climatic conditions lead to thermal selection of allozymes in Indian populations of D. melanogaster .     

Molecular Population Genetics of the Distal Portion of the X Chromosome in Drosophila: Evidence for Genetic Hitchhiking of the Yellow-Achaete Region

We have estimated DNA sequence variation and differentiation within and between Drosophila melanogaster and its sibling species, Drosophila simulans, using six-cutter restriction site variation at yellow-achaete (y-ac), phosphogluconate dehydrogenase (Pgd), and period (per). These three gene regions are of varying distance from the telomere of the X chromosome and range from very low to moderate rates of recombination in D. melanogaster. According to Tajima's test of neutrality, the Pgd region has been influenced by balancing selection in D. melanogaster. This is consistent with previous data suggesting the allozyme polymorphism at this locus is visible to selection. The Hudson, Kreitman, Aguadé test of neutrality reveals a significant departure from neutrality for the y-ac region compared to the per or rosy regions in D. simulans. There is also a significant departure for the y-ac region compared to the Adh 5' flanking region in D. melanogaster. In both species the departure appears to be due to reduced variation at y-ac compared to that expected from divergence between D. simulans and D. melanogaster. We conclude that recent hitchhiking associated with the selective fixation of one or more advantageous mutants in the y-ac region is the best explanation for reduced variation at y-ac.     

The period gene of Drosophila carries species-specific behavioral instructions.

We have analyzed and compared the circadian locomotor activity rhythms of Drosophila melanogaster and D.pseudoobscura. The rhythms of D.pseudoobscura are stronger and the periods shorter than those of D.melanogaster. We have also transformed D.melanogaster flies with a hybrid gene containing the coding region of the D.pseudoobscura period (per) gene. Behavioral assays of flies containing this hybrid gene show that the per protein encoded by the D.pseudoobscura per gene is able to rescue the rhythmic deficiencies of arrhythmic, pero1 D.melanogaster. More important, the rhythms of some of these strains are stronger and the periods shorter than those of D.melanogaster (and those of transformants which carry the equivalent D.melanogaster per gene construct) and hence resemble those of D.pseudoobscura. The results suggest that the primary amino acid sequence of the per gene encodes species-specific behavioral instructions that are detectable when only the per gene is transferred to a different species.     

Estimates of Gene Flow in Drosophila Pseudoobscura Determined from Nucleotide Sequence Analysis of the Alcohol Dehydrogenase Region

The genetic structure of Drosophila pseudoobscura populations was inferred from a nucleotide sequence analysis of a 3.4-kb segment of the alcohol dehydrogenase (Adh) region. A total of 99 isochromosomal strains collected from 13 populations in North and South America were used to determine if any population departed from a neutral model and to estimate levels of gene flow between populations. This study also included the nucleotide sequences from two sibling species, D. persimilis and D. miranda. We estimated the neutral mutation parameter, 4N mu, in synonymous and noncoding sites for 17 subregions of Adh in each of nine populations with sample sizes greater than three. The nucleotide diversity data in the nine populations was tested for departures from an equilibrium neutral model with two statistical tests. The Tajima and the Hudson, Kreitman, Aguade tests showed that each population fails to reject a neutral model. Tests for genetic differentiation between populations fail to show any population substructure among the North American populations of D. pseudoobscura. The nucleotide diversity data is consistent with direct and indirect measures of gene flow that show extensive dispersal between populations of D. pseudoobscura.     

The molecular evolution of the alcohol dehydrogenase and alcohol dehydrogenase-related genes in the Drosophila melanogaster species subgroup

The DNA sequences of the Adh genes of three members of the Drosophila melanogaster species subgroup have been determined. This completes the Adh sequences of the eight species of this subgroup. Two species, D. yakuba and D. teissieri, possess processed Adh pseudogenes. In all of the species of the subgroup, a gene of unknown function, Adhr, is located about 300 bp 3' to Adh. Although this gene is experiencing a higher rate of synonymous substitution than Adh, it is more constrained at the amino acid level. Phylogenetic relationships between all eight members of the melanogaster subgroup have been analyzed using a variety of methods. All analyses suggested that the D. yakuba and D. teissieri pseudogenes have a single common ancestor, rather than evolving independently in each species, and that D. melanogaster is the sister species to D. simulans, D. sechellia, and D. mauritiana. The evolutionary relationships of the latter three species remain equivocal.      

The Drosophila subobscura Adh genomic region contains valuable evolutionary markers.

We have sequenced 4 kb of the genomic region comprising the Adh (Alcohol dehydrogenase) gene of Drosophila subobscura. In agreement with other species which belong to the same subgenus, two structural genes, Adh and Adh-dup, are contained in this region. The main features of these two genes of D. subobscura have been inferred from the sequence data and compared with the homologous region of D. ambigua and D. pseudoobscura. Drosophila subobscura Adh and Adh-dup differ from those of D. ambigua at a corrected estimation of 10.1% and 12.5%, respectively, while from those of D. pseudoobscura they differ by 9.5% and 8.1%, respectively. Our data suggest that Adh and Adh-dup are evolving independently, showing a species-specific pattern. Moreover, particular features of some regions of these genes make them valuable evolutionary hallmarks. For instance, replacement substitutions in the third exon of Adh may indicate the branching of the melanogaster-obscura groups, whereas replacement substitutions in the third exon of the Adh-dup could be used to assess speciation within the obscura group.     

Organization and evolution of the alcohol dehydrogenase gene in Drosophila

The alcohol dehydrogenase (Adh) gene was isolated from Drosophila simulans and D. mauritiana, and the DNA sequence of a 4.6-kb region, containing the structural gene and flanking sequence, was determined for each. These sequences were compared with the Adh region of D. melanogaster to characterize changes that occur in the Drosophila genome during evolution and to identify conserved sequences of functional importance. Drosophila simulans and D. mauritiana Adh are organized in a manner similar to that of D. melanogaster Adh, including the presence of two promoters for the single Adh gene. This study identified conserved flanking elements that, in conjunction with other studies, suggest regions that may be involved in the control of Adh expression. Inter- and intraspecies comparisons revealed differences in the kinds of sequence changes that have accumulated. Sequence divergence in and around the Adh gene was used to assess inter- and intraspecies evolutionary relationships. Finally, there appears to be an unrelated structural gene located directly 3' of the Adh transcribed region.      

Gene Flow and Natural Selection in the Origin of Drosophila Pseudoobscura and Close Relatives

The divergence of Drosophila pseudoobscura and close relatives D. persimilis and D. pseudoobscura bogotana has been studied using comparative DNA sequence data from multiple nuclear loci. New data from the Hsp82 and Adh regions, in conjunction with existing data from Adh and the Period locus, are examined in the light of various models of speciation. The principal finding is that the three loci present very different histories, with Adh indicating large amounts of recent gene flow among the taxa, while little or no gene flow is apparent in the data from the other loci. The data were compared with predictions from several isolation models of divergence. These models include no gene flow, and they were found to be incompatible with the data. Instead the DNA data, taken together with other evidence, seem consistent with divergence models in which natural selection acts against gene flow at some loci more than at others. This family of models includes some sympatric and parapatric speciation models, as well as models of secondary contact and subsequent reinforcement of sexual isolation.     

The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.     

Inferring the history of speciation from multilocus DNA sequence data: the case of Drosophila pseudoobscura and close relatives

The divergence of Drosophila pseudoobscura from its close relatives, D. persimilis and D. pseudoobscura bogotana, was examined using the pattern of DNA sequence variation in a common set of 50 inbred lines at 11 loci from diverse locations in the genome. Drosophila pseudoobscura and D. persimilis show a marked excess of low-frequency variation across loci, consistent with a model of recent population expansion in both species. The different loci vary considerably, both in polymorphism levels and in the levels of polymorphisms that are shared by different species pairs. A major question we address is whether these patterns of shared variation are best explained by gene flow or by persistence since common ancestry. A new test of gene flow, based on patterns of linkage disequilibrium, is developed. The results from these, and other tests, support a model in which D. pseudoobscura and D. persimilis have exchanged genes at some loci. However, the pattern of variation suggests that most gene flow, although occurring after speciation began, was not recent. There is less evidence of gene flow between D. pseudoobscura and D. p. bogotana. The results are compared with recent work on the genomic locations of genes that contribute to reproductive isolation between D. pseudoobscura and D. persimilis. We show that there is a good correspondence between the genomic regions associated with reproductive isolation and the regions that show little or no evidence of gene flow.