Field and laboratory studies reveal interacting effects of stream oxygenation and warming on aquatic ectotherms

Aquatic ecological responses to climatic warming are complicated by interactions between thermal effects and other environmental stressors such as organic pollution and hypoxia. Laboratory experiments have demonstrated how oxygen limitation can set heat tolerance for some aquatic ectotherms, but only at unrealistic lethal temperatures and without field data to assess whether oxygen shortages might also underlie sublethal warming effects. Here, we test whether oxygen availability affects both lethal and nonlethal impacts of warming on two widespread Eurasian mayflies, Ephemera danica, Müller 1764 and Serratella ignita (Poda 1761). Mayfly nymphs are often a dominant component of the invertebrate assemblage in streams, and play a vital role in aquatic and riparian food webs. In the laboratory, lethal impacts of warming were assessed under three oxygen conditions. In the field, effects of oxygen availability on nonlethal impacts of warming were assessed from mayfly occurrence in 42 293 UK stream samples where water temperature and biochemical oxygen demand were measured. Oxygen limitation affected both lethal and sublethal impacts of warming in each species. Hypoxia lowered lethal limits by 5.5 °C (±2.13) and 8.2 °C (±0.62) for E. danica and S. ignita respectively. Field data confirmed the importance of oxygen limitation in warmer waters; poor oxygenation drastically reduced site occupancy, and reductions were especially pronounced under warm water conditions. Consequently, poor oxygenation lowered optimal stream temperatures for both species. The broad concordance shown here between laboratory results and extensive field data suggests that oxygen limitation not only impairs survival at thermal extremes but also restricts species abundance in the field at temperatures well below upper lethal limits. Stream oxygenation could thus control the vulnerability of aquatic ectotherms to global warming. Improving water oxygenation and reducing pollution can provide key facets of climate change adaptation for running waters.     

Hypoxia reduced upper thermal limits causing cellular and nuclear abnormalities of erythrocytes in Nile tilapia,Oreochromis niloticus

Global warming is a threat across the world that leads to estimates of the upper thermal limits of ectothermic species. Increased water temperature up-regulates oxygen consumption and metabolic rates, and alters the physiological processes. In this study, we identified the critical thermal maxima (CTmax) and physiological responses under normoxia and hypoxia in Nile tilapia, Oreochromis niloticus. CTmax was 41.25 °C under hypoxia and 44.50 °C under normoxia. Compared to normoxia, lower values of hemoglobin (Hb) and red blood cells (RBCs) were observed at the CTmax under hypoxia. In contrast, higher values of white blood cells (WBCs) and blood glucose (Glu) levels were observed at the CTmax under hypoxia. Consequently, higher frequencies of micronucleus, cellular and nuclear abnormalities of erythrocytes were observed at the CTmax under hypoxia. These results suggest that high temperature tolerance and subsequent physiology are significantly affected by the oxygen supply in Nile tilapia. As climate vulnerability is intensifying day by day, this data will be helpful in successful management practice for the aquatic environment having low oxygen content.     

Long-term forecast of thermal mortality with climate warming in riverine amphipods

Forecasting long-term consequences of global warming requires knowledge on thermal mortality and how heat stress interacts with other environmental stressors on different timescales. Here, we describe a flexible analytical framework to forecast mortality risks by combining laboratory measurements on tolerance and field temperature records. Our framework incorporates physiological acclimation effects, temporal scale differences and the ecological reality of fluctuations in temperature, and other factors such as oxygen. As a proof of concept, we investigated the heat tolerance of amphipods Dikerogammarus villosus and Echinogammarus trichiatus in the river Waal, the Netherlands. These organisms were acclimated to different temperatures and oxygen levels. By integrating experimental data with high-resolution field data, we derived the daily heat mortality probabilities for each species under different oxygen levels, considering current temperatures as well as 1 and 2°C warming scenarios. By expressing heat stress as a mortality probability rather than a upper critical temperature, these can be used to calculate cumulative annual mortality, allowing the scaling up from individuals to populations. Our findings indicate a substantial increase in annual mortality over the coming decades, driven by projected increases in summer temperatures. Thermal acclimation and adequate oxygenation improved heat tolerance and their effects were magnified on longer timescales. Consequently, acclimation effects appear to be more effective than previously recognized and crucial for persistence under current temperatures. However, even in the best-case scenario, mortality of D. villosus is expected to approach 100% by 2100, while E. trichiatus appears to be less vulnerable with mortality increasing to 60%. Similarly, mortality risks vary spatially: In southern, warmer rivers, riverine animals will need to shift from the main channel toward the cooler head waters to avoid thermal mortality. Overall, this framework generates high-resolution forecasts on how rising temperatures, in combination with other environmental stressors such as hypoxia, impact ecological communities.     

Warming will affect phytoplankton differently: evidence through a mechanistic approach

Although the consequences of global warming in aquatic ecosystems are only beginning to be revealed, a key to forecasting the impact on aquatic communities is an understanding of individual species' vulnerability to increased temperature. Despite their microscopic size, phytoplankton support about half of the global primary production, drive essential biogeochemical cycles and represent the basis of the aquatic food web. At present, it is known that phytoplankton are important targets and, consequently, harbingers of climate change in aquatic systems. Therefore, investigating the capacity of phytoplankton to adapt to the predicted warming has become a relevant issue. However, considering the polyphyletic complexity of the phytoplankton community, different responses to increased temperature are expected. We experimentally tested the effects of warming on 12 species of phytoplankton isolated from a variety of environments by using a mechanistic approach able to assess evolutionary adaptation (the so-called ratchet technique). We found different degrees of tolerance to temperature rises and an interspecific capacity for genetic adaptation. The thermal resistance level reached by each species is discussed in relation to their respective original habitats. Our study additionally provides evidence on the most resistant phytoplankton groups in a future warming scenario.     

Cadmium-dependent oxygen limitation affects temperature tolerance in eastern oysters (Crassostrea virginica Gmelin)

Marine ectotherms, including oysters are exposed to variable environmental conditions in coastal shallow waters and estuaries. In the light of global climate change, additional stressors like pollution might pose higher risk to populations. On the basis of the concept of oxygen- and capacity-limited thermal tolerance in aquatic ectotherms (40), we show that a persistent pollutant, cadmium, can have detrimental effects on oysters (Crassostrea virginica). During acute warming from 20 to 28 degrees C (4 degrees C/48 h) standard metabolic rate (SMR) rose in control and cadmium-exposed (50 microg Cd2+/l) animals, with a consistently higher SMR in Cd-exposed oysters. Additionally, Cd-exposed oysters showed a stronger temperature-dependent decrease in hemolymph oxygen partial pressures. This observation indicates that the effect of temperature on aerobic metabolism was exacerbated due to the additional Cd stress. The oxygen delivery systems could not provide enough oxygen to cover Cd-induced elevated metabolic demands at high temperatures. Interestingly, cardiac performance (measured as the heart rate and hemolymph supply to tissues) rose to a similar extent in control and Cd-exposed oysters with warming indicating that cardiac output was unable to compensate for elevated energy demand in Cd-exposed oysters. Together with the literature data on metal-induced reduction of ventilatory capacity, these findings suggest that synergistic effects of elevated temperatures and cadmium exposure led to oxygen limitation by impaired performance in oxygen supply through ventilation and circulation. Overall, cadmium exposure resulted in progressive hypoxemia in oysters at high temperatures, suggesting that the thermal tolerance window is narrowed in marine ectotherms inhabiting polluted areas compared with pristine environments.     

Linking species thermal tolerance to elevational range shifts in upland dung beetles

Climate warming has been proposed as the main cause of the recent range shifts seen in many species. Although species' thermal tolerances are thought to play a key role in determining responses to climate change, especially in ectotherms, empirical evidence is still limited. We investigate the connection between species' thermal tolerances, elevational range and shifts in the lower elevational limit of dung beetle species (Coleoptera, Aphodiidea) in an upland region in the northwest of England. We measured thermal tolerances in the laboratory, and used current and historical distribution data to test specific hypotheses about the area's three dominant species, particularly the species most likely to suffer from warming: Agollinus lapponum. We found marked differences between species in their minimum and maximum thermal tolerance and in their elevational range and patterns of abundance. Overall, differences in thermal limits among species matched the abundance patterns along the elevation gradient expected if distributions were constrained by climate. Agollinus lapponum abundance increased with elevation and this species showed lower maximum and minimum thermal limits than Acrossus depressus, for which abundance declined with elevation. Consistent with lower tolerance to high temperature, we recorded an uphill retreat of the low elevation limit of A. lapponum (177 m over 57 yr) in line with the increase in summer temperature observed in the region over the same period. Moreover, this species has been replaced at low and mid‐elevations by the other two warm‐tolerant species (A. depressus and Agrilinus ater). Our results provide empirical evidence that species' thermal tolerance constrains elevational ranges and contributes to explain the observed responses to climate warming. A mechanistic understanding of how climate change directly affects species, such as the one presented here, will provide a robust base to inform predictions of how individual species and whole assemblages may change in the future.     

Aquatic hyphomycete strains from metal-contaminated and reference streams might respond differently to future increase in temperature

Aquatic hyphomycetes, a group of polyphyletic fungi, have been reported in streams contaminated with metals. This tolerance to metal contamination however can result in limited performance and limited ability to cope with additional environmental change. The predicted increase in water temperature, as a consequence of global warming, will have an additional effect on many streams. The sensitivity to temperature of strains of three aquatic hyphomycete species isolated from a metal-contaminated stream and an uncontaminated stream was assessed by determining their radial growth and activity (conidial production, oxygen consumption, mycelial biomass accumulation, fine particulate organic matter [FPOM] production, and microbial induced leaf mass loss) at 13 C (present water temperature in autumn) and at 18 C (predicted water temperature under global warming). Growth and reproductive activity generally were depressed for the strains isolated from the metal-contaminated stream when compared with those isolated from the unpolluted stream. These differences however were not translated into differences in FPOM production and leaf-litter mass loss, indicating that the strains isolated from the contaminated stream can decompose leaf litter similar to those of the reference stream. The 5 C increase in temperature stimulated fungal activity and litter decomposition, irrespective of species and strain. This might have strong effect on aquatic food-web and ecosystem functioning under global warming because increases in litter decomposition might lead to food shortage for higher trophic levels. The sensitivity to temperature depended on the response variable, species and strain. FPOM production was the variable most sensitive to temperature across strains and species and that for which temperature sensitivities differed most between strains. Fungal tolerance to metal contamination affects the extent to which its functions are stimulated by an increase in temperature, constituting an additional cost of metal tolerance.     

Improved heat tolerance in air drives the recurrent evolution of air-breathing

The transition to air-breathing by formerly aquatic species has occurred repeatedly and independently in fish, crabs and other animal phyla, but the proximate drivers of this key innovation remain a long-standing puzzle in evolutionary biology. Most studies attribute the onset of air-breathing to the repeated occurrence of aquatic hypoxia; however, this hypothesis leaves the current geographical distribution of the 300 genera of air-breathing crabs unexplained. Here, we show that their occurrence is mainly related to high environmental temperatures in the tropics. We also demonstrate in an amphibious crab that the reduced cost of oxygen supply in air extends aerobic performance to higher temperatures and thus widens the animal''s thermal niche. These findings suggest that high water temperature as a driver consistently explains the numerous times air-breathing has evolved. The data also indicate a central role for oxygen- and capacity-limited thermal tolerance not only in shaping sensitivity to current climate change but also in underpinning the climate-dependent evolution of animals, in this case the evolution of air-breathing.     

Hypoxia tolerance of three native and three alien species of bitterling inhabiting Lake Kasumigaura,Japan

Differences in hypoxia tolerance among three native and three alien bitterling species were examined by means of field surveys and aquarium experiments. I caught fish in minnow traps and measured environmental factors such as dissolved oxygen and current velocity once a month between June and November 2006 at 30 points around Lake Kasumigaura, Japan. In addition, aquarium experiments were carried out to compare aquatic surface respiration thresholds among the bitterling species. Much more individuals of an alien species, Rhodeus ocellatus ocellatus, was caught at the least oxygen range than expected in the field, and showed the highest hypoxia tolerance in the laboratory experiments. Another alien bitterling, Acheilognathus rhombeus, also showed higher hypoxia tolerance than the three native bitterling species. Lake Kasumigaura is well known for eutrophication and water masses with low dissolved oxygen have often been observed there. Differences in hypoxia tolerance among bitterling species might have affected bitterling assemblage structure in Lake Kasumigaura, which is currently characterized by the dominance of R. o. ocellatus.     

Marine species in ambient low‐oxygen regions subject to double jeopardy impacts of climate change

