Phylogenetic analysis of glycolytic enzyme expression

Transformation to ischaemia tolerance of frog brain function corresponds to dynamic changes in mRNA co-expression across metabolic pathways

Neural activity is costly and requires continuous ATP from aerobic metabolism. Brainstem motor function of American bullfrogs normally collapses after minutes of ischaemia, but following hibernation, it becomes ischaemia-tolerant, generating output for up to 2 h without oxygen or glucose delivery. Transforming the brainstem to function during ischaemia involves a switch to anaerobic glycolysis and brain glycogen. We hypothesized that improving neural performance during ischaemia involves a transcriptional program for glycogen and glucose metabolism. Here we measured mRNA copy number of genes along the path from glycogen metabolism to lactate production using real-time quantitative PCR. The expression of individual genes did not reflect enhanced glucose metabolism. However, the number of co-expressed gene pairs increased early into hibernation, and by the end, most genes involved in glycogen metabolism, glucose transport and glycolysis exhibited striking linear co-expression. By contrast, co-expression of genes in the Krebs cycle and electron transport chain decreased throughout hibernation. Our results uncover reorganization of the metabolic transcriptional network associated with a shift to ischaemia tolerance in brain function. We conclude that modifying gene co-expression may be a critical step in synchronizing storage and use of glucose to achieve ischaemia tolerance in active neural circuits.

A unique neurogenomic state emerges after aggressive confrontations in males of the fish Betta splendens

Territorial defense involves frequent aggressive confrontations with competitors, but little is known about how brain-transcriptomic profiles change between individuals competing for territory establishment. Our previous study elucidated that when two fish Betta splendens males interact, transcriptomes across their brains synchronize in a way that reflects a mutual assessment process between them at the gene expression level. Here we aim to evaluate how the brain-transcriptomic profiles of opponents change immediately after shifting their social status (i.e., the winner/loser has emerged) and 30 min after this shift. We showed that changes in the expression of certain genes are unique to different fighting stages and the expression patterns of certain genes are transiently or persistently changed across all fighting stages. These brain transcriptomic responses are in accordance with behavioral changes across the fight. Strikingly, the specificity of the brain-transcriptomic synchronization of a pair during fighting was gradually lost after fighting ceased, leading to the emergence of a basal neurogenomic state in which the changes in gene expression were reduced to minimum and consistent across all individuals. This state shares common characteristics with the hibernation state that animals adopt to minimize their metabolic rates to save energy. Interestingly, expression changes for genes related to metabolism, autism spectrum disorder, and long-term memory still differentiated losers from winners. Together, the fighting system using male B. splendens provides a promising platform for investigating neurogenomic states of aggression in vertebrates.

Thermal adaptation in the honeybee (Apis mellifera) via changes to the structure of malate dehydrogenase

In honeybees there are three alleles of cytosolic malate dehydrogenase gene: F, M and S. Allele frequencies are correlated with environmental temperature, suggesting that the alleles have temperature-dependent fitness benefits. We determined the enzyme activity of each allele across a range of temperatures in vitro The F and S alleles have higher activity and are less sensitive to high temperatures than the M allele, which loses activity after incubation at temperatures found in the thorax of foraging bees in hot climates. Next, we predicted the protein structure of each allele and used molecular dynamics simulations to investigate their molecular flexibility. The M allozyme is more flexible than the S and F allozymes at 50°C, suggesting a plausible explanation for its loss of activity at high temperatures, and has the greatest structural flexibility at 15°C, suggesting that it can retain some enzyme activity at cooler temperatures. MM bees recovered from 2 h of cold narcosis significantly better than all other genotypes. Combined, these results explain clinal variation in malate dehydrogenase allele frequencies in the honeybee at the molecular level.

Evolutionary rate variation among genes involved in galactomannan biosynthesis in Coffea canephora

The endosperm cell walls of mature coffee seeds accumulate large amounts of mannan storage polysaccharides, which serve as nutrient reserve for embryo and contribute to beverage quality. Our study investigated the evolutionary patterns of key galactomannan (GM) biosynthesis genes using d N/d S ratio, synteny, and phylogenetic analysis and detected heterogeneity in rate of evolution among gene copies. Selection ratio index revealed evidence of positive selection in the branch editing gene Coffea canephora alpha (α) galactosidase (Cc-alpha Gal) at Cc11_g15950 copy (ω = 1.12), whereas strong purifying selection on deleterious mutations was observed in the Coffea canephora uridine diphosphate (UDP)-glucose 4'-epimerase (Cc-UG4E) and Coffea canephora mannose-1P guanylytransferase (Cc-MGT) genes controlling the crucial nucleotide carbon sugar building blocks flux in the pathway. Relatively low sequence diversity and strong syntenic linkages were detected in all GM pathway genes except in Cc-alpha Gal, which suggests a correlation between selection pressure and nucleotide diversity or synteny analysis. In addition, phylogenetic analysis revealed independent evolution or expansion of GM pathway genes in different plant species, with no obvious inferable clustering patterns according to either gene family or congruent with evolutionary plants lineages tested due to high dynamic nature and specific biochemical cell wall modification requirements. Altogether, our study shows a significant high rate of evolutionary variation among GM pathway genes in the diploid C. canephora and demonstrates the inherent variation in evolution of gene copies and their potential role in understanding selection rates in a homogenously connected metabolic pathway.

A quantitative framework for characterizing the evolutionary history of mammalian gene expression

The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein–Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.

