Gene expression profiling in the early phases of DMD: a constant molecular signature characterizes DMD muscle from early postnatal life throughout disease progression

Identifying hub genes and dysregulated pathways in Duchenne muscular dystrophy

PURPOSE

The aim of this study was to identify the hub genes and dysregulated pathways in the progression of duchenne muscular dystrophy (DMD) and to unveil detailedly the cellular and molecular mechanisms associated with DMD for developing efficacious treatments in the future.

MATERIAL AND METHODS

Three mRNA microarray datasets (GSE13608, GSE38417 and GSE109178) were downloaded from Gene Expression Omnibus (GEO). The differentially expressed genes (DEGs) between DMD and normal tissues were obtained via R package. Function enrichment analyses were implemented respectively using DAVID online database. The network analysis of protein-protein interaction network (PPI) was conducted using String. Cytoscape and String were used to analyse modules and screen hub genes. The expression of the identified hub genes was confirmed in mdx mice through using qRT-PCR.

RESULTS

In total, 519 DEGs were identified, consisting of 393 upregulated genes and 126 downregulated genes. The enriched functions and pathways of the DEGs mainly involve extracellular matrix organization, collagen fibril organization, interferon-gamma-mediated signaling pathway, muscle contraction, endoplasmic reticulum lumen, MHC class II receptor activity, phagosome, graft-versus-host disease, cardiomyocytes, calcium signaling pathway. Twelve hub genes were discovered and biological process analysis proved that these genes were mainly enriched cell cycle, cell division. The result of qRT-PCR suggested that increase in expression of CD44, ECT2, TYMS, MAGEL2, HLA-DMA, SERPINH1, TNNT2 was confirmed in mdx mice and the downregulation of ASB2 and LEPREL1 was also observed.

CONCLUSION

In conclusion, DEGs and hub genes identified in the current research help us probe the molecular mechanisms underlying the pathogenesis and progression of DMD, and provide candidate targets for diagnosis and treatment of DMD.

Spatiotemporal diversity in molecular and functional abnormalities in the mdx dystrophic brain

Duchenne muscular dystrophy (DMD) is characterized by progressive muscle degeneration and neuropsychiatric abnormalities. Loss of full-length dystrophins is both necessary and sufficient to initiate DMD. These isoforms are expressed in the hippocampus, cerebral cortex (Dp427c), and cerebellar Purkinje cells (Dp427p). However, our understanding of the consequences of their absence, which is crucial for developing targeted interventions, remains inadequate. We combined RNA sequencing with genome-scale metabolic modelling (GSMM), immunodetection, and mitochondrial assays to investigate dystrophic alterations in the brains of the mdx mouse model of DMD. The cerebra and cerebella were analysed separately to discern the roles of Dp427c and Dp427p, respectively. Investigating these regions at 10 days (10d) and 10 weeks (10w) followed the evolution of abnormalities from development to early adulthood. These time points also encompass periods before onset and during muscle inflammation, enabling assessment of the potential damage caused by inflammatory mediators crossing the dystrophic blood-brain barrier. For the first time, we demonstrated that transcriptomic and functional dystrophic alterations are unique to the cerebra and cerebella and vary substantially between 10d and 10w. The common anomalies involved altered numbers of retained introns and spliced exons across mdx transcripts, corresponding with alterations in the mRNA processing pathways. Abnormalities in the cerebra were significantly more pronounced in younger mice. The top enriched pathways included those related to metabolism, mRNA processing, and neuronal development. GSMM indicated dysregulation of glucose metabolism, which corresponded with GLUT1 protein downregulation. The cerebellar dystrophic transcriptome, while significantly altered, showed an opposite trajectory to that of the cerebra, with few changes identified at 10 days. These late defects are specific and indicate an impact on the functional maturation of the cerebella that occurs postnatally. Although no classical neuroinflammation markers or microglial activation were detected at 10 weeks, specific differences indicate that inflammation impacts DMD brains. Importantly, some dystrophic alterations occur late and may therefore be amenable to therapeutic intervention, offering potential avenues for mitigating DMD-related neuropsychiatric defects.

Circulatory CCL2 distinguishes Duchenne muscular dystrophy dogs

To establish a minimally invasive approach to studying body-wide muscle inflammation in the canine Duchenne muscular dystrophy (DMD) model, we evaluated 13 cytokines/chemokines in frozen sera from 90 affected (239 sera) and 73 normal (189 sera) dogs (0.00 to 45.2 months of age). Linear mixed-effects model analysis suggested that ten cytokines/chemokines were significantly elevated in affected dogs, including interleukin (IL)-2, IL-6, IL-7, IL-8, IL-10, IL-15, IL-18, C-C motif chemokine ligand 2 (CCL2), C-X-C motif chemokine ligand 1 (CXCL1) and granulocyte-macrophage colony-stimulating factor (GM-CSF). Further, cytokine/chemokine elevation coincided with the onset of muscle disease. Importantly, only CCL2 showed consistent changes at all ages, with the most pronounced increase occurring between 3 and 9 months. To study the effects of sample storage and type, we compared fresh versus frozen, and serum versus plasma, samples from the same dog. Similar readings were often obtained in fresh and frozen sera. Although plasma readings were significantly lower for many cytokines/chemokines, this did not compromise the robustness of CCL2 as a biomarker. Our study establishes a baseline for using circulatory cytokines/chemokines as biomarkers in canine DMD studies.

Determination of qPCR reference genes suitable for normalizing gene expression in a novel model of Duchenne muscular dystrophy, the D2-mdx mouse

Duchenne muscular dystrophy (DMD) is a X-linked neuromuscular disorder arising from mutations in the dystrophin gene, leading to a progressive muscle wasting and disability. Currently there is no universal therapy, and there is thus a strong interest in preclinical studies for finding novel treatments. The most widely used and characterized mouse model for DMD is the C57BL/10ScSn-Dmdmdx/J (BL10-mdx), but this model exhibits mild pathology and does not replicate key features of human disease. The D2.B10-Dmdmdx/J (D2-mdx) mouse is a more recent model which seems to better mimics the complex human DMD phenotype. However, the D2-mdx mouse remains less extensively characterised than its BL10-mdx counterpart. Quantitative PCR analysis of gene expression is an important tool to monitor disease progression and evaluate therapeutic efficacy, but measurements must be normalised to stably expressed reference genes, which should ideally be determined and validated empirically. We examined gene expression in the gastrocnemius (GC), diaphragm (DIA) and heart in the D2-mdx mouse, the BL10-mdx mouse, and appropriate strain-matched wild-type controls (D2-wt and BL10-wt), from 4 to 52 weeks of age, using a large panel of candidate references (ACTB, AP3D1, CSNK2A2, GAPDH, HPRT1, PAK1IP1, RPL13A, SDHA, and in the heart, also HTATSF1 and HMBS). Data was analyzed using GeNorm, Bestkeeper, deltaCt and Normfinder algorithms to identify stable references under multiple possible scenarios. We show that CSNK2A2, AP3D1 and ACTB represent strong universal reference genes in both GC and DIA, regardless of age, muscle type, strain and genotype, while HTATSF1 and SDHA are optimal for the heart. GAPDH, HPRT1 and RPL13A were conversely revealed to be poor references, showing tissue-, age- or disease-specific changes in expression. Our results illustrate the importance of determining appropriate reference genes for specific comparative scenarios, but also reconfirm that universal panels can nevertheless be identified for normalising gene expression studies in even complex pathological states.

Integrative analysis of whole genome bisulfite and transcriptome sequencing reveals the effect of sodium butyrate on DNA methylation in the differentiation of bovine skeletal muscle satellite cells

Butyric acid as a short-chain fatty acid (SCFA) is one of the key microbial metabolites of ruminants. Numerous studies indicate that butyrate is crucial in muscle growth and development, and plays an important molecular regulatory role mainly by inhibiting histone deacetylation. DNA methylation, a major epigenetic modification, is involved in cell differentiation. Butyrate, in addition to its role in acetylation modifications, can alter the DNA methylation status of cells. However, the impact of butyrate on the DNA methylation of bovine skeletal muscle satellite cells (SMSCs) remains unclear. In this study, we developed a differentiation model of SMSCs and employed RNA sequencing (RNA-seq) alongside whole genome bisulfite sequencing (WGBS) to explore the effects of butyrate treatment on DNA methylation status and its relationship with gene expression. Treatment of SMSCs with sodium butyrate (NaB) at 1.0 mM for 2 days significantly inhibited the expression of DNA methyltransferases (DNMT1, DNMT2, DNMT3A) at the mRNA and protein levels while promoting the expression of demethylases (TET1, TET2, TET3) at mRNA levels. WGBS identified 4292 differentially methylated regions (DMRs), comprising 2294 hypermethylated and 1998 hypomethylated regions. These DMRs were significantly enriched in the MAPK, cAMP, and Wnt signaling pathways, all of which are implicated in myogenesis and development. Combining RNA-seq and WGBS data revealed a total of 130 overlapping genes, including MDFIC, CREBBP, DMD, LTBP2 and KLF4. These genes are predominantly involved in regulating the FoxO, MAPK, PI3K-Akt, and Wnt signaling pathways. This study provides new insights into the effects of butyrate-mediated DNA methylation on SMSC development and enhances our understanding of butyrate as an epigenetic modifier beyond its role in acetylation.