We have learned much about the impacts of warming on the productivity and distribution of marine organisms, but less about the impact of warming combined with other environmental stressors, including oxygen depletion. Also, the combined impact of multiple environmental stressors requires evaluation at the scales most relevant to resource managers. We use the Gulf of St. Lawrence, Canada, characterized by a large permanently hypoxic zone, as a case study. Species distribution models were used to predict the impact of multiple scenarios of warming and oxygen depletion on the local density of three commercially and ecologically important species. Substantial changes are projected within 20–40 years. A eurythermal depleted species already limited to shallow, oxygen‐rich refuge habitat (Atlantic cod) may be relatively uninfluenced by oxygen depletion but increase in density within refuge areas with warming. A more stenothermal, deep‐dwelling species (Greenland halibut) is projected to lose ~55% of its high‐density areas under the combined impacts of warming and oxygen depletion. Another deep‐dwelling, more eurythermal species (Northern shrimp) would lose ~4% of its high‐density areas due to oxygen depletion alone, but these impacts may be buffered by warming, which may increase density by 8% in less hypoxic areas, but decrease density by ~20% in the warmest parts of the region. Due to local climate variability and extreme events, and that our models cannot project changes in species sensitivity to hypoxia with warming, our results should be considered conservative. We present an approach to effectively evaluate the individual and cumulative impacts of multiple environmental stressors on a species‐by‐species basis at the scales most relevant to managers. Our study may provide a basis for work in other low‐oxygen regions and should contribute to a growing literature base in climate science, which will continue to be of support for resource managers as climate change accelerates.     

Unshelled abalone and corrupted urchins: development of marine calcifiers in a changing ocean

The most fragile skeletons produced by benthic marine calcifiers are those that larvae and juveniles make to support their bodies. Ocean warming, acidification, decreased carbonate saturation and their interactive effects are likely to impair skeletogenesis. Failure to produce skeleton in a changing ocean has negative implications for a diversity of marine species. We examined the interactive effects of warming and acidification on an abalone (Haliotis coccoradiata) and a sea urchin (Heliocidaris erythrogramma) reared from fertilization in temperature and pH/pCO(2) treatments in a climatically and regionally relevant setting. Exposure of ectodermal (abalone) and mesodermal (echinoid) calcifying systems to warming (+2°C to 4°C) and acidification (pH 7.6-7.8) resulted in unshelled larvae and abnormal juveniles. Haliotis development was most sensitive with no interaction between stressors. For Heliocidaris, the percentage of normal juveniles decreased in response to both stressors, although a +2°C warming diminished the negative effect of low pH. The number of spines produced decreased with increasing acidification/pCO(2), and the interactive effect between stressors indicated that a +2°C warming reduced the negative effects of low pH. At +4°C, the developmental thermal tolerance was breached. Our results show that projected near-future climate change will have deleterious effects on development with differences in vulnerability in the two species.     

Oxygen supply limits the heat tolerance of lizard embryos

The mechanisms that set the thermal limits to life remain uncertain. Classically, researchers thought that heating kills by disrupting the structures of proteins or membranes, but an alternative hypothesis focuses on the demand for oxygen relative to its supply. We evaluated this alternative hypothesis by comparing the lethal temperature for lizard embryos developing at oxygen concentrations of 10–30%. Embryos exposed to normoxia and hyperoxia survived to higher temperatures than those exposed to hypoxia, suggesting that oxygen limitation sets the thermal maximum. As all animals pass through an embryonic stage where respiratory and cardiovascular systems must develop, oxygen limitation may limit the thermal niches of terrestrial animals as well as aquatic ones.     

Mitochondrial manganese superoxide dismutase from the shrimp Litopenaeus vannamei: Molecular characterization and effect of high temperature,hypoxia and reoxygenation on expression and enzyme activity

Climate warming has been increasing ocean water temperature and decreasing oxygen concentrations, exposing aquatic organisms to environmental stress conditions. The shrimp Litopenaeus vannamei manages to survive these harsh environmental conditions by enhancing their antioxidant defenses, among other strategies. In this study, we report the mitochondrial manganese superoxide dismutase (mMnSOD) nucleotide and deduced amino acid sequences and its gene expression in L. vannamei tissues. The deduced protein has 220 amino acids with a signal peptide of 20 amino acids. Expression of mMnSOD was analyzed in hepatopancreas, gills and muscle, where gills had highest expression in normoxic conditions. In addition, shrimp were subjected to high temperature, hypoxia and reoxygenation to analyze the effect on the expression of mMnSOD and SOD activity in mitochondria. High temperature and hypoxia showed a synergistic effect in the up-regulation on expression of mMnSOD in gills and hepatopancreas. Moreover, induction in SOD activity was found in the mitochondrial fraction from gills of normoxia at high temperature, probably due to an overproduction of reactive oxygen species caused by an elevated metabolic rate due to the stress temperature. These results suggest that the combined stress conditions of hypoxia and high temperature trigger molecularly the antioxidant response in L. vannamei in a higher degree than only one stressor.     

