Absence of dystrophin protein causes cardiac dysfunction in patients with Duchenne muscular dystrophy (DMD). Unlike boys with DMD, the common mouse model of DMD (B10-mdx) does not manifest cardiac deficits until late adulthood. This has limited our understanding of the mechanism and therapeutic approaches to target the pediatric onset of cardiac pathology in DMD. Here we show that the mdx mouse model on the DBA/2 J genetic background (D2-mdx) displays juvenile-onset cardiac degeneration. Molecular and histological analysis revealed that cardiac damage in this model is linked to increased leukocyte chemotactic signaling and an inability to resolve inflammation. These deficiencies result in chronic inflammation and fibrotic conversion of the extracellular matrix (ECM) in the juvenile D2-mdx heart. To address these pathologies, we tested the utility of pro-resolution therapy to clear chronic cardiac inflammation. Use of an N-formyl peptide receptor (FPR) agonist helped physiologically resolve inflammation and mitigate the downstream events that lead to fibrotic degeneration of cardiomyocytes, preventing juvenile onset cardiac muscle loss. These results establish the utility of D2-mdx model to study events associated with pediatric-onset cardiac damage and demonstrates pro-resolution therapy as an alternate to anti-inflammatory therapy for treating degenerative cardiac pathology that leads to cardiomyopathy in DMD.

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.

Biallelic pathogenic variants in the gene encoding nebulin (NEB) are a known cause of congenital myopathy. We present two individuals with congenital myopathy and compound heterozygous variants (NM_001271208.2: c.2079C>A; p.(Cys693Ter) and c.21522+3A>G ) in NEB. Transcriptomic sequencing on patient muscle revealed that the extended splice variant c.21522+3A>G causes exon 144 skipping. Nebulin isoforms containing exon 144 are known to be mutually exclusive with isoforms containing exon 143, and these isoforms are differentially expressed during development and in adult skeletal muscles. Patients MRIs were compared to the known pattern of relative abundance of these two isoforms in muscle. We propose that the pattern of muscle involvement in these patients better fits the distribution of exon 144-containing isoforms in muscle than with previously published MRI findings in NEB-related disease due to other variants. To our knowledge this is the first report hypothesizing disease pathogenesis through the alteration of isoform distributions in muscle.

The application of allele-specific gene editing tools can expand the therapeutic options for dominant genetic conditions, either via gene correction or via allelic gene inactivation in situations where haploinsufficiency is tolerated. Here, we used allele-targeted CRISPR/Cas9 guide RNAs (gRNAs) to introduce inactivating frameshifting indels at a single nucleotide variant in the COL6A1 gene (c.868G>A; G290R), a variant that acts as dominant negative and that is associated with a severe form of congenital muscular dystrophy. We expressed spCas9 along with allele-targeted gRNAs, without providing a repair template, in primary fibroblasts derived from four patients and one control subject. Amplicon deep-sequencing for two gRNAs tested showed that single nucleotide deletions accounted for the majority of indels introduced. While activity of the two gRNAs was greater at the G290R allele, both gRNAs were also active at the wild-type allele. To enhance allele-selectivity, we introduced deliberate additional mismatches to one gRNA. One of these optimized gRNAs showed minimal activity at the WT allele, while generating productive edits and improving collagen VI matrix in cultured patient fibroblasts. This study strengthens the potential of gene editing to treat dominant-negative disorders, but also underscores the challenges in achieving allele selectivity with gRNAs.

Fibro/adipogenic progenitors (FAPs) are skeletal muscle stromal cells that support regeneration of injured myofibers and their maintenance in healthy muscles. FAPs are related to mesenchymal stem cells (MSCs/MeSCs) found in other adult tissues, but there is poor understanding of the extent of similarity between these cells. Using single-cell RNA sequencing (scRNA-seq) datasets from multiple mouse tissues, we have performed comparative transcriptomic analysis. This identified remarkable transcriptional similarity between FAPs and MeSCs, confirmed the suitability of PDGFRα as a reporter for FAPs, and identified extracellular proteolysis as a new FAP function. Using PDGFRα as a cell surface marker, we isolated FAPs from healthy and dysferlinopathic mouse muscles and performed scRNA-seq analysis. This revealed decreased FAP-mediated Wnt signaling as a potential driver of FAP dysfunction in dysferlinopathic muscles. Analysis of FAPs in dysferlin- and dystrophin-deficient muscles identified a relationship between the nature of muscle pathology and alteration in FAP gene expression.

