Large-scale integration of single-cell transcriptomic data captures transitional progenitor states in mouse skeletal muscle regeneration

Current cutting-edge omics techniques on musculoskeletal tissues and diseases

Musculoskeletal disorders, including osteoarthritis, rheumatoid arthritis, osteoporosis, bone fracture, intervertebral disc degeneration, tendinopathy, and myopathy, are prevalent conditions that profoundly impact quality of life and place substantial economic burdens on healthcare systems. Traditional bulk transcriptomics, genomics, proteomics, and metabolomics have played a pivotal role in uncovering disease-associated alterations at the population level. However, these approaches are inherently limited in their ability to resolve cellular heterogeneity or to capture the spatial organization of cells within tissues, thus hindering a comprehensive understanding of the complex cellular and molecular mechanisms underlying these diseases. To address these limitations, advanced single-cell and spatial omics techniques have emerged in recent years, offering unparalleled resolution for investigating cellular diversity, tissue microenvironments, and biomolecular interactions within musculoskeletal tissues. These cutting-edge techniques enable the detailed mapping of the molecular landscapes in diseased tissues, providing transformative insights into pathophysiological processes at both the single-cell and spatial levels. This review presents a comprehensive overview of the latest omics technologies as applied to musculoskeletal research, with a particular focus on their potential to revolutionize our understanding of disease mechanisms. Additionally, we explore the power of multi-omics integration in identifying novel therapeutic targets and highlight key challenges that must be overcome to successfully translate these advancements into clinical applications.

Single-cell and spatial omics: exploring hypothalamic heterogeneity

Elucidating the complex dynamic cellular organization in the hypothalamus is critical for understanding its role in coordinating fundamental body functions. Over the past decade, single-cell and spatial omics technologies have significantly evolved, overcoming initial technical challenges in capturing and analyzing individual cells. These high-throughput omics technologies now offer a remarkable opportunity to comprehend the complex spatiotemporal patterns of transcriptional diversity and cell-type characteristics across the entire hypothalamus. Current single-cell and single-nucleus RNA sequencing methods comprehensively quantify gene expression by exploring distinct phenotypes across various subregions of the hypothalamus. However, single-cell/single-nucleus RNA sequencing requires isolating the cell/nuclei from the tissue, potentially resulting in the loss of spatial information concerning neuronal networks. Spatial transcriptomics methods, by bypassing the cell dissociation, can elucidate the intricate spatial organization of neural networks through their imaging and sequencing technologies. In this review, we highlight the applicative value of single-cell and spatial transcriptomics in exploring the complex molecular-genetic diversity of hypothalamic cell types, driven by recent high-throughput achievements.

Advances in spatial transcriptomics and its application in the musculoskeletal system

While bulk RNA sequencing and single-cell RNA sequencing have shed light on cellular heterogeneity and potential molecular mechanisms in the musculoskeletal system in both physiological and various pathological states, the spatial localization of cells and molecules and intercellular interactions within the tissue context require further elucidation. Spatial transcriptomics has revolutionized biological research by simultaneously capturing gene expression profiles and in situ spatial information of tissues, gradually finding applications in musculoskeletal research. This review provides a summary of recent advances in spatial transcriptomics and its application to the musculoskeletal system. The classification and characteristics of data acquisition techniques in spatial transcriptomics are briefly outlined, with an emphasis on widely-adopted representative technologies and the latest technological breakthroughs, accompanied by a concise workflow for incorporating spatial transcriptomics into musculoskeletal system research. The role of spatial transcriptomics in revealing physiological mechanisms of the musculoskeletal system, particularly during developmental processes, is thoroughly summarized. Furthermore, recent discoveries and achievements of this emerging omics tool in addressing inflammatory, traumatic, degenerative, and tumorous diseases of the musculoskeletal system are compiled. Finally, challenges and potential future directions for spatial transcriptomics, both as a field and in its applications in the musculoskeletal system, are discussed.

snRNA-Seq and Spatial Transcriptome Reveal Cell-Cell Crosstalk Mediated Metabolic Regulation in Porcine Skeletal Muscle

BACKGROUND

Cell-cell crosstalk between myogenic, adipogenic and immune cells in skeletal muscle to regulate energy metabolism and lipid deposition has received considerable attention. The specific mechanisms of interaction between the different cells in skeletal muscle are still unclear.

METHODS

Using integrated analysis of snRNA-seq and spatial transcriptome, the gene expression profile of longissimus dorsi (LD) muscle was compared between adult Taoyuan black (TB, obese, native Chinese breed) and Duroc (lean) pigs.

RESULTS

TB pig had more intramuscular fat (IMF) deposition (3.91%, p = 0.0244) and higher slow myofiber proportion (17.13%, p < 0.0001) compared with Duroc pig (IMF, 2.38%; slow myofiber, 6.92%) at the age of 180 days. We identified eight cell populations in porcine LD muscle. Five subpopulations of myonuclei and 10 subclusters of fibro/adipogenic progenitors (FAPs) were defined by marker genes. CellChat analysis revealed that communication between immune cells and other cells via the BMP and EGF signalling pathway was only observed in Duroc and not in TB pig. Both snRNA-seq and spatial transcriptome pointed out that FAPs are the important source of secretory proteins. A total of 35 upregulated and 23 downregulated differentially expressed genes (DEGs) were annotated as secretory, one upregulated and 36 downregulated secretory DEGs were identified between TB and Duroc pigs in FAPs by snRNA-seq and FAPs-high regions by spatial transcriptome, respectively. The distribution of FAPs was accompanied by the divergent myofiber-type composition. The expression level of slow myofiber marker gene (MYH7) was higher in both FAPs-high and FAPs-low regions of TB compared with Duroc pig (p < 0.0001), and expression level of fast myofiber maker gene (MYH1) was upregulated in FAPs-high region of Duroc compared with FAPs-high region of TB (p < 0.0001) and FAPs-low region of Duroc pig (p = 0.0002). The metabolic differences of myofibers between TB and Duroc pigs were mainly concentrated in energy, lipid and nitrogen metabolism-related pathway (p < 0.05). The significant correlation (R > 0.4, p < 0.05) between secretory and metabolism-related DEGs with spatial aggregation was verified by regression analysis for random region extraction (area of 25 spots, n = 400) from spatial transcriptome, and we speculated that the alteration of secretory proteins forming the microenvironment might regulate myofiber metabolism via target genes such as IRS1, PLPP1 and SLC38A2.

CONCLUSIONS

Our study provides new insights into skeletal muscle microenvironment that contributes to metabolic regulation and new methods and resources to study cell-cell communication in skeletal muscle.

Multi-task benchmarking of spatially resolved gene expression simulation models

BACKGROUND

Computational methods for spatially resolved transcriptomics (SRT) are often developed and assessed using simulated data. The effectiveness of these evaluations relies on the ability of simulation methods to accurately reflect experimental data. However, a systematic evaluation framework for spatial simulators is currently lacking.

RESULTS

