On the Long-term Stability of Clines in Some Metabolic Genes in Drosophila melanogaster.

Very little information exists for long-term changes in genetic variation in natural populations. Here we take the unique opportunity to compare a set of data for SNPs in 15 metabolic genes from eastern US collections of Drosophila melanogaster that span a large latitudinal range and represent two collections separated by 12 to 13 years. We also expand this to a 22-year interval for the Adh gene and approximately 30 years for the G6pd and Pgd genes. During these intervals, five genes showed a statistically significant change in average SNP allele frequency corrected for latitude. While much remains unchanged, we see five genes where latitudinal clines have been lost or gained and two where the slope significantly changes. The long-term frequency shift towards a southern favored Adh S allele reported in Australia populations is not observed in the eastern US over a period of 21 years. There is no general pattern of southern-favored or northern-favored alleles increasing in frequency across the genes. This observation points to the fluid nature of some allelic variation over this time period and the action of selective responses or migration that may be more regional than uniformly imposed across the cline.

New views on the selection acting on genetic polymorphism in central metabolic genes.

Studies of the polymorphism of central metabolic genes as a source of fitness variation in natural populations date back to the discovery of allozymes in the 1960s. The unique features of these genes and their enzymes and our knowledge base greatly facilitates the systems-level study of this group. The expectation that pathway flux control is central to understanding the molecular evolution of genes is discussed, as well as studies that attempt to place gene-specific molecular evolution and polymorphism into a context of pathway and network architecture. There is an increasingly complex picture of the metabolic genes assuming additional roles beyond their textbook anabolic and catabolic reactions. In particular, this review emphasizes the potential role of these genes as part of the energy-sensing machinery. It is underscored that the concentrations of key cellular metabolites are the reflections of cellular energy status and nutritional input. These metabolites are the top-down signaling messengers that set signaling through signaling pathways that are involved in energy economy. I propose that the polymorphisms in central metabolic genes shift metabolite concentrations and in that fashion act as genetic modifiers of the energy-state coupling to the transcriptional networks that affect physiological trade-offs with significant fitness consequences.

Genetic perturbation of key central metabolic genes extends lifespan in Drosophila and affects response to dietary restriction.

There is a connection between nutrient inputs, energy-sensing pathways, lifespan variation and aging. Despite the role of metabolic enzymes in energy homeostasis and their metabolites as nutrient signals, little is known about how their gene expression impacts lifespan. In this report, we use P-element mutagenesis in Drosophila to study the effect on lifespan of reductions in expression of seven central metabolic enzymes, and contrast the effects on normal diet and dietary restriction. The major observation is that for five of seven genes, the reduction of gene expression extends lifespan on one or both diets. Two genes are involved in redox balance, and we observe that lower activity genotypes significantly extend lifespan. The hexokinases also show extension of lifespan with reduced gene activity. Since both affect the ATP/ADP ratio, this connects with the role of AMP-activated protein kinase as an energy sensor in regulating lifespan and mediating caloric restriction. These genes possess significant expression variation in natural populations, and our experimental genotypes span this level of natural activity variation. Our studies link the readout of energy state with the perturbation of the genes of central metabolism and demonstrate their effect on lifespan.

Variation in Drosophila melanogaster central metabolic genes appears driven by natural selection both within and between populations.

In this report, we examine the hypothesis that the drivers of latitudinal selection observed in the eastern US Drosophila melanogaster populations are reiterated within seasons in a temperate orchard population in Pennsylvania, USA. Specifically, we ask whether alleles that are apparently favoured in northern populations are also favoured early in the spring, and decrease in frequency from the spring to autumn with the population expansion. We use SNP data collected for 46 metabolic genes and 128 SNPs representing the central metabolic pathway and examine for the aggregate SNP allele frequencies whether the association of allele change with latitude and that with increasing days of spring-autumn season are reversed. Testing by random permutation, we observe a highly significant negative correlation between these associations that is consistent with this expectation. This correlation is stronger when we confine our analysis to only those alleles that show significant latitudinal changes. This pattern is not caused by association with chromosomal inversions. When data are resampled using SNPs for amino acid change the relationship is not significant but is supported when SNPs associated with cis-expression are only considered. Our results suggest that climate factors driving latitudinal molecular variation in a metabolic pathway are related to those operating on a seasonal level within populations.