Transcriptional changes of genes encoding sarcoplasmic reticulum calcium binding and up-taking proteins in normal and Duchenne muscular dystrophy dogs

BACKGROUND

Cytosolic calcium overload contributes to muscle degradation in Duchenne muscular dystrophy (DMD). The sarcoplasmic reticulum (SR) is the primary calcium storage organelle in muscle. The sarco-endoplasmic reticulum ATPase (SERCA) pumps cytosolic calcium to the SR during muscle relaxation. Calcium is kept in the SR by calcium-binding proteins.

METHODS

Given the importance of the canine DMD model in translational studies, we examined transcriptional changes of SERCA (SERCA1 and SERCA2a), SERCA regulators (phospholamban, sarcolipin, myoregulin, and dwarf open reading frame), and SR calcium-binding proteins (calreticulin, calsequestrin 1, calsequestrin 2, and sarcalumenin) in skeletal muscle (diaphragm and extensor carpi ulnaris) and heart (left ventricle) in normal and affected male dogs by droplet digital PCR before the onset (≤ 2-m-old), at the active stage (8 to 16-m-old), and at the terminal stage (30 to 50-m-old) of the disease. Since many of these proteins are expressed in a fiber type-specific manner, we also evaluated fiber type composition in skeletal muscle.

RESULTS

In affected dog skeletal muscle, SERCA and its regulators were down-regulated at the active stage, but calcium-binding proteins (except for calsequestrin 1) were upregulated at the terminal stage. Surprisingly, nominal differences were detected in the heart. We also examined whether there exists sex-biased expression in 8 to 16-m-old dogs. Multiple transcripts were significantly downregulated in the heart and extensor carpi ulnaris muscle of female dogs. In fiber type analysis, we found significantly more type I fiber in the diaphragm of 8 to 16-m-old affected dogs, and significantly more type II fibers in the extensor carpi ulnaris of 30 to 50-m-old affected dogs. However, no difference was detected between male and female dogs.

CONCLUSIONS

Our study adds new knowledge to the understanding of muscle calcium regulation in normal and dystrophic canines.

Late-Stage Skeletal Muscle Transcriptome in Duchenne muscular dystrophy shows a BMP4-Induced Molecular Signature

Duchenne muscular dystrophy (DMD) is a fatal X-linked recessive disease due to loss-of-function mutations in the DYSTROPHIN gene. DMD-related skeletal muscle wasting is typified by an aberrant immune response involving upregulation of TGFβ family of cytokines. We previously demonstrated that bone morphogenetic protein 4 (BMP4) is increased in DMD and BMP4 stimulation induces a 20-fold upregulation of Smad8 transcription. However, the role of BMP4 in severely affected DMD skeletal muscle is unknown. We hypothesized that transcriptomic signatures in severely affected human DMD skeletal muscle are driven by BMP4 signaling. Transcriptomes from skeletal muscle biopsies of late-stage DMD vs. non-DMD controls and C2C12 muscle cells with or without BMP4 stimulation were generated by RNA-Seq and analyzed for single transcript differential expression as well as by Ingenuity Pathway Analysis and weighted gene co-expression network analyses. A total of 2,328 and 5,291 transcripts in the human muscle and C2C12 muscle cells, respectively, were differentially expressed. We identified an overlapping molecular signature of 1,027 genes dysregulated in DMD muscle that were induced in BMP4-stimulated C2C12 muscle cells. Highly upregulated DMD transcripts that overlapped with BMP4-stimulated C2C12 muscle cells included ADAMTS3, HCAR2, SERPING1, SMAD8 , and UNC13C. The DMD transcriptome was characterized by dysregulation of pathways involving immune function, extracellular matrix remodeling, and metabolic/mitochondrial function. In summary, we define a late-stage DMD skeletal muscle transcriptome that substantially overlaps with the BMP4-induced molecular signature in C2C12 muscle cells. This supports BMP4 as a disease-driving regulator of transcriptomic changes in late-stage DMD skeletal muscle and expands our understanding of the evolution of dystrophic signaling pathways and their associated gene networks that could be explored for therapeutic development.

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.

Ryanodine receptor type 1 content decrease-induced endoplasmic reticulum stress is a hallmark of myopathies

BACKGROUND

Decreased ryanodine receptor type 1 (RyR1) protein levels are a well-described feature of recessive RYR1-related myopathies. The aim of the present study was twofold: (1) to determine whether RyR1 content is also decreased in other myopathies and (2) to investigate the mechanisms by which decreased RyR1 protein triggers muscular disorders.

METHODS

We used publicly available datasets, muscles from human inflammatory and mitochondrial myopathies, an inducible muscle-specific RYR1 recessive mouse model and RyR1 knockdown in C2C12 muscle cells to measure RyR1 content and endoplasmic reticulum (ER) stress markers. Proteomics, lipidomics, molecular biology and transmission electron microscopy approaches were used to decipher the alterations associated with the reduction of RyR1 protein levels.

RESULTS

RYR1 transcripts were reduced in muscle samples of patients suffering from necrotizing myopathy (P = 0.026), inclusion body myopathy (P = 0.003), polymyositis (P < 0.001) and juvenile dermatomyositis (P < 0.001) and in muscle samples of myotonic dystrophy type 2 (P < 0.001), presymptomatic (P < 0.001) and symptomatic (P < 0.001) Duchenne muscular dystrophy, Becker muscular dystrophy (P = 0.004) and limb-girdle muscular dystrophy type 2A (P = 0.004). RyR1 protein content was also significantly decreased in inflammatory myopathy (-75%, P < 0.001) and mitochondrial myopathy (-71%, P < 0.001) muscles. Proteomics data showed that depletion of RyR1 protein in C2C12 myoblasts leads to myotubes recapitulating the common molecular alterations observed in myopathies. Mechanistically, RyR1 protein depletion reduces ER-mitochondria contact length (-26%, P < 0.001), Ca2+ transfer to mitochondria (-48%, P = 0.002) and the mitophagy gene Parkinson protein 2 transcripts (P = 0.037) and induces mitochondrial accumulation (+99%, P = 0.005) and dysfunction (P < 0.001). This was associated to the accumulation of deleterious sphingolipid species. Our data showed increased levels of the ER stress marker chaperone-binding protein/glucose regulated protein 78, GRP78-Bip, in RyR1 knockdown myotubes (+45%, P = 0.046), in mouse RyR1 recessive muscles (+58%, P = 0.001) and in human inflammatory (+96%, P = 0.006) and mitochondrial (+64%, P = 0.049) myopathy muscles. This was accompanied by increased protein levels of the pro-apoptotic protein CCAAT-enhancer-binding protein homologous protein, CHOP-DDIT3, in RyR1 knockdown myotubes (+27%, P < 0.001), mouse RyR1 recessive muscles (+63%, P = 0.009), human inflammatory (+50%, P = 0.038) and mitochondrial (+51%, P = 0.035) myopathy muscles. In publicly available datasets, the decrease in RYR1 content in myopathies was also associated to increased ER stress markers and RYR1 transcript levels are inversely correlated with ER stress markers in the control population.

CONCLUSIONS

Decreased RyR1 is commonly observed in myopathies and associated to ER stress in vitro, in mouse muscle and in human myopathy muscles, suggesting a potent role of RyR1 depletion-induced ER stress in the pathogenesis of myopathies.

A single-cell atlas of bovine skeletal muscle reveals mechanisms regulating intramuscular adipogenesis and fibrogenesis

BACKGROUND

Intramuscular fat (IMF) and intramuscular connective tissue (IMC) are often seen in human myopathies and are central to beef quality. The mechanisms regulating their accumulation remain poorly understood. Here, we explored the possibility of using beef cattle as a novel model for mechanistic studies of intramuscular adipogenesis and fibrogenesis.

