Single nuclei transcriptomics of muscle reveals intra-muscular cell dynamics linked to dystrophin loss and rescue

Beyond the bulk: overview and novel insights into the dynamics of muscle satellite cells during muscle regeneration

Skeletal muscle possesses remarkable regenerative capabilities, fully recovering within a month following severe acute damage. Central to this process are muscle satellite cells (MuSCs), a resident population of somatic stem cells capable of self-renewal and differentiation. Despite the highly predictable course of muscle regeneration, evaluating this process has been challenging due to the heterogeneous nature of myogenic precursors and the limited insight provided by traditional markers with overlapping expression patterns. Notably, recent advancements in single-cell technologies, such as single-cell (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq), have revolutionized muscle research. These approaches allow for comprehensive profiling of individual cells, unveiling dynamic heterogeneity among myogenic precursors and their contributions to regeneration. Through single-cell transcriptome analyses, researchers gain valuable insights into cellular diversity and functional dynamics of MuSCs post-injury. This review aims to consolidate classical and new insights into the heterogeneity of myogenic precursors, including the latest discoveries from novel single-cell technologies.

Single-Nucleus RNA Sequencing Unravels Early Mechanisms of Human Becker Muscular Dystrophy

OBJECTIVE

The transcriptional heterogeneity at a single-nucleus level in human Becker muscular dystrophy (BMD) dystrophic muscle has not been explored. Here, we aimed to understand the transcriptional heterogeneity associated with myonuclei, as well as other mononucleated cell types that underly BMD pathogenesis by performing single-nucleus RNA sequencing.

METHODS

We profiled single-nucleus transcriptional profiles of skeletal muscle samples from 7 BMD patients and 3 normal controls.

RESULTS

A total of 17,216 nuclei (12,879 from BMD patients and 4,337 from controls) were classified into 13 known cell types, including 9 myogenic lineages and 4 non-myogenic lineages, and 1 unclassified nuclear type according to their cell identities. Among them, type IIx myonuclei were the first to degenerate in response to dystrophin reduction. Differential expression analysis revealed that the fibro-adipogenic progenitors (FAPs) population had the largest transcriptional changes among all cell types. Sub-clustering analysis identified a significantly compositional increase in the activated FAPs (aFAPs) subpopulation in BMD muscles. Pseudotime analysis, regulon inference, and deconvolution analysis of bulk RNA-sequencing data derived from 29 BMD patients revealed that the aFAPs subpopulation, a distinctive and previously unrecognized mononuclear subtype, was profibrogenic and expanded in BMD patients. Muscle quantitative real-time polymerase chain reaction and immunofluorescence analysis confirmed that the mRNA and protein levels of the aFAPs markers including LUM, DCN, and COL1A1 in BMD patients were significantly higher than those in controls, respectively.

INTERPRETATION

Our results provide insights into the transcriptional diversity of human BMD muscle at a single-nucleus resolution and new potential targets for anti-fibrosis therapies in BMD. ANN NEUROL 2024.

Caveolin and NOS in the Development of Muscular Dystrophy

Caveolin is a structural protein within caveolae that may be involved in transmembrane molecular transport and/or various intercellular interactions within cells. Specific mutations of caveolin-3 in muscle fibers are well known to cause limb-girdle muscular dystrophy. Altered expression of caveolin-3 has also been detected in Duchenne muscular dystrophy, which may be a part of the pathological process leading to muscle weakness. Interestingly, it has been shown that the renovation of nitric oxide synthase (NOS) in sarcolemma with muscular dystrophy could improve muscle health, suggesting that NOS may be involved in the pathology of muscular dystrophy. Here, we summarize the notable function of caveolin and/or NOS in skeletal muscle fibers and discuss their involvement in the pathology as well as possible tactics for the innovative treatment of muscular dystrophies.

Functional specialisation and coordination of myonuclei

Myofibres serve as the functional unit for locomotion, with the sarcomere as fundamental subunit. Running the entire length of this structure are hundreds of myonuclei, located at the periphery of the myofibre, juxtaposed to the plasma membrane. Myonuclear specialisation and clustering at the centre and ends of the fibre are known to be essential for muscle contraction, yet the molecular basis of this regionalisation has remained unclear. While the 'myonuclear domain hypothesis' helped explain how myonuclei can independently govern large cytoplasmic territories, novel technologies have provided granularity on the diverse transcriptional programs running simultaneously within the syncytia and added a new perspective on how myonuclei communicate. Building upon this, we explore the critical cellular and molecular sources of transcriptional and functional heterogeneity within myofibres, discussing the impact of intrinsic and extrinsic factors on myonuclear programs. This knowledge provides new insights for understanding muscle development, repair, and disease, but also opens avenues for the development of novel and precise therapeutic approaches.

