Abstract
Myotonic dystrophy type 1 (DM1) is a progressive multisystemic disease caused by a CTG repeat expansion in the DMPK gene. The toxic mutant mRNA sequesters MBNL proteins, disrupting global RNA metabolism. Although alternative splicing in DM1 skeletal muscle pathology has been extensively studied, early-stage transcriptomic changes remained uncharacterized. To gain deeper and contextual insight into DM1 transcriptome, we performed the first Weighted Gene Co-expression Network Analysis (WGCNA) on skeletal muscle RNA sequencing data from the widely used DM1 mouse model HSA(LR) (~250 CTG repeats). We identified 532 core genes using data from 16-week-old mice, an age before the onset of muscle weakness. Additional differential expression analysis across multiple HSA(LR) datasets revealed 42 common up-regulated coding and non-coding genes. Within identified core genes, the pathway gene-pair signature analysis enabled contextual selection of functionally related genes involved in maintaining proteostasis, including endoplasmic reticulum (ER) protein processing, the ubiquitin-proteasome system (UPS), macroautophagy and mitophagy, and muscle contraction. The enrichment of ER protein processing with prevailing core genes related to ER-associated degradation suggests adaptive chaperone and UPS activation, while core genes such as Ambra1, Mfn2, and Usp30 indicate adaptations in mitochondrial quality control. Coordinated early alterations in processes maintaining protein homeostasis, critical for muscle mass and function, possibly reflect a response to cellular stress due to repeat expansion and appears before muscle weakness development. Although the study relies exclusively on transcriptomic analyses, it offers a comprehensive, hypothesis-generating perspective that pinpoints candidate pathways, preceding muscle weakness, for future mechanistic validation.