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

Hypoxia Tolerance of Two Killifish Species

Hypoxia tolerance in aquatic ectotherms involves a suite of behavioral and physiological responses at the organismal, tissue, and cellular levels. The current study evaluated two closely related killifish species (Fundulus heteroclitus, Fundulus majalis) to evaluate responses to acute moderate and acute severe hypoxia. Routine metabolic rate and loss of equilibrium were assessed, followed by analysis in skeletal muscle of markers of oxidative damage to proteins (2,4-DNPH), lipids (4-HNE), and DNA (8-OHdG), hypoxia signaling (HIF1α, HIF2α), cellular energy state (p-AMPK: AMPK), and protein degradation (Ubiquitin, LC3B, Calpain 2, Hsp70). Both species had a similar reduction in metabolic rate at low PO2. However, F. heteroclitus was the more hypoxia-tolerant species based on a lower PO2 at which there was loss of equilibrium, perhaps due in part to a lower oxygen demand at all oxygen tensions. Despite the differences in hypoxia tolerance between the species, skeletal muscle molecular markers were largely insensitive to hypoxia, and there were few differences in responses between the species. Thus, the metabolic depression observed at the whole animal level appears to limit perturbations in skeletal muscle in both species during the hypoxia treatments.

Interactive effects of intrinsic and extrinsic factors on metabolic rate

Metabolism energizes all biological processes, and its tempo may importantly influence the ecological success and evolutionary fitness of organisms. Therefore, understanding the broad variation in metabolic rate that exists across the living world is a fundamental challenge in biology. To further the development of a more reliable and holistic picture of the causes of this variation, we review several examples of how various intrinsic (biological) and extrinsic (environmental) factors (including body size, cell size, activity level, temperature, predation and other diverse genetic, cellular, morphological, physiological, behavioural and ecological influences) can interactively affect metabolic rate in synergistic or antagonistic ways. Most of the interactive effects that have been documented involve body size, temperature or both, but future research may reveal additional 'hub factors'. Our review highlights the complex, intimate inter-relationships between physiology and ecology, knowledge of which can shed light on various problems in both disciplines, including variation in physiological adaptations, life histories, ecological niches and various organism-environment interactions in ecosystems. We also discuss theoretical and practical implications of interactive effects on metabolic rate and provide suggestions for future research, including holistic system analyses at various hierarchical levels of organization that focus on interactive proximate (functional) and ultimate (evolutionary) causal networks. This article is part of the theme issue 'The evolutionary significance of variation in metabolic rates'.

Plasticity, repeatability and phenotypic correlations of aerobic metabolic traits in a small estuarine fish

Standard metabolic rate (SMR), maximum metabolic rate (MMR), absolute aerobic scope (AAS) and critical oxygen tension (P_crit) were determined for the Gulf killifish, Fundulus grandis, an ecologically dominant estuarine fish, acclimated to lowered salinity, elevated temperature and lowered oxygen concentration. Acclimation to low salinity resulted in a small, but significant, elevation of P_crit (suggesting lower tolerance of hypoxia); acclimation to elevated temperature increased SMR, MMR, AAS and P_crit; acclimation to low oxygen led to a small increase in SMR, but substantial decreases in MMR, AAS and P_crit Variation in these metabolic traits among individuals was consistent and repeatable when measured during multiple control exposures over 7 months. Trait repeatability was unaffected by acclimation condition, suggesting that repeatability of these traits is not context dependent. There were significant phenotypic correlations between specific metabolic traits: SMR was positively correlated with MMR and P_crit; MMR was positively correlated with AAS; and AAS was negatively correlated with P_crit In general, within-individual variation contributed more than among-individual variation to these phenotypic correlations. The effects of acclimation on these traits demonstrate that aerobic metabolism is plastic and influenced by the conditions experienced by these fish in the dynamic habitats in which they occur; however, the repeatability of these traits and the correlations among them suggest that these traits change in ways that maintain the rank order of performance among individuals across a range of environmental variation.

Plasticity, repeatability and phenotypic correlations of aerobic metabolic traits in a small estuarine fish