A small system--high-resolution study of metabolic adaptation in the central metabolic pathway to temperate climates in Drosophila melanogaster.

In this article, we couple the geographic variation in 127 single-nucleotide polymorphism (SNP) frequencies in genes of 46 enzymes of central metabolism with their associated cis-expression variation to predict latitudinal or climatic-driven gene expression changes in the metabolic architecture of Drosophila melanogaster. Forty-two percent of the SNPs in 65% of the genes show statistically significant clines in frequency with latitude across the 20 local population samples collected from southern Florida to Ontario. A number of SNPs in the screened genes are also associated with significant expression variation within the Raleigh population from North Carolina. A principal component analysis of the full variance-covariance matrix of latitudinal changes in SNP-associated standardized gene expression allows us to identify those major genes in the pathway and its associated branches that are likely targets of natural selection. When embedded in a central metabolic context, we show that these apparent targets are concentrated in the genes of the upper glycolytic pathway and pentose shunt, those controlling glycerol shuttle activity, and finally those enzymes associated with the utilization of glutamate and pyruvate. These metabolites possess high connectivity and thus may be the points where flux balance can be best shifted. We also propose that these points are conserved points associated with coupling energy homeostasis and energy sensing in mammals. We speculate that the modulation of gene expression at specific points in central metabolism that are associated with shifting flux balance or possibly energy-state sensing plays a role in adaptation to climatic variation.

The intensity of selection acting on the couch potato gene--spatial-temporal variation in a diapause cline.

Cosmopolitan populations of Drosophila melanogaster have co-opted a form of reproductive diapause to overwinter in northern populations. Polymorphism in the couch potato gene has been implicated in genetic variation for this diapause trait. Using a collection of 20 populations from Florida to Canada and 11 collections from 3 years in a Pennsylvania orchard, we estimated the allele frequencies for 15 single nucleotide polymorphisms (SNPs) in the couch potato gene. These include the specific polymorphism associated with diapause inducability. We find that the SNP polymorphism, 48034(A/T), is correlated with latitude and its frequencies are predicted by the incidence of diapause trait. We find that the clinal patterns for cpo SNPs sampled in 1997 are similar to the same SNPs sampled in 2009-2010. SNPs that show apparent associations with cpo expression are also clinal with the low-expression allele increasing in frequency, as would be predicted from functional knockout studies of cpo. Finally, we see a significant pattern where the frequency of the diapause-causing allele drops in frequency during the summer season, consistent with the drop in the incidence of the diapause trait. The selection required to drive this response is large, roughly 24% to 59% per generation depending on the degree of dominance.

Molecular population genetics and selection in the glycolytic pathway.

In this review, I discuss the evidence for differential natural selection acting across enzymes in the glycolytic pathway in Drosophila. Across the genome, genes evolve at very different rates and possess markedly varying levels of molecular polymorphism, codon bias and expression variation. Discovering the underlying causes of this variation has been a challenge in evolutionary biology. It has been proposed that both the intrinsic properties of enzymes and their pathway position have direct effects on their molecular evolution, and with the genomic era the study of adaptation has been taken to the level of pathways and networks of genes and their products. Of special interest have been the energy-producing pathways. Using both population genetic and experimental approaches, our laboratory has been engaged in a study of molecular variation across the glycolytic pathway in Drosophila melanogaster and its close relatives. We have observed a pervasive pattern in which genes at the top of the pathway, especially around the intersection at glucose 6-phosphate, show evidence for both contemporary selection, in the form of latitudinal allele clines, and inter-specific selection, in the form of elevated levels of amino acid substitutions between species. To further explore this question, future work will require corroboration in other species, expansion into tangential pathways, and experimental work to better characterize metabolic control through the pathway and to examine the pleiotropic effects of these genes on other traits and fitness components.

