Abstract Hibernation is a state of extraordinary metabolic plasticity. The pathways of amino acid metabolism as they relate to nitrogen homeostasis in hibernating mammals in vivo is unknown. Here we show, using pulse isotopic tracing, evidence of increased myofibrillar (skeletal muscle) protein breakdown and suppressed whole body production (WBP) of metabolites in vivo throughout deep torpor. As WBP of metabolites is suppressed, amino acids with nitrogenous side chains accumulate over torpor while production of urea cycle intermediates do not. Using ^15N stable isotope methodology in arctic ground squirrels (Urocitellus parryii), we provide evidence that free nitrogen is buffered and recycled into essential amino acids (EAA), non-essential amino acids (NEAA), and the gamma-glutamyl system during the interbout arousal period of hibernation. In the absence of nutrient intake or physical activity, our data illustrate the orchestration of metabolic pathways that sustain the provision of essential and non-essential amino acids and prevent ammonia toxicity during hibernation. Introduction: Hibernation is a state of extreme metabolic and physiologic plasticity that occurs in conjunction with physical inactivity and prolonged fasting.^[45]1 Specifically, hibernating Arctic Ground Squirrels (AGS) fast for up to eight months without any exogenous nitrogen intake, yet they are able to maintain physiological function and skeletal muscle.^[46]2–[47]4 Urea nitrogen salvage with nitrogen recycling into amino acids has been hypothesized for decades as a major mechanism to support anabolic processes while preserving protein content in hibernation.^[48]5–[49]8 Less focus has been devoted to free nitrogen. The nature of protein degradation, and how hibernating ground squirrels maintain nitrogen homeostasis and utilize free nitrogen to preserve amino acid stores, has remained unknown due to technical limitations of measuring in vivo fluxes in deep torpor. Hibernation in AGS is composed of torpor bouts (approximately fourteen to twenty days long) where the whole body metabolic rate is reduced to 1-2 percent of basal metabolic rate and core body temperature (T[b]) can dip below 0°C.^[50]3,[51]9 Torpor bouts are interrupted by naturally occurring interbout arousals, with a duration of 12-24 hours.^[52]10 As torpor progresses, amino acids with nitrogenous side chains accumulate in tissues of ground squirrels, signifying increased nitrogen load.^[53]11,[54]12 Multiple tissues, ranging from skeletal muscle, kidney, liver, small intestine, lung and brain, contribute to systemic nitrogen homeostasis, ammonia detoxification and sequestration of nitrogen into amino acids.^[55]13–[56]15 There is evidence anabolism occurs in liver and small intestine during hibernation, which is dependent on recycling of nitrogen into protein.^[57]7 Several other studies indicate a capacity of hibernators to recycle nitrogen, but most have focused on microbial urea nitrogen salvage.^[58]6–[59]8 Nitrogen recycling into non-essential and essential amino acids in hibernating bears (n=2) was measured in a preliminary study using a ^15N[2] urea tracer.^[60]16 This tracer formed ^15N ammonia in plasma, indicating that metabolic pathways facilitating recycling could use either urea or free ammonia.^[61]16 Previous studies have demonstrated that the levels of circulating ammonia remain constant throughout torpor but it is unknown how ammonia is utilized in nitrogen salvage.^[62]12,[63]17 While hypothesized to derive from breakdown of labile tissues,^[64]18 the source of circulating free amino acids and free nitrogen in hibernation is unknown. No study to our knowledge has measured in vivo production rates of protein degradation markers or amino acids in deep torpor. Further, the metabolic pathways and tissues that facilitate nitrogen buffering and recycling of amino acids are unknown. We tested the hypothesis that breakdown of skeletal muscle during torpor contributes to free nitrogen load and that this load is buffered through transamination reactions involving multiple tissues. Our approach included a pulse infusion of a tracer cocktail containing amino acids and metabolites in deep torpor and summer ([65]Fig. 1).^[66]19 In addition we performed metabolomic analysis on blood sampled from indwelling arterial cannula throughout entrance, torpor and arousal and we traced the fate of ^15N following infusion of ^15N ammonium acetate during arousal from torpor in AGS. Figure 1. Figure 1. [67]Open in a new tab Schematic for separate experimental procedures. Experiment #1 infused amino acid isotopes and sampled blood without inducing arousal, body temperature (T[b]) <4°C represents torpor, T[b] >34°C represents interbout arousal. Experiments were repeated in the same animals during summer euthermia in the post-absorptive state. Composition of isotope pulse infusion in deep torpor and summer euthermic animals during experiment #1 reported in nmol/g body weight, n=14 AGS, mean ±SEM. Experiment #2 sampled blood in undisturbed hibernators from entrance into hibernation through a full torpor bout. Experiment #3 induced arousal by handling animals, infused ^15N ammonium acetate and collected blood and tissues. Results: Skeletal muscle breakdown continues albeit at a slow rate in deep torpor Tau-methylhistidine (3-MH) is a skeletal muscle degradation marker and trans-4-hydroxy-L-proline (hPRO) is a collagen degradation marker.^[68]20,[69]21 We found decreases in both 3-MH and hPRO whole body production (WBP) in torpor compared to summer euthermia ([70]Fig. 2a, table insert, p<0.001, t-tests), but the relative reduction of 3-MH was significantly less than hPRO ([71]Fig. 2a, p<0.001, t-test). Specifically, hPRO WBP is reduced to 1% of euthermic rates, which reflects the overall 98% reduction in whole body metabolism during AGS torpor,^[72]9 while 3-MH is reduced to only 22% of euthermic levels. In normal physiology, collagen turnover is far higher than skeletal muscle turnover,^[73]21,[74]22 as is seen in our euthermic summer animals where hPRO WBP is 8.27 (nmol/g body weight/hr) compared to 3-MH WBP of 1.83 (nmol/g body weight/hr). Conversely, in deep torpor 3-MH WBP (0.41 nmol/g body weight/hr) is approximately three times higher than hPRO WBP (0.13 nmol/g body weight/hr). Figure 2. Figure 2. [75]Open in a new tab Skeletal muscle breakdown is on-going and whole body production (WBP) of all metabolites are depressed in deep torpor. a. In deep torpor, tau-methylhistidine (3-MH), a marker of skeletal muscle breakdown, is less repressed than hydroxy L-proline (hPRO), a marker of collagen breakdown (*p<0.001, WBP ratio hPRO vs 3-MH, two tailed t-test, boxplot whiskers represent min-max value, centre line represents the median and box boundaries represent 25^th and 75^th percentiles). Table insert: 3-MH and hPRO both significantly decline during torpor (n=9 AGS) compared to summer (n=5 AGS) (*p<0.001, two tailed t-test). b. Whole body production (WBP) of all metabolites are depressed in torpor compared to summer euthermia (n=5 AGS for all metabolites excluding tyrosine n=3 AGS and isoleucine n=4 AGS). WBP of glutamate in summer was too fast to generate a reliable decay curve and is not reported. (*p<0.05, **p<0.001 euthermic vs. deep torpor, paired two tailed t-test, FDR corrected; phenylalanine, leucine, valine, citrulline, glutamine, taurine p=0.000, tyrosine p=0.007, isoleucine p=0.023, ornithine p=0.001, glycine p=0.053, tryptophan p=0.011). c. WBP of amino acids is suppressed to varying degrees (torpor:summer WBP ratio, n=5 AGS). d. Linear regression analysis for tau-methylhisitidine WBP compared to metabolite WBP in deep torpor. Phenylalanine, glutamine and citrulline WBP are significantly correlated to tau-methylhistidine WBP in torpor (p<0.05, linear regression, n=9 AGS). Data shown are mean ±SEM The reduced suppression of 3-MH WBP compared to hPRO WBP indicates that myofibrillar (skeletal muscle) breakdown occurs throughout torpor, albeit at a low rate. Skeletal muscle breakdown may be a source of amino acids in torpor WBP of all amino acids indicate profound suppression of metabolic processes in deep torpor compared to summer ([76]Fig. 2b, p<0.05, paired t-tests, FDR corrected), however, the degree of suppression varied between amino acids ([77]Fig. 2c). Since our findings show that skeletal muscle degradation is ongoing during torpor, we next asked if amino acid production correlated with skeletal muscle breakdown. Regression analysis showed significant linear relationships between WBP rates of 3-MH in deep torpor and WBP of glutamine, phenylalanine and citrulline ([78]Figure 2d). The linear relationship between rates of WBP of 3-MH and WBP of phenylalanine and glutamine in torpor argue that these amino acids originate from skeletal muscle during torpor. Further, given that glutamine is a precursor to citrulline synthesis in the gut the correlation between citrulline and 3-MH WBP suggests that small intestine metabolism may be linked to glutamine release from skeletal muscle during deep torpor. Nitrogen carriers accumulate over deep torpor with a lack of kidney metabolic recycling Previous metabolic profiling in other hibernating species shows that plasma levels of certain amino acids deplete or accumulate as torpor progresses.^[79]11,[80]23–[81]25 We observed an increase in plasma glutamine, alanine and glutamate as torpor progressed over approximately two weeks ([82]Fig. 3, n=9, glutamate p<0.034, glutamine p<0.001, alanine p<0.001, repeated measures ANOVA). The accumulation of glutamine in liver and brain over torpor in other species support our finding of increased nitrogen load as torpor progresses.^[83]11,[84]12 While some previous plasma metabolomic studies do not mention a significant increase in glutamine or glutamate between entrance and late torpor, it is important to note that we evaluated these changes with high resolution, longitudinal sampling within animals.^[85]23,[86]25 Figure 3. Figure 3. [87]Open in a new tab Circulating nitrogen metabolite pools increase as torpor progresses while urea cycle intermediates do not increase. Data shown are relative abundance normalized to a sample collected at entrance into torpor at a T[b] of 11-12°C. The x-axis indicates duration of the torpor bout expressed as a percent of the total length of the bout in days. (p values compare specific time points in the torpor bout vs early torpor (10- 15%) posthoc LSD two-sided, n=9 AGS). Creatinine, creatine, trans-4-hydroxy-L-proline, allantoin, hypoxanthine, xanthine significantly increase over torpor (p values represent specific time points in the torpor bout vs early torpor (10- 15%) posthoc LSD two-sided, n=9 AGS). Urea cycle intermediates citrulline, ornithine and arginine do not increase over torpor (ns, repeated measures ANOVA). Measurements are from individual animals sampled over a single, undisturbed torpor bout. Data shown as mean ±SEM We also observed that urea cycle intermediates (arginine, ornithine, citrulline) do not increase over torpor ([88]Fig. 3, n=9), consistent with prior evidence that the urea cycle is suppressed during torpor.^[89]12,[90]17,[91]26 Ammonia, present in plasma throughout torpor, is naturally released during amino acid catabolism following protein degradation.^[92]12,[93]17,[94]27 Ammonia is typically incorporated into urea in the liver, but another major pathway for sequestering free nitrogen is via glutamine synthetase (GS) to form glutamine as the major nitrogen carrying amino acid, which occurs in multiple tissue types and skeletal muscle.^[95]13–[96]15,[97]27,[98]28 Thus, our results support a model whereby free nitrogen released during protein breakdown and catabolism accumulates in glutamine, glutamate, and alanine (products of transamination). Multiple waste byproducts of amino acid and purine catabolism, such as creatine, creatinine, and hPRO or hypoxanthine, allantoin and xanthine also increased, respectively ([99]Fig. 3, n=9, p<0.05 for each of these metabolites, repeated measures ANOVA). We interpret the accumulation of metabolites, such as creatinine, as evidence of amino acid catabolism throughout torpor as well as reduction in the kidney glomerular filtration rate (GFR), also supported by prior studies ([100]Fig. 3).^[101]26,[102]29 Accumulation of collagen breakdown products such as hPRO illustrates that production outpaces reabsorption or recycling during torpor. In total, these results suggest low-grade catabolism paired with reduced kidney function and accumulating nitrogen load with flux into non-essential amino acids instead of urea cycle intermediates. Alternatively, urea cycle intermediates may not accumulate due to loss via microbial urease enzyme activity.^[103]8 We next asked if AGS lessen the nitrogen load through metabolic recycling and anabolic metabolism once these processes resume during the interbout arousal period. Nitrogen fluxes into metabolites during arousal from torpor To understand how free nitrogen liberated from protein breakdown is incorporated into amino acids and metabolic pathways we next measured enrichment of ^15N derived from ^15N ammonium acetate infused during interbout arousal. We traced nitrogen incorporated into multiple metabolites in plasma and tissues ([104]Fig. 4). Overall ^15N ammonia incorporated into non-essential amino acids (NEAA) (glutamine, alanine, glutamate, arginine, asparagine, proline, serine, citrulline and ornithine) as well as essential amino acids (EAA) (leucine-isoleucine, valine, threonine, histidine) and the gamma-glutamyl intermediate 5-oxoproline ([105]Fig. 4). All tissues incorporated free ^15N into metabolic pathways ([106]Fig. 4). ^15N incorporation from ammonium acetate into EAA and NEAA supports de novo synthesis. Figure 4. Figure 4. [107]Open in a new tab Free ^15N ammonia is recycled into nonessential amino acids (blue), essential amino acids (red) and 5-oxoproline (turquoise) during arousal from torpor. Data show percent of ^15N incorporation during arousal from torpor following 72 mg/kg ^15N ammonium acetate pulse infusion, n=6, mean ±SEM. Calculations for percent ^15N incorporation: (^15N metabolite peak area)/(^15N metabolite peak area + ^14N metabolite peak area)*100 following natural abundance correction. Notably, we found significant incorporation of ^15N into glutamate in plasma (44.6±4.66% of total Glu), liver (13.5±3.7%), kidney (7.9±1.8%), skeletal muscle (0.7 ±0.2%), small intestine (1.4±0.3%) and lung (2.3±0.6%) (reported here as mean±SEM and denoted as M+1, [108]Fig. 5). The majority of our proposed pathways for ammonia ^15N incorporation during arousal from torpor stem first from synthesis of ^15N glutamate indicating a central role of this metabolite in hibernation nitrogen homeostasis. Previous in vitro studies suggest the reversibly phosphorylated glutamate dehydrogenase (GDH), which can incorporate free nitrogen into glutamate, catalyzes primarily the reverse reaction, glutamate oxidation, during hibernation.^[109]30,[110]31 GDH, however, has complex, strong allosteric regulators^[111]32 and our data shows a clear ability of ^15N to be incorporated into glutamate, indicating in vivo activity of this enzyme to generate glutamate in AGS hibernation. ^15N Enrichment of urea cycle intermediate ornithine was observed in the small intestine, but the percent of ^15N incorporation is less than other metabolites in this tissue ([112]Fig. 5). Figure 5. Figure 5. [113]Open in a new tab Proposed pathway for 15N ammonia recycling into amino acids depends on non-essential amino acid incorporation during arousal from torpor. Unlabeled metabolites (14N), denoted as M+0, measured in the pulse experiments are shown in stacked bar graphs (72 mg/kg 15N ammonium acetate pulse infusion n=6). Bar graphs show percent of 15N incorporation (M+1) into metabolites (72 mg/kg 15N ammonium acetate pulse infusion, n=6). Percent 15N incorporation (denoted as M+1) calculated as: ((15N metabolite peak area/(15N metabolite peak area + 14N metabolite peak area))*100 following natural abundance correction. Suggested location of 15N nitrogen is labeled by red circles, but analysis cannot determine specific location in metabolite. Carbon is labeled by blue circles, unlabeled nitrogen by green circles. Data shown mean ±SEM ^15N Incorporation into glutamate and arginosuccinate suggest that glutamate is a major contributor to ^15N-argininosuccinate and ^15N-arginine synthesis through ^15N-aspartate. ^15N-Aspartate was observed in liver (12.3±3.2%), plasma (29.3±6.3%), kidney (7.5±1.2%) and lung (2.1±0.6%) ([114]Fig. 5). ^15N-Argininosuccinate was observed in liver (4.9±1.4%) and ^15N-arginine was observed in plasma (0.4±0.1%) and kidney (0.2±0.1%) ([115]Fig. 5). Many of these synthesized NEAAs are direct protein precursors and exert protein anabolic effects, which represent metabolic pools available to stimulate protein synthesis during arousal to counteract skeletal muscle breakdown during torpor.^[116]27,[117]33 Recent metabolic modeling predicted increased NEAA synthesis in grizzly bears.^[118]34 These findings also support previous work suggesting that amino acid anabolism is upregulated^[119]4 and nitrogen is redeposited in liver and small intestine during AGS hibernation.^[120]7 ^15N Incorporated into leucine-isoleucine, threonine, and histidine suggests that ^15N-glutamate contributes to the synthesis of branched chain amino acids (BCAA) and essential amino acids (EAA). De novo EAA synthesis in hibernation has been proposed in the hibernating bear,^[121]16 a notion that has been challenged by subsequent studies suggesting that de novo EAA synthesis might not be required due to readily available EAA resources in labile proteins that can be catabolized over the season.^[122]35,[123]36 Incorporation of ^15N into leucine-isoleucine in kidney was 4.6±1.0% ([124]Fig. 5). ^15N-Glutamate is abundant in kidneys and glutamate donates nitrogen to branched-chain amino acid aminotransferases (BCAT), which catalyze nitrogen incorporation/removal between BCAAs and their α-keto acid analogues which are typically rich in kidney.^[125]37 A previous study in dogs revealed renal ability to reincorporate nitrogen into leucine’s keto-acid, supporting this capability in mammals.^[126]38 Exceptionally reduced renal glomerular filtration rate (GFR) in torpor is reversed with the return of kidney function during interbout arousal.^[127]26,[128]29,[129]39 Further, proteomic analysis of ground squirrel kidneys show that BCAA catabolism proteins decrease in hibernation^[130]40 and isotopic leucine tracer studies in black bears indicate very low levels of oxidized leucine in hibernation.^[131]41 Taken together with these studies, our data indicates that renal hibernation BCAA metabolism is shifted toward anabolism rather than catabolism. Further, ^15N incorporation into leucine correlated significantly, or near significantly, with T[b] in kidney, lung, plasma and skeletal muscle illustrating the temperature dependence of this process ([132]Extended Fig. 2). Incorporation of ^15N glutamine also varied with temperature in some tissues ([133]Extended Fig. 4–[134]5). ^15N Incorporation into glutamate only correlated with T[b] in muscle ([135]Extended Fig. 3). ^15N Incorporation into valine ([136]Extended Fig. 6) did not correlate with temperature in any tissue studied. This evidence suggests that rising T[b] facilitates nitrogen incorporation into some BCAA in multiple tissues. ^15N-Glutamate and ^15N-aspartate are also precursors to ^15N-threonine. ^15N Incorporation into threonine was observed in plasma (1.7±0.5%), liver (1.6±0.2%), kidney (0.5±0.2%), small intestine (1.6±0.2%) and lung (1.0±0.2%) ([137]Fig. 5). We observed the highest levels of ^15N incorporation into glutamate and aspartate, precursors to threonine, in the plasma and liver ([138]Fig. 5). We also saw that the highest levels of ^15N incorporation into threonine occur in the small intestine, plasma and liver ([139]Fig. 5). The small intestine serving as the site for threonine synthesis via ^15N-aspartate is supported by previous studies indicating that threonine synthesis in mammals relies on gut microbiota.^[140]42 Sections of the small intestine deliver amino acids to the hepatic portal vein, but our experiments cannot determine specific exchanges in the gut-liver axis. Important to note is that many microbial processes are still dependent on NEAA substrate availability, suggesting that metabolic activity of multiple tissues could play supporting roles in free nitrogen recycling for these metabolic pathways in hibernation.^[141]42 Glutamate and glutamine, along with ATP, are also known nitrogen donors during histidine synthesis.^[142]43 ^15N Incorporation into histidine was observed in plasma (1.4±0.3%), liver (0.2±0.6%) and lung (2.1±0.7%) ([143]Fig. 5). We propose ^15N glutamine or ^15N glutamate are likely sources of ^15N histidine via microbial synthesis; we did not observe ^15N incorporation into ATP. Given ATP’s labile nature, however, lack of detection cannot rule out potential incorporation of ^15N into ATP and lack of detection may be attributable to technical limitations. Interestingly, ^15N histidine mainly appears in lung and plasma and not small intestine ([144]Fig. 5). Most microbial EAA studies focus on the ileum.^[145]42,[146]44 One explanation for our lack of ^15N histidine labeling could be because we did not distinguish between the jejunum and the ileum raising the possibility that our tissue sampling protocol missed a site of small intestine histidine synthesis. The unexpected lung ^15N incorporation into histidine parallels the overall broad ^15N incorporation into multiple metabolites in lung ([147]Fig. 4). This either suggests lung is an important tissue for nitrogen homeostasis in hibernation or that amino acids carrying recycled nitrogen are important for protein synthesis in lung tissue. While known to play a role in amino acid metabolism, lung metabolism in hibernation is a sparsely investigated field to the best of our knowledge. The high amount of ^15N incorporation into histidine and other amino acids in lung suggest that further exploration is worthwhile. ^15N Incorporated into glutamine and glutamate also showed evidence of transfer into 5-oxoproline through the gamma-glutamyl cycle. We observed incorporation of ^15N into glutamine in plasma (14.1±3.5%), liver (15.1±3.9%), kidney (10.9±4.0%), skeletal muscle (1.9±0.4%), small intestine (3.0±1.6%), and lung (8.4±2.9%) ([148]Fig. 5) and incorporation of ^15N into 5-oxoproline in plasma (6.4 ±1.4%), liver (9.0±2.2%), kidney (5.3±1.5%), skeletal muscle (0.4±0.1%), small intestine (1.3±0.4%) and lung (2.7±0.7%) ([149]Fig. 5). Our data is supported by previous findings of significantly increased gamma-glutamylated amino acids during arousal from torpor, showing an active γ-GT (gamma-glutamyltransferase) system, and strong upregulation of antioxidant systems in hibernation.^[150]23,[151]45 In support of this pathway, data additionally shows that the absolute concentrations of naturally occurring ^14N-5-oxoproline significantly correlate with ^14N-glutamine and ^14N-glutamate levels ([152]Supp. Tables 1a and [153]1b). We were surprised by the high abundance of naturally occurring ^14N 5-oxoproline in plasma (393.34±33.42 µM), liver (6.93±0.76 µM), kidney (3.72±0.39 µM), skeletal muscle (11.15±3.78 µM), small intestine (1.35±0.16 µM) and lung (2.33±0.21 µM) ([154]Supp. Table 1a). Although 5-oxoproline is a common metabolite across different life forms, in humans increases in 5-oxoproline are documented in pathologies.^[155]46–[156]48 The gamma-glutamyl enzymes are present in most tissues and essential for glutathione function. Previous studies indicate that glutathione can act as a reserve for glutamate with some hypothesizing 5-oxoproline could do the same and that γ-GT is an alternative amino acid transport system.^[157]49–[158]51 Given the elevated absolute concentration of 5-oxoproline in all tissues and its correlation to glutamine and glutamate, 5-oxoproline may represent a substantial pool for nitrogen sequestration in hibernation which could further support antioxidant function and amino acid pools. Ketoacids increase during natural arousal from torpor During torpor, WBP of branched chain ketoacids and metabolites (ketoisocaproic acid (KIC), ketoisovalerate acid (KIV), keto-beta-methylvalerate (KMV) and beta-hydroxy-beta-methylbutyrate (HMB)) are suppressed ([159]Fig. 6a, p<0.002, WBP torpor vs. WBP summer euthermic, paired two-tailed t-test, n=5). Although it was beyond the scope of this study to capture WBP of these ketoacids during arousal, we next asked if the circulating levels of KIC and KIV and their ketone substrate (acetoacetate) increased in natural arousal from torpor to support nitrogen recycling. As arousal progresses, KIC, KIV and acetoacetate increase as T[b] increases ([160]Fig. 6b, p<0.001, repeated ANOVA, n=6). Interestingly, in AGS cecal content, branched chain fatty acids are higher in arousal compared to summer and torpor and acetate is higher during hibernation, suggesting microbial metabolism may support production of BCAA ketoacids from branched chain fatty acids.^[161]52 At 16°C leucine/isoleucine decreases, potentially indicating an uptake into tissues ([162]Fig. 6b, p<0.001, repeated ANOVA, n=6). Pathway enrichment analysis identified amino acid metabolism and BCAA biosynthesis as metabolically significant over the course of arousal ([163]Fig. 7, p<0.001 ANOVA, FDR-adjusted p value). The increase of BCAA ketoacids and their ketone parent (acetoacetate) in natural arousal from torpor supports nitrogen recycling by increasing availability of substrate for BCAA synthesis. Figure 6: Branched chained ketoacid WBP is suppressed in torpor, but ketoacids and ketones increase during arousal from torpor. Figure 6: [164]Open in a new tab a.) Ketoisocaproic acid (KIC), ketoisovalerate (KIV), keto-beta-methylvalerate (KMV) and betahydroxy-beta-methylbutyrate (HMB) whole body production (WBP) are suppressed in torpor compared to summer euthermia (**p<0.002, WBP torpor vs. WBP summer euthermic, paired two tailed t-test, n=5 AGS). b.) KIC, KIV, 3- hydroxy-3-methyl-2-oxobutanoic acid and acetoacetate increase over arousal (p values compare specific time points vs early arousal (2.5°C), posthoc LSD two-sided, n=6 AGS). Valine does not significantly increase over arousal (ns, repeated measures ANOVA) and leucine-isoleucine decreases (*p<0.05 vs early arousal (2.5°C), posthoc LSD, n=6 AGS). Measurements are from individual animals sampled over a natural arousal. Measurements are normalized to the relative abundance of sample collected during entrance into hibernation at T[b] 11-12°C. Data shown as mean ±SEM Figure 7: Figure 7: [165]Open in a new tab Pathway enrichment analysis shows evidence for amino acid metabolism and branched chain amino acid biosynthesis in arousal from torpor, and points to the prevalence of transamination reactions involving glutamate, alanine and aspartate during torpor. a.) Arousal pathway enrichment analysis compared early arousal (2.5°C) and full arousal (35°C) using MetaboAnalyst 4.0 (n=6) b.) Torpor pathway enrichment analysis compared early torpor (10-15% of torpor bout completed) to late torpor (88-100% of torpor bout completed) using MetaboAnalyst 4.0 (n=9) Discussion: The source and fate of free nitrogen during hibernation has until now been poorly understood. During interbout arousal, free nitrogen is recycled and buffered by transamination reactions into glutamate and glutamine and then incorporated into essential amino acids (EAA), non-essential amino acids (NEAAs) and the gamma-glutamyl intermediate 5-oxoproline. Nitrogen recycling has been proposed for decades as a potential mechanism of nitrogen homeostasis in hibernation. Our work, for the first time, demonstrates metabolic pathways facilitating nitrogen recycling in vivo and shows that free nitrogen is incorporated into metabolites. These metabolites may support protein synthesis and antioxidant capacity during arousal from torpor. Previous studies hypothesized that free EAAs, which support skeletal muscle and whole body protein synthesis, originate from collagen breakdown in labile protein sources during hibernation.^[166]4,[167]18,[168]35,[169]36 We show that, although suppressed, skeletal muscle breakdown is greater than collagen breakdown in torpor suggesting that skeletal muscle breakdown contributes to circulating amino acids and increased nitrogen load in torpor to a greater degree than collagen breakdown. Nitrogen recycling has been proposed for decades as a potential mechanism of metabolic homeostasis and physiological support in hibernation.^[170]6,[171]7,[172]16 Our work demonstrates unequivocally nitrogen recycling and de novo EAA synthesis and highlights pathways that support this metabolic mechanism. Our results suggest that inter-organ cooperation in metabolism supports synthesis of metabolites. These metabolites are essential components of tissue protein synthesis and have anabolic signaling properties and other physiological roles.^[173]27 Our findings highlighting recycling of nitrogen into NEAAs to support EAA synthesis expand on previous models predicting increased NEAA synthesis in hibernating bears.^[174]34 15N Incorporation into 5-oxoproline also represents a novel pathway for recycling and buffering free nitrogen, while supporting previous studies suggesting 5-oxoproline and glutathione could act as pools for glutamate.^[175]49,[176]50 In the greater context, antioxidants are known to be upregulated in hibernation and protect skeletal muscle from injury.^[177]45,[178]53 Sequestration of nitrogen into the gamma-glutamyl cycle would support glutathione synthesis and antioxidant defense. Hibernators preserve their muscle mass and functionality despite extreme inactivity and prolonged fasting.^[179]2,[180]4,[181]18,[182]54–[183]57 While previous research has debated the extent of muscle preservation and the mechanisms responsible, we have documented the retention of lean tissue mass in AGS using dual energy x-ray absorptiometry during prolonged inactivity or forced disuse (neurectomy).^[184]2 Similarly, preservation of lean tissue mass determined by naturally occurring tissue nitrogen isotope ratios (δ15N) was noted in AGS housed throughout the hibernation season at mild ambient temperature of 2°C.^[185]7 By contrast, when housed at a colder ambient temperature (−10°C) such as experienced in the wild, there was a reduction in lean tissue mass in proportion to the total number of days of the hibernation season.^[186]7,[187]9 Thus, milder ambient temperatures reveal a resistance to disuse atrophy not seen in non-hibernating species.^[188]2,[189]55,[190]56,[191]58 Our analysis of liver and small intestine support Lee’s interpretation of protein synthesis in these tissues during the hibernation season and suggest protein synthesis may occur in tissues such as lung, kidney and skeletal muscle to a lesser extent. Our findings also support previous gene expression studies in AGS muscle that indicate protein anabolism persists throughout the hibernation season while proteasome, autophagy and atrogene gene expression are not changed or are downregulated.^[192]4,[193]59 Metabolic pathways facilitating ^15N recycling could directly support muscle protein synthesis during interbout arousal.^[194]1,[195]4,[196]18,[197]55,[198]58,[199]60 While our data clearly shows enrichment of ^15N in these key energy metabolites, there are some important caveats to consider. We observed lower ^15N incorporation into EAA in skeletal muscle than in kidney, lung or liver. One reason for this could be regional temperature differences during rewarming.^[200]61,[201]62 Tissues were collected 2 hours after IV bolus infusion during arousal; the head and thoracic region of the body warms more quickly than distal regions partially due to proximity of brown adipose tissue deposits and distal vasoconstriction.^[202]61 Further, in AGS arousal, blood flow to the extremities significantly lags behind brain and organ blood flow, meaning that metabolic exchanges with circulation occur more quickly for the brain and visceral tissues than for the quadriceps.^[203]62 If experiments had taken place during full arousal (T[b]>34°C) for a longer duration, we might have seen higher incorporation of labeled ^15N in all tissues and peripheral skeletal muscle that lags in rewarming. Indeed, we observed temperature dependent incorporation of ^15N into amino acids in several tissues. We designed the ^15N ammonium acetate infusion to standardize the time between infusion and sample collection. As a result, we collected tissues across a wide range of T[b] due to variation in rates of rewarming. This allowed us to define linear relationships between T[b] and degree of ^15N incorporation in some tissues. The novel insight from this approach is that many of these mechanisms incorporate ^15N into metabolites at temperatures as low as 5°C (albeit, at reduced rates). Ketoacids, however, only increase as temperatures reach 16°C. Therefore, ^15N incorporation into leucine/isoleucine and large-scale synthesis of EAA may be temperature dependent ([204]Extended Fig. 2). This is consistent with previous findings that translation has a temperature threshold of 18°C in 13-lined ground squirrels.^[205]63 The scope of this study did not include measures of 3-MH or hPRO WBP rates during arousal. Arousal from torpor has a period of intense shivering and rapid physiological and metabolic changes. While we saw high rates of 3-MH production in deep torpor, it remains unknown if this trend continues or reverses during the arousal period. Additional metabolic flux analysis was not possible during interbout arousal because isotope from one infusion could have confounded an additional pulse study in the same animal during the hibernation season. While most of the total mammalian 3-MH pools exist in skeletal muscle, 3-MH is also present in actin in all cell types. Prior studies have pointed out that increased rates of 3-MH production could be partially caused by high rates of turnover in these actin pools.^[206]22 We interpret 3-MH as a marker of skeletal protein degradation rather than actin turnover due to the direct correlation of 3-MH to catabolism in leg muscle.^[207]20 In summary, our results support a model of skeletal muscle breakdown and nitrogen accumulation in torpor with ammonia salvage pathways funneled away from the urea cycle and into, EAA, NEAA, and the gamma-glutamyl pathway during arousal. Many of these amino acids are direct precursors and signaling metabolites necessary for protein synthesis. This work has translational significance as a model of metabolic resilience during extreme inactivity and nutrient deficiency. AGS utilize specific metabolic mechanisms to sustain nitrogen stores and deliver critical amino acids which are essential for protein synthesis and organ function. In the absence of mechanical loading, the maintenance oriented mechanisms responsible for the preservation of musculoskeletal system and physical function in hibernation are highly relevant to our understanding of clinical perturbations such as sarcopenia, cachexia and even spaceflight.^[208]64 Now recognized as therapeutic targets for accelerated muscle atrophy in humans, these protective mechanisms may also be positively influenced by the preservation of molecular mechanisms such as miRNAs that dampen the influence of myostatin and ubiquitin ligase in hibernating species.^[209]65,[210]66 Future work should incorporate these active metabolic mechanisms with known molecular mechanisms that augment protein maintenance in hibernation to find novel therapeutic targets. Methods: Trapping: Arctic Ground Squirrels (AGS, Urocitellus parryii) were trapped in early July in the northern foothills of the Brooks Range, Alaska, 40 miles south of Toolik Field Station (68°38 N, 149°38 W; elevation 809 m) and were transported to Fairbanks, Alaska under permit by the Alaska Department of Fish and Game (permit numbers 18-111 and 17-156). Ethic Overview: All procedures were performed in accordance with University of Alaska Fairbanks Institutional Animal Care and use Committee (IACUC, protocols 491809 and 996835). Husbandry: AGS were housed in ambient temperature (T[a]) between 16-18°C at 16:8-hour light/dark cycle until August 15^th, when they were moved to cold chambers with T[a] of 2°C at a 4:20 hour light/dark cycle. Animals were fed Mazuri Rodent chow (#5663 Mazuri, PMI Nutrition International, Richmond, IN, USA) and provided water during the euthermic period. Animals were housed individually in 12”×19”×12” stainless steel wire mesh hanging cages with cotton nests over ammonia absorbing corn cob litter. Once animals exhibited robust hibernation behavior, food was withdrawn; animals were placed in polycarbonate cages (8.5”×17”×8.5”) with shavings, cotton bedding and gel hydration packets. Polycarbonate cages were placed on receivers linked to a data collection system for core body temperature (Data Sciences International, St. Paul, MN, USA). Surgery: AGS were instrumented with chronic femoral arterial and venous cannulas (3 Fr cannula, Instech Laboratories Inc, Plymouth Meeting, PA, USA) and either TA-F40 or CTA-F40 core body temperature loggers (Data Sciences International, St. Paul, MN, USA) in July and August. Cannulae patency was maintained with a heparin/glycerol locking solution (1:1). Chemicals: A mixture of amino acid isotopes (Cambridge Isotope Laboratories, Woburn, MA USA) was compounded professionally in a compounding pharmacy (Temple, TX). For specific mixture composition see [211]Fig. 1. ^15N Ammonium acetate (98%, Sigma-Aldrich, St. Louis, MO) was diluted in saline and sterilized by 0.2 µm filtration. Definition of hibernation state: Torpor is defined as the time when core body temperature is below 4°C. The first day of torpor is defined by a T[b] of 4°C reached after a period of stable euthermic T[b] (34-36°C). Entrance into torpor is defined during the period of decreasing T[b] when T[b] is within 10-11°C. We recorded body weight during the euthermic period preceding experiments. Deep Torpor and Summer: Design for Whole Body Production Measurement with Pulse Isotope Infusions: We administered mixed pulse isotope infusions ([212]Fig. 1) in February during deep torpor on the third day of a torpor bout ([213]Fig. 1, female n=5, male n=4). Infusions did not affect torpor (torpid T[b] did not vary more than 0.12°C during experimental procedures). We flushed AGS cannula with sterile saline on the first and second day of torpor to habituate animals to cannula handling. If flushing induced arousal, animals were excluded from experimental procedures. Experiments began in the morning on the third day of torpor. Heparinized saline locking solution was removed and blood was sampled from the cannula prior to pulse infusion of a mixed amino acid isotope cocktail (0 hour). Blood was then sampled at precisely 10 min, 20 min, 30 min, 1 hour, 2 hour, 3 hour, 4 hour, 12 hour and 24 hour post infusion. Control experiments were repeated in summer (June) with the same animals under post-absorptive conditions after an overnight fast (female n=2, male n=3, cannuale patency was lost on a few animals over the winter season resulting in the decreased summer sample size). Blood sampled during the summer season was sampled immediately prior to infusion (0 hour), and at 5 min, 10 min, 15 min, 20 min, 30 min, 1 hour, 1.5 hour and 2 hour. Blood was sampled via the arterial line and collected into a 1cc syringe via a pin port and immediately transferred to a non-heparinized micro-hematocrit tube (Kimble) and sealed with putty. Blood was centrifuged for 2 min at 4°C and plasma was immediately placed on dry ice. Samples were transferred to a −80°C freezer and stored until shipped to Texas A&M for analysis. Core body temperature (Tb) was recorded at 10-minute intervals throughout the hibernation season. Deep Torpor and Summer Quantification of rates of Whole Body Metabolite Production: Deep torpor and summer metabolites and enrichments of tracers were measured batch-wise by stable tracer-dilution LC-MS/MS as previously described.^[214]19 While amino acid and metabolite production in vivo is often measured by primed-constant infusions of stable isotope tracers, torpor and arousal from torpor are non-steady-state physiological conditions with unknown metabolite production rates and metabolite pool sizes.^[215]19,[216]67 Pulse isotope tracer methodology is the preferred technique to measure non-steady state conditions by measuring decay curves with noncompartmental modeling.^[217]19,[218]67,[219]68 In brief, deproteinized plasma supernatants were diluted in reagents required for a 9-fluorenylmethoxycarbonyl (Fmoc) reaction. Reactions were stopped after 15 minutes. We injected 0.2–1.7 μL into the 24 μl/min pumped flow (microLC 200, Eksigent, part of AB Sciex) of 73:23:4 water-acetonitrile-isopropanol containing 11 mM ammonium acetate, onto a 100 × 0.5 mm column packed by Eksigent with 2.7 μm 90 Å HALO C18 fused-core silica beads (Advanced Materials Technology). Fmoc-amino acids eluted over 3.2 min; bis-Fmoc-amino acids (ornithine, and tyrosine) eluted by 4.5 min. Ions were conveyed by heated electrospray (−3.8 kV IS; 150°C TEM; 14 GAS1; 25 GAS2) to a 5500 QTRAP (Sciex) for multiple reaction monitoring (MRM) of the loss of the Fmoc moiety or moieties, regenerating the amino acid anion. Collision energy (CE) was de-optimized as needed for analytes that would otherwise risk detector saturation. Peaks were automatically identified and integrated by the SignalFinder1 algorithm in MultiQuant v. 3.0 (Sciex) and integrations were inspected. Manual reintegration (redrawing of baselines and peak boundaries) was disallowed, but manual reidentification was allowed as needed. Peak areas were exported to Excel for calculation of area ratios for 1) concentration measurements using a tracer for the internal standard and 2) enrichment measurements of an exogenous tracer (or tracer metabolite) to the corresponding amino acid. Area ratios of calibrators were regressed against calibrator known tracer-to-tracee ratios (TTR) by ordinary least squares (Prism 8, GraphPad). By reverse prediction, each plasma sample area ratio was converted to a TTR. Previous experiments have documented that the pulse isotope infusion does not affect the tracee concentration.^[220]19 Calculations of Whole Body Production: The decay of TTR over time was normalized with body weight and amount of administered tracer (Prism 8, Graphpad). Decay of TTR amino acid isotopes were fitted to the equation y=a*exp(−k1*x)+b*exp(−k2*x), where k1 and k2 are decay rate constants of tracer and tracee, and area under the curve (AUC) was calculated from the integral of the two exponential curves.^[221]68 Rate of appearance (Ra) for each amino acid is calculated by: dose of metabolite in the pulse infusion/decay of metabolite tracer (AUC). Whole body rate of appearance (Ra) of amino acids is a proxy for whole body production (WBP). Metabolomic Deep Torpor Profiling Experimental Procedure: Animals were habituated to cannula flushing with sterile saline. Plasma was sampled from naïve hibernating AGS without disturbing hibernation over a full torpor bout ([222]Fig. 1, female n=5, male n=4). Sampling occurred during entrance to torpor (T[b] 11-12°C) and sequentially through a full torpor bout. Six animals were also sampled through a natural arousal at four time points during re-warming: 2.5°C, 4°C, 16°C and 35°C. A typical torpor bout is a 14 to 20-day period while a typical arousal is a 12-24 hour period, though timing varies depending on ambient temperature and season.^[223]10 Relative abundance of metabolites was normalized to entrance levels (T[b] 11-12°C). Duration of torpor was expressed as a percentage of the total length of the torpor bout in days. Longitudinal Global Metabolomic Quantification: Plasma samples (20 µl) were extracted using ice cold methanol:acetonitrile:water (5:3:2) at 1:25 dilution as described.^[224]69 Extracts were analyzed on a Thermo Vanquish ultra-high performance liquid chromatograph (UHPLC) coupled online to a Thermo Q Exactive mass spectrometer using a 5 min C18 gradient method in positive and negative modes as described previously.^[225]24,[226]70 Peak picking and metabolite assignment were performed using Maven (Princeton University) against the KEGG database ([227]https://www.genome.jp/kegg),^[228]71,[229]72 confirmed against chemical formula determination from isotopic patterns and accurate mass, and validated against experimental retention times for >650 standard compounds (Sigma Aldrich; MLSMS, IROATech, Bolton, MA, USA).^[230]73 Compartmental Nitrogen (^15N) Flux in Arousal from Torpor Experimental Procedure: ^15N-Ammonium acetate experiments occurred on the third day of a torpor bout in February in a separate animal group from previous tracer studies. Arousal was induced by handling between 0730 and 0830 in the morning. Importantly, while re-warming rates vary slightly between handling-induced vs. spontaneous arousal, inducing arousal was not expected to influence outcome of the study. As Williams et al. (2011) has noted, differences in metabolic energetic demand on AGS (in a different circumstance, hibernating at two different temperatures), did not alter transcript levels in major metabolic tissues such as BAT and liver.^[231]74 A ^15N ammonium acetate venous bolus of either 72 mg/kg (female n=3, male n=3) or 360 mg/kg (female n=3, male n=2) at a volume of 1 ml/kg was infused 1.5 hours after inducing arousal ([232]Fig. 1). Blood was sampled two hours post infusion under a surgical plane of anesthesia (induced at 5% isoflurane with 100% medical grade oxygen). Animals were euthanized (by decapitation) directly afterwards and tissues were rapidly sampled and frozen on dry ice. Because the time between bolus and sample collection was standardized to two hours post infusion, rectal temperature of AGS taken at time of tissue collection varied from 4.9-24.4°C. Tissue samples included liver (left lobe), kidney (cortex), skeletal muscle (right quadriceps femoris), lung (right), plasma, and small intestine (cleaned from mesentery). Tissue samples were extracted in ice-cold lysis/extraction buffer (described above) at a concentration of 15 mg/mL in the presence of glass beads. Samples were agitated via bead beater for 3 min at 4°C then vortexed for 30 minutes and centrifuged at 18,200 g for 10 minutes at 4°C. Plasma samples (20 µl) were extracted as described above. ^15N Tracer Metabolite Quantification: Plasma and tissue extracts were analyzed by UHPLC-MS as described above (see longitudinal metabolomics sampling). Animals in this experiment had no previous tracer infusions. Major metabolites for nitrogen metabolism were inspected for ^15N incorporation including purine, amino acid, creatine, spermidine, urea cycle and 5-oxoproline metabolism. Isotopologue distributions were corrected for natural abundance and results were plotted in GraphPad Prism 8.0. Statistical Analysis: Results were expressed as mean +/− standard error of the mean (S.E.M). Significance of difference was assessed by repeated measures ANOVA, post hoc LSD and two-tailed student’s t-test as specified using SPSS (version 25, IBM SPSS Statistics, Armornk, NY). Measurements that violated sphericity were reported with greenhouse-geisser (G-G) corrections and t-tests that violated Levene’s test for equality were performed without assuming equal variances. Pathway enrichment analysis was performed using MetaboAnalyst 4.0. GraphPad Prism (version 8, GraphPad Software Inc., La Jolla, CA) was used for graphic presentation. Treatment groups were matched for sex and hibernation phase timing. Treatment was not blinded and no data was excluded. Sample size was determined by availability of wild-caught animals, but guided by a priori power analysis based on published literature.^[233]11 A priori power analysis (G*Power software) was performed to estimate sample size needed to yield 80% power with 0.05 α error probability for a two tailed t-test. A priori power analysis indicated that a sample size of eight animals with an expected difference in means of 0.20 and standard deviations of 0.15 and 0.11 (based on metabolite data taken from thirteen lined ground squirrel late torpor versus entrance) would have power of 0.81. In our studies, a sample size of nine was chosen to account for attrition. Extended Data Extended Figure 1. Extended Figure 1. [234]Open in a new tab Tracer to tracee ratio (TTR) decay curves show slow decay in torpor (red, n=9) compared to summer euthermic AGS (blue, n=5). Decay of TTR amino acid isotopes were fitted to the equation y=a*exp(−k1*×)+b*exp(−k2*×), and area under the curve (AUC) was calculated from the integral of the two exponential curves. Rate of appearance (Ra) for each amino acid is calculated by: dose of metabolite in the pulse infusion/AUC. Whole body rate of appearance (Ra) of amino acids is a proxy for whole body production (WBP). Extended Figure 2: Extended Figure 2: [235]Open in a new tab Linear regression analysis indicates ^15N incorporation into leucine/isoleucine is correlated to core body temperature during arousal from torpor (72 mg/kg ^15N ammonium acetate pulse infusion n=6, 360 mg/kg ^15N ammonium acetate pulse infusion n=5). In tissues where nitrogen incorporation was not observed in more than one animal, regression analysis was not preformed. Extended Figure 3: Extended Figure 3: [236]Open in a new tab Linear regression analysis indicates ^15N incorporation glutamate in skeletal muscle is correlated to core body temperature during arousal from torpor (72 mg/kg ^15N ammonium acetate pulse infusion n=6, 360 mg/kg ^15N ammonium acetate pulse infusion n=5). In tissues where nitrogen incorporation was not observed in more than one animal, regression analysis was not preformed. Extended Figure 4: Extended Figure 4: [237]Open in a new tab Linear regression analysis indicates ^15N incorporation into glutamine is correlated to core body temperature during arousal from torpor in kidney, plasma and liver (72 mg/kg ^15N ammonium acetate pulse infusion n=6, 360 mg/kg ^15N ammonium acetate pulse infusion n=5). In tissues where nitrogen incorporation was not observed in more than one animal, regression analysis was not preformed. Extended Figure 5: Extended Figure 5: [238]Open in a new tab Linear regression analysis indicates ^15N incorporation into Glutamine M+2 is correlated to core body temperature in kidney and liver during arousal from torpor (72 mg/kg ^15N ammonium acetate pulse infusion n=6, 360 mg/kg ^15N ammonium acetate pulse infusion n=5). In tissues where nitrogen incorporation was not observed in more than one animal, regression analysis was not preformed. Extended Figure 6: Extended Figure 6: [239]Open in a new tab Linear regression in analysis indicates ^15N incorporation into valine cannot be verified to rely on core body temperature during arousal from torpor. In tissues where nitrogen incorporation was not observed in more than one animal, regression analysis was not preformed Extended Figure 7. Extended Figure 7. [240]Open in a new tab Free ^15N ammonia is recycled into nonessential amino acids (blue), essential amino acids (red) and 5-oxoproline (turquoise) during arousal from torpor (360 mg/kg ^15N ammonium acetate pulse infusion, n=5, mean ±SEM). Percent ^15N incorporation calculated as: (^15N metabolite peak area/(^15N metabolite peak area + ^14N metabolite peak area))*100 following natural abundance correction. Supplementary Material 1638169_Supp_Tab1 [241]NIHMS1638169-supplement-1638169_Supp_Tab1.pdf^ (56.1KB, pdf) Acknowledgements