Proteomic and Transcriptomic Changes in Hibernating Grizzly Bears Reveal Metabolic and Signaling Pathways that Protect against Muscle Atrophy

Molecular regulation of reversible cardiac remodeling: lessons from species with extreme physiological adaptations

Some vertebrates evolved to have a remarkable capacity for anatomical and physiological plasticity in response to environmental challenges. One example of such plasticity can be found in the ambush-hunting snakes of the genus Python, which exhibit reversible cardiac growth with feeding. The predation strategy employed by pythons is associated with months-long fasts that are arrested by ingestion of large prey. Consequently, digestion compels a dramatic increase in metabolic rate and hypertrophy of multiple organs, including the heart. In this Review, we summarize the post-prandial cardiac adaptations in pythons at the whole-heart, cellular and molecular scales. We highlight circulating factors and cellular signaling pathways that are altered during digestion to affect cardiac form and function and propose possible mechanisms that may drive the post-digestion regression of cardiac mass. Adaptive physiological cardiac hypertrophy has also been observed in other vertebrates, including in fish acclimated to cold water, birds flying at high altitudes and exercising mammals. To reveal potential evolutionarily conserved features, we summarize the molecular signatures of reversible cardiac remodeling identified in these species and compare them with those of pythons. Finally, we offer a perspective on the potential of biomimetics targeting the natural biology of pythons as therapeutics for human heart disease.

O-GlcNAcylation of SERCA protects skeletal muscle in hibernating Spermophilus dauricus from disuse atrophy

Long-term inactivity of skeletal muscle results in muscular disuse atrophy; however, hibernating animals do not experience muscular disuse atrophy during the hibernation period. The molecular mechanism underlining the anti-atrophy effect in these animals is unclear. O-linked N acetyl-β-D-glucosaminylation (O-GlcNAcylation) and its effect on cell signaling pathways are important mechanisms underlying muscular disuse atrophy; thus, in this study, we investigated O-GlcNAcylation changes during hibernation in Spermophilus dauricus to explore the role of O-GlcNAcylation in the muscle disuse atrophy resistance of hibernating animals. The results showed that during hibernation, the muscle fiber cross-sectional area and ratio of muscle fiber did not change, and the morphological structure of the muscle remained intact, with normal contractile function. The level of O-GlcNAcylation decreased during hibernation, but quickly returned to normal in the periodic arousal stage. The O-GlcNAcylation level of sarcoplasmic/endoplasmic reticulum calcium ATPase 1 (SERCA1) decreased, whereas its activity increased. The decrease in O-GlcNAcylation of SERCA could result in the decreased binding of phospholamban to SERCA1, thus decreasing its inhibition to SERCA1 activity. This in turn can inhibit muscle cell calcium overload, maintain muscle cell calcium homeostasis, and stabilize the calpain proteolytic pathway, ultimately inhibiting skeletal muscle atrophy. Our results demonstrate that periodic arousal along with returning O-GlcNAcylation level to normal are important mechanisms in preventing disuse atrophy of skeletal muscle during hibernation.

Remodeling of skeletal muscle myosin metabolic states in hibernating mammals

Hibernation is a period of metabolic suppression utilized by many small and large mammal species to survive during winter periods. As the underlying cellular and molecular mechanisms remain incompletely understood, our study aimed to determine whether skeletal muscle myosin and its metabolic efficiency undergo alterations during hibernation to optimize energy utilization. We isolated muscle fibers from small hibernators, Ictidomys tridecemlineatus and Eliomys quercinus and larger hibernators, Ursus arctos and Ursus americanus. We then conducted loaded Mant-ATP chase experiments alongside X-ray diffraction to measure resting myosin dynamics and its ATP demand. In parallel, we performed multiple proteomics analyses. Our results showed a preservation of myosin structure in U. arctos and U. americanus during hibernation, whilst in I. tridecemlineatus and E. quercinus, changes in myosin metabolic states during torpor unexpectedly led to higher levels in energy expenditure of type II, fast-twitch muscle fibers at ambient lab temperatures (20 °C). Upon repeating loaded Mant-ATP chase experiments at 8 °C (near the body temperature of torpid animals), we found that myosin ATP consumption in type II muscle fibers was reduced by 77-107% during torpor compared to active periods. Additionally, we observed Myh2 hyper-phosphorylation during torpor in I. tridecemilineatus, which was predicted to stabilize the myosin molecule. This may act as a potential molecular mechanism mitigating myosin-associated increases in skeletal muscle energy expenditure during periods of torpor in response to cold exposure. Altogether, we demonstrate that resting myosin is altered in hibernating mammals, contributing to significant changes to the ATP consumption of skeletal muscle. Additionally, we observe that it is further altered in response to cold exposure and highlight myosin as a potentially contributor to skeletal muscle non-shivering thermogenesis.

