scholarly journals Identification of key genes and pathways associated with cholangiocarcinoma development based on weighted gene correlation network analysis

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7968 ◽  
Author(s):  
Jingwei Liu ◽  
Weixin Liu ◽  
Hao Li ◽  
Qiuping Deng ◽  
Meiqi Yang ◽  
...  

Background As the most frequently occurred tumor in biliary tract, cholangiocarcinoma (CCA) is mainly characterized by its late diagnosis and poor outcome. It is therefore urgent to identify specific genes and pathways associated with its progression and prognosis. Materials and Methods The differentially expressed genes in The Cancer Genome Atlas were analyzed to build the co-expression network by Weighted gene co-expression network analysis (WGCNA). Gene ontology (GO) as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted for the selected genes. Module–clinical trait relationships were analyzed to explore the association with clinicopathological parameters. Log-rank tests and cox regression were used to identify the prognosis-related genes. Results The most related modules with CCA development were tan module containing 181 genes and salmon module with 148 genes. GO analysis suggested enrichment terms of digestion, hormone transport and secretion, epithelial cell proliferation, signal release, fibroblast activation, response to acid chemical, wnt, Nicotinamide adenine dinucleotide phosphate metabolism. KEGG analysis demonstrated 15 significantly altered pathways including glutathione metabolism, wnt, central carbon metabolism, mTOR, pancreatic secretion, protein digestion, axon guidance, retinol metabolism, insulin secretion, salivary secretion, fat digestion. Key genes of SOX2, KIT, PRSS56, WNT9A, SLC4A4, PRRG4, PANX2, PIR, RASSF8, MFSD4A, INS, RNF39, IL1R2, CST1, and PPP3CA might be potential prognostic markers for CCA, of which RNF39 and PRSS56 also showed significant correlation with clinical stage. Discussion Differentially expressed genes and key modules contributing to CCA development were identified by WGCNA. Our results offer novel insights into the characteristics in the etiology, prognosis, and treatment of CCA.

2021 ◽  
Vol 2021 ◽  
pp. 1-15
Author(s):  
Hao Guo ◽  
Jing Zhou ◽  
Yanjun Zhang ◽  
Zhi Wang ◽  
Likun Liu ◽  
...  

Background. Hypoxia closely relates to malignant progression and appears to be prognostic for outcome in hepatocellular carcinoma (HCC). Our research is aimed at mining the hypoxic-related genes (HRGs) and constructing a prognostic predictor (PP) model on clinical prognosis in HCC patients. Methods. RNA-sequencing data about HRGs and clinical data of patients with HCC were obtained from The Cancer Genome Atlas (TCGA) database portal. Differentially expressed HRGs between HCC and para-carcinoma tissue samples were obtained by applying the Wilcox analysis in R statistical software. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were used for gene functional enrichment analyses. Then, the patients who were asked to follow up for at least one month were enrolled in the following study. Cox proportional risk regression model was applied to obtain key HRGs which related to overall survival (OS) in HCC. PP was constructed and defined, and the accuracy of PP was validated by constructing the signature in a training set and validation set. Connectivity map (CMap) was used to find potential drugs, and gene set cancer analysis (GSCA) was also performed to explore the underlying molecular mechanisms. Results. Thirty-seven differentially expressed HRGs were obtained. It contained 28 upregulated and 9 downregulated genes. After the univariate Cox regression model analysis, we obtained 27 prognosis-related HRGs. Of these, 25 genes were risk factors for cancer, and 2 genes were protective factors. The PP was composed by 12 key genes (HDLBP, SAP30, PFKP, DPYSL4, SLC2A1, HMOX1, PGK1, ERO1A, LDHA, ENO2, SLC6A6, and TPI1). GSCA results showed the overall activity of these 12 key genes in 10 cancer-related pathways. Besides, CMap identified deferoxamine, crotamiton, talampicillin, and lycorine might have effects with HCC. Conclusions. This study firstly reported 12 prognostic HRGs and constructed the model of the PP. This comprehensive research of multiple databases helps us gain insight into the biological properties of HCC and provides deferoxamine, crotamiton, talampicillin, and lycorine as potential drugs to fight against HCC.


2021 ◽  
Vol 17 ◽  
Author(s):  
Hui Zhang ◽  
Qidong Liu ◽  
Xiaoru Sun ◽  
Yaru Xu ◽  
Yiling Fang ◽  
...  

