Transcriptional profiling reveals extraordinary diversity among skeletal muscle tissues

Gene expression differences associated with intrinsic hindfoot muscle loss in the jerboa, Jaculus jaculus

Vertebrate animals that run or jump across sparsely vegetated habitats, such as horses and jerboas, have reduced the number of distal limb bones, and many have lost most or all distal limb muscle. We previously showed that nascent muscles are present in the jerboa hindfoot at birth and that these myofibers are rapidly and completely lost soon after by a process that shares features with pathological skeletal muscle atrophy. Here, we apply an intra- and interspecies differential RNA-Seq approach, comparing jerboa and mouse muscles, to identify gene expression differences associated with the initiation and progression of jerboa hindfoot muscle loss. We show evidence for reduced hepatocyte growth factor and fibroblast growth factor signaling and an imbalance in nitric oxide signaling; all are pathways that are necessary for skeletal muscle development and regeneration. We also find evidence for phagosome formation, which hints at how myofibers may be removed by autophagy or by nonprofessional phagocytes without evidence for cell death or immune cell activation. Last, we show significant overlap between genes associated with jerboa hindfoot muscle loss and genes that are differentially expressed in a variety of human muscle pathologies and rodent models of muscle loss disorders. All together, these data provide molecular insight into the process of evolutionary and developmental muscle loss in jerboa hindfeet.

Maternal Hypermethylated Genes Contribute to Intrauterine Growth Retardation of Piglets in Rongchang Pigs

The placenta is a crucial determinant of fetal survival, growth, and development. Deficiency in placental development directly causes intrauterine growth retardation (IUGR). IUGR can lead to fetal growth restriction and an increase in the mortality rate. The genetic mechanisms underlying IUGR development, however, remain unclear. In the present study, we integrated whole-genome DNA methylation and transcriptomic analyses to determine distinct gene expression patterns in various placental tissues to identify pivotal genes that are implicated with IUGR development. By performing RNA-sequencing analysis, 1487 differentially expressed genes (DEGs), with 737 upregulated and 750 downregulated genes, were identified in IUGR pigs (H_IUGR) compared with that in normal birth weight pigs (N_IUGR) (p < 0.05); furthermore, 77 miRNAs, 1331 lncRNAs, and 61 circRNAs were differentially expressed. The protein-protein interaction network analysis revealed that among these DEGs, the genes GNGT1, ANXA1, and CDC20 related to cellular developmental processes and blood vessel development were the key genes associated with the development of IUGR. A total of 495,870 differentially methylated regions were identified between the N_IUGR and H_IUGR groups, which included 25,053 differentially methylated genes (DMEs); moreover, the overall methylation level was higher in the H_IUGR group than in the N_IUGR group. Combined analysis showed an inverse correlation between methylation levels and gene expression. A total of 1375 genes involved in developmental processes, tissue development, and immune system regulation exhibited methylation differences in gene expression levels in the promoter regions and gene ontology regions. Five genes, namely, ANXA1, ADM, NRP2, SHH, and SMAD1, with high methylation levels were identified as potential contributors to IUGR development. These findings provide valuable insights that DNA methylation plays a crucial role in the epigenetic regulation of gene expression and mammalian development and that DNA-hypermethylated genes contribute to IUGR development in Rongchang pigs.

Temporal regulation of the Mediator complex during muscle proliferation, differentiation, regeneration, aging, and disease

Genesis of skeletal muscle relies on the differentiation and fusion of mono-nucleated muscle progenitor cells into the multi-nucleated muscle fiber syncytium. The temporally-controlled cellular and morphogenetic changes underlying this process are initiated by a series of highly coordinated transcription programs. At the core, the myogenic differentiation cascade is driven by muscle-specific transcription factors, i.e., the Myogenic Regulatory Factors (MRFs). Despite extensive knowledge on the function of individual MRFs, very little is known about how they are coordinated. Ultimately, highly specific coordination of these transcription programs is critical for their masterfully timed transitions, which in turn facilitates the intricate generation of skeletal muscle fibers from a naïve pool of progenitor cells. The Mediator complex links basal transcriptional machinery and transcription factors to regulate transcription and could be the integral component that coordinates transcription factor function during muscle differentiation, growth, and maturation. In this study, we systematically deciphered the changes in Mediator complex subunit expression in skeletal muscle development, regeneration, aging, and disease. We incorporated our in vitro and in vivo experimental results with analysis of publicly available RNA-seq and single nuclei RNA-seq datasets and uncovered the regulation of Mediator subunits in different physiological and temporal contexts. Our experimental results revealed that Mediator subunit expression during myogenesis is highly dynamic. We also discovered unique temporal patterns of Mediator expression in muscle stem cells after injury and during the early regeneration period, suggesting that Mediator subunits may have unique contributions to directing muscle stem cell fate. Although we observed few changes in Mediator subunit expression in aging muscles compared to younger muscles, we uncovered extensive heterogeneity of Mediator subunit expression in dystrophic muscle nuclei, characteristic of chronic muscle degeneration and regeneration cycles. Taken together, our study provides a glimpse of the complex regulation of Mediator subunit expression in the skeletal muscle cell lineage and serves as a springboard for mechanistic studies into the function of individual Mediator subunits in skeletal muscle.

Emerging Therapeutic Strategies in Sarcopenia: An Updated Review on Pathogenesis and Treatment Advances

Sarcopenia is a prevalent degenerative skeletal muscle condition in the elderly population, posing a tremendous burden on diseased individuals and healthcare systems worldwide. Conventionally, sarcopenia is currently managed through nutritional interventions, physical therapy, and lifestyle modification, with no pharmaceutical agents being approved for specific use in this disease. As the pathogenesis of sarcopenia is still poorly understood and there is no treatment recognized as universally effective, recent research efforts have been directed at better comprehending this illness and diversifying treatment strategies. In this respect, this paper overviews the new advances in sarcopenia treatment in correlation with its underlying mechanisms. Specifically, this review creates an updated framework for sarcopenia, describing its etiology, pathogenesis, risk factors, and conventional treatments, further discussing emerging therapeutic approaches like new drug formulations, drug delivery systems, stem cell therapies, and tissue-engineered scaffolds in more detail.

Mustn1 is a smooth muscle cell-secreted microprotein that modulates skeletal muscle extracellular matrix composition

OBJECTIVE

Skeletal muscle plasticity and remodeling are critical for adapting tissue function to use, disuse, and regeneration. The aim of this study was to identify genes and molecular pathways that regulate the transition from atrophy to compensatory hypertrophy or recovery from injury. Here, we have used a mouse model of hindlimb unloading and reloading, which causes skeletal muscle atrophy, and compensatory regeneration and hypertrophy, respectively.

METHODS

We analyzed mouse skeletal muscle at the transition from hindlimb unloading to reloading for changes in transcriptome and extracellular fluid proteome. We then used qRT-PCR, immunohistochemistry, and bulk and single-cell RNA sequencing data to determine Mustn1 gene and protein expression, including changes in gene expression in mouse and human skeletal muscle with different challenges such as exercise and muscle injury. We generated Mustn1-deficient genetic mouse models and characterized them in vivo and ex vivo with regard to muscle function and whole-body metabolism. We isolated smooth muscle cells and functionally characterized them, and performed transcriptomics and proteomics analysis of skeletal muscle and aorta of Mustn1-deficient mice.

RESULTS

We show that Mustn1 (Musculoskeletal embryonic nuclear protein 1, also known as Mustang) is highly expressed in skeletal muscle during the early stages of hindlimb reloading. Mustn1 expression is transiently elevated in mouse and human skeletal muscle in response to intense exercise, resistance exercise, or injury. We find that Mustn1 expression is highest in smooth muscle-rich tissues, followed by skeletal muscle fibers. Muscle from heterozygous Mustn1-deficient mice exhibit differences in gene expression related to extracellular matrix and cell adhesion, compared to wild-type littermates. Mustn1-deficient mice have normal muscle and aorta function and whole-body glucose metabolism. We show that Mustn1 is secreted from smooth muscle cells, and that it is present in arterioles of the muscle microvasculature and in muscle extracellular fluid, particularly during the hindlimb reloading phase. Proteomics analysis of muscle from Mustn1-deficient mice confirms differences in extracellular matrix composition, and female mice display higher collagen content after chemically induced muscle injury compared to wild-type littermates.

CONCLUSIONS

We show that, in addition to its previously reported intracellular localization, Mustn1 is a microprotein secreted from smooth muscle cells into the muscle extracellular space. We explore its role in muscle ECM deposition and remodeling in homeostasis and upon muscle injury. The role of Mustn1 in fibrosis and immune infiltration upon muscle injury and dystrophies remains to be investigated, as does its potential for therapeutic interventions.

