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

Estimating maximum oxygen uptake of fishes during swimming and following exhaustive chase - different results, biological bases and applications

The maximum rate at which animals take up oxygen from their environment (ṀO2,max) is a crucial aspect of their physiology and ecology. In fishes, ṀO2,max is commonly quantified by measuring oxygen uptake either during incremental swimming tests or during recovery from an exhaustive chase. In this Commentary, we compile recent studies that apply both techniques to the same fish and show that the two methods typically yield different mean estimates of ṀO2,max for a group of individuals. Furthermore, within a group of fish, estimates of ṀO2,max determined during swimming are poorly correlated with estimates determined during recovery from chasing (i.e. an individual's ṀO2,max is not repeatable across methods). One explanation for the lack of agreement is that these methods measure different physiological states, each with their own behavioural, anatomical and biochemical determinants. We propose that these methods are not directly interchangeable but, rather, each is suited to address different questions in fish biology. We suggest that researchers select the method that reflects the biological contexts of their study, and we advocate for the use of accurate terminology that acknowledges the technique used to elevate ṀO2 (e.g. peak ṀO2,swim or peak ṀO2,recovery). If the study's objective is to estimate the 'true' ṀO2,max of an individual or species, we recommend that pilot studies compare methods, preferably using repeated-measures designs. We hope that these recommendations contribute new insights into the causes and consequences of variation in ṀO2,max within and among fish species.

Method dependency of maximum oxygen uptake rate and its repeatability in the Gulf killifish, Fundulus grandis

The maximum rate at which fish can take up oxygen from their environment to fuel aerobic metabolism is an important feature of their physiology and ecology. Methods to quantify maximum oxygen uptake rate (ṀO2), therefore, should reliably and reproducibly estimate the highest possible ṀO2 by an individual or species under a given set of conditions (peak ṀO2). This study determined peak ṀO2 and its repeatability in Gulf killifish, Fundulus grandis, subjected to three methods to elevate metabolism: swimming at increasing water speeds, during recovery after an exhaustive chase, and after ingestion of a large meal. Estimates of peak ṀO2 during swimming and after an exhaustive chase were repeatable across two trials, whereas peak ṀO2 after feeding was not. Peak ṀO2 determined by the three methods was significantly different from one another, being highest during swimming, lowest after an exhaustive chase, and intermediate after feeding. In addition, peak ṀO2 during recovery from an exhaustive chase depended on the length of time of recovery: in nearly 60% of the trials, values within the first hour of the chase were lower than those measured later. A novel and important finding was that an individual's peak ṀO2 was not repeatable when compared across methods. Therefore, the peak ṀO2 estimated for a group of fish, as well as the ranking of individual ṀO2 within that group, depends on the method used to elevate aerobic metabolism.

Interindividual variation in maximum aerobic metabolism varies with gill morphology and myocardial bioenergetics

This study asked whether interindividual variation in maximum and standard aerobic metabolic rates of the Gulf killifish, Fundulus grandis, correlate with gill morphology and cardiac mitochondrial bioenergetics, traits reflecting critical steps in the O2 transport cascade from the environment to the tissues. Maximum metabolic rate (MMR) was positively related to body mass, total gill filament length, and myocardial oxygen consumption during maximum oxidative phosphorylation (multiple R2=0.836). Standard metabolic rate (SMR) was positively related to body mass, total gill filament length, and myocardial oxygen consumption during maximum electron transport system activity (multiple R2=0.717). After controlling for body mass, individuals with longer gill filaments, summed over all gill arches, or greater cardiac respiratory capacity had higher whole-animal metabolic rates. The overall model fit and the explanatory power of individual predictor variables were better for MMR than for SMR, suggesting that gill morphology and myocardial bioenergetics are more important in determining active rather than resting metabolism. After accounting for body mass, heart ventricle mass was not related to variation in MMR or SMR, indicating that the quality of the heart (i.e., the capacity for mitochondrial metabolism) was more influential than heart size. Finally, the myocardial oxygen consumption required to offset the dissipation of the transmembrane proton gradient in the absence of ATP synthesis was not correlated with either MMR or SMR. The results support the idea that interindividual variation in aerobic metabolism, particularly maximum metabolic rate, is associated with variation in specific steps in the O2 transport cascade.

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

Physiological and behavioural strategies of aquatic animals living in fluctuating environments

Shallow or near-shore environments, such as ponds, estuaries and intertidal zones, are among the most physiologically challenging of all aquatic settings. Animals inhabiting these environments experience conditions that fluctuate markedly over relatively short temporal and spatial scales. Living in these habitats requires the ability to tolerate the physiological disturbances incurred by these environmental fluctuations. This tolerance is achieved through a suite of physiological and behavioural responses that allow animals to maintain homeostasis, including the ability to dynamically modulate their physiology through reversible phenotypic plasticity. However, maintaining the plasticity to adjust to some stresses in a dynamic environment may trade off with the capacity to deal with other stressors. This paper will explore studies on select fishes and invertebrates exposed to fluctuations in dissolved oxygen, salinity and pH. We assess the physiological mechanisms these species employ to achieve homeostasis, with a focus on the plasticity of their responses, and consider the resulting physiological trade-offs in function. Finally, we discuss additional factors that may influence organismal responses to fluctuating environments, such as the presence of multiple stressors, including parasites. We echo recent calls from experimental biologists to consider physiological responses to life in naturally fluctuating environments, not only because they are interesting in their own right but also because they can reveal mechanisms that may be crucial for living with increasing environmental instability as a consequence of climate change.

Interindividual plasticity in metabolic and thermal tolerance traits from populations subjected to recent anthropogenic heating

To better understand temperature's role in the interaction between local evolutionary adaptation and physiological plasticity, we investigated acclimation effects on metabolic performance and thermal tolerance among natural Fundulus heteroclitus (small estuarine fish) populations from different thermal environments. Fundulus heteroclitus populations experience large daily and seasonal temperature variations, as well as local mean temperature differences across their large geographical cline. In this study, we use three populations: one locally heated (32°C) by thermal effluence (TE) from the Oyster Creek Nuclear Generating Station, NJ, and two nearby reference populations that do not experience local heating (28°C). After acclimation to 12 or 28°C, we quantified whole-animal metabolic (WAM) rate, critical thermal maximum (CTmax) and substrate-specific cardiac metabolic rate (CaM, substrates: glucose, fatty acids, lactate plus ketones plus ethanol, and endogenous (i.e. no added substrates)) in approximately 160 individuals from these three populations. Populations showed few significant differences due to large interindividual variation within populations. In general, for WAM and CTmax, the interindividual variation in acclimation response (log2 ratio 28/12°C) was a function of performance at 12°C and order of acclimation (12-28°C versus 28-12°C). CTmax and WAM were greater at 28°C than 12°C, although WAM had a small change (2.32-fold) compared with the expectation for a 16°C increase in temperature (expect 3- to 4.4-fold). By contrast, for CaM, the rates when acclimatized and assayed at 12 or 28°C were nearly identical. The small differences in CaM between 12 and 28°C temperature were partially explained by cardiac remodeling where individuals acclimatized to 12°C had larger hearts than individuals acclimatized to 28°C. Correlation among physiological traits was dependent on acclimation temperature. For example, WAM was negatively correlated with CTmax at 12°C but positively correlated at 28°C. Additionally, glucose substrate supported higher CaM than fatty acid, and fatty acid supported higher CaM than lactate, ketones and alcohol (LKA) or endogenous. However, these responses were highly variable with some individuals using much more FA than glucose. These findings suggest interindividual variation in physiological responses to temperature acclimation and indicate that additional research investigating interindividual may be relevant for global climate change responses in many species.

Oxygen supply capacity breathes new life into critical oxygen partial pressure (Pcrit)

The critical oxygen partial pressure (Pcrit), typically defined as the PO2 below which an animal's metabolic rate (MR) is unsustainable, is widely interpreted as a measure of hypoxia tolerance. Here, Pcrit is defined as the PO2 at which physiological oxygen supply (α0) reaches its maximum capacity (α; µmol O2 g-1 h-1 kPa-1). α is a species- and temperature-specific constant describing the oxygen dependency of the maximum metabolic rate (MMR=PO2×α) or, equivalently, the MR dependence of Pcrit (Pcrit=MR/α). We describe the α-method, in which the MR is monitored as oxygen declines and, for each measurement period, is divided by the corresponding PO2 to provide the concurrent oxygen supply (α0=MR/PO2). The highest α0 value (or, more conservatively, the mean of the three highest values) is designated as α. The same value of α is reached at Pcrit for any MR regardless of previous or subsequent metabolic activity. The MR need not be constant (regulated), standardized or exhibit a clear breakpoint at Pcrit for accurate determination of α. The α-method has several advantages over Pcrit determination and non-linear analyses, including: (1) less ambiguity and greater accuracy, (2) fewer constraints in respirometry methodology and analysis, and (3) greater predictive power and ecological and physiological insight. Across the species evaluated here, α values are correlated with MR, but not Pcrit. Rather than an index of hypoxia tolerance, Pcrit is a reflection of α, which evolves to support maximum energy demands and aerobic scope at the prevailing temperature and oxygen level.