Advancing cancer cachexia diagnosis with -omics technology and exercise as molecular medicine

Muscle atrophy exacerbates disease outcomes and increases mortality, whereas the preservation of skeletal muscle mass and function play pivotal roles in ensuring long-term health and overall quality-of-life. Muscle atrophy represents a significant clinical challenge, involving the continued loss of muscle mass and strength, which frequently accompany the development of numerous types of cancer. Cancer cachexia is a highly prevalent multifactorial syndrome, and although cachexia is one of the main causes of cancer-related deaths, there are still no approved management strategies for the disease. The etiology of this condition is based on the upregulation of systemic inflammation factors and catabolic stimuli, resulting in the inhibition of protein synthesis and enhancement of protein degradation. Numerous necessary cellular processes are disrupted by cachectic pathology, which mediate intracellular signalling pathways resulting in the net loss of muscle and organelles. However, the exact underpinning molecular mechanisms of how these changes are orchestrated are incompletely understood. Much work is still required, but structured exercise has the capacity to counteract numerous detrimental effects linked to cancer cachexia. Primarily through the stimulation of muscle protein synthesis, enhancement of mitochondrial function, and the release of myokines. As a result, muscle mass and strength increase, leading to improved mobility, and quality-of-life. This review summarises existing knowledge of the complex molecular networks that regulate cancer cachexia and exercise, highlighting the molecular interplay between the two for potential therapeutic intervention. Finally, the utility of mass spectrometry-based proteomics is considered as a way of establishing early diagnostic biomarkers of cachectic patients.

Remodelling of Skeletal Muscle Myosin Metabolic States in Hibernating Mammals

Hibernation is a period of metabolic suppression utilized by many small and large mammal species to survive during winter periods. As the underlying cellular and molecular mechanisms remain incompletely understood, our study aimed to determine whether skeletal muscle myosin and its metabolic efficiency undergo alterations during hibernation to optimize energy utilization. We isolated muscle fibers from small hibernators, Ictidomys tridecemlineatus and Eliomys quercinus and larger hibernators, Ursus arctos and Ursus americanus. We then conducted loaded Mant-ATP chase experiments alongside X-ray diffraction to measure resting myosin dynamics and its ATP demand. In parallel, we performed multiple proteomics analyses. Our results showed a preservation of myosin structure in U. arctos and U. americanus during hibernation, whilst in I. tridecemlineatus and E. quercinus, changes in myosin metabolic states during torpor unexpectedly led to higher levels in energy expenditure of type II, fast-twitch muscle fibers at ambient lab temperatures (20°C). Upon repeating loaded Mant-ATP chase experiments at 8°C (near the body temperature of torpid animals), we found that myosin ATP consumption in type II muscle fibers was reduced by 77-107% during torpor compared to active periods. Additionally, we observed Myh2 hyper-phosphorylation during torpor in I. tridecemilineatus, which was predicted to stabilize the myosin molecule. This may act as a potential molecular mechanism mitigating myosin-associated increases in skeletal muscle energy expenditure during periods of torpor in response to cold exposure. Altogether, we demonstrate that resting myosin is altered in hibernating mammals, contributing to significant changes to the ATP consumption of skeletal muscle. Additionally, we observe that it is further altered in response to cold exposure and highlight myosin as a potentially contributor to skeletal muscle non-shivering thermogenesis.

Higher Expression of Denervation-responsive Genes is Negatively Associated with Muscle Volume and Performance Traits in the Study of Muscle, Mobility and Aging (SOMMA)