Evolved genetic and phenotypic differences due to mitochondrial-nuclear interactions

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only pathway with both nuclear and mitochondrial encoded proteins. The importance of the interactions between these two genomes has recently received more attention because of their potential evolutionary effects and how they may affect human health and disease. In many different organisms, healthy nuclear and mitochondrial genome hybrids between species or among distant populations within a species affect fitness and OxPhos functions. However, what is less understood is whether these interactions impact individuals within a single natural population. The significance of this impact depends on the strength of selection for mito-nuclear interactions. We examined whether mito-nuclear interactions alter allele frequencies for ~11,000 nuclear SNPs within a single, natural Fundulus heteroclitus population containing two divergent mitochondrial haplotypes (mt-haplotypes). Between the two mt-haplotypes, there are significant nuclear allele frequency differences for 349 SNPs with a p-value of 1% (236 with 10% FDR). Unlike the rest of the genome, these 349 outlier SNPs form two groups associated with each mt-haplotype, with a minority of individuals having mixed ancestry. We use this mixed ancestry in combination with mt-haplotype as a polygenic factor to explain a significant fraction of the individual OxPhos variation. These data suggest that mito-nuclear interactions affect cardiac OxPhos function. The 349 outlier SNPs occur in genes involved in regulating metabolic processes but are not directly associated with the 79 nuclear OxPhos proteins. Therefore, we postulate that the evolution of mito-nuclear interactions affects OxPhos function by acting upstream of OxPhos.

Steep, coincident, and concordant clines in mitochondrial and nuclear-encoded genes in a hybrid zone between subspecies of Atlantic killifish, Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus-specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non-neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non-neutral nuclear clines, SNPs in two nuclear-encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin-binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non-unimodal, suggesting selection against advanced-generation hybrids, possibly due to effects on processes involved in oxidative metabolism.

Gene by environmental interactions affecting oxidative phosphorylation and thermal sensitivity

The oxidative phosphorylation (OxPhos) pathway is responsible for most aerobic ATP production and is the only metabolic pathway with proteins encoded by both nuclear and mitochondrial genomes. In studies examining mitonuclear interactions among distant populations within a species or across species, the interactions between these two genomes can affect metabolism, growth, and fitness, depending on the environment. However, there is little data on whether these interactions impact natural populations within a single species. In an admixed Fundulus heteroclitus population with northern and southern mitochondrial haplotypes, there are significant differences in allele frequencies associated with mitochondrial haplotype. In this study, we investigate how mitochondrial haplotype and any associated nuclear differences affect six OxPhos parameters within a population. The data demonstrate significant OxPhos functional differences between the two mitochondrial genotypes. These differences are most apparent when individuals are acclimated to high temperatures with the southern mitochondrial genotype having a large acute response and the northern mitochondrial genotype having little, if any acute response. Furthermore, acute temperature effects and the relative contribution of Complex I and II depend on acclimation temperature: when individuals are acclimated to 12°C, the relative contribution of Complex I increases with higher acute temperatures, whereas at 28°C acclimation, the relative contribution of Complex I is unaffected by acute temperature change. These data demonstrate a complex gene by environmental interaction affecting the OxPhos pathway.

Acclimation and acute temperature effects on population differences in oxidative phosphorylation

Temperature changes affect metabolism on acute, acclamatory, and evolutionary time scales. To better understand temperature's affect on metabolism at these different time scales, we quantified cardiac oxidative phosphorylation (OxPhos) in three Fundulus taxa acclimated to 12 and 28°C and measured at three acute temperatures (12, 20, and 28°C). The Fundulus taxa (northern Maine and southern Georgia F. heteroclitus, and a sister taxa, F. grandis) were used to identify evolved changes in OxPhos. Cardiac OxPhos metabolism was quantified by measuring six traits: state 3 (ADP and substrate-dependent mitochondrial respiration); E state (uncoupled mitochondrial activity); complex I, II, and IV activities; and LEAK ratio. Acute temperature affected all OxPhos traits. Acclimation only significantly affected state 3 and LEAK ratio. Populations were significantly different for state 3. In addition to direct effects, there were significant interactions between acclimation and population for complex I and between population and acute temperature for state 3. Further analyses suggest that acclimation alters the acute temperature response for state 3, E state, and complexes I and II: at the low acclimation temperature, the acute response was dampened at low assay temperatures, and at the high acclimation temperature, the acute response was dampened at high assay temperatures. Closer examination of the data also suggests that differences in state 3 respiration and complex I activity between populations were greatest between fish acclimated to low temperatures when assayed at high temperatures, suggesting that differences between the populations become more apparent at the edges of their thermal range.

Proteomic changes in the brain of the western painted turtle (Chrysemys picta bellii) during exposure to anoxia

During anoxia, overall protein synthesis is almost undetectable in the brain of the western painted turtle. The aim of this investigation was to address the question of whether there are alterations to specific proteins by comparing the normoxic and anoxic brain proteomes. Reductions in creatine kinase, hexokinase, glyceraldehyde-3-phosphate dehydrogenase, and pyruvate kinase reflected the reduced production of adenosine triphosphate (ATP) during anoxia while the reduction in transitional endoplasmic reticulum ATPase reflected the conservation of ATP or possibly a decrease in intracellular Ca(2+). In terms of neural protection programed cell death 6 interacting protein (PDCD6IP; a protein associated with apoptosis), dihydropyrimidinase-like protein, t-complex protein, and guanine nucleotide protein G(o) subunit alpha (Go alpha; proteins associated with neural degradation and impaired cognitive function) also declined. A decline in actin, gelsolin, and PDCD6IP, together with an increase in tubulin, also provided evidence for the induction of a neurological repair response. Although these proteomic alterations show some similarities with the crucian carp (another anoxia-tolerant species), there are species-specific responses, which supports the theory of no single strategy for anoxia tolerance. These findings also suggest the anoxic turtle brain could be an etiological model for investigating mammalian hypoxic damage and clinical neurological disorders.

Phylogenetic analyses: comparing species to infer adaptations and physiological mechanisms