Here, we present SpatialSimBench, a comprehensive evaluation framework that assesses 13 simulation methods using ten distinct STR datasets. We introduce simAdaptor, a tool that extends single-cell simulators by incorporating spatial variables, enabling them to simulate spatial data. SimAdaptor ensures SpatialSimBench is backwards compatible, facilitating direct comparisons between spatially aware simulators and existing non-spatial single-cell simulators through the adaption. Using SpatialSimBench, we demonstrate the feasibility of leveraging existing single-cell simulators for SRT data and highlight performance differences among methods. Additionally, we evaluate the simulation methods based on a total of 35 metrics across data property estimation, various downstream analyses, and scalability. In total, we generated 4550 results from 13 simulation methods, ten spatial datasets, and 35 metrics.

CONCLUSIONS

Our findings reveal that model estimation can be influenced by distribution assumptions and dataset characteristics. In summary, our evaluation framework provides guidelines for selecting appropriate methods for specific scenarios and informs future method development.

Inhibiting EZH2 complements steroid effects in Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a devastating X-linked disorder caused by dystrophin gene mutations. Despite recent advances in understanding the disease etiology and applying emerging treatment methodologies, glucocorticoid derivatives remain the only general therapeutic option that can slow disease development. However, the precise molecular mechanism of glucocorticoid action remains unclear, and there is still need for additional remedies to complement the treatment. Here, using single-nucleus RNA sequencing and spatial transcriptome analyses of human and mouse muscles, we investigated pathogenic features in patients with DMD and palliative effects of glucocorticoids. Our approach further illuminated the importance of proliferating satellite cells and revealed increased activity of a signal transduction pathway involving EZH2 in the patient cells. Subsequent administration of EZH2 inhibitors to Dmd mutant mice resulted in improved muscle phenotype through maintaining the immune-suppressing effect but overriding the muscle weakness and fibrogenic effects exerted by glucocorticoids. Our analysis reveals pathogenic mechanisms that can be readily targeted by extant therapeutic options for DMD.

Multi-modal refinement of the human heart atlas during the first gestational trimester

Forty first-trimester human hearts were studied to lay groundwork for further studies of the mechanisms underlying congenital heart defects. We first sampled 49,227 cardiac nuclei from three fetuses at 8.6, 9.0, and 10.7 post-conceptional weeks (pcw) for single-nucleus RNA sequencing, enabling the distinction of six classes comprising 21 cell types. Improved resolution led to the identification of previously unappreciated cardiomyocyte populations and minority autonomic and lymphatic endothelial transcriptomes, among others. After integration with 5-7 pcw heart single-cell RNA-sequencing data, we identified a human cardiomyofibroblast progenitor preceding the diversification of cardiomyocyte and stromal lineages. Spatial transcriptomic analysis (six Visium sections from two additional hearts) was aided by deconvolution, and key spatial markers validated on sectioned and whole hearts in two- and three-dimensional space and over time. Altogether, anatomical-positional features, including innervation, conduction and subdomains of the atrioventricular septum, translate latent molecular identity into specialized cardiac functions. This atlas adds unprecedented spatial and temporal resolution to the characterization of human-specific aspects of early heart formation.

Increased ectodysplasin-A2-receptor EDA2R is a ubiquitous hallmark of aging and mediates parainflammatory responses

Intensive efforts have been made to identify features that could serve as biomarkers of aging. Yet, drug-based interventions aimed at lessening the detrimental effects of getting older are lacking. This is largely attributable to tissue-specificity, sex-related differences, and to the difficulty of identifying actionable targets, which continues to pose a significant challenge. Here, we implement a bioinformatics approach revealing that aging-associated increase of the transmembrane Ectodysplasin-A2-Receptor is a prominent tissue-independent alteration occurring in humans and other species, and is particularly pronounced in models of accelerated aging. We show that strengthening of the Ectodysplasin-A2-Receptor signalling axis in myogenic precursors and differentiated myotubes suffices to trigger potent parainflammatory responses, mirroring aspects of aging-driven sarcopenia. Intriguingly, obesity, insulin-resistance, and aging-related comorbidities, such as type-2-diabetes, result in heightened levels of the Ectodysplasin-A2 ligand. Our findings suggest that targeting the Ectodysplasin-A2 surface receptor represents a promising pharmacological strategy to mitigate the development of aging-associated phenotypes.

A primer on global molecular responses to exercise in skeletal muscle: Omics in focus

Advances in skeletal muscle omics has expanded our understanding of exercise-induced adaptations at the molecular level. Over the past 2 decades, transcriptome studies in muscle have detailed acute and chronic responses to resistance, endurance, and concurrent exercise, focusing on variables such as training status, nutrition, age, sex, and metabolic health profile. Multi-omics approaches, such as the integration of transcriptomic and epigenetic data, along with emerging ribosomal RNA sequencing advancements, have further provided insights into how skeletal muscle adapts to exercise across the lifespan. Downstream of the transcriptome, proteomic and phosphoproteomic studies have identified novel regulators of exercise adaptations, while single-cell/nucleus and spatial sequencing technologies promise to evolve our understanding of cellular specialization and communication in and around skeletal muscle cells. This narrative review highlights (a) the historical foundations of exercise omics in skeletal muscle, (b) current research at 3 layers of the omics cascade (DNA, RNA, and protein), and (c) applications of single-cell omics and spatial sequencing technologies to study skeletal muscle adaptation to exercise. Further elaboration of muscle's global molecular footprint using multi-omics methods will help researchers and practitioners develop more effective and targeted approaches to improve skeletal muscle health as well as athletic performance.

Ginkgolide B increases healthspan and lifespan of female mice

Various anti-aging interventions show promise in extending lifespan, but many are ineffective or even harmful to healthspan. Ginkgolide B (GB), derived from Ginkgo biloba, reduces aging-related morbidities such as osteoporosis, yet its effects on healthspan and longevity have not been fully understood. In this study, we found that continuous oral administration of GB to female mice beginning at 20 months of age extended median survival and median lifespan by 30% and 8.5%, respectively. GB treatment also decreased tumor incidence; enhanced muscle quality, physical performance and metabolism; and reduced systemic inflammation and senescence. Single-nucleus RNA sequencing of skeletal muscle tissue showed that GB ameliorated aging-associated changes in cell type composition, signaling pathways and intercellular communication. GB reduced aging-induced Runx1+ type 2B myonuclei through the upregulation of miR-27b-3p, which suppresses Runx1 expression. Using functional analyses, we found that Runx1 promoted senescence and cell death in muscle cells. Collectively, these findings suggest the translational potential of GB to extend healthspan and lifespan and to promote healthy aging.

Sarcoglycans are enriched at the neuromuscular junction in a nerve-dependent manner