With aging skeletal muscle fibers undergo repeating cycles of denervation and reinnervation. In approximately the 8 th decade of life reinnervation no longer keeps pace, resulting in the accumulation of persistently denervated muscle fibers that in turn cause an acceleration of muscle dysfunction. The significance of denervation in important clinical outcomes with aging is poorly studied. The Study of Muscle, Mobility and Aging (SOMMA) is a large cohort study with the primary objective to assess how aging muscle biology impacts clinically important traits. Using transcriptomics data from vastus lateralis muscle biopsies in 575 participants we have selected 49 denervation-responsive genes to provide insights to the burden of denervation in SOMMA, to test the hypothesis that greater expression of denervation-responsive genes negatively associates with SOMMA participant traits that included time to walk 400 meters, fitness (VO 2peak ), maximal mitochondrial respiration, muscle mass and volume, and leg muscle strength and power. Consistent with our hypothesis, increased transcript levels of: a calcium-dependent intercellular adhesion glycoprotein (CDH15), acetylcholine receptor subunits (Chrna1, Chrnd, Chrne), a glycoprotein promoting reinnervation (NCAM1), a transcription factor regulating aspects of muscle organization (RUNX1), and a sodium channel (SCN5A) were each negatively associated with at least 3 of these traits. VO 2peak and maximal respiration had the strongest negative associations with 15 and 19 denervation-responsive genes, respectively. In conclusion, the abundance of denervation-responsive gene transcripts is a significant determinant of muscle and mobility outcomes in aging humans, supporting the imperative to identify new treatment strategies to restore innervation in advanced age.

Transcriptome Profiling Reveals Enhanced Mitochondrial Activity as a Cold Adaptive Strategy to Hypothermia in Zebrafish Muscle

The utilisation of synthetic torpor for interplanetary travel once seemed farfetched. However, mounting evidence points to torpor-induced protective benefits from the main hazards of space travel, namely, exposure to radiation and microgravity. To determine the radio-protective effects of an induced torpor-like state we exploited the ectothermic nature of the Danio rerio (zebrafish) in reducing their body temperatures to replicate the hypothermic states seen during natural torpor. We also administered melatonin as a sedative to reduce physical activity. Zebrafish were then exposed to low-dose radiation (0.3 Gy) to simulate radiation exposure on long-term space missions. Transcriptomic analysis found that radiation exposure led to an upregulation of inflammatory and immune signatures and a differentiation and regeneration phenotype driven by STAT3 and MYOD1 transcription factors. In addition, DNA repair processes were downregulated in the muscle two days' post-irradiation. The effects of hypothermia led to an increase in mitochondrial translation including genes involved in oxidative phosphorylation and a downregulation of extracellular matrix and developmental genes. Upon radiation exposure, increases in endoplasmic reticulum stress genes were observed in a torpor+radiation group with downregulation of immune-related and ECM genes. Exposing hypothermic zebrafish to radiation also resulted in a downregulation of ECM and developmental genes however, immune/inflammatory related pathways were downregulated in contrast to that observed in the radiation only group. A cross-species comparison was performed with the muscle of hibernating Ursus arctos horribilis (brown bear) to define shared mechanisms of cold tolerance. Shared responses show an upregulation of protein translation and metabolism of amino acids, as well as a hypoxia response with the shared downregulation of glycolysis, ECM, and developmental genes.

Serum plays an important role in reprogramming the seasonal transcriptional profile of brown bear adipocytes

Understanding how metabolic reprogramming happens in cells will aid the progress in the treatment of a variety of metabolic disorders. Brown bears undergo seasonal shifts in insulin sensitivity, including reversible insulin resistance in hibernation. We performed RNA-sequencing on brown bear adipocytes and proteomics on serum to identify changes possibly responsible for reversible insulin resistance. We observed dramatic transcriptional changes, which depended on both the cell and serum season of origin. Despite large changes in adipocyte gene expression, only changes in eight circulating proteins were identified as related to the seasonal shifts in insulin sensitivity, including some that have not previously been associated with glucose homeostasis. The identified serum proteins may be sufficient for shifting hibernation adipocytes to an active-like state.

Revelations About Aging and Disease from Unconventional Vertebrate Model Organisms

Aging is a major risk factor for multiple diseases. Understanding the underlying mechanisms of aging would help to delay and prevent age-associated diseases. Short-lived model organisms have been extensively used to study the mechanisms of aging. However, these short-lived species may be missing the longevity mechanisms that are needed to extend the lifespan of an already long-lived species such as humans. Unconventional long-lived animal species are an excellent resource to uncover novel mechanisms of longevity and disease resistance. Here, we review mechanisms that evolved in nonmodel vertebrate species to counteract age-associated diseases. Some antiaging mechanisms are conserved across species; however, various nonmodel species also evolved unique mechanisms to delay aging and prevent disease. This variety of antiaging mechanisms has evolved due to the remarkably diverse habitats and behaviors of these species. We propose that exploring a wider range of unconventional vertebrates will provide important resources to study antiaging mechanisms that are potentially applicable to humans.