Comparisons among species have been a standard tool in animal physiology to understand how organisms function and adapt to their surrounding environment. During the last two decades, conceptual and methodological advances from different fields, including evolutionary biology and systematics, have revolutionized the way comparative analyses are performed, resulting in the advent of modern phylogenetic statistical methods. This development stems from the realization that conventional analytical methods assume that observations are statistically independent, which is not the case for comparative data because species often resemble each other due to shared ancestry. By taking evolutionary history explicitly into consideration, phylogenetic statistical methods can account for the confounding effects of shared ancestry in interspecific comparisons, improving the reliability of standard approaches such as regressions or correlations in comparative analyses. Importantly, these methods have also enabled researchers to address entirely new evolutionary questions, such as the historical sequence of events that resulted in current patterns of form and function, which can only be studied with a phylogenetic perspective. Here, we provide an overview of phylogenetic approaches and their importance for studying the evolution of physiological processes and mechanisms. We discuss the conceptual framework underlying these methods, and explain when and how phylogenetic information should be employed. We then outline the difficulties and limitations inherent to comparative approaches and discuss potential problems researchers may encounter when designing a comparative study. These issues are illustrated with examples from the literature in which the incorporation of phylogenetic information has been useful, or even crucial, for inferences on how species evolve and adapt to their surrounding environment.

Nature's inordinate fondness for metabolic enzymes: why metabolic enzyme loci are so frequently targets of selection

Metabolic enzyme loci were some of the first genes accessible for molecular evolution and ecology research. New technologies now make the whole genome, transcriptome or proteome readily accessible, allowing unbiased scans for loci exhibiting significant differences in allele frequency or expression level and associated with phenotypes and/or responses to natural selection. With surprising frequency and in many cases in proportions greater than chance relative to other genes, glycolysis and TCA cycle enzyme loci appear among the genes with significant associations in these studies. Hence, there is an ongoing need to understand the basis for fitness effects of metabolic enzyme polymorphisms. Allele-specific effects on the binding affinity and catalytic rate of individual enzymes are well known, but often of uncertain significance because metabolic control theory and in vivo studies indicate that many individual metabolic enzymes do not affect pathway flux rate. I review research, so far little used in evolutionary biology, showing that metabolic enzyme substrates affect signalling pathways that regulate cell and organismal biology, and that these enzymes have moonlighting functions. To date there is little knowledge of how alleles in natural populations affect these phenotypes. I discuss an example in which alleles of a TCA enzyme locus associate with differences in a signalling pathway and development, organismal performance, and ecological dynamics. Ultimately, understanding how metabolic enzyme polymorphisms map to phenotypes and fitness remains a compelling and ongoing need for gaining robust knowledge of ecological and evolutionary processes.

Metabolic characteristics and response to high altitude in Phrynocephalus erythrurus (Lacertilia: Agamidae), a lizard dwell at altitudes higher than any other living lizards in the world

Metabolic response to high altitude remains poorly explored in reptiles. In the present study, the metabolic characteristics of Phrynocephaluserythrurus (Lacertilia: Agamidae), which inhabits high altitudes (4500 m) and Phrynocephalusprzewalskii (Lacertilia: Agamidae), which inhabits low altitudes, were analysed to explore the metabolic regulatory strategies for lizards living at high-altitude environments. The results indicated that the mitochondrial respiratory rates of P. erythrurus were significantly lower than those of P. przewalskii, and that proton leak accounts for 74~79% of state 4 and 7~8% of state3 in P. erythrurus vs. 43~48% of state 4 and 24~26% of state3 in P. przewalskii. Lactate dehydrogenase (LDH) activity in P. erythrurus was lower than in P. przewalskii, indicating that at high altitude the former does not, relatively, have a greater reliance on anaerobic metabolism. A higher activity related to β-hydroxyacyl coenzyme A dehydrogenase (HOAD) and the HOAD/citrate synthase (CS) ratio suggested there was a possible higher utilization of fat in P. erythrurus. The lower expression of PGC-1α and PPAR-γ in P. erythrurus suggested their expression was not influenced by cold and low PO2 at high altitude. These distinct characteristics of P. erythrurus are considered to be necessary strategies in metabolic regulation for living at high altitude and may effectively compensate for the negative influence of cold and low PO2.

Interindividual variation in Complex I activity in Fundulus heteroclitus along a steep thermocline

The first enzyme in the oxidative phosphorylation pathway is Complex I (EC 1.6.5.3). Complex I is a large heteromeric enzyme complex with 45 protein subunits that translocates H(+) ions across the mitochondrial inner membrane. Among northern and southern populations of the teleost fish Fundulus heteroclitus, Complex I subunits have fixed amino acid substitutions. Additionally, there are differences in oxidative phosphorylation activity among populations of F. heteroclitus. To investigate whether these differences are related to Complex I, enzyme activity was measured in 121 individuals from five populations of F. heteroclitus and its sister species Fundulus grandis acclimated to a constant 20°C temperature. Within each population, Complex I activity is highly variable among individuals of F. heteroclitus (coefficient of variation percentage among individuals has a mean of 90% in the five F. heteroclitus populations), and the mean Complex I activity among populations is significantly different at the latitudinal extremes of the range. Importantly, Complex I activity is more similar between F. heteroclitus from the southernmost population and its sister species F. grandis than to the northern populations of F. heteroclitus, suggesting important evolutionary differences. Unexpectedly, the activity is nearly fourfold higher in southern populations than northern populations. Mitochondrial density appears to compensate partially for decreased activity in northern individuals; activity per wet weight is only twofold higher in southern populations. We suggest that some of the variation in Complex I activity is genetically based and thus is being influenced by directional selection. However, this conclusion presents a conundrum: there should not be so much variation in Complex I activity within a population if this variation is biologically important.

Hypoxia tolerance in sculpins is associated with high anaerobic enzyme activity in brain but not in liver or muscle

