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

Laura Cussonneau, Christian Boyer, Charlotte Brun, Christiane Deval, Emmanuelle Loizon, Emmanuelle Meugnier, Elise Gueret, Emeric Dubois, Daniel Taillandier, Cécile Polge, et al.

► To cite this version:

Laura Cussonneau, Christian Boyer, Charlotte Brun, Christiane Deval, Emmanuelle Loizon, et al. Concurrent BMP Signaling Maintenance and TGF-β Signaling Inhibition Is a Hallmark of Natural Resistance to Muscle Atrophy in the Hibernating Bear. Cells, 2021, 10 (8), pp.1873. 10.3390/cells10081873. hal-03314541

HAL Id: hal-03314541

https://hal.inrae.fr/hal-03314541

Submitted on 8 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

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

Laura Cussonneau<fnref n="1" />,*, Charlotte Boyer<fnref n="1" />, Charlotte Brun<fnref n="2" />, Christiane Deval<fnref n="1" />, Emmanuelle Loizon<fnref n="3" />, Emmanuelle Meugnier<fnref n="3" />, Elise Gueret 4,5, Daniel Taillandier<fnref n="1" />, Cécile Polge<fnref n="1" />, Daniel Béchet<fnref n="1" />, Guillemette Gauquelin-Koch<fnref n="6" />, Alina L. Evans<fnref n="7" />, Jon M. Arnemo<fnref n="7" />,<fnref n="8" />, Jon E. Swenson<fnref n="9" />, Stéphane Blanc<fnref n="2" />, Chantal Simon<fnref n="3" />, Etienne Lefai<fnref n="1" />, Fabrice Bertile<fnref n="2" /> and Lydie Combaret<fnref n="1" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Cussonneau, L.; Boyer, C.; Brun, C.; Deval, C.; Loizon, E.; Meugnier, E.; Gueret, E.; Dubois, E.; Taillandier, D.; Polge, C.; et al.
Concurrent BMP Signaling Maintenance and TGF-β Signaling Inhibition Is a Hallmark of Natural Resistance to Muscle Atrophy in the Hibernating Bear. *Cells* **2021**, *10*, 1873.
https://doi.org/10.3390/cells10081873

Academic Editors: Maria Pennuto
and Marco Pirazzini

Received: 16 June 2021

Accepted: 20 July 2021

Published: 23 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** 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.

**Keywords:** brown bear hibernation; mouse unloading; muscle atrophy; physical inactivity; RNA sequencing; TGF-β/BMP signaling

<footnote n="1">¹ INRAE, Unité de Nutrition Humaine, Université Clermont Auvergne, UMR 1019, F-63000 Clermont-Ferrand, France; christian.boyer@inrae.fr (C.B.); christiane.deval@inrae.fr (C.D.); daniel.taillandier@inrae.fr (D.T.); cecile.polge@inrae.fr (C.P.); daniel.bechet@inrae.fr (D.B.); etienne.lefai@inrae.fr (E.L.)</footnote>

<footnote n="2">² Université de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France; charlotte.brun@etu.unistra.fr (C.B.); stephane.blanc@iphc.cnrs.fr (S.B.); fbertile@unistra.fr (F.B.)</footnote>

<footnote n="3">³ CarMen Laboratory, INSERM 1060, INRAE 1397, University of Lyon, F-69600 Oullins, France; emmanuelle.loizon@inserm.fr (E.L.); emmanuelle.meugnier@univ-lyon1.fr (E.M.); chantal.simon@univ-lyon1.fr (C.S.)</footnote>

<footnote n="4">⁴ Institut de Génomique Fonctionnelle (IGF), University Montpellier, CNRS, INSERM, 34094 Montpellier, France; elise.gueret@mgx.cnrs.fr (E.G.); emeric.dubois@mgx.cnrs.fr (E.D.)</footnote>

<footnote n="5">⁵ Montpellier GenomiX, France Génomique, 34095 Montpellier, France</footnote>

<footnote n="6">⁶ Centre National d'Etudes Spatiales, CNES, 75001 Paris, France; guillemette.gauquellinkoch@cnes.fr</footnote>

<footnote n="7">⁷ Department of Forestry and Wildlife Management, Inland Norway University of Applied Sciences, Campus Evenstad, NO-2480 Koppang, Norway; alina-evans@inn.no (A.L.E.); jon.arnemo@inn.no (J.M.A.)</footnote>

<footnote n="8">⁸ Department of Wildlife, Fish, and Environmental Studies, Swedish University of Agricultural Sciences, SE-901 83 Umeå, Sweden</footnote>

<footnote n="9">⁹ Faculty of Environmental Sciences and Natural Resource Management, Norwegian University of Life Sciences, NO-1432 Ås, Norway; jon.swenson@nmbu.no</footnote>

<footnote>* Correspondence: Laura.cussonneau@inrae.fr (L.C.); Lydie.combaret@inrae.fr (L.C.); Tel.: +(33)4-7362-4824 (Lydie Combaret)</footnote>
# 1. Introduction

Muscle atrophy is defined as a loss of muscle mass and strength and is associated with adverse health outcomes, such as an autonomy decline and an increase in morbidity and mortality in many catabolic conditions (e.g., cancer cachexia, heart, and kidney failure, fasting, sepsis, injury, aging, or physical inactivity, etc.) [1-5]. Given the increase in sedentary behavior and improvement in life expectancy, and with to date still no proven therapeutic or preventive treatment to date, muscle atrophy remains a major public health issue (World Health Organization data) [6]. Skeletal muscle tissue represents an important reservoir of amino acids, which are mobilized during catabolic situations to preserve vital functions, resulting in an imbalance of contractile protein turnover (i.e., proteolysis exceeding protein synthesis) [7,8]. Catabolic stimuli (e.g., oxidative stress, endoplasmic reticulum disturbances, nutrient shortage, mitochondrial disruptions, etc.) activate a complex network of intracellular modulators, which in turn lead to the activation of the ubiquitin-proteasome system (UPS) and autophagy [9,10]. These two main proteolytic systems in muscle tissue involve a set of genes, i.e., atrogenes, whose expression at the mRNA levels is commonly altered during atrophy [11]. Cascades of events and players of muscle atrophy are well described and conserved in mammals [12,13] and include the transforming growth factor-β (TGF-β) superfamily, with TGF-β signaling acting as a negative regulator, and Bone Morphogenetic Proteins (BMP) signaling as a positive regulator of muscle mass [14]. The TGF-β pathway mediates muscle atrophy through cytoplasmic and nuclear signaling molecules SMAD2/3, mainly leading to the expression of the atrogenes TRIM63 (MuRF1) and FBXO32 (atrogin-1) [15]. Constitutive expression of SMAD3 triggers muscle wasting, and inhibition of SMAD2 and SMAD3 is sufficient to induce muscle growth in vivo [16-19]. Conversely, the BMP pathway mediates muscle mass maintenance through cytoplasmic and nuclear signaling molecules SMAD1/5/9, promoting a negative transcriptional regulation for a ubiquitin ligase required for muscle wasting, FBXO30 (MUSA1) [20], and increasing the expression of BMP receptors activity in muscles induced hypertrophy through Smad1/5-mediated activation of mTOR signaling [21].

Hibernating bears (Ursidae family) are naturally resistant to muscle atrophy when facing the two major atrophic inducers, prolonged fasting and physical inactivity up to 5-7 months [22,23]. Conversely and during a shorter period, a loss of muscle mass and volume prevails in rodent models [24-29] and humans [4,5]. As in rodent models, the muscles of the active bear are sensitive to disuse after denervation, whereas the muscles of the hibernating brown bear (*Ursus arctos*) are resistant [30]. The hibernating brown bear, therefore, appears as a suitable model to study the underlying mechanisms of muscle mass maintenance [31]. How it withstands muscle loss in conditions where muscle atrophy is expected in non-hibernating mammals remains to be fully elucidated. However, several hypotheses can be raised; our recent analysis of muscle proteome in the hibernating brown bear revealed the maintenance of glycolysis and a turning down of ATP turnover [32]. In addition, we reported (i) a myogenic microRNA signature prone to promoting muscle regeneration and suppressing ubiquitin ligase expression in bear muscle during winter [33], as well as (ii) limited levels of oxidative stress [34]. To unravel the molecular basis of muscle maintenance at the mRNA level, the bear muscle transcriptome has already been explored using cDNA microarrays [35] or RNA sequencing [36,37]. These two transcriptomic studies suggested an overall reduction in energy and protein metabolism, consistent with metabolic suppression and lower energy demand in skeletal muscle during hibernation. Although these studies focused on the changes in bear muscle transcriptome during the hibernating versus active period, our study aimed to compare them with those that occur in the muscle transcriptome of disuse-induced atrophy in a mouse model. The rationale for such a comparative analysis of two contrasted situations of muscle atrophy or maintenance lies in the identification of potential new candidates, beyond already reported metabolic factors [35-37], that may help the hibernating bear resist atrophy, thereby providing new targets for fighting muscle atrophy in humans. Among the transcription factors involved in the regulation of the differentially expressed genes highlighted in bear muscle between
the hibernation and active periods, eight were involved in the regulation of the TGF-β superfamily. We therefore subsequently focused on an in-depth analysis of the TGF-β and BMP intracellular pathways.

**2. Materials and Methods**

2.1. *Animal Experiments*

2.1.1. Bear Sample Collection

Biopsies from the vastus lateralis muscle were collected from 17 free-ranging brown
bears, 2–3 years old (Ursus arctos; 11 females and 6 males), from Dalarna and Gävleborg
counties, Sweden, from 2014 to 2019 (Table S6). The samples were immediately frozen
on dry ice until storage at –80 °C. In a given year, the same bears were captured during
winter hibernation (February) and recaptured during their active period (June). The study
was approved by the Swedish Ethical Committee on Animal Experiment (applications Dnr
C3/2016 and Dnr C18/2015), the Swedish Environmental Protection Agency (NV-0741-18),
and the Swedish Board of Agriculture (Dnr 5.2.18–3060/17). All procedures complied
with Swedish laws and regulations. Capture, anesthesia, and sampling were carried out
according to an established biomedical protocol [38].

2.1.2. Mouse Model of Hindlimb Unloading