Background: The pathophysiology of Alzheimer's disease (AD) is still not fully studied. Objective: This study aimed to explore the differently expressed key genes in AD and build a predictive model of diagnosis and treatment. Methods: Gene expression data of the entorhinal cortex of AD, asymptomatic AD, and control samples from the GEO database were analyzed to explore the relevant pathways and key genes in the progression of AD. Differentially expressed genes between AD and the other two groups in the module were selected to identify biological mechanisms in AD through KEGG and PPI network analysis in Metascape. Furthermore, genes with a high connectivity degree by PPI network analysis were selected to build a predictive model using different machine learning algorithms. Besides, model performance was tested with five-fold cross-validation to select the best fitting model. Results: A total of 20 co-expression gene clusters were identified after the network was constructed. Module 1 (in black) and module 2 (in royal blue) were most positively and negatively correlated with AD, respectively. Total 565 genes in module 1 and 215 genes in module 2, respectively, overlapped in two differentially expressed genes lists. They were enriched in the G protein-coupled receptor signaling pathway, immune-related processes, and so on. 11 genes were screened by using lasso logistic regression, and they were considered to play an important role in predicting AD samples. The model built by the support vector machine algorithm with 11 genes showed the best performance. Conclusion: This result shed light on the diagnosis and treatment of AD.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Yuchu Ye ◽  
Jingyi Wang ◽  
Faya Liang ◽  
Pan Song ◽  
Xiaoqing Yan ◽  
...  

Abstract Background The cause and underlying molecular mechanisms of head and neck squamous cell carcinoma (HNSCC) are unclear. Our study aims to identify the key genes associated with HNSCC and reveal potential biomarkers. Methods In this study, the expression profile dataset GSE83519 of the Gene Expression Omnibus database and the RNA sequencing dataset of HNSCC of The Cancer Genome Atlas were included for analysis. Sixteen differentially expressed genes were screened from these two datasets using R software. Gene Expression Profiling Interactive Analysis 2 (GEPIA2) was then adopted for survival analysis, and finally, three key genes related to the overall survival of HNSCC patients were identified. Furthermore, we verified these three genes using the Oncomine database and from real-time PCR and immunohistochemistry results from HNSCC tissues. Results The expression data of 44 samples from GSE83519 and 545 samples from TCGA-HNSC were collected. Using bioinformatics, the two databases were integrated, and 16 DEGs were screened out. Gene Ontology (GO) enrichment analysis showed that the biological functions of DEGs focused primarily on the apical plasma membrane and regulation of anoikis. Kyoto Encyclopedia of Genes and Genomes (KEGG) signalling pathway analysis showed that these DEGs were mainly involved in drug metabolism-cytochrome P450 and serotonergic synapses. Survival analysis identified three key genes, CEACAM5, CEACAM6 and CLCA4, that were closely related to HNSCC prognosis. The Oncomine database, qRT–PCR and IHC verified that all 3 key genes were downregulated in most HNSCC tissues compared to adjacent normal tissues. Conclusions This study indicates that integrated bioinformatics analyses play an important role in screening for differentially expressed genes and pathways in HNSCC, helping us better understand the biomarkers and molecular mechanism of HNSCC.


Author(s):  
Zhongxiao Lin ◽  
Min Wen ◽  
Enxing Yu ◽  
Xiao Lin ◽  
Hua Wang ◽  
...  

The tumor microenvironment (TME) plays an important role in the growth and invasion of glioma. This study aimed to analyze the composition of the immune microenvironment in glioma samples and analyze the important differentially expressed genes to identify novel immune-targeted therapy for glioma. We downloaded transcriptomic data of 669 glioma samples from The Cancer Genome Atlas database. CIBERSORT and ESTIMATE methods were used to calculate the proportion of tumor-infiltrating immune cells and ratio of immune and stromal components in the TME. The differentially expressed genes (DEGs) were screened by comparing the genes expressed by both stromal and immune cells. Annexin A1 (ANXA1) was determined to be an important prognostic indicator through the common overlap of univariate Cox regression analysis and protein–protein interaction network analysis. The proportion of tumor-infiltrating immune cells, calculated by CIBERSORT algorithm, had a significant difference in distribution among the high and low ANXA1 expression groups, indicating that ANXA1 could be an important immune marker of TME. Furthermore, ANXA1 level was positively correlated with the histopathological factors and negatively related to the survival of glioma patients based on the analysis of multiple databases. Finally, in vitro experiments verified that antagonizing ANXA1 expression promoted cell apoptosis and inhibited the invasion and migration capacities of glioma cells. Therefore, ANXA1 due to its immune-related functions, can be an important prognostic indicator and immune microenvironmental marker for gliomas. Further studies are warranted to confirm ANXA1 as a potential immunotherapeutic target for gliomas.