We assessed hypoxia tolerance in 11 species of fish from the superfamily Cottoidea (commonly called sculpins) that are known to differ in their critical O(2) tensions (P(crit)) and examined whether hypoxia tolerance correlated with larger substrate stores and higher maximal activity of enzymes associated with anaerobic adenosine triphosphate production (especially glycolysis). Among the sculpins studied, there was large variation in time to loss of equilibrium (LOE(50)) at 6.4 ± 0.1 torr, with values ranging between 25 and 538 min, and the variation in LOE(50) was correlated with P(crit). Our measures of time to LOE(50) and P(crit) were regressed against maximal enzyme activities of lactate dehydrogenase (LDH), pyruvate kinase (PK), creatine phosphokinase (CPK), and citrate synthase (CS) as well as the concentrations of glycogen, glucose, and creatine phosphate in the brain, liver, and white muscle. In the brain, there was a phylogenetically independent relationship between P(crit) and tissue LDH, PK, CPK, and CS activities expressed relative to tissue mass. Hypoxia-tolerant sculpins (those with low P(crit) values) had higher levels of brain LDH, PK, CPK, and CS than did hypoxia-sensitive sculpins. Similarly, LOE(50) regressed against brain LDH, PK, and CPK activities expressed relative to tissue mass, with the more hypoxia-tolerant species (i.e., those with higher LOE(50)) having higher enzyme activities. However, when the phylogenetic relationship among our sculpins was taken into account, only the relationship between hypoxia tolerance and LDH activity remained significant. When enzyme activities were expressed relative to total soluble protein in the tissue, the only relationships that remained were between brain LDH activity and P(crit) and LOE(50). In liver and white muscle, there were no relationships between the measures of hypoxia tolerance and enzyme activity or metabolite content. Overall, our analysis suggests that hypoxia-tolerant sculpins maintain higher maximal activities of some of the enzymes involved in anaerobic metabolism in the brain, and this may be an adaptation to hypoxia.

Determining environmental causes of biological effects: the need for a mechanistic physiological dimension in conservation biology

The emerging field of Conservation Physiology links environmental change and ecological success by the application of physiological theory, approaches and tools to elucidate and address conservation problems. Human activity has changed the natural environment to a point where the viability of many ecosystems is now under threat. There are already many descriptions of how changes in biological patterns are correlated with environmental changes. The next important step is to determine the causative relationship between environmental variability and biological systems. Physiology provides the mechanistic link between environmental change and ecological patterns. Physiological research, therefore, should be integrated into conservation to predict the biological consequences of human activity, and to identify those species or populations that are most vulnerable.

The effect of short-term hypoxic exposure on metabolic gene expression

The long-term effect of hypoxia is to decrease both the production and use of ATP and thus decrease the reliance on mitochondrial oxidative energy production. Yet, recent studies include more immediate affects of hypoxia on gene expression and these data suggest the maintenance of mitochondrial function. To better understand the short-term physiological response to hypoxia, we quantified metabolic mRNA expression in the heart ventricles and livers of the teleost fish Fundulus grandis exposed to partial oxygen pressure of 2.8 kPa (-13.5% air saturation).Twenty-eight individuals from a single population were exposed to hypoxia for 0, 4, 8, 12, 24, 48, and 96  hr. Liver and cardiac tissues were sampled from the same individuals at 0-48  hr. At 96  hr, only cardiac tissue was assayed. Gene expression was significantly different (ANOVA, P < 0.05) for 17 of 226 metabolic genes (7.5%) in cardiac tissue and for 20 of 256 (7.8%) metabolic genes in hepatic tissue. For the two tissues examined in this study, the maximum response occurred at different times. For cardiac tissue, using Dunnett's post hoc test, most of these significant differences occurred at 96  hr of exposure. For liver, all but one significant difference occurred at 4  hr. Surprisingly, too many (relative to random expectations) of the genes with significant increase in mRNA are involved in the oxidative phosphorylation pathway: 44% of the significant genes at 96  hr in the heart and 33% of the significant genes at 4  hr in the liver are involved in the oxidative phosphorylation pathway. These data indicate that there are tissue-specific differences in the timing of the response to hypoxia, yet both cardiac and hepatic tissues have increases in mRNA that code for enzyme in the oxidative phosphorylation pathway. If these changes in mRNA produce a similar change in protein, then these data suggest that the initial response to hypoxia involves an increase in the oxidative pathway potentially as a mechanism to maintain ATP production.

Behavioural, morphological, and metabolic maturation of newly emerged adult workers of the bumblebee, Bombus impatiens

Newly emerged adult holometabolous insects must still complete considerable morphological, metabolic, and neural maturation. Despite this, adults have frequently been documented to fly prior to finishing maturation and attaining peak physiological capacity. In some species, flight is limited by the unfurling of the wing, while in other species it may be limited by biochemical capacity. We charted maturation trajectories of adult bumblebee workers (Bombus impatiens) for both morphological and flight muscle metabolic capacities, and compared these to the first age at flight. Workers began regular flights as soon as two days after emergence. The unfurling of the wings was completed well before first flights and before any other studied factor, suggesting this did not initially limit flight. Wing beat frequencies, measured as a struggling response to grasping the hindlegs, were about 90% mature by two days old, and did not significantly change after three days. Conversely, by the initiation of flight, the mean enzyme maturation was only 63% completed relative to adult enzyme capacity, though specific enzyme profiles ranged from 42% to 73%. Maximum ADP-stimulated mitochondrial respiratory activity on pyruvate and proline matured along a similar time frame to glycolytic capacity, reaching its maximum three days after emergence. Bumblebees, as other adult insects, thus begin flights prior to fully maturing.

Switch between life history strategies due to changes in glycolytic enzyme gene dosage in Saccharomyces cerevisiae

Adaptation is the process whereby a population or species becomes better fitted to its habitat through modifications of various life history traits which can be positively or negatively correlated. The molecular factors underlying these covariations remain to be elucidated. Using Saccharomyces cerevisiae as a model system, we have investigated the effects on life history traits of varying the dosage of genes involved in the transformation of resources into energy. Changing gene dosage for each of three glycolytic enzyme genes (hexokinase 2, phosphoglucose isomerase, and fructose-1,6-bisphosphate aldolase) resulted in variation in enzyme activities, glucose consumption rate, and life history traits (growth rate, carrying capacity, and cell size). However, the range of effects depended on which enzyme was expressed differently. Most interestingly, these changes revealed a genetic trade-off between carrying capacity and cell size, supporting the discovery of two extreme life history strategies already described in yeast populations: the "ants," which have lower glycolytic gene dosage, take up glucose slowly, and have a small cell size but reach a high carrying capacity, and the "grasshoppers," which have higher glycolytic gene dosage, consume glucose more rapidly, and allocate it to a larger cell size but reach a lower carrying capacity. These results demonstrate antagonist pleiotropy for glycolytic genes and show that altered dosage of a single gene drives a switch between two life history strategies in yeast.

