A reference single-cell transcriptomic atlas of human skeletal muscle tissue reveals bifurcated muscle stem cell populations

Single-cell RNA sequencing analysis provides novel insights into the role of apoptosis-related genes in muscle aging

OBJECTIVE

This study employed a comprehensive single-cell analysis approach to explore the role of cell apoptosis-related genes in muscle aging.

METHODS

The single-cell RNA sequencing data from the GSE143704 dataset were used to identify distinct cell clusters and assess gene expression patterns related to apoptosis activation. The "limma" package was used to identify hub genes, after which we performed Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to identify relevant pathways. Additionally, Gene Set Enrichment Analysis(GSEA) and Gene Set Variation Analysis (GSVA) were used to uncover relevant biological pathways. The Receiver Operating Characteristic Curve (ROC) was used to evaluate the diagnostic value of the hub genes. Single-sample Gene Set Enrichment Analysis (ssGSEA) was used to analyze the immune cell infiltration levels.

RESULTS

Single-cell sequencing data from muscle aging patients allowed the identification of various cell types, including epithelial cells, adipocytes, and tissue-resident macrophages. By conducting a differential expression analysis that intersected active and nonactive apoptosis, as well as comparing elderly and young samples, a total of 22 hub genes were identified (p < 0.05). The 22 hub genes have discriminative ability as potential biomarkers for diagnosing muscle aging. The enrichment analysis indicated that these genes were closely associated with diverse pathways, including "response to UV-B" and "extracellular matrix organization" (p < 0.05). Furthermore, GSEA and GSVA indicated that multiple pathways emerged-for example, the "complement and coagulation cascades", "proteasome", "insulin signaling pathway", and "MAPK signaling pathway". Additionally, the analysis of immune cell infiltration revealed positive correlations between most of the hub genes and immune cells.

CONCLUSION

Our study identified 22 apoptosis-related genes involved in muscle aging and indicated their potential diagnostic value. These findings offer a novel perspective on the pathogenesis of muscle aging and present potential targets for therapeutic interventions.

In vitro antibody-mediated SARS-CoV-2 infection suppression through human ACE2 receptor blockade

Vaccines and antibodies that specifically target or neutralize components of the SARS-CoV-2 virus are effective in prevention and treatment of human patients with SARS-CoV-2 infection. However, vaccines and SARS-CoV-2 neutralization antibodies target a subset of epitopes of viral proteins, and the fast evolution of the SARS-CoV-2 virus and the continuing emergence of SARS-CoV-2 variants confer SARS-CoV-2 immune escape from these therapies. ACE2 is the human cell receptor that serves as the entry point for SARS-CoV-2 into human cells and thus is the gatekeeper for SARS-CoV-2 infection of humans. We report here the development of 4G8C11, an anti-human ACE2 receptor monoclonal antibody that recognizes ACE2 on human cell surfaces. We determined that 4G8C11 blocks SARS-CoV-2 and variant infection of ACE2+ human cells. Furthermore, 4G8C11 has minimal effects on ACE2 receptor activity. 4G8C11 is therefore a monoclonal antibody for ACE2 receptor detection and potentially an effective immunotherapeutic agent for SARS-CoV-2 and variants.

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.

Vcam1+ Fibro-adipogenic Progenitors Mark Fatty Infiltration in Chronic Limb Threatening Ischemia

Skeletal muscle health and function is a critical determinant of clinical outcomes in patients with peripheral arterial disease (PAD). Herein, we identify fatty infiltration, the ectopic deposition of adipocytes in skeletal muscle, as a histological hallmark of end-stage PAD, also known as chronic limb threatening ischemia (CLTI). Leveraging single cell transcriptome mapping in mouse models of PAD, we identify a pro-adipogenic mesenchymal stromal cell population marked by expression of Vcam1 (termed Vcam1+ FAPs) that expands in the ischemic limb. Mechanistically, we identify Sfrp1 and Nr3c1 as regulators of Vcam1+ FAP adipogenic differentiation. Loss of Sfrp1 and Nr3c1 impair Vcam1+ FAP differentiation into adipocytes in vitro. Finally, we show that Vcam1+ FAPs are enriched in human CLTI patients. Collectively, our results identify a pro-adipogenic FAP subpopulation in CLTI patients and provide a potential therapeutic target for muscle regeneration in PAD.

Single cell RNA sequencing of human FAPs reveals different functional stages in Duchenne muscular dystrophy

Background: 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. Muscle degeneration involves a complex interplay between multiple cell lineages spatially located within areas of damage, termed the degenerative niche, including inflammatory cells, satellite cells (SCs) and fibro-adipogenic precursor cells (FAPs). FAPs are mesenchymal stem cell which have a pivotal role in muscle homeostasis as they can either promote muscle regeneration or contribute to muscle degeneration by expanding fibrotic and fatty tissue. Although it has been described that FAPs could have a different behavior in DMD patients than in healthy controls, the molecular pathways regulating their function as well as their gene expression profile are unknown. Methods: We used single-cell RNA sequencing (scRNAseq) with 10X Genomics and Illumina technology to elucidate the differences in the transcriptional profile of isolated FAPs from healthy and DMD patients. Results: Gene signatures in FAPs from both groups revealed transcriptional differences. Seurat analysis categorized cell clusters as proliferative FAPs, regulatory FAPs, inflammatory FAPs, and myofibroblasts. Differentially expressed genes (DEGs) between healthy and DMD FAPs included upregulated genes CHI3L1, EFEMP1, MFAP5, and TGFBR2 in DMD. Functional analysis highlighted distinctions in system development, wound healing, and cytoskeletal organization in control FAPs, while extracellular organization, degradation, and collagen degradation were upregulated in DMD FAPs. Validation of DEGs in additional samples (n = 9) using qPCR reinforced the specific impact of pathological settings on FAP heterogeneity, reflecting their distinct contribution to fibro or fatty degeneration in vivo. Conclusion: Using the single-cell RNA seq from human samples provide new opportunities to study cellular coordination to further understand the regulation of muscle homeostasis and degeneration that occurs in muscular dystrophies.

Transcriptomic gene signatures measure satellite cell activity in muscular dystrophies

The routine need for myonuclear turnover in skeletal muscle, together with more sporadic demands for hypertrophy and repair, are performed by resident muscle stem cells called satellite cells. Muscular dystrophies are characterized by muscle wasting, stimulating chronic repair/regeneration by satellite cells. Here, we derived and validated transcriptomic signatures for satellite cells, myoblasts/myocytes, and myonuclei using publicly available murine single cell RNA-Sequencing data. Our signatures distinguished disease from control in transcriptomic data from several muscular dystrophies including facioscapulohumeral muscular dystrophy (FSHD), Duchenne muscular dystrophy, and myotonic dystrophy type I. For FSHD, the expression of our gene signatures correlated with direct counts of satellite cells on muscle sections, as well as with increasing clinical and pathological severity. Thus, our gene signatures enable the investigation of myogenesis in bulk transcriptomic data from muscle biopsies. They also facilitate study of muscle regeneration in transcriptomic data from human muscle across health and disease.

Dual inhibition of P38 MAPK and JNK pathways preserves stemness markers and alleviates premature activation of muscle stem cells during isolation

BACKGROUND

