High altitude environments challenge small mammals with persistent low ambient temperatures that require high rates of aerobic heat production in face of low O2 availability. An important component of thermogenic capacity in rodents is non-shivering thermogenesis (NST) mediated by uncoupled mitochondrial respiration in brown adipose tissue (BAT). NST is plastic, and capacity for heat production increases with cold acclimation. However, in lowland native rodents, hypoxia inhibits NST in BAT. We hypothesize that highland deer mice (Peromyscus maniculatus) overcome the hypoxic inhibition of NST through changes in BAT mitochondrial function. We tested this hypothesis using lab born and raised highland and lowland deer mice, and a lowland congeneric (Peromyscus leucopus), acclimated to either warm normoxia (25°C, 760 mmHg) or cold hypoxia (5°C, 430 mmHg). We determined the effects of acclimation and ancestry on whole-animal rates of NST, the mass of interscapular BAT (iBAT), and uncoupling protein (UCP)-1 protein expression. To identify changes in mitochondrial function, we conducted high-resolution respirometry on isolated iBAT mitochondria using substrates and inhibitors targeted to UCP-1. We found that rates of NST increased with cold hypoxia acclimation but only in highland deer mice. There was no effect of cold hypoxia acclimation on iBAT mass in any group, but highland deer mice showed increases in UCP-1 expression and UCP-1-stimulated mitochondrial respiration in response to these stressors. Our results suggest that highland deer mice have evolved to increase the capacity for NST in response to chronic cold hypoxia, driven in part by changes in iBAT mitochondrial function.

The generation of metabolic heat (thermogenesis) is an important mechanism used by endotherms to maintain a high body temperature in face of varying environmental temperatures. In small mammals, non-shivering thermogenesis (NST) is a significant component of whole-animal thermogenic capacity (Cannon and Nedergaard, 2004; Janský, 1973; Himms-Hagen, 1985; Heldmaier, 1971; Heldmaier et al., 1989). In fact, NST can contribute more than 50% of cold-induced O2,max (thermogenic capacity) in some rodents (Van Sant and Hammond, 2008). The major site of NST in eutherian mammals is brown adipose tissue (BAT), where energy from uncoupled mitochondrial respiration, mediated by uncoupling protein 1 (UCP-1) is almost entirely released as heat (Cannon and Nedergaard, 2004). Heat production by NST occurs as the result of activation of UCP-1 by free fatty acids released in response to sympathetic activation of BAT upon cold exposure (Cannon and Nedergaard, 2004; Fedorenko et al., 2012).

Thermogenic capacity is a plastic trait that is influenced by multiple environmental factors. Prolonged exposure to cold is known to increase overall thermogenic capacity in rodents (Cannon and Nedergaard, 2011), which helps in the preparation for seasonal variation in ambient temperatures (e.g. Böckler and Heldmaier, 1983; Klaus et al., 1988). A high thermogenic capacity is especially important for those species that are winter active as it correlates with increased activity levels in the cold (Sears et al., 2006). In rodents, seasonal variation in thermogenic capacity is primarily driven by changes in NST, through remodeling of BAT (Klaus et al., 1988; Van Sant and Hammond, 2008). Studies on wild white-footed mice and bank voles show that NST capacity increases in preparation for winter (Klaus et al., 1988; Wickler, 1980). This plasticity of NST has also been observed in cold-acclimated lab mice (e.g. Beaudry and McClelland, 2010; Lacy et al., 1978; Oufara et al., 1987). The mechanisms that underlie increases in the capacity for NST are thought to involve either changes in: (i) the size of BAT depots, (ii) brown adipocyte phenotype, and/or (iii) mitochondrial function. For example, cold acclimation in rodents has been shown to increase BAT mass by both hypertrophy (Klingenspor et al., 1996) and hyperplasia (Bukowiecki et al., 1982), but also by increasing tissue oxidative capacity through greater mitochondrial abundance and/or the mitochondrial expression of UCP-1 (Beaudry and McClelland, 2010; Klingenspor et al., 1996; Smith and Roberts, 1964). The respiratory capacity of isolated BAT mitochondria may also be a target of plasticity, through increases in UCP-1-induced respiration and capacity for substrate catabolism (e.g. β-oxidation of fatty acids), enhancing brown adipocyte heat production (Ballinger et al., 2016; Heim et al., 2017; Yu et al., 2002).