Sarcoglycanopathies are heterogeneous proximo-distal diseases presenting severe muscle alterations. Although there are 6 different sarcoglycan isoforms, sarcoglycanopathies are caused exclusively by mutations in genes coding for one of the four sarcoglycan transmembrane proteins (alpha, beta, gamma and delta) forming the sarcoglycan complex (SGC) in skeletal and cardiac muscle. Little is known about the different roles of the SGC beyond the dystrophin glycoprotein complex (DGC) structural role. Here, we show that SGC proteins are enriched at the post-synaptic membrane of neuromuscular junctions (NMJs). Using a mouse model lacking the beta-sarcoglycan subunit, we describe for the first time that the loss of the SGC in the NMJ area results in alterations of pre- and postsynaptic membrane, as well as a significant reduction of membrane potential. Moreover, using different denervated wild-type mouse models, we demonstrate that nerve presence precedes the sarcoglycan enrichment at NMJ, suggesting a nerve-dependent sarcoglycan expression. Altogether, our findings suggest that pathological decline should no longer be understood only in terms of sarcolemma damage but also in terms of sarcoglycans' participation in the NMJ. Henceforth, our work paves the way for the identification of new mechanisms involving sarcoglycans and new approaches for the treatment of sarcoglycanopathies.

A panoramic view of cell population dynamics in mammalian aging

To elucidate aging-associated cellular population dynamics, we present PanSci, a single-cell transcriptome atlas profiling >20 million cells from 623 mouse tissues across different life stages, sexes, and genotypes. This comprehensive dataset reveals >3000 different cellular states and >200 aging-associated cell populations. Our panoramic analysis uncovered organ-, lineage-, and sex-specific shifts in cellular dynamics during life-span progression. Moreover, we identify both systematic and organ-specific alterations in immune cell populations associated with aging. We further explored the regulatory roles of the immune system on aging and pinpointed specific age-related cell population expansions that are lymphocyte dependent. Our "cell-omics" strategy enhances comprehension of cellular aging and lays the groundwork for exploring the complex cellular regulatory networks in aging and aging-associated diseases.

Cross-species analyses of thymic mimetic cells reveal evolutionarily ancient origins and both conserved and species-specific elements

Thymic mimetic cells are molecular hybrids between medullary-thymic-epithelial cells (mTECs) and diverse peripheral cell types. They are involved in eliminating autoreactive T cells and can perform supplementary functions reflective of their peripheral-cell counterparts. Current knowledge about mimetic cells derives largely from mouse models. To provide the high resolution that proved revelatory for mice, we performed single-cell RNA sequencing on purified mimetic-cell compartments from human pediatric donors. The single-cell profiles of individual donors were surprisingly similar, with diversification of neuroendocrine subtypes and expansion of the muscle subtype relative to mice. Informatic and imaging studies on the muscle-mTEC population highlighted a maturation trajectory suggestive of skeletal-muscle differentiation, some striated structures, and occasional cellular groupings reminiscent of neuromuscular junctions. We also profiled thymic mimetic cells from zebrafish. Integration of data from the three species identified species-specific adaptations but substantial interspecies conservation, highlighting the evolutionarily ancient nature of mimetic mTECs. Our findings provide a landscape view of human mimetic cells, with anticipated relevance in autoimmunity.

Diacylglycerol kinase δ is required for skeletal muscle development and regeneration

Diacylglycerol kinase δ (DGKδ) phosphorylates diacylglycerol to produce phosphatidic acid. Previously, we demonstrated that down-regulation of DGKδ suppresses the myogenic differentiation of C2C12 myoblasts. However, the myogenic roles of DGKδ in vivo remain unclear. In the present study, we generated DGKδ-conditional knockout mice under the control of the myogenic factor 5 (Myf5) gene promoter, which regulates myogenesis and brown adipogenesis. The knockout mice showed a significant body weight reduction and apparent mass decrease in skeletal muscle, including the tibialis anterior (TA) muscle. Moreover, the thickness of a portion of the myofibers was reduced in DGKδ-deficient TA muscles. However, DGKδ deficiency did not substantially affect brown adipogenesis, suggesting that Myf5-driven DGKδ deficiency mainly affects muscle development. Notably, skeletal muscle injury induced by a cardiotoxin highly up-regulated DGKδ protein expression, and the DGKδ deficiency significantly reduced the thickness of myofibers, the expression levels of myogenic differentiation markers such as embryonic myosin heavy chain and myogenin, and the number of newly formed myofibers containing multiple central nuclei during muscle regeneration. DGKδ was strongly expressed in myogenin-positive satellite cells around the injured myofibers and centronucleated myofibers. These results indicate that DGKδ has important roles in muscle regeneration in activated satellite cells. Moreover, the conditional knockout mice fed with a high-fat diet showed increased fat mass and glucose intolerance. Taken together, these results demonstrate that DGKδ plays crucial roles in skeletal muscle development, regeneration, and function.

Methylome-proteome integration after late-life voluntary exercise training reveals regulation and target information for improved skeletal muscle health

Exercise is a potent stimulus for combatting skeletal muscle ageing. To study the effects of exercise on muscle in a preclinical setting, we developed a combined endurance-resistance training stimulus for mice called progressive weighted wheel running (PoWeR). PoWeR improves molecular, biochemical, cellular and functional characteristics of skeletal muscle and promotes aspects of partial epigenetic reprogramming when performed late in life (22-24 months of age). In this investigation, we leveraged pan-mammalian DNA methylome arrays and tandem mass-spectrometry proteomics in skeletal muscle to provide detailed information on late-life PoWeR adaptations in female mice relative to age-matched sedentary controls (n = 7-10 per group). Differential CpG methylation at conserved promoter sites was related to transcriptional regulation genes as well as Nr4a3, Hes1 and Hox genes after PoWeR. Using a holistic method of -omics integration called binding and expression target analysis (BETA), methylome changes were associated with upregulated proteins related to global and mitochondrial translation after PoWeR (P = 0.03). Specifically, BETA implicated methylation control of ribosomal, mitoribosomal, and mitochondrial complex I protein abundance after training. DNA methylation may also influence LACTB, MIB1 and UBR4 protein induction with exercise - all are mechanistically linked to muscle health. Computational cistrome analysis predicted several transcription factors including MYC as regulators of the exercise trained methylome-proteome landscape, corroborating prior late-life PoWeR transcriptome data. Correlating the proteome to muscle mass and fatigue resistance revealed positive relationships with VPS13A and NPL levels, respectively. Our findings expose differential epigenetic and proteomic adaptations associated with translational regulation after PoWeR that could influence skeletal muscle mass and function in aged mice. KEY POINTS: Late-life combined endurance-resistance exercise training from 22-24 months of age in mice is shown to improve molecular, biochemical, cellular and in vivo functional characteristics of skeletal muscle and promote aspects of partial epigenetic reprogramming and epigenetic age mitigation. Integration of DNA CpG 36k methylation arrays using conserved sites (which also contain methylation ageing clock sites) with exploratory proteomics in skeletal muscle extends our prior work and reveals coordinated and widespread regulation of ribosomal, translation initiation, mitochondrial ribosomal (mitoribosomal) and complex I proteins after combined voluntary exercise training in a sizeable cohort of female mice (n = 7-10 per group and analysis). Multi-omics integration predicted epigenetic regulation of serine β-lactamase-like protein (LACTB - linked to tumour resistance in muscle), mind bomb 1 (MIB1 - linked to satellite cell and type 2 fibre maintenance) and ubiquitin protein ligase E3 component N-recognin 4 (UBR4 - linked to muscle protein quality control) after training. Computational cistrome analysis identified MYC as a regulator of the late-life training proteome, in agreement with prior transcriptional analyses. Vacuolar protein sorting 13 homolog A (VPS13A) was positively correlated to muscle mass, and the glycoprotein/glycolipid associated sialylation enzyme N-acetylneuraminate pyruvate lyase (NPL) was associated to in vivo muscle fatigue resistance.

