Identification of Immune-related Features Involved in Duchenne Muscular Dystrophy: A Bidirectional Transcriptome and Proteome-driven Analysis

Background: We aimed to investigate the biological mechanism and feature genes of Duchenne muscular dystrophy (DMD) by multi-omics and experimental verification strategy.

Methods: We integrated the transcriptomic and proteomic methods to find the differentially expressed mRNAs (DEMs) and proteins (DEPs) between DMD and Control groups. Weighted gene co-expression network analysis (WGCNA) was then used to identify modules of highly correlated genes and hub genes. In the following steps, the immune and stromal cells infiltrations were accomplished by xCELL algorithm. Furthermore, TF and miRNA prediction were performed with Networkanalyst. ELISA, western blot and external datasets were performed to verify the key proteins/mRNAs in DMD patient and mouse. Finally, a nomogram model was established based on the potential biomarkers.

Results: 4515 DEMs and 56 DEPs were obtained from the transcriptomic and proteomic study respectively. 14 common genes were identified, which is enriched in muscle contraction and inflammation-related pathways. Meanwhile, we observed 33 significant differences in the infiltration of cells in DMD. Afterwards, a total of 22 miRNAs and 23 TF genes interacted with the common genes, including TFAP2C, MAX, MYC, NFKB1, RELA, hsa-miR-1255a, hsa-miR-130a, hsa-miR-130b, hsa-miR-152, and hsa-miR-17. In addition, three genes (ATP6AP2, CTSS, and VIM) showed excellent diagnostic performance on discriminating DMD in GSE1004, GSE3307, GSE6011 and GSE38417 datasets (all AUC > 0.8), which is validated in patients (10 DMD vs. 10 controls), DMD with exon 55 mutations, mdx mouse, and nomogram model.

Conclusion: Taken together, ATP6AP2, CTSS, and VIM play important roles in the inflammatory response in DMD, which may serve as diagnostic biomarkers and therapeutic targets.