Population proteomics: quantitative variation within and among populations in cardiac protein expression

Population analysis of gene expression is typically achieved by quantifying levels of mRNA; however, gene expression is also a function of protein translation and turnover. Therefore, a complete understanding of population variation in gene expression requires quantitative knowledge of protein expression within and among natural populations. We used two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) to quantitatively compare expression of heart ventricle proteins among 18 individuals in three populations of the teleost fish Fundulus. Among populations, expressions between orthologous proteins and mRNAs were generally positively correlated. Additionally, similar to the pattern of cardiac mRNA expression for the same populations, we found considerable variation in protein expression both within and among populations: Of 408 protein features in 2D gels, 34% are significantly different (P < 0.01) among individuals within a population, 9% differ between populations, and 12% have a pattern of expression that suggests they have evolved by natural selection. Although similar to mRNA expression, the frequency of significant differences among populations is larger for proteins. Similar to mRNA expressions, expressions of most proteins are correlated to the expressions of many other proteins. However, the correlations among proteins are more extensive than the correlation for similar RNAs. These correlations suggest a greater coordinate regulation of protein than mRNA expression. The larger frequency of significant differences among populations and the greater frequency of correlated expression among proteins versus among RNAs suggest that the molecular mechanisms affecting protein expression enhance the differences among populations, and these regulatory steps could be a source of variation for adaptation.

Genomic approaches with natural fish populations

Natural populations v. inbred stocks provide a much richer resource for identifying the effects of nucleotide substitutions because natural populations have greater polymorphism. Additionally, natural populations offer an advantage over most common research organisms because they are subject to natural selection, and analyses of these adaptations can be used to identify biologically important changes. Among fishes, these analyses are enhanced by having a wide diversity of species (>28 000 species, more than any other group of vertebrates) living in a huge range of environments (from below freezing to > 46 degrees C, in fresh water to salinities >40 ppt.). Moreover, fishes exhibit many different life-history and reproductive strategies and have many different phenotypes and social structures. Although fishes provide numerous advantages over other vertebrate models, there is still a dearth of available genomic tools for fishes. Fishes make up approximately half of all known vertebrate species, yet <0.2% of fish species have significant genomic resources. Nonetheless, genomic approaches with fishes have provided some of the first measures of individual variation in gene expression and insights into environmental and ecological adaptations. Thus, genomic approaches with natural fish populations have the potential to revolutionize fundamental studies of diverse fish species that offer myriad ecological and evolutionary questions.

Adaptation versus allometry: population and body mass effects on hypoxic metabolism in Fundulus grandis

Hypoxia has significant effects on organisms, from metabolic reduction to death, and could be an important evolutionary force affecting the variation among populations within a species. To determine intraspecific variation in hypoxic metabolism and the effect of body mass, we examine rates of oxygen consumption (M(O2)) at seven oxygen concentrations among seven populations of Fundulus grandis that inhabit a mosaic of habitats with different frequencies and intensities of hypoxia. For M(O2), there is a significant interaction (P< 0.05) between body mass and oxygen concentrations: log(10) body mass: log(10) M(O2) slopes were steeper at intermediate oxygen partial pressures (Po(2)) than either normoxic or lowest Po(2) (ANCOVA, P<0.001). Additionally, the PO(2crit) (Po(2) where M(O2) can no longer be maintained) was a negative function of body mass (P < 0.04). At the lowest Po(2) (1.8 kPa), there was a significant difference in M(O2) among populations: one of the populations from environments more frequently stressed by hypoxia has greater M(O2) at the lowest oxygen concentrations. With few differences among populations, the most important effects were how body mass affected M(O2) at intermediate Po(2) and the negative relationship between body mass and PO(2crit). These findings suggest that an increase in body size is a useful strategy to minimize the effect of hypoxia.

Directional evolution of the slope of the metabolic rate-temperature relationship is correlated with climate

Abstract The evolution of metabolic rate-temperature (MR-T) reaction norms is of fundamental importance to physiological ecology. Metabolic cold adaptation (MCA) predicts that populations or species from cooler environments will have either a higher metabolic rate at a common temperature or steeper MR-T relationships, indicating greater sensitivity of respiratory metabolism to temperature. Support for MCA has been found in some insect species by comparing species or populations differing in latitude. However, the generality of these findings are contentious, with most studies either unable to account for phenotypic plasticity or the evolutionary relatedness of species or populations. Hence, the importance of MCA is vigorously debated from both evolutionary and ecological perspectives. Furthermore, few species, particularly from tropical environments, have been shown to differ in MR-T sensitivity along altitudinal temperature gradients. Here, using four populations of tsetse flies (Glossina pallidipes, Diptera: Glossinidae) from thermally distinct geographic regions, we test the hypothesis that there is evolved variation in MR-T relationships to cold climates. We found that a high-altitude equatorial population from a cool habitat has a steeper MR-T reaction norm. By contrast, other populations from warmer environments in East Africa do not differ with respect to their MR-T reaction norms. Squared-change parsimony analyses, based on the combined mitochondrial 16S rDNA ribosomal subunit and cytochrome c oxidase subunit I (COI), support the hypothesis of adaptive differentiation of MR-T reaction norms in the cool-climate population. Seasonal adjustments or laboratory-temperature-induced phenotypic plasticity changed the intercept of the reaction norm rather than the slope, and thus the observed intraspecific variation in slopes of MR-T reaction norms could not be accounted for by phenotypic plasticity. These results therefore suggest evolutionary adaptation of MR-T reaction norms to cool climates (<22 degrees C) in tsetse and provide novel support for MCA within an insect species.

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Comparative characterization of a temperature responsive gene (lactate dehydrogenase-B, ldh-b) in two congeneric tropical fish, Lates calcarifer and Lates niloticus