Exploring the therapeutic potential of fibroadipogenic progenitors in muscle disease

Skeletal muscle relies on its inherent self-repair ability to withstand continuous mechanical damage. Myofiber-intrinsic processes facilitate the repair of damage to sarcolemma and sarcomeres, but it is the coordinated interaction between muscle-resident satellite and stromal cells that are crucial in the regeneration of muscles to replace the lost muscle fibers. Fibroadipogenic progenitors (FAPs), are muscle-resident mesenchymal cells that are notable for their role in creating the dynamic stromal niche required to support long-term muscle homeostasis and regeneration. While FAP-mediated extracellular matrix formation and the establishment of a homeostatic muscle niche are essential for maintaining muscle health, excessive accumulation of FAPs and their aberrant differentiation leads to the fibrofatty degeneration that is a hallmark of myopathies and muscular dystrophies. Recent advancements, including single-cell RNA sequencing and in vivo analysis of FAPs, are providing deeper insights into the functions and specialization of FAPs, shedding light on their roles in both health and disease. This review will explore the above insights, discussing how FAP dysregulation contributes to muscle diseases. It will offer a concise overview of potential therapeutic interventions targeting FAPs to restore disrupted interactions among FAPs and muscle-resident cells, ultimately addressing degenerative muscle loss in neuromuscular diseases.

Single-nucleus transcriptomic profiling of the diaphragm during mechanical ventilation

Mechanical ventilation contributes to diaphragm atrophy and muscle weakness, which is referred to as ventilator-induced diaphragmatic dysfunction (VIDD). The pathogenesis of VIDD has not been fully understood until recently. The aim of this study was to investigate the effects of 24 h of mechanical ventilation on fibro-adipogenic progenitor (FAP) proliferation, endothelial-mesenchymal transition (EndMT), and immune cell infiltration driving diaphragm fibrosis in a rabbit model. The rabbits were anaesthetized and randomly divided into two groups (n = 3 each group): a control group and an experimental group. Diaphragm nuclei for sequencing were prepared by dissociating and filtering muscle tissue. 10X Genomics Platform for single-nucleus RNA sequencing (snRNA-seq) was used to profile the cells. Normalization and clustering were performed by Seurat, and clusters were manually annotated as different cell types. In this study, we performed differentially expressed genes (DEGs) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, pseudotime analysis and high dimensional weighted gene coexpression network analysis (hdWGCNA) to identify the key genes and signaling pathways related to the pathogenesis of VIDD. We further performed quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting to verify the results of snRNA-seq. The snRNA-seq results showed that acute postmechanical ventilation diaphragm cell changes included an increase in the proportion of fibroblasts and a decrease in the proportion of myofibres. The DEGs, KEGG, hdWGCNA and pseudotime analyses demonstrated that fibro-adipogenic progenitor (FAP) proliferation, endothelial-mesenchymal transition (EndMT) and immune cell infiltration are the three main processes involved in early stage of fibrosis development, among which Pdgfd, Sema3a, Cxcr2, are the corresponding regulatory genes. Glycolysis and the gene Pfkfb3 are also important metabolic factors for fibrosis formation. Negr1 and Mef2c are involved in phrenic nerve ending loss and diaphragm fibre atrophy. The qRT-PCR data showed that the mRNA levels of the genes Pdgfd, Cxcr2, Pfkfb3 and Negr1 were significantly greater in the experimental group than in the control group (P < 0.01), and the expression levels of Sema3a and Mef2c were significantly lower (P < 0.01). Despite limitations, including the lack of functional evaluations to confirm ventilator-induced diaphragm dysfunction (VIDD) and the absence of data validating diaphragm unloading during ventilation, our findings suggest that FAP proliferation and immune cell infiltration may play a role in the early stage of driving diaphragm fibrosis during mechanical ventilation. However, future studies are needed to confirm these findings and investigate the potential mechanisms underlying them.

Fibroblast diversification is an embryonic process dependent on muscle contraction

Fibroblasts, the most common cell type found in connective tissues, play major roles in development, homeostasis, regeneration, and disease. Although specific fibroblast subpopulations have been associated with different biological processes, the mechanisms and unique activities underlying their diversity have not been thoroughly examined. Here, we set out to dissect the variation in skeletal-muscle-resident fibroblasts (mrFibroblasts) during development. Our results demonstrate that mrFibroblasts diversify following the transition from embryonic to fetal myogenesis prior to birth. We find that mrFibroblasts segregate into two major subpopulations occupying distinct niches, with interstitial fibroblasts residing between the muscle fibers and delineating fibroblasts sheathing the muscle. We further show that these subpopulations entail distinct cellular dynamics and transcriptomes. Notably, we find that mrFibroblast subpopulations exert distinct regulatory roles on myoblast proliferation and differentiation. Finally, we demonstrate that this diversification depends on muscle contraction. Altogether, these findings establish that mrFibroblasts diversify in a spatiotemporal embryonic process into distinct cell types, entailing different characteristics and roles.

The expanding roles of myonuclei in adult skeletal muscle health and function

Skeletal muscle cells (myofibers) require multiple nuclei to support a cytoplasmic volume that is larger than other mononuclear cell types. It is dogmatic that mammalian resident myonuclei rely on stem cells (specifically satellite cells) for adding new DNA to muscle fibers to facilitate cytoplasmic expansion that occurs during muscle growth. In this review, we discuss the relationship between cell size and supporting genetic material. We present evidence that myonuclei may undergo DNA synthesis as a strategy to increase genetic material in myofibers independent from satellite cells. We then describe the details of our experiments that demonstrated that mammalian myonuclei can replicate DNA in vivo. Finally, we present our findings in the context of expanding knowledge about myonuclear heterogeneity, myonuclear mobility and shape. We also address why myonuclear replication is potentially important and provide future directions for remaining unknowns. Myonuclear DNA replication, coupled with new discoveries about myonuclear transcription, morphology, and behavior in response to stress, may provide opportunities to leverage previously unappreciated skeletal muscle biological processes for therapeutic targets that support muscle mass, function, and plasticity.

A statistical approach for systematic identification of transition cells from scRNA-seq data

Decoding cellular state transitions is crucial for understanding complex biological processes in development and disease. While recent advancements in single-cell RNA sequencing (scRNA-seq) offer insights into cellular trajectories, existing tools primarily study expressional rather than regulatory state shifts. We present CellTran, a statistical approach utilizing paired-gene expression correlations to detect transition cells from scRNA-seq data without explicitly resolving gene regulatory networks. Applying our approach to various contexts, including tissue regeneration, embryonic development, preinvasive lesions, and humoral responses post-vaccination, reveals transition cells and their distinct gene expression profiles. Our study sheds light on the underlying molecular mechanisms driving cellular state transitions, enhancing our ability to identify therapeutic targets for disease interventions.