METHODS

Skeletal muscle single-cell RNAseq was performed on three cattle breeds, including Wagyu (high IMF), Brahman (abundant IMC but scarce IMF), and Wagyu/Brahman cross. Sophisticated bioinformatics analyses, including clustering analysis, gene set enrichment analyses, gene regulatory network construction, RNA velocity, pseudotime analysis, and cell-cell communication analysis, were performed to elucidate heterogeneities and differentiation processes of individual cell types and differences between cattle breeds. Experiments were conducted to validate the function and specificity of identified key regulatory and marker genes. Integrated analysis with multiple published human and non-human primate datasets was performed to identify common mechanisms.

RESULTS

A total of 32 708 cells and 21 clusters were identified, including fibro/adipogenic progenitor (FAP) and other resident and infiltrating cell types. We identified an endomysial adipogenic FAP subpopulation enriched for COL4A1 and CFD (log2FC = 3.19 and 1.92, respectively; P < 0.0001) and a perimysial fibrogenic FAP subpopulation enriched for COL1A1 and POSTN (log2FC = 1.83 and 0.87, respectively; P < 0.0001), both of which were likely derived from an unspecified subpopulation. Further analysis revealed more progressed adipogenic programming of Wagyu FAPs and more advanced fibrogenic programming of Brahman FAPs. Mechanistically, NAB2 drives CFD expression, which in turn promotes adipogenesis. CFD expression in FAPs of young cattle before the onset of intramuscular adipogenesis was predictive of IMF contents in adulthood (R2  = 0.885, P < 0.01). Similar adipogenic and fibrogenic FAPs were identified in humans and monkeys. In aged humans with metabolic syndrome and progressed Duchenne muscular dystrophy (DMD) patients, increased CFD expression was observed (P < 0.05 and P < 0.0001, respectively), which was positively correlated with adipogenic marker expression, including ADIPOQ (R2  = 0.303, P < 0.01; and R2  = 0.348, P < 0.01, respectively). The specificity of Postn/POSTN as a fibrogenic FAP marker was validated using a lineage-tracing mouse line. POSTN expression was elevated in Brahman FAPs (P < 0.0001) and DMD patients (P < 0.01) but not in aged humans. Strong interactions between vascular cells and FAPs were also identified.

CONCLUSIONS

Our study demonstrates the feasibility of beef cattle as a model for studying IMF and IMC. We illustrate the FAP programming during intramuscular adipogenesis and fibrogenesis and reveal the reliability of CFD as a predictor and biomarker of IMF accumulation in cattle and humans.

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

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

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.

Spatial transcriptomics reveal markers of histopathological changes in Duchenne muscular dystrophy mouse models

Duchenne muscular dystrophy is caused by mutations in the DMD gene, leading to lack of dystrophin. Chronic muscle damage eventually leads to histological alterations in skeletal muscles. The identification of genes and cell types driving tissue remodeling is a key step to developing effective therapies. Here we use spatial transcriptomics in two Duchenne muscular dystrophy mouse models differing in disease severity to identify gene expression signatures underlying skeletal muscle pathology and to directly link gene expression to muscle histology. We perform deconvolution analysis to identify cell types contributing to histological alterations. We show increased expression of specific genes in areas of muscle regeneration (Myl4, Sparc, Hspg2), fibrosis (Vim, Fn1, Thbs4) and calcification (Bgn, Ctsk, Spp1). These findings are confirmed by smFISH. Finally, we use differentiation dynamic analysis in the D2-mdx muscle to identify muscle fibers in the present state that are predicted to become affected in the future state.

Comprehensive analysis of m6A regulators characterized by the immune microenvironment in Duchenne muscular dystrophy

BACKGROUND

Duchenne muscular dystrophy (DMD) is an X-linked, incurable, degenerative neuromuscular disease that is exacerbated by secondary inflammation. N⁶-methyladenosine (m⁶A), the most common base modification of RNA, has pleiotropic immunomodulatory effects in many diseases. However, the role of m⁶A modification in the immune microenvironment of DMD remains elusive.

METHODS

Our study retrospectively analyzed the expression data of 56 muscle tissues from DMD patients and 26 from non-muscular dystrophy individuals. Based on single sample gene set enrichment analysis, immune cells infiltration was identified and the result was validated by flow cytometry analysis and immunohistochemical staining. Then, we described the features of genetic variation in 26 m⁶A regulators and explored their relationship with the immune mircoenvironment of DMD patients through a series of bioinformatical analysis. At last, we determined subtypes of DMD patients by unsupervised clustering analysis and characterized the molecular and immune characteristics in different subgroups.

RESULTS

DMD patients have a sophisticated immune microenvironment that is significantly different from non-DMD controls. Numerous m⁶A regulators were aberrantly expressed in the muscle tissues of DMD and inversely related to most muscle-infiltrating immune cell types and immune response-related signaling pathways. A diagnostic model involving seven m⁶A regulators was established using LASSO. Furthermore, we determined three m⁶A modification patterns (cluster A/B/C) with distinct immune microenvironmental characteristics.

CONCLUSION

In summary, our study demonstrated that m⁶A regulators are intimately linked to the immune microenvironment of muscle tissues in DMD. These findings may facilitate a better understanding of the immunomodulatory mechanisms in DMD and provide novel strategies for the treatment.

Efficacy of Rectal Systemic Administration of Mesenchymal Stem Cells to Injury Sites via the CXCL12/CXCR4 Axis to Promote Regeneration in a Rabbit Skeletal Muscle Injury Model

Mesenchymal stem cells (MSCs) have been transplanted directly into lesions or injected intravenously. The administration of MSCs using these delivery methods requires specialized knowledge, techniques, and facilities. Here, we describe intrarectal systemic administration of MSCs, a simple, non-invasive route for homing to the injury sites to promote the regeneration of skeletal muscle injuries. Using a cardiotoxin (CTX)-induced rabbit skeletal muscle injury model, homing to the site of muscle injury was confirmed by intrarectal administration of MSCs; the time required for homing after intrarectal administration was approximately 5 days. In addition, the C-X-C chemokine ligand 12 (CXCL12)/C-X-C chemokine receptor-4 (CXCR4) axis was found to be involved in the homing process. Histopathological examinations showed that skeletal muscle regeneration was promoted in the MSCs-administered group compared to the CTX-only group. Myosin heavy polypeptide 3 (Myh3) expression, an indicator of early muscle regeneration, was detected earlier in the intrarectal MSCs group compared to the CTX-only group. These findings indicate that intrarectal administration of MSCs is effective in homing to the injured area, where they promote injury repair. Since intrarectal administration is a simple and non-invasive delivery route, these findings may be valuable in future research on stem cell therapy.

Duchenne muscular dystrophy: disease mechanism and therapeutic strategies

Duchenne muscular dystrophy (DMD) is a severe, progressive, and ultimately fatal disease of skeletal muscle wasting, respiratory insufficiency, and cardiomyopathy. The identification of the dystrophin gene as central to DMD pathogenesis has led to the understanding of the muscle membrane and the proteins involved in membrane stability as the focal point of the disease. The lessons learned from decades of research in human genetics, biochemistry, and physiology have culminated in establishing the myriad functionalities of dystrophin in striated muscle biology. Here, we review the pathophysiological basis of DMD and discuss recent progress toward the development of therapeutic strategies for DMD that are currently close to or are in human clinical trials. The first section of the review focuses on DMD and the mechanisms contributing to membrane instability, inflammation, and fibrosis. The second section discusses therapeutic strategies currently used to treat DMD. This includes a focus on outlining the strengths and limitations of approaches directed at correcting the genetic defect through dystrophin gene replacement, modification, repair, and/or a range of dystrophin-independent approaches. The final section highlights the different therapeutic strategies for DMD currently in clinical trials.

The Role of P2X7 Purinoceptors in the Pathogenesis and Treatment of Muscular Dystrophies