Standard metabolic rate (SMR), maximum metabolic rate (MMR), absolute aerobic scope (AAS) and critical oxygen tension (Pcrit) were determined for the Gulf killifish, Fundulus grandis, an ecologically dominant estuarine fish, acclimated to lowered salinity, elevated temperature and lowered oxygen concentration. Acclimation to low salinity resulted in a small, but significant, elevation of Pcrit (suggesting lower tolerance of hypoxia); acclimation to elevated temperature increased SMR, MMR, AAS and Pcrit; acclimation to low oxygen led to a small increase in SMR, but substantial decreases in MMR, AAS and Pcrit Variation in these metabolic traits among individuals was consistent and repeatable when measured during multiple control exposures over 7 months. Trait repeatability was unaffected by acclimation condition, suggesting that repeatability of these traits is not context dependent. There were significant phenotypic correlations between specific metabolic traits: SMR was positively correlated with MMR and Pcrit; MMR was positively correlated with AAS; and AAS was negatively correlated with Pcrit In general, within-individual variation contributed more than among-individual variation to these phenotypic correlations. The effects of acclimation on these traits demonstrate that aerobic metabolism is plastic and influenced by the conditions experienced by these fish in the dynamic habitats in which they occur; however, the repeatability of these traits and the correlations among them suggest that these traits change in ways that maintain the rank order of performance among individuals across a range of environmental variation.

Hypoxia and Sprint Swimming Performance of Juvenile Striped Bass, Morone saxatilis

Annual hypoxia in the Chesapeake Bay has expanded to the point where Darwinian fitness of juvenile striped bass (Morone saxatilis) may depend on their ability to perform in low-oxygen environments. The locomotion they use in predator/prey dynamics relies primarily on white (type II) muscle that is powered by anaerobic metabolic pathways and has generally been thought to be immune to aquatic hypoxia. We tested the sprint performance of 15 juvenile striped bass twice under acute hypoxia (20% air saturation [AS]) 5 wk apart and once under normoxia (>85% AS) in between. Average sprint performance was lower under the first hypoxia exposure than in normoxia and increased in the second hypoxia test relative to the first. The rank order of individual sprint performance was significantly repeatable when comparing the two hypoxia tests but not when compared with sprint performance measured under normoxic conditions. The maximum sprint performance of each individual was also significantly repeatable within a given day. Thus, sprint performance of striped bass is reduced under hypoxia, is phenotypically plastic, and improves with repetitive hypoxia exposures but is unrelated to relative sprint performance under normoxia. Since energy to fuel a sprint comes from existing ATP and creatine phosphate stores, the decline in sprint performance probably reflects reduced function of a part of the reflex chain leading from detection of aversive stimuli to activation of the muscle used to power the escape response.

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.

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.

Distinct physiological strategies are used to cope with constant hypoxia and intermittent hypoxia in killifish (Fundulus heteroclitus)

Many fish encounter hypoxia on a daily cycle, but the physiological effects of intermittent hypoxia are poorly understood. We investigated whether acclimation to constant (sustained) hypoxia or to intermittent diel cycles of nocturnal hypoxia (12 h normoxia:12 h hypoxia) had distinct effects on hypoxia tolerance or on several determinants of O2 transport and O2 utilization in estuarine killifish. Adult killifish were acclimated to normoxia, constant hypoxia, or intermittent hypoxia for 7 or 28 days in brackish water (4 ppt). Acclimation to both hypoxia patterns led to comparable reductions in critical O2 tension and resting O2 consumption rate, but only constant hypoxia reduced the O2 tension at loss of equilibrium. Constant (but not intermittent) hypoxia decreased filament length and the proportion of seawater-type mitochondrion-rich cells in the gills (which may reduce ion loss and the associated costs of active ion uptake), increased blood haemoglobin content, and reduced the abundance of oxidative fibres in the swimming muscle. In contrast, only intermittent hypoxia augmented the oxidative and gluconeogenic enzyme activities in the liver and increased the capillarity of glycolytic muscle, each of which should facilitate recovery between hypoxia bouts. Neither exposure pattern affected muscle myoglobin content or the activities of metabolic enzymes in the brain or heart, but intermittent hypoxia increased brain mass. We conclude that the pattern of hypoxia exposure has an important influence on the mechanisms of acclimation, and that the optimal strategies used to cope with intermittent hypoxia may be distinct from those for coping with constant hypoxia.

Relationship between fish size and metabolic rate in the oxyconforming inanga Galaxias maculatus reveals size-dependent strategies to withstand hypoxia

The relationship between metabolic rate and body size in animals is unlikely to be a constant but is instead shaped by a variety of intrinsic (i.e., physiological) and extrinsic (i.e., environmental) factors. This study examined the effect of environmental oxygen tension on oxygen consumption as a function of body mass in the galaxiid fish, inanga (Galaxias maculatus). As an oxyconformer, this fish lacks overt intrinsic regulation of oxygen consumption, eliminating this as a factor affecting the scaling relationship at different oxygen tensions. The relationship between oxygen consumption rate and body size was best described by a power function, with an exponent of 0.82, higher than the theoretical values of 0.66 or 0.75. The value of this exponent was significantly altered by environmental P(O2), first increasing as P(O2) decreased and then declining at the lowest P(O2) tested. These data suggest that the scaling exponent is species specific and regulated by extrinsic factors. Furthermore, the external P(O2) at which fish lost equilibrium was related to fish size, an effect explained by the scaling of anaerobic capacity with fish mass. Therefore, although bigger fish were forced to depress aerobic metabolism more rapidly than small fish when exposed to progressive hypoxia, they were better able to enact anaerobic metabolism, potentially extending their survival in hypoxia.

Low oxygen tolerance of different life stages of temperate freshwater fish species

Data on low dissolved oxygen (DO₂) tolerance of freshwater fish species of north-western Europe were used to create species sensitivity distributions (SSD). Lowest observed effect concentrations (LOEC) and 100% lethal concentrations (LC₁₀₀) data were collected from the scientific literature. Comparisons were made among life stages as well as between native and exotic species. In addition, lethal DO₂ concentrations were compared to oxygen concentrations corresponding to maximum tolerable water temperatures of the same species. Fish eggs and embryos were the least tolerant. Juveniles had a significantly lower mean LOEC than adults, but there was no difference in mean LC₁₀₀ between the two groups. The difference in lethal oxygen concentrations between adults and juveniles was largest for three salmonids, although it remains uncertain if this was a result of smoltification. There were no significant differences between native and exotic species; however, data on exotics are limited. DO₂ concentrations converted from maximum tolerable water temperatures were 3·9 times higher than the measured lethal DO₂ concentrations, which may reflect changes in respiration rates (Q₁₀) and may also relate to the simplicity of the model used.

The relationship between phenotypic and environmental variation: do physiological responses reduce interindividual differences?

What is the effect of a variable environment on phenotypic variation? Does the physiological response to a new environment increase or decrease the differences among individuals? We provide a speculative hypothesis suggesting that the induction of a physiological response to environmental change minimizes phenotypic differences among individuals in outbred genetically variable populations. Although this suggestion runs counter to the general idea that environmental variation induces phenotypic variation, we provide evidence that this is not always the case. One explanation for this counterintuitive hypothesis is that in a variable environment, the physiological mechanism that maintains homeostasis changes the concentrations of active transcription factors (TFs). This change in TFs reduces the effectiveness of nucleotide polymorphisms in TF binding sites and thus reduces the variation among individuals in mRNA expression and in the phenotypes affected by these mRNA transcripts. Thus, there are fewer differences among individuals in a variable environment compared with the variation observed in a constant environment. Our conjecture is that the physiological mechanisms that maintain homeostasis in response to environmental variation canalize phenotypic variation. If our hypothesis is correct, then the physiological canalization of gene expression in a variable environment hides genetic variation and thereby reduces the evolutionary costs of polymorphism. This hypothesis provides a new perspective on the mechanisms by which high levels of genetic variation can persist in real-world populations.

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.

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

Body mass dependence of glycogen stores in the anoxia-tolerant crucian carp (Carassius carassius L.)

Glycogen is a vital energy substrate for anaerobic organisms, and the size of glycogen stores can be a limiting factor for anoxia tolerance of animals. To this end, glycogen stores in 12 different tissues of the crucian carp (Carassius carassius L.), an anoxia-tolerant fish species, were examined. Glycogen content of different tissues was 2-10 times higher in winter (0.68-18.20% of tissue wet weight) than in summer (0.12-4.23%). In scale, bone and brain glycogen stores were strongly dependent on body mass (range between 0.6 and 785 g), small fish having significantly more glycogen than large fish (p < 0.05). In fin and skin, size dependence was evident in winter, but not in summer, while in other tissues (ventricle, atrium, intestine, liver, muscle, and spleen), no size dependence was found. The liver was much bigger in small than large fish (p < 0.001), and there was a prominent enlargement of the liver in winter irrespective of fish size. As a consequence, the whole body glycogen reserves, measured as a sum of glycogen from different tissues, varied from 6.1% of the body mass in the 1-g fish to 2.0% in the 800-g fish. Since anaerobic metabolic rate scales down with body size, the whole body glycogen reserves could provide energy for approximately 79 and 88 days of anoxia in small and large fish, respectively. There was, however, a drastic difference in tissue distribution of glycogen between large and small fish: in the small fish, the liver was the major glycogen store (68% of the stores), while in the large fish, the white myotomal muscle was the principal deposit of glycogen (57%). Since muscle glycogen is considered to be unavailable for blood glucose regulation, its usefulness in anoxia tolerance of the large crucian carp might be limited, although not excluded. Therefore, mobilization of muscle glycogen under anoxia needs to be rigorously tested.

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.