Quantifying interactions within the NADP(H) enzyme network in Drosophila melanogaster.

In this report, we use synthetic, activity-variant alleles in Drosophila melanogaster to quantify interactions across the enzyme network that reduces nicotinamide adenine dinucleotide phosphate (NADP) to NADPH. We examine the effects of large-scale variation in isocitrate dehydrogenase (IDH) or glucose-6-phosphate dehydrogenase (G6PD) activity in a single genetic background and of smaller-scale variation in IDH, G6PD, and malic enzyme across 10 different genetic backgrounds. We find significant interactions among all three enzymes in adults; changes in the activity of any one source of a reduced cofactor generally result in changes in the other two, although the magnitude and directionality of change differs depending on the gene and the genetic background. Observed interactions are presumably through cellular mechanisms that maintain a homeostatic balance of NADPH/NADP, and the magnitude of change in response to modification of one source of reduced cofactor likely reflects the relative contribution of that enzyme to the cofactor pool. Our results suggest that malic enzyme makes the largest single contribution to the NADPH pool, consistent with the results from earlier experiments in larval D. melanogaster using naturally occurring alleles. The interactions between all three enzymes indicate functional interdependence and underscore the importance of examining enzymes as components of a network.

Direct evidence that genetic variation in glycerol-3-phosphate and malate dehydrogenase genes (Gpdh and Mdh1) affects adult ethanol tolerance in Drosophila melanogaster.

Many studies of alcohol adaptation in Drosophila melanogaster have focused on the Adh polymorphism, yet the metabolic elimination of alcohol should involve many enzymes and pathways. Here we evaluate the effects of glycerol-3-phosphate dehydrogenase (Gpdh) and cytosolic malate dehydrogenase (Mdh1) genotype activity on adult tolerance to ethanol. We have created a set of P-element-excision-derived Gpdh, Mdh1, and Adh alleles that generate a range of activity phenotypes from full to zero activity. Comparisons of paired Gpdh genotypes possessing 10 and 60% normal activity and 66 and 100% normal activity show significant effects where higher activity increases tolerance. Mdh1 null allele homozygotes show reductions in tolerance. We use piggyBac FLP-FRT site-specific recombination to create deletions and duplications of Gpdh. Duplications show an increase of 50% in activity and an increase of adult tolerance to ethanol exposure. These studies show that the molecular polymorphism associated with GPDH activity could be maintained in natural populations by selection related to adaptation to alcohols. Finally, we examine the interactions between activity genotypes for Gpdh, Mdh1, and Adh. We find no significant interlocus interactions. Observations on Mdh1 in both Gpdh and Adh backgrounds demonstrate significant increases in ethanol tolerance with partial reductions (50%) in cytosolic MDH activity. This observation strongly suggests the operation of pyruvate-malate and, in particular, pyruvate-citrate cycling in adaptation to alcohol exposure. We propose that an understanding of the evolution of tolerance to alcohols will require a system-level approach, rather than a focus on single enzymes.

An amino acid polymorphism in the couch potato gene forms the basis for climatic adaptation in Drosophila melanogaster.

Diapause is the classic adaptation to seasonality in arthropods, and its expression can result in extreme lifespan extension as well as enhanced resistance to environmental challenges. Little is known about the underlying evolutionary genetic architecture of diapause in any organism. Drosophila melanogaster exhibits a reproductive diapause that is variable within and among populations; the incidence of diapause increases with more temperate climates and has significant pleiotropic effects on a number of life history traits. Using quantitative trait mapping, we identified the RNA-binding protein encoding gene couch potato (cpo) as a major genetic locus determining diapause phenotype in D. melanogaster and independently confirmed this ability to impact diapause expression through genetic complementation mapping. By sequencing this gene in samples from natural populations we demonstrated through linkage association that variation for the diapause phenotype is caused by a single Lys/Ile substitution in one of the six cpo transcripts. Complementation analyses confirmed that the identified amino acid variants are functionally distinct with respect to diapause expression, and the polymorphism also shows geographic variation that closely mirrors the known latitudinal cline in diapause incidence. Our results suggest that a naturally occurring amino acid polymorphism results in the variable expression of a diapause syndrome that is associated with the seasonal persistence of this model organism in temperate habitats.