Integrated genomic and proteomic analyses identify stimulus-dependent molecular changes associated with distinct modes of skeletal muscle atrophy

Skeletal muscle atrophy is a debilitating condition that occurs with aging and disease, but the underlying mechanisms are incompletely understood. Previous work determined that common transcriptional changes occur in muscle during atrophy induced by different stimuli. However, whether this holds true at the proteome level remains largely unexplored. Here, we find that, contrary to this earlier model, distinct atrophic stimuli (corticosteroids, cancer cachexia, and aging) induce largely different mRNA and protein changes during muscle atrophy in mice. Moreover, there is widespread transcriptome-proteome disconnect. Consequently, atrophy markers (atrogenes) identified in earlier microarray-based studies do not emerge from proteomics as generally induced by atrophy. Rather, we identify proteins that are distinctly modulated by different types of atrophy (herein defined as “atroproteins”) such as the myokine CCN1/Cyr61, which regulates myofiber type switching during sarcopenia. Altogether, these integrated analyses indicate that different catabolic stimuli induce muscle atrophy via largely distinct mechanisms.

Skeletal muscle wasting is caused by many catabolic stimuli, which were thought to act via shared mechanisms. Hunt et al. now show that distinct catabolic stimuli induce muscle wasting via largely different molecular changes. The authors identify atrophy-associated proteins (“atroproteins”) that may represent diagnostic biomarkers and/or therapeutic targets.

Concurrent BMP Signaling Maintenance and TGF-β Signaling Inhibition Is a Hallmark of Natural Resistance to Muscle Atrophy in the Hibernating Bear

Muscle atrophy arises from a multiplicity of physio-pathological situations and has very detrimental consequences for the whole body. Although knowledge of muscle atrophy mechanisms keeps growing, there is still no proven treatment to date. This study aimed at identifying new drivers for muscle atrophy resistance. We selected an innovative approach that compares muscle transcriptome between an original model of natural resistance to muscle atrophy, the hibernating brown bear, and a classical model of induced atrophy, the unloaded mouse. Using RNA sequencing, we identified 4415 differentially expressed genes, including 1746 up- and 2369 down-regulated genes, in bear muscles between the active versus hibernating period. We focused on the Transforming Growth Factor (TGF)-β and the Bone Morphogenetic Protein (BMP) pathways, respectively, involved in muscle mass loss and maintenance. TGF-β- and BMP-related genes were overall down- and up-regulated in the non-atrophied muscles of the hibernating bear, respectively, and the opposite occurred for the atrophied muscles of the unloaded mouse. This was further substantiated at the protein level. Our data suggest TGF-β/BMP balance is crucial for muscle mass maintenance during long-term physical inactivity in the hibernating bear. Thus, concurrent activation of the BMP pathway may potentiate TGF-β inhibiting therapies already targeted to prevent muscle atrophy.

Landscape condition influences energetics, reproduction, and stress biomarkers in grizzly bears

Environmental change has been shown to influence mammalian distribution, habitat use, and behavior; however, few studies have investigated the impact on physiological function. This study aimed to determine the influence of landscape condition on the expression of target proteins related to energetics, reproduction, and stress in grizzly bears. We hypothesized that changes in landscape condition explains protein expression. Skin biopsies were collected from free-ranging grizzly bears in Alberta, Canada from 2013-2019 (n = 86 individuals). We used an information theoretic approach to develop 11 a priori candidate generalized linear mixed models to explain protein expression. We compared models using Akaike Information Criteria (AICc) weights and averaged models with ΔAICc < 2 for each protein. Food resources, represented by increased distance to coal mines and decreased crown closure, positively influenced energetic proteins (adiponectin and alpha-1-acid glycoprotein). Proteins related to reproduction (ceruloplasmin and serpin B5) were positively associated with increased wetland and upland food resources in addition to movement, but negatively associated with increased distance to roads. One stress related protein, complement C3, was positively influenced by increased percent conifer. Given the need to detect emerging threats to wildlife, we suggest the assessment of physiological function will lead to improved monitoring of species in rapidly changing landscapes.