The characterization of candidate loci is a critical step in obtaining insight into adaptation and acclimation of organisms. In this study of two non-model tropical (to sub-tropical) congeneric perciformes (Lates calcarifer and Lates niloticus) we characterized both coding and non-coding regions of lactate dehydrogenase-B (ldh-b), a locus which exhibits temperature-adaptive differences among temperate and sub-tropical populations of the North American killifish Fundulus heteroclitus. Ldh-b was 5,004 and 3,527 bp in length in L. calcarifer and L. niloticus, respectively, with coding regions comprising 1,005 bp in both species. A high level of sequence homology existed between species for both coding and non-coding regions of ldh-b (> 97% homology), corresponding to a 98.5% amino acid sequence homology. All six known functional sites within the encoded protein sequence (LDH-B) were conserved between the two Lates species. Ten simple sequence repeat (SSR) motifs (mono-, di-, tri- and tetranucleotide) and thirty putative microRNA elements (miRNAs) were identified within introns 1, 2, 5 and 6 of both Lates species. Five single nucleotide polymorphisms (SNPs) were also identified within miRNA containing intron regions. Such SNPs are implicated in several complex human conditions and/or diseases (as demonstrated by extensive genome-wide association studies). This novel characterization serves as a platform to further examine how non-model species may respond to changes in their native temperatures, which are expected to increase by up to 6 degrees C over the next century.

Proteomic changes in the crucian carp brain during exposure to anoxia

During exposure to anoxia, the crucian carp brain is able to maintain normal overall protein synthesis rates. However, it is not known if there are alterations in the synthesis or expression of specific proteins. This investigation addresses this issue by comparing the normoxic and anoxic brain proteome. Nine proteins were found to be reduced by anoxia. Reductions in the glycolytic pathway proteins creatine kinase, fructose biphosphate aldolase, glyceraldehyde-3-phosphate dehydrogenase, triosephosphate isomerase and lactate dehydrogenase reflect the reduced production and requirement for adenosine tri-phosphate during anoxia. In terms of neural protection, voltage-dependent anion channel, a protein associated with neuronal apoptosis, was reduced, along with gefiltin, a protein associated with the subsequent need for neuronal repair. Additionally the expression of proteins associated with neural degeneration and impaired cognitive function also declined; dihydropyrimidinase-like protein-3 and vesicle amine transport protein-1. One protein was found to be increased by anoxia; pre-proependymin, the precursor to ependymin. Ependymin fulfils multiple roles in neural plasticity, memory formation and learning, neuron growth and regeneration, and is able to reverse the possibility of apoptosis, thus further protecting the anoxic brain.

Swimming performance and energetics as a function of temperature in killifish Fundulus heteroclitus

Populations of the common killifish Fundulus heteroclitus are found along a latitudinal temperature gradient in habitats with high thermal variability. The objectives of this study were to assess the effects of temperature and population of origin on killifish swimming performance (assessed as critical swimming speed, U(crit)). Acclimated fish from northern and southern killifish populations demonstrated a wide zone (from 7 degrees to 33 degrees C) over which U(crit) showed little change with temperature, with performance declining significantly only at lower temperatures. Although we observed significant differences in swimming performance between a northern and a southern population of killifish in one experiment, with northern fish having an approximately 1.5-fold-greater U(crit) than southern fish across all acclimation temperatures, we were unable to replicate this finding in other populations or collection years, and performance was consistently high across all populations and at both low (7 degrees C) and high (23 degrees C) acclimation temperatures. The poor swimming performance of southern killifish from a single collection year was correlated with low muscle [glycogen] rather than with other indicators of fuel stores or body condition. Killifish acclimated to 18 degrees C and acutely challenged at temperatures of 5 degrees , 18 degrees , 25 degrees , or 34 degrees C showed modest thermal sensitivity of U(crit) between 18 degrees and 34 degrees C, with performance declining substantially at 5 degrees C. Thus, much of the zone of relative thermal insensitivity of swimming performance is intrinsic in this species rather than acquired as a result of acclimation. These data suggest that killifish are broadly tolerant of changing temperatures, whether acute or chronic, and demonstrate little evidence of local adaptation in endurance swimming performance in populations from different thermal habitats.

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus. This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark-recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Marine macrophysiology: Studying physiological variation across large spatial scales in marine systems

