Abstract
Collagen VI-related dystrophies (COL6RD) are a group of rare muscle disorders caused by mutations in specific genes responsible for type VI collagen production. It affects muscles, joints, and connective tissues, leading to weakness, joint problems, and structural issues. Currently, there is no effective treatment for COL6RD; its management typically addresses symptoms and complications. Therefore, it is essential to decipher the disease's molecular mechanisms, identify drug targets, and develop effective treatment strategies to treat COL6RD. In this study, we employed differential gene expression analysis, weighted gene co-expression network analysis, and genome-scale metabolic modeling to investigate gene expression patterns in COL6RD patients, uncovering key genes, significant metabolites, and disease-related pathophysiological pathways. First, we performed differential gene expression and weighted gene co-expression network analyses, which led to the identification of 12 genes (CHCHD10, MRPS24, TRIP10, RNF123, MRPS15, NDUFB4, COX10, FUNDC2, MDH2, RPL3L, NDUFB11, PARVB) as potential hub genes involved in the disease. Second, we utilized a drug repurposing strategy to identify pharmaceutical candidates that could potentially modulate these genes and be effective in the treatment. Next, we utilized context-specific genome-scale metabolic models to compare metabolic variations between healthy individuals and COL6RD patients. Finally, we conducted reporter metabolite analysis to identify reporter metabolites (e.g., phosphatidates, nicotinate ribonucleotide, ubiquinol, ferricytochrome C). In summary, our analysis revealed critical genes and pathways associated with COL6RD and identified potential targets, reporter metabolites, and candidate drugs for therapeutic interventions.