Body Protein Sparing in Hibernators: A Source for Biomedical Innovation

Proteins are not only the major structural components of living cells but also ensure essential physiological functions within the organism. Any change in protein abundance and/or structure is at risk for the proper body functioning and/or survival of organisms. Death following starvation is attributed to a loss of about half of total body proteins, and body protein loss induced by muscle disuse is responsible for major metabolic disorders in immobilized patients, and sedentary or elderly people. Basic knowledge of the molecular and cellular mechanisms that control proteostasis is continuously growing. Yet, finding and developing efficient treatments to limit body/muscle protein loss in humans remain a medical challenge, physical exercise and nutritional programs managing to only partially compensate for it. This is notably a major challenge for the treatment of obesity, where therapies should promote fat loss while preserving body proteins. In this context, hibernating species preserve their lean body mass, including muscles, despite total physical inactivity and low energy consumption during torpor, a state of drastic reduction in metabolic rate associated with a more or less pronounced hypothermia. The present review introduces metabolic, physiological, and behavioral adaptations, e.g., energetics, body temperature, and nutrition, of the torpor or hibernation phenotype from small to large mammals. Hibernating strategies could be linked to allometry aspects, the need for periodic rewarming from torpor, and/or the ability of animals to fast for more or less time, thus determining the capacity of individuals to save proteins. Both fat- and food-storing hibernators rely mostly on their body fat reserves during the torpid state, while minimizing body protein utilization. A number of them may also replenish lost proteins during arousals by consuming food. The review takes stock of the physiological, molecular, and cellular mechanisms that promote body protein and muscle sparing during the inactive state of hibernation. Finally, the review outlines how the detailed understanding of these mechanisms at play in various hibernators is expected to provide innovative solutions to fight human muscle atrophy, to better help the management of obese patients, or to improve the ex vivo preservation of organs.

Nitrogen recycling buffers against ammonia toxicity from skeletal muscle breakdown in hibernating arctic ground squirrels

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 are unknown. Here we show, using pulse isotopic tracing, evidence of increased myofibrillar (skeletal muscle) protein breakdown and suppressed whole-body production of metabolites in vivo throughout deep torpor. As whole-body production of metabolites is suppressed, amino acids with nitrogenous side chains accumulate during torpor, while urea cycle intermediates do not. Using ¹⁵N stable isotope methodology in arctic ground squirrels (Urocitellus parryii), we provide evidence that free nitrogen is buffered and recycled into essential amino acids, non-essential amino acids and the gamma-glutamyl system during the inter-bout 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.

Specific shifts in the endocannabinoid system in hibernating brown bears

In small hibernators, global downregulation of the endocannabinoid system (ECS), which is involved in modulating neuronal signaling, feeding behavior, energy metabolism, and circannual rhythms, has been reported to possibly drive physiological adaptation to the hibernating state. In hibernating brown bears (Ursus arctos), we hypothesized that beyond an overall suppression of the ECS, seasonal shift in endocannabinoids compounds could be linked to bear’s peculiar features that include hibernation without arousal episodes and capacity to react to external disturbance. We explored circulating lipids in serum and the ECS in plasma and metabolically active tissues in free-ranging subadult Scandinavian brown bears when both active and hibernating. In winter bear serum, in addition to a 2-fold increase in total fatty acid concentration, we found significant changes in relative proportions of circulating fatty acids, such as a 2-fold increase in docosahexaenoic acid C22:6 n-3 and a decrease in arachidonic acid C20:4 n-6. In adipose and muscle tissues of hibernating bears, we found significant lower concentrations of 2-arachidonoylglycerol (2-AG), a major ligand of cannabinoid receptors 1 (CB1) and 2 (CB2). Lower mRNA level for genes encoding CB1 and CB2 were also found in winter muscle and adipose tissue, respectively. The observed reduction in ECS tone may promote fatty acid mobilization from body fat stores, and favor carbohydrate metabolism in skeletal muscle of hibernating bears. Additionally, high circulating level of the endocannabinoid-like compound N-oleoylethanolamide (OEA) in winter could favor lipolysis and fatty acid oxidation in peripheral tissues. We also speculated on a role of OEA in the conservation of an anorexigenic signal and in the maintenance of torpor during hibernation, while sustaining the capacity of bears to sense stimuli from the environment.