Biallelic pathogenic variants in the genes encoding the dolichol-phosphate mannose synthase subunits (DPM) which produce mannosyl donors for glycosylphosphatidylinositols, N-glycan and protein O- and C-mannosylation, are rare causes of congenital disorders of glycosylation. Pathogenic variants in DPM1 and DPM2 are associated with muscle-eye-brain (MEB) disease, whereas DPM3 variants have mostly been reported in patients with isolated muscle disease-dystroglycanopathy. Thus far, only one affected individual with compound heterozygous DPM3 variants presenting with myopathy, mild intellectual disability, seizures, and nonspecific white matter abnormalities (WMA) around the lateral ventricles has been described. Here we present five affected individuals from four unrelated families with global developmental delay/intellectual disability ranging from mild to severe, microcephaly, seizures, WMA, muscle weakness and variable cardiomyopathy. Exome sequencing of the probands revealed an ultra-rare homozygous pathogenic missense DPM3 variant NM_018973.4:c.221A>G, p.(Tyr74Cys) which segregated with the phenotype in all families. Haplotype analysis indicated that the variant arose independently in three families. Functional analysis did not reveal any alteration in the N-glycosylation pathway caused by the variant; however, this does not exclude its pathogenicity in the function of the DPM complex and related cellular pathways. This report provides supporting evidence that, besides DPM1 and DPM2, defects in DPM3 can also lead to a muscle and brain phenotype.

Collagen VI-related dystrophies (COL6-RDs) are disorders of the myomatrix that manifest with a spectrum of severity. Over 400 variants have been described in the COL6A1, COL6A2, and COL6A3 genes in association with COL6-RD, although many have unclear pathogenicity. The recent discovery of a recurrent deep intronic splice-inducing variant (COL6A1 c.930+189C>T) that is actionable with antisense therapy prompted us to systematically explore the impact of splice-modulating variants in the COL6 genes. We queried our internal database of 225 molecularly confirmed and/or suspected COL6-RD patients. We found 130 patients (57.8%) carrying variants affecting, or potentially affecting, splicing, totaling 88 unique variants (61 dominant, 27 recessive; 40 in COL6A1, 29 in COL6A2, and 19 in COL6A3). To further investigate the effect of these variants on splicing, we selected 60 patients (representing 61 unique variants) for RNA-sequencing of muscle biopsy specimen (n=3) or of cultured dermal fibroblasts (n=57). The most frequent mechanism for dominant variants was in-frame exon skipping. However, unexpectedly in select cases, variants were leaky, causing only moderate levels of exon skipping predicted to translate into lesser clinical impact. In addition to the COL6A1 c.930+189C>T variant, we found two variants leading to in-frame insertions (COL6A2 c.900+1G>A; COL6A3 c.6309+10C>G), suggesting that this mechanism might be more frequent than previously thought. A recurrent recessive variant was COL6A2 c.1970-9G>A (n=6) was also characterized and found to create a novel acceptor splice site and to insert seven nucleotides in the mRNA. This study helps to further clarify the pathogenicity of COL6 variants and resulting genotype-phenotype correlations and to identify variants that may be candidates for therapeutic splice modulation. As individualized RNA-targeting therapies are becoming more feasible, the appropriate recognition of "actionable" variants is critical. Collagen VI-related dystrophies (COL6-RDs) are disorders of the myomatrix that manifest with a spectrum of severity. Over 400 variants have been described in the COL6A1, COL6A2, and COL6A3 genes in association with COL6-RD, although many have unclear pathogenicity. The recent discovery of a recurrent deep intronic splice-inducing variant (COL6A1 c.930+189C>T) that is actionable with antisense therapy prompted us to systematically explore the impact of splice-modulating variants in the COL6 genes. We queried our internal database of 225 molecularly confirmed and/or suspected COL6-RD patients. We found 130 patients (57.8%) carrying variants affecting, or potentially affecting, splicing, totaling 88 unique variants (61 dominant, 27 recessive; 40 in COL6A1, 29 in COL6A2, and 19 in COL6A3). To further investigate the effect of these variants on splicing, we selected 60 patients (representing 61 unique variants) for RNA-sequencing of muscle biopsy specimen (n=3) or of cultured dermal fibroblasts (n=57). The most frequent mechanism for dominant variants was in-frame exon skipping. However, unexpectedly in select cases, variants were leaky, causing only moderate levels of exon skipping predicted to translate into lesser clinical impact. In addition to the COL6A1 c.930+189C>T variant, we found two variants leading to in-frame insertions (COL6A2 c.900+1G>A; COL6A3 c.6309+10C>G), suggesting that this mechanism might be more frequent than previously thought. A recurrent recessive variant was COL6A2 c.1970-9G>A (n=6) was also characterized and found to create a novel acceptor splice site and to insert seven nucleotides in the mRNA. This study helps to further clarify the pathogenicity of COL6 variants and resulting genotype-phenotype correlations and to identify variants that may be candidates for therapeutic splice modulation. As individualized RNA-targeting therapies are becoming more feasible, the appropriate recognition of "actionable" variants is critical.