Sperm polymorphism within the sea urchin Strongylocentrotus droebachiensis: divergence between Pacific and Atlantic oceans.

The rapid evolution of traits related to fertilization such as sperm morphology may be pivotal in the evolution of reproductive barriers and speciation. The sea urchin Strongylocentrotus droebachiensis has a circumarctic distribution and shows substantial genetic subdivision between northeastern Atlantic populations and northwestern Atlantic and Pacific populations. Using transmission electron microscopy, we show here that sperm shape, size, and ultrastructure differ markedly among populations of S. droebachiensis from different oceans and reflect patterns of genetic divergence. Sperm nuclei from northwestern Atlantic and Pacific populations were longer and narrower than those from the northeastern Atlantic. We additionally demonstrate population-level differences in the amount and location of filamentous actin (F-actin) prior to the occurrence of the acrosome reaction. Sperm from Pacific and northwest Atlantic populations differed from that of all other echinoids examined in that intact sperm contains a partly preformed acrosomal process, a structure more closely resembling the acrosomal rod seen in some molluscs. Immunofluorescent studies using anti-bindin antibodies and the F-actin-specific stain phalloidin confirmed these findings. Divergence of reproductive traits such as sperm morphology may be related to divergence in gamete compatibility and genetic divergence, and could represent the first stages of speciation in free-spawning marine invertebrates.

Flux control and excess capacity in the enzymes of glycolysis and their relationship to flight metabolism in Drosophila melanogaster.

An important question in evolutionary and physiological genetics is how the control of flux-base phenotypes is distributed across the enzymes in a pathway. This control is often related to enzyme-specific levels of activity that are reported to be in excess of that required for demand. In glycolysis, metabolic control is frequently considered vested in classical regulatory enzymes, each strongly displaced from equilibrium. Yet the contribution of individual steps to control is unclear. To assess enzyme-specific control in the glycolytic pathway, we used P-element excision-derived mutagenesis in Drosophila melanogaster to generate full and partial knockouts of seven metabolic genes and to measure tethered flight performance. For most enzymes, we find that reduction to half of the normal activity has no measurable impact on wing beat frequency. The enzymes catalyzing near-equilibrium reactions, phosphoglucose isomerase, phosphoglucomutase, and triosephosphate isomerase fail to show any decline in flight performance even when activity levels are reduced to 17% or less. At reduced activities, the classic regulatory enzymes, hexokinase and glycogen phosphorylase, show significant drops in flight performance and are nearer to saturation. Our results show that flight performance is canalized or robust to the activity variation found in natural populations. Furthermore, enzymes catalyzing near-equilibrium reactions show strong genetic dominance down to low levels of activity. This implies considerable excess enzyme capacity for these enzymes.

Triglyceride pools, flight and activity variation at the Gpdh locus in Drosophila melanogaster.

We have created a set of P-element excision-derived Gpdh alleles that generate a range of GPDH activity phenotypes ranging from zero to full activity. By placing these synthetic alleles in isogenic backgrounds, we characterize the effects of minor and major activity variation on two different aspects of Gpdh function: the standing triglyceride pool and glycerol-3-phosphate shuttle-assisted flight. We observe small but statistically significant reductions in triglyceride content for adult Gpdh genotypes possessing 33-80% reductions from normal activity. These small differences scale to a notable proportion of the observed genetic variation in triglyceride content in natural populations. Using a tethered fly assay to assess flight metabolism, we observed that genotypes with 100 and 66% activity exhibited no significant difference in wing-beat frequency (WBF), while activity reductions from 60 to 10% showed statistically significant reductions of approximately 7% in WBF. These studies show that the molecular polymorphism associated with GPDH activity could be maintained in natural populations by selection in the triglyceride pool.

Geographic variation in diapause incidence, life-history traits, and climatic adaptation in Drosophila melanogaster.