Extraocular muscle stem cells exhibit distinct cellular properties associated with non-muscle molecular signatures

Skeletal muscle stem cells (MuSCs) are recognised as functionally heterogeneous. Cranial MuSCs are reported to have greater proliferative and regenerative capacity when compared with those in the limb. A comprehensive understanding of the mechanisms underlying this functional heterogeneity is lacking. Here, we have used clonal analysis, live imaging and single cell transcriptomic analysis to identify crucial features that distinguish extraocular muscle (EOM) from limb muscle stem cell populations. A MyogeninntdTom reporter showed that the increased proliferation capacity of EOM MuSCs correlates with deferred differentiation and lower expression of the myogenic commitment gene Myod. Unexpectedly, EOM MuSCs activated in vitro expressed a large array of extracellular matrix components typical of mesenchymal non-muscle cells. Computational analysis underscored a distinct co-regulatory module, which is absent in limb MuSCs, as driver of these features. The EOM transcription factor network, with Foxc1 as key player, appears to be hardwired to EOM identity as it persists during growth, disease and in vitro after several passages. Our findings shed light on how high-performing MuSCs regulate myogenic commitment by remodelling their local environment and adopting properties not generally associated with myogenic cells.

Determining the feasibility of characterising cellular senescence in human skeletal muscle and exploring associations with muscle morphology and physical function at different ages: findings from the MASS_Lifecourse Study

Cellular senescence may be associated with morphological changes in skeletal muscle and changes in physical function with age although there have been few human studies. We aimed to determine the feasibility of characterising cellular senescence in skeletal muscle and explored sex-specific associations between markers of cellular senescence, muscle morphology, and physical function in participants from the MASS_Lifecourse Study. Senescence markers (p16, TAF (Telomere-Associated DNA Damage Foci), HMGB1 (High Mobility Group Box 1), and Lamin B1) and morphological characteristics (fibre size, number, fibrosis, and centrally nucleated fibres) were assessed in muscle biopsies from 40 men and women (age range 47-84) using spatially-resolved methods (immunohistochemistry, immunofluorescence, and RNA and fluorescence in situ hybridisation). The associations between senescence, morphology, and physical function (muscle strength, mass, and physical performance) at different ages were explored. We found that most senescence markers and morphological characteristics were weakly associated with age in men but more strongly, although non-significantly, associated with age in women. Associations between senescence markers, morphology, and physical function were also stronger in women for HMGB1 and grip strength (r = 0.52); TAF, BMI, and muscle mass (r > 0.4); Lamin B1 and fibrosis (r =  - 0.5); fibre size and muscle mass (r ≥ 0.4); and gait speed (r =  - 0.5). However, these associations were non-significant. In conclusion, we have demonstrated that it is feasible to characterise cellular senescence in human skeletal muscle and to explore associations with morphology and physical function in women and men of different ages. The findings require replication in larger studies.

Mechanosensory encoding in ex vivo muscle-nerve preparations

Our objective was to evaluate an ex vivo muscle-nerve preparation used to study mechanosensory signalling by low threshold mechanosensory receptors (LTMRs). Specifically, we aimed to assess how well the ex vivo preparation represents in vivo firing behaviours of the three major LTMR subtypes of muscle primary sensory afferents, namely type Ia and II muscle spindle (MS) afferents and type Ib tendon organ afferents. Using published procedures for ex vivo study of LTMRs in mouse hindlimb muscles, we replicated earlier reports on afferent firing in response to conventional stretch paradigms applied to non-contracting, that is passive, muscle. Relative to in vivo studies, stretch-evoked firing for confirmed MS afferents in the ex vivo preparation was markedly reduced in firing rate and deficient in encoding dynamic features of muscle stretch. These deficiencies precluded conventional means of discriminating type Ia and II afferents. Muscle afferents, including confirmed Ib afferents were often indistinguishable based on their similar firing responses to the same physiologically relevant stretch paradigms. These observations raise uncertainty about conclusions drawn from earlier ex vivo studies that either attribute findings to specific afferent types or suggest an absence of treatment effects on dynamic firing. However, we found that replacing the recording solution with bicarbonate buffer resulted in afferent firing rates and profiles more like those seen in vivo. Improving representation of the distinctive sensory encoding properties in ex vivo muscle-nerve preparations will promote accuracy in assigning molecular markers and mechanisms to heterogeneous types of muscle mechanosensory neurons.

Tree-based differential testing using inferential uncertainty for RNA-Seq

Identifying differentially expressed transcripts poses a crucial yet challenging problem in transcriptomics. Substantial uncertainty is associated with the abundance estimates of certain transcripts which, if ignored, can lead to the exaggeration of false positives and, if included, may lead to reduced power. For a given set of RNA-Seq samples, TreeTerminus arranges transcripts in a hierarchical tree structure that encodes different layers of resolution for interpretation of the abundance of transcriptional groups, with uncertainty generally decreasing as one ascends the tree from the leaves. We introduce trenDi, which utilizes the tree structure from TreeTerminus for differential testing. The candidate nodes are determined in a data-driven manner to maximize the signal that can be extracted from the data while controlling for the uncertainty associated with estimating the transcript abundances. The identified candidate nodes can include transcripts and inner nodes, with no two nodes having an ancestor/descendant relationship. We evaluated our method on both simulated and experimental datasets, comparing its performance with other tree-based differential methods as well as with uncertainty-aware differential transcript/gene expression methods. Our method detects inner nodes that show a strong signal for differential expression, which would have been overlooked when analyzing the transcripts alone.

Fusion-negative rhabdomyosarcoma 3D organoids to predict effective drug combinations: A proof-of-concept on cell death inducers

Rhabdomyosarcoma (RMS) is the main form of pediatric soft-tissue sarcoma. Its cure rate has not notably improved in the last 20 years following relapse, and the lack of reliable preclinical models has hampered the design of new therapies. This is particularly true for highly heterogeneous fusion-negative RMS (FNRMS). Although methods have been proposed to establish FNRMS organoids, their efficiency remains limited to date, both in terms of derivation rate and ability to accurately mimic the original tumor. Here, we present the development of a next-generation 3D organoid model derived from relapsed adult and pediatric FNRMS. This model preserves the molecular features of the patients' tumors and is expandable for several months in 3D, reinforcing its interest to drug combination screening with longitudinal efficacy monitoring. As a proof-of-concept, we demonstrate its preclinical relevance by reevaluating the therapeutic opportunities of targeting apoptosis in FNRMS from a streamlined approach based on transcriptomic data exploitation.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

High-throughput transcriptome analyses from ASPIRO, a phase 1/2/3 study of gene replacement therapy for X-linked myotubular myopathy

X-linked myotubular myopathy (XLMTM) is a severe congenital disease characterized by profound muscle weakness, respiratory failure, and early death. No approved therapy for XLMTM is currently available. Adeno-associated virus (AAV)-mediated gene replacement therapy has shown promise as an investigational therapeutic strategy. We aimed to characterize the transcriptomic changes in muscle biopsies of individuals with XLMTM who received resamirigene bilparvovec (AT132; rAAV8-Des-hMTM1) in the ASPIRO clinical trial and to identify potential biomarkers that correlate with therapeutic outcome. We leveraged RNA-sequencing data from the muscle biopsies of 15 study participants and applied differential expression analysis, gene co-expression analysis, and machine learning to characterize the transcriptomic changes at baseline (pre-dose) and at 24 and 48 weeks after resamirigene bilparvovec dosing. As expected, MTM1 expression levels were significantly increased after dosing (p < 0.0001). Differential expression analysis identified upregulated genes after dosing that were enriched in several pathways, including lipid metabolism and inflammatory response pathways, and downregulated genes were enriched in cell-cell adhesion and muscle development pathways. Genes involved in inflammatory and immune pathways were differentially expressed between participants exhibiting ventilator support reduction of either greater or less than 6 h/day after gene therapy compared to pre-dosing. Co-expression analysis identified similarly regulated genes, which were grouped into modules. Finally, the machine learning model identified five genes, including MTM1, as potential RNA biomarkers to monitor the progress of AAV gene replacement therapy. These findings further extend our understanding of AAV-mediated gene therapy in individuals with XLMTM at the transcriptomic level.

Suppressed Akt/GSK-3β/β-catenin signaling contributes to excessive adipogenesis of fibro-adipogenic progenitors after rotator cuff tears

