Oxygen limitation may affect the temperature and size dependence of metabolism in aquatic ectotherms
- aDivision of Biological Sciences, University of Montana, Missoula, MT 59812;
- bDepartment of Biology, McGill University, Montreal, QC H3A 1B1, Canada;
- cDepartment of Animal Ecology and Physiology, Institute for Water and Wetland Research, Radboud University Nijmegen, 6500 GL Nijmegen, the Netherlands;
- dDepartment of Environmental Science, Institute for Water and Wetland Research, Radboud University Nijmegen, 6500 GL Nijmegen, the Netherlands
See allHide authors and affiliations
Edited by James H. Brown, University of New Mexico, Morro Bay, CA, and approved October 22, 2020 (received for review February 20, 2020)

Significance
Organismal responses to climate change are mediated through its effects on physiology and metabolism. In aquatic environments, both water temperature and oxygen availability may modulate these responses by altering the aerobic metabolism fueling physiological performance. However, ecological models aimed at predicting how environmental factors shape aerobic metabolism disregard the role of oxygen supply. Here, we expand on these models by explicitly incorporating oxygen supply. Our results show that warmer water increases oxygen demand relative to supply, and the resulting reduction in aerobic scope appears to be stronger in larger individuals. Smaller aerobic scopes in warming water imply that climate change will reduce energy budgets needed to support the activities of aquatic animals and their physiological performance in the future.
Abstract
Both oxygen and temperature are fundamental factors determining metabolic performance, fitness, ecological niches, and responses of many aquatic organisms to climate change. Despite the importance of physical and physiological constraints on oxygen supply affecting aerobic metabolism of aquatic ectotherms, ecological theories such as the metabolic theory of ecology have focused on the effects of temperature rather than oxygen. This gap currently impedes mechanistic models from accurately predicting metabolic rates (i.e., oxygen consumption rates) of aquatic organisms and restricts predictions to resting metabolism, which is less affected by oxygen limitation. Here, we expand on models of metabolic scaling by accounting for the role of oxygen availability and temperature on both resting and active metabolic rates. Our model predicts that oxygen limitation is more likely to constrain metabolism in larger, warmer, and active fish. Consequently, active metabolic rates are less responsive to temperature than are resting metabolic rates, and metabolism scales to body size with a smaller exponent whenever temperatures or activity levels are higher. Results from a metaanalysis of fish metabolic rates are consistent with our model predictions. The observed interactive effects of temperature, oxygen availability, and body size predict that global warming will limit the aerobic scope of aquatic ectotherms and may place a greater metabolic burden on larger individuals, impairing their physiological performance in the future. Our model reconciles the metabolic theory with empirical observations of oxygen limitation and provides a formal, quantitative framework for predicting both resting and active metabolic rate and hence aerobic scope of aquatic ectotherms.
Aerobic metabolism fueling the physiological performance of ectotherms depends strongly on environmental temperature and can be limited by insufficient supply of oxygen (1, 2). Both temperature and oxygen availability are thus key environmental variables shaping the fitness of ectotherms and hence their geographical distributions (3⇓–5). The limiting effect of oxygen may be especially important for aquatic ectotherms, for which obtaining oxygen is more challenging (6⇓⇓⇓–10). However, mechanistic theories such as the metabolic theory of ecology (MTE), which have been widely adopted by ecologists to investigate organismal responses to climate change (11), have placed the thermal dependence of metabolic rates at the center of their synthesis (12, 13). While MTE successfully models the factors determining resting metabolism, predictions for active metabolic rates are far less accurate, likely because higher metabolic rates are progressively more constrained by oxygen supply. Here we extend models of the temperature dependence of metabolic scaling to include mechanistic descriptions of how oxygen uptake may constrain metabolic rates (i.e., measured as the rate of oxygen consumption) of active, warm, large-bodied aquatic ectotherms. We propose that such models will be crucial for making accurate predictions under future climate scenarios.
The temperature dependence of metabolism in many ecological theories relies on the Boltzmann–Arrhenius kinetics,
General frameworks of the MTE and the oxygen limitation hypothesis. While MTE focuses on the exponential relationship between RMRs and temperature, the oxygen limitation hypothesis emphasizes that metabolism is constrained by oxygen availability, especially at high temperatures or during intense activity (MMR). Therefore, MTE cannot be used to predict MMR (dashed line) and the aerobic scope (i.e., the difference between MMR and RMR, shaded areas), which according to the oxygen limitation hypothesis, represents the energy available for fitness-enhancing activities such as reproduction, feeding, or fighting.
It is well known that metabolic rates can be limited by insufficient oxygen, especially during intense activity or in hypoxia (1, 16), and that oxygen is linked to both how temperature affects metabolism (2, 17) and how metabolic rate scales with body size (18⇓⇓⇓–22). At high temperatures or under intense physical activity, the oxygen supply can no longer meet the increase in oxygen demand, at which point the metabolic rate approaches a limit (Fig. 1 and refs. 2, 5, 17, 23, 24). The difference between this upper limit, or maximum metabolic rate, and the resting metabolic rate is the aerobic scope, which represents the rate of energy available for fitness-enhancing activities like movement, feeding, digestion, growth, and reproduction. According to the oxygen limitation hypothesis, the aerobic scope is functionally related to physiological performance and fitness of aquatic ectotherms (2).
Although recent models have investigated the effects of hypoxia on physiological performance (25⇓–27), we still lack a mechanistic, quantitative framework that derives oxygen supply and demand from first principles. Here, we build upon the MTE equation of temperature and mass dependence (12) to show that oxygen limitation emerges from fundamental biophysical principles. We model the effect of oxygen limitation by introducing Michaelis–Menten kinetics for oxygen consumption. Although originally developed to describe saturation effects in individual chemical reactions, the Michaelis–Menten equation also captures how more complex, multicomponent systems, like metabolism, respond to limitation in supply (10, 28, 29). We use our model to investigate how oxygen limitation influences the temperature dependence and mass scaling of resting and maximum metabolic rates, and we show that the model successfully predicts patterns of resting and active metabolic rate across 286 species of fishes. More broadly, this analysis shows that predictions from the MTE represent a special case of a more general model where the affinity for oxygen is infinitely high.
Methods
Model Description.
The metabolic reactions that comprise whole-organism metabolic rate depend on temperature and the concentration of substrates such as carbohydrates and oxygen. In enzyme-catalyzed chemical reactions, the dependence on substrate concentration is given by the Hill kinetic equation,
where ρO2,gill is the concentration of oxygen at the gill surface (mg m−3). Eq. 1 incorporates an additional constant δ representing the increase in metabolic demand above the minimum (e.g., due to activities such as locomotion, feeding, or fighting). δ = 1 thus represents the minimum metabolic requirements (resting metabolic rate, RMR), whereas δ > 1 simulates increasing metabolic demand (up to maximum metabolic rate, MMR). Thus, the parameter δ is equivalent to the factorial aerobic scope (FAS = MMR/RMR), which provides a means to investigate the effect of oxygen limitation at different activity levels.
The key challenge is to model the oxygen concentration at the gill surface, ρO2,gill, as a function of environmental variables together with the animal’s morphology and physiology. In water-breathing ectotherms, oxygen intake depends on the rate of convective oxygen transfer, which involves both simple diffusion governed by gradients in partial pressure of oxygen, and the motion of water powered by branchial ventilation (6, 7, 30⇓–32). To model this convective process, we need to apply an analog of Fick’s law that also incorporates the advective motion of water at the gill surface:
where
Now we have the ingredients to model oxygen consumption by metabolism,
where ρO2,gill* is the steady-state concentration of oxygen at the respiratory surface. The left side of Eq. 3 contains the relationship between metabolic rate and both temperature (i.e., the Boltzmann–Arrhenius equation) and oxygen availability (i.e., the Michalis–Menten kinetic equation). The right side describes the rate of oxygen intake from the environment, modeled here via boundary-layer effects on external oxygen transport. A KM of 0 thus indicates that the affinity for oxygen is infinitely high and that there is no oxygen limitation. Under this condition, metabolic rate increases exponentially with temperature following the Arrhenius equation, consistent with how this relationship is expressed in the MTE. Higher KM values make metabolism increasingly dependent on the rate of oxygen intake given by the right side of the equation. This equation also contains several variables that the animal can potentially alter to increase oxygen intake rates, for example, increasing the velocity of the current of water at the gill surface, v (increasing ventilation rates) or the surface area of the respiratory organs, A (increasing exposure of respiratory organs). We analyzed the effect of each parameter by exploring Eq. 3 numerically (SI Appendix, Numerical analyses).
Model Predictions.
We first calibrated the model by comparing predicted and observed metabolic rates of eight species of teleost fishes to obtain realistic values of model parameters b0, E, KM, and δ (SI Appendix, Model calibration). We then used the model to analyze how incorporating oxygen limitation changes both 1) the thermal dependence of metabolism and 2) the mass scaling of metabolism. Following MTE, we used Arrhenius plots, which present metabolic rates as a function of inverse temperature, (kT) −1, where k is the Boltzmann constant (eV K−1) and T is absolute temperature (K). The slope of this relationship describes the sensitivity of metabolism to changes in temperature (more negative slope means greater dependence on temperature) and, according to MTE, this slope is determined by the average activation energy of the enzymatic reactions of metabolism. In our model, this slope also depends on the availability of and affinity for oxygen and thereby on temperature and activity level (SI Appendix, Numerical analyses). Specifically, at higher temperatures, more intense activity, or both, steady-state oxygen uptake rates are more and more constrained by oxygen availability, reducing the slope of the relationship of metabolism vs. inverse temperature (Fig. 2A). This predicted effect of oxygen limitation on the slope indicates that we should expect a significant interaction between temperature and activity level on metabolic rate.
Biophysical model predictions for temperature dependence (A) and mass scaling (B) of resting, RMR, and MMR. RMRs are predicted to exhibit greater temperature dependence than MMRs, which are more constrained by oxygen limitation. The metabolic scaling exponent (the slope of the log-MR vs. log-mass relationship) is predicted to vary with temperature, especially for MMR, where greater oxygen limitation with increasing temperature and larger body size gives rise to a shallower scaling relationship.
The predicted scaling of metabolic rates also depends on temperature and oxygen availability. Warmer water increases metabolism and thereby the intercept of the scaling relationship. However, the model predicts that the limits to oxygen supply are more readily experienced in larger fish, as indicated in a declining slope of the scaling relationship. Such a shallower scaling relationship is especially pronounced under intense activity and in warm waters (Fig. 2B). Specifically, the mass-scaling slope falls from its initial value, set here to b = 0.75, and approaches ∼0.60 with increasing demand for oxygen. These specific values of the mass-scaling slope in relation to temperature also depend on how other traits involved in oxygen transport scale with body size (30). We performed an additional sensitivity analysis to examine how the mass scaling of gill surface (A), velocity of the gill water flow (v), and activity level (δ) influence the predicted relationship between metabolism and body mass (SI Appendix, Sensitivity analysis). This predicted effect of temperature on the slope of the allometric scaling means that we should expect an interaction between temperature and mass on metabolic rate, and that this interaction should be more evident on maximum than on resting metabolic rates.
Data Analyses.
We tested model predictions for temperature dependence and mass scaling of RMR and MMR using empirical data from 286 species of fishes. These data were taken from White et al. (35) and supplemented with data from FishBase (36). The full database includes both intra- and interspecific data comprising 3,291 RMR measurements (241 species) and 1,212 MMR measurements (95 species).
We analyzed the effects of temperature, body mass, and activity level (either RMR or MMR) on oxygen consumption rate. We first included the three independent variables and analyzed their interactions to test the first prediction of the model (metabolic rates are determined by an interaction between activity level and temperature). We then split the data and performed separate analyses for RMR and MMR to test the second prediction of the model (the interaction between temperature and mass is more important for maximum than resting metabolic rates). We used mixed-effect models including the species identity as a random factor, allowing for both random intercepts and random slopes for each species. We also controlled for the phylogenetic relatedness between species using phylogenetic information in a subset of 103 fish species (SI Appendix, Phylogenetic analysis). We used the R package MCMCglmm (37) to perform both phylogenetically controlled and ordinary mixed-effect models in R v3.5.3 (38).
Results
Temperature Dependence.
Fish metabolic rates (both resting and maximum) were negatively and linearly related to inverse temperature (1/kT), which means that increasing water temperature increases oxygen demand (Table 1). The slope of this relationship, however, varied between resting and maximum metabolic rates; resting metabolic increased with temperature more strongly than maximum metabolic rates (Fig. 3A). This interaction term suggests that oxygen is more likely to be limiting for maximum metabolic rates, reducing their sensitivity to temperature. This result supports the first prediction of our model (Fig. 3A).
Results of the comparative mixed-effects models relating metabolic rates with inverse temperature (1/kT), log-body mass (mass), and activity level (resting, RMR, and MMR)
Empirical patterns for temperature dependence (A) and mass scaling (B) of resting, RMR, and MMR of fishes. Lines represent fits to empirical data (note that analyses were performed using linear mixed models that control for species identity).
Mass Scaling.
Although both RMR and MMR scaled with body mass, the scaling exponent varied with temperature, as indicated by the significant interaction between mass and temperature on MMR (Table 1). Specifically, increasing water temperature reduced the slope of the mass-scaling relationship and this effect was more evident for MMRs (Table 1). This result supports the second prediction of our model that oxygen limitation has a greater impact on large animals at higher water temperature or under intense activity (Fig. 3B). The results of both temperature dependence and the mass scaling were consistent with those obtained from phylogenetic models (SI Appendix, Phylogenetic analysis).
Discussion
The model provides quantitative predictions to evaluate how oxygen availability may shape the temperature dependence and mass scaling of metabolism in aquatic ectotherms. The model’s key predictions are that MMR is less sensitive to temperature variation than is RMR, and that metabolism scales to body size with a smaller exponent whenever temperatures are warmer or activity levels are higher (Fig. 2). Predictions on both temperature and size dependence were supported by a metaanalysis of metabolic rates in 286 species of fishes (Fig. 3). We propose that a parsimonious explanation for the observed interactions between body size, water temperature, and activity is the influence of physical processes shaping oxygen uptake dynamics in water. By integrating oxygen limitation (Michaelis–Menten kinetics), temperature dependence (Arrhenius kinetics), and oxygen transport dynamics, our model provides a mechanistic framework that integrates effects of temperature dependence, mass scaling, and activity on the aerobic metabolism of aquatic ectotherms.
Including both temperature and substrate dependence in metabolic models provides a more general view and better predictions than theories focused on temperature dependence alone (28). The first attempt to include oxygen limitation in the MTE was the metabolic index by Deutsch et al. (3), which focused on the net effects of oxygen supply and demand rather than on modeling each explicitly. Their metabolic index highlighted the need to incorporate both water temperature and oxygenation as determinants of species geographical ranges, as animals can lose aerobic scope due to warming and deoxygenation (3, 27). Dynamic energy budget theory (DEB) offers another framework to model the effects of metabolic fluxes of mass and energy on self-maintenance, growth, and reproduction (39, 40). Recent models have integrated the effect of hypoxia on metabolic fluxes (26) within the context of DEB theory (25), relying on indices and empirical formulas to describe the relationship between oxygen consumption and availability, rather than deriving them from underlying biophysical and chemical processes. Here we take a first-principles approach to modeling rates of oxygen supply and consumption. By linking metabolism to environmental factors and the physical properties of water that determine oxygen supply, such an approach offers the possibility of integrating temperature and oxygen dependence more mechanistically into theoretical frameworks such as MTE or DEB.
Oxygen moves from environment to mitochondria via a cascade that includes external transport, absorption across respiratory surfaces, distribution by internal transport systems, and diffusion within tissues and cells. Constraints on oxygen supply can arise at any or all of these levels (20, 22). While the MTE focuses on the internal transport network to predict the scaling exponent (14), the role of external factors such as temperature and notably oxygen availability have received less attention (40). Here we show that modeling external transport explicitly has two major implications for understanding metabolic scaling under oxygen limitation. First, this approach takes into consideration key compensatory mechanisms used by aquatic ectotherms to mitigate constraints on oxygen supply, such as increasing gill surface area (A), and the velocity of water at the gill surface via active branchial ventilation (v), or reducing the activity level by adopting a lifestyle with lower oxygen demands (δ). It has been argued that oxygen supply capacity can evolve to meet demand at the current oxygen partial pressure (41), and the increase in gill surface area with mass may be critical to reducing the impacts of oxygen limitation (42, 43). Here, we find that, when considering ranges of variation of gill surface, velocity of gill water flow, and activity level that are representative for fish, size constraints on oxygen supply remain evident at higher temperatures and activity levels (SI Appendix, Sensitivity analysis), and the empirical results are consistent with the model prediction of increased susceptibility to oxygen limitation in larger, active, hot fish (Fig. 3). Although fish exhibit a variety of adaptations to enhance oxygen supply or reduce oxygen demand, they are apparently not able to fully mitigate such oxygen constraints. A likely reason is that these adaptations themselves carry costs and involve trade-offs with other life-history traits that ultimately limit the effectiveness of selection. Although our observational results are consistent with predictions based on well-established physical principles on mass transfer, further experimental research is necessary to provide definitive support for oxygen limitation as the primary mechanism driving these patterns.
The second key implication of modeling external transport results from the fact that the oxygen-transfer coefficient depends on variables such as water temperature, current velocity, flow regime, and body size, making metabolic scaling exponents dependent on environmental variables and the activity of the animal. Fry et al. (1) first noticed that metabolic scaling exponents depend on environmental factors such as salinity and the activity of the animal. The metabolic-level boundaries hypothesis by Glazier (18) later postulated that both intrinsic and extrinsic factors can affect metabolic scaling and predicted a negative relationship between the intercept and the slope of the scaling relationship (19, 44). Killen et al. (19), found support for this prediction for resting metabolic rates: as the intercept of the scaling relationship increased with temperature, the mass-scaling exponent decreased. Here, we propose that stronger constraints on oxygen supply at warmer temperatures and larger body sizes are the main mechanism underpinning this negative relationship, and further show that the effect is more evident in maximum than in resting metabolic rates (Fig. 3B). Because our model incorporates both internal (parameter b from MTE) and external constraints on oxygen supply, it provides a comprehensive framework for quantitatively predicting scaling exponents under oxygen limitation. The Arrhenius exponential equation adopted by the MTE describes the temperature dependence of metabolic rates when oxygen is not limiting. In contrast, our model better captures responses to higher temperatures, higher activity levels, and environments with low oxygen availability, representing an advance in the field of metabolic scaling over previous studies which investigated differences in scaling between RMR and MMR but did not consider temperature effects on oxygen limitation (18, 19, 44).
Because the diffusivity of oxygen is much lower in water than in air, the effects of oxygen limitation on both temperature dependence and scaling exponents of metabolic rates should be more apparent in aquatic than in terrestrial taxa. Thus, for example, Gillooly et al. (12) showed that the Arrhenius slopes for fishes and amphibians are less negative than in other taxa but did not provide an explanation. According to our model, these shallower slopes result from greater risks of oxygen limitation in aquatic and semiaquatic than in terrestrial organisms. A second-order extrapolation of our model is that, in aquatic ectotherms, larger individuals likely experience greater oxygen limitation whenever temperatures are warmer, or activity levels are higher. Therefore, large animals may have some aerobic scope left to fuel activity, growth, and reproduction in cold waters but not in warmer waters––a prediction that is consistent with the temperature-size rule, i.e., the tendency of many ectotherms to grow to smaller final sizes when reared in warmer temperatures (45). Thus, Forster et al. (46) found that the temperature-size rule is more pronounced in aquatic than terrestrial species. Furthermore, Horne et al. (47) found that body sizes of aquatic species generally decrease more than those of terrestrial organisms with both warming and decreasing latitude. Finally, the evidence supporting the role of oxygen limitation of thermal tolerance limits is stronger for water- than for air-breathing arthropods (8, 9). Overall, these observations support the conclusion that oxygen limitation is a major factor driving patterns of body-size variations and physiological traits in aquatic but not in terrestrial organisms. This result does not contradict the idea that oxygen supply capacity in animals has evolved to meet demand at the current oxygen partial pressure (41), but it suggests that many aquatic ectotherms can readily reach their maximum supply capacity during bursts of activity. Warmer water temperatures will make demand for oxygen increasingly outpace supply (7, 31), which will curtail performance, especially at larger body sizes.
The biological impacts of climate warming are mediated primarily through its effects on organismal physiology (11, 48), and our model has general implications for predicting these impacts. First, rising water temperatures are predicted to raise resting metabolic rates more than MMRs, which may reduce the overall aerobic scopes of aquatic ectotherms. Smaller aerobic scopes imply that less energy remains for supporting physiological performance and fitness-related activities (17, 23, 49). Second, large individuals are likely more susceptible to oxygen limitation and therefore they could experience a greater reduction in aerobic scope with warming. This prediction supports the idea that oxygen limitation is at least partly responsible for driving the reduction of body size in response to global warming (32, 45, 49). Finally, risks of oxygen limitation will be greater in aquatic ecosystems. Consequently, oxygen limitation of aerobic metabolism may be a fundamental mechanism driving the response of aquatic ectotherms to climate warming.
Data Availability.
All study data are included in the article and SI Appendix.
Acknowledgments
We thank Craig White for insightful comments on the study and the editor and three anonymous reviewers for their effort toward improving our manuscript. J.G.R. was funded by an European Commission’s Marie Curie Individual Fellowship (H2020-MSCA-IF-2018, 843094).
Footnotes
- ↵1To whom correspondence may be addressed. Email: jg.rubalcaba{at}gmail.com.
Author contributions: J.G.R. and H.A.W. designed research; J.G.R., W.C.E.P.V., and H.A.W. performed research; J.G.R. contributed new reagents/analytic tools; J.G.R. and W.C.E.P.V. analyzed data; J.G.R., W.C.E.P.V., A.J.H., B.S., and H.A.W. wrote the paper; J.G.R. designed theoretical model; and W.C.E.P.V., A.J.H., and B.S. collected data.
The authors declare no competing interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.2003292117/-/DCSupplemental.
- Copyright © 2020 the Author(s). Published by PNAS.
This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).
References
- ↵
- W. S. Hoar,
- D. J. Randall
- F. E. J. Fry
- ↵
- ↵
- C. Deutsch,
- A. Ferrel,
- B. Seibel,
- H. O. Pörtner,
- R. B. Huey
- ↵
- W. C. E. P. Verberk,
- I. Durance,
- I. P. Vaughan,
- S. J. Ormerod
- ↵
- R. Ern
- ↵
- P. Dejours
- ↵
- ↵
- ↵
- ↵
- N. Galic,
- T. Hawkins,
- V. E. Forbes
- ↵
- ↵
- J. F. Gillooly,
- J. H. Brown,
- G. B. West,
- V. M. Savage,
- E. L. Charnov
- ↵
- ↵
- G. B. West,
- J. H. Brown,
- B. J. Enquist
- ↵
- J. R. Banavar et al
- ↵
- ↵
- H. O. Pörtner,
- R. Knust
- ↵
- ↵
- ↵
- J. F. Gillooly,
- J. P. Gomez,
- E. V. Mavrodiev,
- Y. Rong,
- E. S. McLamore
- ↵
- ↵
- ↵
- ↵
- P. M. Schulte
- ↵
- Y. Thomas et al
- ↵
- P. Neubauer,
- K. H. Andersen
- ↵
- C. Deutsch,
- J. L. Penn,
- B. A. Seibel
- ↵
- ↵
- C. M. Wood
- ↵
- ↵
- ↵
- W. C. E. P. Verberk,
- D. Atkinson
- ↵
- H. E. García,
- L. I. Gordon
- ↵
- D. Randall,
- H. Lin,
- P. A. Wright
- ↵
- C. R. White et al
- ↵
- R. Froese,
- D. Pauly
- ↵
- ↵
- R Core Team
- ↵
- S. A. L. M. Kooijman
- ↵
- ↵
- B. A. Seibel,
- C. Deutsch
- ↵
- S. Lefevre,
- D. J. McKenzie,
- G. E. Nilsson
- ↵
- D. Pauly,
- W. W. L. Cheung
- ↵
- D. Glazier
- ↵
- ↵
- J. Forster,
- A. G. Hirst,
- D. Atkinson
- ↵
- ↵
- C. A. Deutsch et al
- ↵
- H. O. Pörtner,
- A. P. Farrell
Citation Manager Formats
Article Classifications
- Biological Sciences
- Ecology