Muscular dystrophies are inherited neuromuscular diseases, resulting in progressive disability and often affecting life expectancy. The most severe, common types are Duchenne muscular dystrophy (DMD) and Limb-girdle sarcoglycanopathy, which cause advancing muscle weakness and wasting. These diseases share a common pathomechanism where, due to the loss of the anchoring dystrophin (DMD, dystrophinopathy) or due to mutations in sarcoglycan-encoding genes (LGMDR3 to LGMDR6), the α-sarcoglycan ecto-ATPase activity is lost. This disturbs important purinergic signaling: An acute muscle injury causes the release of large quantities of ATP, which acts as a damage-associated molecular pattern (DAMP). DAMPs trigger inflammation that clears dead tissues and initiates regeneration that eventually restores normal muscle function. However, in DMD and LGMD, the loss of ecto-ATPase activity, that normally curtails this extracellular ATP (eATP)-evoked stimulation, causes exceedingly high eATP levels. Thus, in dystrophic muscles, the acute inflammation becomes chronic and damaging. The very high eATP over-activates P2X7 purinoceptors, not only maintaining the inflammation but also tuning the potentially compensatory P2X7 up-regulation in dystrophic muscle cells into a cell-damaging mechanism exacerbating the pathology. Thus, the P2X7 receptor in dystrophic muscles is a specific therapeutic target. Accordingly, the P2X7 blockade alleviated dystrophic damage in mouse models of dystrophinopathy and sarcoglycanopathy. Therefore, the existing P2X7 blockers should be considered for the treatment of these highly debilitating diseases. This review aims to present the current understanding of the eATP-P2X7 purinoceptor axis in the pathogenesis and treatment of muscular dystrophies.

The P2X7 purinoceptor in pathogenesis and treatment of dystrophino- and sarcoglycanopathies

Dystrophinopathy and sarcoglycanopathies are incurable diseases caused by mutations in the genes encoding dystrophin or members of the dystrophin associated protein complex (DAPC). Restoration of the missing dystrophin or sarcoglycans via genetic approaches is complicated by the downsides of personalised medicines and immune responses against re-expressed proteins. Thus, the targeting of disease mechanisms downstream from the mutant protein has a strong translational potential. Acute muscle damage causes release of large quantities of ATP, which activates P2X7 purinoceptors, resulting in inflammation that clears dead tissues and triggers regeneration. However, in dystrophic muscles, loss of α-sarcoglycan ecto-ATPase activity further elevates extracellular ATP (eATP) levels, exacerbating the pathology. Moreover, seemingly compensatory P2X7 upregulation in dystrophic muscle cells, combined with high eATP leads to further damage. Accordingly, P2X7 blockade alleviated dystrophic damage in mouse models of both dystrophinopathy and sarcoglycanopathy. Existing P2X7 blockers could be re-purposed for the treatment of these highly debilitating diseases.

Lipin1 plays complementary roles in myofibre stability and regeneration in dystrophic muscles

Duchenne muscular dystrophy (DMD) is a severe muscle wasting disorder caused by dystrophin mutations, leading to the loss of sarcolemmal integrity, and resulting in progressive myofibre necrosis and impaired muscle function. Our previous studies suggest that lipin1 is important for skeletal muscle regeneration and myofibre integrity. Additionally, we discovered that mRNA expression levels of lipin1 were significantly reduced in skeletal muscle of DMD patients and the mdx mouse model. To understand the role of lipin1 in dystrophic muscle, we generated dystrophin/lipin1 double knockout (DKO) mice, and compared the limb muscle pathology and function of wild-type B10, muscle-specific lipin1 deficient (lipin1Myf5cKO ), mdx and DKO mice. We found that further knockout of lipin1 in dystrophic muscle exhibited a more severe phenotype characterized by increased necroptosis, fibrosis and exacerbated membrane damage in DKO compared to mdx mice. In barium chloride-induced muscle injury, both lipin1Myf5cKO and DKO showed prolonged regeneration at day 14 post-injection, suggesting that lipin1 is critical for muscle regeneration. In situ contractile function assays showed that lipin1 deficiency in dystrophic muscle led to reduced specific force production. Using a cell culture system, we found that lipin1 deficiency led to elevated expression levels of necroptotic markers and medium creatine kinase, which could be a result of sarcolemmal damage. Most importantly, restoration of lipin1 inhibited the elevation of necroptotic markers in differentiated primary lipin1-deficient myoblasts. Overall, our data suggests that lipin1 plays complementary roles in myofibre stability and muscle function in dystrophic muscles, and overexpression of lipin1 may serve as a potential therapeutic strategy for dystrophic muscles. KEY POINTS: We identified that lipin1 mRNA expression levels are significantly reduced in skeletal muscles of Duchenne muscular dystrophy patients and mdx mice. We found that further depletion of lipin1 in skeletal muscles of mdx mice induces more severe dystrophic phenotypes, including enhanced myofibre sarcolemma damage, muscle necroptosis, inflammation, fibrosis and reduced specific force production. Lipin1 deficiency leads to elevated expression levels of necroptotic markers, whereas restoration of lipin1 inhibits their expression. Our results suggest that lipin1 is functionally complementary to dystrophin in muscle membrane integrity and muscle regeneration.

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.

Identification of hub genes related to Duchenne muscular dystrophy by weighted gene co-expression network analysis

BACKGROUND

The study was aimed to analyze the potential gene modules and hub genes of Duchenne muscular dystrophy (DMD) by weighted gene co-expression network analysis.

METHODS

Based on the muscular dystrophy tissue expression profiling microarray GSE13608 from gene expression omnibus, gene co-expression modules were analyzed using weighted gene co-expression network analysis, gene modules related to DMD were screened, gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses were performed, and signature genes in the modules were screened. The protein-protein interaction network was constructed through Cytoscape, and hub genes were identified. The expression of hub genes in DMD versus normal muscle tissue was calculated in GSE6011.

RESULTS

12 co-expressed gene modules were identified, among which black module was significantly related to DMD. The characteristic genes in the module were enriched in the regulation of immune effector processes, immune response mediated by immunoglobulin, immune response mediated by B cells, etc. SERPING1, F13A1, C1S, C1R, and HLA-DPA1 were considered as hub genes in protein-protein interaction network. Analysis of GSE6011 shows that expression of SERPING1, F13A1, C1S, C1R, and HLA-DPA1 in tissues of DMD patients were higher than normal.

CONCLUSION

SERPING1, F13A1, C1S, C1R, and HLA-DPA1 may participate in the development of DMD by regulating innate immunity and inflammation, and they are expected to be a potential biomarker and novel therapeutic targets for DMD.

Ketogenic diet containing medium-chain triglyceride ameliorates transcriptome disruption in skeletal muscles of rat models of duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a myopathy characterized by progressive muscle weakness caused by a mutation in the dystrophin gene on the X chromosome. We recently showed that a medium-chain triglyceride-containing ketogenic diet (MCTKD) improves skeletal muscle myopathy in a CRISPR/Cas9 gene-edited rat model of DMD. We examined the effects of the MCTKD on transcription profiles in skeletal muscles of the model rats to assess the underlying mechanism of the MCTKD-induced improvement in DMD. DMD rats were fed MCTKD or normal diet (ND) from weaning to 9 months, and wild-type rats were fed with the ND, then tibialis anterior muscles were sampled for mRNA-seq analysis. Pearson correlation heatmaps revealed a one-node transition in the expression profile between DMD and wild-type rats. A total of 10,440, 11,555 and 11,348 genes were expressed in the skeletal muscles of wild-type and ND-fed DMD rats the MCTKD-fed DMD rats, respectively. The MCTKD reduced the number of DMD-specific mRNAs from 1624 to 1350 and increased the number of mRNAs in common with wild-type rats from 9931 to 9998. Among 2660 genes were differentially expressed in response to MCTKD intake, the mRNA expression of 1411 and 1249 of them was respectively increased and decreased. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses suggested that the MCTKD significantly suppressed the mRNA expression of genes associated with extracellular matrix organization and inflammation. This suggestion was consistent with our previous findings that the MCTKD significantly suppressed fibrosis and inflammation in DMD rats. In contrast, the MCTKD significantly increased the mRNA expression of genes associated with oxidative phosphorylation and ATP production pathways, suggesting altered energy metabolism. The decreased and increased mRNA expression of Sln and Atp2a1 respectively suggested that Sarco/endoplasmic reticulum Ca²⁺-ATPase activation is involved in the MCTKD-induced improvement of skeletal muscle myopathy in DMD rats. This is the first report to examine transcription profiles in the skeletal muscle of CRISPR/Cas9 gene-edited DMD model rats and the effect of MCTKD feeding on it.

•

We evaluated the effects of an MCTKD on the global transcriptome of DMD rats.

•

DMD rats are suitable models of human DMD for assessing transcriptome changes.

•

MCTKD suppressed fibrosis and inflammatory pathways at the transcriptional level.

•

MCTKD upregulated oxidative phosphorylation and ATP production pathways.

•

MCTKD might activate SERCA at the transcriptional level.