Adult skeletal muscle contains resident muscle stem cells (MuSC) with high myogenic and engraftment potentials, making them suitable for cell therapy and regenerative medicine approaches. However, purification process of MuSC remains a major hurdle to their use in the clinic. Indeed, muscle tissue enzymatic dissociation triggers a massive activation of stress signaling pathways, among which P38 and JNK MAPK, associated with a premature loss of MuSC quiescence. While the role of these pathways in the myogenic progression of MuSC is well established, the extent to which their dissociation-induced activation affects the functionality of these cells remains unexplored.

METHODS

We assessed the effect of P38 and JNK MAPK induction on stemness marker expression and MuSC activation state during isolation by pharmacological approaches. MuSC functionality was evaluated by in vitro assays and in vivo transplantation experiments. We performed a comparative analysis of the transcriptome of human MuSC purified with pharmacological inhibitors of P38 and JNK MAPK (SB202190 and SP600125, respectively) versus available RNAseq resources.

RESULTS

We monitored PAX7 protein levels in murine MuSC during muscle dissociation and demonstrated a two-step decline partly dependent on P38 and JNK MAPK activities. We showed that simultaneous inhibition of these pathways throughout the MuSC isolation process preserves the expression of stemness markers and limits their premature activation, leading to improved survival and amplification in vitro as well as increased engraftment in vivo. Through a comparative RNAseq analysis of freshly isolated human MuSC, we provide evidence that our findings in murine MuSC could be relevant to human MuSC. Based on these findings, we implemented a purification strategy, significantly improving the recovery yields of human MuSC.

CONCLUSION

Our study highlights the pharmacological limitation of P38 and JNK MAPK activities as a suitable strategy to qualitatively and quantitatively ameliorate human MuSC purification process, which could be of great interest for cell-based therapies.

Single-nucleus RNA sequencing and lipidomics reveal characteristics of transcriptional and lipid composition in porcine longissimus dorsi muscle

BACKGROUND

Global per capita meat consumption continues to rise, especially pork. Meat quality is influenced by the content of intramuscular fat (IMF) as a key factor. The longissimus dorsi muscle of Dahe pigs (DHM, IMF: 7.98% ± 1.96%) and Dahe black pigs (DHBM, IMF: 3.30% ± 0.64%) was studied to explore cellular heterogeneity and differentially expressed genes (DEGs) associated with IMF deposition using single-nucleus RNA sequencing (snRNA-seq). The lipid composition was then analyzed using non-targeted lipidomics.

RESULTS

A total of seven cell subpopulations were identified, including myocytes, fibroblast/fibro/adipogenic progenitors (FAPs), satellite cells, endothelial cells, macrophages, pericytes, and adipocytes. Among them, FAPs and adipocytes were more focused because they could be associated with lipid deposition. 1623 DEGs in the FAPs subpopulation of DHBM were up-regulated compared with DHM, while 1535 were down-regulated. These DEGs enriched in the glycolysis/gluconeogenesis pathway. 109 DEGs were up-regulated and 806 were down-regulated in the adipocyte subpopulation of DHBM compared with DHM, which were mainly enriched in the PPAR signaling pathway and fatty acid (FA) biosynthesis. The expression level of PPARG, ABP4, LEP, and ACSL1 genes in DHM was higher than that in DHBM. Lipidomics reveals porcine lipid composition characteristics of muscle tissue. A total of 41 lipid classes and 2699 lipid species were identified in DHM and DHBM groups. The top ten relative peak areas of lipid classes in DHM and DHBM were triglyceride (TG), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), diglyceride (DG), cardiolipin (CL), ceramides (Cer), Simple Glc series (Hex1Cer), sphingomyelin (phSM), and phosphatidylinositol (PI). The relative peak areas of 35 lipid species in DHM were lower than DHBM, and 28 lipid species that were higher. There was a significant increase in the TG fatty acyl chains C6:0, C17:0, and C11:4, and a significant decrease in C16:0, C18:1, C18:2, and C22:4 in DHBM (p < 0.05).

CONCLUSIONS

C16:0 FA may downregulate the expression level of PPARG gene, which leads to the downregulation of fat metabolism-related genes such as ACSL, PLIN2, and FABP4 in DHBM compared with DHM. This may be the reason that the lipid deposition ability of Dahe pigs is stronger than that of Dahe black pigs, which need further investigation.

Fibro-adipogenic progenitors in physiological adipogenesis and intermuscular adipose tissue remodeling

Excessive accumulation of intermuscular adipose tissue (IMAT) is a common pathological feature in various metabolic and health conditions and can cause muscle atrophy, reduced function, inflammation, insulin resistance, cardiovascular issues, and unhealthy aging. Although IMAT results from fat accumulation in muscle, the mechanisms underlying its onset, development, cellular components, and functions remain unclear. IMAT levels are influenced by several factors, such as changes in the tissue environment, muscle type and origin, extent and duration of trauma, and persistent activation of fibro-adipogenic progenitors (FAPs). FAPs are a diverse and transcriptionally heterogeneous population of stromal cells essential for tissue maintenance, neuromuscular stability, and tissue regeneration. However, in cases of chronic inflammation and pathological conditions, FAPs expand and differentiate into adipocytes, resulting in the development of abnormal and ectopic IMAT. This review discusses the role of FAPs in adipogenesis and how they remodel IMAT. It highlights evidence supporting FAPs and FAP-derived adipocytes as constituents of IMAT, emphasizing their significance in adipose tissue maintenance and development, as well as their involvement in metabolic disorders, chronic pathologies and diseases. We also investigated the intricate molecular pathways and cell interactions governing FAP behavior, adipogenesis, and IMAT accumulation in chronic diseases and muscle deconditioning. Finally, we hypothesize that impaired cellular metabolic flexibility in dysfunctional muscles impacts FAPs, leading to IMAT. A deeper understanding of the biology of IMAT accumulation and the mechanisms regulating FAP behavior and fate are essential for the development of new therapeutic strategies for several debilitating conditions.

Identifying the miRNA-gene networks contributes to exploring paravertebral muscle degeneration's underlying pathogenesis and therapy strategy

Low back pain (LBP) is a worldwide problem with public health. Paravertebral muscle degeneration (PMD) is believed to be associated with LBP. Increasing evidence has demonstrated that microRNA (miRNA)-mRNA signaling networks have been implicated in the pathophysiology of diseases. Research suggests that cell death, oxidative stress, inflammatory and immune response, and extracellular matrix (ECM) metabolism are the pathogenesis of PMD; however, the miRNA-mRNA mediated the pathological process of PMD remains elusive. RNA sequencing (RNA-seq) and single cell RNA-seq (scRNA-seq) are invaluable tools for uncovering the functional biology underlying these miRNA and gene expression changes. Using scRNA-seq, we show that multiple immunocytes are presented during PMD, revealing that they may have been implicated with PMD. Additionally, using RNA-seq, we identified 76 differentially expressed genes (DEGs) and 106 differentially expressed miRNAs (DEMs), among which IL-24 and CCDC63 were the top upregulated and downregulated genes in PMD. Comprehensive bioinformatics analyses, including Venn diagrams, differential expression, functional enrichment, and protein-protein interaction analysis, were then conducted to identify six ferroptosis-related DEGs, two oxidative stress-related DEGs, eleven immunity-related DEGs, five ECM-related DEGs, among which AKR1C2/AKR1C3/SIRT1/ALB/IL-24 belong to inflammatory genes. Furthermore, 67 DEMs were predicted to be upstream miRNAs of 25 key DEGs by merging RNA-seq, TargetScan, and mirDIP databases. Finally, a miRNA-gene network was constructed using Cytoscape software and an alluvial plot. ROC curve analysis unveiled multiple key DEGs with the high clinical diagnostic value, providing novel approaches for diagnosing and treating PMD diseases.

Apelin regulates skeletal muscle adaptation to exercise in a high-intensity interval training model

Plasma apelin levels are reduced in aging and muscle wasting conditions. We aimed to investigate the significance of apelin signaling in cardiac and skeletal muscle responses to physiological stress. Apelin knockout (KO) and wild-type (WT) mice were subjected to high-intensity interval training (HIIT) by treadmill running. The effects of apelin on energy metabolism were studied in primary mouse skeletal muscle myotubes and cardiomyocytes. Apelin increased mitochondrial ATP production and mitochondrial coupling efficiency in myotubes and promoted the expression of mitochondrial genes both in primary myotubes and cardiomyocytes. HIIT induced mild concentric cardiac hypertrophy in WT mice, whereas eccentric growth was observed in the left ventricles of apelin KO mice. HIIT did not affect myofiber size in skeletal muscles of WT mice but decreased the myofiber size in apelin KO mice. The decrease in myofiber size resulted from a fiber type switch toward smaller slow-twitch type I fibers. The increased proportion of slow-twitch type I fibers in apelin KO mice was associated with upregulation of myosin heavy chain slow isoform expression, accompanied with upregulated expression of genes related to fatty acid transport and downregulated expression of genes related to glucose metabolism. Mechanistically, skeletal muscles of apelin KO mice showed defective induction of insulin-like growth factor-1 signaling in response to HIIT. In conclusion, apelin is required for proper skeletal and cardiac muscle adaptation to high-intensity exercise. Promoting apelinergic signaling may have benefits in aging- or disease-related muscle wasting conditions.NEW & NOTEWORTHY Apelin levels decline with age. This study demonstrates that in trained mice, apelin deficiency results in a switch from fast type II myofibers to slow oxidative type I myofibers. This is associated with a concomitant change in gene expression profile toward fatty acid utilization, indicating an aged-muscle phenotype in exercised apelin-deficient mice. These data are of importance in the design of exercise programs for aging individuals and could offer therapeutic target to maintain muscle mass.

Human skeletal muscle aging atlas

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Revealing myopathy spectrum: integrating transcriptional and clinical features of human skeletal muscles with varying health conditions

Myopathy refers to a large group of heterogeneous, rare muscle diseases. Bulk RNA-sequencing has been utilized for the diagnosis and research of these diseases for many years. However, the existing valuable sequencing data often lack integration and clinical interpretation. In this study, we integrated bulk RNA-sequencing data from 1221 human skeletal muscles (292 with myopathies, 929 controls) from both databases and our local samples. By applying a method similar to single-cell analysis, we revealed a general spectrum of muscle diseases, ranging from healthy to mild disease, moderate muscle wasting, and severe muscle disease. This spectrum was further partly validated in three specific myopathies (97 muscles) through clinical features including trinucleotide repeat expansion, magnetic resonance imaging fat fraction, pathology, and clinical severity scores. This spectrum helped us identify 234 genuinely healthy muscles as unprecedented controls, providing a new perspective for deciphering the hallmark genes and pathways among different myopathies. The newly identified featured genes of general myopathy, inclusion body myositis, and titinopathy were highly expressed in our local muscles, as validated by quantitative polymerase chain reaction.

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.

Cellular and Molecular Mechanisms of Heterotopic Ossification in Fibrodysplasia Ossificans Progressiva

Fibrodysplasia ossificans progressiva (FOP) is a debilitating genetic disorder characterized by recurrent episodes of heterotopic ossification (HO) formation in muscles, tendons, and ligaments. FOP is caused by a missense mutation in the ACVR1 gene (activin A receptor type I), an important signaling receptor involved in endochondral ossification. The ACVR1R206H mutation induces increased downstream canonical SMAD-signaling and drives tissue-resident progenitor cells with osteogenic potential to participate in endochondral HO formation. In this article, we review aberrant ACVR1R206H signaling and the cells that give rise to HO in FOP. FOP mouse models and lineage tracing analyses have been used to provide strong evidence for tissue-resident mesenchymal cells as cellular contributors to HO. We assess how the underlying mutation in FOP disrupts muscle-specific dynamics during homeostasis and repair, with a focus on muscle-resident mesenchymal cells known as fibro-adipogenic progenitors (FAPs). Accumulating research points to FAPs as a prominent HO progenitor population, with ACVR1R206H FAPs not only aberrantly differentiating into chondro-osteogenic lineages but creating a permissive environment for bone formation at the expense of muscle regeneration. We will further discuss the emerging role of ACVR1R206H FAPs in muscle regeneration and therapeutic targeting of these cells to reduce HO formation in FOP.

Ex vivo two-photon imaging of whole-mount skeletal muscles to visualize stem cell behavior

Quiescent skeletal muscle stem cells (MuSCs) are morphologically and functionally heterogeneous and exhibit different lengths of cellular extensions, which we call protrusions. Here, we present a protocol for ex vivo two-photon imaging of MuSCs in their native environment. We describe steps for muscle dissection, fixation, embedding, imaging, and analysis of datasets. This protocol allows the examination of MuSC morphology and protrusions at the single-cell level as well as stem cell numbers. For complete details on the use and execution of this protocol, please refer to Ma et al. (2022).1.

A transcriptomics-based drug repositioning approach to identify drugs with similar activities for the treatment of muscle pathologies in spinal muscular atrophy (SMA) models

Spinal muscular atrophy (SMA) is a genetic neuromuscular disorder caused by the reduction of survival of motor neuron (SMN) protein levels. Although three SMN-augmentation therapies are clinically approved that significantly slow down disease progression, they are unfortunately not cures. Thus, complementary SMN-independent therapies that can target key SMA pathologies and that can support the clinically approved SMN-dependent drugs are the forefront of therapeutic development. We have previously demonstrated that prednisolone, a synthetic glucocorticoid (GC) improved muscle health and survival in severe Smn-/-;SMN2 and intermediate Smn2B/- SMA mice. However, long-term administration of prednisolone can promote myopathy. We thus wanted to identify genes and pathways targeted by prednisolone in skeletal muscle to discover clinically approved drugs that are predicted to emulate prednisolone's activities. Using an RNA-sequencing, bioinformatics, and drug repositioning pipeline on skeletal muscle from symptomatic prednisolone-treated and untreated Smn-/-; SMN2 SMA and Smn+/-; SMN2 healthy mice, we identified molecular targets linked to prednisolone's ameliorative effects and a list of 580 drug candidates with similar predicted activities. Two of these candidates, metformin and oxandrolone, were further investigated in SMA cellular and animal models, which highlighted that these compounds do not have the same ameliorative effects on SMA phenotypes as prednisolone; however, a number of other important drug targets remain. Overall, our work further supports the usefulness of prednisolone's potential as a second-generation therapy for SMA, identifies a list of potential SMA drug treatments and highlights improvements for future transcriptomic-based drug repositioning studies in SMA.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Effect of aging on the human myometrium at single-cell resolution

Age-associated myometrial dysfunction can prompt complications during pregnancy and labor, which is one of the factors contributing to the 7.8-fold increase in maternal mortality in women over 40. Using single-cell/single-nucleus RNA sequencing and spatial transcriptomics, we have constructed a cellular atlas of the aging myometrium from 186,120 cells across twenty perimenopausal and postmenopausal women. We identify 23 myometrial cell subpopulations, including contractile and venous capillary cells as well as immune-modulated fibroblasts. Myometrial aging leads to fewer contractile capillary cells, a reduced level of ion channel expression in smooth muscle cells, and impaired gene expression in endothelial, smooth muscle, fibroblast, perivascular, and immune cells. We observe altered myometrial cell-to-cell communication as an aging hallmark, which associated with the loss of 25 signaling pathways, including those related to angiogenesis, tissue repair, contractility, immunity, and nervous system regulation. These insights may contribute to a better understanding of the complications faced by older individuals during pregnancy and labor.

Quality Control of Stem Cell-Based Cultured Meat According to Specific Differentiation Abilities

The demand for stem cell-based cultured meat as an alternative protein source is increasing in response to global food scarcity. However, the definition of quality controls, including appropriate growth factors and cell characteristics, remains incomplete. Cluster of differentiation (CD) 29 is ubiquitously expressed in bovine muscle tissue and is a marker of progenitor cells in cultured meat. However, CD29+ cells are naturally heterogeneous, and this quality control issue must be resolved. In this study, the aim was to identify the subpopulation of the CD29+ cell population with potential utility in cultured meat production. The CD29+ cell population exhibited heterogeneity, discernible through the CD44 and CD344 markers. CD29+CD44-CD344- cells displayed the ability for long-term culture, demonstrating high adipogenic potential and substantial lipid droplet accumulation, even within 3D cultures. Conversely, CD29+CD44+ cells exhibited rapid proliferation but were not viable for prolonged culture. Using cells suitable for adipocyte and muscle differentiation, we successfully designed meat buds, especially those rich in fat. Collectively, the identification and comprehension of distinct cell populations within bovine tissues contribute to quality control predictions in meat production. They also aid in establishing a stable and reliable cultured meat production technique.

Distinct human stem cell subpopulations drive adipogenesis and fibrosis in musculoskeletal injury

Fibroadipogenic progenitors (FAPs) maintain healthy skeletal muscle in homeostasis but drive muscle degeneration in chronic injuries by promoting adipogenesis and fibrosis. To uncover how these stem cells switch from a pro-regenerative to pro-degenerative role we perform single-cell mRNA sequencing of human FAPs from healthy and injured human muscles across a spectrum of injury, focusing on rotator cuff tears. We identify multiple subpopulations with progenitor, adipogenic, or fibrogenic gene signatures. We utilize full spectrum flow cytometry to identify distinct FAP subpopulations based on highly multiplexed protein expression. Injury severity increases adipogenic commitment of FAP subpopulations and is driven by the downregulation of DLK1. Treatment of FAPs both in vitro and in vivo with DLK1 reduces adipogenesis and fatty infiltration, suggesting that during injury, reduced DLK1 within a subpopulation of FAPs may drive degeneration. This work highlights how stem cells perform varied functions depending on tissue context, by dynamically regulating subpopulation fate commitment, which can be targeted improve patient outcomes after injury.

Mass spectrometry-based proteomics approaches to interrogate skeletal muscle adaptations to exercise

Acute exercise and chronic exercise training elicit beneficial whole-body changes in physiology that ultimately depend on profound alterations to the dynamics of tissue-specific proteins. Since the work accomplished during exercise owes predominantly to skeletal muscle, it has received the majority of interest from exercise scientists that attempt to unravel adaptive mechanisms accounting for salutary metabolic effects and performance improvements that arise from training. Contemporary scientists are also beginning to use mass spectrometry-based proteomics, which is emerging as a powerful approach to interrogate the muscle protein signature in a more comprehensive manner. Collectively, these technologies facilitate the analysis of skeletal muscle protein dynamics from several viewpoints, including changes to intracellular proteins (expression proteomics), secreted proteins (secretomics), post-translational modifications as well as fiber-, cell-, and organelle-specific changes. This review aims to highlight recent literature that has leveraged new workflows and advances in mass spectrometry-based proteomics to further our understanding of training-related changes in skeletal muscle. We call attention to untapped areas in skeletal muscle proteomics research relating to exercise training and metabolism, as well as basic points of contention when applying mass spectrometry-based analyses, particularly in the study of human biology. We further encourage researchers to couple the hypothesis-generating and descriptive nature of omics data with functional analyses that propel our understanding of the complex adaptive responses in skeletal muscle that occur with acute and chronic exercise.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Skeletal muscle regeneration failure in ischemic-damaged limbs is associated with pro-inflammatory macrophages and premature differentiation of satellite cells

BACKGROUND

Chronic limb-threatening ischemia (CLTI), a severe manifestation of peripheral arterial disease (PAD), is associated with a 1-year limb amputation rate of approximately 15-20% and substantial mortality. A key feature of CLTI is the compromised regenerative ability of skeletal muscle; however, the mechanisms responsible for this impairment are not yet fully understood. In this study, we aim to delineate pathological changes at both the cellular and transcriptomic levels, as well as in cell-cell signaling pathways, associated with compromised muscle regeneration in limb ischemia in both human tissue samples and murine models of CLTI.

METHODS

We performed single-cell transcriptome analysis of ischemic and non-ischemic muscle from the same CLTI patients and from a murine model of CLTI. In both datasets, we analyzed gene expression changes in macrophage and muscle satellite cell (MuSC) populations as well as differential cell-cell signaling interactions and differentiation trajectories.

RESULTS

Single-cell transcriptomic profiling and immunofluorescence analysis of CLTI patient skeletal muscle demonstrated that ischemic-damaged tissue displays a pro-inflammatory macrophage signature. Comparable results were observed in a murine CLTI model. Moreover, integrated analyses of both human and murine datasets revealed premature differentiation of MuSCs to be a key feature of failed muscle regeneration in the ischemic limb. Furthermore, in silico inferences of intercellular communication and in vitro assays highlight the importance of macrophage-MuSC signaling in ischemia induced muscle injuries.

CONCLUSIONS

Collectively, our research provides the first single-cell transcriptome atlases of skeletal muscle from CLTI patients and a murine CLTI model, emphasizing the crucial role of macrophages and inflammation in regulating muscle regeneration in CLTI through interactions with MuSCs.

[What if the origin of FAPs was contributing to their heterogeneity in muscle?]

Fibro-adipogenic progenitors (FAPs) are resident mesenchymal stromal cells (MSCs) of skeletal muscle. They play a crucial role in muscle homeostasis and regeneration through their paracrine activity. Recent technological advances in single-cell RNA sequencing have allowed the characterization of the heterogeneity within this cell population. In this article, we will present the different subpopulations of FAPs under basal, injury, or degenerative conditions, as well as their associated functions in mice and humans. We will then discuss the potential extramuscular origin of a post-injury FAP population. Indeed, our recent work demonstrates that MSCs from adipose tissue, infiltrating the muscle, could contribute to FAP heterogeneity.

A roadmap for delivering a human musculoskeletal cell atlas

Advances in single-cell technologies have transformed the ability to identify the individual cell types present within tissues and organs. The musculoskeletal bionetwork, part of the wider Human Cell Atlas project, aims to create a detailed map of the healthy musculoskeletal system at a single-cell resolution throughout tissue development and across the human lifespan, with complementary generation of data from diseased tissues. Given the prevalence of musculoskeletal disorders, this detailed reference dataset will be critical to understanding normal musculoskeletal function in growth, homeostasis and ageing. The endeavour will also help to identify the cellular basis for disease and lay the foundations for novel therapeutic approaches to treating diseases of the joints, soft tissues and bone. Here, we present a Roadmap delineating the critical steps required to construct the first draft of a human musculoskeletal cell atlas. We describe the key challenges involved in mapping the extracellular matrix-rich, but cell-poor, tissues of the musculoskeletal system, outline early milestones that have been achieved and describe the vision and directions for a comprehensive musculoskeletal cell atlas. By embracing cutting-edge technologies, integrating diverse datasets and fostering international collaborations, this endeavour has the potential to drive transformative changes in musculoskeletal medicine.

Senescent fibro-adipogenic progenitors are potential drivers of pathology in inclusion body myositis

Inclusion body myositis (IBM) is unique across the spectrum of idiopathic inflammatory myopathies (IIM) due to its distinct clinical presentation and refractoriness to current treatment approaches. One explanation for this resistance may be the engagement of cell-autonomous mechanisms that sustain or promote disease progression of IBM independent of inflammatory activity. In this study, we focused on senescence of tissue-resident cells as potential driver of disease. For this purpose, we compared IBM patients to non-diseased controls and immune-mediated necrotizing myopathy patients. Histopathological analysis suggested that cellular senescence is a prominent feature of IBM, primarily affecting non-myogenic cells. In-depth analysis by single nuclei RNA sequencing allowed for the deconvolution and study of muscle-resident cell populations. Among these, we identified a specific cluster of fibro-adipogenic progenitors (FAPs) that demonstrated key hallmarks of senescence, including a pro-inflammatory secretome, expression of p21, increased β-galactosidase activity, and engagement of senescence pathways. FAP function is required for muscle cell health with changes to their phenotype potentially proving detrimental. In this respect, the transcriptomic landscape of IBM was also characterized by changes to the myogenic compartment demonstrating a pronounced loss of type 2A myofibers and a rarefication of acetylcholine receptor expressing myofibers. IBM muscle cells also engaged a specific pro-inflammatory phenotype defined by intracellular complement activity and the expression of immunogenic surface molecules. Skeletal muscle cell dysfunction may be linked to FAP senescence by a change in the collagen composition of the latter. Senescent FAPs lose collagen type XV expression, which is required to support myofibers' structural integrity and neuromuscular junction formation in vitro. Taken together, this study demonstrates an altered phenotypical landscape of muscle-resident cells and that FAPs, and not myofibers, are the primary senescent cell type in IBM.

An efficient protocol for studying human pluripotent stem cell-derived myotube senescence

Sarcopenia, an age-related skeletal muscle condition characterized by a progressive decline in muscle mass and function, is linked to increased vulnerability, a higher likelihood of falls, and higher mortality rates in older individuals. A comprehensive understanding of the intricate mechanisms driving skeletal muscle aging is of great significance in both scientific and clinical fields. Consequently, myotube models that facilitate studying regulatory mechanisms underlying skeletal muscle aging are important tools required to advance intervention strategies against skeletal muscle aging and associated disorders. Here, we provide a detailed protocol to generate human pluripotent stem cells-derived myotubes and describe their applications in aging studies, as well as a troubleshooting for potential problems. Overall, this protocol serves as a valuable methodological reference for exploring the role and mechanism of genes involved in skeletal muscle aging.

Single-Nuclei RNA-Sequencing of the Gastrocnemius Muscle in Peripheral Artery Disease

BACKGROUND

Lower extremity peripheral artery disease (PAD) is a growing epidemic with limited effective treatment options. Here, we provide a single-nuclei atlas of PAD limb muscle to facilitate a better understanding of the composition of cells and transcriptional differences that comprise the diseased limb muscle.

METHODS

We obtained gastrocnemius muscle specimens from 20 patients with PAD and 12 non-PAD controls. Nuclei were isolated and single-nuclei RNA-sequencing was performed. The composition of nuclei was characterized by iterative clustering via principal component analysis, differential expression analysis, and the use of known marker genes. Bioinformatics analysis was performed to determine differences in gene expression between PAD and non-PAD nuclei, as well as subsequent analysis of intercellular signaling networks. Additional histological analyses of muscle specimens accompany the single-nuclei RNA-sequencing atlas.

RESULTS

Single-nuclei RNA-sequencing analysis indicated a fiber type shift with patients with PAD having fewer type I (slow/oxidative) and more type II (fast/glycolytic) myonuclei compared with non-PAD, which was confirmed using immunostaining of muscle specimens. Myonuclei from PAD displayed global upregulation of genes involved in stress response, autophagy, hypoxia, and atrophy. Subclustering of myonuclei also identified populations that were unique to PAD muscle characterized by metabolic dysregulation. PAD muscles also displayed unique transcriptional profiles and increased diversity of transcriptomes in muscle stem cells, regenerating myonuclei, and fibro-adipogenic progenitor cells. Analysis of intercellular communication networks revealed fibro-adipogenic progenitors as a major signaling hub in PAD muscle, as well as deficiencies in angiogenic and bone morphogenetic protein signaling which may contribute to poor limb function in PAD.

CONCLUSIONS

This reference single-nuclei RNA-sequencing atlas provides a comprehensive analysis of the cell composition, transcriptional signature, and intercellular communication pathways that are altered in the PAD condition.

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.

White Adipose Tissue Heterogeneity in the Single-Cell Era: From Mice and Humans to Cattle

Adipose tissue is a major modulator of metabolic function by regulating energy storage and by acting as an endocrine organ through the secretion of adipokines. With the advantage of next-generation sequencing-based single-cell technologies, adipose tissue has been studied at single-cell resolution, thus providing unbiased insight into its molecular composition. Recent single-cell RNA sequencing studies in human and mouse models have dissected the transcriptional cellular heterogeneity of subcutaneous (SAT), visceral (VAT), and intramuscular (IMAT) white adipose tissue depots and revealed unique populations of adipose tissue progenitor cells, mature adipocytes, immune cell, vascular cells, and mesothelial cells that play direct roles on adipose tissue function and the development of metabolic disorders. In livestock species, especially in bovine, significant gaps of knowledge remain in elucidating the roles of adipose tissue cell types and depots on driving the pathogenesis of metabolic disorders and the distinct fat deposition in VAT, SAT, and IMAT in meat animals. This review summarizes the current knowledge on the transcriptional and functional cellular diversity of white adipose tissue revealed by single-cell approaches and highlights the depot-specific function of adipose tissue in different mammalian species, with a particular focus on recent findings and future implications in cattle.

Single-cell analysis of bovine muscle-derived cell types for cultured meat production

Cultured meat technologies leverage the proliferation and differentiation of animal-derived stem cells ex vivo to produce edible tissues for human consumption in a sustainable fashion. However, skeletal muscle is a dynamic and highly complex tissue, involving the interplay of numerous mono- and multinucleated cells, including muscle fibers, satellite cells (SCs) and fibro-adipogenic progenitors (FAPs), and recreation of the tissue in vitro thus requires the characterization and manipulation of a broad range of cell types. Here, we use a single-cell RNA sequencing approach to characterize cellular heterogeneity within bovine muscle and muscle-derived cell cultures over time. Using this data, we identify numerous distinct cell types, and develop robust protocols for the easy purification and proliferation of several of these populations. We note overgrowth of undesirable cell types within heterogeneous proliferative cultures as a barrier to efficient cultured meat production, and use transcriptomics to identify conditions that favor the growth of SCs in the context of serum-free medium. Combining RNA velocities computed in silico with time-resolved flow cytometric analysis, we characterize dynamic subpopulations and transitions between active, quiescent, and committed states of SCs, and demonstrate methods for modulation of these states during long-term proliferative cultures. This work provides an important reference for advancing our knowledge of bovine skeletal muscle biology, and its application in the development of cultured meat technologies.

