scholarly journals Integrated weighted gene coexpression network analysis identifies Frizzled 2 (FZD2) as a key gene in invasive malignant pleomorphic adenoma

2022 ◽  
Vol 20 (1) ◽  
Author(s):  
Zhenyuan Han ◽  
Huiping Ren ◽  
Jingjing Sun ◽  
Lihui Jin ◽  
Qin Wang ◽  
...  

Abstract Background Invasive malignant pleomorphic adenoma (IMPA) is a highly malignant neoplasm of the oral salivary glands with a poor prognosis and a considerable risk of recurrence. Many disease-causing genes of IMPA have been identified in recent decades (e.g., P53, PCNA and HMGA2), but many of these genes remain to be explored. Weighted gene coexpression network analysis (WGCNA) is a newly emerged algorithm that can cluster genes and form modules based on similar gene expression patterns. This study constructed a gene coexpression network of IMPA via WGCNA and then carried out multifaceted analysis to identify novel disease-causing genes. Methods RNA sequencing (RNA-seq) was performed for 10 pairs of IMPA and normal tissues to acquire the gene expression profiles. Differentially expressed genes (DEGs) were screened out with the cutoff criteria of |log2 Fold change (FC)|> 1 and adjusted p value  < 0.05. Then, WGCNA was applied to systematically identify the hidden diagnostic hub genes of IMPA. Results In this research, a total of 1970 DEGs were screened out in IMPA tissues, including 1056 upregulated DEGs and 914 downregulated DEGs. Functional enrichment analysis was performed for identified DEGs and revealed an enrichment of tumor-associated GO terms and KEGG pathways. We used WGCNA to identify gene module most relevant with the histological grade of IMPA. The gene FZD2 was then recognized as the hub gene of the selected module with the highest module membership (MM) value and intramodule connectivity in protein–protein interaction (PPI) network. According to immunohistochemistry (IHC) staining, the expression level of FZD2 was higher in low-grade IMPA than in high-grade IMPA. Conclusion FZD2 shows an expression dynamic that is negatively correlated with the clinical malignancy of IMPA and it plays a central role in the transcription network of IMPA. Thus, FZD2 serves as a promising histological indicator for the precise prediction of IMPA histological stages.

2021 ◽  
Vol 2021 ◽  
pp. 1-22
Author(s):  
Aoran Yang ◽  
Xinhuan Wang ◽  
Yaofeng Hu ◽  
Chao Shang ◽  
Yang Hong

The function of glutamate ionotropic receptor NMDA type subunit 1 (GRIN1) in neurodegenerative diseases has been widely reported; however, its role in the occurrence of glioma remains less explored. We obtained clinical data and transcriptome data from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). Hub gene’s expression differential analysis and survival analysis were conducted by browsing the Gene Expression Profiling Interactive Analysis (GEPIA) database, Human Protein Atlas database, and LOGpc database. We conducted a variation analysis of datasets obtained from GEO and TCGA and performed a weighted gene coexpression network analysis (WGCNA) using the R programming language (3.6.3). Kaplan-Meier (KM) analysis was used to calculate the prognostic value of GRIN1. Finally, we conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Using STRING, we constructed a protein–protein interaction (PPI) network. Cytoscape software, a prerequisite of visualizing core genes, was installed, and CytoHubba detected the 10 most tumor-related core genes. We identified 185 differentially expressed genes (DEGs). GO and KEGG enrichment analyses illustrated that the identified DEGs are imperative in different biological functions and ascertained the potential pathways in which the DEGs may be enriched. The overall survival calculated by KM analysis showed that patients with lower expression of GRIN1 had worse prognoses than patients with higher expression of GRIN1 ( p = 0.004 ). The GEPIA and LOGpc databases were used to verify the expression difference of GRIN1 among GBM, LGG, and normal brain tissues. Ultimately, immunohistochemical assay results showed that GRIN1 was detected in normal tissue and not in the tumor specimens. Our results highlight a potential target for glioma treatment and will further our understanding of the molecular mechanisms underlying the treatment of glioma.


2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Mi Zhou ◽  
Ruru Guo ◽  
Yong-Fei Wang ◽  
Wanling Yang ◽  
Rongxiu Li ◽  
...  