Serum inflammatory cytokines as disease biomarkers in the DE50-MD dog model of Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a fatal muscle-wasting disease, caused by mutations in the dystrophin gene, characterised by cycles of muscle degeneration, inflammation and regeneration. Recently, there has been renewed interest specifically in drugs that ameliorate muscle inflammation in DMD patients. The DE50-MD dog is a model of DMD that closely mimics the human DMD phenotype. We quantified inflammatory proteins in serum from wild-type (WT) and DE50-MD dogs aged 3-18 months to identify biomarkers for future pre-clinical trials. Significantly higher concentrations of C-C motif chemokine ligand 2 (CCL2), granulocyte-macrophage colony-stimulating factor (GM-CSF or CSF2), keratinocyte chemotactic-like (KC-like, homologous to mouse CXCL1), TNFα (or TNF), and interleukins IL2, IL6, IL7, IL8 (CXCL8), IL10, IL15 and IL18 were detected in DE50-MD serum compared to WT serum. Of these, CCL2 best differentiated the two genotypes. The relative level of CCL2 mRNA was greater in the vastus lateralis muscle of DE50-MD dogs than in that of WT dogs, and CCL2 was expressed both within and at the periphery of damaged myofibres. Serum CCL2 concentration was significantly associated with acid phosphatase staining in vastus lateralis biopsy samples in DE50-MD dogs. In conclusion, the serum cytokine profile suggests that inflammation is a feature of the DE50-MD phenotype. Quantification of serum CCL2 in particular is a useful non-invasive biomarker of the DE50-MD phenotype.

Muscle cell-derived cytokines in skeletal muscle regeneration

Regeneration of the mammalian adult skeletal muscle is a well-orchestrated process regulated by multiple proteins and signalling pathways. Cytokines constitute a major class of regulators of skeletal myogenesis. It is well established that infiltrating immune cells at the site of muscle injury secrete cytokines, which play critical roles in the myofibre repair and regeneration process. In the past 10-15 years, skeletal muscle itself has emerged as a prolific producer of cytokines. Much attention in the field has been focused on the endocrine effects of muscle-secreted cytokines (myokines) on metabolic regulation. However, ample evidence suggests that muscle-derived cytokines also regulate myogenic differentiation and muscle regeneration in an autocrine manner. In this review, we survey cytokines that meet two criteria: (a) evidence of expression by muscle cells; (b) evidence demonstrating a myogenic function. Dozens of cytokines representing several major classes make up this group, and together they regulate all steps of the myogenic process. How such a large array of cytokines coordinate their signalling to form a regulatory network is a fascinating, pressing question. Functional studies that can distinguish the source of the cytokines in vivo are also much needed in order to facilitate exploration of their full therapeutic potential.

Protocol to identify functional doppelgängers and verify biomedical gene expression data using doppelgangerIdentifier

Functional doppelgängers (FDs) are independently derived sample pairs that confound machine learning model (ML) performance when assorted across training and validation sets. Here, we detail the use of doppelgangerIdentifier (DI), providing software installation, data preparation, doppelgänger identification, and functional testing steps. We demonstrate examples with biomedical gene expression data. We also provide guidelines for the selection of user-defined function arguments. For complete details on the use and execution of this protocol, please refer to Wang et al. (2022).

Dlk1-Dio3 cluster miRNAs regulate mitochondrial functions in the dystrophic muscle in Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is a severe muscle disease caused by impaired expression of dystrophin. Whereas mitochondrial dysfunction is thought to play an important role in DMD, the mechanism of this dysfunction remains to be clarified. Here we demonstrate that in DMD and other muscular dystrophies, a large number of Dlk1-Dio3 clustered miRNAs (DD-miRNAs) are coordinately up-regulated in regenerating myofibers and in the serum. To characterize the biological effect of this dysregulation, 14 DD-miRNAs were simultaneously overexpressed in vivo in mouse muscle. Transcriptomic analysis revealed highly similar changes between the muscle ectopically overexpressing 14 DD-miRNAs and the mdx diaphragm, with naturally up-regulated DD-miRNAs. Among the commonly dysregulated pathway we found repressed mitochondrial metabolism, and oxidative phosphorylation (OxPhos) in particular. Knocking down the DD-miRNAs in iPS-derived skeletal myotubes resulted in increased OxPhos activities. The data suggest that (1) DD-miRNAs are important mediators of dystrophic changes in DMD muscle, (2) mitochondrial metabolism and OxPhos in particular are targeted in DMD by coordinately up-regulated DD-miRNAs. These findings provide insight into the mechanism of mitochondrial dysfunction in muscular dystrophy.

Loss of full-length dystrophin expression results in major cell-autonomous abnormalities in proliferating myoblasts

Duchenne muscular dystrophy (DMD) affects myofibers and muscle stem cells, causing progressive muscle degeneration and repair defects. It was unknown whether dystrophic myoblasts-the effector cells of muscle growth and regeneration-are affected. Using transcriptomic, genome-scale metabolic modelling and functional analyses, we demonstrate, for the first time, convergent abnormalities in primary mouse and human dystrophic myoblasts. In Dmdmdx myoblasts lacking full-length dystrophin, the expression of 170 genes was significantly altered. Myod1 and key genes controlled by MyoD (Myog, Mymk, Mymx, epigenetic regulators, ECM interactors, calcium signalling and fibrosis genes) were significantly downregulated. Gene ontology analysis indicated enrichment in genes involved in muscle development and function. Functionally, we found increased myoblast proliferation, reduced chemotaxis and accelerated differentiation, which are all essential for myoregeneration. The defects were caused by the loss of expression of full-length dystrophin, as similar and not exacerbated alterations were observed in dystrophin-null Dmdmdx-βgeo myoblasts. Corresponding abnormalities were identified in human DMD primary myoblasts and a dystrophic mouse muscle cell line, confirming the cross-species and cell-autonomous nature of these defects. The genome-scale metabolic analysis in human DMD myoblasts showed alterations in the rate of glycolysis/gluconeogenesis, leukotriene metabolism, and mitochondrial beta-oxidation of various fatty acids. These results reveal the disease continuum: DMD defects in satellite cells, the myoblast dysfunction affecting muscle regeneration, which is insufficient to counteract muscle loss due to myofiber instability. Contrary to the established belief, our data demonstrate that DMD abnormalities occur in myoblasts, making these cells a novel therapeutic target for the treatment of this lethal disease.

NR1D1 controls skeletal muscle calcium homeostasis through myoregulin repression

The sarcoplasmic reticulum (SR) plays an important role in calcium homeostasis. SR calcium mishandling is described in pathological conditions such as myopathies. Here, we investigated whether the nuclear receptor subfamily 1 group D member (NR1D1, also called REV-ERBα) regulates skeletal muscle SR calcium homeostasis. Our data demonstrate that NR1D1 deficiency in mice impairs SERCA-dependent SR calcium uptake. NR1D1 acts on calcium homeostasis by repressing the SERCA inhibitor myoregulin through direct binding to its promoter. Restoration of myoregulin counteracts the effects of NR1D1 overexpression on SR calcium content. Interestingly, myoblasts from Duchenne myopathy patients display lower NR1D1 expression, whereas pharmacological NR1D1 activation ameliorates SR calcium homeostasis, and improves muscle structure and function in dystrophic mdx/Utr+/- mice. Our findings demonstrate that NR1D1 regulates muscle SR calcium homeostasis, pointing to its therapeutic interest for mitigating myopathy.

P2X7 Purinoceptor Affects Ectopic Calcification of Dystrophic Muscles

Ectopic calcification (EC) of myofibers is a pathological feature of muscle damage in Duchenne muscular dystrophy (DMD). Mineralisation of muscle tissue occurs concomitantly with macrophage infiltration, suggesting a link between ectopic mineral deposition and inflammation. One potential link is the P2X7 purinoceptor, a key trigger of inflammation, which is expressed on macrophages but also up-regulated in dystrophic muscle cells. To investigate the role of P2X7 in dystrophic calcification, we utilised the Dmd mdx-βgeo dystrophin-null mouse model of DMD crossed with a global P2X7 knockout (P2rx7 -/- ) or with our novel P2X7 knockin-knockout mouse (P2x7 KiKo ), which expresses P2X7 in macrophages but not muscle cells. Total loss of P2X7 increased EC, indicating that P2X7 overexpression is a protective mechanism against dystrophic mineralisation. Given that muscle-specific P2X7 ablation did not affect dystrophic EC, this underlined the role of P2X7 receptor expression on the inflammatory cells. Serum phosphate reflected dystrophic calcification, with the highest serum phosphate levels found in genotypes with the most ectopic mineral. To further investigate the underlying mechanisms, we measured phosphate release from cells in vitro, and found that dystrophic myoblasts released less phosphate than non-dystrophic cells. Treatment with P2X7 antagonists increased phosphate release from both dystrophic and control myoblasts indicating that muscle cells are a potential source of secreted phosphate while macrophages protect against ectopic mineralisation. Treatment of cells with high phosphate media engendered mineral deposition, which was decreased in the presence of the P2X7 agonist BzATP, particularly in cultures of dystrophic cells, further supporting a protective role for P2X7 against ectopic mineralisation in dystrophic muscle.