In vitro-generated human muscle reserve cells are heterogeneous for Pax7 with distinct molecular states and metabolic profiles

BACKGROUND

The capacity of skeletal muscles to regenerate relies on Pax7+ muscle stem cells (MuSC). While in vitro-amplified MuSC are activated and lose part of their regenerative capacity, in vitro-generated human muscle reserve cells (MuRC) are very similar to quiescent MuSC with properties required for their use in cell-based therapies.

METHODS

In the present study, we investigated the heterogeneity of human MuRC and characterized their molecular signature and metabolic profile.

RESULTS

We observed that Notch signaling is active and essential for the generation of quiescent human Pax7+ MuRC in vitro. We also revealed, by immunofluorescence and flow cytometry, two distinct subpopulations of MuRC distinguished by their relative Pax7 expression. After 48 h in differentiation medium (DM), the Pax7High subpopulation represented 35% of the total MuRC pool and this percentage increased to 61% after 96 h in DM. Transcriptomic analysis revealed that Pax7High MuRC were less primed for myogenic differentiation as compared to Pax7Low MuRC and displayed a metabolic shift from glycolysis toward fatty acid oxidation. The bioenergetic profile of human MuRC displayed a 1.5-fold decrease in glycolysis, basal respiration and ATP-linked respiration as compared to myoblasts. We also observed that AMPKα1 expression was significantly upregulated in human MuRC that correlated with an increased phosphorylation of acetyl-CoA carboxylase (ACC). Finally, we showed that fatty acid uptake was increased in MuRC as compared to myoblasts, whereas no changes were observed for glucose uptake.

CONCLUSIONS

Overall, these data reveal that the quiescent MuRC pool is heterogeneous for Pax7 with a Pax7High subpopulation being in a deeper quiescent state, less committed to differentiation and displaying a reduced metabolic activity. Altogether, our data suggest that human Pax7High MuRC may constitute an appropriate stem cell source for potential therapeutic applications in skeletal muscle diseases.

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.

A single-cell transcriptome atlas of pig skin characterizes anatomical positional heterogeneity

Different anatomical locations of the body skin show differences in their gene expression patterns depending on different origins, and the inherent heterogeneous information can be maintained in adults. However, highly resolvable cellular specialization is less well characterized in different anatomical regions of the skin. Pig is regarded as an excellent model animal for human skin research in view of its similar physiology to human. In this study, single-cell RNA sequencing was performed on pig skin tissues from six different anatomical regions of Chenghua (CH) pigs, with a superior skin thickness trait, and the back site of large white (LW) pigs. We obtained 233,715 cells, representing seven cell types, among which we primarily characterized the heterogeneity of the top three cell types, including smooth muscle cells (SMCs), endothelial cells (ECs), and fibroblasts (FBs). Then, we further identified several subtypes of SMCs, ECs, and FBs, and discovered the expression patterns of site-specific genes involved in some important pathways such as the immune response and extracellular matrix (ECM) synthesis in different anatomical regions. By comparing differentially expressed genes of skin FBs among different anatomical regions, we considered TNN, COL11A1, and INHBA as candidate genes for facilitating ECM accumulation. These findings of heterogeneity in the main skin cell types from different anatomical sites will contribute to a better understanding of inherent skin information and place the potential focus on skin generation, transmission, and transplantation, paving the foundation for human skin priming.

Single-cell RNA sequencing in skeletal muscle developmental biology

Skeletal muscle is the most extensive tissue in mammals, and they perform several functions; it is derived from paraxial mesodermal somites and undergoes hyperplasia and hypertrophy to form multinucleated, contractile, and functional muscle fibers. Skeletal muscle is a complex heterogeneous tissue composed of various cell types that establish communication strategies to exchange biological information; therefore, characterizing the cellular heterogeneity and transcriptional signatures of skeletal muscle is central to understanding its ontogeny's details. Studies of skeletal myogenesis have focused primarily on myogenic cells' proliferation, differentiation, migration, and fusion and ignored the intricate network of cells with specific biological functions. The rapid development of single-cell sequencing technology has recently enabled the exploration of skeletal muscle cell types and molecular events during development. This review summarizes the progress in single-cell RNA sequencing and its applications in skeletal myogenesis, which will provide insights into skeletal muscle pathophysiology.

Biomarkers of aging

Aging biomarkers are a combination of biological parameters to (i) assess age-related changes, (ii) track the physiological aging process, and (iii) predict the transition into a pathological status. Although a broad spectrum of aging biomarkers has been developed, their potential uses and limitations remain poorly characterized. An immediate goal of biomarkers is to help us answer the following three fundamental questions in aging research: How old are we? Why do we get old? And how can we age slower? This review aims to address this need. Here, we summarize our current knowledge of biomarkers developed for cellular, organ, and organismal levels of aging, comprising six pillars: physiological characteristics, medical imaging, histological features, cellular alterations, molecular changes, and secretory factors. To fulfill all these requisites, we propose that aging biomarkers should qualify for being specific, systemic, and clinically relevant.

microRNA-501 controls myogenin+/CD74+ myogenic progenitor cells during muscle regeneration

OBJECTIVE

Skeletal muscle regeneration is markedly impaired during aging. How adult muscle stem cells contribute to this decrease in regenerative capacity is incompletely understood. We investigated mechanisms of age-related changes in myogenic progenitor cells using the tissue-specific microRNA 501.

METHODS