Sensitivity to thermal extremes in Australian Drosophila implies similar impacts of climate change on the distribution of widespread and tropical species

Climatic factors influence the distribution of ectotherms, raising the possibility that distributions of many species will shift rapidly under climate change and/or that species will become locally extinct. Recent studies have compared performance curves of species from different climate zones and suggested that tropical species may be more susceptible to climate change than those from temperate environments. However, in other comparisons involving responses to thermal extremes it has been suggested that mid‐latitude populations are more susceptible. Using a group of 10 closely related Drosophila species with known tropical or widespread distribution, we undertake a detailed investigation of their growth performance curves and their tolerance to thermal extremes. Thermal sensitivity of life history traits (fecundity, developmental success, and developmental time) and adult heat resistance were similar in tropical and widespread species groups, while widespread species had higher adult cold tolerance under all acclimation regimes. Laboratory measurements of either population growth capacity or acute tolerance to heat and cold extremes were compared to daily air temperature under current (2002–2007) and future (2100) conditions to investigate if these traits could explain current distributions and, therefore, also forecast future effects of climate change. Life history traits examining the thermal sensitivity of population growth proved to be a poor predictor of current species distributions. In contrast, we validate that adult tolerance to thermal extremes provides a good correlate of current distributions. Thus, in their current distribution range, most of the examined species experience heat exposure close to, but rarely above, the functional heat resistance limit. Similarly, adult functional cold resistance proved a good predictor of species distribution in cooler climates. When using the species’ functional tolerance limits under a global warming scenario, we find that both tropical and widespread Drosophila species will face a similar proportional reduction in distribution range under future warming.     

The effects of changing climate on faunal depth distributions determine winners and losers

Changing climate is predicted to impact all depths of the global oceans, yet projections of range shifts in marine faunal distributions in response to changing climate seldom evaluate potential shifts in depth distribution. Marine ectotherms' thermal tolerance is limited by their ability to maintain aerobic metabolism (oxygen‐ and capacity‐limited tolerance), and is functionally associated with their hypoxia tolerance. Shallow‐water (<200 m depth) marine invertebrates and fishes demonstrate limited tolerance of increasing hydrostatic pressure (pressure exerted by the overlying mass of water), and hyperbaric (increased pressure) tolerance is proposed to depend on the ability to maintain aerobic metabolism, too. Here, we report significant correlation between the hypoxia thresholds and the hyperbaric thresholds of taxonomic groups of shallow‐water fauna, suggesting that pressure tolerance is indeed oxygen limited. Consequently, it appears that the combined effects of temperature, pressure and oxygen concentration constrain the fundamental ecological niches (FENs) of marine invertebrates and fishes. Including depth in a conceptual model of oxygen‐ and capacity‐limited FENs' responses to ocean warming and deoxygenation confirms previous predictions made based solely on consideration of the latitudinal effects of ocean warming (e.g. Cheung et al., 2009), that polar taxa are most vulnerable to the effects of climate change, with Arctic fauna experiencing the greatest FEN contraction. In contrast, the inclusion of depth in the conceptual model reveals for the first time that temperate fauna as well as tropical fauna may experience substantial FEN expansion with ocean warming and deoxygenation, rather than FEN maintenance or contraction suggested by solely considering latitudinal range shifts.     

Sub‐lethal effects on fish provide insight into a biologically‐relevant threshold of hypoxia

Hypoxia (low dissolved oxygen) is a mounting concern for aquatic ecosystems as its prevalence increases with rising anthropogenic nutrient inputs. Hypoxia is most commonly defined as 2.0 mg l^–1 of dissolved oxygen, although this level varies widely across studies and agency regulations. Such definitions may be too conservative, as ecologically‐relevant non‐lethal effects (e.g. consumption and growth) of hypoxia on important aquatic species, such as fish, often occur at oxygen levels much higher than 2.0 mg l^–1. In addition, many mechanisms that regulate hypoxia tolerance in fish have been proposed, including temperature, habitat, location in the water column, and body size, but there is ongoing debate over which mechanisms are most important. Using a structured meta‐analysis of published studies, we showed consistent, significant negative effects on fish growth and consumption below 4.5 mg l^–1. While the total amount of variation explained was generally low, below 4.5 mg l^–1 of dissolved oxygen, phylogenetic relationships accounted for most of the explained variation in fish growth. Ecological factors including body size, location in the water column (pelagic, demersal, or benthopelagic), habitat (freshwater, marine, or diadromous), and temperature explained very little of the effect of hypoxia on fish growth and explained only a moderate level of variation in consumption. Our results suggest a dramatically higher threshold for sub‐lethal effects of hypoxia on fish than oxygen levels generally set for regulation purposes, and provide little support for accepted ecological mechanisms thought to influence hypoxia tolerance.     

