Assessment of Weighted Gene Co-Expression Network Analysis to Explore Key Pathways and Novel Biomarkers in Muscular Dystrophy.
Ontology highlight
ABSTRACT: Purpose
This study aimed to explore the key molecular pathways involved in Duchenne muscular dystrophy (DMD) and Becker muscular dystrophy (BMD) and thereby identify hub genes to be potentially used as novel biomarkers using a bioinformatics approach.Methods
Raw GSE109178 data were collected from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was conducted on the top 50% of altered genes. The key modules associated with the clinical features of DMD and BMD were identified. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the DAVID website. A protein-protein interaction (PPI) network was constructed using the STRING website. MCODE, together with the Cytohubba plug-ins of Cytoscape, screened out the potential hub genes, which were subsequently verified via receiver operating characteristic (ROC) curves in other datasets.Results
Among the 11 modules obtained, the black module was predominantly associated with pathology and DMD, whereas the light-green module was primarily related to age and BMD. Functional enrichment assessments indicated that the genes in the black module were primarily clustered in "immune response" and "phagosome," whereas the ones in the light-green module were chiefly enriched in "protein polyubiquitination". Eleven essential genes were eventually identified, including VCAM1, TYROBP, CD44, ITGB2, CSF1R, LCP2, C3AR1, CCL2, and ITGAM for DMD, along with UBA5 and UBR2 for BMD.Conclusion
Overall, our findings may be useful for investigating the mechanisms underlying DMD and BMD. In addition, the hub genes discovered might serve as novel molecular markers correlated with dystrophinopathies.
SUBMITTER: Xu X
PROVIDER: S-EPMC8053709 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundThe study was aimed to analyze the potential gene modules and hub genes of Duchenne muscular dystrophy (DMD) by weighted gene co-expression network analysis.MethodsBased 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.Results12 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.ConclusionSERPING1, 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.
| S-EPMC9803489 | biostudies-literature
Project description:BACKGROUND:Weighted co-expression network analysis (WGCNA) is a powerful systems biology method to describe the correlation of gene expression based on the microarray database, which can be used to facilitate the discovery of therapeutic targets or candidate biomarkers in diseases. OBJECTIVE:To explore the key genes in the development of Alzheimer's disease (AD) by using WGCNA. METHODS:The whole gene expression data GSE1297 from AD and control human hippocampus was obtained from the GEO database in NCBI. Co-expressed genes were clustered into different modules. Modules of interest were identified through calculating the correlation coefficient between the module and phenotypic traits. GO and pathway enrichment analyses were conducted, and the central players (key hub genes) within the modules of interest were identified through network analysis. The expression of the identified key genes was confirmed in AD transgenic mice through using qRT-PCR. RESULTS:Two modules were found to be associated with AD clinical severity, which functioning mainly in mineral absorption, NF-?B signaling, and cGMP-PKG signaling pathways. Through analysis of the two modules, we found that metallothionein (MT), Notch2, MSX1, ADD3, and RAB31 were highly correlated with AD phenotype. Increase in expression of these genes was confirmed in aged AD transgenic mice. CONCLUSION:WGCNA analysis can be used to analyze and predict the key genes in AD. MT1, MT2, MSX1, NOTCH2, ADD3, and RAB31 are identified to be the most relevant genes, which may be potential targets for AD therapy.
| S-EPMC6218130 | biostudies-literature
Project description:BackgroundThe 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).MethodsThe 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.ResultsOne 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.ConclusionIn 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).
| S-EPMC8684282 | biostudies-literature
Project description:Bronchopulmonary dysplasia (BPD) is a complex disorder resulting from interactions between genes and the environment. The accurate molecular etiology of BPD remains largely unclear. This study aimed to identify key BPD-associated genes and pathways functionally enriched using weighted gene co-expression network analysis (WGCNA). We analyzed microarray data of 62 pre-term patients with BPD and 38 pre-term patients without BPD from Gene Expression Omnibus (GEO). WGCNA was used to construct a gene expression network, and genes were classified into definite modules. In addition, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of BPD-related hub genes were performed. Firstly, we constructed a weighted gene co-expression network, and genes were divided into 10 modules. Among the modules, the yellow module was related to BPD progression and severity and included the following hub genes: MMP25, MMP9, SIRPA, CKAP4, SLCO4C1, and SLC2A3; and the red module included some co-expression molecules that displayed a continuous decline in expression with BPD progression and included the following hub genes: LEF1, ITK, CD6, RASGRP1, IL7R, SKAP1, CD3E, and ICOS. GO and KEGG analyses showed that high expression of inflammatory response-related genes and low expression of T cell receptor activation-related genes are significantly correlated with BPD progression. The present WGCNA-based study thus provides an overall perspective of BPD and lays the foundation for identifying potential pathways and hub genes that contribute to the development of BPD.
| S-EPMC7509191 | biostudies-literature
Project description:Colorectal cancer (CRC) is one of the most common malignant diseases in the world. Although mechanistic studies have been conducted on the pathogenesis of CRC, the molecular mechanism of CRC tumorigenesis remains unclear. In the present study, the weighted gene co?expression network analysis was performed for the Gene Expression Omnibus (GEO) dataset GSE87211, in order to analyze the key modules involved in the pathogenesis of CRC. Next, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed on the key module genes to analyze the functional pathways involved. The hub genes were screened using the Cytoscape platform and verified by a second GEO dataset, GSE21510. Finally, 10 hub genes were identified in 2 key modules (the green and brown modules) as the genes most significantly associated with the tumorigenesis of CRC. The 5 hub genes from the green module included collagen type I ?1 chain, collagen type XII ?1 chain, collagen triple helix repeat containing 1, inhibin subunit ?a (INHBA) and chromobox 2 (CBX2), while the 5 hub genes from the brown module included bestrophin 2 (BEST2), carbonic anhydrase 2, glucagon, solute carrier family 4 member 4 and gliomedin. The 2 key modules with the 10 hub genes identified may regulate the occurrence and development of CRC through the extracellular matrix pathway, PI3K?Akt and chemokine signaling pathways, thus providing a reference for understanding the complex mechanism of tumorigenesis in CRC. Of note, few studies have reported the pathogenesis of CRC with the 3 identified hub genes, INHBA, CBX2 and BEST2. Further investigation of the molecular mechanism of these genes in CRC is recommended.
| S-EPMC7185270 | biostudies-literature
Project description:ObjectiveTo reveal the molecular mechanism underlying the pathogenesis of HCM and find new effective therapeutic strategies using a systematic biological approach.MethodsThe WGCNA algorithm was applied to building the co-expression network of HCM samples. A sample cluster analysis was performed using the hclust tool and a co-expression module was constructed. The WGCNA algorithm was used to study the interactive connection between co-expression modules and draw a heat map to show the strength of interactions between modules. The genetic information of the respective modules was mapped to the associated GO terms and KEGG pathways, and the Hub Genes with the highest connectivity in each module were identified. The Wilcoxon test was used to verify the expression level of hub genes between HCM and normal samples, and the "pROC" R package was used to verify the possibility of hub genes as biomarkers. Finally, the potential functions of hub genes were analyzed by GSEA software.ResultsSeven co-expression modules were constructed using sample clustering analysis. GO and KEGG enrichment analysis judged that the turquoise module is an important module. The hub genes of each module are RPL35A for module Black, FH for module Blue, PREI3 for module Brown, CREB1 for module Green, LOC641848 for module Pink, MYH7 for module Turquoise and MYL6 for module Yellow. The results of the differential expression analysis indicate that MYH7 and FH are considered true hub genes. In addition, the ROC curves revealed their high diagnostic value as biomarkers for HCM. Finally, in the results of the GSEA analysis, MYH7 and FH highly expressed genes were enriched with the "proteasome" and a "PPAR signaling pathway," respectively.ConclusionsThe MYH7 and FH genes may be the true hub genes of HCM. Their respective enriched pathways, namely the "proteasome" and the "PPAR signaling pathway," may play an important role in the development of HCM.
| S-EPMC7585681 | biostudies-literature
Project description:PurposeAt present, there is a lack of accurate early diagnostic markers for ischemic stroke.MethodsBy using dimensionality reduction cluster analysis, differential expression analysis, weighted co-expression network analysis, protein-protein interaction network analysis, cell heterogeneity and key pathogenic genes were identified in ischemic stroke. Immunomicroenvironment analysis was used to explore the immune landscape and immune associations of key genes in ischemic stroke. The analysis platform we use is R software (version 4.0.5). PCR experiments were used to verify the expression of key genes.ResultsSingle cell sequencing data in ischemic stroke can be annotated as fibroblast cells, pre-B cell CD34, neutrophils cells, bone marrow (BM), keratinocytes, macrophage, neurons and mesenchymal stem cells (MSC). By the intersection of differential expression analysis and WGCNA analysis, 385 genes were obtained. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these genes were highly correlated with multiple functions and pathways. Protein-protein interaction network analysis revealed that MRPS11 and MRPS12 were key genes, both of which were down-regulated in ischemic stroke. The Pseudo-time series analysis found that the expression of MRPS12 decreased gradually with the differentiation of pre-B cell CD34 cells in ischemic stroke, suggesting that the downregulation of MRPS12 expression may play an important role in ischemic stroke. At last, PCR showed that MRPS11 and MRPS12 were significantly down-regulated in peripheral blood of patients with ischemic stroke.ConclusionsOur study provides a reference for the study of pathogenesis and key targets of ischemic stroke.
| S-EPMC10373980 | biostudies-literature
Project description:Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.
| S-EPMC10008864 | biostudies-literature
Project description:Background: Endometriosis is a common gynecological disorder with high rates of infertility and pelvic pain. However, its pathogenesis and diagnostic biomarkers remain unclear. This study aimed to elucidate potential hub genes and key pathways associated with endometriosis in ectopic endometrium (EC) and eutopic endometrium (EU). Material and Method: EC and EU-associated microarray datasets were obtained from the gene expression omnibus (GEO) database. Gene set enrichment analysis was performed to obtain further biological insight into the EU and EC-associated genes. Weighted gene co-expression network analysis (WGCNA) was performed to find clinically significant modules of highly-correlated genes. The hub genes that belong to both the weighted gene co-expression network and protein-protein interaction (PPI) network were identified using a Venn diagram. Results: We obtained EC and EU-associated microarray datasets GSE7305 and GSE120103. Genes in the EC were mainly enriched in the immune response and immune cell trafficking, and genes in the EU were mainly enriched in stress response and steroid hormone biosynthesis. PPI networks and weighted gene co-expression networks were constructed. An EC-associated blue module and an EU-associated magenta module were identified, and their function annotations revealed that hormone receptor signaling or inflammatory microenvironments may promote EU passing through the oviducts and migrating to the ovarian surfaces, and adhesion and immune correlated genes may induce the successful ectopic implantation of the endometrium (EC). Twelve hub genes in the EC and sixteen hub genes in the EU were recognized and further validated in independent datasets. Conclusion: Our study identified, for the first time, the hub genes and enrichment pathways in the EC and EU using WGCNA, which may provide a comprehensive understanding of the pathogenesis of endometriosis and have important clinical implications for the treatment and diagnosis of endometriosis.
| S-EPMC8436105 | biostudies-literature
Project description:BackgroundPsoriasis is a chronic inflammatory skin disease, the pathogenesis of which is more complicated and often requires long-term treatment. In particular, moderate to severe psoriasis usually requires systemic treatment. Psoriasis is also associated with many diseases, such as cardiometabolic diseases, malignant tumors, infections, and mood disorders. Psoriasis can appear at any age, and lead to a substantial burden for individuals and society. At present, psoriasis is still a treatable, but incurable, disease. Previous studies have found that microRNAs (miRNAs) play an important regulatory role in the progression of various diseases. Currently, miRNAs studies in psoriasis and dermatology are relatively new. Therefore, the identification of key miRNAs in psoriasis is helpful to elucidate the molecular mechanism of psoriasis.AimTo identify key molecular markers and signaling pathways to provide potential basis for the treatment and management of psoriasis.MethodsThe miRNA and mRNA data were obtained from the Gene Expression Omnibus database. Then, differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs) were screened out by limma R package. Subsequently, DEmRNAs were analyzed for Gene Ontology and Kyoto Encyclopedia of Genes and Genomics functional enrichment. The "WGCNA" R package was used to analyze the co-expression network of all miRNAs. In addition, we constructed miRNA-mRNA regulatory networks based on identified hub miRNAs. Finally, in vitro validation was performed. All experimental procedures were approved by the ethics committee of Chinese PLA General Hospital (S2021-012-01).ResultsA total of 639 DEmRNAs and 84 DEmiRNAs were identified. DEmRNAs screening criteria were adjusted P (adj. P) value < 0.01 and |logFoldChange| (|logFC|) > 1. DEmiRNAs screening criteria were adj. P value < 0.01 and |logFC| > 1.5. KEGG functional analysis demonstrated that DEmRNAs were significantly enriched in immune-related biological functions, for example, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction, and chemokine signaling pathway. In weighted gene co-expression network analysis, turquoise module was the hub module. Moreover, 10 hub miRNAs were identified. Among these 10 hub miRNAs, only 8 hub miRNAs predicted the corresponding target mRNAs. 97 negatively regulated miRNA-mRNA pairs were involved in the miRNA-mRNA regulatory network, for example, hsa-miR-21-5p-claudin 8 (CLDN8), hsa-miR-30a-3p-interleukin-1B (IL-1B), and hsa-miR-181a-5p/hsa-miR-30c-2-3p-C-X-C motif chemokine ligand 9 (CXCL9). Real-time polymerase chain reaction results showed that IL-1B and CXCL9 were up-regulated and CLDN8 was down-regulated in psoriasis with statistically significant differences.ConclusionThe identification of potential key molecular markers and signaling pathways provides potential research directions for further understanding the molecular mechanisms of psoriasis. This may also provide new research ideas for the prevention and treatment of psoriasis in the future.
| S-EPMC9254198 | biostudies-literature