A new approach toward understanding marine ecosystems has emerged through the integration of ecological physiology and macroecology. This multidisciplinary approach, titled here marine macrophysiology, facilitates unique insight into the foundation of macro-scale ecological patterns, such as biogeographic distributions, via examination of functional attributes of marine organisms across large spatial scales. For example, these broad-scale physiological inquiries confer the ability to directly assess the abundant-center hypothesis (aka Brown's principle) which proposes that species have decreased performance toward their range edges. By extension, the marine macrophysiological perspective also stands to clarify our understanding of more complex macro-scale phenomena such as biological invasions, the design of marine protected areas, and species' responses to global climate change. In this article, we review recent marine macrophysiology research and offer insights into future directions for this emerging field.

The biological importance of measuring individual variation

Functional genomics research using Fundulus heteroclitus has focused on variation among individuals because of the evolutionary importance and value of Fundulus in explaining the human condition (why individual humans are different and are affected differently by stress,disease and drugs). Among different populations and species of Fundulus, there are evolutionarily adaptive differences in gene expression. This natural variation in gene expression seems to affect cardiac metabolism because up to 81% of the variation in glucose utilization observed in isolated heart ventricles is related to specific patterns of gene expression. The surprising result from this research is that among different groups of individuals, the expression of mRNA from different metabolic pathways explains substrate-specific metabolism. For example, variation in oxidative phosphorylation mRNAs explains glucose metabolism for one group of individuals but expression of glucose metabolism genes explains this metabolism in a different group of individuals. This variation among individuals has important implications for studies using inbred strains:conclusions based on one individual or one strain will not necessarily reflect a generalized conclusion for a population or species. Finally, there are surprisingly strong positive and negative correlations among metabolic genes,both within and between pathways. These data suggest that measures of mRNA expression are meaningful, yet there is a complexity in how gene expression is related to physiological processes.

A falsification of the thermal specialization paradigm: compensation for elevated temperatures in Antarctic fishes

Specialization to a particular environment is one of the main factors used to explain species distributions. Antarctic fishes are often cited as a classic example to illustrate the specialization process and are regarded as the archetypal stenotherms. Here we show that the Antarctic fish Pagothenia borchgrevinki has retained the capacity to compensate for chronic temperature change. By displaying astounding plasticity in cardiovascular response and metabolic control, the fishes maintained locomotory performance at elevated temperatures. Our falsification of the specialization paradigm indicates that the effect of climate change on species distribution and extinction may be overestimated by current models of global warming.

Lifetime performance in foraging honeybees: behaviour and physiology

Honeybees, Apis mellifera, gradually increase their rate of forage uptake as they gain foraging experience. This increase in foraging performance has been proposed to occur as a result of learning; however, factors affecting flight ability such as changes in physiological components of flight metabolism could also contribute to this pattern. Thus, the purpose of this study was to assess the contribution of physiological changes to the increase in honeybee foraging performance. We investigated aspects of honeybee flight muscle biochemistry throughout the adult life, from non-foraging hive bees, through young and mature foragers, to old foragers near the end of their lifespan. Two-dimensional gel proteomic analysis on honeybee thorax muscle revealed an increase in several proteins from hive bees to mature foragers including troponin T 10a, aldolase and superoxide dismutase. By contrast, the activities (V(max)) of enzymes involved in aerobic performance, phosphofructokinase, hexokinase, pyruvate kinase and cytochrome c oxidase, did not increase in the flight muscles of hive bees, young foragers, mature foragers and old foragers. However, citrate synthase activity was found to increase with foraging experience. Hence, our results suggest plasticity in both structural and metabolic components of flight muscles with foraging experience.

Variation within and among species in gene expression: raw material for evolution

Heritable variation in regulatory or coding regions is the raw material for evolutionary processes. The advent of microarrays has recently promoted examination of the extent of variation in gene expression within and among taxa and examination of the evolutionary processes affecting variation. This review examines these issues. We find: (i) microarray-based measures of gene expression are precise given appropriate experimental design; (ii) there is large inter-individual variation, which is composed of a minor nongenetic component and a large heritable component; (iii) variation among populations and species appears to be affected primarily by neutral drift and stabilizing selection, and to a lesser degree by directional selection; and (iv) neutral evolutionary divergence in gene expression becomes nonlinear with greater divergence times due to functional constraint. Evolutionary analyses of gene expression reviewed here provide unique insights into partitioning of regulatory variation in nature. However, common limitations of these studies include the tendency to assume a linear relationship between expression divergence and species divergence, and failure to test explicit hypotheses that involve the ecological context of evolutionary divergence.

Neutral and adaptive variation in gene expression

Variation among populations in gene expression should be related to the accumulation of random-neutral changes and evolution by natural selection. The following evolutionary analysis has general applicability to biological and medical science because it accounts for genetic relatedness and identifies patterns of expression variation that are affected by natural selection. To identify genes evolving by natural selection, we allocate the maximum among-population variation to genetic distance and then examine the remaining variation relative to a hypothesized important ecological parameter (temperature). These analyses measure the expression of metabolic genes in common-gardened populations of the fish Fundulus heteroclitus whose habitat is distributed along a steep thermal gradient. Although much of the variation in gene expression fits a null model of neutral drift, the variation in expression for 22% of the genes that regress with habitat temperature was far greater than could be accounted for by genetic distance alone. The most parsimonious explanation for among-population variation for these genes is evolution by natural selection. In addition, many metabolic genes have patterns of variation incongruent with neutral evolution: They have too much or too little variation. These patterns of biological variation in expression may reflect important physiological or ecological functions.

Evolution of enzymes in a series is driven by dissimilar functional demands

That distinct enzyme activities in an unbranched metabolic pathway are evolutionarily tuned to a single functional requirement is a pervasive assumption. Here we test this assumption by examining the activities of two consecutively acting enzymes in human erythrocytes with an approach to quantitative evolutionary design that avoids the above-mentioned assumption. We previously found that avoidance of NADPH depletion during the pulses of oxidative load to which erythrocytes are normally exposed is the main functional requirement mediating selection for high glucose-6-phosphate dehydrogenase activity. In the present study, we find that, in contrast, the maintenance of oxidized glutathione at low concentrations is the main functional requirement mediating selection for high glutathione reductase activity. The results in this case show that, contrary to the assumption of a single functional requirement, natural selection for the normal activities of the distinct enzymes in the pathway is mediated by different requirements. On the other hand, the results agree with the more general principles that underlie our approach. Namely, that (i) the values of biochemical parameters evolve so as to fulfill the various performance requirements that are relevant to achieve high fitness, and (ii) these performance requirements can be inferred from quantitative systems theory considerations, informed by knowledge of specific aspects of the biochemistry, physiology, genetics, and ecology of the organism.

Steady-state effects of temperature acclimation on the transcriptome of the rainbow trout heart

Cold-acclimated (CA) phenotype of trout heart was induced by 4-wk acclimation at 4°C and was characterized by 32.7% increase in relative heart mass and 49.8% increase in ventricular myocyte size compared with warm-acclimated (WA; 18°C) fish ( P < 0.001). Effect of temperature acclimation on transcriptome of the rainbow trout heart was examined using species-specific microarray chips containing 1,380 genes. After 4 wk of temperature acclimation, 8.8% (122) of the genes were differently expressed in CA and WA hearts, and most of them (82%) were upregulated in the cold ( P < 0.01). Transcripts of genes engaged in protein synthesis and intermediary metabolism were most strongly upregulated, whereas genes contributing to the connective tissue matrix were clearly repressed. Extensive upregulation of the genes coding for ribosomal proteins and translation elongation and initiation factors suggest that the protein synthesis machinery of the trout heart is enhanced in the cold and is an essential part of the compensatory mechanism causing and maintaining the hypertrophy of cardiac myocytes. The prominent depression of collagen genes may be indicative of a reduced contribution of extracellular matrix to the remodeling of the CA fish heart. Temperature-related changes in transcripts of metabolic enzymes suggest that at mRNA level, glycolytic energy production from carbohydrates is compensated in the heart of CA rainbow trout, while metabolic compensation is absent in mitochondria. In addition, the analysis revealed three candidate genes: muscle LIM protein, atrial natriuretic peptide B, and myosin light chain 2, which might be central for induction and maintenance of the hypertrophic phenotype of the CA trout heart. These findings indicate that extensive modification of gene expression is needed to maintain the temperature-specific phenotype of the fish heart.

Fish as models for environmental genomics

Fish offer important advantages for defining the organism-environment interface and responses to natural or anthropogenic stressors. Genomic approaches using fish promise increased investigative power, and have already provided insights into the mechanisms that underlie short-term and long-term environmental adaptations. The range of fish species for which genomic resources are available is increasing, but will require significant further expansion for the optimal application of fish environmental genomics.

FunnyBase: a systems level functional annotation of Fundulus ESTs for the analysis of gene expression

BACKGROUND

While studies of non-model organisms are critical for many research areas, such as evolution, development, and environmental biology, they present particular challenges for both experimental and computational genomic level research. Resources such as mass-produced microarrays and the computational tools linking these data to functional annotation at the system and pathway level are rarely available for non-model species. This type of "systems-level" analysis is critical to the understanding of patterns of gene expression that underlie biological processes.

RESULTS

We describe a bioinformatics pipeline known as FunnyBase that has been used to store, annotate, and analyze 40,363 expressed sequence tags (ESTs) from the heart and liver of the fish, Fundulus heteroclitus. Primary annotations based on sequence similarity are linked to networks of systematic annotation in Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) and can be queried and computationally utilized in downstream analyses. Steps are taken to ensure that the annotation is self-consistent and that the structure of GO is used to identify higher level functions that may not be annotated directly. An integrated framework for cDNA library production, sequencing, quality control, expression data generation, and systems-level analysis is presented and utilized. In a case study, a set of genes, that had statistically significant regression between gene expression levels and environmental temperature along the Atlantic Coast, shows a statistically significant (P < 0.001) enrichment in genes associated with amine metabolism.