Muscular fatty infiltration is a common and troublesome pathology after rotator cuff tears (RCT), which mainly derives from fibro-adipogenic progenitors (FAPs). Compared to the RCT, fatty infiltration is not so severe in Achilles tendon tears (ATT). The knowledge of why fatty infiltration is more likely to occur after RCT is limited. In this study, more severe fatty infiltration was verified in supraspinatus than gastrocnemius muscles after tendon injury. Additionally, we revealed higher adipogenic differentiation ability of RCT-FAPs in vitro. Activation of Akt significantly stimulated GSK-3β/β-catenin signaling and thus decreased PPARγ expression and adipogenesis of RCT-FAPs, while the inhibition effect was attenuated by β-catenin inhibitor. Furthermore, Wnt signaling activator BML-284 limited adipogenesis of RCT-FAPs, alleviated muscular fatty infiltration, and improved parameters in gait analysis and treadmill test for RCT model. In conclusion, our study demonstrated that suppressed Akt/GSK-3β/β-catenin signaling increased PPARγ expression and thus contributed to excessive adipogenesis in RCT-FAPs. Modulation of Akt/GSK-3β/β-catenin signaling ameliorated excessive fatty infiltration of rotator cuff muscles and improved shoulder function after RCT.

Transcriptomic analysis of paired healthy human skeletal muscles to identify modulators of disease severity in DMD

Muscle damage and fibro-fatty replacement of skeletal muscles is a main pathologic feature of Duchenne muscular dystrophy (DMD) with more proximal muscles affected earlier and more distal affected later in the disease course, suggesting that different skeletal muscle groups possess distinctive characteristics that influence their susceptibility to disease. To explore transcriptomic factors driving differential gene expression and modulating DMD skeletal muscle severity, we characterized the transcriptome of vastus lateralis (VL), a more proximal and susceptible muscle, relative to tibialis anterior (TA), a more distal and protected muscle, in 15 healthy individuals using bulk RNA sequencing to identify gene expression differences that may mediate their relative susceptibility to damage with loss of dystrophin. Matching single nuclei RNA sequencing data was generated for 3 of the healthy individuals, to infer cell composition in the bulk RNA sequencing dataset and to improve mapping of differentially expressed genes to their cell source of expression. A total of 3,410 differentially expressed genes were identified and mapped to cell type using single nuclei RNA sequencing of muscle, including long non-coding RNAs and protein coding genes. There was an enrichment of genes involved in calcium release from the sarcoplasmic reticulum, particularly in the myofibers and these myofiber genes were higher in the VL. There was an enrichment of genes in "Collagen-Containing Extracellular Matrix" expressed by fibroblasts, endothelial, smooth muscle and pericytes, with most genes higher in the TA, as well as genes in "Regulation Of Apoptotic Process" expressed across all cell types. Previously reported genetic modifiers were also enriched within the differentially expressed genes. We also identify 6 genes with differential isoform usage between the VL and TA. Lastly, we integrate our findings with DMD RNA sequencing data from the TA, and identify "Collagen-Containing Extracellular Matrix" and "Negative Regulation Of Apoptotic Process" as differentially expressed between DMD compared to healthy. Collectively, these findings propose novel candidate mechanisms that may mediate differential muscle susceptibility in muscular dystrophies and provide new insight into potential therapeutic targets.

When Size Really Matters: The Eccentricities of Dystrophin Transcription and the Hazards of Quantifying mRNA from Very Long Genes

At 2.3 megabases in length, the dystrophin gene is enormous: transcription of a single mRNA requires approximately 16 h. Principally expressed in skeletal muscle, the dystrophin protein product protects the muscle sarcolemma against contraction-induced injury, and dystrophin deficiency results in the fatal muscle-wasting disease, Duchenne muscular dystrophy. This gene is thus of key clinical interest, and therapeutic strategies aimed at eliciting dystrophin restoration require quantitative analysis of its expression. Approaches for quantifying dystrophin at the protein level are well-established, however study at the mRNA level warrants closer scrutiny: measured expression values differ in a sequence-dependent fashion, with significant consequences for data interpretation. In this manuscript, we discuss these nuances of expression and present evidence to support a transcriptional model whereby the long transcription time is coupled to a short mature mRNA half-life, with dystrophin transcripts being predominantly nascent as a consequence. We explore the effects of such a model on cellular transcriptional dynamics and then discuss key implications for the study of dystrophin gene expression, focusing on both conventional (qPCR) and next-gen (RNAseq) approaches.

TreeTerminus -creating transcript trees using inferential replicate counts

A certain degree of uncertainty is always associated with the transcript abundance estimates. The uncertainty may make many downstream analyses, such as differential testing, difficult for certain transcripts. Conversely, gene-level analysis, though less ambiguous, is often too coarse-grained. We introduce TreeTerminus, a data-driven approach for grouping transcripts into a tree structure where leaves represent individual transcripts and internal nodes represent an aggregation of a transcript set. TreeTerminus constructs trees such that, on average, the inferential uncertainty decreases as we ascend the tree topology. The tree provides the flexibility to analyze data at nodes that are at different levels of resolution in the tree and can be tuned depending on the analysis of interest. We evaluated TreeTerminus on two simulated and two experimental datasets and observed an improved performance compared to transcripts (leaves) and other methods under several different metrics.

RhoA/ROCK signalling activated by ARHGEF3 promotes muscle weakness via autophagy in dystrophic mdx mice

BACKGROUND

Duchenne muscular dystrophy (DMD), caused by dystrophin deficiency, leads to progressive and fatal muscle weakness through yet-to-be-fully deciphered molecular perturbations. Emerging evidence implicates RhoA/Rho-associated protein kinase (ROCK) signalling in DMD pathology, yet its direct role in DMD muscle function, and related mechanisms, are unknown.

METHODS

Three-dimensionally engineered dystrophin-deficient mdx skeletal muscles and mdx mice were used to test the role of ROCK in DMD muscle function in vitro and in situ, respectively. The role of ARHGEF3, one of the RhoA guanine nucleotide exchange factors (GEFs), in RhoA/ROCK signalling and DMD pathology was examined by generating Arhgef3 knockout mdx mice. The role of RhoA/ROCK signalling in mediating the function of ARHGEF3 was determined by evaluating the effects of wild-type or GEF-inactive ARHGEF3 overexpression with ROCK inhibitor treatment. To gain more mechanistic insights, autophagy flux and the role of autophagy were assessed in various conditions with chloroquine.

RESULTS

Inhibition of ROCK with Y-27632 improved muscle force production in 3D-engineered mdx muscles (+25% from three independent experiments, P < 0.05) and in mice (+25%, P < 0.001). Unlike suggested by previous studies, this improvement was independent of muscle differentiation or quantity and instead related to increased muscle quality. We found that ARHGEF3 was elevated and responsible for RhoA/ROCK activation in mdx muscles, and that depleting ARHGEF3 in mdx mice restored muscle quality (up to +36%, P < 0.01) and morphology without affecting regeneration. Conversely, overexpressing ARHGEF3 further compromised mdx muscle quality (-13% vs. empty vector control, P < 0.01) in GEF activity- and ROCK-dependent manner. Notably, ARHGEF3/ROCK inhibition exerted the effects by rescuing autophagy which is commonly impaired in dystrophic muscles.

CONCLUSIONS

Our findings uncover a new pathological mechanism of muscle weakness in DMD involving the ARHGEF3-ROCK-autophagy pathway and the therapeutic potential of targeting ARHGEF3 in DMD.

A transcriptome atlas of leg muscles from healthy human volunteers reveals molecular and cellular signatures associated with muscle location