Dystrophin deficiency impairs cell junction formation during embryonic myogenesis from pluripotent stem cells

Mutations in the DMD gene lead to Duchenne muscular dystrophy (DMD), a severe neuromuscular disorder affecting young boys as they acquire motor functions. DMD is typically diagnosed at 2-4 years of age, but the absence of dystrophin has negative impacts on skeletal muscles before overt symptoms appear in patients, which poses a serious challenge in current standards of care. Here, we investigated the consequences of dystrophin deficiency during skeletal muscle development. We used single-cell transcriptome profiling to characterize the myogenic trajectory of human pluripotent stem cells and showed that DMD cells bifurcate to an alternative branch when they reach the somite stage. Dystrophin deficiency was linked to marked dysregulations of cell junction proteins involved in the cell state transitions characteristic of embryonic somitogenesis. Altogether, this work demonstrates that in vitro, dystrophin deficiency has deleterious effects on cell-cell communication during myogenic development, which should be considered in future therapeutic strategies for DMD.

A Novel MAO-B/SSAO Inhibitor Improves Multiple Aspects of Dystrophic Phenotype in mdx Mice

Duchenne muscular dystrophy (DMD) is one of the most frequent and severe childhood muscle diseases. Its pathophysiology is multifaceted and still incompletely understood, but we and others have previously shown that oxidative stress plays an important role. In particular, we have demonstrated that inhibition of mitochondrial monoamine oxidases could improve some functional and biohumoral markers of the pathology. In the present study we report the use of dystrophic mdx mice to evaluate the efficacy of a dual monoamine oxidase B (MAO-B)/semicarbazide-sensitive amine oxidase (SSAO) inhibitor, PXS-5131, in reducing inflammation and fibrosis and improving muscle function. We found that a one-month treatment starting at three months of age was able to decrease reactive oxygen species (ROS) production, fibrosis, and inflammatory infiltrate in the tibialis anterior (TA) and diaphragm muscles. Importantly, we also observed a marked improvement in the capacity of the gastrocnemius muscle to maintain its force when challenged with eccentric contractions. Upon performing a bulk RNA-seq analysis, PXS-5131 treatment affected the expression of genes involved in inflammatory processes and tissue remodeling. We also studied the effect of prolonged treatment in older dystrophic mice, and found that a three-month administration of PXS-5131 was able to greatly reduce the progression of fibrosis not only in the diaphragm but also in the heart. Taken together, these results suggest that PXS-5131 is an effective inhibitor of fibrosis and inflammation in dystrophic muscles, a finding that could open a new therapeutic avenue for DMD patients.

Dystrophin deficiency impairs cell junction formation during embryonic myogenesis

Mutations in the DMD gene lead to Duchenne muscular dystrophy, a severe X-linked neuromuscular disorder that manifests itself as young boys acquire motor functions. DMD is typically diagnosed at 2 to 4 years of age, but the absence of dystrophin negatively impacts muscle structure and function before overt symptoms appear in patients, which poses a serious challenge in the optimization of standards of care. In this report, we investigated the early consequences of dystrophin deficiency during skeletal muscle development. We used single-cell transcriptome profiling to characterize the myogenic trajectory of human pluripotent stem cells and showed that DMD cells bifurcate to an alternative branch when they reach the somite stage. Here, dystrophin deficiency was linked to marked dysregulations of cell junction protein families involved in the cell state transitions characteristic of embryonic somitogenesis. Altogether, this work demonstrates that in vitro, dystrophin deficiency has deleterious effects on cell-cell communication during myogenic development, which should be considered in future therapeutic strategies for DMD.

Decoding the transcriptome of Duchenne muscular dystrophy to the single nuclei level reveals clinical-genetic correlations