Aside from myofibers, numerous mononucleated cells reside in the skeletal muscle. These include the mesenchymal cells called fibro-adipogenic progenitors (FAPs), that support muscle development and regeneration in adult muscles. Recent evidence shows that defects in FAP function contributes to chronic muscle diseases and targeting FAPs offers avenues for treating these diseases.

Age-related loss of muscle mass and strength is widely attributed to limitation in the capacity of muscle resident satellite cells to perform their myogenic function. This idea contains two notions that have not been comprehensively evaluated by experiment. First, it entails the idea that we damage and lose substantial amounts of muscle in the course of our normal daily activities. Second, it suggests that mechanisms of muscle repair are in some way exhausted, thus limiting muscle regeneration. A third potential option is that the aged environment becomes inimical to the conduct of muscle regeneration. In the present study, we used our established model of human muscle xenografting to test whether muscle samples taken from cadavers, of a range of ages, maintained their myogenic potential after being transplanted into immunodeficient mice. We find no measurable difference in regeneration across the range of ages investigated up to 78 years of age. Moreover, we report that satellite cells maintained their myogenic capacity even when muscles were grafted 11 days postmortem in our model. We conclude that the loss of muscle mass with increasing age is not attributable to any intrinsic loss of myogenicity and is most likely a reflection of progressive and detrimental changes in the muscle microenvironment such as to disfavor the myogenic function of these cells.

The urea cycle protects the central nervous system from ammonia toxicity by converting ammonia to urea. N-acetylglutamate synthase (NAGS) catalyzes formation of N-acetylglutamate, an essential allosteric activator of carbamylphosphate synthetase 1. Enzymatic activity of mammalian NAGS doubles in the presence of L-arginine, but the physiological significance of NAGS activation by L-arginine has been unknown. The NAGS knockout ( Nags −/− ) mouse is an animal model of inducible hyperammonemia, which develops hyperammonemia without N-carbamylglutamate and L-citrulline supplementation (NCG + Cit). We used adeno associated virus (AAV) based gene transfer to correct NAGS deficiency in the Nags −/− mice, established the dose of the vector needed to rescue Nags −/− mice from hyperammonemia and measured expression levels of Nags mRNA and NAGS protein in the livers of rescued animals. This methodology was used to investigate the effect of L-arginine on ureagenesis in vivo by treating Nags −/− mice with AAV vectors encoding either wild-type or E354A mutant mouse NAGS (mNAGS), which is not activated by L-arginine. The Nags −/− mice expressing E354A mNAGS were viable but had elevated plasma ammonia concentration despite similar levels of the E354A and wild-type mNAGS proteins. The corresponding mutation in human NAGS (NP_694551.1:p.E360D) that abolishes binding and activation by L-arginine was identified in a patient with NAGS deficiency. Our results show that NAGS deficiency can be rescued by gene therapy, and suggest that L-arginine binding to the NAGS enzyme is essential for normal ureagenesis.

To define the transcriptomic changes responsible for the histologic alterations in skeletal muscle and their progression in collagen VI-related muscular dystrophy (COL6-RD). COL6-RD patient muscle biopsies were stratified into three groups based on the overall level of pathologic severity considering degrees of fibrosis, muscle fiber atrophy, and fatty replacement of muscle tissue. Using microarray and RNA-Seq, we then performed global gene expression profiling on the same muscle biopsies and compared their transcriptome with age- and sex-matched controls. COL6-RD muscle biopsy transcriptomes as a group revealed prominent upregulation of muscle extracellular matrix component genes and the downregulation of skeletal muscle and mitochondrion-specific genes. Upregulation of the TGFβ pathway was the most conspicuous change across all biopsies and was fully evident even in the mildest/earliest histological group. There was no difference in the overall transcriptional signature between the different histologic groups but polyserial analysis identified relative changes along with COL6-RD histological severity. Overall, our study establishes the prominent dysregulation of extracellular matrix genes, TGFβ signaling, and its downstream cellular pathways at the transcriptomic level in COL6-RD muscle.