Does simultaneous UV-B exposure enhance the lethal and sub-lethal effects of aquatic hypoxia on developing anuran embryos and larvae?

Recent catastrophic global amphibian declines have been partially linked to increases in UV-B radiation as a consequence of stratospheric ozone depletion. Previous studies have shown that in the presence of other environmental stressors including aquatic pH and temperature and the presence of contaminants or pathogens, the lethal effects of UV-B on amphibian larvae are enhanced due to interactions between the stressors. Little is known about the interactions between UV-B and aquatic hypoxia, a common and significant natural stressor of amphibian larvae. We examined the potential effects of UV-B and aquatic hypoxia in combination on embryonic survival, developmental rate, body mass and locomotor performance of embryos and larvae of the striped marsh frog, Limnodynastes peronii. We found that while both UV-B and hypoxia independently had substantial negative effects on the developing embryos of L. peronii, they did not interact in a multiplicative or antagonistic manner. The effects of the stressors in combination were as might be predicted based on the knowledge of their independent actions alone (i.e. an additive effect). In all cases developing embryos exposed to both UV-B and hypoxia were more severely affected than those exposed to either UV-B or hypoxia alone. The results of this study show the importance of examining both the direct actions of individual stressors and how these may be influenced by the presence of other environmental factors.     

Testing local and global stressor impacts on a coastal foundation species using an ecologically realistic framework

Despite the abundance of literature on organismal responses to multiple environmental stressors, most studies have not matched the timing of experimental manipulations with the temporal pattern of stressors in nature. We test the interactive effects of diel‐cycling hypoxia with both warming and decreased salinities using ecologically realistic exposures. Surprisingly, we found no evidence of negative synergistic effects on Olympia oyster growth; rather, we found only additive and opposing effects of hypoxia (detrimental) and warming (beneficial). We suspect that diel‐cycling provided a temporal refuge that allowed physiological compensation. We also tested for latent effects of warming and hypoxia to low‐salinity tolerance using a seasonal delay between stressor events. However, we did not find a latent effect, rather a threshold survival response to low salinity that was independent of early life‐history exposure to warming or hypoxia. The absence of synergism is likely the result of stressor treatments that mirror the natural timing of environmental stressors. We provide environmental context for laboratory experimental data by examining field time series environmental data from four North American west coast estuaries and find heterogeneous environmental signals that characterize each estuary, suggesting that the potential stressor exposure to oysters will drastically differ over moderate spatial scales. This heterogeneity implies that efforts to conserve and restore oysters will require an adaptive approach that incorporates knowledge of local conditions. We conclude that studies of multiple environmental stressors can be greatly improved by integrating ecologically realistic exposure and timing of stressors found in nature with organismal life‐history traits.     

Effects of oxygen on responses to heating in two lizard species sampled along an elevational gradient

Thermal tolerance is an important variable in predictive models about the effects of global climate change on species distributions, yet the physiological mechanisms responsible for reduced performance at high temperatures in air-breathing vertebrates are not clear. We conducted an experiment to examine how oxygen affects three variables exhibited by ectotherms as they heat—gaping threshold, panting threshold, and loss of righting response (the latter indicating the critical thermal maximum)—in two lizard species along an elevational (and therefore environmental oxygen partial pressure) gradient. Oxygen partial pressure did not impact these variables in either species. We also exposed lizards at each elevation to severely hypoxic gas to evaluate their responses to hypoxia. Severely low oxygen partial pressure treatments significantly reduced the gaping threshold, panting threshold, and critical thermal maximum. Further, under these extreme hypoxic conditions, these variables were strongly and positively related to partial pressure of oxygen. In an elevation where both species overlapped, the thermal tolerance of the high elevation species was less affected by hypoxia than that of the low elevation species, suggesting the high elevation species may be adapted to lower oxygen partial pressures. In the high elevation species, female lizards had higher thermal tolerance than males. Our data suggest that oxygen impacts the thermal tolerance of lizards, but only under severely hypoxic conditions, possibly as a result of hypoxia-induced anapyrexia.