Comprehensive multi-cohort transcriptional meta-analysis of muscle diseases identifies a signature of disease severity

Muscle diseases share common pathological features suggesting common underlying mechanisms. We hypothesized there is a common set of genes dysregulated across muscle diseases compared to healthy muscle and that these genes correlate with severity of muscle disease. We performed meta-analysis of transcriptional profiles of muscle biopsies from human muscle diseases and healthy controls. Studies obtained from public microarray repositories fulfilling quality criteria were divided into six categories: (i) immobility, (ii) inflammatory myopathies, (iii) intensive care unit (ICU) acquired weakness (ICUAW), (iv) congenital muscle diseases, (v) chronic systemic diseases, (vi) motor neuron disease. Patient cohorts were separated in discovery and validation cohorts retaining roughly equal proportions of samples for the disease categories. To remove bias towards a specific muscle disease category we repeated the meta-analysis five times by removing data sets corresponding to one muscle disease class at a time in a "leave-one-disease-out" analysis. We used 636 muscle tissue samples from 30 independent cohorts to identify a 52 gene signature (36 up-regulated and 16 down-regulated genes). We validated the discriminatory power of this signature in 657 muscle biopsies from 12 additional patient cohorts encompassing five categories of muscle diseases with an area under the receiver operating characteristic curve of 0.91, 83% sensitivity, and 85.3% specificity. The expression score of the gene signature inversely correlated with quadriceps muscle mass (r = -0.50, p-value = 0.011) in ICUAW and shoulder abduction strength (r = -0.77, p-value = 0.014) in amyotrophic lateral sclerosis (ALS). The signature also positively correlated with histologic assessment of muscle atrophy in ALS (r = 0.88, p-value = 1.62 × 10^-3) and fibrosis in muscular dystrophy (Jonckheere trend test p-value = 4.45 × 10^-9). Our results identify a conserved transcriptional signature associated with clinical and histologic muscle disease severity. Several genes in this conserved signature have not been previously associated with muscle disease severity.

Doppelgänger Spotting in Biomedical Gene Expression Data

Doppelgänger effects (DEs) occur when samples exhibit chance similarities such that, when split across training and validation sets, inflates the trained machine learning (ML) model performance. This inflationary effect causes misleading confidence on the deployability of the model. Thus, so far, there are no tools for doppelgänger identification or standard practices to manage their confounding implications. We present doppelgangerIdentifier, a software suite for doppelgänger identification and verification. Applying doppelgangerIdentifier across a multitude of diseases and data types, we show the pervasive nature of DEs in biomedical gene expression data. We also provide guidelines toward proper doppelgänger identification by exploring the ramifications of lingering batch effects from batch imbalances on the sensitivity of our doppelgänger identification algorithm. We suggest doppelgänger verification as a useful procedure to establish baselines for model evaluation that may inform on whether feature selection and ML on the data set may yield meaningful insights.

•

Doppelgänger effects inflate the machine learning performance

•

Doppelgänger effects exist in RNA-Seq and microarray gene expression data

•

Developed doppelgangerIdentifier, a software to identify and verify doppelgängers

•

Provide guidelines for proper doppelgänger identification

Skeletal muscle-specific overexpression of miR-486 limits mammary tumor-induced skeletal muscle functional limitations

miR-486 is a myogenic microRNA, and its reduced skeletal muscle expression is observed in muscular dystrophy. Transgenic overexpression of miR-486 using muscle creatine kinase promoter (MCK-miR-486) partially rescues muscular dystrophy phenotype. We had previously demonstrated reduced circulating and skeletal muscle miR-486 levels with accompanying skeletal muscle defects in mammary tumor models. To determine whether skeletal muscle miR-486 is functionally similar in dystrophies and cancer, we performed functional limitations and biochemical studies of skeletal muscles of MMTV-Neu mice that mimic HER2+ breast cancer and MMTV-PyMT mice that mimic luminal subtype B breast cancer and these mice crossed to MCK-miR-486 mice. miR-486 significantly prevented tumor-induced reduction in muscle contraction force, grip strength, and rotarod performance in MMTV-Neu mice. In this model, miR-486 reversed cancer-induced skeletal muscle changes, including loss of p53, phospho-AKT, and phospho-laminin alpha 2 (LAMA2) and gain of hnRNPA0 and SRSF10 phosphorylation. LAMA2 is a part of the dystrophin-associated glycoprotein complex, and its loss of function causes congenital muscular dystrophy. Complementing these beneficial effects on muscle, miR-486 indirectly reduced tumor growth and improved survival, which is likely due to systemic effects of miR-486 on production of pro-inflammatory cytokines such as IL-6. Thus, similar to dystrophy, miR-486 has the potential to reverse skeletal muscle defects and cancer burden.

Identification of Auxiliary Biomarkers and Description of the Immune Microenvironmental Characteristics in Duchenne Muscular Dystrophy by Bioinformatical Analysis and Experiment

Background

Duchenne muscular dystrophy (DMD) is a genetic muscle disorder characterized by progressive muscle wasting associated with persistent inflammation. In this study, we aimed to identify auxiliary biomarkers and further characterize the immune microenvironment in DMD.

Methods

Differentially expressed genes (DEGs) were identified between DMD and normal muscle tissues based on Gene Expression Omnibus (GEO) datasets. Bioinformatical analysis was used to screen and identify potential diagnostic signatures of DMD which were further validated by real-time quantitative reverse transcription PCR (RT-qPCR). We also performed single-sample gene-set enrichment analysis (ssGSEA) to characterize the proportion of tissue-infiltrating immune cells to determine the inflammatory state of DMD.

Results

In total, 182 downregulated genes and 263 upregulated genes were identified in DMD. C3, SPP1, TMSB10, TYROBP were regarded as adjunct biomarkers and successfully validated by RT-qPCR. The infiltration of macrophages, CD4+, and CD8+ T cells was significantly higher (p < 0.05) in DMD compared with normal muscle tissues, while the infiltration of activated B cells, CD56dim natural killer cells, and type 17 T helper (Th17) cells was lower. In addition, the four biomarkers (C3, SPP1, TMSB10, TYROBP) were strongly associated with immune cells and immune-related pathways in DMD muscle tissues.

Conclusion

Analyses demonstrated C3, SPP1, TMSB10, and TYROBP may serve as biomarkers and enhance our understanding of immune responses in DMD. The infiltration of immune cells into the muscle microenvironment might exert a critical impact on the development and occurrence of DMD.

An Investigation of How Normalisation and Local Modelling Techniques Confound Machine Learning Performance In a Mental Health Study

Machine learning (ML) is increasingly deployed on biomedical studies for biomarker development (feature selection) and diagnostic/prognostic technologies (classification). While different ML techniques produce different feature sets and classification performances, less understood is how upstream data processing methods (e.g., normalisation) impact downstream analyses. Using a clinical mental health dataset, we investigated the impact of different normalisation techniques on classification model performance. Gene Fuzzy Scoring (GFS), an in-house developed normalisation technique, is compared against widely used normalisation methods such as global quantile normalisation, class-specific quantile normalisation and surrogate variable analysis. We report that choice of normalisation technique has strong influence on feature selection. with GFS outperforming other techniques. Although GFS parameters are tuneable, good classification model performance (ROC-AUC > 0.90) is observed regardless of the GFS parameter settings. We also contrasted our results against local modelling, which is meant to improve the resolution and meaningfulness of classification models built on heterogeneous data. Local models, when derived from non-biologically meaningful subpopulations, perform worse than global models. A deep dive however, revealed that the factors driving cluster formation has little to do with the phenotype-of-interest. This finding is critical, as local models are often seen as a superior means of clinical data modelling. We advise against such naivete. Additionally, we have developed a combinatorial reasoning approach using both global and local paradigms: This helped reveal potential data quality issues or underlying factors causing data heterogeneity that are often overlooked. It also assists to explain the model as well as provides directions for further improvement.

Beneficial impacts of neuromuscular electrical stimulation on muscle structure and function in the zebrafish model of Duchenne muscular dystrophy