In Drosophila melanogaster, exposure of females to low temperature and shortened photoperiod can induce the expression of reproductive quiescence or diapause. Diapause expression is highly variable within and among natural populations and has significant effects on life-history profiles, including patterns of longevity, fecundity, and stress resistance. We hypothesized that if diapause expression is associated with overwintering mechanisms and adaptation to temperate environments, the frequency of diapause incidence would exhibit a latitudinal cline among natural populations. Because stress resistance and reproductive traits are also clinal in this species, we also examined how patterns of fecundity and longevity varied with geography and how stress resistance and associated traits differed constitutively between diapause and nondiapause lines. Diapause incidence was shown to vary predictably with latitude, ranging from 35% to 90% among natural populations in the eastern United States Survivorship under starvation stress differed between diapause and nondiapause lines; diapause phenotypes were also distinct for total body triglyceride content and the developmental distribution of oocytes in the ovary following stress exposure. Patterns of longevity, fecundity, and ovariole number also varied with geography. The data suggest that, for North American populations, diapause expression is functionally associated with overwintering mechanisms and may be an integral life-history component in natural populations.

Natural and synthetic alleles provide complementary insights into the nature of selection acting on the Men polymorphism of Drosophila melanogaster.

Two malic enzyme alleles, Men(113A) and Men(113G), occur at approximately equal frequency in North American populations of Drosophila melanogaster, while only Men(113A) occurs in African populations. We investigated the population genetics, biochemical characteristics, and selective potential of these alleles. Comparable levels of nucleotide polymorphism in both alleles suggest that the Men(113G) allele is not recently derived, but we find no evidence in the DNA sequence data for selection maintaining the polymorphism. Interestingly, the alleles differ in both V(max) and K(m) for the substrate malate. Triglyceride concentration and isocitrate dehydrogenase (IDH) and glucose-6-phosphate dehydrogenase (G6PD) activities are negatively correlated with the in vivo activities of the Men alleles. We examined the causality of the observed correlations using P-element excision-derived knockout alleles of the Men gene and found significant changes in the maximum activities of both IDH and G6PD, but not in triglyceride concentration, suggesting compensatory interactions between MEN, IDH, and G6PD. Additionally, we found significantly higher than expected levels of MEN activity in knockout heterozygotes, which we attribute to transvection effects. The distinct differences in biochemistry and physiology between the naturally occurring alleles and between the engineered alleles suggest the potential for selection on the Men locus.

The structure and population genetics of the breakpoints associated with the cosmopolitan chromosomal inversion In(3R)Payne in Drosophila melanogaster.

We report here the breakpoint structure and sequences of the Drosophila melanogaster cosmopolitan chromosomal inversion In(3R)P. Combining in situ hybridization to polytene chromosomes and long-range PCR, we have identified and sequenced the distal and proximal breakpoints. The breakpoints are not simple cut-and-paste structures; gene fragments and small duplications of DNA are associated with both breaks. The distal breakpoint breaks the tolkin (tok) gene and the proximal breakpoint breaks CG31279 and the tolloid (tld) gene. Functional copies of all three genes are found at the opposite breakpoints. We sequenced a representative sample of standard (St) and In(3R)P karyotypes for a 2-kb portion of the tok gene, as well as the same 2 kb from the pseudogene tok fragment found at the distal breakpoint of In(3R)P chromosomes. The tok gene in St arrangements possesses levels of polymorphism typical of D. melanogaster genes. The functional tok gene associated with In(3R)P shows little polymorphism. Numerous single-base changes, as well as deletions and duplications, are associated with the truncated copy of tok. The overall pattern of polymorphism is consistent with a recent origin of In(3R)P, on the order of Ne generations. The identification of these breakpoint sequences permits a simple PCR-based screen for In(3R)P.

Single-locus latitudinal clines and their relationship to temperate adaptation in metabolic genes and derived alleles in Drosophila melanogaster.