Skeletal muscles support the stability and mobility of the skeleton but differ in biomechanical properties and physiological functions. The intrinsic factors that regulate muscle-specific characteristics are poorly understood. To study these, we constructed a large atlas of RNA-seq profiles from six leg muscles and two locations from one muscle, using biopsies from 20 healthy young males. We identified differential expression patterns and cellular composition across the seven tissues using three bioinformatics approaches confirmed by large-scale newly developed quantitative immune-histology procedures. With all three procedures, the muscle samples clustered into three groups congruent with their anatomical location. Concomitant with genes marking oxidative metabolism, genes marking fast- or slow-twitch myofibers differed between the three groups. The groups of muscles with higher expression of slow-twitch genes were enriched in endothelial cells and showed higher capillary content. In addition, expression profiles of Homeobox (HOX) transcription factors differed between the three groups and were confirmed by spatial RNA hybridization. We created an open-source graphical interface to explore and visualize the leg muscle atlas (https://tabbassidaloii.shinyapps.io/muscleAtlasShinyApp/). Our study reveals the molecular specialization of human leg muscles, and provides a novel resource to study muscle-specific molecular features, which could be linked with (patho)physiological processes.

Distinctive chaperonopathy in skeletal muscle associated with the dominant variant in DNAJB4

DnaJ homolog, subfamily B, member 4, a member of the heat shock protein 40 chaperones encoded by DNAJB4, is highly expressed in myofibers. We identified a heterozygous c.270 T > A (p.F90L) variant in DNAJB4 in a family with a dominantly inherited distal myopathy, in which affected members have specific features on muscle pathology represented by the presence of cytoplasmic inclusions and the accumulation of desmin, p62, HSP70, and DNAJB4 predominantly in type 1 fibers. Both Dnajb4F90L knockin and knockout mice developed muscle weakness and recapitulated the patient muscle pathology in the soleus muscle, where DNAJB4 has the highest expression. These data indicate that the identified variant is causative, resulting in defective chaperone function and selective muscle degeneration in specific muscle fibers. This study demonstrates the importance of DNAJB4 in skeletal muscle proteostasis by identifying the associated chaperonopathy.

Defining the age-dependent and tissue-specific circadian transcriptome in male mice

Cellular circadian clocks direct a daily transcriptional program that supports homeostasis and resilience. Emerging evidence has demonstrated age-associated changes in circadian functions. To define age-dependent changes at the systems level, we profile the circadian transcriptome in the hypothalamus, lung, heart, kidney, skeletal muscle, and adrenal gland in three age groups. We find age-dependent and tissue-specific clock output changes. Aging reduces the number of rhythmically expressed genes (REGs), indicative of weakened circadian control. REGs are enriched for the hallmarks of aging, adding another dimension to our understanding of aging. Analyzing differential gene expression within a tissue at four different times of day identifies distinct clusters of differentially expressed genes (DEGs). Increased variability of gene expression across the day is a common feature of aged tissues. This analysis extends the landscape for understanding aging and highlights the impact of aging on circadian clock function and temporal changes in gene expression.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Exogenous Oncostatin M induces Cardiac Dysfunction, Musculoskeletal Atrophy, and Fibrosis

Musculoskeletal diseases such as muscular dystrophy, cachexia, osteoarthritis, and rheumatoid arthritis impair overall physical health and reduce survival. Patients suffer from pain, dysfunction, and dysmobility due to inflammation and fibrosis in bones, muscles, and joints, both locally and systemically. The Interleukin-6 (IL-6) family of cytokines, most notably IL-6, is implicated in musculoskeletal disorders and cachexia. Here we show elevated circulating levels of OSM in murine pancreatic cancer cachexia and evaluate the effects of the IL-6 family member, Oncostatin M (OSM), on muscle and bone using adeno-associated virus (AAV) mediated over-expression of murine OSM in wildtype and IL-6 deficient mice. Initial studies with high titer AAV-OSM injection yielded high circulating OSM and IL-6, thrombocytosis, inflammation, and 60% mortality without muscle loss within 4 days. Subsequently, to mimic OSM levels in cachexia, a lower titer of AAV-OSM was used in wildtype and Il6 null mice, observing effects out to 4 weeks and 12 weeks. AAV-OSM caused muscle atrophy and fibrosis in the gastrocnemius, tibialis anterior, and quadriceps of the injected limb, but these effects were not observed on the non-injected side. In contrast, OSM induced both local and distant trabecular bone loss as shown by reduced bone volume, trabecular number, and thickness, and increased trabecular separation. OSM caused cardiac dysfunction including reduced ejection fraction and reduced fractional shortening. RNA-sequencing of cardiac muscle revealed upregulation of genes related to inflammation and fibrosis. None of these effects were different in IL-6 knockout mice. Thus, OSM induces local muscle atrophy, systemic bone loss, tissue fibrosis, and cardiac dysfunction independently of IL-6, suggesting a role for OSM in musculoskeletal conditions with these characteristics, including cancer cachexia.

Phosphoproteomics of three exercise modalities identifies canonical signaling and C18ORF25 as an AMPK substrate regulating skeletal muscle function

Exercise induces signaling networks to improve muscle function and confer health benefits. To identify divergent and common signaling networks during and after different exercise modalities, we performed a phosphoproteomic analysis of human skeletal muscle from a cross-over intervention of endurance, sprint, and resistance exercise. This identified 5,486 phosphosites regulated during or after at least one type of exercise modality and only 420 core phosphosites common to all exercise. One of these core phosphosites was S67 on the uncharacterized protein C18ORF25, which we validated as an AMPK substrate. Mice lacking C18ORF25 have reduced skeletal muscle fiber size, exercise capacity, and muscle contractile function, and this was associated with reduced phosphorylation of contractile and Ca2+ handling proteins. Expression of C18ORF25 S66/67D phospho-mimetic reversed the decreased muscle force production. This work defines the divergent and canonical exercise phosphoproteome across different modalities and identifies C18ORF25 as a regulator of exercise signaling and muscle function.

Gene expression profiles of specific chicken skeletal muscles

The chicken provides large amounts of protein for the human diet and is also used as a model organism for biomedical research. Increasing meat production is an important goal in the poultry industry and skeletal muscles have highly diverse origins, shapes, metabolic features, and physical functions. Previous gene expression atlases have largely ignored the differences among diverse types of skeletal muscles; therefore, comprehensive transcriptional maps of all skeletal muscles are needed to improve meat production traits. In this study, we sequenced 58 samples from 10 different skeletal muscles of 42-day-old White Plymouth Rock chickens. We also measured myofiber diameter and generated myofiber-type datasets of these 10 tissues. We generated 418.4 Gb high-quality bulk RNA-Seq data from four or six biological replicates of each skeletal muscle (four replicates from extraocular samples) (approximately 7.4 Gb per sample). This dataset provides valuable information for understanding the muscle fiber characteristics of White Plymouth Rock chickens. Furthermore, our data can be used as a model for heterogeneity analysis between tissues with similar properties.

Measurement(s)	Gene expression profiles of specific chicken skeletal muscles • cross section area of chicken skeletal muscles and type I muscle fiber quantity
Technology Type(s)	RNA sequence • Histological Procedure
Factor Type(s)	different skeletal muscles
Sample Characteristic - Organism	Gallus gallus
Sample Characteristic - Environment	farm

Systematic identification of aberrant non-coding RNAs and their mediated modules in rotator cuff tears

Background: Rotator cuff tears (RCT) is the most common cause of shoulder dysfunction, however, its molecular mechanisms remain unclear. Non-coding RNAs(ncRNAs), such as long ncRNA (lncRNA), microRNA (miRNA) and circular RNA (circRNA), are involved in a variety of diseases, but little is known about their roles in RCT. Therefore, the purpose of this study is to identify dysregulated ncRNAs and understand how they influence RCT.

Methods: We performed RNA sequencing and miRNA sequencing on five pairs of torn supraspinatus muscles and matched unharmed subscapularis muscles to identify RNAs dysregulated in RCT patients. To better comprehend the fundamental biological processes, we carried out enrichment analysis of these dysregulated mRNAs or the co-expressed genes of dysregulated ncRNAs. According to the competing endogenous RNA (ceRNA) theory, we finally established ceRNA networks to explore the relationship among dysregulated RNAs in RCT.

Results: A total of 151 mRNAs, 38 miRNAs, 20 lncRNAs and 90 circRNAs were differentially expressed between torn supraspinatus muscles and matched unharmed subscapularis muscles, respectively. We found that these dysregulated mRNAs, the target mRNAs of these dysregulated miRNAs or the co-expressed mRNAs of these dysregulated ncRNAs were enriched in muscle structure development, actin-mediated cell contraction and actin binding. Then we constructed and analyzed the ceRNA network and found that the largest module in the ceRNA network was associated with vasculature development. Based on the topological properties of the largest module, we identified several important ncRNAs including hsa_circ_0000722, hsa-miR-129-5p and hsa-miR-30c-5p, whose interacting mRNAs related to muscle diseases, fat and inflammation.

Conclusion: This study presented a systematic dissection of the expression profile of mRNAs and ncRNAs in RCT patients and revealed some important ncRNAs which may contribute to the development of RCT. Such results could provide new insights for further research on RCT.

MeDIP-seq and RNA-seq analysis during porcine testis development reveals functional DMRs at the promoter of LDHC

Testis development requires tight regulation of gene expression programmed by epigenetic modifiers. However, their mechanism remains to be elucidated. Here, we investigated the genome-wide DNA methylation landscape in the Duroc and Meishan boar testes using methylated DNA immunoprecipitation sequencing (MeDIP-seq). We identified over 1100 promoter differential methylation genes (DMGs) before and after puberty, most of which are associated with testis development. Furthermore, we discovered that the expression of lactate dehydrogenase C (LDHC) gene during testis development is regulated by DNA methylation. The promoter of LDHC in pre-pubertal testes is substantially methylated, whereas considerably demethylated in post-pubertal testes. Artificial demethylation with the demethylating agent 5-Aza-CdR induced LDHC expression in immature Sertoli cells (SCs). Mechanistically, we confirmed the transcription factor SP1 was recruited to bind in hypomethylated differentially methylated regions (DMRs) in LDHC promoter, which upregulated the expression of LDHC. Functionally, we demonstrated that LDHC gene was activated in mature SCs (mSCs) and its overexpression significantly increases lactate secretion in SCs. In conclusion, our results highlight the function and regulation of dynamic DNA methylation in testis development.

A potent physiological method to magnify and sustain soleus oxidative metabolism improves glucose and lipid regulation

Slow oxidative muscle, most notably the soleus, is inherently well equipped with the molecular machinery for regulating blood-borne substrates. However, the entire human musculature accounts for only ∼15% of the body’s oxidative metabolism of glucose at the resting energy expenditure, despite being the body’s largest lean tissue mass. We found the human soleus muscle could raise local oxidative metabolism to high levels for hours without fatigue, during a type of soleus-dominant activity while sitting, even in unfit volunteers. Muscle biopsies revealed there was minimal glycogen use. Magnifying the otherwise negligible local energy expenditure with isolated contractions improved systemic VLDL-triglyceride and glucose homeostasis by a large magnitude, e.g., 52% less postprandial glucose excursion (∼50 mg/dL less between ∼1 and 2 h) with 60% less hyperinsulinemia. Targeting a small oxidative muscle mass (∼1% body mass) with local contractile activity is a potent method for improving systemic metabolic regulation while prolonging the benefits of oxidative metabolism.

•

We developed a method to capitalize upon the unique phenotype of the soleus

•

“A high quality versus large quantity perspective” for muscle activation

•

Singular movement targeting the 1 kg soleus easily sustains oxidative metabolism

•

This method provides a distinct muscular activity stimulus for metabolic control

Potassium homeostasis: sensors, mediators, and targets

Transmembrane potassium (K) gradients are key determinants of membrane potential that can modulate action potentials, control muscle contractility, and influence ion channel and transporter activity. Daily K intake is normally equal to the amount of K in the entire extracellular fluid (ECF) creating a critical challenge - how to maintain ECF [K] and membrane potential in a narrow range during feast and famine. Adaptations to maintain ECF [K] include sensing the K intake, sensing ECF [K] vs. desired set-point and activating mediators that regulate K distribution between ECF and ICF, and regulate renal K excretion. In this focused review, we discuss the basis of these adaptions, including (1) potential mechanisms for rapid feedforward signaling to kidney and muscle after a meal (before a rise in ECF [K]), (2) how skeletal muscles sense and respond to changes in ECF [K], (3) effects of K on aldosterone biosynthesis, and (4) how the kidney responds to changes in ECF [K] to modify K excretion. The concepts of sexual dimorphisms in renal K handling adaptation are introduced, and the molecular mechanisms that can account for the benefits of a K-rich diet to maintain cardiovascular health are discussed. Although the big picture of K homeostasis is becoming more clear, we also highlight significant pieces of the puzzle that remain to be solved, including knowledge gaps in our understanding of initiating signals, sensors and their connection to homeostatic adjustments of ECF [K].

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

Comprehensive multi-cohort transcriptional meta-analysis of muscle diseases identifies a signature of disease severity

Muscle diseases share common pathological features suggesting common underlying mechanisms. We hypothesized there is a common set of genes dysregulated across muscle diseases compared to healthy muscle and that these genes correlate with severity of muscle disease. We performed meta-analysis of transcriptional profiles of muscle biopsies from human muscle diseases and healthy controls. Studies obtained from public microarray repositories fulfilling quality criteria were divided into six categories: (i) immobility, (ii) inflammatory myopathies, (iii) intensive care unit (ICU) acquired weakness (ICUAW), (iv) congenital muscle diseases, (v) chronic systemic diseases, (vi) motor neuron disease. Patient cohorts were separated in discovery and validation cohorts retaining roughly equal proportions of samples for the disease categories. To remove bias towards a specific muscle disease category we repeated the meta-analysis five times by removing data sets corresponding to one muscle disease class at a time in a "leave-one-disease-out" analysis. We used 636 muscle tissue samples from 30 independent cohorts to identify a 52 gene signature (36 up-regulated and 16 down-regulated genes). We validated the discriminatory power of this signature in 657 muscle biopsies from 12 additional patient cohorts encompassing five categories of muscle diseases with an area under the receiver operating characteristic curve of 0.91, 83% sensitivity, and 85.3% specificity. The expression score of the gene signature inversely correlated with quadriceps muscle mass (r = -0.50, p-value = 0.011) in ICUAW and shoulder abduction strength (r = -0.77, p-value = 0.014) in amyotrophic lateral sclerosis (ALS). The signature also positively correlated with histologic assessment of muscle atrophy in ALS (r = 0.88, p-value = 1.62 × 10^-3) and fibrosis in muscular dystrophy (Jonckheere trend test p-value = 4.45 × 10^-9). Our results identify a conserved transcriptional signature associated with clinical and histologic muscle disease severity. Several genes in this conserved signature have not been previously associated with muscle disease severity.

Relationship between gene expression networks and muscle contractile physiology differences in Anolis lizards

Muscles facilitate most animal behavior, from eating to fleeing. However, to generate the variation in behavior necessary for survival, different muscles must perform differently; for instance, sprinting requires multiple rapid muscle contractions, whereas biting may require fewer contractions but greater force. Here, we use a transcriptomic approach to identify genes associated with variation in muscle contractile physiology among different muscles from the same individual. We measured differential gene expression between a leg and jaw muscle of Anolis lizards known to differ in muscle contractile physiology and performance. For each individual, one muscle was used to measure muscle contractile physiology, including contractile velocity (V_max and V₄₀), specific tension, power ratio, and twitch time, whereas the contralateral muscle was used to extract RNA for transcriptomic sequencing. Using the transcriptomic data, we found clear clustering of muscle type. Expression of genes clustered in gene ontology (GO) terms related to muscle contraction and extracellular matrix was, on average, negatively correlated with V_max and slower twitch times but positively correlated to power ratio and V₄₀. Conversely, genes related to the GO terms related to aerobic respiration were downregulated in muscles with higher power ratio and V₄₀, and over-expressed as twitch time decreased. Determining the molecular mechanisms that underlie variation in muscle contractile physiology can begin to explain how organisms are able to optimize behavior under variable conditions. Future studies pursuing the effects of differential gene expression across muscle types in different environments might inform researchers about how differences develop across species, populations, and individuals varying in ecological history.

Specialized androgen synthesis in skeletal muscles that actuate elaborate social displays

Androgens mediate the expression of many reproductive behaviors, including the elaborate displays used to navigate courtship and territorial interactions. In some vertebrates, males can produce androgen-dependent sexual behavior even when levels of testosterone (T) is low in the bloodstream. One idea is that select tissues make their own androgens from scratch to support behavioral performance. We first study this phenomenon in the skeletal muscles that actuate elaborate sociosexual displays in downy woodpeckers and two songbirds. We show that the woodpecker display muscle maintains elevated T when the testes are regressed in the non-breeding season. Both the display muscles of woodpeckers, as well as the display muscles in the avian vocal organ (syrinx or SYR) of songbirds, express all transporters and enzymes necessary to convert cholesterol into bioactive androgens locally. In a final analysis, we broaden our study by looking for these same transporters and enzymes in mammalian muscles that operate at different speeds. Using RNA-seq data, we find that the capacity for de novo synthesis is only present in "superfast" extraocular muscle. Together, our results suggest that skeletal muscle specialized to generate extraordinary twitch-times and/or extremely rapid contractile speeds may depend on androgenic hormones produced locally within the muscle itself. Our study therefore uncovers an important new dimension of androgenic regulation of behavior.

Unraveling the Molecular Basis of the Dystrophic Process in Limb-Girdle Muscular Dystrophy LGMD-R12 by Differential Gene Expression Profiles in Diseased and Healthy Muscles

Limb-girdle muscular dystrophy R12 (LGMD-R12) is caused by two mutations in anoctamin-5 (ANO5). Our aim was to identify genes and pathways that underlie LGMD-R12 and explain differences in the molecular predisposition and susceptibility between three thigh muscles that are severely (semimembranosus), moderately (vastus lateralis) or mildly (rectus femoris) affected in this disease. We performed transcriptomics on these three muscles in 16 male LGMD-R12 patients and 15 age-matched male controls. Our results showed that LGMD-R12 dystrophic muscle is associated with the expression of genes indicative of fibroblast and adipocyte replacement, such as fibroadipogenic progenitors and immune cell infiltration, while muscle protein synthesis and metabolism were downregulated. Muscle degeneration was associated with an increase in genes involved in muscle injury and inflammation, and muscle repair/regeneration. Baseline differences between muscles in healthy individuals indicated that muscles that are the most affected by LGMD-R12 have the lowest expression of transcription factor networks involved in muscle (re)generation and satellite stem cell activation. Instead, they show relative high levels of fetal/embryonic myosins, all together indicating that muscles differ in their baseline regenerative potential. To conclude, we profiled the gene expression landscape in LGMD-R12, identified baseline differences in expression levels between differently affected muscles and characterized disease-associated changes.

Daily running enhances molecular and physiological circadian rhythms in skeletal muscle

Objective

Exercise is a critical component of a healthy lifestyle and a key strategy for the prevention and management of metabolic disease. Identifying molecular mechanisms underlying adaptation in response to chronic physical activity is of critical interest in metabolic physiology. Circadian rhythms broadly modulate metabolism, including muscle substrate utilization and exercise capacity. Here, we define the molecular and physiological changes induced across the daily cycle by voluntary low intensity daily exercise.

Methods

Wildtype C57BL6/J male and female mice were housed with or without access to a running wheel for six weeks. Maximum running speed was measured at four different zeitgeber times (ZTs, hours after lights on) using either electrical or manual stimulation to motivate continued running on a motorized treadmill. RNA isolated from plantaris muscles at six ZTs was sequenced to establish the impact of daily activity on genome-wide transcription. Patterns of gene expression were analyzed using Gene Set Enrichment Analysis (GSEA) and Detection of Differential Rhythmicity (DODR). Blood glucose, lactate, and ketones, and muscle and liver glycogen were measured before and after exercise.

Results

We demonstrate that the use of mild electrical shocks to motivate running negatively impacts maximum running speed in mice, and describe a manual method to motivate running in rodent exercise studies. Using this method, we show that time of day influences the increase in exercise capacity afforded by six weeks of voluntary wheel running: when maximum running speed is measured at the beginning of the nighttime active period in mice, there is no measurable benefit from a history of daily voluntary running, while maximum increase in performance occurs at the end of the night. We show that daily voluntary exercise dramatically remodels the murine muscle circadian transcriptome. Finally, we describe daily rhythms in carbohydrate metabolism associated with the time-dependent response to moderate daily exercise in mice.

Conclusions

Collectively, these data indicate that chronic nighttime physical activity dramatically remodels daily rhythms of murine muscle gene expression, which in turn support daily fluctuations in exercise performance.

•

Daily voluntary running dramatically remodels the mouse muscle circadian transcriptome.

•

Daily voluntary running maximally increases mouse running speed in the late active period.

•

Muscle and liver glycogen content exhibit robust daily rhythms in laboratory mice.

•

Use of mild electric shocks to motivate running in mice impairs maximum running speed.

Bridging the muscle genome to phenome across multiple biological scales

Muscle is highly hierarchically organized, with functions shaped by genetically controlled expression of protein ensembles with different isoform profiles at the sarcomere scale. However, it remains unclear how isoform profiles shape whole-muscle performance. We compared two mouse hindlimb muscles, the slow, relatively parallel-fibered soleus and the faster, more pennate-fibered tibialis anterior (TA), across scales: from gene regulation, isoform expression and translation speed, to force-length-velocity-power for intact muscles. Expression of myosin heavy-chain (MHC) isoforms directly corresponded with contraction velocity. The fast-twitch TA with fast MHC isoforms had faster unloaded velocities (actin sliding velocity, Vactin; peak fiber velocity, Vmax) than the slow-twitch soleus. For the soleus, Vactin was biased towards Vactin for purely slow MHC I, despite this muscle's even fast and slow MHC isoform composition. Our multi-scale results clearly identified a consistent and significant dampening in fiber shortening velocities for both muscles, underscoring an indirect correlation between Vactin and fiber Vmax that may be influenced by differences in fiber architecture, along with internal loading due to both passive and active effects. These influences correlate with the increased peak force and power in the slightly more pennate TA, leading to a broader length range of near-optimal force production. Conversely, a greater force-velocity curvature in the near-parallel fibered soleus highlights the fine-tuning by molecular-scale influences including myosin heavy and light chain expression along with whole-muscle characteristics. Our results demonstrate that the individual gene, protein and whole-fiber characteristics do not directly reflect overall muscle performance but that intricate fine-tuning across scales shapes specialized muscle function.

Exercise reduces intramuscular stress and counteracts muscle weakness in mice with breast cancer

BACKGROUND

Patients with breast cancer exhibit muscle weakness, which is associated with increased mortality risk and reduced quality of life. Muscle weakness is experienced even in the absence of loss of muscle mass in breast cancer patients, indicating intrinsic muscle dysfunction. Physical activity is correlated with reduced cancer mortality and disease recurrence. However, the molecular processes underlying breast cancer-induced muscle weakness and the beneficial effect of exercise are largely unknown.

METHODS

Eight-week-old breast cancer (MMTV-PyMT, PyMT) and control (WT) mice had access to active or inactive in-cage voluntary running wheels for 4 weeks. Mice were also subjected to a treadmill test. Muscle force was measured ex vivo. Tumour markers were determined with immunohistochemistry. Mitochondrial biogenesis and function were assessed with transcriptional analyses of PGC-1α, the electron transport chain (ETC) and antioxidants superoxide dismutase (Sod) and catalase (Cat), combined with activity measurements of SOD, citrate synthase (CS) and β-hydroxyacyl-CoA-dehydrogenase (βHAD). Serum and intramuscular stress levels were evaluated by enzymatic assays, immunoblotting, and transcriptional analyses of, for example, tumour necrosis factor-α (TNF-α) and p38 mitogen-activated protein kinase (MAPK) signalling.

RESULTS

PyMT mice endured shorter time and distance during the treadmill test (~30%, P < 0.05) and ex vivo force measurements revealed ~25% weaker slow-twitch soleus muscle (P < 0.001). This was independent of cancer-induced alteration of muscle size or fibre type. Inflammatory stressors in serum and muscle, including TNF-α and p38 MAPK, were higher in PyMT than in WT mice (P < 0.05). Cancer-induced decreases in ETC (P < 0.05, P < 0.01) and antioxidant gene expression were observed (P < 0.05). The exercise intervention counteracted the cancer-induced muscle weakness and was accompanied by a less aggressive, differentiated tumour phenotype, determined by increased CK8 and reduced CK14 expression (P < 0.05). In PyMT mice, the exercise intervention led to higher CS activity (P = 0.23), enhanced β-HAD and SOD activities (P < 0.05), and reduced levels of intramuscular stressors together with a normalization of the expression signature of TNFα-targets and ETC genes (P < 0.05, P < 0.01). At the same time, the exercise-induced PGC-1α expression, and CS and β-HAD activity was blunted in muscle from the PyMT mice as compared with WT mice, indicative that breast cancer interfere with transcriptional programming of mitochondria and that the molecular adaptation to exercise differs between healthy mice and those afflicted by disease.

CONCLUSIONS

Four-week voluntary wheel running counteracted muscle weakness in PyMT mice which was accompanied by reduced intrinsic stress and improved mitochondrial and antioxidant profiles and activities that aligned with muscles of healthy mice.

RNA Sequencing of Whole Blood Defines the Signature of High Intensity Exercise at Altitude in Elite Speed Skaters

Although high altitude training has been increasingly popular among endurance athletes, the molecular and cellular bases of this adaptation remain poorly understood. We aimed to define the underlying physiological changes and screen for potential biomarkers of adaptation using transcriptional profiling of whole blood. Seven elite female speed skaters were profiled on the 18th day of high-altitude adaptation. Whole blood RNA-seq before and after an intense 1 h skating bout was used to measure gene expression changes associated with exercise. In order to identify the genes specifically regulated at high altitudes, we have leveraged the data from eight previously published microarray datasets studying blood expression changes after exercise at sea level. Using cell type-specific signatures, we were able to deconvolute changes of cell type abundance from individual gene expression changes. Among these were PHOSPHO1, with a known role in erythropoiesis, and MARC1 with a role in endogenic NO metabolism. We find that platelet and erythrocyte counts uniquely respond to altitude exercise, while changes in neutrophils represent a more generic marker of intense exercise. Publicly available data from both single cell atlases and exercise-related blood profiling dramatically increases the value of whole blood RNA-seq for the dynamic evaluation of physiological changes in an athlete’s body.

Tongue Muscle for the Analysis of Head Muscle Regeneration Dynamics

Tongue muscle damage impairs speaking and eating, thereby degrading overall health and quality of life. Skeletal muscles of the body are diverse in embryonic origin, anatomic location, and gene expression profiles. Responses to disease, atrophy, aging, or drugs vary among different muscles. Currently, most muscle studies are focused on limb muscles and the tongue is neglected. The regenerative ability of tongue muscle remains unknown, and thus there is need for tongue muscle research models. Here, we present a comprehensive characterization of the spatiotemporal dynamics in a mouse model of tongue muscle regeneration and establish a method for the isolation of primary tongue-derived satellite cells. We compare and contrast our observations with the tibialis anterior (TA) limb muscle. Acute injury was induced by intramuscular injection of cardiotoxin, a cytolytic agent, and examined at multiple timepoints. Initially, necrotic myofibers with fragmented sarcoplasm became infiltrated with inflammatory cells. Concomitantly, satellite cells expanded rapidly. Seven days postinjury, regenerated myofibers with centralized nuclei appeared. Full regeneration, as well as an absence of fibrosis, was evident 21 d postinjury. Primary tongue-derived satellite cells were isolated by enzymatic separation of tongue epithelium from mesenchyme followed by magnetic-activated cell sorting. We observed that tongue displays an efficient regenerative response similar to TA but with slightly faster kinetics. In vitro, tongue-derived satellite cells differentiated robustly into mature myotubes with spontaneous contractile behavior and myogenic marker expression. Comparison of gene expression signatures between tongue and TA-derived satellite cells revealed differences in the expression of positional-identity genes, including the HOX family. In conclusion, we have established a model for tongue regeneration useful for investigations of orofacial muscle biology. Furthermore, we showed that tongue is a viable source of satellite cells with unique properties and inherited positional memory.

Lipidomic Analyses Reveal Specific Alterations of Phosphatidylcholine in Dystrophic Mdx Muscle

In Duchenne muscular dystrophy (DMD), lack of dystrophin increases the permeability of myofiber plasma membranes to ions and larger macromolecules, disrupting calcium signaling and leading to progressive muscle wasting. Although the biological origin and meaning are unclear, alterations of phosphatidylcholine (PC) are reported in affected skeletal muscles of patients with DMD that may include higher levels of fatty acid (FA) 18:1 chains and lower levels of FA 18:2 chains, possibly reflected in relatively high levels of PC 34:1 (with 16:0_18:1 chain sets) and low levels of PC 34:2 (with 16:0_18:2 chain sets). Similar PC alterations have been reported to occur in the mdx mouse model of DMD. However, altered ratios of PC 34:1 to PC 34:2 have been variably reported, and we also observed that PC 34:2 levels were nearly equally elevated as PC 34:1 in the affected mdx muscles. We hypothesized that experimental factors that often varied between studies; including muscle types sampled, mouse ages, and mouse diets; may strongly impact the PC alterations detected in dystrophic muscle of mdx mice, especially the PC 34:1 to PC 34:2 ratios. In order to test our hypothesis, we performed comprehensive lipidomic analyses of PC and phosphatidylethanolamine (PE) in several muscles (extensor digitorum longus, gastrocnemius, and soleus) and determined the mdx-specific alterations. The alterations in PC 34:1 and PC 34:2 were closely monitored from the neonate period to the adult, and also in mice raised on several diets that varied in their fats. PC 34:1 was naturally high in neonate’s muscle and decreased until age ∼3-weeks (disease onset age), and thereafter remained low in WT muscles but was higher in regenerated mdx muscles. Among the muscle types, soleus showed a distinctive phospholipid pattern with early and diminished mdx alterations. Diet was a major factor to impact PC 34:1/PC 34:2 ratios because mdx-specific alterations of PC 34:2 but not PC 34:1 were strictly dependent on diet. Our study identifies high PC 34:1 as a consistent biochemical feature of regenerated mdx-muscle and indicates nutritional approaches are also effective to modify the phospholipid compositions.

The skeletal muscle circadian clock regulates titin splicing through RBM20

Circadian rhythms are maintained by a cell-autonomous, transcriptional-translational feedback loop known as the molecular clock. While previous research suggests a role of the molecular clock in regulating skeletal muscle structure and function, no mechanisms have connected the molecular clock to sarcomere filaments. Utilizing inducible, skeletal muscle specific, Bmal1 knockout (iMSBmal1^-/-) mice, we showed that knocking out skeletal muscle clock function alters titin isoform expression using RNAseq, liquid chromatography-mass spectrometry, and sodium dodecyl sulfate-vertical agarose gel electrophoresis. This alteration in titin's spring length resulted in sarcomere length heterogeneity. We demonstrate the direct link between altered titin splicing and sarcomere length in vitro using U7 snRNPs that truncate the region of titin altered in iMSBmal1^-/- muscle. We identified a mechanism whereby the skeletal muscle clock regulates titin isoform expression through transcriptional regulation of Rbm20, a potent splicing regulator of titin. Lastly, we used an environmental model of circadian rhythm disruption and identified significant downregulation of Rbm20 expression. Our findings demonstrate the importance of the skeletal muscle circadian clock in maintaining titin isoform through regulation of RBM20 expression. Because circadian rhythm disruption is a feature of many chronic diseases, our results highlight a novel pathway that could be targeted to maintain skeletal muscle structure and function in a range of pathologies.

Epigenetics of Mitochondria-Associated Genes in Striated Muscle

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3' end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson's disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene's four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

Six1 promotes skeletal muscle thyroid hormone response through regulation of the MCT10 transporter

BACKGROUND

The Six1 transcription factor is implicated in controlling the development of several tissue types, notably skeletal muscle. Six1 also contributes to muscle metabolism and its activity is associated with the fast-twitch, glycolytic phenotype. Six1 regulates the expression of certain genes of the fast muscle program by directly stimulating their transcription or indirectly acting through a long non-coding RNA. We hypothesized that additional mechanisms of action of Six1 might be at play.

METHODS

A combined analysis of gene expression profiling and genome-wide location analysis data was performed. Results were validated using in vivo RNA interference loss-of-function assays followed by measurement of gene expression by RT-PCR and transcriptional reporter assays.

RESULTS

The Slc16a10 gene, encoding the thyroid hormone transmembrane transporter MCT10, was identified as a gene with a transcriptional enhancer directly bound by Six1 and requiring Six1 activity for full expression in adult mouse tibialis anterior, a predominantly fast-twitch muscle. Of the various thyroid hormone transporters, MCT10 mRNA was found to be the most abundant in skeletal muscle, and to have a stronger expression in fast-twitch compared to slow-twitch muscle groups. Loss-of-function of MCT10 in the tibialis anterior recapitulated the effect of Six1 on the expression of fast-twitch muscle genes and led to lower activity of a thyroid hormone receptor-dependent reporter gene.

CONCLUSIONS

These results shed light on the molecular mechanisms controlling the tissue expression profile of MCT10 and identify modulation of the thyroid hormone signaling pathway as an additional mechanism by which Six1 influences skeletal muscle metabolism.

Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle

Muscle cells have different phenotypes adapted to different usage, and can be grossly divided into fast/glycolytic and slow/oxidative types. While most muscles contain a mixture of such fiber types, we aimed at providing a genome-wide analysis of the epigenetic landscape by ChIP-Seq in two muscle extremes, the fast/glycolytic extensor digitorum longus (EDL) and slow/oxidative soleus muscles. Muscle is a heterogeneous tissue where up to 60% of the nuclei can be of a different origin. Since cellular homogeneity is critical in epigenome-wide association studies we developed a new method for purifying skeletal muscle nuclei from whole tissue, based on the nuclear envelope protein Pericentriolar material 1 (PCM1) being a specific marker for myonuclei. Using antibody labelling and a magnetic-assisted sorting approach, we were able to sort out myonuclei with 95% purity in muscles from mice, rats and humans. The sorting eliminated influence from the other cell types in the tissue and improved the myo-specific signal. A genome-wide comparison of the epigenetic landscape in EDL and soleus reflected the differences in the functional properties of the two muscles, and revealed distinct regulatory programs involving distal enhancers, including a glycolytic super-enhancer in the EDL. The two muscles were also regulated by different sets of transcription factors; e.g. in soleus, binding sites for MEF2C, NFATC2 and PPARA were enriched, while in EDL MYOD1 and SIX1 binding sites were found to be overrepresented. In addition, more novel transcription factors for muscle regulation such as members of the MAF family, ZFX and ZBTB14 were identified.

Complex tissues like skeletal muscle contain a variety of cells which confound the analysis of each cell type when based on homogenates, thus only about half of the cell nuclei in muscles reside inside the muscle cells. We here describe a labelling and sorting technique that allowed us to study the epigenetic landscape in purified muscle cell nuclei leaving the other cell types out. Differences between a fast/glycolytic and a slow/oxidative muscle were studied. While all skeletal muscle fibers have a similar make up and basic function, they differ in their physiology and the way they are used. Thus, some fibers are fast contracting but fatigable, and are used for short lasting explosive tasks such as sprinting. Other fibers are slow and are used for more prolonged tasks such as standing or long distance running. Since fiber type correlate with metabolic profile these features can also be related to metabolic diseases. We here show that the epigenetic landscape differed in gene loci corresponding to the differences in functional properties, and revealed that the two types are enriched in different gene regulatory networks. Exercise can alter muscle phenotype, and the epigenetic landscape might be related to how plastic different properties are.

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.

Gene Expression Profiling of Skeletal Muscles

Next-generation sequencing provides an opportunity for an in-depth biocomputational analysis to identify gene expression patterns between soleus and tibialis anterior, two well-characterized skeletal muscles, and analyze their gene expression profiling. RNA read counts were analyzed for differential gene expression using the R package edgeR. Differentially expressed genes were filtered using a false discovery rate of less than 0.05 c, a fold-change value of more than twenty, and an association with overrepresented pathways based on the Reactome pathway over-representation analysis tool. Most of the differentially expressed genes associated with soleus are coded for components of lipid metabolism and unique contractile elements. Differentially expressed genes associated with tibialis anterior encoded mostly for glucose and glycogen metabolic pathway regulatory enzymes and calcium-sensitive contractile components. These gene expression distinctions partly explain the genetic basis for skeletal muscle specialization, and they may help to explain skeletal muscle susceptibility to disease and drugs and further refine tissue engineering approaches.

Troponin Variants as Markers of Skeletal Muscle Health and Diseases

Ca2 +-regulated contractility is a key determinant of the quality of muscles. The sarcomeric myofilament proteins are essential players in the contraction of striated muscles. The troponin complex in the actin thin filaments plays a central role in the Ca2+-regulation of muscle contraction and relaxation. Among the three subunits of troponin, the Ca2+-binding subunit troponin C (TnC) is a member of the calmodulin super family whereas troponin I (TnI, the inhibitory subunit) and troponin T (TnT, the tropomyosin-binding and thin filament anchoring subunit) are striated muscle-specific regulatory proteins. Muscle type-specific isoforms of troponin subunits are expressed in fast and slow twitch fibers and are regulated during development and aging, and in adaptation to exercise or disuse. TnT also evolved with various alternative splice forms as an added capacity of muscle functional diversity. Mutations of troponin subunits cause myopathies. Owing to their physiological and pathological importance, troponin variants can be used as specific markers to define muscle quality. In this focused review, we will explore the use of troponin variants as markers for the fiber contents, developmental and differentiation states, contractile functions, and physiological or pathophysiological adaptations of skeletal muscle. As protein structure defines function, profile of troponin variants illustrates how changes at the myofilament level confer functional qualities at the fiber level. Moreover, understanding of the role of troponin modifications and mutants in determining muscle contractility in age-related decline of muscle function and in myopathies informs an approach to improve human health.

In Vitro, In Vivo, and In Silico Methods for Assessment of Muscle Size and Muscle Growth Regulation

Trauma, burn injury, sepsis, and ischemia lead to acute and chronic loss of skeletal muscle mass and function. Healthy muscle is essential for eating, posture, respiration, reproduction, and mobility, as well as for appropriate function of the senses including taste, vision, and hearing. Beyond providing support and contraction, skeletal muscle also exerts essential roles in temperature regulation, metabolism, and overall health. As the primary reservoir for amino acids, skeletal muscle regulates whole-body protein and glucose metabolism by providing substrate for protein synthesis and supporting hepatic gluconeogenesis during illness and starvation. Overall, greater muscle mass is linked to greater insulin sensitivity and glucose disposal, strength, power, and longevity. In contrast, low muscle mass correlates with dysmetabolism, dysmobility, and poor survival. Muscle mass is highly plastic, appropriate to its role as reservoir, and subject to striking genetic control. Defining mechanisms of muscle growth regulation holds significant promise to find interventions that promote health and diminish morbidity and mortality after trauma, sepsis, inflammation, and other systemic insults. In this invited review, we summarize techniques and methods to assess and manipulate muscle size and muscle mass in experimental systems, including cell culture and rodent models. These approaches have utility for studies of myopenia, sarcopenia, cachexia, and acute muscle growth or atrophy in the setting of health or injury.

A pig BodyMap transcriptome reveals diverse tissue physiologies and evolutionary dynamics of transcription

A comprehensive transcriptomic survey of pigs can provide a mechanistic understanding of tissue specialization processes underlying economically valuable traits and accelerate their use as a biomedical model. Here we characterize four transcript types (lncRNAs, TUCPs, miRNAs, and circRNAs) and protein-coding genes in 31 adult pig tissues and two cell lines. We uncover the transcriptomic variability among 47 skeletal muscles, and six adipose depots linked to their different origins, metabolism, cell composition, physical activity, and mitochondrial pathways. We perform comparative analysis of the transcriptomes of seven tissues from pigs and nine other vertebrates to reveal that evolutionary divergence in transcription potentially contributes to lineage-specific biology. Long-range promoter–enhancer interaction analysis in subcutaneous adipose tissues across species suggests evolutionarily stable transcription patterns likely attributable to redundant enhancers buffering gene expression patterns against perturbations, thereby conferring robustness during speciation. This study can facilitate adoption of the pig as a biomedical model for human biology and disease and uncovers the molecular bases of valuable traits.

A comprehensive transcriptomic survey of the pig could enable mechanistic understanding of tissue specialization and accelerate its use as a biomedical model. Here the authors characterize four distinct transcript types in 31 adult pig tissues to dissect their distinct structural and transcriptional features and uncover transcriptomic variability related to tissue physiology.

Myofiber stretch induces tensile and shear deformation of muscle stem cells in their native niche

Muscle stem cells (MuSCs) are requisite for skeletal muscle regeneration and homeostasis. Proper functioning of MuSCs, including activation, proliferation, and fate decision, is determined by an orchestrated series of events and communication between MuSCs and their niche. A multitude of biochemical stimuli are known to regulate MuSC fate and function. However, in addition to biochemical factors, it is conceivable that MuSCs are subjected to mechanical forces during muscle stretch-shortening cycles because of myofascial connections between MuSCs and myofibers. MuSCs respond to mechanical forces in vitro, but it remains to be proven whether physical forces are also exerted on MuSCs in their native niche and whether they contribute to the functioning and fate of MuSCs. MuSC deformation in their native niche resulting from mechanical loading of ex vivo myofiber bundles was visualized utilizing mT/mG double-fluorescent Cre-reporter mouse and multiphoton microscopy. MuSCs were subjected to 1 h pulsating fluid shear stress (PFSS) with a peak shear stress rate of 6.5 Pa/s. After PFSS treatment, nitric oxide, messenger RNA (mRNA) expression levels of genes involved in regulation of MuSC proliferation and differentiation, ERK 1/2, p38, and AKT activation were determined. Ex vivo stretching of extensor digitorum longus and soleus myofiber bundles caused compression as well as tensile and shear deformation of MuSCs in their niche. MuSCs responded to PFSS in vitro with increased nitric oxide production and an upward trend in iNOS mRNA levels. PFSS enhanced gene expression of c-Fos, Cdk4, and IL-6, whereas expression of Wnt1, MyoD, Myog, Wnt5a, COX2, Rspo1, Vangl2, Wnt10b, and MGF remained unchanged. ERK 1/2 and p38 MAPK signaling were also upregulated after PFSS treatment. We conclude that MuSCs in their native niche are subjected to force-induced deformations due to myofiber stretch-shortening. Moreover, MuSCs are mechanoresponsive, as evidenced by PFSS-mediated expression of factors by MuSCs known to promote proliferation.