Systemic juvenile idiopathic arthritis (sJIA) is a severe autoinflammatory disorder with a still not clearly defined molecular mechanism. To better understand the disease, we used scattered datasets from public domains and performed a weighted gene coexpression network analysis (WGCNA) to identify key modules and hub genes underlying sJIA pathogenesis. Two gene expression datasets, GSE7753 and GSE13501, were used to construct the WGCNA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to the genes and hub genes in the sJIA modules. Cytoscape was used to screen and visualize the hub genes. We further compared the hub genes with the genome-wide association study (GWAS) genes and used a consensus WGCNA to verify that our conclusions were conservative and reproducible across multiple independent datasets. A total of 5,414 genes were obtained for WGCNA, from which highly correlated genes were divided into 17 modules. The red module demonstrated the highest correlation with the sJIA module ( r = 0.8 , p = 3 e − 29 ), whereas the green-yellow module was found to be closely related to the non-sJIA module ( r = 0.62 , p = 1 e − 14 ). Functional enrichment analysis demonstrated that the red module was mostly enriched in the activation of immune responses, infection, nucleosomes, and erythrocytes, and the green-yellow module was mostly enriched in immune responses and inflammation. Additionally, the hub genes in the red module were highly enriched in erythrocyte differentiation, including ALAS2, AHSP, TRIM10, TRIM58, and KLF1. The hub genes from the green-yellow module were mainly associated with immune responses, as exemplified by the genes KLRB1, KLRF1, CD160, and KIRs. We identified sJIA-related modules and several hub genes that might be associated with the development of sJIA. Particularly, the modules may help understand the mechanisms of sJIA, and the hub genes may become biomarkers and therapeutic targets of sJIA in the future.


2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Biao Yang ◽  
Shuxun Wei ◽  
Yan-Bin Ma ◽  
Sheng-Hua Chu

Meningiomas are the most common primary intracranial tumor in adults. However, to date, systemic coexpression analyses for meningiomas fail to explain its pathogenesis. The aim of the present study was to construct coexpression modules and identify potential biomarkers associated with meningioma progression. Weighted gene coexpression network analysis (WGCNA) was performed based on GSE43290, and module preservation was tested by GSE74385. Functional annotations were performed to analyze biological significance. Hub genes were selected for efficacy evaluations and correlation analyses using two independent cohorts. A total of 14 coexpression modules were identified, and module lightcyan was significantly associated with WHO grades. Functional enrichment analyses of module lightcyan were associated with tumor pathogenesis. The top 10 hub genes were extracted. Ten biomarkers, particularly AHCYL2, FGL2, and KCNMA1, were significantly related to grades and prognosis of meningioma. These findings not only construct coexpression modules leading to the better understanding of its pathogenesis but also provide potential biomarkers that represent specific on tumor grades and identify recurrence, predicting prognosis and progression of meningiomas.


2019 ◽  
Vol 8 (8) ◽  
pp. 1160 ◽  
Author(s):  
Wang ◽  
Li ◽  
Cai ◽  
Sheu ◽  
Tsai ◽  
...  

Breast cancer is one of the most common malignancies. However, the molecular mechanisms underlying its pathogenesis remain to be elucidated. The present study aimed to identify the potential prognostic marker genes associated with the progression of breast cancer. Weighted gene coexpression network analysis was used to construct free-scale gene coexpression networks, evaluate the associations between the gene sets and clinical features, and identify candidate biomarkers. The gene expression profiles of GSE48213 were selected from the Gene Expression Omnibus database. RNA-seq data and clinical information on breast cancer from The Cancer Genome Atlas were used for validation. Four modules were identified from the gene coexpression network, one of which was found to be significantly associated with patient survival time. The expression status of 28 genes formed the black module (basal); 18 genes, dark red module (claudin-low); nine genes, brown module (luminal), and seven genes, midnight blue module (nonmalignant). These modules were clustered into two groups according to significant difference in survival time between the groups. Therefore, based on betweenness centrality, we identified TXN and ANXA2 in the nonmalignant module, TPM4 and LOXL2 in the luminal module, TPRN and ADCY6 in the claudin-low module, and TUBA1C and CMIP in the basal module as the genes with the highest betweenness, suggesting that they play a central role in information transfer in the network. In the present study, eight candidate biomarkers were identified for further basic and advanced understanding of the molecular pathogenesis of breast cancer by using co-expression network analysis.


2021 ◽  
Vol 2021 ◽  
pp. 1-12
Author(s):  
Wei Liu ◽  
Yongquan Shi ◽  
Tao Cheng ◽  
Ruixue Jia ◽  
Ming-Zhong Sun ◽  
...  