Our objective was to compare the muscle transcriptome of the hibernating bear with that in a rodent model of long-term physical inactivity. We chose the 10-day mouse model of unloading as an established disuse-atrophy model, with atrophic pathways still activated [28,39]. All experiments were conducted with the approval of the regional ethics committee (agreement n° D6334515) following the European Directive 2010/63/EU on the protection of vertebrate animals used for experimental and scientific purposes. This study was performed with 12 C57BL6/J adult male mice purchased from Janvier Labs (Le Genest-Saint-Isle, France). They were housed by pairs upon arrival in a polycarbonate cage in a controlled room (22 ± 2 °C, 60% ± 5% humidity, 12 h light/dark cycle, light period starting at 8 h), fed ad libitum, and given free access to water.

After 10 days of acclimatization, the mice were either kept unsuspended (Control, *n* = 6) or subjected to hindlimb unloading through tail suspension (Unloaded, *n* = 6) for 10 days. Custom tail suspension cages were adapted from previous studies [40,41]. The cages (43 × 29 × 24 cm) have an overhead frame to which two suspension systems are fixed in parallel on the width of the cage. These suspension systems are widely spaced in a cage so that mice can always be housed in pairs without touching each other. Unloaded mice had a metal ring attached near the base of their tails using surgical adhesive tape. This ring was then attached to a swivel that allowed a 360-degree rotation and was attached on a rail that covered the upper width of the cage. The height of the swivel has been adjusted to keep the mouse at an angle of about 30° from the head so that the hindlimbs could not touch the ground or the walls. During the 10 days of unloading, mice showed only a very small body weight loss (<7%) that occurred within the first 3 days, with no change in food intake. At the end of the experiment, soleus muscles were rapidly dissected out and immediately frozen in liquid nitrogen and stored at –80 °C until analyses. As for the data used from the RNA sequencing of Zhang et al. [42], the soleus muscle atrophied by 37% (from 7.24 ± 0.27 mg in control mice to 4.58 ± 0.31 mg in unloaded mice, *p* < 0.001 according to the unpaired Student's *t*-test).

2.2. RNA Sequencing of Brown Bear Muscle

2.2.1. RNA Isolation

Total RNA from bear muscles was isolated as described [43]. Briefly, muscle RNA
from six bears (paired samples collected in summer and winter in a given year for the same
individual) was extracted using TRIzol reagent (Invitrogen, Courtaboeuf, France).
### 2.2.2. Illumina RNA Sequencing, Data Assembly, Statistical Analysis

We constructed RNA-Seq libraries with the Truseq stranded mRNA sample preparation kit from Illumina and sequenced them in two lanes on an Illumina HiSeq2500 (single-end, 50 bp, six libraries per lane). Image analyses and base calling were performed using the HiSeq Control Software (v2.2.70, Illumina, San Diego, CA, USA) and Real-Time Analysis component (v1.18.66.4, Illumina, San Diego, CA, USA). Demultiplexing was performed using Illumina's conversion software (bcl2fastq 2.20). The quality of the raw data were assessed using FastQC from the Babraham Institute and the Illumina software SAV (Sequencing Analysis Viewer, Illumina, San Diego, CA, USA). A splice junction mapper, TopHat 2.1.1 [44] (using Bowtie 2.3.5.1 [45], Johns Hopkins University, MD, USA), was used to align the RNA-Seq reads to the *Ursus arctos* genome (GCA_003584765.1 ASM358476v1 assembly downloaded from NCBI) with a set of gene model annotations (GCF_003584765.1 ASM358476v1_genomic.gff downloaded on 17 June 2019, from NCBI). Final read alignments having more than three mismatches were discarded. Samtools (v1.9) (http://samtools.sourceforge.net) was used to sort the alignment files. Then, the gene quantification was performed with Featurecounts 1.6.2 (http://subread.sourceforge.net/) [46]. As the data were from a strand-specific assay, the read had to be mapped to the opposite strand of the gene (-s 2 option). Before statistical analysis, genes with less than 30 reads (cumulating over all the analyzed samples) were filtered out. Differentially expressed genes were identified using the Bioconductor(https://bioconductor.org/) [47] package DESeq2 1.26.0 [48] (R version 3.6.1 https://www.r-project.org/). Data were normalized using the DESeq2(https://bioconductor.org/packages/release/bioc/html/DESeq2.html, accessed on 16 June 2021) normalization method. Genes with an adjusted *p*-value below 5% (according to the Benjamini–Hochberg procedure that controls the FDR) were declared differentially expressed.

### 2.2.3. Functional and Pathway Enrichment Analysis

Hierarchical clustering of bear transcriptomic data (log-transformed) was performed using Cluster v3.0 software (University of Tokyo, Tokyo, Japan) from the 13531 reads [49]. Parameters were set as follows: median centering and normalization of genes for adjusting data and centroid linkage clustering for both genes and arrays. Dendrograms were generated and viewed using the Java Treeview v1.3.3 program (Alok Saldanha, Stanford University, Stanford, CA, USA) [50]. To identify the differentially expressed genes (DEGs), we selected a fold change (FC) Winter/Summer >1.3 or <1.0.77 and an adjusted *p*-value <1.01e-05 as cut-off standards, for the up- and down-regulated genes, respectively. Visualization of functional enrichment was performed using Metascape [51], a web-based portal for visualizing the inference of enriched biological pathways among the DEGs. For the given DEGs gene list, pathway and process enrichment analysis has been carried out with the following ontology sources: KEGG Pathway, GO Biological Processes, Reactome Gene Sets, Canonical Pathways, CORUM, TRRUST, DisGeNET, PaGenBase, Transcription Factor Targets, WikiPathways, PANTHER Pathway, and COVID. All genes in the genome have been used as the enrichment background. Terms with a *p*-value < 0.01, a minimum count of 3, and an enrichment factor >1.5 (the enrichment factor is the ratio between the observed counts and the counts expected by chance) are collected and grouped into clusters based on their membership similarities. More specifically, *p*-values are calculated based on the accumulative hypergeometric distribution, and *q*-values are calculated using the Benjamini–Hochberg procedure to account for multiple testings. Kappa scores are used as the similarity metric when performing hierarchical clustering on the enriched terms, and sub-trees with a similarity of >0.3 are considered a cluster. The most statistically significant term within a cluster is chosen to represent the cluster. The 10 top-score enrichment terms from that analysis are shown in Figure 1b,c, and the 10 top-score enrichment transcription factors regulating the DEGs are shown in Figure 1d. The heat map representing the expression changes of the TGF-β/BMP gene sets in bear versus mouse muscles was made using the PHeatmap package (R 1.4.1106, University of Tartu, Tartu, Estonia). Briefly,
the gene hierarchical clustering is based on Euclidean distance calculated from the log2FC value (Winter/Summer and Unloaded/Control for bear and mouse muscles, respectively).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Deep changes in brown bear muscle transcriptome during hibernation. (a) Heatmap from brown bear muscle (vastus lateralis) transcripts (n = 6 bears/season, the same individuals were sampled and analyzed in summer and winter); red indicates high and green indicates low expression level of the 13531 genes. (b) Graph representing the 10 top-score of significantly enriched</figcaption></figure>
terms in winter compared to summer, from the 1746 up-regulated differentially expressed genes (DEGs) (FC W/S > 1.3 | and padj < 10.01 |) or (c) from the 2369 down-regulated DEGs (FC W/S < 10.77 | and padj < 10.01 |; Table S1). The color code for the functional cluster is indicated in the respective graphs, and the bold numbers in the different bars represent the numbers of DEGs found in the enriched terms. (d) Graph representing the 10 top-score of transcription factors involved in DEGs regulation from "Formation of a pool of free 40S subunits" and "Extracellular matrix organization" enriched terms. The bold TFs are involved in TGF-β superfamily regulation.

## 2.3. Transcriptomic Data Assembly and Statistical Analysis of Mouse Muscle

We used transcriptomic data from an already published study [42]. Briefly, in this study, C57BL6/J adult male mice were either kept unsuspended (Control, *n* = 4) or subjected to hindlimb unloading through tail suspension (Unloaded, *n* = 4) for 10 days. The fastq files of eight soleus muscles were downloaded from GEO (GSE102284). A splice junction map- per, TopHat 2.1.1(Johns Hopkins University, MD, USA) [44], was used to align the RNA-Seq reads to the mouse genome (UCSC mm10) with a set of gene model annotations (genes.gtf downloaded from UCSC on 29 October 2019; GeneIDs come from the NCBI gene2refseq file). Final read alignments having more than three mismatches were discarded. Samtools (v1.9, http://www.htslib.org/) was used to sort the alignment files. Then, the gene quan- tification was performed with Featurecounts 2.0.0 (http://subread.sourceforge.net/) [46]. As the data were from a strand-specific assay, the read had to be mapped to the opposite strand of the gene (-s 2 option). Before statistical analysis, genes with less than 20 reads (cumulating all the analyzed samples) were filtered out. Differentially expressed genes were identified using the Bioconductor [47] package DESeq2 1.26.0 [48] as previously described (cf. 2.2.2).

## 2.4. Western Blot

Vastus lateralis muscles from eleven bears (paired samples collected in summer and winter in a given year for the same individual; Table S6) and soleus muscles from 10-days control or unloaded mice (*n* = 6/group) (~30 mg) were used. Samples were homoge- nized using a polytron in 1 mL of an ice-cold buffer (10 mM Tris pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% Triton X-100, 0.5% Igepal CA630) containing inhibitors of proteases (Protease Inhibitor Cocktail) and phosphatases (1 mM Na3VO3, 10 mM NaF) (Sigma, Saint-Quentin-Fallavier, France). The homogenates were stirred for 1 h at 4 °C and then centrifuged at 10,000 g for 15 min at 4 °C. The resulting supernatants were then stored at −80 °C until use. The concentration of proteins was determined using the Bradford Protein Assay Kit (Biorad, Marnes-la-Coquette, France). Proteins were then diluted in Laemmli buffer and stored at −80 °C until use. Protein extracts were subjected to SDS-PAGE (sodium dodecyl sulfate-polyacrylamide gel electrophoresis) using TGX™ FastCast™ 10% Acrylamide gels (Biorad, Marnes-la-Coquette, France) and transferred onto a PVDF membrane (Hybond P, Amersham, England). Blots were blocked for 1 h at room temperature with 5% bovine serum albumin in TBS buffer with 0.1% Tween-20 (TBS-T, pH = 7.8), then washed thrice in TBS-T and incubated (overnight, stirring, 4 °C) with appropriated primary antibodies against SMAD1/5 (PA5-80036, Thermo Fisher, Illkirch, France), SMAD2/3 (#8685, Cell Signaling Technology, Saint-Cyr-L'Ecole, France), SMAD4 (ab230815), CTGF (ab227180), and GDF5 (ab137698) (Abcam, Cambridge, United Kingdom). Blots were then washed and incubated for 1 h with an appropriate secondary horseradish peroxidase-conjugated antibody at room temperature. Signals were detected after incuba- tion with Luminata Crescendo Western HRP substrate (Millipore, Burlington, MA, USA) and visualized using G: BOX ChemiXT4 (XL1) imaging system (Syngene, Frederick, MD, USA). Signals were then quantified using the GeneTools software (Syngene, Cambridge, UK) and normalized against the total amount of proteins determined by TGX signals to correct for uneven loading. Protein data were presented as individual values. The bilateral ratio paired Student's *t*-test was used to compare the muscles of bears during summer and winter (S and W, respectively). For muscles of control and unloaded mice (C and U, respectively), statistical significance was determined using the bilateral unpaired Student's
t-test. Statistical analysis was performed using Prism 8 (GraphPad Prism 9, San Diego, CA, USA).

**3. Results**

**3.1. Deep Changes in Brown Bear Muscle Transcriptome during Hibernation**

The brown bear transcriptome data set revealed that, from 13531 transcripts commonly identified in all individuals, gene expression differed markedly between summer and winter (Figure 1a). We identified 4115 differentially expressed genes (DEGs) between muscles of the active and hibernating bear with mainly down-regulated genes (Table S1). The 10 top-score obtained from an annotation enrichment analysis performed from the down- and up-regulated DEGs highlighted several significant enriched terms regulated differentially in bear muscles between the two seasons (Figure 1b,c). For instance, the protein metabolism functional cluster was the most up-regulated in bear muscles in winter compared to summer, with "Formation of a pool of free 40S subunits", "Ribonucleoprotein complex biogenesis", or "Translation" enriched terms (Figure 1b). In addition, the tissue structure remodeling functional cluster was the most down-regulated in muscles of the hibernating versus active bear, with "Extracellular matrix organization", "Cell-substrate adhesion", or "Supramolecular fiber organization" enriched terms (Figure 1c). We then run a transcriptional regulatory network analysis to identify transcription factors (TFs) involved in the DEGs regulation from the two most differentially enriched terms between the two seasons (e.g., "Formation of a pool of free 40S subunits" and "Extracellular matrix organization") (Figure 1d). We found that SP1 was the TF the most involved in the DEGs regulation, as well as others TFs such as SP3, NFKB1, JUN, RELA, TFAP2A, ETS1, or SMAD4. Interestingly, these TFs cited above are all involved in the regulation or signal transduction of the TGF-β superfamily [52–54]. We therefore decided to focus on that superfamily, including (i) TGF-β signaling, which is a master regulator of the extracellular matrix organization and also involved in muscle mass loss, and (ii) BMP signaling, which has recently been discovered to be involved in muscle mass maintenance [20].

**3.2. Hibernation Induces a Transcriptional Shift from the TGF-β to the BMP Pathway**

From a thorough analysis of the literature [55–57], we have drawn up a list of the actors and regulators of the TGF-β superfamily. We then analyzed precisely how they were regulated at the mRNA level in the muscles of the hibernating bear between the winter and summer seasons.

The expression levels of two main TGF-β ligands, INHBA and MSTN, were dramatically lower (Fold change (FC) = 0.28 and 0.52, respectively) in winter, whereas INHBB expression was higher (FC = 1.85) (Figure 2, Tables S2 and S3). In BMP signaling, the main ligand described in muscle mass maintenance, GDF5, showed higher levels (FC = 2.5) in hibernating bear muscles compared to active muscles. Extracellular actors inhibiting (KCP, DCN, MGP, NOV, or CHRD) or promoting (CCN2 or BMPER) TGF-β and/or BMP signals were mainly down-regulated in winter compared to summer. Receptors from the TGF-β signaling were differentially expressed during hibernation, with ACVR1C and TGBBR2 levels being considerably lower (FC = 0.28 and 0.80, respectively) in winter compared to summer, whereas TGFBR1 and ACVR1B were higher (FC = 1.42 and 1.47, respectively). The GDF5 receptor, BMPR1B, was higher in winter compared to summer (FC = 1.37). The co-receptors that control intensity and specificity of the downstream TGF-β/BMP signaling were mainly down-regulated or unchanged in winter for both pathways, except for MUSK, an important BMP co-receptor in muscle cells, which was up-regulated (FC = 2.39). Overall, actors involved in the initiation of the TGF-β signal were mainly repressed while those driving the BMP signal were increased.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Deep transcriptomic reprogramming of TGF-β and BMP pathways in muscle during brown bear hibernation. Scheme showing brown bear vastus lateralis muscle transcripts involved in TGF-β and BMP signaling and depicting (i) their relationships [55–57] and (ii) the difference in their expression levels between hibernation and activity periods. Red and green boxes indicate, respectively up- and down-regulated genes during hibernation compared to the summer season, and white boxes indicate unchanged genes. Target genes of the TGF and BMP pathways are indicated in italic and are in green when down-regulated, in red when up-regulated, and in black when unchanged. Dashed lines show the SMAD-independent pathway and full lines the canonical signaling pathway. Arrows indicate activation, and ⊥ bars indicate inhibition. (n = 6 bears/season, the same individuals were sampled and analyzed in summer and winter, padj < 10^5`; Tables S2 and S3). ECAg: Extracellular Agonist, ECAn: Extracellular Antagonist, L: Ligand, Co-R: Co-Receptor, RII: Receptor type II, RI: Receptor type I, TA: Transcriptional Activator, TR: Transcriptional Repressor, SBEs: SMAD Binding Element. Created with BioRender.com.</figcaption></figure>

For TGF-β signaling, intracellular inhibitors such as ERBIN, LDLRAD4, EIF3I, STK17B, and PP2CA were up-regulated in winter bear muscles (FC = 1.26, 1.50, 1.37, 1.86, and 1.56, respectively), whereas some of the actors promoting the signal were lower expressed, i.e., DAB2 and TRAP1 (FC = 0.50 and 0.63, respectively). By contrast, for BMP signaling, intracellular inhibitors were mainly lower expressed in muscles of the hibernating bear, such as CTDNEP1, FKBP1A (FC = 0.70 and 0.67, respectively). Regarding the intracellular actors triggering TGF-β/BMP signaling, SMAD3 (TGF-β signaling; FC = 0.76) was less expressed, SMAD1 and SMAD5 (BMP signaling; FC = 1.99 and 1.30, respectively), and SMAD4 (common to TGF-β and BMP signaling; FC = 1.34) were more expressed in muscles of the hibernating bear compared to the active one. Thus, expression changes of the intracellular actors again suggest repression of TGF-β signaling but maintenance of the BMP signaling.

For nuclear components, transcriptional activators were either up- (FOXO3 and SP3, FC = 1.68 and 1.70) or down-regulated (e.g., KAT2B, ATF3 and ETS1, FC = 0.55, 0.46, and 0.54) for TGF-β, whereas mainly up-regulated in winter for BMP with YAP1, ZCCHC12, and HOXC8 (FC = 1.63, 2.21, and 2.46, respectively). Conversely, the transcriptional repressors were mainly up-regulated for the TGF-β pathway, i.e., TRIM33, YAP1, and SIRT1 (FC = 1.60, 1.63, and 1.82, respectively), whereas unchanged or lower expressed, i.e.,
TOB1 (FC = 0.78) for BMP during hibernation. Considering TGF-β target genes, an overall down-regulation was observed in winter compared to summer, as highlighted for the different collagen isoforms or the multiple metalloproteinases COL1A1/2 (FC = 0.04 and 0.07), COL3A1 (FC = 0.06), COL5A2 (FC = 0.28), COL6A1/3 (FC = 0.22 and 0.20), COL14A1 (FC = 0.14), and MMP2/14 (FC = 0.33 and 0.46) (Figure 2 and Table S2). For BMP target genes, it was less contrasted, with either an unchanged (e.g., RUNX2 or ID4), down-regulated (e.g., ID1 or ID2, FC = 0.31 and 0.65), or up-regulated expression (RGS4 or KLF10, FC = 1.72 and 1.32) during hibernation, with RGS4 being a muscle-specific gene. Overall, this supports a general down-regulation of the transcriptional activity that drives TGF-β signaling, with possible maintenance of that for BMP signaling.

Finally, TGF-β and BMP signaling also use a shared SMAD-independent pathway, involving a branch of the MAPK (Mitogen-Activated Protein Kinase) pathway [58]. In this pathway, the expression of TRAF6 and its downstream actor MAP3K7 were higher in muscles of the hibernating bear compared to the active one (FC = 1.78 and 1.60, respectively). Some of the TRAF6 downstream actors, including MEF2A and MEF2C, which are key muscle transcription factors, were as well up-regulated (FC = 1.89 and 2.23) in winter compared to summer bear muscles (Figure 2).

TGF-β and BMP signaling pathways are tightly regulated by several UPS members, such as E3 ubiquitin ligases (E3s) and deubiquitinating enzymes (DUB) [55,59]. TGF-β inhibitors localized from the receptors to the nucleus level were up-regulated in winter, such as for CUL1 (FC = 1.85), NEDD4L (FC = 1.63), and TRIM33 (FC = 1.60), the latter being also a positive regulator of BMP signal transduction (Figure 3 and Table S4). SMURF1, an E3 ligase inhibiting both TGF-β and BMP signaling, was also up-regulated (FC = 1.64) in muscles of the hibernating bear. For UPS activators of the TGF-β signaling pathway alone, several DUBs were up-regulated in winter, e.g., UCHL5, USP11, and USP4 (FC = 1.78, 1.77, and 1.37, respectively), the latter also promoting BMP signaling, and others were down-regulated such as for USP15 and TRAF4 (FC = 0.52 and 0.47). The TGF-β and BMP pathways regulate the transcription of some muscle-specific E3s (TRIM63, FBXO32, and FBXO30). None of them were differentially expressed in winter compared to summer (Figure 3 and Table S4).

Overall, this wide transcriptomic analysis highlighted a winter transcriptional pattern in muscles of the brown bear that was prone to shutting down TGF-β signaling while maintaining or even over-activating the BMP pathway.

### 3.3. Divergent Regulation of TGF-β and BMP Pathways in Atrophy-Resistant Muscles of the Hibernating Brown Bear versus Atrophied Muscles of the Unloaded Mouse

We compared the above-described brown bear muscle transcriptome to published transcriptomic data from a model of long-term physical inactivity in mice induced by 10 days of unloading, where the soleus muscle mass decreased by ~30% [42]. Comparing the two models, the muscle transcriptomic profiles appeared very different for selected gene expression related to TGF-β and BMP signaling pathways (Figure S1 and Tables S2–S4). Of note, gene expressions of TGF-β and BMP ligands were mainly differently regulated, as particularly evidenced for GDF5, which was up-regulated during bear hibernation, but down-regulated during unloading in the mouse model (FC = 0.18), and as well with MSTN strongly down-regulated during bear hibernation but unchanged in unloaded mouse muscles (FC = 1.57) (Figure 4a and Tables S2 and S3). Regarding receptor expressions, the two models responded quite similarly (Figure 4b and Tables S2 and S3), with three notable exceptions. Firstly, the gene expression of ACRV1C did not change in muscles of the unloaded mouse, unlike muscles of the hibernating bear, where a strong down-regulation occurred (FC = 1.11 and 0.28, respectively), what we also observed with the gene expression of TGFBR2 (FC = 1.08 and 0.80, respectively). Secondly, the gene expression of BMPR1B, GDF5 receptor, was up-regulated in muscles of the hibernating bear but was down-regulated in muscles of the unloaded mouse (FC = 0.66) (Figure 4b). The intracellular actors SMAD3, SMAD4, and SMAD1 were regulated similarly by bear hibernation or mouse unloading, although not to the same extent (FC = 0.73, 1.10, and 1.50, respectively).
However, SMAD2 was up-regulated (FC = 1.24), and SMAD5 and SMAD9 remained unchanged only in the muscles of the unloaded mouse (Figure 4c).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** A transcriptomic reprogramming of UPS components involved in TGF-β and BMP regulation prevails in muscle during brown bear hibernation. Scheme showing the E3s/DUBs enzymes regulation of brown bear vastus lateralis muscle transcripts involved in TGF-β and BMP signaling pathways and depicting (i) their relationships [12,60] and (ii) the difference in their expression levels between hibernation and activity periods. Red and green boxes indicate, respectively up- down- regulated genes during hibernation compared to the summer season, and white boxes indicate unchanged genes. Dashed lines show the SMAD-independent pathway and full lines the canonical signaling. Arrows indicate activation and bars inhibition. (n = 6 bears/season, the same individuals were sampled and analyzed in summer and winter, padj < 10^5); Table S4). SBEs: SMAD Binding Elements. Created with BioRender.com.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>(b) Receptors</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(c) Intracellular actors</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** The gene expression pattern of TGF-β and BMP components is different in brown bear muscle resistant to atrophy during hibernation compared to atrophied muscles of the unloaded mouse. Genes expression level in vastus lateralis muscle of</figcaption></figure>
active and hibernating brown bears (n = 6 bears/season, the same individuals were sampled and analyzed in summer and winter, log2FC Winter/Summer, dotted white bars) and in soleus muscle of control and unloaded mice (n = 4 mice per condition, log2FC Unloaded/Control, gray bars). The depicted genes were categorized as TGF-β and BMP (a) ligands, (b) receptors, and (c) intracellular actors. Data are expressed as log2FC ± IfcSE. Statistical significance is shown (* padj < |0.05|; ** padj < |0.01|; *** padj < |0.001|).

Regarding the E3s and DUB enzymes involved in the BMP pathway alone (Figure 5 upper panel), NEDD4 was similarly up-regulated in both models (FC = 1.63 and 1.28, respectively in bear and mouse), whereas FBXO30 (FC = 0.82) was down-regulated only in muscles of the unloaded mouse (Figure 5 and Table S4). The enzymes involved in the regulation of both TGF-β/BMP signaling were mainly up- or down-regulated in muscles of the hibernating bear but were mainly unaffected in unloaded mouse muscles (Figure 5 middle panel, Table S4). Finally, the E3s/DUB involved only in TGF-β signaling regulation were for half of them commonly unchanged or up-regulated, e.g., CUL1, UCHL5 or NEDD4L (FC = 1.17, 1.30 or 2.38, respectively), and for the other half oppositely regulated, e.g., SMURF2, CBLB, or TRIM62 (FC = 0.81, 1.50, and 1.97, respectively) in muscles of the unloaded mouse compared to muscles of the hibernating bear (Figure 5 lower panel and Table S4).

Overall, TGF-β and BMP signaling pathways were differentially regulated between atrophy-resistant muscles of the hibernating bear and atrophy-sensitive muscles of the unloaded mouse.

### 3.4. Hibernation Induces Changes in TGF-β and BMP Pathway Components at the Protein Level

To further compare these models of resistance and vulnerability to atrophy, we explored the protein levels of the SMAD intracellular actors. In the brown bear muscle, we observed a tendency to decrease for SMAD2 protein levels (p = 0.09) in winter compared to summer, whereas SMAD3 levels remained quite similar between the two seasons (Figure 6a,c,d). Protein levels of SMAD4 were higher in muscles of the hibernating bear compared to the active one (Figure 6a,e), but the converse was observed for SMAD1/5 (Figure 6a,f). As for mRNA, these SMAD proteins did not follow the same regulation pattern in muscles of the unloaded mouse, where none changed at the protein level (Figure 6b–f). The protein levels of CCN2, a TGF-β target gene that is also an extracellular activator of the TGF-β pathway and an inhibitor of the BMP pathway, were strongly lower in winter bear muscles but were unaffected in unloaded mouse muscles (Figure 6a,b,g). Finally, the protein levels for GDF5, a BMP ligand, remained unchanged in both muscles from the hibernating bear and the unloaded mouse (Figure 6a,b,h).
<figure><img src="image_10.png" /><figcaption>**Figure 5.** The gene expression pattern of muscle E3s/DUB enzymes involved in TGF-β and BMP pathway regulation differed in the atrophy-resistant brown bear muscle during hibernation compared to the atrophied muscle of the unloaded mouse. Genes expression level (a) in vastus lateralis muscle of active and hibernating brown bears (n = 6 bears/season, the same individuals were sampled and analyzed in summer and winter, log2FC Winter/Summer, dotted white bars), and (b) in soleus muscle of control and unloaded mice (n = 4 mice per condition, log2FC Unloaded/Control, gray bars). Data are expressed as log2FC ± IfcSE. Statistical significance is shown (* padj < |0.05|; ** padj < |0.01|; *** padj < |0.001|).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 6.** *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 6.** Hibernation induces changes in TGF-β and BMP components at protein level. Protein levels for total SMAD2/3, SMAD1/5, SMAD4, CCN2, and GDF5 were assessed by Western blots (a) in the vastus lateralis muscle of brown bears during summer (S) and winter (W) and (b) in the soleus muscle of control (C) and unloaded mice (U). Representative western blots are shown for three couples of bears and mice. (c–h). Data are presented as individuals’ values with mean bars (n = 11 bears/season, the same individuals were sampled and analyzed in summer and winter, and n = 6 mice per condition). Gray and black dots are for muscles of bears, in summer and winter, respectively, and gray and black triangles are for control and unloaded muscle of mice, respectively.</figcaption></figure>
**4. Discussion**

Although basic knowledge regarding the underlying mechanisms of muscle atrophy is continuously growing, essentially from rodent models and clinical studies in humans, there are still no efficient therapeutic strategies for its prevention and treatment. To explore new avenues, we compared a model of muscle atrophy resistance, the hibernating brown bear [23], and a mouse model of disuse-induced muscle atrophy [42]. Therefore, we analyzed the bear muscle transcriptome and identified sweeping changes in gene expression between the summer-active period and the winter-hibernating period, and we compared them with transcriptomic data from muscles of the unloaded mouse.

Whereas the loss of muscle mass during inactivity is often associated with a decrease in muscle protein synthesis [42,61–63], we reported here that genes implicated in protein metabolism were mainly up-regulated in muscles of the hibernating bear. This is consistent with previous studies that linked muscle atrophy resistance during hibernation to induction of protein translation [35,36,64]. Under activation of FOXOs transcription factors, atrogenes involved in both autophagy and UPS pathways are enhanced in rodents [65–67]. These atrogenes (e.g., MAP1LC3A, FBXO32, and ZFAND5) were indeed up-regulated or maintained in the unloaded mouse model but down-regulated or unchanged in muscles of the hibernating bear (Table S5). In agreement with previous studies, our data confirmed that proteolytic actors were up-regulated in disuse-induced muscle atrophy in rodents [11] but not in the bear model of muscle atrophy resistance [36]. Despite discrepancies between bear and mouse models (e.g., fed status, torpor vs. hindlimb muscle disuse), both display similar pathways controlling skeletal muscle mass and protein balance. For instance, the TGF-β signaling pathway has been reported to be evolutionarily conserved among several species, from *Caenorhabditis elegans* and *Drosophila melanogaster* to *Mus musculus* [68]. In addition, both models are responsive to denervation-induced atrophy when in active conditions [20,30]. We performed our analyses on the fast-twitch vastus lateralis muscle for the brown bear and the slow-twitch soleus muscle for the mouse. Thus, we cannot exclude the possibility that the differences recorded here between the models may partly have resulted from the specific nature of these muscles concerning their metabolic and contractile properties. However, it is noteworthy that, although fast-twitch muscles are not as sensitive to physical inactivity as slow-twitch muscles, both types of muscles atrophied in classical models of long-term physical inactivity in rodents [25,29,39,69,70] and humans [71–73]. Despite this difference in metabolic and contractile properties, a general down-regulation of genes involved in extracellular matrix (ECM) structure organization was observed in atrophied muscles from the unloaded mouse and atrophy-resistant muscles from the hibernating brown bear (Figure 1 and Table S2). This ECM structure remodeling is a common feature of other atrophic models reported during dystrophic diseases [74] or muscle disuse in rodents or humans [26,69,75–77].

TGF-β is currently of major interest within the field of skeletal muscle biology. Indeed, the gene disruption of two of its ligands Myostatin (MSTN) or Activin A (INHBA), and the inhibition of their shared receptor ActRIIB (ACVR2B) promote a profound muscle hypertrophy phenotype in various conditions and species [78–81]. On the contrary, over-expression of MSTN or INHBA leads to the recruitment and phosphorylation of SMAD2-3, triggering an atrophic transcriptional program [16,17,82–84]. Here, we reported that the expression levels of both ligands MSTN and INHBA were lower only in muscles of the hibernating bear. In addition, ACVR2B and SMAD2 mRNA levels were up-regulated in muscles of the unloaded mouse whereas maintained in muscles of the hibernating bear. At the protein level, SMAD2 was maintained in the muscles of the unloaded mouse and showed a tendency to decrease in muscles of the hibernating bear. However, despite extensive efforts and numerous antibodies tested, we could not characterize the SMADs phosphorylation status in bear muscles, which thus remains to be defined.

One TGF-β target gene, CCN2 (also known as CTGF), is an ECM protein associated with fibrotic activity that is up-regulated in several muscle chronic disorders (i.e., Duchenne muscular dystrophy or the amyotrophic lateral sclerosis) [85]. CCN2 is one of the main
pro-fibrotic cytokines acting downstream of TGF-β signaling and can amplify its effects through enhancement of TGF-β ligand-receptor binding [85–90]. Inhibition of CCN2 gene expression reduced fibrosis and improved muscle and locomotor performance in a rodent model of amyotrophic lateral sclerosis [86]. During unloading, mRNA and protein levels of CCN2 were unchanged in the muscles of the unloaded mouse but were strongly lower in the muscles of the hibernating brown bear. Taken together, our data strongly suggest that TGF-β signaling is overall inhibited only in muscles resistant to atrophy during bear hibernation.

The BMP pathway is a potent inducer of bone and cartilage formation [91]. BMP signaling was also recently discovered as a regulator of muscle mass, as its inhibition abolished the hypertrophic phenotype of the *MSTN*-KO mouse [20], and an increase in its receptor activity induced important muscle hypertrophy [21]. Regarding BMP ligands, GDF5 is essential to muscle mass maintenance, binding preferentially to the type I receptor BMPR1B. GDF5 expression was strongly induced in denervated mouse muscles, and its inhibition worsened muscle atrophy, suggesting a role in counteracting denervation-induced atrophy [20]. We report here that the gene expression for both *GDF5* and *BMPR1B* was strongly down-regulated in mouse muscles during unloading, whereas up-regulated in muscles of the hibernating bear. However, GDF5 protein levels were stable in the muscles of both the unloaded mouse and the hibernating bear. We recently demonstrated that circulating components of the hibernating bear serum were able to induce trans-species effects on human myotubes, notably inhibition of protein degradation. Therefore, we hypothesized that those components could be involved in the maintenance of muscle mass and strength in the hibernating bear [92]. Along with muscle mass maintenance, bear hibernation is also associated with bone mass maintenance [93], and thus GDF5 may constitute a possible target toward muscle and bone protection during long periods of physical inactivity and/or fasting. Unfortunately, the exploration of GDF5 concentration in bear serum was hampered by species cross-reactivity concerns with commercially available Elisa kits and thus remains to be addressed.

It has been proposed that the BMP and TGF-β common actor SMAD4 mainly engages with the TGF-β pathway, but switches to the BMP pathway when TGF-β transduction is reduced, and thus could be the limiting factor between these two signalings [20]. Moreover, denervation-induced muscle atrophy was exacerbated in the SMAD4 deficient mouse [20], and muscle mass was increased in humans with a mutation-associated gain of function in the SMAD4 gene [94]. We observed here higher mRNA and protein levels for SMAD4 in muscles of the hibernating bear, but only at mRNA levels in muscles of the unloaded mouse. This is concomitant with the overall down-regulation highlighted for the TGF-β signaling in the hibernating bear, which was not observed for the unloaded mouse. Thus, the inhibition of the TGF-β signaling in muscles of the hibernating bear may have released SMAD4 from TGF-β to BMP pathway to maintain muscle mass in a long period of disuse (Graphical abstract).

TGF-β and BMP share a SMAD-independent pathway that activates the E3 ubiquitin ligase TRAF6 [58]. In addition to its pro-atrophic role [95], TRAF6 is also required for myogenic differentiation and muscle regeneration via the MEF2 axis [96]. MEF2 is a conserved family of transcription factors involved in the control of muscle gene expression [97]. A recent muscle transcriptome analysis highlighted an inhibition of MEF2 transcription factors during human bed rest leading to skeletal muscle alterations [98]. Here, only *MEF2A* was up-regulated during unloading in mouse muscles, whereas *TRAF6*, *MEF2A*, and *MEF2C* were up-regulated in the hibernating bear muscles. We already observed a myogenic microRNA signature mediated by MEF2A signaling in the muscles of the hibernating bear, promoting mechanisms of muscle regeneration, suppression of ubiquitin ligases, and resistance to muscle atrophy [33]. Further studies are required to address whether this MEF2 signature could be under the control of the TGF-β and/or the BMP pathway through TRAF6.
**5. Conclusions**

Resistance to muscle atrophy in hibernating brown bears has so far been linked to a reduction in protein and energy metabolism. Here, we show for the first time that the TGF-β pathway is down-regulated whereas the BMP pathway is concomitantly sustained or even up-regulated in atrophy-resistant muscles of the hibernating brown bear. Thus, beyond strengthening the previous hypothesis of a hypometabolism enabling this natural resistance to muscle atrophy, our study provides new insights regarding the underlying mechanisms.

The originality of the current work lies in the choice to study the mechanisms involved in resistance to atrophy, and not solely, as in many studies, the mechanisms involved in the onset of atrophy. Our comparison of resistance and sensitivity to muscle atrophy animal models suggested the balance between the TGF-β and the BMP pathways as critical for preventing skeletal muscle atrophy over a long period of disuse. Many targeted therapies to counteract muscle atrophy already focus on TGF-β inhibition [99]. Our data open the way for further studies and clinical trials to test the effects of strategies to switch on (or sustain) the BMP pathway in combination with TGF-β inhibition to prevent disuse-induced muscle atrophy.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/cells10081873/s1, Figure S1, Tables S1–S6, and Blots quantification.

**Author Contributions:** Individual contributions of authors is as followed: Conceptualization, L.C. (Laura Cussonneau), E.L. (Etienne Lefai), F.B. and L.C. (Lydie Combaret); Methodology, L.C. (Laura Cussonneau), C.B. (Christian Boyer), C.B. (Charlotte Brun), C.D., E.L. (Etienne Lefai), E.M., A.L.E., J.M.A., J.E.S., E.L. (Emmanuelle Loizon), F.B. and L.C. (Lydie Combaret); Analysis, L.C. (Laura Cussonneau), E.G., E.D., E.L. (Etienne Lefai), F.B. and L.C. (Lydie Combaret); Writing—original draft preparation, L.C. (Laura Cussonneau), E.L. (Etienne Lefai), F.B. and L.C. (Lydie Combaret); Writing—review and editing, all authors; Project administration, E.L. (Etienne Lefai), F.B. and L.C. (Lydie Combaret); Funding acquisition, E.L. (Etienne Lefai), F.B. and L.C. (Lydie Combaret). All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the Center National d'Etudes Spatiales (CNES, #865, #974, #1006, #1905), the iSITE Challenge 3 Mobility program (Université Clermont Auvergne), and the CNRS and Strasbourg University (H2E project; MyoBears project of the PEPS ExoMod program). MGX acknowledges financial support from France Génomique National infrastructure, funded as part of "Investissement d'Avenir" program managed by the Agence Nationale pour la Recherche (contract ANR-10-INBS-09). L.C. (Laura Cussonneau) was supported by a grant from the Institut National de la Recherche Agronomique et Environnement and Clermont Métropole. C.B. (Christian Boyerand) and C.B. (Charlotte Brun) were supported by grants from the Ministère Français de l'Enseignement Supérieur, de la Recherche et de l'Innovation. The long-term funding of Scandinavian Brown Bear Research Project (SBBRP) has come primarily from the Swedish Environmental Protection Agency, the Norwegian Environment Agency, the Austrian Science Fund, and the Swedish Association for Hunting and Wildlife Management.

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki and of the European Directive 2010/63/EU and approved by the Institutional Review Board (or Ethics Committee) of (1) the Swedish Ethical Committee on Animal Experiment (applications Dnr C3/2016 and Dnr C18/2015), the Swedish Environmental Protection Agency (NV-0741-18), and the Swedish Board of Agriculture (Dnr 5.2.18–3060/17), and (2) the C2E2A (Comité d'Ethique pour l'Expérimentation Animale Auvergne) (D6334515-08/2018).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The analyzed transcriptomic mouse data that support the findings of this study are openly available in the GEO repository database at https://doi.org/10.1093/gerona/gly051, accessed on 16 June 2021, reference number (GSE102284). The generated transcriptomic bear data that support the findings of this study are openly available in the GEO repository database at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi, accessed on 16 June 2021, reference number
(GSE144856). Other data that support the findings of this study are available in the supplementary material of this article.

**Acknowledgments:** The authors thank the field capture team (D. Ahlqvist, A. Friebe, H. Nordin, H. Blomgren, and S. Persson) and the IEN (INRAE Clermont-Ferrand-Theix, France) for excellent assistance with mouse care (A. Cissoire, Y. Delorme, and M. Djelloul-Mazouz), and L. Parry for help during animal slaughtering. The authors are also grateful to T. Brioche (Unité Mixte de Recherche INRAe, Dynamique Musculaire et Métabolisme, Université de Montpellier, France) for his help in implementing the unloading model in mouse and F. Rocher for helpful discussions. This is paper N°317 of the Scandinavian Brown Bear Research Project.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Jamshidpour, B.; Bahrpeyma, F.; Khatami, M.-R. The Effect of Aerobic and Resistance Exercise Training on the Health Related Quality of Life, Physical Function, and Muscle Strength among Hemodialysis Patients with Type 2 Diabetes. *J. Bodyw. Mov. Ther.* **2019**, *24*, S1360859219303134. [CrossRef] [PubMed]

2. Leenders, M.; Verdijk, L.B.; van der Hoeven, L.; Adam, J.J.; van Kranenburg, J.; Nilwik, R.; van Loon, L.J.C. Patients With Type 2 Diabetes Show a Greater Decline in Muscle Mass, Muscle Strength, and Functional Capacity With Aging. *J. Am. Med. Dir. Assoc.* **2013**, *14*, 585–592. [CrossRef]

3. DeFronzo, R.A.; Tripathy, D. Skeletal Muscle Insulin Resistance Is the Primary Defect in Type 2 Diabetes. *Diabetes Care* **2009**, *32*, S157–S163. [CrossRef] [PubMed]

4. Alkner, B.A.; Tesch, P.A. Knee Extensor and Plantar Flexor Muscle Size and Function Following 90 Days of Bed Rest with or without Resistance Exercise. *Eur. J. Appl. Physiol.* **2004**, *93*, 294–305. [CrossRef] [PubMed]

5. Trappe, S.; Trappe, T.; Gallagher, P.; Harber, M.; Alkner, B.; Tesch, P. Human Single Muscle Fibre Function with 84 Day Bed-Rest and Resistance Exercise. *J. Physiol.* **2004**, *557*, 501–513. [CrossRef] [PubMed]

6. Thyfault, J.P.; Du, M.; Kraus, W.E.; Levine, J.A.; Booth, F.W. Physiology of Sedentary Behavior and Its Relationship to Health Outcomes. *Med. Sci. Sports Exerc.* **2015**, *47*, 1301–1305. [CrossRef] [PubMed]

7. Argilés, J.M.; Campos, N.; Lopez-Pedrosa, J.M.; Rueda, R.; Rodriguez-Mañas, L. Skeletal Muscle Regulates Metabolism via Interorgan Crosstalk: Roles in Health and Disease. *J. Am. Med. Dir. Assoc.* **2016**, *17*, 789–796. [CrossRef] [PubMed]

8. Powers, S.K.; Lynch, G.S.; Murphy, K.T.; Reid, M.B.; Zijdewind, I. Disease-Induced Skeletal Muscle Atrophy and Fatigue. *Med. Sci. Sports Exerc.* **2016**, *48*, 2307–2319. [CrossRef]

9. Bonaldo, P.; Sandri, M. Cellular and Molecular Mechanisms of Muscle Atrophy. *Dis. Model. Mech.* **2013**, *6*, 25–39. [CrossRef] [PubMed]

10. Lecker, S.H.; Goldberg, A.L.; Mitch, W.E. Protein Degradation by the Ubiquitin–Proteasome Pathway in Normal and Disease States. *J. Am. Soc. Nephrol.* **2006**, *17*, 1807–1819. [CrossRef]

11. Taillandier, D.; Polge, C. Skeletal Muscle Atrogenes: From Rodent Models to Human Pathologies. *Biochimie* **2019**, *166*, 251–269. [CrossRef] [PubMed]

12. Peris-Moreno, D.; Cussonneau, L.; Combaret, L.; Polge, C.; Taillandier, D. Ubiquitin Ligases at the Heart of Skeletal Muscle Atrophy Control. *Molecules* **2021**, *26*, 407. [CrossRef]

13. Vainshtein, A.; Sandri, M. Signaling Pathways That Control Muscle Mass. *Int. J. Mol. Sci.* **2020**, *21*, 4759. [CrossRef]

14. Sartori, R.; Gregorevic, P.; Sandri, M. TGFβ and BMP Signaling in Skeletal Muscle: Potential Significance for Muscle-Related Disease. *Trends Endocrinol. Metab.* **2014**, *25*, 464–471. [CrossRef] [PubMed]

15. Lokireddy, S.; Mouly, V.; Butler-Browne, G.; Gluckman, P.D.; Sharma, M.; Kambadur, R.; McFarlane, C. Myostatin Promotes the Wasting of Human Myoblast Cultures through Promoting Ubiquitin-Proteasome Pathway-Mediated Loss of Sarcomeric Proteins. *Am. J. Physiol.-Cell Physiol.* **2011**, *301*, C1316–C1324. [CrossRef] [PubMed]

16. Sartori, R.; Milan, G.; Patron, M.; Mammucari, C.; Blaauw, B.; Abraham, R.; Sandri, M. Smad2 and 3 Transcription Factors Control Muscle Mass in Adulthood. *Am. J. Physiol.-Cell Physiol.* **2009**, *296*, C1248–C1257. [CrossRef]

17. Tando, T.; Hirayama, A.; Furukawa, M.; Sato, Y.; Kobayashi, T.; Funayama, A.; Kanaji, A.; Hao, W.; Watanabe, R.; Morita, M.; et al. Smad2/3 Proteins Are Required for Immobilization-Induced Skeletal Muscle Atrophy. *J. Biol. Chem.* **2016**, *291*, 12184–12194. [CrossRef]

18. Goodman, C.A.; McNally, R.M.; Hoffmann, F.M.; Hornberger, T.A. Smad3 Induces Atrogin-1, Inhibits MTOR and Protein Synthesis, and Promotes Muscle Atrophy In Vivo. *Mol. Endocrinol.* **2013**, *27*, 1946–1957. [CrossRef]

19. Carlson, M.E.; Hsu, M.; Conboy, I.M. Imbalance between PSmad3 and Notch Induces CDK Inhibitors in Old Muscle Stem Cells. *Nature* **2008**, *454*, 528–532. [CrossRef] [PubMed]

20. Sartori, R.; Schirwis, E.; Blaauw, B.; Bortolanza, S.; Zhao, J.; Enzo, E.; Stantzou, A.; Mouisel, E.; Toniolo, L.; Ferry, A.; et al. BMP Signaling Controls Muscle Mass. *Nat. Genet.* **2013**, *45*, 1309–1318. [CrossRef]

21. Winbanks, C.E.; Chen, J.L.; Qian, H.; Liu, Y.; Bernardo, B.C.; Beyer, C.; Watt, K.I.; Thomson, R.E.; Connor, T.; Turner, B.J.; et al. The Bone Morphogenetic Protein Axis Is a Positive Regulator of Skeletal Muscle Mass. *J. Cell Biol.* **2013**, *203*, 345–357. [CrossRef] [PubMed]
22. Tinker, D.B.; Harlow, H.J.; Beck, T.D.I. Protein Use and Muscle-Fiber Changes in Free-Ranging, Hibernating Black Bears. *Physiol. Zool.* **1998**, *71*, 414–424. [CrossRef]

23. Lohuis, T.D.; Harlow, H.J.; Beck, T.D.I. Hibernating Black Bears (Ursus Americanus) Experience Skeletal Muscle Protein Balance during Winter Anorexia. *Comp. Biochem. Physiol. B Biochem. Mol. Biol.* **2007**, *147*, 20–28. [CrossRef]

24. Deval, C.; Capel, F.; Laillet, B.; Polge, C.; Béchet, D.; Taillandier, D.; Attaix, D.; Combaret, L. Docosahexaenoic Acid-Supplementation Prior to Fasting Prevents Muscle Atrophy in Mice: Docosahexaenoic Acid Limits Muscle Wasting in Fasted Mice. *J. Cachexia Sarcopenia Muscle* **2016**, *7*, 587–603. [CrossRef] [PubMed]

25. Vazeille, E.; Codran, A.; Claustre, A.; Averous, J.; Listrat, A.; Béchet, D.; Taillandier, D.; Dardevet, D.; Attaix, D.; Combaret, L. The Ubiquitin-Proteasome and the Mitochondria-Associated Apoptotic Pathways Are Sequentially Downregulated during Recovery after Immobilization-Induced Muscle Atrophy. *Am. J. Physiol. Endocrinol. Metab.* **2008**, *295*, E1181–E1190. [CrossRef] [PubMed]

26. Slimani, L.; Vazeille, E.; Deval, C.; Meunier, B.; Polge, C.; Dardevet, D.; Béchet, D.; Taillandier, D.; Micol, D.; Listrat, A.; et al. The Delayed Recovery of the Remobilized Rat Tibialis Anterior Muscle Reflects a Defect in Proliferative and Terminal Differentiation That Impairs Early Regenerative Processes. *J. Cachexia Sarcopenia Muscle* **2015**, *6*, 73–83. [CrossRef]

27. Gao, Y.; Arfat, Y.; Wang, H.; Goswami, N. Muscle Atrophy Induced by Mechanical Unloading: Mechanisms and Potential Countermeasures. *Front. Physiol.* **2018**, *9*, 235. [CrossRef] [PubMed]

28. Cui, Q.; Yang, H.; Gu, Y.; Zong, C.; Chen, X.; Lin, Y.; Sun, H.; Shen, Y.; Zhu, J. RNA Sequencing (RNA-Seq) Analysis of Gene Expression Provides New Insights into Hindlimb Unloading-Induced Skeletal Muscle Atrophy. *Ann. Transl. Med.* **2020**, *8*, 1595. [CrossRef] [PubMed]

29. Deval, C.; Calonne, J.; Coudy-Gandilhon, C.; Vazeille, E.; Bechet, D.; Polge, C.; Taillandier, D.; Attaix, D.; Combaret, L. Mitophagy and Mitochondria Biogenesis Are Differentially Induced in Rat Skeletal Muscles during Immobilization and/or Remobilization. *Int. J. Mol. Sci.* **2020**, *21*, 3691. [CrossRef] [PubMed]

30. Lin, D.C.; Hershey, J.D.; Mattoon, J.S.; Robbins, C.T. Skeletal Muscles of Hibernating Brown Bears Are Unusually Resistant to Effects of Denervation. *J. Exp. Biol.* **2012**, *215*, 2081–2087. [CrossRef]

31. Bertile, F.; Habold, C.; Le Maho, Y.; Giroud, S. Body Protein Sparing in Hibernators: A Source for Biomedical Innovation. *Front. Physiol.* **2021**, *12*, 634953. [CrossRef] [PubMed]

32. Chazarin, B.; Storey, K.B.; Ziemianin, A.; Chanon, S.; Plumel, M.; Chery, I.; Durand, C.; Evans, A.L.; Arnemo, J.M.; Zedrosser, A.; et al. Metabolic Reprogramming Involving Glycolysis in the Hibernating Brown Bear Skeletal Muscle. *Front. Zool.* **2019**, *16*, 12. [CrossRef] [PubMed]

33. Luu, B.E.; Lefai, E.; Giroud, S.; Swenson, J.E.; Chazarin, B.; Gauquelin-Koch, G.; Arnemo, J.M.; Evans, A.L.; Bertile, F.; Storey, K.B. MicroRNAs Facilitate Skeletal Muscle Maintenance and Metabolic Suppression in Hibernating Brown Bears. *J. Cell. Physiol.* **2019**, *235*, 3984–3993. [CrossRef] [PubMed]

34. Chazarin, B.; Ziemianin, A.; Evans, A.L.; Meugnier, E.; Loizon, E.; Chery, I.; Arnemo, J.M.; Swenson, J.E.; Gauquelin-Koch, G.; Simon, C.; et al. Limited Oxidative Stress Favors Resistance to Skeletal Muscle Atrophy in Hibernating Brown Bears (Ursus Arctos). *Antioxidants* **2019**, *8*, 334. [CrossRef] [PubMed]

35. Fedorov, V.B.; Goropashnaya, A.V.; Stewart, N.C.; Tøien, Ø.; Chang, C.; Wang, H.; Yan, J.; Showe, L.C.; Showe, M.K.; Barnes, B.M. Comparative Functional Genomics of Adaptation to Muscular Disuse in Hibernating Mammals. *Mol. Ecol.* **2014**, *23*, 5524–5537. [CrossRef]

36. Jansen, H.T.; Trojahn, S.; Saxton, M.W.; Quackenbush, C.R.; Evans Hutzenbiler, B.D.; Nelson, O.L.; Cornejo, O.E.; Robbins, C.T.; Kelley, J.L. Hibernation Induces Widespread Transcriptional Remodeling in Metabolic Tissues of the Grizzly Bear. *Commun. Biol.* **2019**, *2*, 336. [CrossRef]

37. Mugahid, D.A.; Sengul, T.G.; You, X.; Wang, Y.; Steil, L.; Bergmann, N.; Radke, M.H.; Ofenbauer, A.; Gesell-Salazar, M.; Balogh, A.; et al. Proteomic and Transcriptomic Changes in Hibernating Grizzly Bears Reveal Metabolic and Signaling Pathways That Protect against Muscle Atrophy. *Sci. Rep.* **2019**, *9*, 19976. [CrossRef] [PubMed]

38. Arnemo, J.M.; Evans, A.L. Biomedical Protocols for Free-Ranging Brown Bears, *Wolves, Wolverines and Lynx*; Inland Norway University of Applied Sciences: Elverum, Norway, 2017.

39. Chacon-Cabrera, A.; Gea, J.; Barreiro, E. Short- and Long-Term Hindlimb Immobilization and Reloading: Profile of Epigenetic Events in Gastrocnemius. *J. Cell. Physiol.* **2017**, *232*, 1415–1427. [CrossRef]

40. Anderson, M.J.; Diko, S.; Baehr, L.M.; Baar, K.; Bodine, S.C.; Christiansen, B.A. Contribution of Mechanical Unloading to Trabecular Bone Loss Following Non-Invasive Knee Injury in Mice. *J. Orthop. Res.* **2016**, *34*, 1680–1687. [CrossRef] [PubMed]

41. Pagano, A.F.; Demangel, R.; Brioche, T.; Jublanc, E.; Bertrand-Gaday, C.; Candau, R.; Dechesne, C.A.; Dani, C.; Bonnieu, A.; Py, G.; et al. Muscle Regeneration with Intermuscular Adipose Tissue (IMAT) Accumulation Is Modulated by Mechanical Constraints. *PLOS ONE* **2015**, *10*, e0144230. [CrossRef] [PubMed]

42. Zhang, X.; Trevino, M.B.; Wang, M.; Gardell, S.J.; Ayala, J.E.; Han, X.; Kelly, D.P.; Goodpaster, B.H.; Vega, R.B.; Coen, P.M. Impaired Mitochondrial Energetics Characterize Poor Early Recovery of Muscle Mass Following Hind Limb Unloading in Old Mice. *J. Gerontol. Ser. A* **2018**, *73*, 1313–1322. [CrossRef]

43. Boyer, C.; Cussonneau, L.; Brun, C.; Deval, C.; Pais de Barros, J.-P.; Chanon, S.; Bernoud-Hubac, N.; Daira, P.; Evans, A.L.; Arnemo, J.M.; et al. Specific Shifts in the Endocannabinoid System in Hibernating Brown Bears. *Front. Zool.* **2020**, *17*, 35. [CrossRef]

44. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate Alignment of Transcriptomes in the Presence of Insertions, Deletions and Gene Fusions. *Genome Biol.* **2013**, *14*, R36. [CrossRef] [PubMed]
45. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and Memory-Efficient Alignment of Short DNA Sequences to the Human Genome. *Genome Biol.* **2009**, *10*, R25. [CrossRef]

46. Liao, Y.; Smyth, G.K.; Shi, W. FeatureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features. *Bioinformatics* **2014**, *30*, 923–930. [CrossRef] [PubMed]

47. Gentleman, R.C.; Carey, V.J.; Bates, D.M.; Bolstad, B.; Dettling, M.; Dudoit, S.; Ellis, B.; Gautier, L.; Ge, Y.; Gentry, J.; et al. Bioconductor: Open Software Development for Computational Biology and Bioinformatics. *Genome Biol.* **2004**, *5*, 16.

48. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef] [PubMed]

49. De Hoon, M.J.L.; Imoto, S.; Nolan, J.; Miyano, S. Open Source Clustering Software. *Bioinformatics* **2004**, *20*, 1453–1454. [CrossRef] [PubMed]

50. Saldanha, A.J. Java Treeview–Extensible Visualization of Microarray Data. *Bioinformatics* **2004**, *20*, 3246–3248. [CrossRef]

51. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. *Nat. Commun.* **2019**, *10*, 1523. [CrossRef]

52. Koinuma, D.; Tsutsumi, S.; Kamimura, N.; Taniguchi, H.; Miyazawa, K.; Sunamura, M.; Imamura, T.; Miyazono, K.; Aburatani, H. Chromatin Immunoprecipitation on Microarray Analysis of Smad2/3 Binding Sites Reveals Roles of ETS1 and TFAP2A in Transforming Growth Factor β Signaling. *Mol. Cell. Biol.* **2009**, *29*, 172–186. [CrossRef]

53. Luo, K. Signaling Cross Talk between TGF-β/Smad and Other Signaling Pathways. *Cold Spring Harb. Perspect. Biol.* **2017**, *9*, a022137. [CrossRef]

54. Córdova, G.; Rochard, A.; Riquelme-Guzmán, C.; Cofré, C.; Scherman, D.; Bigey, P.; Brandan, E. SMAD3 and SP1/SP3 Transcription Factors Collaborate to Regulate Connective Tissue Growth Factor Gene Expression in Myoblasts in Response to Transforming Growth Factor β. *J. Cell. Biochem.* **2015**, *116*, 1880–1887. [CrossRef] [PubMed]

55. Hata, A.; Chen, Y.-G. TGF-β Signaling from Receptors to Smads. *Cold Spring Harb. Perspect. Biol.* **2016**, *8*, a022061. [CrossRef] [PubMed]

56. Itoh, S.; ten Dijke, P. Negative Regulation of TGF-β Receptor/Smad Signal Transduction. *Curr. Opin. Cell Biol.* **2007**, *19*, 176–184. [CrossRef]

57. Weiss, A.; Attisano, L. The TGFbeta Superfamily Signaling Pathway: TGFbeta Superfamily Signaling Pathway. *Wiley Interdiscip. Rev. Dev. Biol.* **2013**, *2*, 47–63. [CrossRef] [PubMed]

58. Zhang, Y.E. Non-Smad Signaling Pathways of the TGF-β Family. *Cold Spring Harb. Perspect. Biol.* **2017**, *9*, a022129. [CrossRef] [PubMed]

59. Miyazono, K.; Kamiya, Y.; Morikawa, M. Bone Morphogenetic Protein Receptors and Signal Transduction. *J. Biochem. (Tokyo)* **2010**, *147*, 35–51. [CrossRef]

60. Inui, N.; Sakai, S.; Kitagawa, M. Molecular Pathogenesis of Pulmonary Fibrosis, with Focus on Pathways Related to TGF-β and the Ubiquitin-Proteasome Pathway. *Int. J. Mol. Sci.* **2021**, *22*, 6107. [CrossRef]

61. Ferrando, A.A.; Lane, H.W.; Stuart, C.A.; Davis-Street, J.; Wolfe, R.R. Prolonged Bed Rest Decreases Skeletal Muscle and Whole Body Protein Synthesis. *Am. J. Physiol.-Endocrinol. Metab.* **1996**, *270*, E627–E633. [CrossRef]

62. Ferrando, A.A.; Tipton, K.D.; Bamman, M.M.; Wolfe, R.R. Resistance Exercise Maintains Skeletal Muscle Protein Synthesis during Bed Rest. *J. Appl. Physiol.* **1997**, *82*, 807–810. [CrossRef]

63. Idris, I.; Atherton, P.J. Human Skeletal Muscle Disuse Atrophy: Effects on Muscle Protein Synthesis, Breakdown, and Insulin Resistance—A Qualitative Review. *Front. Physiol.* **2016**, *7*, 10.

64. Fedorov, V.B.; Goropashnaya, A.V.; Tøien, Ø.; Stewart, N.C.; Gracey, A.Y.; Chang, C.; Qin, S.; Pertea, G.; Quackenbush, J.; Showe, L.C.; et al. Elevated Expression of Protein Biosynthesis Genes in Liver and Muscle of Hibernating Black Bears (*Ursus Americanus*). *Physiol. Genomics* **2009**, *37*, 108–118. [CrossRef]

65. Zhao, J.; Brault, J.J.; Schild, A.; Cao, P.; Sandri, M.; Schiaffino, S.; Lecker, S.H.; Goldberg, A.L. FoxO3 Coordinately Activates Protein Degradation by the Autophagic/Lysosomal and Proteasomal Pathways in Atrophying Muscle Cells. *Cell Metab.* **2007**, *6*, 472–483. [CrossRef]

66. Mammucari, C.; Milan, G.; Romanello, V.; Masiero, E.; Rudolf, R.; Del Piccolo, P.; Burden, S.J.; Di Lisi, R.; Sandri, C.; Zhao, J.; et al. FoxO3 Controls Autophagy in Skeletal Muscle In Vivo. *Cell Metab.* **2007**, *6*, 458–471. [CrossRef] [PubMed]

67. Sandri, M.; Sandri, C.; Gilbert, A.; Skurk, C.; Calabria, E.; Picard, A.; Walsh, K.; Schiaffino, S.; Lecker, S.H.; Goldberg, A.L. Foxo Transcription Factors Induce the Atrophy-Related Ubiquitin Ligase Atrogin-1 and Cause Skeletal Muscle Atrophy. *Cell* **2004**, *117*, 399–412. [CrossRef]

68. Fabian, D.K.; Fuentealba, M.; Dönertas, H.M.; Partridge, L.; Thornton, J.M. Functional Conservation in Genes and Pathways Linking Ageing and Immunity. *Immun. Ageing* **2021**, *18*, 23. [CrossRef]

69. Slimani, L.; Micol, D.; Amat, J.; Delcros, G.; Meunier, B.; Taillandier, D.; Polge, C.; Béchet, D.; Dardevet, D.; Picard, B.; et al. The Worsening of Tibialis Anterior Muscle Atrophy during Recovery Post-Immobilization Correlates with Enhanced Connective Tissue Area, Proteolysis, and Apoptosis. *Am. J. Physiol.-Endocrinol. Metab.* **2012**, *303*, E1335–E1347. [CrossRef]

70. Pigna, E.; Renzini, A.; Greco, E.; Simonazzi, E.; Fulle, S.; Mancinelli, R.; Moresi, V.; Adamo, S. HDAC4 Preserves Skeletal Muscle Structure Following Long-Term Denervation by Mediating Distinct Cellular Responses. *Skelet. Muscle* **2018**, *8*, 6. [CrossRef]
71. Leermakers, P.A.; Kneppers, A.E.M.; Schols, A.M.W.J.; Kelders, M.C.J.M.; Theije, C.C.; Verdijk, L.B.; Loon, L.J.C.; Langen, R.C.J.; Gosker, H.R. Skeletal Muscle Unloading Results in Increased Mitophagy and Decreased Mitochondrial Biogenesis Regulation. *Muscle Nerve* **2019**, *60*, 769–778. [CrossRef]

72. Dirks, M.L.; Wall, B.T.; van de Valk, B.; Holloway, T.M.; Holloway, G.P.; Chabowski, A.; Goossens, G.H.; van Loon, L.J.C. One Week of Bed Rest Leads to Substantial Muscle Atrophy and Induces Whole-Body Insulin Resistance in the Absence of Skeletal Muscle Lipid Accumulation. *Diabetes* **2016**, *65*, 2862–2875. [CrossRef] [PubMed]

73. Brocca, L.; Longa, E.; Cannavino, J.; Seynnes, O.; de Vito, G.; McPhee, J.; Narici, M.; Pellegrino, M.A.; Bottinelli, R. Human Skeletal Muscle Fibre Contractile Properties and Proteomic Profile: Adaptations to 3 Weeks of Unilateral Lower Limb Suspension and Active Recovery: Impact of ULLS and Active Rehabilitation on Human Muscle Fibres. *J. Physiol.* **2015**, *593*, 5361–5385. [CrossRef] [PubMed]

74. Gillies, A.R.; Lieber, R.L. Structure and Function of the Skeletal Muscle Extracellular Matrix: Skeletal Muscle ECM. *Muscle Nerve* **2011**, *44*, 318–331. [CrossRef]

75. Martinez, D.A.; Vailas, A.C.; Vanderby, R.; Grindeland, R.E. Temporal Extracellular Matrix Adaptations in Ligament during Wound Healing and Hindlimb Unloading. *Am. J. Physiol.-Regul. Integr. Comp. Physiol.* **2007**, *293*, R1552–R1560. [CrossRef]

76. Miller, T.A.; Lesniewski, L.A.; Muller-Delp, J.M.; Majors, A.K.; Scalise, D.; Delp, M.D. Hindlimb Unloading Induces a Collagen Isoform Shift in the Soleus Muscle of the Rat. *Am. J. Physiol.-Regul. Integr. Comp. Physiol.* **2001**, *281*, R1710–R1717. [CrossRef]

77. Urso, M.L.; Scrimgeour, A.G.; Chen, Y.-W.; Thompson, P.D.; Clarkson, P.M. Analysis of Human Skeletal Muscle after 48 h Immobilization Reveals Alterations in mRNA and Protein for Extracellular Matrix Components. *J. Appl. Physiol.* **2006**, *101*, 13. [CrossRef]

78. Elkina, Y.; von Haehling, S.; Anker, S.D.; Springer, J. The Role of Myostatin in Muscle Wasting: An Overview. *J. Cachexia Sarcopenia Muscle* **2011**, *2*, 143–151. [CrossRef]

79. Lee, S.-J.; Reed, L.A.; Davies, M.V.; Girgenrath, S.; Goad, M.E.P.; Tomkinson, K.N.; Wright, J.F.; Barker, C.; Ehrmantraut, G.; Holmstrom, J.; et al. Regulation of Muscle Growth by Multiple Ligands Signaling through Activin Type II Receptors. *Proc. Natl. Acad. Sci. USA* **2005**, *102*, 18117–18122. [CrossRef]

80. Chen, J.L.; Walton, K.L.; Hagg, A.; Colgan, T.D.; Johnson, K.; Qian, H.; Gregorevic, P.; Harrison, C.A. Specific Targeting of TGF-β Family Ligands Demonstrates Distinct Roles in the Regulation of Muscle Mass in Health and Disease. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 201620013. [CrossRef]

81. Latres, E.; Mastaitis, J.; Fury, W.; Miloscio, L.; Trejos, J.; Pangilinan, J.; Okamoto, H.; Cavino, K.; Na, E.; Papatheodorou, A.; et al. Activin A More Prominently Regulates Muscle Mass in Primates than Does GDF8. *Nat. Commun.* **2017**, *8*, 15153. [CrossRef]

82. Zimmers, T.A. Induction of Cachexia in Mice by Systemically Administered Myostatin. *Science* **2002**, *296*, 1486–1488. [CrossRef]

83. Langley, B.; Thomas, M.; Bishop, A.; Sharma, M.; Gilmour, S.; Kambadur, R. Myostatin Inhibits Myoblast Differentiation by Down-Regulating MyoD Expression. *J. Biol. Chem.* **2002**, *277*, 49831–49840. [CrossRef]

84. Ding, H.; Zhang, G.; Sin, K.W.T.; Liu, Z.; Lin, R.-K.; Li, M.; Li, Y.-P. Activin A Induces Skeletal Muscle Catabolism via P38β Mitogen-Activated Protein Kinase: Activin A Induces Skeletal Muscle Catabolism. *J. Cachexia Sarcopenia Muscle* **2017**, *8*, 202–212. [CrossRef]

85. Sun, G.; Haginoya, K.; Wu, Y.; Chiba, Y.; Nakanishi, T.; Onuma, A.; Sato, Y.; Takigawa, M.; Iinuma, K.; Tsuchiya, S. Connective Tissue Growth Factor Is Overexpressed in Muscles of Human Muscular Dystrophy. *J. Neurol. Sci.* **2008**, *267*, 48–56. [CrossRef] [PubMed]

86. Gonzalez, D.; Rebolledo, D.L.; Correa, L.M.; Court, F.A.; Cerpa, W.; Lipson, K.E.; van Zundert, B.; Brandan, E. The Inhibition of CTGF/CCN2 Activity Improves Muscle and Locomotor Function in a Murine ALS Model. *Hum. Mol. Genet.* **2018**, *27*, 2913–2926. [CrossRef]

87. Morales, M.G.; Gutierrez, J.; Cabello-Verrugio, C.; Cabrera, D.; Lipson, K.E.; Goldschmeding, R.; Brandan, E. Reducing CTGF/CCN2 Slows down Mdx Muscle Dystrophy and Improves Cell Therapy. *Hum. Mol. Genet.* **2013**, *22*, 4938–4951. [CrossRef]

88. Abreu, J.G.; Ketpura, N.I.; Reversade, B.; De Robertis, E.M. Connective-Tissue Growth Factor (CTGF) Modulates Cell Signalling by BMP and TGF-β. *Nat. Cell Biol.* **2002**, *4*, 599–604. [CrossRef]

89. Vial, C.; Zúñiga, L.M.; Cabello-Verrugio, C.; Cañón, P.; Fadic, R.; Brandan, E. Skeletal Muscle Cells Express the Profibrotic Cytokine Connective Tissue Growth Factor (CTGF/CCN2), Which Induces Their Dedifferentiation. *J. Cell. Physiol.* **2008**, *215*, 410–421. [CrossRef] [PubMed]

90. Hillege, M.; Galli Caro, R.; Offringa, C.; de Wit, G.; Jaspers, R.; Hoogaars, W. TGF-β Regulates Collagen Type I Expression in Myoblasts and Myotubes via Transient Ctgf and Fgf-2 Expression. *Cells* **2020**, *9*, 375. [CrossRef] [PubMed]

91. Salazar, V.S.; Gamer, L.W.; Rosen, V. BMP Signalling in Skeletal Development, Disease and Repair. *Nat. Rev. Endocrinol.* **2016**, *12*, 203–221. [CrossRef] [PubMed]

92. Chanon, S.; Chazarin, B.; Toubhans, B.; Durand, C.; Chery, I.; Robert, M.; Vieille-Marchiset, A.; Swenson, J.E.; Zedrosser, A.; Evans, A.L.; et al. Proteolysis Inhibition by Hibernating Bear Serum Leads to Increased Protein Content in Human Muscle Cells. *Sci. Rep.* **2018**, *8*, 5525. [CrossRef] [PubMed]

93. Donahue, S.W.; McGee, M.E.; Harvey, K.B.; Vaughan, M.R.; Robbins, C.T. Hibernating Bears as a Model for Preventing Disuse Osteoporosis. *J. Biomech.* **2006**, *39*, 1480–1488. [CrossRef] [PubMed]
94. Le Goff, C.; Mahaut, C.; Abhyankar, A.; Le Goff, W.; Serre, V.; Afenjar, A.; Destrée, A.; di Rocco, M.; Héron, D.; Jacquemont, S.; et al. Mutations at a Single Codon in Mad Homology 2 Domain of SMAD4 Cause Myhre Syndrome. *Nat. Genet.* **2012**, *44*, 85–88.
[CrossRef]

95. Paul, P.K.; Bhatnagar, S.; Mishra, V.; Srivastava, S.; Darnay, B.G.; Choi, Y.; Kumar, A. The E3 Ubiquitin Ligase TRAF6 Intercedes in Starvation-Induced Skeletal Muscle Atrophy through Multiple Mechanisms. *Mol. Cell. Biol.* **2012**, *32*, 1248–1259.
[CrossRef]
[PubMed]

96. Xiao, F.; Wang, H.; Fu, X.; Li, Y.; Wu, Z. TRAF6 Promotes Myogenic Differentiation via the TAK1/P38 Mitogen-Activated Protein Kinase and Akt Pathways. *PLoS ONE* **2012**, *7*, e34081.
[CrossRef]

97. Taylor, M.V.; Hughes, S.M. Mef2 and the Skeletal Muscle Differentiation Program. *Semin. Cell Dev. Biol.* **2017**, *72*, 33–44.
[CrossRef]
[PubMed]

98. Rullman, E.; Fernandez-Gonzalo, R.; Mekjavić, I.B.; Gustafsson, T.; Eiken, O. MEF2 as Upstream Regulator of the Transcriptome Signature in Human Skeletal Muscle during Unloading. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* **2018**, *315*, R799–R809.
[CrossRef]
[PubMed]

99. Abrigo, J.; Simon, F.; Cabrera, D.; Cordova, G.; Trollet, C.; Cabello-Verrugio, C. Central Role of Transforming Growth Factor Type Beta 1 in Skeletal Muscle Dysfunctions: An Update on Therapeutic Strategies. *Curr. Protein Pept. Sci.* **2018**, *19*, 1189–1200.
[CrossRef]