High altitude regions are generally cooler year-round than those of low altitude at the same latitude (Hayes, 1989). Since low oxygen availability is known to restrict aerobic metabolism, the high alpine climate presents a significant thermoregulatory challenge for small mammals inhabiting these regions. Indeed, deer mice [Peromyscus maniculatus (Wagner 1945)] native to high altitude have higher daily energy requirements than their low altitude conspecifics, presumably because of the high demand for thermogenesis (Hayes, 1989). This high aerobic demand has likely contributed to directional selection for enhanced aerobic capacity seen in high altitude deer mice (Hayes and O'Connor, 1999). As a result, high altitude deer mice have evolved a higher thermogenic capacity compared with both their low altitude conspecifics and a strictly low altitude congeneric [white-footed mice, Peromyscus leucopus (Rafinesque 1818)] (Cheviron et al., 2012, 2014). This higher capacity for aerobic thermogenesis in highland deer mice is associated with constitutively greater aerobic capacity of skeletal muscle (Cheviron et al., 2014; Lui et al., 2015; Scott et al., 2015), presumably to enhance the capacity for shivering. However, NST is a major component of thermogenesis in deer mice (Van Sant and Hammond, 2008), suggesting that BAT function also contributes to the greater thermogenic capacity in highlanders. Indeed, Velotta et al. (2016) found that wild highland deer mice had higher NST than lowland conspecifics when measured at their native altitudes. These population differences in NST may represent genetically fixed traits or may be the product of phenotypic plasticity induced by chronic exposure to cold hypoxia in the high alpine. Comparisons of NST and interscapular BAT (iBAT) transcriptomics between wild deer mice from low and high altitudes, and their first-generation laboratory born and raised descendants, suggest phenotypic plasticity (developmental and/or adult) plays an important role in the evolved higher thermogenic capacity seen at high altitude (Velotta et al., 2016).

In low altitude natives, hypoxia has antagonistic effects on cold-induced changes in NST and BAT remodeling (Beaudry and McClelland, 2010). This response is distinct from the effect of either cold or hypoxia acclimation alone (Beaudry and McClelland, 2010). While cold acclimation has a stimulatory effect on BAT growth and metabolic remodeling, hypoxia acclimation tends to have an inhibitory effect. Although the effects of hypoxia acclimation on NST have not been extensively studied, it has been shown to increase cytochrome c oxidase activity and stimulate UCP-1 mRNA expression in BAT but decreases the capacity for NST in CD-1 strain lab mice (Beaudry and McClelland, 2010). In contrast, when combined with cold, hypoxia led to a blunting of the cold-induced increases in NST (Beaudry and McClelland, 2010). Whether these antagonistic influences of combined cold hypoxia on NST and BAT function persist in highland natives in response to the combined stressors experienced in their native environment is presently unknown.

In this study, we used first generation lab born and raised lowland and highland deer mice (P. maniculatus), and the strictly lowland white-footed mice (P. leucopus) acclimated to warm normoxia control conditions and to cold hypoxia simulating high altitude. We tested the hypothesis that deer mice native to high altitude have evolved a robust response to cold hypoxia acclimation that results in an increase in NST due to BAT remodeling. We predicted that cold hypoxia acclimation will lead to higher rates of NST, greater BAT mass, and an augmented capacity for brown adipocytes and their mitochondria to increase heat production in highlanders. In contrast, we predict cold hypoxia will have antagonistic effects in lowlanders resulting in no net change in NST, BAT mass or its capacity for heat production.

Animal model

Captive breeding populations were established from wild-caught mice, as previously described (Cheviron et al., 2012; Lau et al., 2017; Tate et al., 2017). Lowland deer mice (P. maniculatus nebrascensis) and white-footed mice (P. leucopus) were captured at Nine-mile Prairie, NE, USA (320 m a.s.l.). Highland deer mice (P. maniculatus rufinus) were captured at the summit of Mount Evans, CO, USA at 4348 m a.s.l. Wild-caught mice were transported to McMaster University in Hamilton, ON, Canada (∼90 m a.s.l.) and housed in common laboratory conditions for breeding within their respective populations to the first generation. Mice were held under standard temperature (25°C) at sea level with a 12 h:12 h light:dark cycle and had access to standard mouse chow and water ad libitum.

Acclimation conditions

Mice of all three groups, which included adult mice from a mix of different families and sexes, were exposed to two acclimation conditions for 6–8 weeks: warm normoxia (WN, control conditions) and cold hypoxia (CH). Mice in the WN acclimation were held at atmospheric pressure of ∼101 kPa and at a room temperature of 25°C, while CH mice were held in hypobaric chambers (McClelland et al., 1998) at 60 kPa housed in a 5°C cold room, simulating conditions at an altitude of 4300 m a.s.l. The cages were cleaned once per week, which required the CH mice to be temporarily exposed to normobaric conditions (<1 h per week). All experiments and protocols approved by the McMaster University Animal Research Ethics Board in accordance with the Canadian Council for Animal Care.

Non-shivering thermogenesis

To assess rates of NST, mice were first removed from their cages and immediately placed in metabolic chambers (475 ml) maintained within their thermoneutral zone at 28°C, using a temperature control cabinet (Sable Systems, Las Vegas, NV, USA). Mice were allowed to adjust to the chambers for 45 min while metabolic rate was monitored by flow-through respirometry in normoxia, as previously described (Robertson and McClelland, 2019). Briefly, air dried and stripped of CO2 (Drierite, W. A. Hammond, Xenia, OH; soda lime and Ascarite, Fisher Scientific, Pittsburgh, PA) was pushed through the chamber at a rate of 1000 ml min−1, and a subsample of the excurrent air was dried before being drawn through CO2 and O2 analyzers (Sable Systems). Mice were then removed from the metabolic chamber and injected subcutaneously with a standardized dose of norepinephrine (NE) based on body mass (Wunder and Gettinger, 1996) in a total of 500 µl saline. We found in preliminary experiments that any increase in metabolism due to injection of saline alone returned to baseline within 10 min. Thus, we quantified NST as the 5 min of the maximal NE-induced O2 while animals were not active (verified using a webcam) and after the effect of saline was complete. We also calculated the increase in O2 above basal metabolic rate (BMR) as NST–BMR (Robertson and McClelland, 2019). BMR for each group was estimated from their average mass using the allometric relationship of body mass and BMR from a separate set of 179 individuals across all three populations in our captive breeding colony, as the lowest O2 in the post-absorptive state, awake and in thermoneutral conditions during 3 h of measurement (N. Wall and G.B.M., unpublished data). To calculate O2 we used the following equations from Withers (1977):
(1)
(2)
where CO2 is production of CO2, i is the flow rate of incurrent air [at standard temperature and pressure (stp)] into the metabolic chamber, and F is the fractional concentration of O2 or CO2 in the incurrent (i) or excurrent (e) air.

Western blotting

Mice were euthanized by an overdose of anesthetic (isoflurane), followed by cervical dislocation to sample iBAT, a major BAT depot in mice (de Jong et al., 2015), which was stored at −80°C before analysis. Western blotting was used to quantify the protein expression of UCP-1, as previously described (Robertson et al., 2019). Unless described otherwise, all reagents were purchased from Millipore Sigma (St. Louis, MO). Frozen iBAT samples were powdered under liquid N2 and homogenized in ice-cold RIPA buffer containing (in mM) 150 NaCl, 1.0% Triton X-100, 0.5% deoxycholic acid, 0.1% SDS, in 50 Tris-HCl, pH 8.0. Total homogenate protein was quantified via Bradford assay (Bradford, 1976) and 15 μg per sample was denatured by heating to 95°C for 5 min in Laemmli sample buffer (Bio-Rad, Mississauga, ON, Canada) with β-mercaptoethanol. Proteins were separated on 12% SDS PAGE gels and then transferred to PVDF membrane. Membranes were probed with primary antibody against UCP-1 (UCP-11-A, Alpha Diagnostics International Inc., San Antonio, TX) at a dilution of 1:500, followed by HRP-conjugated secondary antibody at 1:10,000 (goat anti-rabbit, Santa Cruz Biotechnology, Santa Cruz, CA). Band density was normalized to total lane protein determined by membrane staining with Coomassie Blue, using the Image Lab software package (Bio-Rad). Owing to low availability of P. leucopus in our captive colony, these measurements were only carried out on lowland and highland P. maniculatus.

Mitochondrial isolation

BAT mitochondria were isolated as described previously (Cannon and Nedergaard, 2008; McFarlane et al., 2017; Muleme et al., 2006). For each mitochondrial preparation, iBAT was blunt dissected from one to two individual mice of the same family, cleaned of white adipose tissue and muscle, pooled, and immediately placed in ice-cold homogenization buffer (HB: 250 mmol l−1sucrose, 1 mmol l−1 EGTA, 10 mmol l−1 Hepes, pH 7.4). Tissue was minced and then mechanically homogenized in a glass mortar with a Teflon pestle (6 strokes at 100 rpm) in ice-cold HB. The homogenate was filtered through two layers of cheesecloth and the filtrate was centrifuged at 4°C for 10 min at 8700 g. The pellet was resuspended in HB and centrifuged again at 800 g for 10 min. The supernatant was transferred to a clean tube and was centrifuged again at 8700 g for 10 min. The pellet was resuspended in HB and centrifuged at 8700 g for 10 min twice more. The final pellet was then resuspended in ∼300 µl of HB and kept on ice until further use. Total mitochondrial protein content was quantified by Bradford assay in triplicate.

Mitochondrial respiration

Respiration of isolated iBAT mitochondria was determined by high-resolution respirometry using Clark-type polarographic oxygen electrodes (Oxygraph-2k; Oroboros, Innsbruck Austria) in a total volume of 2 ml respiration buffer [110 mmol l−1 sucrose, 60 mmol l−1 K-lactobionate, 20 mmol l−1 Hepes, 20 mmol l−1 taurine, 10 mmol l−1 KH2PO4, 3 mmol l−1 MgCl2, 0.5 mmol l−1 EGTA and 0.1% (w/v) fatty-acid free BSA, pH 7.1]. Calibration was performed using air-saturated buffer and an oxygen depleted buffer (using a yeast suspension). Measurements were made with constant stirring at 750 rpm and at 37°C. Background respiration through complex I, in the absence of exogenous fatty acid stimulation, was determined using 50 µg of isolated BAT mitochondria in the presence of 1 mmol l−1 pyruvate and 1 mmol l−1 malate. UCP-1-stimulated uncoupled respiration was induced by the addition of 1 mmol l−1 malate and 1 mmol l−1 octanoyl-carnitine (8:0) or 25 µmol l−1 palmitoyl-carnitine (16:0) [free fatty acids (FFAs) of different chain lengths], at a concentration determined in preliminary experiments. The addition of 2 mmol l−1 GDP was used to inhibit UCP-1-mediated respiration (Shabalina et al., 2004), and thus GDP-sensitive respiration (accounting for any residual respiration after GDP addition) is reported as respiration directly attributed to UCP-1 stimulation. Respiration due to oxidative phosphorylation was determined by adding 1.25 mmol l−1 ADP, which was then inhibited with the addition of 1.25 µmol l−1 oligomycin. All respiration measurements were conducted in triplicate, except for measurements with the addition of FFAs, which were measured once per mitochondrial preparation. Isolated mitochondria not used for respiration were stored at −80°C for enzyme assays.

Enzyme assays

Following one round of freeze–thaw, isolated iBAT mitochondria were assayed to determine the apparent Vmax of citrate synthase (CS) and β-hydroxyacyl-CoA dehydrogenase (HOAD). Each assay was performed in triplicate, and an assay without added substrate was performed to correct for any background activity. Assay conditions were as follows for CS: 0.5 mmol l−1 oxaloacetate, 0.15 mmol l−1 acetyl-CoA, 0.05% Triton X-100 and 0.15 mmol l−1 5,5′-dithiobis-(2-nitrobenzoic acid) (DTNB) in 0.04 mmol l−1 Tris-HCl (pH 8.0). For HOAD: 0.15 mmol l−1 acetoacetyl-CoA, 0.15 mmol l−1 NADH in 0.1 mmol l−1 triethanolamine (pH 7.0). Assays were carried out in 96-well format at 37°C using a Spectromax Plus 384 microplate reader (Molecular Devices).

Statistical analysis

Data were analyzed by a two-way analysis of variance (ANOVA), an analysis of covariance (ANCOVA) with body mass as a covariate and population and acclimation as main effects, or a linear mixed model with mass as a covariate and family and/or sex as random variables, in R (https://www.r-project.org/). Multiple comparisons were performed using Tukey HSD or Holm–Šídák post hoc tests. Data are presented as means±s.e.m.

Non-shivering thermogenesis (NST)

We assessed the capacity for NST and the phenotypic plasticity of this important performance variable in lowland deer mice and white-footed mice, compared to highland deer mice with and without acclimation to CH conditions. We defined NST as peak rates of O2 after injection of norepinephrine and found it to be similar between all groups of mice acclimated to control WN conditions (P>0.05). Acclimation to CH led to an increase in NST but only in highlanders, where it was 1.3-fold higher than in warm normoxic conditions (P=0.03). Rates of NST above predicted BMR were also similar across groups kept in WN. There was an effect of CH but only in highlanders (Fig. 1), where NST-BMR increased by 1.7-fold from 1.15±0.21 ml min−1 to 1.92±0.26 ml min−1, a rate significantly higher than that observed in P. leucopus (1.35±0.25 ml min−1) under the same conditions (P=0.04).

Fig. 1.

Norepinephrine-induced non-shivering thermogenesis (NST) and NST-basal metabolic rate of lowland native white-footed mice (Peromyscus leucopus) and deer mice (Peromyscus maniculatus) compared with highland native deer mice acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia. (A) Norepinephrine-induced NST. (B) Norepinephrine-induced NST–BMR (predicted from average body mass for each group based on BMR=5.229×Mb0.689). Sample sizes (mass and predicted BMR) for warm normoxia and cold hypoxia, respectively, were N=9 (20.7±1.5 g, 0.70 ml min−1) and N=4 (16.6±0.7 g, 0.60 ml min−1) for lowland deer mice; N=6 (29.3±1.7, 0.89 ml min−1) and N=6 (26.5±1.9 g, 0.83 ml min−1) for white-footed mice; and N=6 (24.0±3.6 g, 0.78 ml min−1) and N=5 (17.4±0.8 g, 0.62 ml min−1) for highland deer mice. Using mass as a covariate, NST showed a significant population×acclimation interaction (F2,28=4.676, P=0.018). Similarly, for NST–BMR when mass was used as a covariate, the effect of acclimation approached significance (F1,28=3.77, P=0.062), and there was a significant population×acclimation interaction (F2,28=5.32, P=0.011). Bars that do not share a letter are significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

Fig. 1.

Norepinephrine-induced non-shivering thermogenesis (NST) and NST-basal metabolic rate of lowland native white-footed mice (Peromyscus leucopus) and deer mice (Peromyscus maniculatus) compared with highland native deer mice acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia. (A) Norepinephrine-induced NST. (B) Norepinephrine-induced NST–BMR (predicted from average body mass for each group based on BMR=5.229×Mb0.689). Sample sizes (mass and predicted BMR) for warm normoxia and cold hypoxia, respectively, were N=9 (20.7±1.5 g, 0.70 ml min−1) and N=4 (16.6±0.7 g, 0.60 ml min−1) for lowland deer mice; N=6 (29.3±1.7, 0.89 ml min−1) and N=6 (26.5±1.9 g, 0.83 ml min−1) for white-footed mice; and N=6 (24.0±3.6 g, 0.78 ml min−1) and N=5 (17.4±0.8 g, 0.62 ml min−1) for highland deer mice. Using mass as a covariate, NST showed a significant population×acclimation interaction (F2,28=4.676, P=0.018). Similarly, for NST–BMR when mass was used as a covariate, the effect of acclimation approached significance (F1,28=3.77, P=0.062), and there was a significant population×acclimation interaction (F2,28=5.32, P=0.011). Bars that do not share a letter are significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

iBAT mass and UCP-1 protein expression

To investigate the underlying traits responsible for changes in NST with acclimation, we assessed iBAT mass and expression of UCP-1. Changes in iBAT mass and/or UCP-1 expression have both been shown to correlate with changes in NST in other studies (McDevitt and Andrews, 1997; Beaudry and McClelland, 2010). Since iBAT mass varies as a function of body mass in other mammals (Aleksiuk and Frohlinger, 1971; Król and Speakman, 2019; Smith and Roberts, 1964), we first determined this relationship for deer mice and white-footed mice. We found that iBAT mass varied directly with body mass across populations and with acclimation conditions (Fig. 2A). After accounting for any differences in body mass, we found that iBAT mass tended to be lower in white-footed mice than lowland deer mice, but this difference did not reach statistical significance with multiple comparisons test (P=0.06). However, there was no significant effect of CH acclimation on iBAT mass in any of the experimental groups. In fact, all mice maintained consistent iBAT mass across the two acclimation conditions (Fig. 2B).

Fig. 2.

Mass of interscapular brown adipose tissue (iBAT) in lowland P. leucopus, and in lowland and highland P. maniculatus. (A) Mass (Mb) of iBAT varied with body mass as BATmass=0.0052×Mb0.997. (B) Mass of iBAT in lowland and highland native mice, all acclimated to 6–8 weeks of either warm normoxia (control conditions) or cold hypoxia. Samples sizes (and average body mass) for warm normoxia and cold hypoxia, respectively, were N=9 (20.7±1.8 g) and N=4 (17.3±1.9 g) for lowland deer mice; N=10 (26.2±1.8 g) and N=12 (26.2±1.4 g) for white-footed mice; and N=16 (21.0±1.3 g) and N=12 (19.3±1.0 g) for highland deer mice. Using mass as a covariate, there was a significant effect of population (F2,55=4.802, P=0.012) but not of acclimation (F1,55=2.589, P=0.113) nor a significant interaction (F2,55=2.53, P=0.089). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values in B are presented as means±s.e.m.

Fig. 2.

Mass of interscapular brown adipose tissue (iBAT) in lowland P. leucopus, and in lowland and highland P. maniculatus. (A) Mass (Mb) of iBAT varied with body mass as BATmass=0.0052×Mb0.997. (B) Mass of iBAT in lowland and highland native mice, all acclimated to 6–8 weeks of either warm normoxia (control conditions) or cold hypoxia. Samples sizes (and average body mass) for warm normoxia and cold hypoxia, respectively, were N=9 (20.7±1.8 g) and N=4 (17.3±1.9 g) for lowland deer mice; N=10 (26.2±1.8 g) and N=12 (26.2±1.4 g) for white-footed mice; and N=16 (21.0±1.3 g) and N=12 (19.3±1.0 g) for highland deer mice. Using mass as a covariate, there was a significant effect of population (F2,55=4.802, P=0.012) but not of acclimation (F1,55=2.589, P=0.113) nor a significant interaction (F2,55=2.53, P=0.089). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values in B are presented as means±s.e.m.

In contrast, CH acclimation led to an increased expression of UCP-1 protein (Fig. 3), but only in highland deer mice. Compared with highlanders kept in WN conditions, those acclimated to CH showed a 1.9-fold increase in UCP-1 expression (P<0.01), to levels seen in CH-acclimated lowlanders.

Fig. 3.

Relative protein expression of uncoupling protein (UCP)-1 in brown adipose tissue of lowland compared with highland native deer mice acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia. Representative western blot (A) and plotted protein abundance (B) of UCP-1. Sample sizes for warm normoxia and cold hypoxia, respectively, were N=5 and N=6 for lowland deer mice, and N=5 and N=4 for highland deer mice. There was a significant effect of acclimation (F1,16=8.435, P=0.010) and a significant population×acclimation interaction (F1,16=5.525, P=0.032). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

Fig. 3.

Relative protein expression of uncoupling protein (UCP)-1 in brown adipose tissue of lowland compared with highland native deer mice acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia. Representative western blot (A) and plotted protein abundance (B) of UCP-1. Sample sizes for warm normoxia and cold hypoxia, respectively, were N=5 and N=6 for lowland deer mice, and N=5 and N=4 for highland deer mice. There was a significant effect of acclimation (F1,16=8.435, P=0.010) and a significant population×acclimation interaction (F1,16=5.525, P=0.032). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

Isolated mitochondria respiration and enzyme activities

To understand if the metabolic properties of mitochondria might contribute to differences in NST, we isolated mitochondria from iBAT and determined basal and UCP-1-stimulated rates of respiration. When we assessed respiration using pyruvate and malate (in the absence of exogenous ADP or FFAs), WN acclimated lowland deer mice showed higher rates of this background respiration than both white-footed mice (P=0.03) and highland deer mice (P=0.01). In general, lowland deer mice also had higher rates of background respiration than all other groups after CH acclimation (Table 1). The exception was white-footed mice, which showed an increase in background respiration with CH acclimation (P<0.001) to rates that were equivalent to both lowland and highland deer mice. When respiration was further stimulated using exogenous FFAs, which are known to activate UCP-1, rates increased in all WN acclimated mice, but respiration was higher in lowland deer mice compared with either white-footed mice (P=0.006) or highland deer mice (P=0.016). In contrast, only highland deer mice showed an increase in FFA-simulated respiration after CH acclimation, although the increase was only marginally significant compared with levels in the control (P=0.054). Nevertheless, all CH acclimated deer mice and white footed mice had similar FFA-stimulated respiration rates. FFA-stimulated respiration was strongly inhibited by the addition of GDP, by 4- to 5-fold in all groups. Respiration that was directly attributed to UCP-1 (i.e. GDP-sensitive palmitoyl-carnitine respiration) was compared between WN acclimated mice and showed no significant difference between populations. However, both highland deer mice (P=0.062) and white-footed mice (P=0.07) showed a strong trend to increase respiration, by 1.6-fold and 1.8-fold, respectively, after CH acclimation (Fig. 4). Indeed, after CH acclimation, highland mice had higher rates of respiration than white-footed mice in warm normoxia (P=0.018). Regardless, all groups showed similar rates of mitochondrial respiration after CH acclimation. We also found that acyl chain length had little impact on FFA-stimulated respiration, with similar rates of UCP-1-dependent respiration with octanoyl-carnitine to those with palmitoyl-carnitine. When oxidative phosphorylation was induced by the addition of ADP, rates of respiration were relatively low, and only ∼20–40% the rate of UCP-1-stimulated respiration (Table 1).

Fig. 4.

UCP-1-stimulated respiration in isolated mitochondria of interscapular brown adipose tissue of P. leucopus and P. maniculatus. UCP-1-stimulated respiration, measured as palmitoyl-carnitine-stimulated respiration that was GDP sensitive (in nmol O2 min−1 mg−1) in mitochondria isolated from interscapular brown adipose tissue from mice acclimated to 6–8 weeks of either warm normoxia (control conditions) or cold hypoxia. Sample sizes (each sample represents BAT pooled from N=2 mice of the same family) for warm normoxia and cold hypoxia, respectively, were N=8 and N=6 for lowland deer mice, N=5 and N=10 for white-footed mice, and N=8 and N=8 for highland deer mice. There was a significant effect of acclimation (F1,39=17.17, P<0.001) and an effect of population that approached statistical significance (F2,39=2.99, P=0.062). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

Fig. 4.

UCP-1-stimulated respiration in isolated mitochondria of interscapular brown adipose tissue of P. leucopus and P. maniculatus. UCP-1-stimulated respiration, measured as palmitoyl-carnitine-stimulated respiration that was GDP sensitive (in nmol O2 min−1 mg−1) in mitochondria isolated from interscapular brown adipose tissue from mice acclimated to 6–8 weeks of either warm normoxia (control conditions) or cold hypoxia. Sample sizes (each sample represents BAT pooled from N=2 mice of the same family) for warm normoxia and cold hypoxia, respectively, were N=8 and N=6 for lowland deer mice, N=5 and N=10 for white-footed mice, and N=8 and N=8 for highland deer mice. There was a significant effect of acclimation (F1,39=17.17, P<0.001) and an effect of population that approached statistical significance (F2,39=2.99, P=0.062). Bars that do not share a letter were significantly different from each other as determined by multiple comparisons (Tukey). Values presented as means±s.e.m.

Table 1.

Respiration rates (in nmol O2 min−1 mg−1 protein) of isolated mitochondria from brown adipose tissue of first generation lab born and raised lowland native deer mice (P. maniculatus) and white-footed mice (P. leucopus) compared with highland native deer mice, all acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia

Respiration rates (in nmol O2 min−1 mg−1 protein) of isolated mitochondria from brown adipose tissue of first generation lab born and raised lowland native deer mice (P. maniculatus) and white-footed mice (P. leucopus) compared with highland native deer mice, all acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia
Respiration rates (in nmol O2 min−1 mg−1 protein) of isolated mitochondria from brown adipose tissue of first generation lab born and raised lowland native deer mice (P. maniculatus) and white-footed mice (P. leucopus) compared with highland native deer mice, all acclimated to 6-8 weeks of either warm normoxia (control conditions) or cold hypoxia

Differences in mitochondrial respiration were not related to mitochondrial-specific activities (apparent Vmax) of CS or HOAD. In fact, CS activity tended to decline with CH acclimation (significant acclimation effect), likely driven by reductions in activity in the two LA species. In contrast, mitochondrial HOAD activity showed little change in activity with acclimation (Fig. 5).

Fig. 5.

Apparent Vmax of citrate synthase (CS) and β-hydroxyacyl-CoA dehydrogenase (HOAD) of isolated mitochondria from P. leucopus and P. maniculatus. Enzyme activity of (A) CS and (B) HOAD in isolated mitochondria. Sample sizes for warm normoxia and cold hypoxia, respectively, were N=3 and N=3 for lowland deer mice, N=9 and N=5 for white-footed mice, and N=7 and N=10 for highland deer mice. For HOAD there was an effect of acclimation that approached statistical significance (F1,28=3.726, P=0.064). For CS there was a significant effect of acclimation (*F1,28=3.298, P=0.008), and a nearly significant population×acclimation interaction (F2,28=2.838, P=0.076). Values presented as means±s.e.m.

Fig. 5.

Apparent Vmax of citrate synthase (CS) and β-hydroxyacyl-CoA dehydrogenase (HOAD) of isolated mitochondria from P. leucopus and P. maniculatus. Enzyme activity of (A) CS and (B) HOAD in isolated mitochondria. Sample sizes for warm normoxia and cold hypoxia, respectively, were N=3 and N=3 for lowland deer mice, N=9 and N=5 for white-footed mice, and N=7 and N=10 for highland deer mice. For HOAD there was an effect of acclimation that approached statistical significance (F1,28=3.726, P=0.064). For CS there was a significant effect of acclimation (*F1,28=3.298, P=0.008), and a nearly significant population×acclimation interaction (F2,28=2.838, P=0.076). Values presented as means±s.e.m.

The main objective of this study was to determine if highland native deer mice have evolved an enhanced plasticity in NST with acclimation to cold hypoxia (CH) conditions that simulate their native altitude. Our data support this hypothesis, with highland mice showing a significant increase in NST after acclimation. In contrast, lowland conspecifics and congenerics showed no change in NST when acclimated to the same conditions. Variation in NST with CH acclimation was not driven by changes in iBAT mass, which remained consistent across experimental groups and acclimation conditions. However, acclimation to CH in highland deer mice led to an increase in the expression of UCP-1 in iBAT and in UCP-1-stimulated respiration of isolated mitochondria, changes not seen in lowlanders.

We found that high altitude native mice increased NST by 32% with acclimation to conditions that simulate an elevation of 4300 m, while low altitude native mice showed little change in this trait in response to the same conditions. Indeed, the response to CH acclimation was even greater for NST-BMR, which increased by 67% from WN in highlanders. This suggests that environment-induced plasticity plays an important role in elevated NST seen in wild highland deer mice measured in their native environment, when compared with lowland mice measured at low altitude (Velotta et al., 2016). This population difference likely reflects the influence of the persistent low temperatures experienced by the highland mice, even in summer (Hayes, 1989; Velotta et al., 2016). Indeed, descendants of lowland and highland deer mice born and raised in common laboratory conditions had equivalent NST capacities prior to acclimation (Fig. 1; Velotta et al., 2016). It is well recognized that cold acclimation increases NST capacity by ∼40% in lowland mammals, including laboratory mice and deer mice (Beaudry and McClelland, 2010; Oufara et al., 1987; Van Sant and Hammond, 2008). However, the addition of chronic hypoxia exposure can blunt the effect of cold on thermogenic capacity and NST (Beaudry and McClelland, 2010). Our data show that highland deer mice have overcome the antagonistic effect that hypoxia has on the plasticity of BAT-based heat production, with a distinct response to the combination of stressors experienced at high elevations. Here, we show that iBAT also undergoes phenotypic plasticity in response to CH with increased rates of NST that would enable highland mice to cope with the cold conditions at altitude.

Increases in NST capacity seen in highland mice with CH were not associated with an increase in iBAT mass, in contrast to changes in the size of this tissue that other species experience with seasonal changes in temperature (e.g. Apodemus flavicollis, Klaus et al., 1988). Acclimation to cold alone is known to induce BAT growth in rodents (Bukowiecki et al., 1982; Nedergaard et al., 1980; Smith and Roberts, 1964), but when combined with hypoxia we found any potential growth was inhibited in deer mice and white-footed mice (Fig. 2). This is consistent with data from laboratory mice, where acclimation to equivalent CH conditions did not influence absolute iBAT mass (Beaudry and McClelland, 2010). Limiting BAT growth may be beneficial in hypoxic environments to minimize the metabolic costs associated with maintaining a large and metabolically active tissue. Instead, what may underlie the CH-stimulated increase in NST are changes in the capacity of brown adipocytes for heat production. Consistent with this idea, wild highland deer mice show 1.5-fold higher Ucp1 gene expression in iBAT compared with lab-reared descendants (Velotta et al., 2016). Furthermore, we found that lab-reared highland deer mice increased UCP-1 protein expression and had greater UCP-1-stimulated mitochondrial respiration after CH acclimation. The benefit of increased UCP-1 expression is that it provides a higher capacity for uncoupled respiratory capacity, without affecting basal proton leak (Shabalina et al., 2010). Interestingly, basal respiration, presumably to support proton leak independent of UCP-1 stimulation, tended to be lower in highland deer mice compared with lowland deer mice (Table 1). Therefore, iBAT thermogenic capacity is increased to support higher NST rates, with limited effect on energy expenditure when mice are not actively thermoregulating. These changes, along with increased expression of genes involved in BAT vascularization (Velotta et al., 2016) and an enhanced systemic blood flow (Tate et al., 2017, 2020), contribute to the increased capacity for NST seen in CH acclimated highlanders.

The environmentally induced plasticity of NST at high altitude likely contributes to the high overall thermogenic capacity seen in highland deer mice (Cheviron et al., 2012, 2014). NST is known to contribute >50% of total thermogenic capacity in rodents (e.g. Phyllotis xanthopugus, Nespolo et al., 1999; P. maniculatus, Van Sant and Hammond, 2008). The rates of NST for WN acclimated deer mice and white- footed mice contributed 35–57% of previously reported thermogenic capacity determined in hypoxia (Tate et al., 2020). After CH acclimation, lowland mice generally showed a decline in the contribution of NST to total thermogenesis, but highland deer mice showed an increase in the relative contribution of NST. The adaptive significance of increased NST reliance with CH acclimation in highland deer mice is presently unclear. While NST may be a more efficient method of heat production than shivering, highland deer mice have a more oxidative muscle phenotype relative to lowlanders (Lui et al., 2015), presumably for high rates of shivering. However, CH acclimation does not influence the aerobic phenotype of gastrocnemius muscle in highland deer mice (Mahalingam et al., 2020), suggesting that the capacity for shivering remains unchanged. While high alpine environments are colder than low altitude sites of the same latitude, they still experience seasonal and daily temperature variation. Thus, the evolved plasticity of NST in response to CH may be essential for further increases in thermogenic capacity during the coldest parts of the year at high altitude.

Similarly to other rodent species, our lowland native Peromyscus spp. showed no increase in NST with CH (e.g. Beaudry and McClelland, 2010), suggesting this may be a conserved response for lowland native mice. The mechanistic underpinnings of the antagonistic effect of hypoxia on cold in lowlanders is unclear, although there are several possibilities. The hypothalamus integrates ambient and body temperature information, activating both shivering and NST when required. Hypoxia tends to lower the body temperature setpoint in the hypothalamus of rodents, making them less responsive to cold temperatures (Tattersall and Milsom, 2009). Hypoxia has also been shown to inhibit adrenergic signaling in lowland rats (Voelkel et al., 1981), which is a requisite for BAT recruitment (Cannon and Nedergaard, 2004). Together, these data may explain the resistance of NST and iBAT phenotype to CH acclimation in lowland natives if BAT activation is generally inhibited. Previous studies on adrenergic signaling in hypoxia have focused primarily on cardiac function (León-Velarde et al., 2001; Voelkel et al., 1981) and the effects of hypoxia on adrenergic signaling in BAT during hypoxia have not been explored. Moreover, active BAT requires appropriate O2 and fuel delivery (Bartelt et al., 2011; Heim and Hull, 1966). Highland deer mice maintain higher arterial O2 saturation than lowlanders in hypoxia (Tate et al., 2017), suggesting that highlanders can deliver sufficient O2 to activated BAT in hypoxia while lowlanders may not. Together, these traits likely allow for a robust environmentally induced increase in NST and thermogenic capacity necessary to survive the harsh high alpine conditions.

Conclusions

Our results suggest that evolved plasticity of thermogenic capacity is important for survival in adult deer mice at high altitude, consistent with previous reports (Cheviron et al., 2014; Tate et al., 2020; Velotta et al., 2016). Here, we show that highland deer mice have an enhanced plasticity to cold hypoxia that increases the capacity for NST. Skeletal muscles of highland deer mice show constitutively higher aerobic capacities compared with lowland mice, and this aerobic phenotype is resistant to hypoxia (Lui et al., 2015; Mahalingam et al., 2017) and CH (Mahalingam et al., 2020) acclimation. This suggests that plasticity of whole animal thermogenic capacity in highland mice is principally the result of changes in NST rather than shivering. Our data also suggest that it is an evolved plasticity of iBAT function that contributes to enhanced rates of NST-based heat production at high altitude.

In this study, we have found that highland deer mice increase NST capacity with CH acclimation by functionally remodeling iBAT mitochondria. However, the convergence in iBAT phenotype between mice from different altitudes after CH acclimation suggests that iBAT phenotype cannot completely explain changes in whole-animal thermogenesis. Instead, differences may lie at the integrated whole-tissue or whole-animal level. A greater capacity for tissue O2 delivery and blood O2 carrying capacity of highland deer mice in hypoxia may underly some of these mechanisms. The phenotypic plasticity of BAT and NST observed in high altitude deer mice contributes to high thermogenic capacity, which allows them to be active in high altitude conditions and invade an ecological niche inaccessible to other mice. Future areas of study to determine the proximate mechanism behind this enhanced plasticity include comparisons of thermogenic function in intact brown adipocytes between acclimations, fuel provisioning to BAT mitochondria in vivo and contributions to NST from other depots of BAT and other tissues, such as beige fat.

The authors thank N. Dawson for technical assistance with enzyme assays and G. Scott for use of his Oroboros Oxygraph-2k.

Author contributions

Conceptualization: G.B.M.; Methodology: S.Z.C., C.E.R., S.M.; Formal analysis: S.Z.C., G.B.M.; Resources: G.B.M.; Data curation: S.Z.C., C.E.R., S.M.; Writing - original draft: S.Z.C., G.B.M.; Writing - review & editing: S.Z.C., C.E.R., S.M., G.B.M.; Supervision: G.B.M.; Funding acquisition: G.B.M.

Funding

This work was funded by a Natural Science and Engineering Research Council of Canada (NSERC) Discovery Grant, an NSERC Discovery Accelerator Supplement (G.B.M.), an NSERC Canada Graduate Scholarship (C.E.R.) and an Ontario Graduate Scholarship (S.M.).

Data availability

Data are available from figshare: https://doi.org/10.6084/m9.figshare.14357297.v1

Aleksiuk
,
M.
and
Frohlinger
,
A.
(
1971
).
Seasonal metabolic organization in the muskrat (Ondatra zibethica). I. Changes in growth, thyroid activity, brown adipose tissue, and organ weights in nature
.
Can. J. Zool.
49
,
1143
-
1154
.
Ballinger
,
M. A.
,
Hess
,
C.
,
Napolitano
,
M. W.
,
Bjork
,
J. A.
and
Andrews
,
M. T.
(
2016
).
Seasonal changes in brown adipose tissue mitochondria in a mammalian hibernator: from gene expression to function
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
311
,
R325
-
R336
.
Bartelt
,
A.
,
Bruns
,
O. T.
,
Reimer
,
R.
,
Hohenberg
,
H.
,
Ittrich
,
H.
,
Peldschus
,
K.
,
Kaul
,
M. G.
,
Tromsdorf
,
U. I.
,
Weller
,
H.
,
Waurisch
,
C.
et al. 
(
2011
).
Brown adipose tissue activity controls triglyceride clearance
.
Nat. Med.
17
,
200
-
205
.
Beaudry
,
J. L.
and
McClelland
,
G. B.
(
2010
).
Thermogenesis in CD-1 mice after combined chronic hypoxia and cold acclimation
.
Comp. Biochem. Physiol. B Biochem. Mol. Biol.
157
,
301
-
309
.
Böckler
,
H.
and
Heldmaier
,
G.
(
1983
).
Interaction of shivering and non-shivering thermogenesis during cold exposure in seasonally-acclimatized Djungarian hamsters (Phodopus sungorus)
.
J. Therm. Biol.
8
,
97
-
98
.
Bradford
,
M. M.
(
1976
).
A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding
.
Anal. Biochem.
72
,
248
-
254
.
Bukowiecki
,
L.
,
Collet
,
A. J.
and
Follea
,
N.
,
Guay
,
G.
and
Jahjah
,
L.
(
1982
).
Brown adipose tissue hyperplasia: a fundamental mechanism of adaptation to cold and hyperphagia
.
Am. J. Physiol. Endocrinol. Metab.
242
,
E353
-
E359
.
Cannon
,
B.
and
Nedergaard
,
J.
(
2004
).
Brown adipose tissue: function and physiological significance
.
Physiol. Rev.
84
,
277
-
359
.
Cannon
,
B.
and
Nedergaard
,
J.
(
2008
).
Studies of thermogenesis and mitochondrial function in adipose tissues
. In
Adipose Tissue Protocols
(ed.
K.
Yang
), pp.
109
-
121
.
Totowa
,
USA
:
Humana Press
.
Cannon
,
B.
and
Nedergaard
,
J.
(
2011
).
Nonshivering thermogenesis and its adequate measurement in metabolic studies
.
J. Exp. Biol.
214
,
242
-
253
.
Cheviron
,
Z. A.
,
Bachman
,
G. C.
,
Connaty
,
A. D.
,
McClelland
,
G. B.
and
Storz
,
J. F.
(
2012
).
Regulatory changes contribute to the adaptive enhancement of thermogenic capacity in high-altitude deer mice
.
Proc. Natl. Acad. Sci. USA
109
,
8635
-
8640
.
Cheviron
,
Z. A.
,
Connaty
,
A. D.
,
McClelland
,
G. B.
and
Storz
,
J. F.
(
2014
).
Functional genomics of adaptation to hypoxic cold-stress in high-altitude deer mice: transcriptomic plasticity and thermogenic performance
.
Evolution
68
,
48
-
62
.
de Jong
,
J. M. A.
,
Larsson
,
O.
,
Cannon
,
B.
and
Nedergaard
,
J.
(
2015
).
A stringent validation of mouse adipose tissue identity markers
.
Am. J. Physiol. Endocrinol. Metab.
308
,
E1085
-
E1105
.
Fedorenko
,
A.
,
Lishko
,
P. V.
and
Kirichok
,
Y.
(
2012
).
Mechanism of fatty-acid-dependent UCP1 uncoupling in brown fat mitochondria
.
Cell
151
,
400
-
413
.
Hayes
,
J. P.
(
1989
).
Field and maximal metabolic rates of deer mice (Peromyscus maniculatus) at low and high altitudes
.
Physiol. Zool.
62
,
732
-
744
.
Hayes
,
J. P.
and
O'Connor
,
C. S.
(
1999
).
Natural selection on thermogenic capacity of high-altitude deer mice
.
Evolution
53
,
1280
-
1287
.
Heim
,
T.
and
Hull
,
D.
(
1966
).
The blood flow and oxygen consumption of brown adipose tissue in the new-born rabbit
.
J. Physiol.
186
,
42
-
55
.
Heim
,
A. B.
,
Chung
,
D.
,
Florant
,
G. L.
and
Chicco
,
A. J.
(
2017
).
Tissue-specific seasonal changes in mitochondrial function of a mammalian hibernator
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
313
,
R180
-
R190
.
Heldmaier
,
G.
(
1971
).
Zitterfreie Wärmebildung und Körpergröße bei Säugetieren
.
Z. Vergl. Physiologie
73
,
222
-
247
.
Heldmaier
,
G.
,
Klaus
,
S.
,
Wiesinger
,
U. F.
and
Wenzel
,
M.
(
1989
).
Cold acclimation and thermogenesis
. In
Living in the Cold II
(ed.
A.
Malan
and
B.
Canguilhem
), pp.
347
-
358
.
Cobham
,
UK
:
John Libbey Eurotext Ltd
.
Himms-Hagen
,
J.
(
1985
).
Brown adipose tissue metabolism and thermogenesis
.
Annu. Rev. Nutr.
5
,
69
-
94
.
Janský
,
L.
(
1973
).
Non-shivering thermogenesis and its thermoregulatory significance
.
Biol. Rev.
48
,
85
-
132
.
Klaus
,
S.
,
Heldmaier
,
G.
and
Ricquier
,
D.
(
1988
).
Seasonal acclimation of bank voles and wood mice: nonshivering thermogenesis and thermogenic properties of brown adipose tissue mitochondria
.
J. Comp. Physiol. B
158
,
157
-
164
.
Klingenspor
,
M.
,
Ebbinghaus
,
C.
,
Hülshorst
,
G.
,
Stöhr
,
S.
,
Spiegelhalter
,
F.
,
Haas
,
K.
and
Heldmaier
,
G.
(
1996
).
Multiple regulatory steps are involved in the control of lipoprotein lipase activity in brown adipose tissue
.
J. Lipid Res.
37
,
1685
-
1695
.
Król
,
E.
and
Speakman
,
J. R.
(
2019
).
Switching off the furnace: brown adipose tissue and lactation
.
Mol. Asp. Med.
68
,
18
-
41
.
Lacy
,
R. C.
,
Lynch
,
C. B.
and
Lynch
,
G. R.
(
1978
).
Developmental and adult acclimation effects of ambient temperature on temperature regulation of mice selected for high and low levels of nest-building
.
J. Comp. Physiol.
123
,
185
-
192
.
Lau
,
D. S.
,
Connaty
,
A. D.
,
Mahalingam
,
S.
,
Wall
,
N.
,
Cheviron
,
Z. A.
,
Storz
,
J. F.
,
Scott
,
G. R.
and
McClelland
,
G. B.
(
2017
).
Acclimation to hypoxia increases carbohydrate use during exercise in high-altitude deer mice
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
312
,
R400
-
R411
.
León-Velarde
,
F.
,
Bourin
,
M.-C.
,
Germack
,
R.
,
Mohammadi
,
K.
,
Crozatier
,
B.
and
Richalet
,
J.-P.
(
2001
).
Differential alterations in cardiac adrenergic signaling in chronic hypoxia or norepinephrine infusion
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
280
,
R274
-
R281
.
Lui
,
M. A.
,
Mahalingam
,
S.
,
Patel
,
P.
,
Connaty
,
A. D.
,
Ivy
,
C. M.
,
Cheviron
,
Z. A.
,
Storz
,
J. F.
,
McClelland
,
G. B.
and
Scott
,
G. R.
(
2015
).
High-altitude ancestry and hypoxia acclimation have distinct effects on exercise capacity and muscle phenotype in deer mice
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
308
,
R779
-
R791
.
Mahalingam
,
S.
,
McClelland
,
G. B.
and
Scott
,
G. R.
(
2017
).
Evolved changes in the intracellular distribution and physiology of muscle mitochondria in high–altitude native deer mice
.
J. Physiol.
595
,
4785
-
4801
.
Mahalingam
,
S.
,
Cheviron
,
Z. A.
,
Storz
,
J. F.
,
McClelland
,
G. B.
and
Scott
,
G. R.
(
2020
).
Chronic cold exposure induces mitochondrial plasticity in deer mice native to high altitudes
.
J. Physiol.
598
,
5411
-
5426
.
McClelland
,
G. B.
,
Hochachka
,
P. W.
and
Weber
,
J.-M.
(
1998
).
Carbohydrate utilization during exercise after high-altitude acclimation: a new perspective
.
Proc. Natl. Acad. Sci. USA
95
,
10288
-
10293
.
McDevitt
,
R. M.
and
Andrews
,
J. F.
(
1997
).
Seasonal variation in brown adipose tissue mass and lipid droplet size of Sorex minutus, the pygmy shrew; the relationship between morphology and metabolic rate
.
J. Therm. Biol.
22
,
127
-
135
.
McFarlane
,
S. V.
,
Mathers
,
K. E.
and
Staples
,
J. F.
(
2017
).
Reversible temperature-dependent differences in brown adipose tissue respiration during torpor in a mammalian hibernator
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
312
,
R434
-
R442
.
Muleme
,
H. M.
,
Walpole
,
A. C.
and
Staples
,
J. F.
(
2006
).
Mitochondrial metabolism in hibernation: metabolic suppression, temperature effects, and substrate preferences
.
Physiol. Biochem. Zool.
79
,
474
-
483
.
Nedergaard
,
J.
,
Alexson
,
S.
and
Cannon
,
B.
(
1980
).
Cold adaptation in the rat: increased brown fat peroxisomal beta-oxidation relative to maximal mitochondrial oxidative capacity
.
Am. J. Physiol. Cell Physiol.
239
,
C208
-
C216
.
Nespolo
,
R. F.
,
Opazo
,
J. C.
,
Rosenmann
,
M.
and
Bozinovic
,
F.
(
1999
).
Thermal acclimation, maximum metabolic rate, and nonshivering thermogenesis of Phyllotis xanthopygus (Rodentia) in the Andes mountains
.
J. Mammal.
80
,
742
-
748
.
Oufara
,
S.
,
Barre
,
H.
,
Rouanet
,
J. L.
and
Chatonnet
,
J.
(
1987
).
Adaptation to extreme ambient temperatures in cold-acclimated gerbils and mice
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
253
,
R39
-
R45
.
Robertson
,
C. E.
and
McClelland
,
G. B.
(
2019
).
Developmental delay in shivering limits thermogenic capacity in juvenile high-altitude deer mice (Peromyscus maniculatus)
.
J. Exp. Biol.
222
,
jeb210963
.
Robertson
,
C. E.
,
Tattersall
,
G. J.
and
McClelland
,
G. B.
(
2019
).
Development of homeothermic endothermy is delayed in high-altitude native deer mice (Peromyscus maniculatus)
.
Proc. R. Soc. B Biol. Sci.
286
,
20190841
.
Scott
,
G. R.
,
Elogio
,
T. S.
,
Lui
,
M. A.
,
Storz
,
J. F.
and
Cheviron
,
Z. A.
(
2015
).
Adaptive modifications of muscle phenotype in high-altitude deer mice are associated with evolved changes in gene regulation
.
Mol. Biol. Evol.
32
,
1962
-
1976
.
Sears
,
M. W.
,
Hayes
,
J. P.
,
O'Conner
,
C. S.
,
Geluso
,
K.
and
Sedinger
,
J. S.
(
2006
).
Individual variation in thermogenic capacity affects above-ground activity of high-altitude deer mice
.
Funct. Ecol.
20
,
97
-
104
.
Shabalina
,
I. G.
,
Jacobsson
,
A.
,
Cannon
,
B.
and
Nedergaard
,
J.
(
2004
).
Native UCP1 displays simple competitive kinetics between the regulators purine nucleotides and fatty acids
.
J. Biol. Chem.
279
,
38236
-
38248
.
Shabalina
,
I. G.
,
Ost
,
M.
,
Petrovic
,
N.
,
Vrbacky
,
M.
,
Nedergaard
,
J.
and
Cannon
,
B.
(
2010
).
Uncoupling protein-1 is not leaky
.
Biochim. Biophys. Acta (BBA) Bioenerg.
1797
,
773
-
784
.
Smith
,
R. E.
and
Roberts
,
J. C.
(
1964
).
Thermogenesis of brown adipose tissue in cold-acclimated rats
.
Am. J. Physiol. Legacy Content
206
,
143
-
148
.
Tate
,
K. B.
,
Ivy
,
C. M.
,
Velotta
,
J. P.
,
Storz
,
J. F.
,
McClelland
,
G. B.
,
Cheviron
,
Z. A.
and
Scott
,
G. R.
(
2017
).
Circulatory mechanisms underlying adaptive increases in thermogenic capacity in high-altitude deer mice
.
J. Exp. Biol.
220
,
3616
-
3620
.
Tate
,
K. B.
,
Wearing
,
O. H.
,
Ivy
,
C. M.
,
Cheviron
,
Z. A.
,
Storz
,
J. F.
,
McClelland
,
G. B.
and
Scott
,
G. R.
(
2020
).
Coordinated changes across the O2 transport pathway underlie adaptive increases in thermogenic capacity in high-altitude deer mice
.
Proc. Biol. Sci. B Biol. Sci.
287
,
20192750
.
Tattersall
,
G. J.
and
Milsom
,
W. K.
(
2009
).
Hypoxia reduces the hypothalamic thermogenic threshold and thermosensitivity
.
J. Physiol.
587
,
5259
-
5274
.
Van Sant
,
M. J.
and
Hammond
,
K. A.
(
2008
).
Contribution of shivering and nonshivering thermogenesis to thermogenic capacity for the deer mouse (Peromyscus maniculatus)
.
Physiol. Biochem. Zool.
81
,
605
-
611
.
Velotta
,
J. P.
,
Jones
,
J.
,
Wolf
,
C. J.
and
Cheviron
,
Z. A.
(
2016
).
Transcriptomic plasticity in brown adipose tissue contributes to an enhanced capacity for nonshivering thermogenesis in deer mice
.
Mol. Ecol.
25
,
2870
-
2886
.
Voelkel
,
N. F.
,
Hegstrand
,
L.
,
Reeves
,
J. T.
,
McMurty
,
I. F.
and
Molinoff
,
P. B.
(
1981
).
Effects of hypoxia on density of beta-adrenergic receptors
.
J. Appl. Physiol.
50
,
363
-
366
.
Wickler
,
S. J.
(
1980
).
Maximal thermogenic capacity and body temperatures of white-footed mice (Peromyscus) in summer and winter
.
Physiol. Zool.
53
,
338
-
346
.
Withers
,
P. C.
(
1977
).
Measurement of VO2, VCO2, and evaporative water loss with a flow-through mask
.
J. Appl. Physiol.
42
,
120
-
123
.
Wunder
,
B. A.
and
Gettinger
,
R. D.
(
1996
).
Effects of body mass and temperature acclimation on the nonshivering thermogenic response of small mammals
. In
Adaptations to the Cold: Tenth International Hibernation Symposium
(ed.
F.
Geiser
,
A. J.
Hulbert
and
S. C.
Nicol
), pp.
131
-
139
.
Armidale
,
AUS
:
University of New England Press
.
Yu
,
X. X.
,
Lewin
,
D. A.
,
Forrest
,
W.
and
Adams
,
S. H.
(
2002
).
Cold elicits the simultaneous induction of fatty acid synthesis and β-oxidation in murine brown adipose tissue: prediction from differential gene expression and confirmation in vivo
.
FASEB J.
16
,
155
-
168
.

Competing interests

The authors declare no competing or financial interests.