Neuromuscular electrical stimulation (NMES) allows activation of muscle fibers in the absence of voluntary force generation. NMES could have the potential to promote muscle homeostasis in the context of muscle disease, but the impacts of NMES on diseased muscle are not well understood. We used the zebrafish Duchenne muscular dystrophy (dmd) mutant and a longitudinal design to elucidate the consequences of NMES on muscle health. We designed four neuromuscular stimulation paradigms loosely based on weightlifting regimens. Each paradigm differentially affected neuromuscular structure, function, and survival. Only endurance neuromuscular stimulation (eNMES) improved all outcome measures. We found that eNMES improves muscle and neuromuscular junction morphology, swimming, and survival. Heme oxygenase and integrin alpha7 are required for eNMES-mediated improvement. Our data indicate that neuromuscular stimulation can be beneficial, suggesting that the right type of activity may benefit patients with muscle disease.

Treating Duchenne Muscular Dystrophy: The Promise of Stem Cells, Artificial Intelligence, and Multi-Omics

Muscular dystrophies are chronic and debilitating disorders caused by progressive muscle wasting. Duchenne muscular dystrophy (DMD) is the most common type. DMD is a well-characterized genetic disorder caused by the absence of dystrophin. Although some therapies exist to treat the symptoms and there are ongoing efforts to correct the underlying molecular defect, patients with muscular dystrophies would greatly benefit from new therapies that target the specific pathways contributing directly to the muscle disorders. Three new advances are poised to change the landscape of therapies for muscular dystrophies such as DMD. First, the advent of human induced pluripotent stem cells (iPSCs) allows researchers to design effective treatment strategies that make up for the gaps missed by conventional “one size fits all” strategies. By characterizing tissue alterations with single-cell resolution and having molecular profiles for therapeutic treatments for a variety of cell types, clinical researchers can design multi-pronged interventions to not just delay degenerative processes, but regenerate healthy tissues. Second, artificial intelligence (AI) will play a significant role in developing future therapies by allowing the aggregation and synthesis of large and disparate datasets to help reveal underlying molecular mechanisms. Third, disease models using a high volume of multi-omics data gathered from diverse sources carry valuable information about converging and diverging pathways. Using these new tools, the results of previous and emerging studies will catalyze precision medicine-based drug development that can tackle devastating disorders such as DMD.

Inhibition of sphingolipid de novo synthesis counteracts muscular dystrophy

Duchenne muscular dystrophy (DMD), the most common muscular dystrophy, is a severe muscle disorder, causing muscle weakness, loss of independence, and premature death. Here, we establish the link between sphingolipids and muscular dystrophy. Transcripts of sphingolipid de novo biosynthesis pathway are up-regulated in skeletal muscle of patients with DMD and other muscular dystrophies, which is accompanied by accumulation of metabolites of the sphingolipid pathway in muscle and plasma. Pharmacological inhibition of sphingolipid synthesis by myriocin in the mdx mouse model of DMD ameliorated the loss in muscle function while reducing inflammation, improving Ca²⁺ homeostasis, preventing fibrosis of the skeletal muscle, heart, and diaphragm, and restoring the balance between M1 and M2 macrophages. Myriocin alleviated the DMD phenotype more than glucocorticoids. Our study identifies inhibition of sphingolipid synthesis, targeting multiple pathogenetic pathways simultaneously, as a strong candidate for treatment of muscular dystrophies.

Skeletal Ryanodine Receptors Are Involved in Impaired Myogenic Differentiation in Duchenne Muscular Dystrophy Patients

Duchenne muscular dystrophy (DMD) is characterized by progressive muscle wasting following repeated muscle damage and inadequate regeneration. Impaired myogenesis and differentiation play a major role in DMD as well as intracellular calcium (Ca2+) mishandling. Ca2+ release from the sarcoplasmic reticulum is mostly mediated by the type 1 ryanodine receptor (RYR1) that is required for skeletal muscle differentiation in animals. The study objective was to determine whether altered RYR1-mediated Ca2+ release contributes to myogenic differentiation impairment in DMD patients. The comparison of primary cultured myoblasts from six boys with DMD and five healthy controls highlighted delayed myoblast differentiation in DMD. Silencing RYR1 expression using specific si-RNA in a healthy control induced a similar delayed differentiation. In DMD myotubes, resting intracellular Ca2+ concentration was increased, but RYR1-mediated Ca2+ release was not changed compared with control myotubes. Incubation with the RYR-calstabin interaction stabilizer S107 decreased resting Ca2+ concentration in DMD myotubes to control values and improved calstabin1 binding to the RYR1 complex. S107 also improved myogenic differentiation in DMD. Furthermore, intracellular Ca2+ concentration was correlated with endomysial fibrosis, which is the only myopathologic parameter associated with poor motor outcome in patients with DMD. This suggested a potential relationship between RYR1 dysfunction and motor impairment. Our study highlights RYR1-mediated Ca2+ leakage in human DMD myotubes and its key role in myogenic differentiation impairment. RYR1 stabilization may be an interesting adjunctive therapeutic strategy in DMD.

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.

Transcriptome analysis of collagen VI-related muscular dystrophy muscle biopsies

Objective

To define the transcriptomic changes responsible for the histologic alterations in skeletal muscle and their progression in collagen VI‐related muscular dystrophy (COL6‐RD).

Methods

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.

Results

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.

Interpretation

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.

Disrupted Calcium Homeostasis in Duchenne Muscular Dystrophy: A Common Mechanism behind Diverse Consequences

Duchenne muscular dystrophy (DMD) leads to disability and death in young men. This disease is caused by mutations in the DMD gene encoding diverse isoforms of dystrophin. Loss of full-length dystrophins is both necessary and sufficient for causing degeneration and wasting of striated muscles, neuropsychological impairment, and bone deformities. Among this spectrum of defects, abnormalities of calcium homeostasis are the common dystrophic feature. Given the fundamental role of Ca2+ in all cells, this biochemical alteration might be underlying all the DMD abnormalities. However, its mechanism is not completely understood. While abnormally elevated resting cytosolic Ca2+ concentration is found in all dystrophic cells, the aberrant mechanisms leading to that outcome have cell-specific components. We probe the diverse aspects of calcium response in various affected tissues. In skeletal muscles, cardiomyocytes, and neurons, dystrophin appears to serve as a scaffold for proteins engaged in calcium homeostasis, while its interactions with actin cytoskeleton influence endoplasmic reticulum organisation and motility. However, in myoblasts, lymphocytes, endotheliocytes, and mesenchymal and myogenic cells, calcium abnormalities cannot be clearly attributed to the loss of interaction between dystrophin and the calcium toolbox proteins. Nevertheless, DMD gene mutations in these cells lead to significant defects and the calcium anomalies are a symptom of the early developmental phase of this pathology. As the impaired calcium homeostasis appears to underpin multiple DMD abnormalities, understanding this alteration may lead to the development of new therapies. In fact, it appears possible to mitigate the impact of the abnormal calcium homeostasis and the dystrophic phenotype in the total absence of dystrophin. This opens new treatment avenues for this incurable disease.

Inflammation in Duchenne Muscular Dystrophy-Exploring the Role of Neutrophils in Muscle Damage and Regeneration

Duchenne muscular dystrophy (DMD) is a severe and progressive, X-linked, neuromuscular disorder caused by mutations in the dystrophin gene. In DMD, the lack of functional dystrophin protein makes the muscle membrane fragile, leaving the muscle fibers prone to damage during contraction. Muscle degeneration in DMD patients is closely associated with a prolonged inflammatory response, and while this is important to stimulate regeneration, inflammation is also thought to exacerbate muscle damage. Neutrophils are one of the first immune cells to be recruited to the damaged muscle and are the first line of defense during tissue injury or infection. Neutrophils can promote inflammation by releasing pro-inflammatory cytokines and compounds, including myeloperoxidase (MPO) and neutrophil elastase (NE), that lead to oxidative stress and are thought to have a role in prolonging inflammation in DMD. In this review, we provide an overview of the roles of the innate immune response, with particular focus on mechanisms used by neutrophils to exacerbate muscle damage and impair regeneration in DMD.