Duchenne muscular dystrophy is a genetic disease produced by mutations in the dystrophin gene characterized by early onset muscle weakness leading to severe and irreversible disability. The cellular and molecular consequences of the lack of dystrophin in humans are only partially known, which is crucial for the development of new therapies aiming to slow or stop the progression of the disease. Here we have analyzed quadriceps muscle biopsies of seven DMD patients aged 2 to 4 years old and five age and gender matched controls using single nuclei RNA sequencing (snRNAseq) and correlated the results obtained with clinical data. SnRNAseq identified significant differences in the proportion of cell population present in the muscle samples, including an increase in the number of regenerative fibers, satellite cells, and fibro-adipogenic progenitor cells (FAPs) and a decrease in the number of slow fibers and smooth muscle cells. Muscle samples from the younger patients with stable mild weakness were characterized by an increase in regenerative fibers, while older patients with moderate and progressive weakness were characterized by loss of muscle fibers and an increase in FAPs. An analysis of the gene expression profile in muscle fibers identified a strong regenerative signature in DMD samples characterized by the upregulation of genes involved in myogenesis and muscle hypertrophy. In the case of FAPs, we observed upregulation of genes involved in the extracellular matrix regeneration but also several signaling pathways. Indeed, further analysis of the potential intercellular communication profile showed a dysregulation of the communication profile in DMD samples identifying FAPs as a key regulator of cell signaling in DMD muscle samples. In conclusion, our study has identified significant differences at the cellular and molecular levels in the different cell populations present in skeletal muscle samples of patients with DMD compared to controls.

The metabolomic plasma profile of patients with Duchenne muscular dystrophy: providing new evidence for its pathogenesis

BACKGROUND

Duchenne muscular dystrophy (DMD) is a fatal genetic muscle-wasting disease that affects 1 in 5000 male births with no current cure. Despite great progress has been made in the research of DMD, its underlying pathological mechanism based on the metabolomics is still worthy of further study. Therefore, it is necessary to gain a deeper understanding of the mechanisms or pathogenesis underlying DMD, which may reveal potential therapeutic targets and/or biomarkers.

RESULTS

Plasma samples from 42 patients with DMD from a natural history study and 40 age-matched healthy volunteers were subjected to a liquid chromatography-mass spectrometry-based non-targeted metabolomics approach. Acquired metabolic data were evaluated by principal component analysis, partial least squares-discriminant analysis, and metabolic pathway analysis to explore distinctive metabolic patterns in patients with DMD. Differentially expressed metabolites were identified using publicly available and integrated databases. By comparing the DMD and healthy control groups, 25 differential metabolites were detected, including amino acids, unsaturated fatty acids, carnitine, lipids, and metabolites related to the gut microbiota. Correspondingly, linoleic acid metabolism, D-glutamine and D-glutamate metabolism, glycerophospholipid metabolism, and alanine, aspartate, and glutamate metabolism were significantly altered in patients with DMD, compared with those of healthy volunteers.

CONCLUSIONS

Our study demonstrated the abnormal metabolism of amino acids, energy, and lipids in patients with DMD, consistent with pathological features, such as recurrent muscle necrosis and regeneration, interstitial fibrosis, and fat replacement. Additionally, we found that metabolites of intestinal flora were disordered in DMD patients, providing support for treatment of intestinal microbia disturbance in DMD diseases. Our study provides a new research strategy for understanding the pathogenesis of DMD.

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.

Muscle denervation promotes functional interactions between glial and mesenchymal cells through NGFR and NGF

We performed scRNA-seq/snATAC-seq of skeletal muscles post sciatic nerve transection to delineate cell type-specific patterns of gene expression/chromatin accessibility at different time points post-denervation. Unlike myotrauma, denervation selectively activates glial cells and Thy1/CD90-expressing mesenchymal cells. Glial cells expressed Ngf receptor (Ngfr) and were located near neuromuscular junctions (NMJs), close to Thy1/CD90-expressing cells, which provided the main cellular source of NGF post-denervation. Functional communication between these cells was mediated by NGF/NGFR, as either recombinant NGF or co-culture with Thy1/CD90-expressing cells could increase glial cell number ex vivo. Pseudo-time analysis in glial cells revealed an initial bifurcation into processes related to either cellular de-differentiation/commitment to specialized cell types (e.g., Schwann cells), or failure to promote nerve regeneration, leading to extracellular matrix remodeling toward fibrosis. Thus, interactions between denervation-activated Thy1/CD90-expressing and glial cells represent an early abortive process toward NMJs repair, ensued by the conversion of denervated muscles into an environment hostile for NMJ repair.

Trained immunity as a potential target for therapeutic immunomodulation in Duchenne muscular dystrophy

Dysregulated inflammation involving innate immune cells, particularly of the monocyte/macrophage lineage, is a key contributor to the pathogenesis of Duchenne muscular dystrophy (DMD). Trained immunity is an evolutionarily ancient protective mechanism against infection, in which epigenetic and metabolic alterations confer non-specific hyperresponsiveness of innate immune cells to various stimuli. Recent work in an animal model of DMD (mdx mice) has shown that macrophages exhibit cardinal features of trained immunity, including the presence of innate immune system "memory". The latter is reflected by epigenetic changes and durable transmissibility of the trained phenotype to healthy non-dystrophic mice by bone marrow transplantation. Mechanistically, it is suggested that a Toll-like receptor (TLR) 4-regulated, memory-like capacity of innate immunity is induced at the level of the bone marrow by factors released from the damaged muscles, leading to exaggerated upregulation of both pro- and anti-inflammatory genes. Here we propose a conceptual framework for the involvement of trained immunity in DMD pathogenesis and its potential to serve as a new therapeutic target.

Single-nucleus and bulk RNA sequencing reveal cellular and transcriptional mechanisms underlying lipid dynamics in high marbled pork

Pork is the most consumed meat in the world, and its quality is associated with human health. Intramuscular fat (IMF) deposition (also called marbling) is a key factor positively correlated with various quality traits and lipo-nutritional values of meat. However, the cell dynamics and transcriptional programs underlying lipid deposition in highly marbled meat are still unclear. Here, we used Laiwu pigs with high (HLW) or low (LLW) IMF contents to explore the cellular and transcriptional mechanisms underlying lipid deposition in highly-marbled pork by single-nucleus RNA sequencing (snRNA-seq) and bulk RNA sequencing. The HLW group had higher IMF contents but less drip loss than the LLW group. Lipidomics results revelled the changes of overall lipid classes composition (e.g., glycerolipids including triglycerides, diglycerides, and monoglycerides; sphingolipids including ceramides and monohexose ceramide significantly increased) between HLW and LLW groups. SnRNA-seq revealed nine distinct cell clusters, and the HLW group had a higher percentage of adipocytes (1.40% vs. 0.17%) than the LLW group. We identified 3 subpopulations of adipocytes, including PDE4D⁺/PDE7B⁺ (in HLW and LLW), DGAT2⁺/SCD⁺ (mostly in HLW) and FABP5⁺/SIAH1⁺ cells (mostly in HLW). Moreover, we showed that fibro/adipogenic progenitors could differentiate into IMF cells and contribute to 43.35% of adipocytes in mice. In addition, RNA-seq revealed different genes involved in lipid metabolism and fatty acid elongation. Our study provides new insights into the cellular and molecular signatures of marbling formation; such knowledge may facilitate the development of new strategies to increase IMF deposition and the lipo-nutritional quality of high marbled pork.

A protocol for single nucleus RNA-seq from frozen skeletal muscle

Single-cell technologies are a method of choice to obtain vast amounts of cell-specific transcriptional information under physiological and diseased states. Myogenic cells are resistant to single-cell RNA sequencing because of their large, multinucleated nature. Here, we report a novel, reliable, and cost-effective method to analyze frozen human skeletal muscle by single-nucleus RNA sequencing. This method yields all expected cell types for human skeletal muscle and works on tissue frozen for long periods of time and with significant pathological changes. Our method is ideal for studying banked samples with the intention of studying human muscle disease.

Childhood muscular dystrophies

Infancy- and childhood-onset muscular dystrophies are associated with a characteristic distribution and progression of motor dysfunction. The underlying causes of progressive childhood muscular dystrophies are heterogeneous involving diverse genetic pathways and genes that encode proteins of the plasma membrane, extracellular matrix, sarcomere, and nuclear membrane components. The prototypical clinicopathological features in an affected child may be adequate to fully distinguish it from other likely diagnoses based on four common features: (1) weakness and wasting of pelvic-femoral and scapular muscles with involvement of heart muscle; (2) elevation of serum muscle enzymes in particular serum creatine kinase; (3) necrosis and regeneration of myofibers; and (4) molecular neurogenetic assessment particularly utilizing next-generation sequencing of the genome of the likeliest candidates genes in an index case or family proband. A number of different animal models of therapeutic strategies have been developed for gene transfer therapy, but so far these techniques have not yet entered clinical practice. Treatment remains for the most part symptomatic with the goal of ameliorating locomotor and cardiorespiratory manifestations of the disease.