We report a study in Drosophila melanogaster of latitudinal clines for 23 SNPs embedded in 13 genes (Pgi, Gapdh1, UGPase, Pglym78, Pglym87, Eno, Men, Gdh, Sod, Pgk, Mdh1, TreS, Treh) representing various metabolic enzymes. Our samples are from 10 populations spanning latitude from southern Florida to northern Vermont. Three new clines with latitude were detected. These are the amino acid polymorphisms in the NAD-dependent glutamate dehydrogenase (Gdh) and trehalase (Treh) genes, and a silent site polymorphism in the UDP-glucose pyrophosphorylase gene (UGPase). The result, when combined with the overall incidence and pattern of reports for six other genes (Adh, Gpdh, Pgm, G6pd, 6Pgd, Hex-C), presents a picture of latitudinal clines in metabolic genes prevalent around the branch point of competing pathways. For six of the seven amino acid polymorphisms showing significant latitudinal clines in North America, the derived allele is the one increasing with latitude, suggesting temperate adaptation. This is consistent with a model of an Afrotropical ancestral species adapting to temperate climates through selection favoring new mutations.

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.

Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila.

This study focuses on the population genetics of alcohol dehydrogenase (Adh) in cactophilic Drosophila. Drosophila mojavensis and D. arizonae utilize cactus hosts, and each host contains a characteristic mixture of alcohol compounds. In these Drosophila species there are two functional Adh loci, an adult form (Adh-2) and a larval and ovarian form (Adh-1). Overall, the greater level of variation segregating in D. arizonae than in D. mojavensis suggests a larger population size for D. arizonae. There are markedly different patterns of variation between the paralogs across both species. A 16-bp intron haplotype segregates in both species at Adh-2, apparently the product of an ancient gene conversion event between the paralogs, which suggests that there is selection for the maintenance of the intron structure possibly for the maintenance of pre-mRNA structure. We observe a pattern of variation consistent with adaptive protein evolution in the D. mojavensis lineage at Adh-1, suggesting that the cactus host shift that occurred in the divergence of D. mojavensis from D. arizonae had an effect on the evolution of the larval expressed paralog. Contrary to previous work we estimate a recent time for both the divergence of D. mojavensis and D. arizonae (2.4 +/- 0.7 MY) and the age of the gene duplication (3.95 +/- 0.45 MY).

INSIGHT INTO SPECIATION FROM HISTORICAL DEMOGRAPHY IN THE PHYTOPHAGOUS BEETLE GENUS OPHRAELLA.

Speciation in phytophagous insects is commonly associated with shifts in host use. Using a phylogenetic framework to identify recently diverged taxa that have undergone a radical host shift, this study focuses on how reconstruction of the historical demography of a species, in conjunction with branching patterns between species, provides insight into mode of speciation. Analyses of mitochondrial cytochrome oxidase I sequences indicate that the leaf beetle Ophraella communa exhibits significant population structure, as shown by patterns of genealogical relationships among mitochondrial haplotypes and high FST -values. However, the absence of regional localization of old clades of haplotypes, negative Tajima's D, and unimodal rather than bimodal frequency distribution of the number of pairwise differences between sequences suggests an absence of long-term barriers to gene flow. Furthermore, we found no evidence of isolation by distance. This pattern of genetic variation is consistent with episodes of gene flow on a large geographic scale, perhaps owing to Pleistocene changes in climate. Ophraella communa and its sister species O. bilineata diverged during the early Pleistocene. The evidence of dynamic population structure in O. communa, potentially including episodic but massive gene flow, suggests that reproductive isolation evolved quite rapidly on a localized geographic scale, because speciation would probably have been reversed by gene flow if the evolution of reproductive isolation had been prolonged. That is, gene flow occasioned by range shifts during the Pleistocene would likely have interrupted speciation unless it occurred very rapidly. Sequence diversity implies a large effective population size (> 106 ) in both O. communa and O. bilineata. However, a model based on a drastic bottleneck did not have a lower likelihood than a model with no bottleneck, simply because the time since speciation has been great enough for coalescence to a single ancestor that existed after the speciation event. Sequence diversity in itself, without reference to the time since speciation, cannot provide evidence on the demography of speciation.