C-X-C motif chemokine ligand 12: a potential therapeutic target in Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is an X-linked recessive disease caused by a mutant dystrophin protein. DMD patients undergo gradual progressive paralysis until death. Chronic glucocorticoid therapy remains one of the main treatments for DMD, despite the significant side effects. However, its mechanisms of action remain largely unknown. We used bioinformatics tools to identify pathogenic genes involved in DMD and glucocorticoid target genes. Two gene expression profiles containing data from DMD patients and healthy controls (GSE38417 and GSE109178) were downloaded for further analysis. Differentially expressed genes (DEGs) between DMD patients and controls were identified using GEO2R, and glucocorticoid target genes were predicted from the Pharmacogenetics and Pharmacogenomics Knowledge Base. Surprisingly, only one gene, CXCL12 (C-X-C motif chemokine ligand 12), was both a glucocorticoid target and a DEG. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, Gene Ontology term enrichment analysis, and gene set enrichment analysis were performed. A protein-protein interaction network was constructed and hub genes identified using the Search Tool for the Retrieval of Interacting Genes (STRING) database and Cytoscape. Enriched pathways involving the DEGs, including CXCL12, were associated with the immune response and inflammation. Levels of CXCL12 and its receptor CXCR4 (C-X-C motif chemokine receptor 4) were increased in X-linked muscular dystrophy (mdx) mice (DMD models) but became significantly reduced after prednisone treatment. Metformin also reduced the expression of CXCL12 and CXCR4 in mdx mice. In conclusion, the CXCL12-CXCR4 pathway may be a potential target for DMD therapy.

Multiomic Approaches to Uncover the Complexities of Dystrophin-Associated Cardiomyopathy

Despite major progress in treating skeletal muscle disease associated with dystrophinopathies, cardiomyopathy is emerging as a major cause of death in people carrying dystrophin gene mutations that remain without a targeted cure even with new treatment directions and advances in modelling abilities. The reasons for the stunted progress in ameliorating dystrophin-associated cardiomyopathy (DAC) can be explained by the difficulties in detecting pathophysiological mechanisms which can also be efficiently targeted within the heart in the widest patient population. New perspectives are clearly required to effectively address the unanswered questions concerning the identification of authentic and effectual readouts of DAC occurrence and severity. A potential way forward to achieve further therapy breakthroughs lies in combining multiomic analysis with advanced preclinical precision models. This review presents the fundamental discoveries made using relevant models of DAC and how omics approaches have been incorporated to date.

Validation of Chemokine Biomarkers in Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is a progressive muscle disease involving complex skeletal muscle pathogenesis. The pathogenesis is triggered by sarcolemma instability due to the lack of dystrophin protein expression, leading to Ca²⁺ influx, muscle fiber apoptosis, inflammation, muscle necrosis, and fibrosis. Our lab recently used two high-throughput multiplexing techniques (e.g., SomaScan^® aptamer assay and tandem mass tag-(TMT) approach) and identified a series of serum protein biomarkers tied to different pathobiochemical pathways. In this study, we focused on validating the circulating levels of three proinflammatory chemokines (CCL2, CXCL10, and CCL18) that are believed to be involved in an early stage of muscle pathogenesis. We used highly specific and reproducible MSD ELISA assays and examined the association of these chemokines with DMD pathogenesis, age, disease severity, and response to glucocorticoid treatment. As expected, we confirmed that these three chemokines were significantly elevated in serum and muscle samples of DMD patients relative to age-matched healthy controls (p-value < 0.05, CCL18 was not significantly altered in muscle samples). These three chemokines were not significantly elevated in Becker muscular dystrophy (BMD) patients, a milder form of dystrophinopathy, when compared in a one-way ANOVA to a control group but remained significantly elevated in the age-matched DMD group (p < 0.05). CCL2 and CCL18 but not CXCL10 declined with age in DMD patients, whereas all three chemokines remained unchanged with age in BMD and controls. Only CCL2 showed significant association with time to climb four steps in the DMD group (r = 0.48, p = 0.038) and neared significant association with patients' reported outcome in the BMD group (r = 0.39, p = 0.058). Furthermore, CCL2 was found to be elevated in a serum of the mdx mouse model of DMD, relative to wild-type mouse model. This study suggests that CCL2 might be a suitable candidate biomarker for follow-up studies to demonstrate its physiological significance and clinical utility in DMD.

Urolithin A improves muscle function by inducing mitophagy in muscular dystrophy

Duchenne muscular dystrophy (DMD) is the most common muscular dystrophy, and despite advances in genetic and pharmacological disease-modifying treatments, its management remains a major challenge. Mitochondrial dysfunction contributes to DMD, yet the mechanisms by which this occurs remain elusive. Our data in experimental models and patients with DMD show that reduced expression of genes involved in mitochondrial autophagy, or mitophagy, contributes to mitochondrial dysfunction. Mitophagy markers were reduced in skeletal muscle and in muscle stem cells (MuSCs) of a mouse model of DMD. Administration of the mitophagy activator urolithin A (UA) rescued mitophagy in DMD worms and mice and in primary myoblasts from patients with DMD, increased skeletal muscle respiratory capacity, and improved MuSCs' regenerative ability, resulting in the recovery of muscle function and increased survival in DMD mouse models. These data indicate that restoration of mitophagy alleviates symptoms of DMD and suggest that UA may have potential therapeutic applications for muscular dystrophies.

Longitudinal metabolomic analysis of plasma enables modeling disease progression in Duchenne muscular dystrophy mouse models

Duchenne muscular dystrophy is a severe pediatric neuromuscular disorder caused by the lack of dystrophin. Identification of biomarkers is needed to support and accelerate drug development. Alterations of metabolites levels in muscle and plasma have been reported in pre-clinical and clinical cross-sectional comparisons. We present here a 7-month longitudinal study comparing plasma metabolomic data in wild-type and mdx mice. A mass spectrometry approach was used to study metabolites in up to five time points per mouse at 6, 12, 18, 24 and 30 weeks of age, providing an unprecedented in depth view of disease trajectories. A total of 106 metabolites were studied. We report a signature of 31 metabolites able to discriminate between healthy and disease at various stages of the disease, covering the acute phase of muscle degeneration and regeneration up to the deteriorating phase. We show how metabolites related to energy production and chachexia (e.g. glutamine) are affected in mdx mice plasma over time. We further show how the signature is connected to molecular targets of nutraceuticals and pharmaceutical compounds currently in development as well as to the nitric oxide synthase pathway (e.g. arginine and citrulline). Finally, we evaluate the signature in a second longitudinal study in three independent mouse models carrying 0, 1 or 2 functional copies of the dystrophin paralog utrophin. In conclusion, we report an in-depth metabolomic signature covering previously identified associations and new associations, which enables drug developers to peripherally assess the effect of drugs on the metabolic status of dystrophic mice.

A Network Medicine Approach for Drug Repurposing in Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is a progressive hereditary muscular disease caused by a lack of dystrophin, leading to membrane instability, cell damage, and inflammatory response. However, gene-editing alone is not enough to restore the healthy phenotype and additional treatments are required. In the present study, we have first conducted a meta-analysis of three microarray datasets, GSE38417, GSE3307, and GSE6011, to identify the differentially expressed genes (DEGs) between healthy donors and DMD patients. We have then integrated this analysis with the knowledge obtained from DisGeNET and DIAMOnD, a well-known algorithm for drug–gene association discoveries in the human interactome. The data obtained allowed us to identify novel possible target genes and were used to predict potential therapeutical options that could reverse the pathological condition.

Dynamic and thermodynamic models of adaptation

The concept of biological adaptation was closely connected to some mathematical, engineering and physical ideas from the very beginning. Cannon in his "The wisdom of the body" (1932) systematically used the engineering vision of regulation. In 1938, Selye enriched this approach by the notion of adaptation energy. This term causes much debate when one takes it literally, as a physical quantity, i.e. a sort of energy. Selye did not use the language of mathematics systematically, but the formalization of his phenomenological theory in the spirit of thermodynamics was simple and led to verifiable predictions. In 1980s, the dynamics of correlation and variance in systems under adaptation to a load of environmental factors were studied and the universal effect in ensembles of systems under a load of similar factors was discovered: in a crisis, as a rule, even before the onset of obvious symptoms of stress, the correlation increases together with variance (and volatility). During 30 years, this effect has been supported by many observations of groups of humans, mice, trees, grassy plants, and on financial time series. In the last ten years, these results were supplemented by many new experiments, from gene networks in cardiology and oncology to dynamics of depression and clinical psychotherapy. Several systems of models were developed: the thermodynamic-like theory of adaptation of ensembles and several families of models of individual adaptation. Historically, the first group of models was based on Selye's concept of adaptation energy and used fitness estimates. Two other groups of models are based on the idea of hidden attractor bifurcation and on the advection-diffusion model for distribution of population in the space of physiological attributes. We explore this world of models and experiments, starting with classic works, with particular attention to the results of the last ten years and open questions.