Cellular, Molecular, and Physiological Adaptations of Hibernation: The Solution to Environmental Challenges

Thriving in times of resource scarcity requires an incredible flexibility of behavioral, physiological, cellular, and molecular functions that must change within a relatively short time. Hibernation is a collection of physiological strategies that allows animals to inhabit inhospitable environments, where they experience extreme thermal challenges and scarcity of food and water. Many different kinds of animals employ hibernation, and there is a spectrum of hibernation phenotypes. Here, we focus on obligatory mammalian hibernators to identify the unique challenges they face and the adaptations that allow hibernators to overcome them. This includes the cellular and molecular strategies used to combat low environmental and body temperatures and lack of food and water. We discuss metabolic, neuronal, and hormonal cues that regulate hibernation and how they are thought to be coordinated by internal clocks. Last, we touch on questions that are left to be addressed in the field of hibernation research. Studies from the last century and more recent work reveal that hibernation is not simply a passive reduction in body temperature and vital parameters but rather an active process seasonally regulated at the molecular, cellular, and organismal levels.

Functional Studies with Primary Cells Provide a System for Genome-to-Phenome Investigations in Marine Mammals

Marine mammals exhibit some of the most dramatic physiological adaptations in their clade and offer unparalleled insights into the mechanisms driving convergent evolution on relatively short time scales. Some of these adaptations, such as extreme tolerance to hypoxia and prolonged food deprivation, are uncommon among most terrestrial mammals and challenge established metabolic principles of supply and demand balance. Non-targeted omics studies are starting to uncover the genetic foundations of such adaptations, but tools for testing functional significance in these animals are currently lacking. Cellular modeling with primary cells represents a powerful approach for elucidating the molecular etiology of physiological adaptation, a critical step in accelerating genome-to-phenome studies in organisms in which transgenesis is impossible (e.g., large-bodied, long-lived, fully aquatic, federally protected species). Gene perturbation studies in primary cells can directly evaluate whether specific mutations, gene loss, or duplication confer functional advantages such as hypoxia or stress tolerance in marine mammals. Here, we summarize how genetic and pharmacological manipulation approaches in primary cells have advanced mechanistic investigations in other non-traditional mammalian species, and highlight the need for such investigations in marine mammals. We also provide key considerations for isolating, culturing, and conducting experiments with marine mammal cells under conditions that mimic in vivo states. We propose that primary cell culture is a critical tool for conducting functional mechanistic studies (e.g., gene knockdown, over-expression, or editing) that can provide the missing link between genome- and organismal-level understanding of physiological adaptations in marine mammals.

Gut transcriptomic changes during hibernation in the greater horseshoe bat (Rhinolophus ferrumequinum)

Background

The gut is the major organ for nutrient absorption and immune response in the body of animals. Although effects of fasting on the gut functions have been extensively studied in model animals (e.g. mice), little is known about the response of the gut to fasting in a natural condition (e.g. hibernation). During hibernation, animals endure the long term of fasting and hypothermia.

Results

Here we generated the first gut transcriptome in a wild hibernating bat (Rhinolophus ferrumequinum). We identified 1614 differentially expressed genes (DEGs) during four physiological states (Torpor, Arousal, Winter Active and Summer Active). Gene co-expression network analysis assigns 926 DEGs into six modules associated with Torpor and Arousal. Our results reveal that in response to the stress of luminal nutrient deficiency during hibernation, the gut helps to reduce food intake by overexpressing genes (e.g. CCK and GPR17) that regulate the sensitivity to insulin and leptin. At the same time, the gut contributes energy supply by overexpressing genes that increase capacity for ketogenesis (HMGCS2) and selective autophagy (TEX264). Furthermore, we identified separate sets of multiple DEGs upregulated in Torpor and Arousal whose functions are involved in innate immunity.

Conclusion

This is the first gut transcriptome of a hibernating mammal. Our study identified candidate genes associated with regulation of food intake and enhance of innate immunity in the gut during hibernation. By comparing with previous studies, we found that two DEGs (CPE and HSPA8) were also significantly elevated during torpor in liver and brain of R. ferrumequinum and several DEGs (e.g. TXNIP and PDK1/4) were commonly upregulated during torpor in multiple tissues of different mammals. Our results support that shared expression changes may underlie the hibernation phenotype by most mammals.