Identifying the hub genes for Duchenne muscular dystrophy and Becker muscular dystrophy by weighted correlation network analysis

Abstract

Background: The goal of this study is to identify the hub genes for Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) via weighted correlation network analysis (WGCNA). Methods: The gene expression profile of vastus lateralis biopsy samples obtained in 17 patients with DMD, 11 patients with BMD and 6 healthy individuals was downloaded from the Gene Expression Omnibus (GEO) database (GSE109178). After obtaining different expressed genes (DEGs) via GEO2R, WGCNA was conducted using R package, modules and genes that highly associated with DMD, BMD, and their age or pathology were screened. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analysis and protein–protein interaction (PPI) network analysis were also conducted. Hub genes and highly correlated clustered genes were identified using Search Tool for the Retrieval of Interacting Genes (STRING) and Cystoscape software. Results: One thousand four hundred seventy DEGs were identified between DMD and control, with 1281 upregulated and 189 downregulated DEGs. Four hundred and twenty DEGs were found between BMD and control, with 157 upregulated and 263 upregulated DEGs. Fourteen modules with different colors were identified for DMD vs control, and 7 modules with different colors were identified for BMD vs control. Ten hub genes were summarized for DMD and BMD respectively, 5 hub genes were summarized for BMD age, 5 and 3 highly correlated clustered genes were summarized for DMD age and BMD pathology, respectively. In addition, 20 GO enrichments were found to be involved in DMD, 3 GO enrichments were found to be involved in BMD, 3 GO enrichments were found to be involved in BMD age. Conclusion: In DMD, several hub genes were identified: C3AR1, TLR7, IRF8, FYB and CD33(immune and inflammation associated genes), TYROBP, PLEK, AIF1(actin reorganization associated genes), LAPTM5 and NT5E(cell death and arterial calcification associated genes, respectively). In BMD, a number of hub genes were identified: LOX, ELN, PLEK, IKZF1, CTSK, THBS2, ADAMTS2, COL5A1(extracellular matrix associated genes), BCL2L1 and CDK2(cell cycle associated genes).