Latitudinal clines for nucleotide polymorphisms in the Esterase 6 gene of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Previous studies have found non-neutral patterns of nucleotide polymorphism in the promoter and coding regions of Est6 in D. melanogaster. Coding region polymorphism peaks around two closely linked replacement differences associated with the EST6-F/EST6-S allozyme polymorphism. The promoter contains two common, highly diverged haplotype groups, P1 and P7, that differentially affect Est6 expression. Allozyme studies have also revealed latitudinal clines in EST6-F and EST6-S frequencies that recur across continents. Here we analyse nucleotide polymorphisms across the promoter and the region of peak coding sequence polymorphism in 10 Australian populations along a 25° latitudinal gradient in order to examine the basis for the allozyme clines. As with the earlier studies, we find an excess of intermediate to high frequency variants in both the P1/P7 region and around the two EST6-F/EST6-S replacements in some populations. The two EST6-F/EST6-S replacement polymorphisms show latitudinal clines whereas the P1 and P7 groups of promoter haplotypes do not. However the strongest clines are for three co-segregating silent site polymorphisms in a 4 bp stretch at the 3′ end of the sequenced region. Monte Carlo simulations show that the clines for those three sites can explain all others in the data but none of the others can explain those three. Thus the allozyme clines may not reflect selection on either the P1/P7 polymorphism or the two replacements previously associated with the EST6-F/EST-S difference.     

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.     

Historical Selection, Amino Acid Polymorphism and Lineage-Specific Divergence at the G6pd Locus in Drosophila Melanogaster and D. Simulans

The nucleotide diversity across 1705 bp of the G6pd gene is studied in 50 Drosophila melanogaster and 12 D. simulans lines. Our earlier report contrasted intraspecific polymorphism and interspecific differences at silent and replacement sites in these species. This report expands the number of European and African lines and examines the pattern of polymorphism with respect to the common A/B allozymes. In D. melanogaster the silent nucleotide diversity varies 2.8-fold across localities. The B allele sequences are two- to fourfold more variable than the derived A allele, and differences between allozymes are twice as among B alleles. There is strong linkage disequilibrium across the G6pd region. In both species the level of silent polymorphism increases from the 5' to 3' ends, while there is no comparable pattern in level of silent site divergence or fixation. The neutral model is not rejected in either species. Using D. yakuba as an outgroup, the D. melanogaster lineage shows a twofold greater rate of silent fixation, but less than half the rate of amino acid replacement. Lineage-specific differences in mutation fixation are inconsistent with neutral expectations and suggest the interaction of species-specific population size differences with both weakly advantageous and deleterious selection.     

Studies of esterase 6 in Drosophila melanogaster. XVIII. Biochemical differences between the slow and fast allozymes

Most natural populations of Drosophila melanogaster are polymorphic for two major electrophoretic variants at the esterase-6 locus. The frequency of the EST 6F allozyme is greatest in populations in warmer latitudes, whereas the EST 6S allozyme is predominant in colder latitudes. Latitudinal clines in electromorph frequencies are found on three continents. Purified preparations of the allozymes have been characterized for their pH optimum, substrate specificity, organophosphate inhibition, alcohol activation, thermal stability, and kinetic parameters. These and previous analyses of the EST 6 allozymes reveal that the two variants have differences in their physical and kinetic properties that may provide a basis for the selective maintenance of the polymorphisms and an explanation of the clinal variation observed in natural populations.      

Nucleotide Polymorphism in the 5' Promoter Region of Esterase 6 in Drosophila Melanogaster and Its Relationship to Enzyme Activity Variation

A 974-bp region immediately 5' of the esterase 6 gene was sequenced in 17 field derived third chromosome isoallelic lines. Twenty-three polymorphisms were identified, only two in the first 400 bp 5' but 16 in a 325-bp region from -494 to -819 bp. This distribution differs from previously published patterns in Drosophila simulans and D. mauritiana, where the first 800 bp are highly conserved. Fourteen common polymorphisms in the 325-bp region above are all in strong linkage disequilibrium with each other. Moreover, most of the haplotypes defined by the total of 23 polymorphisms fall into two groups that differ as a block at all 14 of these latter sites. Sequence differences between the two groups include some restriction sites that were scored in an earlier study of RFLPs and EST6 enzyme phenotypes among 42 isoallelic lines from the same population. By collating the two studies, we show that one haplotype group yields ~15% lower EST6 enzyme activity in adult males than the other. The promoter haplotypes show only weak disequilibrium with the esterase 6 fast/slow allozyme polymorphism, so it seems unlikely that previously reported latitudinal clines in the allozyme frequencies are due to their hitchhiking along with selection on the promoter difference.     

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Intraspecific nuclear DNA variation in Drosophila

We have summarized and analyzed all available nuclear DNA sequence polymorphism studies for three species of Drosophila, D. melanogaster (24 loci), D. simulans (12 loci), and D. pseudoobscura (5 loci). Our major findings are: (1) The average nucleotide heterozygosity ranges from about 0.4% to 2% depending upon species and function of the region, i.e., coding or noncoding. (2) Compared to D. simulans and D. pseudoobscura (which are about equally variable), D. melanogaster displays a low degree of DNA polymorphism. (3) Noncoding introns and 3' and 5' flanking DNA shows less polymorphism than silent sites within coding DNA. (4) X-linked genes are less variable than autosomal genes. (5) Transition (Ts) and transversion (Tv) polymorphisms are about equally frequent in non-coding DNA and at fourfold degenerate sites in coding DNA while Ts polymorphisms outnumber Tv polymorphisms by about 2:1 in total coding DNA. The increased Ts polymorphism in coding regions is likely due to the structure of the genetic code: silent changes are more often Ts's than are replacement substitutions. (6) The proportion of replacement polymorphisms is significantly higher in D. melanogaster than in D. simulans. (7) The level of variation in coding DNA and the adjacent noncoding DNA is significantly correlated indicating regional effects, most notably recombination. (8) Surprisingly, the level of polymorphism at silent coding sites in D. melanogaster is positively correlated with degree of codon usage bias. (9) Three proposed tests of the neutral theory of DNA polymorphisms have been performed on the data: Tajima's test, the HKA test, and the McDonald-Kreitman test. About half of the loci fail to conform to the expectations of neutral theory by one of the tests. We conclude that many variables are affecting levels of DNA polymorphism in Drosophila, from properties of nucleotides to population history and, perhaps, mating structure. No simple, all encompassing explanation satisfactorily accounts for the data.      

The functional impact of Pgm amino acid polymorphism on glycogen content in Drosophila melanogaster

Earlier studies of the common PGM allozymes in Drosophila melanogaster reported no in vitro activity differences. However, our study of nucleotide variation observed that PGM allozymes are a heterogeneous mixture of amino acid polymorphisms. In this study, we analyze 10 PGM protein haplotypes with respect to PGM activity, thermostability, and adult glycogen content. We find a twofold difference in activity among PGM protein haplotypes that is associated with a threefold difference in glycogen content. The latitudinal clines for several Pgm amino acid polymorphisms show that high PGM activity, and apparently higher flux to glycogen synthesis, parallel the low activity clines at G6PD for reduced pentose shunt flux in northern latitudes. This suggests that amino acid polymorphism is under selection at this branch point and may be favored for increased metabolic storage associated with stress resistance and adaptation to temperate regions.     

Contrasting molecular population genetics of four hexokinases in Drosophila melanogaster, D. simulans and D. yakuba

As part of a larger study contrasting patterns of variation in regulatory and nonregulatory enzymes of the central metabolic pathways we have examined the molecular variation in four uncharacterized hexokinase genes unique to muscle, fat body, and testis in Drosophila melanogaster, D. simulans, and D. yakuba. Earlier isoenzyme studies had designated these genes as Hex-A, Hex-C, and Hex-t. There are two tightly linked testes-specific genes designated here as Hex-t1 and Hex-t2. Substantial and concordant differences across species are seen in levels of both amino acid and silent polymorphism. The flight muscle form Hex-A is the most conserved followed by the fat body hexokinase Hex-C and testis-specific hexokinases Hex-t1 and Hex-t2. While constraints acting at the amino acid level are expected, the silent polymorphisms follow this pattern as well. All genes are in regions of normal recombination, therefore hitchhiking and background selection are not likely causes of interlocus differences. In D. melanogaster latitudinal clines are seen for amino acid polymorphisms at the Hex-C and Hex-t2 loci. There is evidence for accelerated amino acid substitution in Hex-t1 that has lost residues known to be associated with glucose and glucose-6-phosphate binding. D. simulans shows substantial linkage phase structuring that suggests historical population subdivision.     

Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans

Surveys of molecular variation in Drosophila melanogaster and Drosophila simulans have suggested that diversity outside of Africa is a subset of that within Africa. It has been argued that reduced levels of diversity in non-African populations reflect a population bottleneck, adaptation to temperate climates, or both. Here, I summarize the available single-nucleotide polymorphism data for both species. A simple "out of Africa" bottleneck scenario is consistent with geographic patterns for loci on the X chromosome but not with loci on the autosomes. Interestingly, there is a trend toward lower nucleotide diversity on the X chromosome relative to autosomes in non-African populations of D. melanogaster, but the opposite trend is seen in African populations. In African populations, autosomal inversion polymorphisms in D. melanogaster may contribute to reduced autosome diversity relative to the X chromosome. To elucidate the role that selection might play in shaping patterns of variability, I present a summary of within- and between-species patterns of synonymous and replacement variation in both species. Overall, D. melanogaster autosomes harbor an excess of amino acid replacement polymorphisms relative to D. simulans. Interestingly, range expansion from Africa appears to have had little effect on synonymous-to-replacement polymorphism ratios.     

Molecular evolution of the Est-6 gene in Drosophila melanogaster: contrasting patterns of DNA variability in adjacent functional regions

We have investigated nucleotide polymorphism at the esterase 6 gene (Est-6) gene, including the complete coding region (1686 bp), as well as the 5'-flanking (1183 bp) and 3'-flanking (193 bp) regions of the gene, in 30 strains of Drosophila melanogaster and in one strain of Drosophila simulans. The level of silent variation is similar in the coding and in the 3'-flanking region, but smaller in the 5'-flanking region. Strong linkage disequilibrium occurs within each region; and also, although less pronounced, between the 5'-flanking region and the rest of the gene, including the 3'-flanking region. We suggest that the pattern of nucleotide polymorphism of Est-6 may be shaped by: (1) directional and balancing selection acting on the promoter and the coding region; and (2) interactions between the two regions that involve variable degrees of hitchhiking. The patterns of linkage disequilibrium, as well as the statistics Z(nS) (Genetics 146 (1997)1197) and B and Q (Genet. Res. 74 (1999) 65), may be interpreted as there being multiple targets of selection within the gene. The previously reported Est-6 allozyme latitudinal clines may be accounted for by the interaction between selective processes in the promoter and coding regions.     

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.     

Nucleotide variation in the triosephosphate isomerase (Tpi) locus of Drosophila melanogaster and Drosophila simulans

DNA sequence variation in a 1.1-kb region including the coding portion of the Tpi locus was examined in 25 homozygous third-chromosome lines of Drosophila melanogaster, nine lines of Drosophila simulans, and one line of Drosophila yakuba. Our data show that the widespread allozyme polymorphism observed in cosmopolitan D. melanogaster is due to a glutamic acid substitution occurring in a phylogenetically conserved lysine that has been identified as part of the "hinged-lid" active site of the enzyme. This observation suggests that the replacement polymorphism may have important functional consequences. One replacement polymorphism was also observed in D. simulans, although its functional relevance is more difficult to assess, since it affects a site that is not strongly conserved. This amino acid change in D. simulans is associated with a single lineage possessing seven unique silent substitutions, which may be indicative of balancing selection or population subdivision. The absence of fixed amino acid differences between D. melanogaster and D. simulans and only a single difference with D. yakuba suggests that triose phosphate isomerase is under strong functional constraint. Silent variation is slightly higher for D. melanogaster than for D. simulans. Finally, we outline the general lack of evidence for old balanced polymorphisms at allozyme loci in D. melanogaster.      

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.     

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 polymorphism in Drosophila melanogaster and D. simulans: what have we learned from recent studies?

We present a review of recent studies of molecular polymorphism in Drosophila melanogaster and D. simulans. The availability of African and non-African samples for these species makes it possible to compare microsatellite and DNA sequence polymorphism between these species, both inside and outside their native regions. There are four main results: (i) variability is larger in African populations; (ii) variation is usually higher on the autosomes, except for African D. melanogaster; (iii) DNA sequence variation is higher on D. simulans than on D. melanogaster autosomes; (iv) the ratio of replacement to silent polymorphisms is higher for D. melanogaster autosomal loci. We summarize the main hypotheses put forward to explain these results.     

The effect of developmental temperature on the genetic architecture underlying size and thermal clines in Drosophila melanogaster and D. simulans from the east coast of Australia

Body size and thermal tolerance clines in Drosophila melanogaster occur along the east coast of Australia. However the extent to which temperature affects the genetic architecture underlying the observed clinal divergence remains unknown. Clinal variation in these traits is associated with cosmopolitan chromosome inversions that cline in D. melanogaster. Whether this association influences the genetic architecture for these traits in D. melanogaster is unclear. Drosophila simulans shows linear clines in body size, but nonlinear clines in cold resistance. Clinally varying inversions are absent in D. simulans. Line-cross and clinal analyses were performed between tropical and temperate populations of D. melanogaster and D. simulans from the east coast of Australia to investigate whether clinal patterns and genetic effects contributing to clinal divergence in wing centroid size, thorax length, wing-to-thorax ratio, cold and heat resistance differed under different developmental temperatures (18 °C, 25 °C, and 29 °C). Developmental temperature influenced the genetic architecture in both species. Similarities between D. melanogaster and D. simulans suggest clinally varying inversion polymorphisms have little influence on the genetic architecture underlying clinal divergence in size in D. melanogaster. Differing genetic architectures across different temperatures highlight the need to consider different environments in future evolutionary and molecular studies of phenotypic divergence.     

Kinetics and thermodynamics of ethanol oxidation catalyzed by genetic variants of the alcohol dehydrogenase from Drosophila melanogaster and D. simulans

Four naturally occurring variants of the alcohol dehydrogenase enzyme (ADH; EC 1.1.1.1) from Drosophila melanogaster and D. simulans, with different primary structures, have been subjected to kinetic studies of ethanol oxidation at five temperatures. Two amino acid replacements in the N-terminal region which distinguish the ADH of D. simulans from the three ADH allozymes of D. melanogaster generate a significantly different activation enthalpy and entropy, and Gibbs free energy change. The one or two amino acid replacements in the C-terminal region between the ADH allozymes of D. melanogaster do not have such clear-cut effects. All four ADH variants show highly negative activation entropies. Sarcosine oxidation by the ADH-71k variant of D. melanogaster has an activation energy barrier similar to that of ethanol oxidation. Three amino acid differences between the ADH of D. simulans and the ADH-F variant of D. melanogaster influence the kappa cat and kappa cat/Kethm constant by a maximum factor of about 2 and 2.5, respectively, over the whole temperature range. Product inhibition patterns suggest a 'rapid equilibrium random' mechanism of ethanol oxidation by the ADH-71k, and the ADH of D. simulans.     

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.