2021 ◽  
Author(s):  
Xiao-Yan Huang ◽  
Wen-Tao Qin ◽  
Qi-Sheng Su ◽  
Cheng-Cheng Qiu ◽  
Ruo-Chuan Liu ◽  
...  

Abstract Objective: Aim of this study was to identify the stemness-related genes in Pancreatic ductal adenocarcinoma (PDAC). Methods: The RNA-seq data of PADC were downloaded from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) databases. mRNA expression base-index (mRNAsi) and epigenetically regulated mRNAsi (EREG-mRNAsi) of PADC were evaluated. The mRNAsi-based gene sets in PADC were identified by Weighted gene co-expression network analysis (WGCNA). Functional Enrichment Analyses were performed with key genes. Kaplan–Meier survival analysis and the Cox proportional hazards model were used to evaluate the prognostic value of the key genes. Prognosis-associated hub genes were applied to establish nomograms. The receiver operating characteristic curves (ROC) and concordance index (C-index) were utilized to assess the discrimination and accuracy of the nomogram. Finally, these results were validated in the Gene Expression Omnibus (GEO) database and immunohistochemistry (IHC). Results: 36 key genes associate with mRNAsi were screened via WGCNA. Next a five-gene signature compromised TPX2, NCAPH, UBE2C, CCNB2, CEP55. Based on the expression of the signature, the PADC patients were classify patients into high- and low-risk groups. Cox regression analysis revealed that the high-risk group was significantly positive with overall survival (OS). Moreover, the nomogram has better sensitivity and specificity for predicting the OS. And the ROC, C-index indicated good performance of the prognostic signature in the TCGA and GEO dataset. Conclusion: Prognostic model associated with cancer stem cells (CSCs) reliably predict OS in PADC. this might be beneficial for the treatment of PADC patients.


BMC Neurology ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Yi Wang ◽  
Guogang Dai ◽  
Ling Jiang ◽  
Shichuan Liao ◽  
Jiao Xia

Abstract Background Although the pathology of sciatica has been studied extensively, the transcriptional changes in the peripheral blood caused by sciatica have not been characterized. This study aimed to characterize the peripheral blood transcriptomic signature for sciatica. Methods We used a microarray to identify differentially expressed genes in the peripheral blood of patients with sciatica compared with that of healthy controls, performed a functional analysis to reveal the peripheral blood transcriptomic signature for sciatica, and conducted a network analysis to identify key genes that contribute to the observed transcriptional changes. The expression levels of these key genes were assessed by qRT-PCR. Results We found that 153 genes were differentially expressed in the peripheral blood of patients with sciatica compared with that of healthy controls, and 131 and 22 of these were upregulated and downregulated, respectively. A functional analysis revealed that these differentially expressed genes (DEGs) were strongly enriched for the inflammatory response or immunity. The network analysis revealed that a group of genes, most of which are related to the inflammatory response, played a key role in the dysregulation of these DEGs. These key genes are Toll-like receptor 4, matrix metallopeptidase 9, myeloperoxidase, cathelicidin antimicrobial peptide, resistin and Toll-like receptor 5, and a qRT-PCR analysis validated the higher transcript levels of these key genes in the peripheral blood of patients with sciatica than in that of healthy controls. Conclusion We revealed inflammatory characteristics that serve as a peripheral blood transcriptomic signature for sciatica and identified genes that are essential for mRNA dysregulation in the peripheral blood of patients with sciatica.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Sheng Zheng ◽  
Zizhen Zhang ◽  
Ning Ding ◽  
Jiawei Sun ◽  
Yifeng Lin ◽  
...  

Abstract Introduction Angiogenesis is a key factor in promoting tumor growth, invasion and metastasis. In this study we aimed to investigate the prognostic value of angiogenesis-related genes (ARGs) in gastric cancer (GC). Methods mRNA sequencing data with clinical information of GC were downloaded from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. The differentially expressed ARGs between normal and tumor tissues were analyzed by limma package, and then prognosis‑associated genes were screened using Cox regression analysis. Nine angiogenesis genes were identified as crucially related to the overall survival (OS) of patients through least absolute shrinkage and selection operator (LASSO) regression. The prognostic model and corresponding nomograms were establish based on 9 ARGs and verified in in both TCGA and GEO GC cohorts respectively. Results Eighty-five differentially expressed ARGs and their enriched pathways were confirmed. Significant enrichment analysis revealed that ARGs-related signaling pathway genes were highly related to tumor angiogenesis development. Kaplan–Meier analysis revealed that patients in the high-risk group had worse OS rates compared with the low-risk group in training cohort and validation cohort. In addition, RS had a good prognostic effect on GC patients with different clinical features, especially those with advanced GC. Besides, the calibration curves verified fine concordance between the nomogram prediction model and actual observation. Conclusions We developed a nine gene signature related to the angiogenesis that can predict overall survival for GC. It’s assumed to be a valuable prognosis model with high efficiency, providing new perspectives in targeted therapy.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 1037.2-1038
Author(s):  
X. Sun ◽  
S. X. Zhang ◽  
S. Song ◽  
T. Kong ◽  
C. Zheng ◽  
...  

Background:Psoriasis is an immune-mediated, genetic disease manifesting in the skin or joints or both, and also has a strong genetic predisposition and autoimmune pathogenic traits1. The hallmark of psoriasis is sustained inflammation that leads to uncontrolled keratinocyte proliferation and dysfunctional differentiation. And it’s also a chronic relapsing disease, which often necessitates a long-term therapy2.Objectives:To investigate the molecular mechanisms of psoriasis and find the potential gene targets for diagnosis and treating psoriasis.Methods:Total 334 gene expression data of patients with psoriasis research (GSE13355 GSE14905 and GSE30999) were obtained from the Gene Expression Omnibus database. After data preprocessing and screening of differentially expressed genes (DEGs) by R software. Online toll Metascape3 was used to analyze Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs. Interactions of proteins encoded by DEGs were discovered by Protein-protein interaction network (PPI) using STRING online software. Cytoscape software was utilized to visualize PPI and the degree of each DEGs was obtained by analyzing the topological structure of the PPI network.Results:A total of 611 DEGs were found to be differentially expressed in psoriasis. GO analysis revealed that up-regulated DEGs were mostly associated with defense and response to external stimulus while down-regulated DEGs were mostly associated with metabolism and synthesis of lipids. KEGG enrichment analysis suggested they were mainly enriched in IL-17 signaling, Toll-like receptor signaling and PPAR signaling pathways, Cytokine-cytokine receptor interaction and lipid metabolism. In addition, top 9 key genes (CXCL10, OASL, IFIT1, IFIT3, RSAD2, MX1, OAS1, IFI44 and OAS2) were identified through Cytoscape.Conclusion:DEGs of psoriasis may play an essential role in disease development and may be potential pathogeneses of psoriasis.References:[1]Boehncke WH, Schon MP. Psoriasis. Lancet 2015;386(9997):983-94. doi: 10.1016/S0140-6736(14)61909-7 [published Online First: 2015/05/31].[2]Zhang YJ, Sun YZ, Gao XH, et al. Integrated bioinformatic analysis of differentially expressed genes and signaling pathways in plaque psoriasis. Mol Med Rep 2019;20(1):225-35. doi: 10.3892/mmr.2019.10241 [published Online First: 2019/05/23].[3]Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019;10(1):1523. doi: 10.1038/s41467-019-09234-6 [published Online First: 2019/04/05].Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared


2021 ◽  
Vol 20 ◽  
pp. 153303382098329
Author(s):  
Yujie Weng ◽  
Wei Liang ◽  
Yucheng Ji ◽  
Zhongxian Li ◽  
Rong Jia ◽  
...  

Human epidermal growth factor 2 (HER2)+ breast cancer is considered the most dangerous type of breast cancers. Herein, we used bioinformatics methods to identify potential key genes in HER2+ breast cancer to enable its diagnosis, treatment, and prognosis prediction. Datasets of HER2+ breast cancer and normal tissue samples retrieved from Gene Expression Omnibus and The Cancer Genome Atlas databases were subjected to analysis for differentially expressed genes using R software. The identified differentially expressed genes were subjected to gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses followed by construction of protein-protein interaction networks using the STRING database to identify key genes. The genes were further validated via survival and differential gene expression analyses. We identified 97 upregulated and 106 downregulated genes that were primarily associated with processes such as mitosis, protein kinase activity, cell cycle, and the p53 signaling pathway. Visualization of the protein-protein interaction network identified 10 key genes ( CCNA2, CDK1, CDC20, CCNB1, DLGAP5, AURKA, BUB1B, RRM2, TPX2, and MAD2L1), all of which were upregulated. Survival analysis using PROGgeneV2 showed that CDC20, CCNA2, DLGAP5, RRM2, and TPX2 are prognosis-related key genes in HER2+ breast cancer. A nomogram showed that high expression of RRM2, DLGAP5, and TPX2 was positively associated with the risk of death. TPX2, which has not previously been reported in HER2+ breast cancer, was associated with breast cancer development, progression, and prognosis and is therefore a potential key gene. It is hoped that this study can provide a new method for the diagnosis and treatment of HER2 + breast cancer.


Sign in / Sign up

Export Citation Format

Share Document