In mouse models, the recovery of liver volume is mainly mediated by the proliferation of hepatocytes after partial hepatectomy that is commonly accompanied with ischemia-reperfusion. The identification of differently expressed genes in liver following partial hepatectomy benefits the better understanding of the molecular mechanisms during liver regeneration (LR) with appliable clinical significance. Briefly, studying different gene expression patterns in liver tissues collected from the mice group that survived through extensive hepatectomy will be of huge critical importance in LR than those collected from the mice group that survived through appropriate hepatectomy. In this study, we performed the weighted gene coexpression network analysis (WGCNA) to address the central candidate genes and to construct the free-scale gene coexpression networks using the identified dynamic different expressive genes in liver specimens from the mice with 85% hepatectomy (20% for seven-day survial rate) and 50% hepatectomy (100% for seven-day survial rate under ischemia-reperfusion condition compared with the sham group control mice). The WGCNA combined with Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses pinpointed out the apparent distinguished importance of three gene expression modules: the blue module for apoptotic process, the turquoise module for lipid metabolism, and the green module for fatty acid metabolic process in LR following extensive hepatectomy. WGCNA analysis and protein-protein interaction (PPI) network construction highlighted FAM175B, OGT, and PDE3B were the potential three hub genes in the previously mentioned three modules. This work may help to provide new clues to the future fundamental study and treatment strategy for LR following liver injury and hepatectomy.


2021 ◽  
Vol 12 ◽  
Author(s):  
Zhen-Qing Zhang ◽  
Wei-Wei Wu ◽  
Jin-Dong Chen ◽  
Guang-Yin Zhang ◽  
Jing-Yu Lin ◽  
...  

Bipolar disorder (BD) is a major and highly heritable mental illness with severe psychosocial impairment, but its etiology and pathogenesis remains unclear. This study aimed to identify the essential pathways and genes involved in BD using weighted gene coexpression network analysis (WGCNA), a bioinformatic method studying the relationships between genes and phenotypes. Using two available BD gene expression datasets (GSE5388, GSE5389), we constructed a gene coexpression network and identified modules related to BD. The analyses of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways were performed to explore functional enrichment of the candidate modules. A protein-protein interaction (PPI) network was further constructed to identify the potential hub genes. Ten coexpression modules were identified from the top 5,000 genes in 77 samples and three modules were significantly associated with BD, which were involved in several biological processes (e.g., the actin filament-based process) and pathways (e.g., MAPK signaling). Four genes (NOTCH1, POMC, NGF, and DRD2) were identified as candidate hub genes by PPI analysis and CytoHubba. Finally, we carried out validation analyses in a separate dataset, GSE12649, and verified NOTCH1 as a hub gene and the involvement of several biological processes such as actin filament-based process and axon development. Taken together, our findings revealed several candidate pathways and genes (NOTCH1) in the pathogenesis of BD and call for further investigation for their potential research values in BD diagnosis and treatment.


2020 ◽  
Vol 2020 ◽  
pp. 1-20
Author(s):  
Xuanchen Zhou ◽  
Xiaoyue Zhen ◽  
Yiqing Liu ◽  
Zhaoyang Cui ◽  
Zhiyong Yue ◽  
...  

Chronic rhinosinusitis with nasal polyps (CRSwNP) is a chronic inflammatory disease with relatively easy recurrence. However, the precise molecular mechanisms of this disease are poorly known. Based on gene sequencing data obtained from the Gene Expression Omnibus (GEO) database, we constructed coexpression networks by weighted gene coexpression network analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by the Database for Annotation, Visualization, and Integrated Discovery (DAVID). The core gene of pathogenesis, CRSwNP, was screened by protein-protein interaction data (PPI) from the HPRD database. Unsupervised clustering was applied to screen hub genes related to the phenotype of CRSwNP. Blue and turquoise modules were found to be most significantly related to the pathogenicity of CRSwNP. Functional enrichment analysis showed that cell proliferation in the blue modules, the apoptotic process in the turquoise module, and the cancer pathway in both modules were mostly significantly correlated with the development of CRSwNP. The noncoding RNAs (long noncoding RNA and microRNA) and the top 10 core genes in each module were found to be associated with the pathogenesis of CRSwNP. A total of nine hub genes were identified to be related to the CRSwNP phenotype. By qRT-PCR analysis, AKT1, CDH1, PIK3R1, CBL, LRP1, MALAT1, and XIST were proven to be associated with the pathogenesis of CRSwNP. AGR2, FAM3D, PIP, DSE, and TMC were identified to be related to the CRSwNP phenotype. Further exploration of these genes will reveal more important information about the mechanisms of CRSwNP.


Sign in / Sign up

Export Citation Format

Share Document