BET1 is required, together with its SNARE complex partners GOSR2, SEC22b, and Syntaxin-5 for fusion of endoplasmic reticulum-derived vesicles with the ER-Golgi intermediate compartment (ERGIC) and the cis-Golgi. Here, we report three individuals, from two families, with severe congenital muscular dystrophy (CMD) and biallelic variants in BET1 (P1 p.(Asp68His)/p.(Ala45Valfs*2); P2 and P3 homozygous p.(Ile51Ser)). Due to aberrant splicing and frameshifting, the variants in P1 result in low BET1 protein levels and impaired ER-to-Golgi transport. Since in silico modeling suggested that p.(Ile51Ser) interferes with binding to interaction partners other than SNARE complex subunits, we set off and identified novel BET1 interaction partners with low affinity for p.(Ile51Ser) BET1 protein compared to wild-type, among them ERGIC-53. The BET1/ERGIC-53 interaction was validated by endogenous co-immunoprecipitation with both proteins colocalizing to the ERGIC compartment. Mislocalization of ERGIC-53 was observed in P1 and P2's derived fibroblasts; while in the p.(Ile51Ser) P2 fibroblasts specifically, mutant BET1 was also mislocalized along with ERGIC-53. Thus, we establish BET1 as a novel CMD/epilepsy gene and confirm the emerging role of ER/Golgi SNAREs in CMD.

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon-junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium ( Nrap −5kb), large (nebulin - 22 kb) and very-large (titin - 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type specific differential expression of these novel transcripts. The improved transcript identification and quantification demonstrated by our approach removes previous impediments to studies aimed at quantitative differential expression of ultra-long transcripts.

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium (Nrap: 5 kb), large (Neb: 22 kb), and very large (Ttn: 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type-specific differential expression of these novel transcripts. The improved transcript identification and quantification shown by our approach removes previous impediments to studies aimed at quantitative differential expression of ultralong transcripts.

MSTO1 encodes a cytosolic mitochondrial fusion protein, misato homolog 1 or MSTO1. While the full genotype–phenotype spectrum remains to be explored, pathogenic variants in MSTO1 have recently been reported in a small number of patients presenting with a phenotype of cerebellar ataxia, congenital muscle involvement with histologic findings ranging from myopathic to dystrophic and pigmentary retinopathy. The proposed underlying pathogenic mechanism of MSTO1-related disease is suggestive of impaired mitochondrial fusion secondary to a loss of function of MSTO1. Disorders of mitochondrial fusion and fission have been shown to also lead to mitochondrial DNA (mtDNA) depletion, linking them to the mtDNA depletion syndromes, a clinically and genetically diverse class of mitochondrial diseases characterized by a reduction of cellular mtDNA content. However, the consequences of pathogenic variants in MSTO1 on mtDNA maintenance remain poorly understood. We present extensive phenotypic and genetic data from 12 independent families, including 15 new patients harbouring a broad array of bi-allelic MSTO1 pathogenic variants, and we provide functional characterization from seven MSTO1-related disease patient fibroblasts. Bi-allelic loss-of-function variants in MSTO1 manifest clinically with a remarkably consistent phenotype of childhood-onset muscular dystrophy, corticospinal tract dysfunction and early-onset non-progressive cerebellar atrophy. MSTO1 protein was not detectable in the cultured fibroblasts of all seven patients evaluated, suggesting that pathogenic variants result in a loss of protein expression and/or affect protein stability. Consistent with impaired mitochondrial fusion, mitochondrial networks in fibroblasts were found to be fragmented. Furthermore, all fibroblasts were found to have depletion of mtDNA ranging from 30 to 70% along with alterations to mtDNA nucleoids. Our data corroborate the role of MSTO1 as a mitochondrial fusion protein and highlight a previously unrecognized link to mtDNA regulation. As impaired mitochondrial fusion is a recognized cause of mtDNA depletion syndromes, this novel link to mtDNA depletion in patient fibroblasts suggests that MSTO1-deficiency should also be considered a mtDNA depletion syndrome. Thus, we provide mechanistic insight into the disease pathogenesis associated with MSTO1 mutations and further define the clinical spectrum and the natural history of MSTO1-related disease.

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan.