Gene expression profiling of Duchenne muscular dystrophy skeletal muscle

The Innovative Role of Nuclear Receptor Interaction Protein in Orchestrating Invadosome Formation for Myoblast Fusion

BACKGROUND

Nuclear receptor interaction protein (NRIP) is versatile and engages with various proteins to execute its diverse biological function. NRIP deficiency was reported to cause small myofibre size in adult muscle regeneration, indicating a crucial role of NRIP in myoblast fusion.

METHODS

The colocalization and interaction of NRIP with actin were investigated by immunofluorescence and immunoprecipitation assay, respectively. The participation of NRIP in myoblast fusion was demonstrated by cell fusion assay and time-lapse microscopy. The NRIP mutants were generated for mechanism study in NRIP-null C2C12 (termed KO19) cells and muscle-specific NRIP knockout (NRIP cKO) mice. A GEO profile database was used to analyse NRIP expression in Duchenne muscular dystrophy (DMD) patients.

RESULTS

In this study, we found that NRIP directly and reciprocally interacted with actin both in vitro and in cells. Immunofluorescence microscopy showed that the endogenous NRIP colocalized with components of invadosome, such as actin, Tks5, and cortactin, at the tips of cells during C2C12 differentiation. The KO19 cells were generated and exhibited a significant deficit in myoblast fusion compared with wild-type C2C12 cells (3.16% vs. 33.67%, p < 0.005). Overexpressed NRIP in KO19 cells could rescue myotube formation compared with control (3.37% vs. 1.00%, p < 0.01). We further confirmed that NRIP directly participated in cell fusion by using a cell-cell fusion assay. We investigated the mechanism of invadosome formation for myoblast fusion, which depends on NRIP-actin interaction, by analysing NRIP mutants in NRIP-null cells. Loss of actin-binding of NRIP reduced invadosome (enrichment ratio, 1.00 vs. 2.54, p < 0.01) and myotube formation (21.82% vs. 35.71%, p < 0.05) in KO19 cells and forced NRIP expression in KO19 cells and muscle-specific NRIP knockout (NRIP cKO) mice increased myofibre size compared with controls (over 1500 μm2, 61.01% vs. 20.57%, p < 0.001). We also found that the NRIP mRNA level was decreased in DMD patients compared with healthy controls (18 072 vs. 28 289, p < 0.001, N = 10 for both groups).

CONCLUSIONS

NRIP is a novel actin-binding protein for invadosome formation to induce myoblast fusion.

Oncostatin M signaling drives cancer-associated skeletal muscle wasting

Progressive weakness and muscle loss are associated with multiple chronic conditions, including muscular dystrophy and cancer. Cancer-associated cachexia, characterized by dramatic weight loss and fatigue, leads to reduced quality of life and poor survival. Inflammatory cytokines have been implicated in muscle atrophy; however, available anticytokine therapies failed to prevent muscle wasting in cancer patients. Here, we show that oncostatin M (OSM) is a potent inducer of muscle atrophy. OSM triggers cellular atrophy in primary myotubes using the JAK/STAT3 pathway. Identification of OSM targets by RNA sequencing reveals the induction of various muscle atrophy-related genes, including Atrogin1. OSM overexpression in mice causes muscle wasting, whereas muscle-specific deletion of the OSM receptor (OSMR) and the neutralization of circulating OSM preserves muscle mass and function in tumor-bearing mice. Our results indicate that activated OSM/OSMR signaling drives muscle atrophy, and the therapeutic targeting of this pathway may be useful in preventing muscle wasting.

A novel biomarker of fibrofatty replacement in dystrophinopathies identified by integrating transcriptome, magnetic resonance imaging, and pathology data

BACKGROUND

We aimed to analyse genome-wide transcriptome differences between Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) patients and identify biomarkers that correlate well with muscle magnetic resonance imaging (MRI) and histological fibrofatty replacement in both patients, which have not been reported.

METHODS

One hundred and one male patients with dystrophinopathies (55 DMD and 46 BMD) were enrolled. Muscle-derived genome-wide RNA-sequencing was performed in 31 DMD patients, 29 BMD patients, and 11 normal controls. Fibrofatty replacement was scored on muscle MRI and histological levels in all patients. A unique pipeline, single-sample gene set enrichment analysis combined with Spearman's rank correlations (ssGSEA-Cor) was developed to identify the most correlated gene signature for fibrofatty replacement. Quantitative real-time PCR (qRT-PCR) analysis, western blot analysis, and single-nucleus RNA-sequencing (snRNA-seq) were performed in the remaining patients to validate the most correlated gene signature.

RESULTS

Comparative transcriptomic analysis revealed that 31 DMD muscles were characterized by a significant increase of inflammation/immune response and extracellular matrix remodelling compared with 29 BMD muscles (P < 0.05). The ssGSEA-Cor pipeline revealed that the gene set of CDKN2A and CDKN2B was the most correlated gene signature for fibrofatty replacement (histological rs  = 0.744, P < 0.001; MRI rs  = 0.718, P < 0.001). Muscle qRT-PCR confirmed that CDKN2A mRNA expression in both 15 DMD (median = 25.007, P < 0.001) and 12 BMD (median = 5.654, P < 0.001) patients were significantly higher than that in controls (median = 1.101), while no significant difference in CDKN2B mRNA expression was found among DMD, BMD, and control groups. In the 27 patients, muscle CDKN2A mRNA expression respectively correlated with muscle MRI (rs  = 0.883, P < 0.001) and histological fibrofatty replacement (rs  = 0.804, P < 0.001) and disease duration (rs  = 0.645, P < 0.001) and North Star Ambulatory Assessment total scores (rs  = -0.698, P < 0.001). Muscle western blot analysis confirmed that both four DMD (median = 2.958, P < 0.05) and four BMD (median = 1.959, P < 0.01) patients had a significantly higher level of CDKN2A protein expression than controls (median = 1.068). The snRNA-seq analysis of two DMD muscles revealed that CDKN2A was mainly expressed in fibro-adipogenic progenitors, satellite cells, and myoblasts.

CONCLUSIONS

We identify CDKN2A expression as a novel biomarker of fibrofatty replacement, which might be a new target for antifibrotic therapy in dystrophinopathies.

Transcriptional dysregulation of autophagy in the muscle of a mouse model of Duchenne muscular dystrophy

It has been reported that autophagic activity is disturbed in the skeletal muscles of dystrophin-deficient mdx mice and patients with Duchenne muscular dystrophy (DMD). Transcriptional regulations of autophagy by FoxO transcription factors (FoxOs) and transcription factor EB (TFEB) play critical roles in adaptation to cellular stress conditions. Here, we investigated whether autophagic activity is dysregulated at the transcription level in dystrophin-deficient muscles. Expression levels of autophagy-related genes were globally decreased in tibialis anterior and soleus muscles of mdx mice compared with those of wild-type mice. DNA microarray data from the NCBI database also showed that genes related to autophagy were globally downregulated in muscles from patients with DMD. These downregulated genes are known as targets of FoxOs and TFEB. Immunostaining showed that nuclear localization of FoxO1 and FoxO3a was decreased in mdx mice. Western blot analyses demonstrated increases in phosphorylation levels of FoxO1 and FoxO3a in mdx mice. Nuclear localization of TFEB was also reduced in mdx mice, which was associated with elevated phosphorylation levels of TFEB. Collectively, the results suggest that autophagy is disturbed in dystrophin-deficient muscles via transcriptional downregulation due to phosphorylation-mediated suppression of FoxOs and TFEB.

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

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

Modeling Duchenne Muscular Dystrophy Cardiomyopathy with Patients' Induced Pluripotent Stem-Cell-Derived Cardiomyocytes

Duchenne muscular dystrophy (DMD) is an X-linked progressive muscle degenerative disease caused by mutations in the dystrophin gene, resulting in death by the end of the third decade of life at the latest. A key aspect of the DMD clinical phenotype is dilated cardiomyopathy, affecting virtually all patients by the end of the second decade of life. Furthermore, despite respiratory complications still being the leading cause of death, with advancements in medical care in recent years, cardiac involvement has become an increasing cause of mortality. Over the years, extensive research has been conducted using different DMD animal models, including the mdx mouse. While these models present certain important similarities to human DMD patients, they also have some differences which pose a challenge to researchers. The development of somatic cell reprograming technology has enabled generation of human induced pluripotent stem cells (hiPSCs) which can be differentiated into different cell types. This technology provides a potentially endless pool of human cells for research. Furthermore, hiPSCs can be generated from patients, thus providing patient-specific cells and enabling research tailored to different mutations. DMD cardiac involvement has been shown in animal models to include changes in gene expression of different proteins, abnormal cellular Ca²⁺ handling, and other aberrations. To gain a better understanding of the disease mechanisms, it is imperative to validate these findings in human cells. Furthermore, with the recent advancements in gene-editing technology, hiPSCs provide a valuable platform for research and development of new therapies including the possibility of regenerative medicine. In this article, we review the DMD cardiac-related research performed so far using human hiPSCs-derived cardiomyocytes (hiPSC-CMs) carrying DMD mutations.

EDA2R-NIK signalling promotes muscle atrophy linked to cancer cachexia

Skeletal muscle atrophy is a hallmark of the cachexia syndrome that is associated with poor survival and reduced quality of life in patients with cancer¹. Muscle atrophy involves excessive protein catabolism and loss of muscle mass and strength². An effective therapy against muscle wasting is currently lacking because mechanisms driving the atrophy process remain incompletely understood. Our gene expression analysis in muscle tissues indicated upregulation of ectodysplasin A2 receptor (EDA2R) in tumour-bearing mice and patients with cachectic cancer. Here we show that activation of EDA2R signalling promotes skeletal muscle atrophy. Stimulation of primary myotubes with the EDA2R ligand EDA-A2 triggered pronounced cellular atrophy by induction of the expression of muscle atrophy-related genes Atrogin1 and MuRF1. EDA-A2-driven myotube atrophy involved activation of the non-canonical NFĸB pathway and was dependent on NFκB-inducing kinase (NIK) activity. Whereas EDA-A2 overexpression promoted muscle wasting in mice, deletion of either EDA2R or muscle NIK protected tumour-bearing mice from loss of muscle mass and function. Tumour-induced oncostatin M (OSM) upregulated muscle EDA2R expression, and muscle-specific oncostatin M receptor (OSMR)-knockout mice were resistant to tumour-induced muscle wasting. Our results demonstrate that EDA2R-NIK signalling mediates cancer-associated muscle atrophy in an OSM-OSMR-dependent manner. Thus, therapeutic targeting of these pathways may be beneficial in prevention of muscle loss.

Downregulation of Dystrophin Expression Occurs across Diverse Tumors, Correlates with the Age of Onset, Staging and Reduced Survival of Patients

Altered dystrophin expression was found in some tumors and recent studies identified a developmental onset of Duchenne muscular dystrophy (DMD). Given that embryogenesis and carcinogenesis share many mechanisms, we analyzed a broad spectrum of tumors to establish whether dystrophin alteration evokes related outcomes. Transcriptomic, proteomic, and mutation datasets from fifty tumor tissues and matching controls (10,894 samples) and 140 corresponding tumor cell lines were analyzed. Interestingly, dystrophin transcripts and protein expression were found widespread across healthy tissues and at housekeeping gene levels. In 80% of tumors, DMD expression was reduced due to transcriptional downregulation and not somatic mutations. The full-length transcript encoding Dp427 was decreased in 68% of tumors, while Dp71 variants showed variability of expression. Notably, low expression of dystrophins was associated with a more advanced stage, older age of onset, and reduced survival across different tumors. Hierarchical clustering analysis of DMD transcripts distinguished malignant from control tissues. Transcriptomes of primary tumors and tumor cell lines with low DMD expression showed enrichment of specific pathways in the differentially expressed genes. Pathways consistently identified: ECM-receptor interaction, calcium signaling, and PI3K-Akt are also altered in DMD muscle. Therefore, the importance of this largest known gene extends beyond its roles identified in DMD, and certainly into oncology.

Application of Droplet Digital PCR Technology in Muscular Dystrophies Research

Although they are considered rare disorders, muscular dystrophies have a strong impact on people’s health. Increased disease severity with age, frequently accompanied by the loss of ability to walk in some people, and the lack of treatment, have directed the researchers towards the development of more effective therapeutic strategies aimed to improve the quality of life and life expectancy, slow down the progression, and delay the onset or convert a severe phenotype into a milder one. Improved understanding of the complex pathology of these diseases together with the tremendous advances in molecular biology technologies has led to personalized therapeutic procedures. Different approaches that are currently under extensive investigation require more efficient, sensitive, and less invasive methods. Due to its remarkable analytical sensitivity, droplet digital PCR has become a promising tool for accurate measurement of biomarkers that monitor disease progression and quantification of various therapeutic efficiency and can be considered a tool for non-invasive prenatal diagnosis and newborn screening. Here, we summarize the recent applications of droplet digital PCR in muscular dystrophy research and discuss the factors that should be considered to get the best performance with this technology.

Skeletal muscle phenotyping of Hippo gene-mutated mice reveals that Lats1 deletion increases the percentage of type I muscle fibers

The Hippo signal transduction network regulates transcription through Yap/Taz-Tead1-4 in many tissues including skeletal muscle. Whilst transgenic mice have been generated for many Hippo genes, the resultant skeletal muscle phenotypes were not always characterized. Here, we aimed to phenotype the hindlimb muscles of Hippo gene-mutated Lats1^−/−, Mst2^−/−, Vgll3^−/−, and Vgll4^+/− mice. This analysis revealed that Lats1^−/− mice have 11% more slow type I fibers than age and sex-matched wild-type controls. Moreover, the mRNA expression of slow Myh7 increased by 50%, and the concentration of type I myosin heavy chain is 80% higher in Lats1^−/− mice than in age and sex-matched wild-type controls. Second, to find out whether exercise-related stimuli affect Lats1, we stimulated C2C12 myotubes with the hypertrophy agent clenbuterol or the energy stress agent AICAR. We found that both stimulated Lats1 expression by 1.2 and 1.3 fold respectively. Third, we re-analyzed published datasets and found that Lats1 mRNA in muscle is 63% higher in muscular dystrophy, increases by 17–77% after cardiotoxin-induced muscle injury, by 41–71% in muscles during overload-induced hypertrophy, and by 19–21% after endurance exercise when compared to respective controls. To conclude, Lats1 contributes to the regulation of muscle fiber type proportions, and its expression is regulated by physiological and pathological situations in skeletal muscle.

Becker muscular dystrophy: case report, review of the literature, and analysis of differentially expressed hub genes

INTRODUCTION

Becker muscular dystrophy (BMD) is a genetic and progressive neuromuscular disease caused by mutations in the dystrophin gene with no available cure. A case report and comprehensive review of BMD cases aim to provide important clues for early diagnosis and implications for clinical practice. Genes and pathways identified from microarray data of muscle samples from patients with BMD help uncover the potential mechanism and provide novel therapeutic targets for dystrophin-deficient muscular dystrophies.

METHODS

We describe a BMD family with a 10-year-old boy as the proband and reviewed BMD cases from PubMed. Datasets from the Gene Expression Omnibus database were downloaded and integrated with the online software.

RESULTS

The systematic review revealed the clinical manifestations and mutation points of the dystrophin gene. Gene ontology analysis showed that extracellular matrix organization and extracellular structure organization with enrichment of upregulated genes coexist in three datasets. We present the first report of TUBA1A involvement in the development of BMD/Duchenne muscular dystrophy (DMD).

DISCUSSION

This study provides important implications for clinical practice, uncovering the potential mechanism of the progress of BMD/DMD, and provided new therapeutic targets.

Perturbation of PI3K/Akt signaling affected autophagy modulation in dystrophin-deficient myoblasts

BACKGROUND

The absence of dystrophin has gave a massive impact on myotube development in Muscular Dystrophy pathogenesis. One of the conserved signaling pathways involved in skeletal muscle differentiation is the PI3K/Akt/mTOR pathway that plays a vital role in autophagy regulation. To further understand and establish targeted therapy in dystrophin-deficient myoblasts, protein expression profiling has been determined which provides information on perturbed autophagy modulation and activation.

METHODS

In this study, a dystrophin-deficient myoblast cell line established from the skeletal muscle of a dystrophic (mdx) mouse was used as a model. The dfd13 (dystrophin-deficient) and C2C12 (non-dystrophic) myoblasts were cultured in low mitogen conditions for 10 days to induce differentiation. The cells were subjected to total protein extraction prior to Western blotting assay technique. Protein sub-fractionation has been conducted to determine protein localization. The live-cell analysis of autophagy assay was done using a flow cytometer.

RESULTS

In our culture system, the dfd13 myoblasts did not achieve terminal differentiation. PTEN expression was profoundly increased in dfd13 myoblasts throughout the differentiation day subsequently indicates perturbation of PI3K/Akt/mTOR regulation. In addition, rictor-mTORC2 was also found inactivated in this event. This occurrence has caused FoxO3 misregulation leads to higher activation of autophagy-related genes in dfd13 myoblasts. Autophagosome formation was increased as LC3B-I/II showed accumulation upon differentiation. However, the ratio of LC3B lipidation and autophagic flux were shown decreased which exhibited dystrophic features.

CONCLUSION

Perturbation of the PTEN-PI3K/Akt pathway triggers excessive autophagosome formation and subsequently reduced autophagic flux within dystrophin-deficient myoblasts where these findings are of importance to understand Duchenne Muscular Dystrophy (DMD) patients. We believe that some manipulation within its regulatory signaling reported in this study could help restore muscle homeostasis and attenuate disease progression. Video Abstract.

Muscle Enriched Lamin Interacting Protein (Mlip) Binds Chromatin and Is Required for Myoblast Differentiation

Muscle-enriched A-type lamin-interacting protein (Mlip) is a recently discovered Amniota gene that encodes proteins of unknown biological function. Here we report Mlip’s direct interaction with chromatin, and it may function as a transcriptional co-factor. Chromatin immunoprecipitations with microarray analysis demonstrated a propensity for Mlip to associate with genomic regions in close proximity to genes that control tissue-specific differentiation. Gel mobility shift assays confirmed that Mlip protein complexes with genomic DNA. Blocking Mlip expression in C2C12 myoblasts down-regulates myogenic regulatory factors (MyoD and MyoG) and subsequently significantly inhibits myogenic differentiation and the formation of myotubes. Collectively our data demonstrate that Mlip is required for C2C12 myoblast differentiation into myotubes. Mlip may exert this role as a transcriptional regulator of a myogenic program that is unique to amniotes.

Single Cell ADNP Predictive of Human Muscle Disorders: Mouse Knockdown Results in Muscle Wasting

Activity-dependent neuroprotective protein (ADNP) mutations are linked with cognitive dysfunctions characterizing the autistic-like ADNP syndrome patients, who also suffer from delayed motor maturation. We thus hypothesized that ADNP is deregulated in versatile myopathies and that local ADNP muscle deficiency results in myopathy, treatable by the ADNP fragment NAP. Here, single-cell transcriptomics identified ADNP as a major constituent of the developing human muscle. ADNP transcript concentrations further predicted multiple human muscle diseases, with concentrations negatively correlated with the ADNP target interacting protein, microtubule end protein 1 (EB1). Reverting back to modeling at the single-cell level of the male mouse transcriptome, Adnp mRNA concentrations age-dependently correlated with motor disease as well as with sexual maturation gene transcripts, while Adnp expressing limb muscle cells significantly decreased with aging. Mouse Adnp heterozygous deficiency exhibited muscle microtubule reduction and myosin light chain (Myl2) deregulation coupled with motor dysfunction. CRISPR knockdown of adult gastrocnemius muscle Adnp in a Cas9 mouse resulted in treadmill (male) and gait (female) dysfunctions that were specifically ameliorated by treatment with the ADNP snippet, microtubule interacting, Myl2-regulating, NAP (CP201). Taken together, our studies provide new hope for personalized diagnosis/therapeutics in versatile myopathies.

PTEN Inhibition Ameliorates Muscle Degeneration and Improves Muscle Function in a Mouse Model of Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is caused by a mutation of the muscle membrane protein dystrophin and characterized by severe degeneration of myofibers, progressive muscle wasting, loss of mobility, and, ultimately, cardiorespiratory failure and premature death. Currently there is no cure for DMD. Herein, we report that skeletal muscle-specific knockout (KO) of the phosphatase and tensin homolog (Pten) gene in an animal model of DMD (mdx mice) alleviates myofiber degeneration and restores muscle function without increasing tumor incidence. Specifically, Pten KO normalizes myofiber size and prevents muscular atrophy, and it improves grip strength and exercise performance in mdx mice. Pten KO also reduces fibrosis and inflammation, and it ameliorates muscle pathology in mdx mice. Unbiased RNA sequencing reveals that Pten KO upregulates extracellular matrix and basement membrane components positively correlated with wound healing and suppresses negative regulators of wound healing and lipid biosynthesis, thus improving the integrity of muscle basement membrane at the ultrastructural level. Importantly, pharmacological inhibition of PTEN similarly ameliorates muscle pathology and improves muscle integrity and function in mdx mice. Our findings provide evidence that PTEN inhibition may represent a potential therapeutic strategy to restore muscle function in DMD.

VGLL3 operates via TEAD1, TEAD3 and TEAD4 to influence myogenesis in skeletal muscle

VGLL proteins are transcriptional co-factors that bind TEAD family transcription factors to regulate events ranging from wing development in fly, to muscle fibre composition and immune function in mice. Here, we characterise Vgll3 in skeletal muscle. We found that mouse Vgll3 was expressed at low levels in healthy muscle but that its levels increased during hypertrophy or regeneration; in humans, VGLL3 was highly expressed in tissues from patients with various muscle diseases, such as in dystrophic muscle and alveolar rhabdomyosarcoma. Interaction proteomics revealed that VGLL3 bound TEAD1, TEAD3 and TEAD4 in myoblasts and/or myotubes. However, there was no interaction with proteins from major regulatory systems such as the Hippo kinase cascade, unlike what is found for the TEAD co-factors YAP (encoded by YAP1) and TAZ (encoded by WWTR1). Vgll3 overexpression reduced the activity of the Hippo negative-feedback loop, affecting expression of muscle-regulating genes including Myf5, Pitx2 and Pitx3, and genes encoding certain Wnts and IGFBPs. VGLL3 mainly repressed gene expression, regulating similar genes to those regulated by YAP and TAZ. siRNA-mediated Vgll3 knockdown suppressed myoblast proliferation, whereas Vgll3 overexpression strongly promoted myogenic differentiation. However, skeletal muscle was overtly normal in Vgll3-null mice, presumably due to feedback signalling and/or redundancy. This work identifies VGLL3 as a transcriptional co-factor operating with the Hippo signal transduction network to control myogenesis.

Summary: VGLL3 interacts with TEAD transcription factors to direct expression of crucial muscle regulatory genes and contribute to the control of skeletal myogenesis.

Tracking disease progression non-invasively in Duchenne and Becker muscular dystrophies

BACKGROUND

Analysis of muscle biopsies allowed to characterize the pathophysiological changes of Duchenne and Becker muscular dystrophies (D/BMD) leading to the clinical phenotype. Muscle tissue is often investigated during interventional dose finding studies to show in situ proof of concept and pharmacodynamics effect of the tested drug. Less invasive readouts are needed to objectively monitor patients' health status, muscle quality, and response to treatment. The identification of serum biomarkers correlating with clinical function and able to anticipate functional scales is particularly needed for personalized patient management and to support drug development programs.

METHODS

A large-scale proteomic approach was used to identify serum biomarkers describing pathophysiological changes (e.g. loss of muscle mass), association with clinical function, prediction of disease milestones, association with in vivo ³¹ P magnetic resonance spectroscopy data and dystrophin levels in muscles. Cross-sectional comparisons were performed to compare DMD patients, BMD patients, and healthy controls. A group of DMD patients was followed up for a median of 4.4 years to allow monitoring of individual disease trajectories based on yearly visits.

RESULTS

Cross-sectional comparison enabled to identify 10 proteins discriminating between healthy controls, DMD and BMD patients. Several proteins (285) were able to separate DMD from healthy, while 121 proteins differentiated between BMD and DMD; only 13 proteins separated BMD and healthy individuals. The concentration of specific proteins in serum was significantly associated with patients' performance (e.g. BMP6 serum levels and elbow flexion) or dystrophin levels (e.g. TIMP2) in BMD patients. Analysis of longitudinal trajectories allowed to identify 427 proteins affected over time indicating loss of muscle mass, replacement of muscle by adipose tissue, and cardiac involvement. Over-representation analysis of longitudinal data allowed to highlight proteins that could be used as pharmacodynamic biomarkers for drugs currently in clinical development.

CONCLUSIONS

Serum proteomic analysis allowed to not only discriminate among DMD, BMD, and healthy subjects, but it enabled to detect significant associations with clinical function, dystrophin levels, and disease progression.

Fractal dimension analysis reveals skeletal muscle disorganization in mdx mice

Duchenne Muscular Dystrophy (DMD) is characterized by muscle extracellular matrix disorganization due to the increased collagen deposition leading to fibrosis that significantly exacerbates disease progression. Fractal dimension analysis is a method that quantifies tissue/cellular disorganization and characterizes complex structures. The first objective of the present study was use fractal analysis to evaluate extracellular matrix disorganization in mdx mice soleus muscle. Next, we mimic a hyper-proliferation of fibrogenic cells by co-culturing NIH3T3 fibroblasts and C2C12 myoblasts to test whether fibroblasts induce disorganization in myoblast arrangement. Here, we show mdx presented high skeletal muscle disorganization as revealed by fractal analysis. Similarly, this method revealed that myoblasts co-cultured with fibroblast also presented cellular arrangement disorganization. We also reanalyzed skeletal muscle microarrays transcriptomic data from mdx and DMD patients that revealed transcripts related to extracellular matrix organization. This analysis also identified Osteoglycin, which was validated as a potential regulator of ECM organization in mdx dystrophic muscles. Our results demonstrate that fractal dimension is useful tool for the analysis of skeletal muscle disorganization in DMD and also reveal a fibroblast-myoblast cross-talk that contributes to "in vitro" myoblast disarrangement.

Cross-sectional serum metabolomic study of multiple forms of muscular dystrophy

Muscular dystrophies are characterized by a progressive loss of muscle tissue and/or muscle function. While metabolic alterations have been described in patients’‐derived muscle biopsies, non‐invasive readouts able to describe these alterations are needed in order to objectively monitor muscle condition and response to treatment targeting metabolic abnormalities. We used a metabolomic approach to study metabolites concentration in serum of patients affected by multiple forms of muscular dystrophy such as Duchenne and Becker muscular dystrophies, limb‐girdle muscular dystrophies type 2A and 2B, myotonic dystrophy type 1 and facioscapulohumeral muscular dystrophy. We show that 15 metabolites involved in energy production, amino acid metabolism, testosterone metabolism and response to treatment with glucocorticoids were differentially expressed between healthy controls and Duchenne patients. Five metabolites were also able to discriminate other forms of muscular dystrophy. In particular, creatinine and the creatine/creatinine ratio were significantly associated with Duchenne patients performance as assessed by the 6‐minute walk test and north star ambulatory assessment. The obtained results provide evidence that metabolomics analysis of serum samples can provide useful information regarding muscle condition and response to treatment, such as to glucocorticoids treatment.

Quantitative proteome profiling of dystrophic dog skeletal muscle reveals a stabilized muscular architecture and protection against oxidative stress after systemic delivery of MuStem cells

Proteomic profiling plays a decisive role in the elucidation of molecular signatures representative of a specific clinical context. MuStem cell based therapy represents a promising approach for clinical applications to cure Duchenne muscular dystrophy (DMD). To expand our previous studies collected in the clinically relevant DMD animal model, we decided to investigate the skeletal muscle proteome 4 months after systemic delivery of allogenic MuStem cells. Quantitative proteomics with isotope-coded protein labeling was used to compile quantitative changes in the protein expression profiles of muscle in transplanted Golden Retriever muscular dystrophy (GRMD) dogs as compared to Golden Retriever muscular dystrophy dogs. A total of 492 proteins were quantified, including 25 that were overrepresented and 46 that were underrepresented after MuStem cell transplantation. Interestingly, this study demonstrates that somatic stem cell therapy impacts on the structural integrity of the muscle fascicle by acting on fibers and its connections with the extracellular matrix. We also show that cell infusion promotes protective mechanisms against oxidative stress and favors the initial phase of muscle repair. This study allows us to identify putative candidates for tissue markers that might be of great value in objectively exploring the clinical benefits resulting from our cell-based therapy for DMD. All MS data have been deposited in the ProteomeXchange with identifier PXD001768 (http://proteomecentral.proteomexchange.org/dataset/PXD001768).

Inositol 1,4,5-trisphosphate (IP3)-dependent Ca2+ signaling mediates delayed myogenesis in Duchenne muscular dystrophy fetal muscle

Duchenne muscular dystrophy (DMD) is a progressive neuromuscular disorder characterized by muscle wasting and premature death. The defective gene is dystrophin, a structural protein, absence of which causes membrane fragility and myofiber necrosis. Several lines of evidence showed that in adult DMD patients dystrophin is involved in signaling pathways that regulate calcium homeostasis and differentiation programs. However, secondary aspects of the disease, such as inflammation and fibrosis development, might represent a bias in the analysis. Because fetal muscle is not influenced by gravity and does not suffer from mechanical load and/or inflammation, we investigated 12-week-old fetal DMD skeletal muscles, highlighting for the first time early alterations in signaling pathways mediated by the absence of dystrophin itself. We found that PLC/IP3/IP3R/Ryr1/Ca2+ signaling is widely active in fetal DMD skeletal muscles and, through the calcium-dependent PKCα protein, exerts a fundamental regulatory role in delaying myogenesis and in myofiber commitment. These data provide new insights into the origin of DMD pathology during muscle development.

Glucocorticoids enhance muscle endurance and ameliorate Duchenne muscular dystrophy through a defined metabolic program

Classic physiology studies dating to the 1930s demonstrate that moderate or transient glucocorticoid (GC) exposure improves muscle performance. The ergogenic properties of GCs are further evidenced by their surreptitious use as doping agents by endurance athletes and poorly understood efficacy in Duchenne muscular dystrophy (DMD), a genetic muscle-wasting disease. A defined molecular basis underlying these performance-enhancing properties of GCs in skeletal muscle remains obscure. Here, we demonstrate that ergogenic effects of GCs are mediated by direct induction of the metabolic transcription factor KLF15, defining a downstream pathway distinct from that resulting in GC-related muscle atrophy. Furthermore, we establish that KLF15 deficiency exacerbates dystrophic severity and muscle GC-KLF15 signaling mediates salutary therapeutic effects in the mdx mouse model of DMD. Thus, although glucocorticoid receptor (GR)-mediated transactivation is often associated with muscle atrophy and other adverse effects of pharmacologic GC administration, our data define a distinct GR-induced gene regulatory pathway that contributes to therapeutic effects of GCs in DMD through proergogenic metabolic programming.

Expression of nuclear factor of activated T cells (NFAT) and downstream muscle-specific proteins in ground squirrel skeletal and heart muscle during hibernation

The thirteen-lined ground squirrel (Ictidomys tridecemlineatus) undergoes remarkable adaptive changes during hibernation. Interestingly, skeletal muscle remodelling occurs during the torpor-arousal cycle of hibernation to prevent net muscle loss despite inactivity. Reversible cardiomyocyte hypertrophy occurs in cardiac muscle, allowing the heart to preserve cardiac output during hibernation, while avoiding chronic maladaptive hypertrophy post-hibernation. We propose that calcium signalling proteins [calcineurin (Cn), calmodulin (CaM), and calpain], the nuclear factor of activated T cell (NFAT) family of transcription factors, and the NFAT targets myoferlin and myomaker contribute significantly to adaptations taking place in skeletal and cardiac muscle during hibernation. Protein-level analyses were performed over several conditions: euthermic room temperature (ER), euthermic cold room (EC), entrance into (EN), early (ET), and late torpor (LT) time points, in addition to early (EA), interbout (IA), and late arousal (LA) time points using immunoblotting and DNA-protein interaction (DPI) enzyme-linked immunosorbent assay (ELISAs). In skeletal and cardiac muscle, NFATc2 protein levels were elevated during torpor. NFATc4 increased throughout the torpor-arousal cycle in both tissues, and NFATc1 showed this trend in cardiac muscle only. NFATc3 showed an elevation in DNA-binding activity but not expression during torpor. Myoferlin protein levels dramatically increased during torpor in both skeletal and cardiac muscle. Myomaker levels also increased significantly in cardiac muscle during torpor. Cardiac Cn levels remained stable, whereas CaM and calpain decreased throughout the torpor-arousal cycle. Activation and/or upregulation of NFATc2, c3, myoferlin, and myomaker at torpor could be part of a stress-response mechanism to preserve skeletal muscle mass, whereas CaM and calpain appear to initiate the rapid reversal of cardiac hypertrophy during arousal through downregulation of the NFAT-Cn pathway.

MicroRNA and mRNA Expression Changes in Steroid Naïve and Steroid Treated DMD Patients

Background:

Duchenne Muscular Dystrophy (DMD) is a recessive X-linked form of muscular dystrophy. Steroid therapy has clinical benefits for DMD patients, but the mechanism remains unclear.

Objective:

This study was designed to identify mRNAs and microRNAs regulated in Duchenne Muscular Dystrophy patients prior to and after steroid therapy.

Methods:

Genome wide transcriptome profiling of whole blood was performed to identify mRNAs and microRNAs regulated in DMD patients.

Results:

The data show many regulated mRNAs and some microRNAs, including some muscle-specific microRNAs (e.g., miR-206), that were significantly altered in blood of young (age 3–10) DMD patients compared to young controls. A total of 95 microRNAs, but no mRNAs, were differentially expressed in older DMD patients compared to matched controls (age 11–20). Steroid treatment reversed expression patterns of several microRNAs (miR-206, miR-181a, miR-4538, miR-4539, miR-606, and miR-454) that were altered in the young DMD patients. As an example, the over-expression of miR-206 in young DMD patients is predicted to down-regulate a set of target genes (e.g., RHGAP31, KHSRP, CORO1B, PTBP1, C7orf58, DLG4, and KLF4) that would worsen motor function. Since steroids decreased miR-206 expression to control levels, this could provide one mechanism by which steroids improve motor function.

Conclusions:

These identified microRNA-mRNA alterations will help better understand the pathophysiology of DMD and the response to steroid treatment.

Circulating Biomarkers for Duchenne Muscular Dystrophy

Duchenne muscular dystrophy is the most common form of muscular dystrophy. Genetic and biochemical research over the years has characterized the cause, pathophysiology and development of the disease providing several potential therapeutic targets and/or biomarkers. High throughput - omic technologies have provided a comprehensive understanding of the changes occurring in dystrophic muscles. Murine and canine animal models have been a valuable source to profile muscles and body fluids, thus providing candidate biomarkers that can be evaluated in patients. This review will illustrate known circulating biomarkers that could track disease progression and response to therapy in patients affected by Duchenne muscular dystrophy. We present an overview of the transcriptomic, proteomic, metabolomics and lipidomic biomarkers described in literature. We show how studies in muscle tissue have led to the identification of serum and urine biomarkers and we highlight the importance of evaluating biomarkers as possible surrogate endpoints to facilitate regulatory processes for new medicinal products.

Differential Gene Expression Profiling of Dystrophic Dog Muscle after MuStem Cell Transplantation

BACKGROUND

Several adult stem cell populations exhibit myogenic regenerative potential, thus representing attractive candidates for therapeutic approaches of neuromuscular diseases such as Duchenne Muscular Dystrophy (DMD). We have recently shown that systemic delivery of MuStem cells, skeletal muscle-resident stem cells isolated in healthy dog, generates the remodelling of muscle tissue and gives rise to striking clinical benefits in Golden Retriever Muscular Dystrophy (GRMD) dog. This global effect, which is observed in the clinically relevant DMD animal model, leads us to question here the molecular pathways that are impacted by MuStem cell transplantation. To address this issue, we compare the global gene expression profile between healthy, GRMD and MuStem cell treated GRMD dog muscle, four months after allogenic MuStem cell transplantation.

RESULTS

In the dystrophic context of the GRMD dog, disease-related deregulation is observed in the case of 282 genes related to various processes such as inflammatory response, regeneration, calcium ion binding, extracellular matrix organization, metabolism and apoptosis regulation. Importantly, we reveal the impact of MuStem cell transplantation on several molecular and cellular pathways based on a selection of 31 genes displaying signals specifically modulated by the treatment. Concomitant with a diffuse dystrophin expression, a histological remodelling and a stabilization of GRMD dog clinical status, we show that cell delivery is associated with an up-regulation of genes reflecting a sustained enhancement of muscle regeneration. We also identify a decreased mRNA expression of a set of genes having metabolic functions associated with lipid homeostasis and energy. Interestingly, ubiquitin-mediated protein degradation is highly enhanced in GRMD dog muscle after systemic delivery of MuStem cells.

CONCLUSIONS

Overall, our results provide the first high-throughput characterization of GRMD dog muscle and throw new light on the complex molecular/cellular effects associated with muscle repair and the clinical efficacy of MuStem cell-based therapy.

The Coiled-Coil Domain Containing 80 (ccdc80) gene regulates gadd45β2 expression in the developing somites of zebrafish as a new player of the hedgehog pathway

The Coiled-Coil Domain Containing 80 (CCDC80) gene has been identified as strongly induced in rat thyroid PC CL3 cells immortalized by the adenoviral E1A gene. In human, CCDC80 is a potential oncosoppressor due to its down-regulation in several tumor cell lines and tissues and it is expressed in almost all tissues. CCDC80 has homologous in mouse, chicken, and zebrafish. We cloned the zebrafish ccdc80 and analyzed its expression and function during embryonic development. The in-silico translated zebrafish protein shares high similarity with its mammalian homologous, with nuclear localization signals and a signal peptide. Gene expression analysis demonstrates that zebrafish ccdc80 is maternally and zygotically expressed throughout the development. In particular, ccdc80 is strongly expressed in the notochord and it is under the regulation of the Hedgehog pathway. In this work we investigated the functional effects of ccdc80-loss-of-function during embryonic development and verified its interaction with gadd45β2 in somitogenesis.

Pharmacologic management of Duchenne muscular dystrophy: target identification and preclinical trials

Duchenne muscular dystrophy (DMD) is an X-linked human disorder in which absence of the protein dystrophin causes degeneration of skeletal and cardiac muscle. For the sake of treatment development, over and above definitive genetic and cell-based therapies, there is considerable interest in drugs that target downstream disease mechanisms. Drug candidates have typically been chosen based on the nature of pathologic lesions and presumed underlying mechanisms and then tested in animal models. Mammalian dystrophinopathies have been characterized in mice (mdx mouse) and dogs (golden retriever muscular dystrophy [GRMD]). Despite promising results in the mdx mouse, some therapies have not shown efficacy in DMD. Although the GRMD model offers a higher hurdle for translation, dogs have primarily been used to test genetic and cellular therapies where there is greater risk. Failed translation of animal studies to DMD raises questions about the propriety of methods and models used to identify drug targets and test efficacy of pharmacologic intervention. The mdx mouse and GRMD dog are genetically homologous to DMD but not necessarily analogous. Subcellular species differences are undoubtedly magnified at the whole-body level in clinical trials. This problem is compounded by disparate cultures in clinical trials and preclinical studies, pointing to a need for greater rigor and transparency in animal experiments. Molecular assays such as mRNA arrays and genome-wide association studies allow identification of genetic drug targets more closely tied to disease pathogenesis. Genes in which polymorphisms have been directly linked to DMD disease progression, as with osteopontin, are particularly attractive targets.

Muscle lim protein isoform negatively regulates striated muscle actin dynamics and differentiation

Muscle lim protein (MLP) has emerged as a critical regulator of striated muscle physiology and pathophysiology. Mutations in cysteine and glycine-rich protein 3 (CSRP3), the gene encoding MLP, have been directly associated with human cardiomyopathies, whereas aberrant expression patterns are reported in human cardiac and skeletal muscle diseases. Increasing evidence suggests that MLP has an important role in both myogenic differentiation and myocyte cytoarchitecture, although the full spectrum of its intracellular roles has not been delineated. We report the discovery of an alternative splice variant of MLP, designated as MLP-b, showing distinct expression in neuromuscular disease and direct roles in actin dynamics and muscle differentiation. This novel isoform originates by alternative splicing of exons 3 and 4. At the protein level, it contains the N-terminus first half LIM domain of MLP and a unique sequence of 22 amino acids. Physiologically, it is expressed during early differentiation, whereas its overexpression reduces C2C12 differentiation and myotube formation. This may be mediated through its inhibition of MLP/cofilin-2-mediated F-actin dynamics. In differentiated striated muscles, MLP-b localizes to the sarcomeres and binds directly to Z-disc components, including α-actinin, T-cap and MLP. The findings of the present study unveil a novel player in muscle physiology and pathophysiology that is implicated in myogenesis as a negative regulator of myotube formation, as well as in differentiated striated muscles as a contributor to sarcomeric integrity.

Statistical insights into major human muscular diseases

Muscular diseases lead to muscle fiber degeneration, impairment of mobility, and in some cases premature death. Many of these muscular diseases are largely idiopathic. The goal of this study was to identify biomarkers based on their functional role and possible mechanisms of pathogenesis, specific to individual muscular disease. We analyzed the muscle transcriptome from five major muscular diseases: acute quadriplegic myopathy (AQM), amyotrophic lateral sclerosis (ALS), mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS), dermatomyositis (DM) and polymyositis (PM) using pairwise statistical comparison to identify uniquely regulated genes in each muscular disease. The genome-wide information encoded in the transcriptome provided biomarkers and functional insights into dysregulation in each muscular disease. The analysis showed that the dysregulation of genes in forward membrane pathway, responsible for transmitting action potential from neural excitation, is unique to AQM, while the dysregulation of myofibril genes, determinant of the mechanical properties of muscle, is unique to ALS, dysregulation of ER protein processing, responsible for correct protein folding, is unique to DM, and upregulation of immune response genes is unique to PM. We have identified biomarkers specific to each muscular disease which can be used for diagnostic purposes.

Comparative Genomics of X-linked Muscular Dystrophies: The Golden Retriever Model

Duchenne muscular dystrophy (DMD) is a devastating disease that dramatically decreases the lifespan and abilities of affected young people. The primary molecular cause of the disease is the absence of functional dystrophin protein, which is critical to proper muscle function. Those with DMD vary in disease presentation and dystrophin mutation; the same causal mutation may be associated with drastically different levels of disease severity. Also contributing to this variation are the influences of additional modifying genes and/or changes in functional elements governing such modifiers. This genetic heterogeneity complicates the efficacy of treatment methods and to date medical interventions are limited to treating symptoms. Animal models of DMD have been instrumental in teasing out the intricacies of DMD disease and hold great promise for advancing knowledge of its variable presentation and treatment. This review addresses the utility of comparative genomics in elucidating the complex background behind phenotypic variation in a canine model of DMD, Golden Retriever muscular dystrophy (GRMD). This knowledge can be exploited in the development of improved, more personalized treatments for DMD patients, such as therapies that can be tailor-matched to the disease course and genomic background of individual patients.

Gene expression profiling identifies molecular pathways associated with collagen VI deficiency and provides novel therapeutic targets

Ullrich congenital muscular dystrophy (UCMD), caused by collagen VI deficiency, is a common congenital muscular dystrophy. At present, the role of collagen VI in muscle and the mechanism of disease are not fully understood. To address this we have applied microarrays to analyse the transcriptome of UCMD muscle and compare it to healthy muscle and other muscular dystrophies. We identified 389 genes which are differentially regulated in UCMD relative to controls. In addition, there were 718 genes differentially expressed between UCMD and dystrophin deficient muscle. In contrast, only 29 genes were altered relative to other congenital muscular dystrophies. Changes in gene expression were confirmed by real-time PCR. The set of regulated genes was analysed by Gene Ontology, KEGG pathways and Ingenuity Pathway analysis to reveal the molecular functions and gene networks associated with collagen VI defects. The most significantly regulated pathways were those involved in muscle regeneration, extracellular matrix remodelling and inflammation. We characterised the immune response in UCMD biopsies as being mainly mediated via M2 macrophages and the complement pathway indicating that anti-inflammatory treatment may be beneficial to UCMD as for other dystrophies. We studied the immunolocalisation of ECM components and found that biglycan, a collagen VI interacting proteoglycan, was reduced in the basal lamina of UCMD patients. We propose that biglycan reduction is secondary to collagen VI loss and that it may be contributing towards UCMD pathophysiology. Consequently, strategies aimed at over-expressing biglycan and restore the link between the muscle cell surface and the extracellular matrix should be considered.

Combining docking site and phosphosite predictions to find new substrates: identification of smoothelin-like-2 (SMTNL2) as a c-Jun N-terminal kinase (JNK) substrate

Specific docking interactions between mitogen-activated protein kinases (MAPKs), their regulators, and their downstream substrates, are crucial for efficient and accurate signal transmission. To identify novel substrates of the c-Jun N-terminal kinase (JNK) family of MAPKs, we searched the human genome for proteins that contained (1), a predicted JNK-docking site (D-site); and (2), a cluster of putative JNK target phosphosites located close to the D-site. Here we describe a novel JNK substrate that emerged from this analysis, the functionally uncharacterized protein smoothelin-like 2 (SMTNL2). SMTNL2 protein bound with high-affinity to multiple MAPKs including JNK1-3 and ERK2; furthermore, the identity of conserved amino acids in the predicted docking site (residues 180-193) was necessary for this high-affinity binding. In addition, purified full-length SMTNL2 protein was phosphorylated by JNK1-3 in vitro, and this required the integrity of the D-site. Using mass spectrometry and mutagenesis, we identified four D-site-dependent phosphoacceptor sites in close proximity to the docking site, at S217, S241, T236 and T239. A short peptide comprised of the SMTNL2 D-site inhibited JNK-mediated phosphorylation of the ATF2 transcription factor, showing that SMTNL2 can compete with other substrates for JNK binding. Moreover, when transfected into HEK293 cells, SMTNL2 was phosphorylated by endogenous JNK in a D-site dependent manner, on the same residues identified in vitro. SMTNL2 protein was expressed in many mammalian tissues, with a notably high expression in skeletal muscle. Consistent with the hypothesis that SMTNL2 has a function in skeletal muscle, SMTNL2 protein expression was strongly induced during the transition from myoblasts to myotubes in differentiating C2C12 cells.

DMD transcript imbalance determines dystrophin levels

Duchenne and Becker muscular dystrophies are caused by out-of-frame and in-frame mutations, respectively, in the dystrophin encoding DMD gene. Molecular therapies targeting the precursor-mRNA are in clinical trials and show promising results. These approaches will depend on the stability and expression levels of dystrophin mRNA in skeletal muscles and heart. We report that the DMD gene is more highly expressed in heart than in skeletal muscles, in mice and humans. The transcript mutated in the mdx mouse model shows a 5' to 3' imbalance compared with that of its wild-type counterpart and reading frame restoration via antisense-mediated exon skipping does not correct this event. We also report significant transcript instability in 22 patients with Becker dystrophy, clarifying the fact that transcript imbalance is not caused by premature nonsense mutations. Finally, we demonstrate that transcript stability, rather than transcriptional rate, is an important determinant of dystrophin protein levels in patients with Becker dystrophy. We suggest that the availability of the complete transcript is a key factor to determine protein abundance and thus will influence the outcome of mRNA-targeting therapies.

Dystrophin-deficient pigs provide new insights into the hierarchy of physiological derangements of dystrophic muscle

Duchenne muscular dystrophy (DMD) is caused by mutations in the X-linked dystrophin (DMD) gene. The absence of dystrophin protein leads to progressive muscle weakness and wasting, disability and death. To establish a tailored large animal model of DMD, we deleted DMD exon 52 in male pig cells by gene targeting and generated offspring by nuclear transfer. DMD pigs exhibit absence of dystrophin in skeletal muscles, increased serum creatine kinase levels, progressive dystrophic changes of skeletal muscles, impaired mobility, muscle weakness and a maximum life span of 3 months due to respiratory impairment. Unlike human DMD patients, some DMD pigs die shortly after birth. To address the accelerated development of muscular dystrophy in DMD pigs when compared with human patients, we performed a genome-wide transcriptome study of biceps femoris muscle specimens from 2-day-old and 3-month-old DMD and age-matched wild-type pigs. The transcriptome changes in 3-month-old DMD pigs were in good concordance with gene expression profiles in human DMD, reflecting the processes of degeneration, regeneration, inflammation, fibrosis and impaired metabolic activity. In contrast, the transcriptome profile of 2-day-old DMD pigs showed similarities with transcriptome changes induced by acute exercise muscle injury. Our studies provide new insights into early changes associated with dystrophin deficiency in a clinically severe animal model of DMD.

The involvement of collagen triple helix repeat containing 1 in muscular dystrophies

Fibrosis is the main complication of muscular dystrophies. We identified collagen triple helix repeat containing 1 (Cthrc1) in skeletal and cardiac muscles of mice, representing Duchenne and congenital muscle dystrophies (DMD and CMD, respectively), and dysferlinopathy. In all of the mice, Cthrc1 was associated with high collagen type I levels; no Cthrc1 or collagen was observed in muscles of control mice. High levels of Cthrc1 were also observed in biopsy specimens from patients with DMD, in whom they were reversibly correlated with that of β-dystroglycan, whereas collagen type I levels were elevated in all patients with DMD. At the muscle sites where collagen and Cthrc1 were adjacent, collagen fibers appeared smaller, suggesting involvement of Cthrc1 in collagen turnover. Halofuginone, an inhibitor of Smad3 phosphorylation downstream of the transforming growth factor-β signaling, reduced Cthrc1 levels in skeletal and cardiac muscles of mice, representing DMD, CMD, and dysferlinopathy. The myofibroblasts infiltrating the dystrophic muscles of the murine models of DMD, CMD, and dysferlinopathy were the source of Cthrc1. Transforming growth factor-β did not affect Cthrc1 levels in the mdx fibroblasts but decreased them in the control fibroblasts, in association with increased migration of mdx fibroblasts and dystrophic muscle invasion by myofibroblasts. To our knowledge, this is the first demonstration of Cthrc1 as a marker of the severity of the disease progression in the dystrophic muscles, and as a possible target for therapy.

Functional classification of skeletal muscle networks. II. Applications to pathophysiology

In our preceding companion paper (Wang Y, Winters J, Subramaniam S. J Appl Physiol. doi: 10.1152/japplphysiol.01514.2011), we used extensive expression profile data on normal human subjects, in combination with legacy knowledge to classify skeletal muscle function into four models, namely excitation-activation, mechanical, metabolic, and signaling-production model families. In this paper, we demonstrate how this classification can be applied to study two well-characterized myopathies: amyotrophic lateral sclerosis (ALS) and Duchenne muscular dystrophy (DMD). Using skeletal muscle profile data from ALS and DMD patients compared with that from normal subjects, normal young in the case of DMD, we delineate molecular mechanisms that are causative and consequential to skeletal muscle dysfunction. In ALS, our analysis establishes the metabolic role and specifically identifies the mechanisms of calcium dysregulation and defects in mitochondrial transport of materials as important for muscle dysfunction. In DMD, we illustrate how impaired mechanical function is strongly coordinated with other three functional networks, resulting in transformation of the skeletal muscle into hybrid forms as a compensatory mechanism. Our functional models also provide, in exquisite detail, the mechanistic role of myriad proteins in these four families in normal and disease function.

C3KO mouse expression analysis: downregulation of the muscular dystrophy Ky protein and alterations in muscle aging

Mutations in CAPN3 gene cause limb-girdle muscular dystrophy type 2A (LGMD2A) characterized by muscle wasting and progressive degeneration of scapular and pelvic musculature. Since CAPN3 knockout mice (C3KO) display features of muscle pathology similar to those features observed in the earliest-stage or preclinical LGMD2A patients, gene expression profiling analysis in C3KO mice was performed to gain insight into mechanisms of disease. Two different comparisons were carried out in order to determine, first, the differential gene expression between wild-type (WT) and C3KO soleus and, second, to identify the transcripts differentially expressed in aging muscles of WT and C3KO mice. The up/downregulation of two genes, important for normal muscle function, was identified in C3KO mice: the Ky gene, encoding a protease implicated in muscle development, and Park2 gene encoding an E3 ubiquitin ligase (parkin). The Ky gene was downregulated in C3KO muscles suggesting that Ky protease may play a complementary role in regulating muscle cytoskeleton homeostasis in response to changes in muscle activity. Park2 was upregulated in the aged WT muscles but not in C3KO muscles. Taking into account the known functions of parkin E3 ligase, it is possible that it plays a role in ubiquitination and degradation of atrophy-specific and damaged proteins that are necessary to avoid cellular toxicity and a cellular stress response in aging muscles.

Transcriptional response to GAA deficiency (Pompe disease) in infantile-onset patients

Pompe disease is a genetic disorder resulting from a deficiency of lysosomal acid alpha-glucosidase (GAA) that manifests as a clinical spectrum with regard to symptom severity and rate of progression. In this study, we used microarrays to examine gene expression from the muscle of two cohorts of infantile-onset Pompe patients to identify transcriptional differences that may contribute to the disease phenotype. We found strong similarities among the gene expression profiles generated from biceps and quadriceps, and identified a number of signaling pathways altered in both cohorts. We also found that infantile-onset Pompe patient muscle had a gene expression pattern characteristic of immature or regenerating muscle, and exhibited many transcriptional markers of inflammation, despite having few overt signs of inflammatory infiltrate. Further, we identified genes exhibiting correlation between expression at baseline and response to therapy. This combined dataset can serve as a foundation for biological discovery and biomarker development to improve the treatment of Pompe disease.

Impaired calcium calmodulin kinase signaling and muscle adaptation response in the absence of calpain 3

Mutations in the non-lysosomal, cysteine protease calpain 3 (CAPN3) result in the disease limb girdle muscular dystrophy type 2A (LGMD2A). CAPN3 is localized to several subcellular compartments, including triads, where it plays a structural, rather than a proteolytic, role. In the absence of CAPN3, several triad components are reduced, including the major Ca(2+) release channel, ryanodine receptor (RyR). Furthermore, Ca(2+) release upon excitation is impaired in the absence of CAPN3. In the present study, we show that Ca-calmodulin protein kinase II (CaMKII) signaling is compromised in CAPN3 knockout (C3KO) mice. The CaMK pathway has been previously implicated in promoting the slow skeletal muscle phenotype. As expected, the decrease in CaMKII signaling that was observed in the absence of CAPN3 is associated with a reduction in the slow versus fast muscle fiber phenotype. We show that muscles of WT mice subjected to exercise training activate the CaMKII signaling pathway and increase expression of the slow form of myosin; however, muscles of C3KO mice do not exhibit these adaptive changes to exercise. These data strongly suggest that skeletal muscle's adaptive response to functional demand is compromised in the absence of CAPN3. In agreement with our mouse studies, RyR levels were also decreased in biopsies from LGMD2A patients. Moreover, we observed a preferential pathological involvement of slow fibers in LGMD2A biopsies. Thus, impaired CaMKII signaling and, as a result, a weakened muscle adaptation response identify a novel mechanism that may underlie LGMD2A and suggest a pharmacological target that should be explored for therapy.

Ccdc80-l1 Is involved in axon pathfinding of zebrafish motoneurons

Axon pathfinding is a subfield of neural development by which neurons send out axons to reach the correct targets. In particular, motoneurons extend their axons toward skeletal muscles, leading to spontaneous motor activity. In this study, we identified the zebrafish Ccdc80 and Ccdc80-like1 (Ccdc80-l1) proteins in silico on the basis of their high aminoacidic sequence identity with the human CCDC80 (Coiled-Coil Domain Containing 80). We focused on ccdc80-l1 gene that is expressed in nervous and non-nervous tissues, in particular in territories correlated with axonal migration, such as adaxial cells and muscle pioneers. Loss of ccdc80-l1 in zebrafish embryos induced motility issues, although somitogenesis and myogenesis were not impaired. Our results strongly suggest that ccdc80-l1 is involved in axon guidance of primary and secondary motoneurons populations, but not in their proper formation. ccdc80-l1 has a differential role as regards the development of ventral and dorsal motoneurons, and this is consistent with the asymmetric distribution of the transcript. The axonal migration defects observed in ccdc80-l1 loss-of-function embryos are similar to the phenotype of several mutants with altered Hedgehog activity. Indeed, we reported that ccdc80-l1 expression is positively regulated by the Hedgehog pathway in adaxial cells and muscle pioneers. These findings strongly indicate ccdc80-l1 as a down-stream effector of the Hedgehog pathway.

Canine models of Duchenne muscular dystrophy and their use in therapeutic strategies

Duchenne muscular dystrophy (DMD) is an X-linked recessive disorder in which the loss of dystrophin causes progressive degeneration of skeletal and cardiac muscle. Potential therapies that carry substantial risk, such as gene- and cell-based approaches, must first be tested in animal models, notably the mdx mouse and several dystrophin-deficient breeds of dogs, including golden retriever muscular dystrophy (GRMD). Affected dogs have a more severe phenotype, in keeping with that of DMD, so may better predict disease pathogenesis and treatment efficacy. Various phenotypic tests have been developed to characterize disease progression in the GRMD model. These biomarkers range from measures of strength and joint contractures to magnetic resonance imaging. Some of these tests are routinely used in clinical veterinary practice, while others require specialized equipment and expertise. By comparing serial measurements from treated and untreated groups, one can document improvement or delayed progression of disease. Potential treatments for DMD may be broadly categorized as molecular, cellular, or pharmacologic. The GRMD model has increasingly been used to assess efficacy of a range of these therapies. A number of these studies have provided largely general proof-of-concept for the treatment under study. Others have demonstrated efficacy using the biomarkers discussed. Importantly, just as symptoms in DMD vary among patients, GRMD dogs display remarkable phenotypic variation. Though confounding statistical analysis in preclinical trials, this variation offers insight regarding the role that modifier genes play in disease pathogenesis. By correlating functional and mRNA profiling results, gene targets for therapy development can be identified.

Regulation of DMD pathology by an ankyrin-encoded miRNA

Background

Duchenne muscular dystrophy (DMD) is an X-linked myopathy resulting from the production of a nonfunctional dystrophin protein. MicroRNA (miRNA) are small 21- to 24-nucleotide RNA that can regulate both individual genes and entire cell signaling pathways. Previously, we identified several mRNA, both muscle-enriched and inflammation-induced, that are dysregulated in the skeletal muscles of DMD patients. One particularly muscle-enriched miRNA, miR-486, is significantly downregulated in dystrophin-deficient mouse and human skeletal muscles. miR-486 is embedded within the ANKYRIN1(ANK1) gene locus, which is transcribed as either a long (erythroid-enriched) or a short (heart muscle- and skeletal muscle-enriched) isoform, depending on the cell and tissue types.

Results

Inhibition of miR-486 in normal muscle myoblasts results in inhibited migration and failure to repair a wound in primary myoblast cell cultures. Conversely, overexpression of miR-486 in primary myoblast cell cultures results in increased proliferation with no changes in cellular apoptosis. Using bioinformatics and miRNA reporter assays, we have identified platelet-derived growth factor receptor β, along with several other downstream targets of the phosphatase and tensin homolog deleted on chromosome 10/AKT (PTEN/AKT) pathway, as being modulated by miR-486. The generation of muscle-specific transgenic mice that overexpress miR-486 revealed that miR-486 alters the cell cycle kinetics of regenerated myofibers in vivo, as these mice had impaired muscle regeneration.

Conclusions

These studies demonstrate a link for miR-486 as a regulator of the PTEN/AKT pathway in dystrophin-deficient muscle and an important factor in the regulation of DMD muscle pathology.

STIM1 as a key regulator for Ca2+ homeostasis in skeletal-muscle development and function

Stromal interaction molecules (STIM) were identified as the endoplasmic-reticulum (ER) Ca2+ sensor controlling store-operated Ca2+ entry (SOCE) and Ca2+-release-activated Ca2+ (CRAC) channels in non-excitable cells. STIM proteins target Orai1-3, tetrameric Ca2+-permeable channels in the plasma membrane. Structure-function analysis revealed the molecular determinants and the key steps in the activation process of Orai by STIM. Recently, STIM1 was found to be expressed at high levels in skeletal muscle controlling muscle function and properties. Novel STIM targets besides Orai channels are emerging.Here, we will focus on the role of STIM1 in skeletal-muscle structure, development and function. The molecular mechanism underpinning skeletal-muscle physiology points toward an essential role for STIM1-controlled SOCE to drive Ca2+/calcineurin/nuclear factor of activated T cells (NFAT)-dependent morphogenetic remodeling programs and to support adequate sarcoplasmic-reticulum (SR) Ca2+-store filling. Also in our hands, STIM1 is transiently up-regulated during the initial phase of in vitro myogenesis of C2C12 cells. The molecular targets of STIM1 in these cells likely involve Orai channels and canonical transient receptor potential (TRPC) channels TRPC1 and TRPC3. The fast kinetics of SOCE activation in skeletal muscle seem to depend on the triad-junction formation, favoring a pre-localization and/or pre-formation of STIM1-protein complexes with the plasma-membrane Ca2+-influx channels. Moreover, Orai1-mediated Ca2+ influx seems to be essential for controlling the resting Ca2+ concentration and for proper SR Ca2+ filling. Hence, Ca2+ influx through STIM1-dependent activation of SOCE from the T-tubule system may recycle extracellular Ca2+ losses during muscle stimulation, thereby maintaining proper filling of the SR Ca2+ stores and muscle function. Importantly, mouse models for dystrophic pathologies, like Duchenne muscular dystrophy, point towards an enhanced Ca2+ influx through Orai1 and/or TRPC channels, leading to Ca2+-dependent apoptosis and muscle degeneration. In addition, human myopathies have been associated with dysfunctional SOCE. Immunodeficient patients harboring loss-of-function Orai1 mutations develop myopathies, while patients suffering from Duchenne muscular dystrophy display alterations in their Ca2+-handling proteins, including STIM proteins. In any case, the molecular determinants responsible for SOCE in human skeletal muscle and for dysregulated SOCE in patients of muscular dystrophy require further examination.

Alignment of gene expression profiles from test samples against a reference database: New method for context-specific interpretation of microarray data