Young and old C57Bl/6 mice were used (3 months or 24 months of age, respectively) with or without global or tissue-specific genetic deletion of miR-501. Muscle regeneration was induced using intramuscular cardiotoxin injection or treadmill exercise and analysed using single cell and bulk RNA sequencing, qRT-PCR and immunofluorescence. Muscle fiber damage was assessed with Evan`s blue dye (EBD). In vitro analysis was performed in primary muscle cells obtained from mice and humans.

RESULTS

Single cell sequencing revealed myogenic progenitor cells in miR-501 knockout mice at day 6 after muscle injury that are characterized by high levels of myogenin and CD74. In control mice these cells were less in number and already downregulated after day 3 of muscle injury. Muscle from knockout mice had reduced myofiber size and reduced myofiber resilience to injury and exercise. miR-501 elicits this effect by regulating sarcomeric gene expression through its target gene estrogen-related receptor gamma (Esrrg). Importantly, in aged skeletal muscle where miR-501 was significantly downregulated and its target Esrrg significantly upregulated, the number of myog+/CD74+ cells during regeneration was upregulated to similar levels as observed in 501 knockout mice. Moreover, myog+/CD74+-aged skeletal muscle exhibited a similar decrease in the size of newly formed myofibers and increased number of necrotic myofibers after injury as observed in mice lacking miR-501.

CONCLUSIONS

miR-501 and Esrrg are regulated in muscle with decreased regenerative capacity and loss of miR-501 is permissive to the appearance of CD74+ myogenic progenitors. Our data uncover a novel link between the metabolic transcription factor Esrrg and sarcomere formation and demonstrate that stem cell heterogeneity in skeletal muscle during aging is under miRNA control. Targeting Esrrg or myog+/CD74+ progenitor cells might improve fiber size and myofiber resilience to exercise in aged skeletal muscle.

A human skeletal muscle stem/myotube model reveals multiple signaling targets of cancer secretome in skeletal muscle

Skeletal muscle dysfunction or reprogramming due to the effects of the cancer secretome is observed in multiple malignancies. Although mouse models are routinely used to study skeletal muscle defects in cancer, because of species specificity of certain cytokines/chemokines in the secretome, a human model system is required. Here, we establish simplified multiple skeletal muscle stem cell lines (hMuSCs), which can be differentiated into myotubes. Using single nuclei ATAC-seq (snATAC-seq) and RNA-seq (snRNA-seq), we document chromatin accessibility and transcriptomic changes associated with the transition of hMuSCs to myotubes. Cancer secretome accelerated stem to myotube differentiation, altered the alternative splicing machinery and increased inflammatory, glucocorticoid receptor, and wound healing pathways in hMuSCs. Additionally, cancer secretome reduced metabolic and survival pathway associated miR-486, AKT, and p53 signaling in hMuSCs. hMuSCs underwent myotube differentiation when engrafted into NSG mice and thus providing a humanized in vivo skeletal muscle model system to study cancer cachexia.

Odd skipped-related 1 controls the pro-regenerative response of fibro-adipogenic progenitors

Skeletal muscle regeneration requires the coordinated interplay of diverse tissue-resident- and infiltrating cells. Fibro-adipogenic progenitors (FAPs) are an interstitial cell population that provides a beneficial microenvironment for muscle stem cells (MuSCs) during muscle regeneration. Here we show that the transcription factor Osr1 is essential for FAPs to communicate with MuSCs and infiltrating macrophages, thus coordinating muscle regeneration. Conditional inactivation of Osr1 impaired muscle regeneration with reduced myofiber growth and formation of excessive fibrotic tissue with reduced stiffness. Osr1-deficient FAPs acquired a fibrogenic identity with altered matrix secretion and cytokine expression resulting in impaired MuSC viability, expansion and differentiation. Immune cell profiling suggested a novel role for Osr1-FAPs in macrophage polarization. In vitro analysis suggested that increased TGFβ signaling and altered matrix deposition by Osr1-deficient FAPs actively suppressed regenerative myogenesis. In conclusion, we show that Osr1 is central to FAP function orchestrating key regenerative events such as inflammation, matrix secretion and myogenesis.

Pro-inflammatory macrophages impair skeletal muscle regeneration in ischemic-damaged limbs by inducing precocious differentiation of satellite cells

Chronic limb-threatening ischemia (CLTI), representing the end-stage of peripheral arterial disease (PAD), is associated with a one-year limb amputation rate of ∼15-20% and significant mortality. A key characteristic of CLTI is the failure of the innate regenerative capacity of skeletal muscle, though the underlying mechanisms remain unclear. Here, single-cell transcriptome analysis of ischemic and non-ischemic muscle from the same CLTI patients demonstrated that ischemic-damaged tissue is enriched with pro-inflammatory macrophages. Comparable results were also observed in a murine CLTI model. Importantly, integrated analyses of both human and murine data revealed premature differentiation of muscle satellite cells (MuSCs) in damaged tissue and indications of defects in intercellular signaling communication between MuSCs and their inflammatory niche. Collectively, our research provides the first single-cell transcriptome atlases of skeletal muscle from CLTI patients and murine models, emphasizing the crucial role of macrophages and inflammation in regulating muscle regeneration in CLTI through interactions with MuSCs.

Single cell analysis reveals satellite cell heterogeneity for proinflammatory chemokine expression

Background: The expression of proinflammatory signals at the site of muscle injury are essential for efficient tissue repair and their dysregulation can lead to inflammatory myopathies. Macrophages, neutrophils, and fibroadipogenic progenitor cells residing in the muscle are significant sources of proinflammatory cytokines and chemokines. However, the inducibility of the myogenic satellite cell population and their contribution to proinflammatory signaling is less understood.

Methods: Mouse satellite cells were isolated and exposed to lipopolysaccharide (LPS) to mimic sterile skeletal muscle injury and changes in the expression of proinflammatory genes was examined by RT-qPCR and single cell RNA sequencing. Expression patterns were validated in skeletal muscle injured with cardiotoxin by RT-qPCR and immunofluorescence.

Results: Satellite cells in culture were able to express Tnfa, Ccl2, and Il6, within 2 h of treatment with LPS. Single cell RNA-Seq revealed seven cell clusters representing the continuum from activation to differentiation. LPS treatment led to a heterogeneous pattern of induction of C-C and C-X-C chemokines (e.g., Ccl2, Ccl5, and Cxcl0) and cytokines (e.g., Tgfb1, Bmp2, Il18, and Il33) associated with innate immune cell recruitment and satellite cell proliferation. One cell cluster was enriched for expression of the antiviral interferon pathway genes under control conditions and LPS treatment. Activation of this pathway in satellite cells was also detectable at the site of cardiotoxin induced muscle injury.

Conclusion: These data demonstrate that satellite cells respond to inflammatory signals and secrete chemokines and cytokines. Further, we identified a previously unrecognized subset of satellite cells that may act as sensors for muscle infection or injury using the antiviral interferon pathway.

Mechanisms of skeletal muscle-tendon development and regeneration/healing as potential therapeutic targets

Skeletal muscle contraction is essential for the movement of our musculoskeletal system. Tendons and ligaments that connect the skeletal muscles to bones in the correct position at the appropriate time during development are also required for movement to occur. Since the musculoskeletal system is essential for maintaining basic bodily functions as well as enabling interactions with the environment, dysfunctions of these tissues due to disease can significantly reduce quality of life. Unfortunately, as people live longer, skeletal muscle and tendon/ligament diseases are becoming more common. Sarcopenia, a disease in which skeletal muscle function declines, and tendinopathy, which involves chronic tendon dysfunction, are particularly troublesome because there have been no significant advances in their treatment. In this review, we will summarize previous reports on the development and regeneration/healing of skeletal muscle and tendon tissues, including a discussion of the molecular and cellular mechanisms involved that may be used as potential therapeutic targets.

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.

MME+ fibro-adipogenic progenitors are the dominant adipogenic population during fatty infiltration in human skeletal muscle

Fatty infiltration, the ectopic deposition of adipose tissue within skeletal muscle, is mediated via the adipogenic differentiation of fibro-adipogenic progenitors (FAPs). We used single-nuclei and single-cell RNA sequencing to characterize FAP heterogeneity in patients with fatty infiltration. We identified an MME⁺ FAP subpopulation which, based on ex vivo characterization as well as transplantation experiments, exhibits high adipogenic potential. MME⁺ FAPs are characterized by low activity of WNT, known to control adipogenic commitment, and are refractory to the inhibitory role of WNT activators. Using preclinical models for muscle damage versus fatty infiltration, we show that many MME⁺ FAPs undergo apoptosis during muscle regeneration and differentiate into adipocytes under pathological conditions, leading to a reduction in their abundance. Finally, we utilized the varying fat infiltration levels in human hip muscles and found less MME⁺ FAPs in fatty infiltrated human muscle. Altogether, we have identified the dominant adipogenic FAP subpopulation in skeletal muscle.

Inflammation balance in skeletal muscle damage and repair

Responding to tissue injury, skeletal muscles undergo the tissue destruction and reconstruction accompanied with inflammation. The immune system recognizes the molecules released from or exposed on the damaged tissue. In the local minor tissue damage, tissue-resident macrophages sequester pro-inflammatory debris to prevent initiation of inflammation. In most cases of the skeletal muscle injury, however, a cascade of inflammation will be initiated through activation of local macrophages and mast cells and recruitment of immune cells from blood circulation to the injured site by recongnization of damage-associated molecular patterns (DAMPs) and activated complement system. During the inflammation, macrophages and neutrophils scavenge the tissue debris to release inflammatory cytokines and the latter stimulates myoblast fusion and vascularization to promote injured muscle repair. On the other hand, an abundance of released inflammatory cytokines and chemokines causes the profound hyper-inflammation and mobilization of immune cells to trigger a vicious cycle and lead to the cytokine storm. The cytokine storm results in the elevation of cytolytic and cytotoxic molecules and reactive oxygen species (ROS) in the damaged muscle to aggravates the tissue injury, including the healthy bystander tissue. Severe inflammation in the skeletal muscle can lead to rhabdomyolysis and cause sepsis-like systemic inflammation response syndrome (SIRS) and remote organ damage. Therefore, understanding more details on the involvement of inflammatory factors and immune cells in the skeletal muscle damage and repair can provide the new precise therapeutic strategies, including attenuation of the muscle damage and promotion of the muscle repair.

Application of Single-Cell and Spatial Omics in Musculoskeletal Disorder Research

Musculoskeletal disorders, including fractures, scoliosis, heterotopic ossification, osteoporosis, osteoarthritis, disc degeneration, and muscular injury, etc., can occur at any stage of human life. Understanding the occurrence and development mechanism of musculoskeletal disorders, as well as the changes in tissues and cells during therapy, might help us find targeted treatment methods. Single-cell techniques provide excellent tools for studying alterations at the cellular level of disorders. However, the application of these techniques in research on musculoskeletal disorders is still limited. This review summarizes the current single-cell and spatial omics used in musculoskeletal disorders. Cell isolation, experimental methods, and feasible experimental designs for single-cell studies of musculoskeletal system diseases have been reviewed based on tissue characteristics. Then, the paper summarizes the latest findings of single-cell studies in musculoskeletal disorders from three aspects: bone and ossification, joint, and muscle and tendon disorders. Recent discoveries about the cell populations involved in these diseases are highlighted. Furthermore, the therapeutic responses of musculoskeletal disorders, especially single-cell changes after the treatments of implants, stem cell therapies, and drugs are described. Finally, the application potential and future development directions of single-cell and spatial omics in research on musculoskeletal diseases are discussed.

Multi-omic integration reveals cell-type-specific regulatory networks of insulin resistance in distinct ancestry populations

Our knowledge of the cell-type-specific mechanisms of insulin resistance remains limited. To dissect the cell-type-specific molecular signatures of insulin resistance, we performed a multiscale gene network analysis of adipose and muscle tissues in African and European ancestry populations. In adipose tissues, a comparative analysis revealed ethnically conserved cell-type signatures and two adipocyte subtype-enriched modules with opposite insulin sensitivity responses. The modules enriched for adipose stem and progenitor cells as well as immune cells showed negative correlations with insulin sensitivity. In muscle tissues, the modules enriched for stem cells and fibro-adipogenic progenitors responded to insulin sensitivity oppositely. The adipocyte and muscle fiber-enriched modules shared cellular-respiration-related genes but had tissue-specific rearrangements of gene regulations in response to insulin sensitivity. Integration of the gene co-expression and causal networks further pinpointed key drivers of insulin resistance. Together, this study revealed the cell-type-specific transcriptomic networks and signaling maps underlying insulin resistance in major glucose-responsive tissues. A record of this paper's transparent peer review process is included in the supplemental information.

Integrative genetic and single cell RNA sequencing analysis provides new clues to the amyotrophic lateral sclerosis neurodegeneration

Introduction

The gradual loss of motor neurons (MNs) in the brain and spinal cord is a hallmark of amyotrophic lateral sclerosis (ALS), but the mechanisms underlying neurodegeneration in ALS are still not fully understood.

Methods

Based on 75 ALS-pathogenicity/susceptibility genes and large-scale single-cell transcriptomes of human/mouse brain/spinal cord/muscle tissues, we performed an expression enrichment analysis to identify cells involved in ALS pathogenesis. Subsequently, we created a strictness measure to estimate the dosage requirement of ALS-related genes in linked cell types.

Results

Remarkably, expression enrichment analysis showed that α- and γ-MNs, respectively, are associated with ALS-susceptibility genes and ALS-pathogenicity genes, revealing differences in biological processes between sporadic and familial ALS. In MNs, ALS-susceptibility genes exhibited high strictness, as well as the ALS-pathogenicity genes with known loss of function mechanism, indicating the main characteristic of ALS-susceptibility genes is dosage-sensitive and the loss of function mechanism of these genes may involve in sporadic ALS. In contrast, ALS-pathogenicity genes with gain of function mechanism exhibited low strictness. The significant difference of strictness between loss of function genes and gain of function genes provided a priori understanding for the pathogenesis of novel genes without an animal model. Besides MNs, we observed no statistical evidence for an association between muscle cells and ALS-related genes. This result may provide insight into the etiology that ALS is not within the domain of neuromuscular diseases. Moreover, we showed several cell types linked to other neurological diseases [i.e., spinocerebellar ataxia (SA), hereditary motor neuropathies (HMN)] and neuromuscular diseases [i.e. hereditary spastic paraplegia (SPG), spinal muscular atrophy (SMA)], including an association between Purkinje cells in brain and SA, an association between α-MNs in spinal cord and SA, an association between smooth muscle cells and SA, an association between oligodendrocyte and HMN, a suggestive association between γ-MNs and HMN, a suggestive association between mature skeletal muscle and HMN, an association between oligodendrocyte in brain and SPG, and no statistical evidence for an association between cell type and SMA.

Discussion

These cellular similarities and differences deepened our understanding of the heterogeneous cellular basis of ALS, SA, HMN, SPG, and SMA.

Single nuclei profiling identifies cell specific markers of skeletal muscle aging, frailty, and senescence

Aging is accompanied by a loss of muscle mass and function, termed sarcopenia, which causes numerous morbidities and economic burdens in human populations. Mechanisms implicated in age-related sarcopenia or frailty include inflammation, muscle stem cell depletion, mitochondrial dysfunction, and loss of motor neurons, but whether there are key drivers of sarcopenia are not yet known. To gain deeper insights into age-related muscle loss, we performed transcriptome profiling on lower limb muscle biopsies from 72 young, elderly, and frail human subjects using bulk RNA-seq (N = 72) and single-nuclei RNA-seq (N = 17). This combined approach revealed changes in gene expression that occur with age and frailty in multiple cell types comprising mature skeletal muscle. Notably, we found increased expression of the genes MYH8 and PDK4, and decreased expression of the gene IGFN1, in aged muscle. We validated several key genes changes in fixed human muscle tissue using digital spatial profiling. We also identified a small population of nuclei that express CDKN1A, present only in aged samples, consistent with p21cip1-driven senescence in this subpopulation. Overall, our findings identify unique cellular subpopulations in aged and sarcopenic skeletal muscle, which will facilitate the development of new therapeutic strategies to combat age-related frailty.