Adaptive evolution of metabolic pathways in Drosophila.

The adaptive significance of enzyme variation has been of central interest in population genetics. Yet, how natural selection operates on enzymes in the larger context of biochemical pathways has not been broadly explored. A basic expectation is that natural selection on metabolic phenotypes will target enzymes that control metabolic flux, but how adaptive variation is distributed among enzymes in metabolic networks is poorly understood. Here, we use population genetic methods to identify enzymes responding to adaptive selection in the pathways of central metabolism in Drosophila melanogaster and Drosophila simulans. We report polymorphism and divergence data for 17 genes that encode enzymes of 5 metabolic pathways that converge at glucose-6-phosphate (G6P). Deviations from neutral expectations were observed at five loci. Of the 10 genes that encode the enzymes of glycolysis, only aldolase (Ald) deviated from neutrality. The other 4 genes that were inconsistent with neutral evolution (glucose-6-phosphate dehydrogenase [G6pd]), phosphoglucomutase [Pgm], trehalose-6-phosphate synthetase [Tps1], and glucose-6phosphatase [G6pase] encode G6P branch point enzymes that catalyze reactions at the entry point to the pentose-phosphate, glycogenic, trehalose synthesis, and gluconeogenic pathways. We reconcile these results with population genetics theory and existing arguments on metabolic regulation and propose that the incidence of adaptive selection in this system is related to the distribution of flux control. The data suggest that adaptive evolution of G6P branch point enzymes may have special significance in metabolic adaptation.

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.

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster.

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.

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.

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.

Adaptive evolution of a candidate gene for aging in Drosophila.

Examination of the phenotypic effects of specific mutations has been extensively used to identify candidate genes affecting traits of interest. However, such analyses do not reveal anything about the evolutionary forces acting at these loci, or whether standing allelic variation contributes to phenotypic variance in natural populations. The Drosophila gene methuselah (mth) has been proposed as having major effects on organismal stress response and longevity phenotype. Here, we examine patterns of polymorphism and divergence at mth in population level samples of Drosophila melanogaster, D. simulans, and D. yakuba. Mth has experienced an unusually high level of adaptive amino acid divergence concentrated in the intra- and extracellular loop domains of the receptor protein, suggesting the historical action of positive selection on those regions of the molecule that modulate signal transduction. Further analysis of single nucleotide polymorphisms (SNPs) in D. melanogaster provided evidence for contemporary and spatially variable selection at the mth locus. In ten surveyed populations, the most common mth haplotype exhibited a 40% cline in frequency that coincided with population level differences in multiple life-history traits including lifespan. This clinal pattern was not associated with any particular SNP in the coding region, indicating that selection is operating at a closely linked site that may be involved in gene expression. Together, these consistently nonneutral patterns of inter- and intraspecific variation suggest adaptive evolution of a signal transduction pathway that may modulate lifespan in nature.

Lack of nucleotide polymorphism in the Y-linked sperm flagellar dynein gene Dhc-Yh3 of Drosophila melanogaster and D. simulans.

We studied levels of intra- and interspecific nucleotide variation associated with a Y-linked gene in five members of the Drosophila melanogaster subgroup. Using published sequence for 348 bp of the Dhc-Yh3 gene, and degenerate PCR primers designed from comparisons of the sea urchin and Chlamydomonas flagellar dynein genes, we recovered a 1738-bp region in D. melanogaster. Analyses of sequence variation in a worldwide collection of 11 lines of D. melanogaster and 10 lines of D. simulans found only a single silent polymorphism in the latter species. The synonymous site divergence per site for Dhc-Yh3 is comparable to values for X and autosomal genes. Assuming a Wright-Fisher population model, the lack of variation is statistically less than expected using appropriately reduced estimates of theta from the X and autosomes. Because the Y chromosome encodes only six known genes, genetic hitchhiking associated with background selection is unlikely to explain this low variation. Conversely, adaptive hitchhiking, as associated with sex-ratio chromosomes, or a large variance in male fertility may reduce the polymorphism on the Y chromosome. Codon bias is very low, as seen for other genes in regions of low recombination.

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.

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.

Contrasting histories of three gene regions associated with In(3L)Payne of Drosophila melanogaster.

In the present report, we studied nucleotide variation in three gene regions of Drosophila melanogaster, spanning > 5 kb and showing different degrees of association with the cosmopolitan inversion In(3L)Payne. The analysis of sequence variation in the regions surrounding the breakpoints and the heat shock 83 (Hsp83) gene locus, located close to the distal breakpoint, revealed the absence of shared polymorphisms and the presence of a number of fixed differences between arrangements, indicating absence of genetic exchange. In contrast, for the esterase-6 gene region, located in the center of the inversion, we observed the presence of shared polymorphisms between arrangements suggesting genetic exchange. In the regions close to the breakpoints, the common St arrangement is 10 times more polymorphic than inverted chromosomes. We propose that the lack of recombination between arrangements in these regions coupled with genetic hitchhiking is the best explanation for the low heterozygosity observed in inverted lines. Using the data for the breakpoints, we estimate that this inversion polymorphism is around 0.36 million yr old. Although it is widely accepted that inversions are examples of balanced polymorphisms, none of the current neutrality tests including our Monte Carlo simulations showed significant departure from neutral expectations.

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.