The 24-hour molecular landscape after exercise in humans reveals MYC is sufficient for muscle growth

A detailed understanding of molecular responses to a hypertrophic stimulus in skeletal muscle leads to therapeutic advances aimed at promoting muscle mass. To decode the molecular factors regulating skeletal muscle mass, we utilized a 24-h time course of human muscle biopsies after a bout of resistance exercise. Our findings indicate: (1) the DNA methylome response at 30 min corresponds to upregulated genes at 3 h, (2) a burst of translation- and transcription-initiation factor-coding transcripts occurs between 3 and 8 h, (3) changes to global protein-coding gene expression peaks at 8 h, (4) ribosome-related genes dominate the mRNA landscape between 8 and 24 h, (5) methylation-regulated MYC is a highly influential transcription factor throughout recovery. To test whether MYC is sufficient for hypertrophy, we periodically pulse MYC in skeletal muscle over 4 weeks. Transient MYC increases muscle mass and fiber size in the soleus of adult mice. We present a temporally resolved resource for understanding molecular adaptations to resistance exercise in muscle ( http://data.myoanalytics.com ) and suggest that controlled MYC doses influence the exercise-related hypertrophic transcriptional landscape.

Transcriptomic analysis of skeletal muscle regeneration across mouse lifespan identifies altered stem cell states

In aging, skeletal muscle regeneration declines due to alterations in both myogenic and non-myogenic cells and their interactions. This regenerative dysfunction is not understood comprehensively or with high spatiotemporal resolution. We collected an integrated atlas of 273,923 single-cell transcriptomes and high-resolution spatial transcriptomic maps from muscles of young, old and geriatric mice (~5, 20 and 26 months old) at multiple time points following myotoxin injury. We identified eight immune cell types that displayed accelerated or delayed dynamics by age. We observed muscle stem cell states and trajectories specific to old and geriatric muscles and evaluated their association with senescence by scoring experimentally derived and curated gene signatures in both single-cell and spatial transcriptomic data. This revealed an elevation of senescent-like muscle stem cell subsets within injury zones uniquely in aged muscles. This Resource provides a holistic portrait of the altered cellular states underlying muscle regenerative decline across mouse lifespan.

SMN Deficiency Induces an Early Non-Atrophic Myopathy with Alterations in the Contractile and Excitatory Coupling Machinery of Skeletal Myofibers in the SMN∆7 Mouse Model of Spinal Muscular Atrophy

Spinal muscular atrophy (SMA) is caused by a deficiency of the ubiquitously expressed survival motor neuron (SMN) protein. The main pathological hallmark of SMA is the degeneration of lower motor neurons (MNs) with subsequent denervation and atrophy of skeletal muscle. However, increasing evidence indicates that low SMN levels not only are detrimental to the central nervous system (CNS) but also directly affect other peripheral tissues and organs, including skeletal muscle. To better understand the potential primary impact of SMN deficiency in muscle, we explored the cellular, ultrastructural, and molecular basis of SMA myopathy in the SMNΔ7 mouse model of severe SMA at an early postnatal period (P0-7) prior to muscle denervation and MN loss (preneurodegenerative [PND] stage). This period contrasts with the neurodegenerative (ND) stage (P8-14), in which MN loss and muscle atrophy occur. At the PND stage, we found that SMN∆7 mice displayed early signs of motor dysfunction with overt myofiber alterations in the absence of atrophy. We provide essential new ultrastructural data on focal and segmental lesions in the myofibrillar contractile apparatus. These lesions were observed in association with specific myonuclear domains and included abnormal accumulations of actin-thin myofilaments, sarcomere disruption, and the formation of minisarcomeres. The sarcoplasmic reticulum and triads also exhibited ultrastructural alterations, suggesting decoupling during the excitation-contraction process. Finally, changes in intermyofibrillar mitochondrial organization and dynamics, indicative of mitochondrial biogenesis overactivation, were also found. Overall, our results demonstrated that SMN deficiency induces early and MN loss-independent alterations in myofibers that essentially contribute to SMA myopathy. This strongly supports the growing body of evidence indicating the existence of intrinsic alterations in the skeletal muscle in SMA and further reinforces the relevance of this peripheral tissue as a key therapeutic target for the disease.

Spatial multi-omics in whole skeletal muscle reveals complex tissue architecture

Myofibers are large multinucleated cells that have long thought to have a rather simple organization. Single-nucleus transcriptomics, spatial transcriptomics and spatial metabolomics analysis have revealed distinct transcription profiles in myonuclei related to myofiber type. However, the use of local tissue collection or dissociation methods have obscured the spatial organization. To elucidate the full tissue architecture, we combine two spatial omics, RNA tomography and mass spectrometry imaging. This enables us to map the spatial transcriptomic, metabolomic and lipidomic organization of the whole murine tibialis anterior muscle. Our findings on heterogeneity in fiber type proportions are validated with multiplexed immunofluorescence staining in tibialis anterior, extensor digitorum longus and soleus. Our results demonstrate unexpectedly strong regionalization of gene expression, metabolic differences and variable myofiber type proportion along the proximal-distal axis. These new insights in whole-tissue level organization reconcile sometimes conflicting results coming from previous studies relying on local sampling methods.

Molecular aspects of the exercise response and training adaptation in skeletal muscle

Skeletal muscle plasticity enables an enormous potential to adapt to various internal and external stimuli and perturbations. Most notably, changes in contractile activity evoke a massive remodeling of biochemical, metabolic and force-generating properties. In recent years, a large number of signals, sensors, regulators and effectors have been implicated in these adaptive processes. Nevertheless, our understanding of the molecular underpinnings of training adaptation remains rudimentary. Specifically, the mechanisms that underlie signal integration, output coordination, functional redundancy and other complex traits of muscle adaptation are unknown. In fact, it is even unclear how stimulus-dependent specification is brought about in endurance or resistance exercise. In this review, we will provide an overview on the events that describe the acute perturbations in single endurance and resistance exercise bouts. Furthermore, we will provide insights into the molecular principles of long-term training adaptation. Finally, current gaps in knowledge will be identified, and strategies for a multi-omic and -cellular analyses of the molecular mechanisms of skeletal muscle plasticity that are engaged in individual, acute exercise bouts and chronic training adaptation discussed.

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.

SpaDiT: diffusion transformer for spatial gene expression prediction using scRNA-seq

The rapid development of spatially resolved transcriptomics (SRT) technologies has provided unprecedented opportunities for exploring the structure of specific organs or tissues. However, these techniques (such as image-based SRT) can achieve single-cell resolution, but can only capture the expression levels of tens to hundreds of genes. Such spatial transcriptomics (ST) data, carrying a large number of undetected genes, have limited its application value. To address the challenge, we develop SpaDiT, a deep learning framework for spatial reconstruction and gene expression prediction using scRNA-seq data. SpaDiT employs scRNA-seq data as an a priori condition and utilizes shared genes between ST and scRNA-seq data as latent representations to construct inputs, thereby facilitating the accurate prediction of gene expression in ST data. SpaDiT enhances the accuracy of spatial gene expression predictions over a variety of spatial transcriptomics datasets. We have demonstrated the effectiveness of SpaDiT by conducting extensive experiments on both seq-based and image-based ST data. We compared SpaDiT with eight highly effective baseline methods and found that our proposed method achieved an 8%-12% improvement in performance across multiple metrics. Source code and all datasets used in this paper are available at https://github.com/wenwenmin/SpaDiT and https://zenodo.org/records/12792074.

Spatiotemporal transcriptomic mapping of regenerative inflammation in skeletal muscle reveals a dynamic multilayered tissue architecture

Tissue regeneration is orchestrated by macrophages that clear damaged cells and promote regenerative inflammation. How macrophages spatially adapt and diversify their functions to support the architectural requirements of actively regenerating tissue remains unknown. In this study, we reconstructed the dynamic trajectories of myeloid cells isolated from acutely injured and early stage dystrophic muscles. We identified divergent subsets of monocytes/macrophages and DCs and validated markers (e.g., glycoprotein NMB [GPNMB]) and transcriptional regulators associated with defined functional states. In dystrophic muscle, specialized repair-associated subsets exhibited distinct macrophage diversity and reduced DC heterogeneity. Integrating spatial transcriptomics analyses with immunofluorescence uncovered the ordered distribution of subpopulations and multilayered regenerative inflammation zones (RIZs) where distinct macrophage subsets are organized in functional zones around damaged myofibers supporting all phases of regeneration. Importantly, intermittent glucocorticoid treatment disrupted the RIZs. Our findings suggest that macrophage subtypes mediated the development of the highly ordered architecture of regenerative tissues, unveiling the principles of the structured yet dynamic nature of regenerative inflammation supporting effective tissue repair.

A comprehensive single-cell RNA transcriptomic analysis identifies a unique SPP1+ macrophages subgroup in aging skeletal muscle

Senescence of skeletal muscle (SkM) has been a primary contributor to senior weakness and disability in recent years. The gradually declining SkM function associated with senescence has recently been connected to an imbalance between damage and repair. Macrophages (Mac) are involved in SkM aging, and different macrophage subgroups hold different biological functions. Through comprehensive single-cell transcriptomic analysis, we first compared the metabolic pathways and biological functions of different types of cells in young (Y) and old (O) mice SkM. Strikingly, the Mac population in mice SkM was also explored, and we identified a unique Mac subgroup in O SkM characterized by highly expressed SPP1 with strong senescence and adipogenesis features. Further work was carried out on the metabolic and biological processes for these Mac subgroups. Besides, we verified that the proportion of the SPP1+ Mac was increased significantly in the quadriceps tissues of O mice, and the senotherapeutic drug combination dasatinib + quercetin (D + Q) could dramatically reduce its proportion. Our study provides novel insight into the potential role of SPP1+ Mac in SkM, which may serve as a senotherapeutic target in SkM aging.

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.

Lipidomics and single-cell RNA sequencing reveal lipid and cell dynamics of porcine glycerol-injured skeletal muscle regeneration model

AIMS

Intramuscular fat (IMF) infiltration and extracellular matrix (ECM) deposition are characteristic features of muscle dysfunction, such as muscular dystrophy and severe muscle injuries. However, the underlying mechanisms of cellular origin, adipocyte formation and fibrosis in skeletal muscle are still unclear.

MAIN METHODS

Pigs were injected with 50 % glycerol (GLY) to induce skeletal muscle injury and regeneration. The acyl chain composition was analyzed by lipidomics, and the cell atlas and molecular signatures were revealed via single-cell RNA sequencing (scRNA-seq). Adipogenesis analysis was performed on fibroblast/fibro-adipogenic progenitors (FAPs) isolated from pigs.

KEY FINDINGS

The porcine GLY-injured skeletal muscle regeneration model was characterized by IMF infiltration and ECM deposition. Skeletal muscle stem cells (MuSCs) and FAP clusters were analyzed to explore the potential mechanisms of adipogenesis and fibrosis, and it was found that the TGF-β signaling pathway might be a key switch that regulates differentiation. Consistently, activation of the TGF-β signaling pathway increased SMAD2/3 phosphorylation and inhibited adipogenesis in FAPs, while inhibition of the TGF-β signaling pathway increased the expression of PPARγ and promoted adipogenesis.

SIGNIFICANCE

GLY-induced muscle injury and regeneration provides comprehensive insights for the development of therapies for human skeletal muscle dysfunction and fatty infiltration-related diseases in which the TGF-β/SMAD signaling pathway might play a primary regulatory role.

RhoA-mediated G12-G13 signaling maintains muscle stem cell quiescence and prevents stem cell loss

Multiple processes control quiescence of muscle stem cells (MuSCs), which is instrumental to guarantee long-term replenishment of the stem cell pool. Here, we describe that the G-proteins G12-G13 integrate signals from different G-protein-coupled receptors (GPCRs) to control MuSC quiescence via activation of RhoA. Comprehensive screening of GPCR ligands identified two MuSC-niche-derived factors, endothelin-3 (ET-3) and neurotensin (NT), which activate G12-G13 signaling in MuSCs. Stimulation with ET-3 or NT prevented MuSC activation, whereas pharmacological inhibition of ET-3 or NT attenuated MuSC quiescence. Inactivation of Gna12-Gna13 or Rhoa but not of Gnaq-Gna11 completely abrogated MuSC quiescence, which depleted the MuSC pool and was associated with accelerated sarcopenia during aging. Expression of constitutively active RhoA prevented exit from quiescence in Gna12-Gna13 mutant MuSCs, inhibiting cell cycle entry and differentiation via Rock and formins without affecting Rac1-dependent MuSC projections, a hallmark of quiescent MuSCs. The study uncovers a critical role of G12-G13 and RhoA signaling for active regulation of MuSC quiescence.

Thrombospondin-1 promotes fibro-adipogenic stromal expansion and contractile dysfunction of the diaphragm in obesity

Pulmonary disorders affect 40%-80% of individuals with obesity. Respiratory muscle dysfunction is linked to these conditions; however, its pathophysiology remains largely undefined. Mice subjected to diet-induced obesity (DIO) develop diaphragm muscle weakness. Increased intradiaphragmatic adiposity and extracellular matrix (ECM) content correlate with reductions in contractile force. Thrombospondin-1 (THBS1) is an obesity-associated matricellular protein linked with muscular damage in genetic myopathies. THBS1 induces proliferation of fibro-adipogenic progenitors (FAPs) - mesenchymal cells that differentiate into adipocytes and fibroblasts. We hypothesized that THBS1 drives FAP-mediated diaphragm remodeling and contractile dysfunction in DIO. We tested this by comparing the effects of dietary challenge on diaphragms of wild-type (WT) and Thbs1-knockout (Thbs1-/-) mice. Bulk and single-cell transcriptomics demonstrated DIO-induced stromal expansion in WT diaphragms. Diaphragm FAPs displayed upregulation of ECM and TGF-β-related expression signatures and augmentation of a Thy1-expressing subpopulation previously linked to type 2 diabetes. Despite similar weight gain, Thbs1-/- mice were protected from these transcriptomic changes and from obesity-induced increases in diaphragm adiposity and ECM deposition. Unlike WT controls, Thbs1-/- diaphragms maintained normal contractile force and motion after DIO challenge. THBS1 is therefore a necessary mediator of diaphragm stromal remodeling and contractile dysfunction in overnutrition and a potential therapeutic target in obesity-associated respiratory dysfunction.

Single-cell spatial transcriptomics reveals a dystrophic trajectory following a developmental bifurcation of myoblast cell fates in facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in "FSHD-hi" and "FSHD-lo" states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

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.

Integration mapping of cardiac fibroblast single-cell transcriptomes elucidates cellular principles of fibrosis in diverse pathologies

Single-cell technology has allowed researchers to probe tissue complexity and dynamics at unprecedented depth in health and disease. However, the generation of high-dimensionality single-cell atlases and virtual three-dimensional tissues requires integrated reference maps that harmonize disparate experimental designs, analytical pipelines, and taxonomies. Here, we present a comprehensive single-cell transcriptome integration map of cardiac fibrosis, which underpins pathophysiology in most cardiovascular diseases. Our findings reveal similarity between cardiac fibroblast (CF) identities and dynamics in ischemic versus pressure overload models of cardiomyopathy. We also describe timelines for commitment of activated CFs to proliferation and myofibrogenesis, profibrotic and antifibrotic polarization of myofibroblasts and matrifibrocytes, and CF conservation across mouse and human healthy and diseased hearts. These insights have the potential to inform knowledge-based therapies.

Spatial transcriptomics in embryonic mouse diaphragm muscle reveals regional gradients and subdomains of developmental gene expression

The murine embryonic diaphragm is a primary model for studying myogenesis and neuro-muscular synaptogenesis, both representing processes regulated by spatially organized genetic programs of myonuclei located in distinct myodomains. However, a spatial gene expression pattern of embryonic mouse diaphragm has not been reported. Here, we provide spatially resolved gene expression data for horizontally sectioned embryonic mouse diaphragms at embryonic days E14.5 and E18.5. These data reveal gene signatures for specific muscle regions with distinct maturity and fiber type composition, as well as for a central neuromuscular junction (NMJ) and a peripheral myotendinous junction (MTJ) compartment. Comparing spatial expression patterns of wild-type mice with those of transgenic mice lacking either the skeletal muscle calcium channel CaV1.1 or β-catenin, reveals curtailed muscle development and dysregulated expression of genes potentially involved in NMJ formation. Altogether, these datasets provide a powerful resource for further studies of muscle development and NMJ formation in the mouse.

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.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

A Panoramic View of Cell Population Dynamics in Mammalian Aging

To elucidate the aging-associated cellular population dynamics throughout the body, here we present PanSci, a single-cell transcriptome atlas profiling over 20 million cells from 623 mouse tissue samples, encompassing a range of organs across different life stages, sexes, and genotypes. This comprehensive dataset allowed us to identify more than 3,000 unique cellular states and catalog over 200 distinct aging-associated cell populations experiencing significant depletion or expansion. Our panoramic analysis uncovered temporally structured, organ- and lineage-specific shifts of cellular dynamics during lifespan progression. Moreover, we investigated aging-associated alterations in immune cell populations, revealing both widespread shifts and organ-specific changes. We further explored the regulatory roles of the immune system on aging and pinpointed specific age-related cell population expansions that are lymphocyte-dependent. The breadth and depth of our 'cell-omics' methodology not only enhance our comprehension of cellular aging but also lay the groundwork for exploring the complex regulatory networks among varied cell types in the context of aging and aging-associated diseases.

Muscle stem cells as immunomodulator during regeneration

The skeletal muscle is well known for its remarkable ability to regenerate after injuries. The regeneration is a complex and dynamic process that involves muscle stem cells (also called muscle satellite cells, MuSCs), fibro-adipogenic progenitors (FAPs), immune cells, and other muscle-resident cell populations. The MuSCs are the myogenic cell populaiton that contribute nuclei directly to the regenerated myofibers, while the other cell types collaboratively establish a microenvironment that facilitates myogenesis of MuSCs. The myogenic process includes activation, proliferation and differentiationof MuSCs, and subsequent fusion their descendent mononuclear myocytes into multinuclear myotubes. While the contributions of FAPs and immune cells to this microenvironment have been well studied, the influence of MuSCs on other cell types remains poorly understood. This review explores recent evidence supporting the potential role of MuSCs as immunomodulators during muscle regeneration, either through cytokine production or ligand-receptor interactions.

Advances and Challenges in Cell Biology for Cultured Meat

Cultured meat is an emerging biotechnology that aims to produce meat from animal cell culture, rather than from the raising and slaughtering of livestock, on environmental and animal welfare grounds. The detailed understanding and accurate manipulation of cell biology are critical to the design of cultured meat bioprocesses. Recent years have seen significant interest in this field, with numerous scientific and commercial breakthroughs. Nevertheless, these technologies remain at a nascent stage, and myriad challenges remain, spanning the entire bioprocess. From a cell biological perspective, these include the identification of suitable starting cell types, tuning of proliferation and differentiation conditions, and optimization of cell-biomaterial interactions to create nutritious, enticing foods. Here, we discuss the key advances and outstanding challenges in cultured meat, with a particular focus on cell biology, and argue that solving the remaining bottlenecks in a cost-effective, scalable fashion will require coordinated, concerted scientific efforts. Success will also require solutions to nonscientific challenges, including regulatory approval, consumer acceptance, and market feasibility. However, if these can be overcome, cultured meat technologies can revolutionize our approach to food.

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.

Decoding aging-dependent regenerative decline across tissues at single-cell resolution

Regeneration across tissues and organs exhibits significant variation throughout the body and undergoes a progressive decline with age. To decode the relationships between aging and regenerative capacity, we conducted a comprehensive single-cell transcriptome analysis of regeneration in eight tissues from young and aged mice. We employed diverse analytical models to study tissue regeneration and unveiled the intricate cellular and molecular mechanisms underlying the attenuated regenerative processes observed in aged tissues. Specifically, we identified compromised stem cell mobility and inadequate angiogenesis as prominent contributors to this age-associated decline in regenerative capacity. Moreover, we discovered a unique subset of Arg1+ macrophages that were activated in young tissues but suppressed in aged regenerating tissues, suggesting their important role in age-related immune response disparities during regeneration. This study provides a comprehensive single-cell resource for identifying potential targets for interventions aimed at enhancing regenerative outcomes in the aging population.

ApoE isoform does not influence skeletal muscle regeneration in adult mice

Introduction: Apolipoprotein E (ApoE) has been shown to be necessary for proper skeletal muscle regeneration. Consistent with this finding, single-cell RNA-sequencing analyses of skeletal muscle stem cells (MuSCs) revealed that Apoe is a top marker of quiescent MuSCs that is downregulated upon activation. The purpose of this study was to determine if muscle regeneration is altered in mice which harbor one of the three common human ApoE isoforms, referred to as ApoE2, E3 and E4. Methods: Histomorphometric analyses were employed to assess muscle regeneration in ApoE2, E3, and E4 mice after 14 days of recovery from barium chloride-induced muscle damage in vivo, and primary MuSCs were isolated to assess proliferation and differentiation of ApoE2, E3, and E4 MuSCs in vitro. Results: There was no difference in the basal skeletal muscle phenotype of ApoE isoforms as evaluated by section area, myofiber cross-sectional area (CSA), and myonuclear and MuSC abundance per fiber. Although there were no differences in fiber-type frequency in the soleus, Type IIa relative frequency was significantly lower in plantaris muscles of ApoE4 mice compared to ApoE3. Moreover, ApoE isoform did not influence muscle regeneration as assessed by fiber frequency, fiber CSA, and myonuclear and MuSC abundance. Finally, there were no differences in the proliferative capacity or myogenic differentiation potential of MuSCs between any ApoE isoform. Discussion: Collectively, these data indicate nominal effects of ApoE isoform on the ability of skeletal muscle to regenerate following injury or the in vitro MuSC phenotype.

Spatial transcriptomic interrogation of the murine bone marrow signaling landscape

Self-renewal and differentiation of skeletal stem and progenitor cells (SSPCs) are tightly regulated processes, with SSPC dysregulation leading to progressive bone disease. While the application of single-cell RNA sequencing (scRNAseq) to the bone field has led to major advancements in our understanding of SSPC heterogeneity, stem cells are tightly regulated by their neighboring cells which comprise the bone marrow niche. However, unbiased interrogation of these cells at the transcriptional level within their native niche environment has been challenging. Here, we combined spatial transcriptomics and scRNAseq using a predictive modeling pipeline derived from multiple deconvolution packages in adult mouse femurs to provide an endogenous, in vivo context of SSPCs within the niche. This combined approach localized SSPC subtypes to specific regions of the bone and identified cellular components and signaling networks utilized within the niche. Furthermore, the use of spatial transcriptomics allowed us to identify spatially restricted activation of metabolic and major morphogenetic signaling gradients derived from the vasculature and bone surfaces that establish microdomains within the marrow cavity. Overall, we demonstrate, for the first time, the feasibility of applying spatial transcriptomics to fully mineralized tissue and present a combined spatial and single-cell transcriptomic approach to define the cellular components of the stem cell niche, identify cell‒cell communication, and ultimately gain a comprehensive understanding of local and global SSPC regulatory networks within calcified tissue.

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.

Sox11 is enriched in myogenic progenitors but dispensable for development and regeneration of the skeletal muscle

Transcription factors (TFs) play key roles in regulating differentiation and function of stem cells, including muscle satellite cells (MuSCs), a resident stem cell population responsible for postnatal regeneration of the skeletal muscle. Sox11 belongs to the Sry-related HMG-box (SOX) family of TFs that play diverse roles in stem cell behavior and tissue specification. Analysis of single-cell RNA-sequencing (scRNA-seq) datasets identify a specific enrichment of Sox11 mRNA in differentiating but not quiescent MuSCs. Consistent with the scRNA-seq data, Sox11 levels increase during differentiation of murine primary myoblasts in vitro. scRNA-seq data comparing muscle regeneration in young and old mice further demonstrate that Sox11 expression is reduced in aged MuSCs. Age-related decline of Sox11 expression is associated with reduced chromatin contacts within the topologically associating domains. Unexpectedly, Myod1Cre-driven deletion of Sox11 in embryonic myoblasts has no effects on muscle development and growth, resulting in apparently healthy muscles that regenerate normally. Pax7CreER- or Rosa26CreER- driven (MuSC-specific or global) deletion of Sox11 in adult mice similarly has no effects on MuSC differentiation or muscle regeneration. These results identify Sox11 as a novel myogenic differentiation marker with reduced expression in quiescent and aged MuSCs, but the specific function of Sox11 in myogenesis remains to be elucidated.

Voyager: exploratory single-cell genomics data analysis with geospatial statistics

Exploratory spatial data analysis (ESDA) can be a powerful approach to understanding single-cell genomics datasets, but it is not yet part of standard data analysis workflows. In particular, geospatial analyses, which have been developed and refined for decades, have yet to be fully adapted and applied to spatial single-cell analysis. We introduce the Voyager platform, which systematically brings the geospatial ESDA tradition to (spatial) -omics, with local, bivariate, and multivariate spatial methods not yet commonly applied to spatial -omics, united by a uniform user interface. Using Voyager, we showcase biological insights that can be derived with its methods, such as biologically relevant negative spatial autocorrelation. Underlying Voyager is the SpatialFeatureExperiment data structure, which combines Simple Feature with SingleCellExperiment and AnnData to represent and operate on geometries bundled with gene expression data. Voyager has comprehensive tutorials demonstrating ESDA built on GitHub Actions to ensure reproducibility and scalability, using data from popular commercial technologies. Voyager is implemented in both R/Bioconductor and Python/PyPI, and features compatibility tests to ensure that both implementations return consistent results.

Thrombospondin-1 promotes fibro-adipogenic stromal expansion and contractile dysfunction of the diaphragm in obesity

Pulmonary disorders impact 40-80% of individuals with obesity. Respiratory muscle dysfunction is linked to these conditions; however, its pathophysiology remains largely undefined. Mice subjected to diet-induced obesity (DIO) develop diaphragmatic weakness. Increased intra-diaphragmatic adiposity and extracellular matrix (ECM) content correlate with reductions in contractile force. Thrombospondin-1 (THBS1) is an obesity-associated matricellular protein linked with muscular damage in genetic myopathies. THBS1 induces proliferation of fibro-adipogenic progenitors (FAPs)-mesenchymal cells that differentiate into adipocytes and fibroblasts. We hypothesized that THBS1 drives FAP-mediated diaphragm remodeling and contractile dysfunction in DIO. We tested this by comparing effects of dietary challenge on diaphragms of wild-type (WT) and Thbs1 knockout ( Thbs1 -/- ) mice. Bulk and single-cell transcriptomics demonstrated DIO-induced stromal expansion in WT diaphragms. Diaphragm FAPs displayed upregulation of ECM and TGFβ-related expression signatures, and augmentation of a Thy1 -expressing sub-population previously linked to type 2 diabetes. Despite similar weight gain, Thbs1 -/- mice were protected from these transcriptomic changes, and from obesity-induced increases in diaphragm adiposity and ECM deposition. Unlike WT controls, Thbs1 -/- diaphragms maintained normal contractile force and motion after DIO challenge. These findings establish THBS1 as a necessary mediator of diaphragm stromal remodeling and contractile dysfunction in overnutrition, and potential therapeutic target in obesity-associated respiratory dysfunction.