Background

Gene expression microarray data have been organized and made available as public databases, but the utilization of such highly heterogeneous reference datasets in the interpretation of data from individual test samples is not as developed as e.g. in the field of nucleotide sequence comparisons. We have created a rapid and powerful approach for the alignment of microarray gene expression profiles (AGEP) from test samples with those contained in a large annotated public reference database and demonstrate here how this can facilitate interpretation of microarray data from individual samples.

Methods

AGEP is based on the calculation of kernel density distributions for the levels of expression of each gene in each reference tissue type and provides a quantitation of the similarity between the test sample and the reference tissue types as well as the identity of the typical and atypical genes in each comparison. As a reference database, we used 1654 samples from 44 normal tissues (extracted from the Genesapiens database).

Results

Using leave-one-out validation, AGEP correctly defined the tissue of origin for 1521 (93.6%) of all the 1654 samples in the original database. Independent validation of 195 external normal tissue samples resulted in 87% accuracy for the exact tissue type and 97% accuracy with related tissue types. AGEP analysis of 10 Duchenne muscular dystrophy (DMD) samples provided quantitative description of the key pathogenetic events, such as the extent of inflammation, in individual samples and pinpointed tissue-specific genes whose expression changed (SAMD4A) in DMD. AGEP analysis of microarray data from adipocytic differentiation of mesenchymal stem cells and from normal myeloid cell types and leukemias provided quantitative characterization of the transcriptomic changes during normal and abnormal cell differentiation.

Conclusions

The AGEP method is a widely applicable method for the rapid comprehensive interpretation of microarray data, as proven here by the definition of tissue- and disease-specific changes in gene expression as well as during cellular differentiation. The capability to quantitatively compare data from individual samples against a large-scale annotated reference database represents a widely applicable paradigm for the analysis of all types of high-throughput data. AGEP enables systematic and quantitative comparison of gene expression data from test samples against a comprehensive collection of different cell/tissue types previously studied by the entire research community.

Meta-analysis of muscle transcriptome data using the MADMuscle database reveals biologically relevant gene patterns

Background

DNA microarray technology has had a great impact on muscle research and microarray gene expression data has been widely used to identify gene signatures characteristic of the studied conditions. With the rapid accumulation of muscle microarray data, it is of great interest to understand how to compare and combine data across multiple studies. Meta-analysis of transcriptome data is a valuable method to achieve it. It enables to highlight conserved gene signatures between multiple independent studies. However, using it is made difficult by the diversity of the available data: different microarray platforms, different gene nomenclature, different species studied, etc.

Description

We have developed a system tool dedicated to muscle transcriptome data. This system comprises a collection of microarray data as well as a query tool. This latter allows the user to extract similar clusters of co-expressed genes from the database, using an input gene list. Common and relevant gene signatures can thus be searched more easily. The dedicated database consists in a large compendium of public data (more than 500 data sets) related to muscle (skeletal and heart). These studies included seven different animal species from invertebrates (Drosophila melanogaster, Caenorhabditis elegans) and vertebrates (Homo sapiens, Mus musculus, Rattus norvegicus, Canis familiaris, Gallus gallus). After a renormalization step, clusters of co-expressed genes were identified in each dataset. The lists of co-expressed genes were annotated using a unified re-annotation procedure. These gene lists were compared to find significant overlaps between studies.

Conclusions

Applied to this large compendium of data sets, meta-analyses demonstrated that conserved patterns between species could be identified. Focusing on a specific pathology (Duchenne Muscular Dystrophy) we validated results across independent studies and revealed robust biomarkers and new pathways of interest. The meta-analyses performed with MADMuscle show the usefulness of this approach. Our method can be applied to all public transcriptome data.

Maged1, a new regulator of skeletal myogenic differentiation and muscle regeneration

BACKGROUND

In normal adult skeletal muscle, cell turnover is very slow. However, after an acute lesion or in chronic pathological conditions, such as primary myopathies, muscle stem cells, called satellite cells, are induced to proliferate, then withdraw definitively from the cell cycle and fuse to reconstitute functional myofibers.

RESULTS

We show that Maged1 is expressed at very low levels in normal adult muscle but is strongly induced after injury, during the early phase of myoblast differentiation. By comparing in vitro differentiation of myoblasts derived from wild-type or Maged1 knockout mice, we observed that Maged1 deficiency results in reduced levels of p21CIP1/WAF1, defective cell cycle exit and impaired myotube maturation. In vivo, this defect results in delayed regeneration of injured muscle.

CONCLUSIONS

These data demonstrate for the first time that Maged1 is an important factor required for proper skeletal myoblast differentiation and muscle healing.

Myoferlin regulation by NFAT in muscle injury, regeneration and repair

Ferlin proteins mediate membrane-fusion events in response to Ca(2+). Myoferlin, a member of the ferlin family, is required for normal muscle development, during which it mediates myoblast fusion. We isolated both damaged and intact myofibers from a mouse model of muscular dystrophy using laser-capture microdissection and found that the levels of myoferlin mRNA and protein were increased in damaged myofibers. To better define the components of the muscle-injury response, we identified a discreet 1543-bp fragment of the myoferlin promoter, containing multiple NFAT-binding sites, and found that this was sufficient to drive high-level myoferlin expression in cells and in vivo. This promoter recapitulated normal myoferlin expression in that it was downregulated in healthy myofibers and was upregulated in response to myofiber damage. Transgenic mice expressing GFP under the control of the myoferlin promoter were generated and GFP expression in this model was used to track muscle damage in vivo after muscle injury and in muscle disease. Myoferlin modulates the response to muscle injury through its activity in both myoblasts and mature myofibers.

Role of thrombospondin 1 in macrophage inflammation in dysferlin myopathy

Muscle inflammation can be a prominent feature in several muscular dystrophies. In dysferlin myopathy, it is mainly composed of macrophages. To understand the origin of inflammation in dysferlin-deficient muscle, we analyzed soluble factors involved in monocyte chemotaxis released by myoblasts and myotubes from control and dysferlinopathy patients using a transwell system. Dysferlin-deficient myotubes released more soluble factors involved in monocyte chemotaxis compared with controls (p < 0.001). Messenger RNA microarray analysis showed a 3.2-fold increase of thrombospondin 1 (TSP-1) expression in dysferlin-deficient myotubes. Retrotranscriptasepolymerase chain reaction analysis, ELISA, and immunohistochemistry confirmed these results. Dysferlin mRNA knockdown with short-interfering RNA in normal myogenic cells resulted in TSP-1 mRNA upregulation and increased chemotaxis. Furthermore, monocyte chemotaxis was decreased when TSP-1 was blocked by specific antibodies. In muscle biopsies from dysferlinopathy patients, TSP-1 expression was increased in muscle fibers but not in biopsies of patientswith other myopathies with inflammation; TSP-1 was seen in some macrophages in all samples analyzed. Taken together, the data demonstrate that dysferlin-deficient muscle upregulates TSP-1 in vivoand in vitro and indicate that endogenous chemotactic factors arecrucial to the sustained inflammatory process observed in dysferlinopathies.

Transcriptomic analysis of dystrophin RNAi knockdown reveals a central role for dystrophin in muscle differentiation and contractile apparatus organization

BACKGROUND

Duchenne muscular dystrophy (DMD) is a fatal muscle wasting disorder caused by mutations in the dystrophin gene. DMD has a complex and as yet incompletely defined molecular pathophysiology hindering development of effective ameliorative approaches. Transcriptomic studies so far conducted on dystrophic cells and tissues suffer from non-specific changes and background noise due to heterogeneous comparisons and secondary pathologies. A study design in which a perfectly matched control cell population is used as reference for transcriptomic studies will give a much more specific insight into the effects of dystrophin deficiency and DMD pathophysiology.

RESULTS

Using RNA interference (RNAi) to knock down dystrophin in myotubes from C57BL10 mice, we created a homogenous model to study the transcriptome of dystrophin-deficient myotubes. We noted significant differences in the global gene expression pattern between these myotubes and their matched control cultures. In particular, categorical analyses of the dysregulated genes demonstrated significant enrichment of molecules associated with the components of muscle cell contractile unit, ion channels, metabolic pathways and kinases. Additionally, some of the dysregulated genes could potentially explain conditions and endophenotypes associated with dystrophin deficiency, such as dysregulation of calcium homeostasis (Pvalb and Casq1), or cardiomyopathy (Obscurin, Tcap). In addition to be validated by qPCR, our data gains another level of validity by affirmatively reproducing several independent studies conducted previously at genes and/or protein levels in vivo and in vitro.

CONCLUSION

Our results suggest that in striated muscles, dystrophin is involved in orchestrating proper development and organization of myofibers as contractile units, depicting a novel pathophysiology for DMD where the absence of dystrophin results in maldeveloped myofibers prone to physical stress and damage. Therefore, it becomes apparent that any gene therapy approaches for DMD should target early stages in muscle development to attain a maximum clinical benefit. With a clear and specific definition of the transcriptome of dystrophin deficiency, manipulation of identified dysregulated molecules downstream of dystrophin may lead to novel ameliorative approaches for DMD.