Amino acid polymorphism and rare electrophoretic variants of G6PD from natural populations of Drosophila melanogaster.

Identifying the amino acid changes responsible for electrophoretic variants is essential to understanding the significance of allozyme polymorphism in adaptation. The amino acid mutations responsible for the common G6PD allozyme polymorphisms in Drosophila melanogaster have been recently described. This study characterizes the amino acid changes associated with 11 rare electrophoretic G6PD variants. The 11 rare electrophoretic variants result from six independent amino acid mutations. The in vivo function of the rare variants was determined in an earlier study and most variants fell into one of two function classes. It is shown here that the function of the rare variants reflects the state of the Pro/Leu mutation responsible for the A/B allozyme polymorphism in each variant. Two mutations destabilize quaternary structure resulting in shifts from tetrameric dimeric alleles, and one of these also results in a variant with in vivo function intermediate to A and B. That mutation is an aspartic-acid-to-asparagine change that is two residues away from the Pro/Leu polymorphism responsible for the A/B dimertetramer quaternary shift. Structure-function relationships based on studies of human G6PD deficiency-associated mutations predict that these last two amino acid changes fall within the protein domain responsible for NADP binding.

Isolation and analysis of the breakpoint sequences of chromosome inversion In(3L)Payne in Drosophila melanogaster.

Chromosomal rearrangements constitute a significant feature of genome evolution, and inversion polymorphisms in Drosophila have been studied intensely for decades. Population geneticists have long recognized that the sequence features associated with inversion breakpoints would reveal much about the mutational origin, uniqueness, and genealogical history of individual inversion polymorphisms, but the cloning of breakpoint sequences is not trivial. With the aid of a method for rapid recovery of DNA clones spanning rearrangement breakpoints, we recover and examine the DNA sequences spanning the breakpoints of the cosmopolitan inversion In(3L)Payne in Drosophila melanogaster. By examining the sequence diversity associated with six standard and seven inverted chromosomes from natural populations, we find that the inversion is monophyletic in origin, the sequences are genetically isolated from recombination at the breakpoints, and there is no association with features such as transposable elements. The inverted sequences show 17-fold less nucleotide polymorphism, but there are eight fixed differences in the region spanning both breakpoints. This suggests that this inversion is not recently derived. Finally, Northern analysis and transcript mapping find that the distal breakpoint has disrupted three transcripts that are normally expressed in the standard arrangement. Incidentally, the method introduced here can be used to isolate breakpoint sequences of arrangements associated with many human diseases.

Evidence for adaptive evolution of the G6pd gene in the Drosophila melanogaster and Drosophila simulans lineages.

Proponents of the neutral theory argue that evolution at the molecular level largely reflects a process of random genetic drift of neutral mutations. Under this theory, levels of interspecific divergence and intraspecific polymorphism are expected to be correlated across classes of nucleotide or amino acid sequences with different degrees of functional constraint, such as synonymous and replacement sites. Nucleotide sites with reduced polymorphism should show comparably reduced levels of interspecific divergence. To examine this hypothesis, we have sequenced 32 and 12 copies of the glucose-6-phosphate dehydrogenase (G6pd) gene in Drosophila melanogaster and Drosophila simulans, respectively. Both species exhibit similar levels of nucleotide polymorphism at synonymous sites. D. melanogaster shows two amino acid polymorphisms, one associated with the cosmopolitan allozyme polymorphism and a second with an allozyme polymorphism endemic to European and North African populations. In contrast, D. simulans shows no replacement polymorphism. While synonymous divergence between species is 10%, which is typical of other genes, there are 21 replacement differences. This level of amino acid sequence divergence, when contrasted with levels of amino acid polymorphism, silent polymorphism, and divergence, is in 10-fold excess over that expected under the neutral model of molecular evolution. We propose that this excess divergence reflects episodes of natural selection on G6pd resulting in fixation of advantageous amino acid mutations in these two recently separated lineages.

Direct measurement of in vivo flux differences between electrophoretic variants of G6PD from Drosophila melanogaster.

Demonstrating that naturally occurring enzyme polymorphisms significantly impact metabolic pathway flux is a fundamental step in examining the possible adaptive significance of such polymorphisms. In earlier studies of the glucose-6-phosphate dehydrogenase (G6PD) polymorphism in Drosophila melanogaster, we used two different methods, exploiting both genotype-dependent interactions with the 6Pgd locus, and conventional steady-state kinetics to examine activity differences between the two common allozymes. In this report we use 1-14C- and 6-14C-labeled glucose to estimate directly genotype-dependent flux differences through the pentose shunt. Our results show that G6pdA genotype possesses statistically lower pentose shunt flux than G6pdB at 25 degrees. We estimate this to be about a 32% reduction, which is consistent with the two former studies. These results reflect a significant responsiveness of pentose shunt flux to activity variation at the G6PD-catalyzed step, and predict that the G6PD allozymes generate a polymorphism for pentose shunt flux.

Accumulation of P elements in minority inversions in natural populations of Drosophila melanogaster.

The accumulation of a transposable element inside chromosomal inversions is examined theoretically by a mathematical model, and empirically by counts of P elements associated with inversion polymorphisms in natural populations of Drosophila melanogaster. The model demonstrates that, if heterozygosity for an inversion effectively reduces element associated production of detrimental chromosome rearrangements, a differential accumulation of elements is expected, with increased copy number inside the minority inversion. Several-fold differential accumulations are possible with certain parameter values. We present data on P element counts for inversion polymorphisms on all five chromosome arms of 157 haploid genomes from two African populations. Our observations show significantly increased numbers of elements within the regions associated with the least common, or minority arrangements, in natural inversion polymorphisms.

Comparison of in vitro and in vivo activities associated with the G6PD allozyme polymorphism in Drosophila melanogaster.

Earlier studies of the A and B allozymes at the G6pd locus show a differential ability of the genotypes to suppress the loss of viability associated with a low activity 6-phosphogluconate dehydrogenase mutation, 6Pgdlo1. This observation indicates a relatively lower activity for the A allozyme genotype, but it is not known if this level of suppression required a large difference in in vivo activity. To clarify this difference an analysis of the biochemical properties of the purified allozymes was carried out, as well as an analysis of the activity level associated with an original low activity P element-derived allele which had partially reverted and lost its suppression ability. G6PD activity and protein level were studied in 47 X chromosome lines from North America. The A genotype averages a 9% lower Vmax. From analysis of the correlation between G6PD activity and protein level it remains unclear whether the allozyme Vmax difference results from dissimilarity in protein level or kcat. At 25 degrees and physiological pH, comparative studies of the steady-state kinetics show the two purified allozyme variants differ significantly in their KM values for glucose-6-phosphate and NADP, and the K1 for NADPH. In aggregate these parameters predict the A genotype possesses a 20% lower in vitro catalytic efficiency. A partial revertant of a P element-derived low activity B variant, was shown to lose the ability to suppress 6Pgdlo1 low viability after acquiring only 60% of normal B activity. This last comparison shows the A genotype activity must be reduced in vivo by at least 40%.

Cloning regions of the Drosophila genome by microdissection of polytene chromosome DNA and PCR with nonspecific primer.

A simple and rapid procedure to isolate clones carrying sequences from a specific region of the polytene chromosome of Drosophila is demonstrated. The procedure involves microdissection of the region of interest, amplification of the DNA by PCR using a primer designed to prime the synthesis nonspecifically, labeling of the amplified DNA using the random primer method, and screening of a standard library with the probe to identify and isolate clones carrying sequences homologous to the dissected region. This procedure has the potential to replace the difficult procedure of microcloning, as well as facilitate chromosome walking.

Restriction-map variation with the yellow-achaete-scute region in five populations of Drosophila melanogaster.

It has been proposed that the degree of recombination for a genomic region will affect the level of both nucleotide heterozygosity and the density of transposable elements. Both features of genomic diversity have been examined in a number of recent reports for regions undergoing relatively normal levels of recombination in Drosophila melanogaster. In this study the genomic variation associated with yellow-achaete-scute loci located at the tip of the X chromosome is examined by six-cutter restriction mapping. In this region, as usual for regions adjacent to telomeres, crossing-over is dramatically reduced, and published studies of visible mutants indicate extremely little restriction-map variation. Eight six-cutter restriction endonucleases were used to locate sequence variation in 14- and 16.5-kb regions in 109 lines sampled from North America, Africa, and Europe. The overall level of heterozygosity is estimated as 0.29%. Nine large insertions, all presumed to be transposable elements, were observed. Base-pair heterozygosity appears to be reduced compared with regions having normal levels of recombination. The estimated heterozygosity is much higher than reported in earlier studies of restriction-map variation among visible mutations in the complex. The incidence of large insertions is not elevated compared with that in other regions of the genome. This suggests that asymmetric synapsis and exchange is not an important mechanism for the elimination of transposable elements.

Restriction-map variation associated with the G6PD polymorphism in natural populations of Drosophila melanogaster.

Restriction-map variation was studied in 126 copies of the G6pd region in X chromosome lines of Drosophila melanogaster from North America, Europe, and Africa. Special attention was focused on the distribution of variation relative to the geographically variable polymorphism for two electrophoretic variants. Nucleotide heterozygosity as determined by eight six-cutter restriction enzymes for the 13-kb region is estimated, on the basis of the worldwide sample, to be 0.065%, which is the lowest value reported for any comparable region in the D. melanogaster genome. Significant linkage disequilibrium between electrophoretic alleles and restriction-site variation is observed for several sites. In contrast to published studies of other genetic regions, there are large insertions that reach significant frequencies and are found across considerable geographic distances. There is a clustering of this variation inside the first large intervening sequence of the G6PD gene.