CONCLUSION

The methods described have application for functional genomics studies, particularly among non-model organisms. The web interface for FunnyBase can be accessed at http://genomics.rsmas.miami.edu/funnybase/super_craw4/. Data and source code are available by request at jpaschall@bioinfobase.umkc.edu.

Natural variation in cardiac metabolism and gene expression in Fundulus heteroclitus

Individual variation in gene expression is important for evolutionary adaptation and susceptibility to diseases and pathologies. In this study, we address the functional importance of this variation by comparing cardiac metabolism to patterns of mRNA expression using microarrays. There is extensive variation in both cardiac metabolism and the expression of metabolic genes among individuals of the teleost fish Fundulus heteroclitus from natural outbred populations raised in a common environment: metabolism differed among individuals by a factor of more than 2, and expression levels of 94% of genes were significantly different (P < 0.01) between individuals in a population. This unexpectedly high variation in metabolic gene expression explains much of the variation in metabolism, suggesting that it is biologically relevant. The patterns of gene expression that are most important in explaining cardiac metabolism differ between groups of individuals. Apparently, the variation in metabolism seems to be related to different patterns of gene expression in the different groups of individuals. The magnitude of differences in gene expression in these groups is not important; large changes in expression have no greater predictive value than small changes. These data suggest that variation in physiological performance is related to the subtle variation in gene expression and that this relationship differs among individuals.

Changes in gene expression as biochemical adaptations to environmental change: a tribute to Peter Hochachka

Changes in gene expression are likely to play a critical role in both acclimation and adaptation to a changing environment. There is a rapidly growing body of literature implicating quantitative changes in gene expression during acclimation to environmental change, but less is known about the role of qualitative changes in gene expression, such as switching between alternative isoforms. Alternative isoforms can arise via gene duplication, alternative splicing, or alternative promoter usage. Organisms that have undergone recent genome duplication events may make use of environment-specific isoforms coded by multiple genes, but their role in other organisms is less well known. However, recent data suggest that isoforms arising from alternative splicing may be an under-appreciated source of physiological variation. The role of changes in gene expression during evolutionary adaptation has received comparatively limited attention, but novel approaches to addressing the adaptive significance of changes in gene expression have been applied to a few cases of differences in gene expression among taxa. Recent advances in genomics, including microarray technology, knock-out and knock-down approaches, and the wealth of data coming from large-scale sequencing projects have provided (and will continue to provide at ever increasing rates) new insights into these classic questions in comparative biochemistry.

Stress and interpopulation variation in glycolytic enzyme activity and expression in a teleost fish Fundulus heteroclitus

Northern populations of Fundulus heteroclitus have twofold greater activity of lactate dehydrogenase-B (LDH-B) than southern populations, but exposure to stress increases LDH-B in southern populations, abolishing this difference. To test whether differences in the activity of other hepatic glycolytic enzymes between populations are sensitive to stress, we injected fish with a pharmacological dose of cortisol in coconut oil (400 microg g(-1)) or exposed them to handling stress and measured the activities of all the glycolytic enzymes. At rest, liver phosphofructokinase (PFK) and aldolase (ALD) activities were greater in southern fish, whereas LDH-B activity was greater in northern fish. No other glycolytic enzymes differed in activity between populations in control fish. Cortisol injection and handling stress decreased PFK and ALD and increased LDH activities in the southern but not the northern population, such that the populations no longer differed in the activity of any enzyme following treatment. Unlike Ldh-B mRNA, Pfk and Ald mRNA levels did not parallel enzyme activity, suggesting complex kinetics or regulation at multiple levels. Plasma cortisol did not differ between populations at rest but was significantly different between populations in treated fish. These data suggest that differences in liver enzyme activity may be related to differences in stress hormone physiology between populations.