scholarly journals Identification of 5 Hub Genes Related to the Early Diagnosis, Tumour Stage, and Poor Outcomes of Hepatitis B Virus-Related Hepatocellular Carcinoma by Bioinformatics Analysis

2021 ◽  
Vol 2021 ◽  
pp. 1-20
Author(s):  
Rui Qiang ◽  
Zitong Zhao ◽  
Lu Tang ◽  
Qian Wang ◽  
Yanhong Wang ◽  
...  

Background. The majority of primary liver cancers in adults worldwide are hepatocellular carcinomas (HCCs, or hepatomas). Thus, a deep understanding of the underlying mechanisms for the pathogenesis and carcinogenesis of HCC at the molecular level could facilitate the development of novel early diagnostic and therapeutic treatments to improve the approaches and prognosis for HCC patients. Our study elucidates the underlying molecular mechanisms of HBV-HCC development and progression and identifies important genes related to the early diagnosis, tumour stage, and poor outcomes of HCC. Methods. GSE55092 and GSE121248 gene expression profiling data were downloaded from the Gene Expression Omnibus (GEO) database. There were 119 HCC samples and 128 nontumour tissue samples. GEO2R was used to screen for differentially expressed genes (DEGs). Volcano plots and Venn diagrams were drawn by using the ggplot2 package in R. A heat map was generated by using Heatmapper. By using the clusterProfiler R package, KEGG and GO enrichment analyses of DEGs were conducted. Through PPI network construction using the STRING database, key hub genes were identified by cytoHubba. Finally, KM survival curves and ROC curves were generated to validate hub gene expression. Results. By GO enrichment analysis, 694 DEGs were enriched in the following GO terms: organic acid catabolic process, carboxylic acid catabolic process, carboxylic acid biosynthetic process, collagen-containing extracellular matrix, blood microparticle, condensed chromosome kinetochore, arachidonic acid epoxygenase activity, arachidonic acid monooxygenase activity, and monooxygenase activity. In the KEGG pathway enrichment analysis, DEGs were enriched in arachidonic acid epoxygenase activity, arachidonic acid monooxygenase activity, and monooxygenase activity. By PPI network construction and analysis of hub genes, we selected the top 10 genes, including CDK1, CCNB2, CDC20, BUB1, BUB1B, CCNB1, NDC80, CENPF, MAD2L1, and NUF2. By using TCGA and THPA databases, we found five genes, CDK1, CDC20, CCNB1, CENPF, and MAD2L1, that were related to the early diagnosis, tumour stage, and poor outcomes of HBV-HCC. Conclusions. Five abnormally expressed hub genes of HBV-HCC are informative for early diagnosis, tumour stage determination, and poor outcome prediction.

2021 ◽  
Author(s):  
Daoquan Liu ◽  
Jianhong Ma ◽  
Sheng Wei ◽  
Jianmin Liu ◽  
Mingzhou Li ◽  
...  

Abstract Background: Bladder cancer (BLCA) is the most popular malignant carcinomas in genitourinary system which has a high incidence and is prone to relapse. However, the molecular mechanism of BLCA remains to be unclear. Moreover, there is still a shortage of effective biomarkers that can predict progression and prognosis of BLCA. The objective of current study is to screen significant genes as biomarkers to forecast the progression and prognosis of BLCA patients.Methods: Gene expression profile downloaded from TCGA database and GEO database was used. Differential gene expression analysis and WGCNA were conducted to identify differential co-expression genes. In addition, GO enrichment analysis and KEGG pathway analysis were used to explore the functions of these genes. Moreover, PPI network, OS and DFS were used to identify survival-related hub genes. Finally, the expression levels of these genes were validated by qRT-PCR and HPA database.Results: About 124 differential co-expression genes were identified. And these genes were mainly enriched in muscle system process and muscle contraction (BP), contractile fiber, myofibril, sarcomere, focal adhesion and cell-substrate junction (CC) and actin binding (MF) in GO enrichment analysis, while enriched in vascular smooth muscle contraction, focal adhesion, cardiac muscle contraction, hypertrophic cardiomyopathy, dilated cardiomyopathy and regulation of actin cytoskeleton in KEGG analysis. Furthermore, five survival-related hub genes (MYH11, ACTA2, CALD1, TPM1, MYLK) were identified via overall OS and DFS. In addition, the expression levels of the five survival-related genes were upregulated with the procession of BLCA, such as grade, stage and TNM stage. Finally, all survival-related hub genes were found to be down-regulated in BLCA via qRT-PCR and HPA database.Conclusions: Our current study verified five new key genes in BLCA, which could help us better understand the pathogenesis of BLCA. And these five hub genes may be involved in the development and progression of BLCA and served as potential biomarkers.


2020 ◽  
Author(s):  
Basavaraj Vastrad ◽  
Chanabasayya Vastrad ◽  
Iranna Kotturshetti

AbstractTriple receptor negative breast cancer (TNBC) is the type of gynecological cancer in the elderly women. This study is aimed to explore molecular mechanism of TNBC via bioinformatics analysis. The gene expression profiles of GSE88715 (including 38 TNBC and 38 normal control) was downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were screened using the limma package in R software. Pathway and gene ontology (GO) enrichment analysis were performed based on various pathway dabases and GO database. Then, InnateDb interactome database, Cytoscape and PEWCC1 were applied to construct the protein-protein interaction (PPI) network and screen hub genes. Similarly, miRNet database, NetworkAnalyst database and Cytoscape were applied to construct the target gene - miRNA network and target gene - TF network, and screen targate genes. Pathway and GO enrichment analysis was further performed for hub genes, gene clusters identified via module analysis and targate genes. The expression of hub genes with prognostic values was validated on the UALCAN, cBio Portal, The Human Protein Atlas, receiver operator characteristic (ROC) curve analysis, RT-PCR analysis and immune infiltration analysis. A total of 949 DEGs were identified in TNBC (469 up regulated genes, and 480 down regulated genes), and they were mainly enriched in the terms of phospholipases, toxoplasmosis, immune response, cell surface, glycolysis, biosynthesis of amino acids, carboxylic acid metabolic process and organic substance catabolic process extracellular space. Hub genes including UBD, HLA-B, MYC and HSP90AB1 were identified via PPI network and modules, which were mainly enriched in immune response, antigen processing and presentation, cell cycle and pathways in cancer. Targate genes including CCDC80, PEG10, HOPX and CCNA2 were identified via target gene - miRNA network and target gene - TF network, which were mainly enriched in extracellular structure organization, validated targets of C-MYC transcriptional activation, ensemble of genes encoding core extracellular matrix including ECM glycoproteins and cell cycle. The top five significantly overexpressed mRNA (ADAM15, BATF, NOTCH3, ITGAX and SDC1) and the top five significantly underexpressed mRNA (RPL4, EEF1G, RPL3, RBMX and ABCC2) were selected for further validation in TNBCpatients and healthy controls. Analysis of the expression of genes in the various databases showed that ADAM15, BATF, NOTCH3, ITGAX, SDC1, RPL4, EEF1G, RPL3, RBMX and ABCC2 expressions have a cancer specific pattern in TNBC. Collectively, ADAM15, BATF, NOTCH3, ITGAX, SDC1, RPL4, EEF1G, RPL3, RBMX and ABCC2 may be useful candidate biomarkers for TNBC diagnosis, prognosis and theraputic targates.


Author(s):  
Xitong Yang ◽  
Pengyu Wang ◽  
Shanquan Yan ◽  
Guangming Wang

AbstractStroke is a sudden cerebrovascular circulatory disorder with high morbidity, disability, mortality, and recurrence rate, but its pathogenesis and key genes are still unclear. In this study, bioinformatics was used to deeply analyze the pathogenesis of stroke and related key genes, so as to study the potential pathogenesis of stroke and provide guidance for clinical treatment. Gene Expression profiles of GSE58294 and GSE16561 were obtained from Gene Expression Omnibus (GEO), the differentially expressed genes (DEGs) were identified between IS and normal control group. The different expression genes (DEGs) between IS and normal control group were screened with the GEO2R online tool. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were performed. Using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and gene set enrichment analysis (GSEA), the function and pathway enrichment analysis of DEGS were performed. Then, a protein–protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. Cytoscape with CytoHubba were used to identify the hub genes. Finally, NetworkAnalyst was used to construct the targeted microRNAs (miRNAs) of the hub genes. A total of 85 DEGs were screened out in this study, including 65 upward genes and 20 downward genes. In addition, 3 KEGG pathways, cytokine − cytokine receptor interaction, hematopoietic cell lineage, B cell receptor signaling pathway, were significantly enriched using a database for labeling, visualization, and synthetic discovery. In combination with the results of the PPI network and CytoHubba, 10 hub genes including CEACAM8, CD19, MMP9, ARG1, CKAP4, CCR7, MGAM, CD79A, CD79B, and CLEC4D were selected. Combined with DEG-miRNAs visualization, 5 miRNAs, including hsa-mir-146a-5p, hsa-mir-7-5p, hsa-mir-335-5p, and hsa-mir-27a- 3p, were predicted as possibly the key miRNAs. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of ischemic stroke, and provide a new strategy for clinical therapy.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Ruoting Lin ◽  
Conor E. Fogarty ◽  
Bowei Ma ◽  
Hejie Li ◽  
Guoying Ni ◽  
...  

Abstract Background Papillary thyroid carcinoma (PTC) is the most common thyroid cancer. While many patients survive, a portion of PTC cases display high aggressiveness and even develop into refractory differentiated thyroid carcinoma. This may be alleviated by developing a novel model to predict the risk of recurrence. Ferroptosis is an iron-dependent form of regulated cell death (RCD) driven by lethal accumulation of lipid peroxides, is regulated by a set of genes and shows a variety of metabolic changes. To elucidate whether ferroptosis occurs in PTC, we analyse the gene expression profiles of the disease and established a new model for the correlation. Methods The thyroid carcinoma (THCA) datasets were downloaded from The Cancer Genome Atlas (TCGA), UCSC Xena and MisgDB, and included 502 tumour samples and 56 normal samples. A total of 60 ferroptosis related genes were summarised from MisgDB database. Gene set enrichment analysis (GSEA) and Gene set variation analysis (GSVA) were used to analyse pathways potentially involving PTC subtypes. Single sample GSEA (ssGSEA) algorithm was used to analyse the proportion of 28 types of immune cells in the tumour immune infiltration microenvironment in THCA and the hclust algorithm was used to conduct immune typing according to the proportion of immune cells. Spearman correlation analysis was performed on the ferroptosis gene expression and the correlation between immune infiltrating cells proportion. We established the WGCNA to identify genes modules that are highly correlated with the microenvironment of immune invasion. DEseq2 algorithm was further used for differential analysis of sequencing data to analyse the functions and pathways potentially involving hub genes. GO and KEGG enrichment analysis was performed using Clusterprofiler to explore the clinical efficacy of hub genes. Univariate Cox analysis was performed for hub genes combined with clinical prognostic data, and the results was included for lasso regression and constructed the risk regression model. ROC curve and survival curve were used for evaluating the model. Univariate Cox analysis and multivariate Cox analysis were performed in combination with the clinical data of THCA and the risk score value, the clinical efficacy of the model was further evaluated. Results We identify two subtypes in PTC based on the expression of ferroptosis related genes, with the proportion of cluster 1 significantly higher than cluster 2 in ferroptosis signature genes that are positively associated. The mutations of Braf and Nras are detected as the major mutations of cluster 1 and 2, respectively. Subsequent analyses of TME immune cells infiltration indicated cluster 1 is remarkably richer than cluster 2. The risk score of THCA is in good performance evaluated by ROC curve and survival curve, in conjunction with univariate Cox analysis and multivariate Cox analysis results based on the clinical data shows that the risk score of the proposed model could be used as an independent prognostic indicator to predict the prognosis of patients with papillary thyroid cancer. Conclusions Our study finds seven crucial genes, including Ac008063.2, Apoe, Bcl3, Acap3, Alox5ap, Atxn2l and B2m, and regulation of apoptosis by parathyroid hormone-related proteins significantly associated with ferroptosis and immune cells in PTC, and we construct the risk score model which can be used as an independent prognostic index to predict the prognosis of patients with PTC.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Baojie Wu ◽  
Shuyi Xi

Abstract Background This study aimed to explore and identify key genes and signaling pathways that contribute to the progression of cervical cancer to improve prognosis. Methods Three gene expression profiles (GSE63514, GSE64217 and GSE138080) were screened and downloaded from the Gene Expression Omnibus database (GEO). Differentially expressed genes (DEGs) were screened using the GEO2R and Venn diagram tools. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Gene set enrichment analysis (GSEA) was performed to analyze the three gene expression profiles. Moreover, a protein–protein interaction (PPI) network of the DEGs was constructed, and functional enrichment analysis was performed. On this basis, hub genes from critical PPI subnetworks were explored with Cytoscape software. The expression of these genes in tumors was verified, and survival analysis of potential prognostic genes from critical subnetworks was conducted. Functional annotation, multiple gene comparison and dimensionality reduction in candidate genes indicated the clinical significance of potential targets. Results A total of 476 DEGs were screened: 253 upregulated genes and 223 downregulated genes. DEGs were enriched in 22 biological processes, 16 cellular components and 9 molecular functions in precancerous lesions and cervical cancer. DEGs were mainly enriched in 10 KEGG pathways. Through intersection analysis and data mining, 3 key KEGG pathways and related core genes were revealed by GSEA. Moreover, a PPI network of 476 DEGs was constructed, hub genes from 12 critical subnetworks were explored, and a total of 14 potential molecular targets were obtained. Conclusions These findings promote the understanding of the molecular mechanism of and clinically related molecular targets for cervical cancer.


2021 ◽  
Author(s):  
Li Guoquan ◽  
Du Junwei ◽  
He Qi ◽  
Fu Xinghao ◽  
Ji Feihong ◽  
...  

Abstract BackgroundHashimoto's thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is a common autoimmune disease, which mainly occurs in women. The early manifestation was hyperthyroidism, however, hypothyroidism may occur if HT was not controlled for a long time. Numerous studies have shown that multiple factors, including genetic, environmental, and autoimmune factors, were involved in the pathogenesis of the disease, but the exact mechanisms were not yet clear. The aim of this study was to identify differentially expressed genes (DEGs) by comprehensive analysis and to provide specific insights into HT. MethodsTwo gene expression profiles (GSE6339, GSE138198) about HT were downloaded from the Gene Expression Omnibus (GEO) database. The DEGs were assessed between the HT and normal groups using the GEO2R. The DEGs were then sent to the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The hub genes were discovered using Cytoscape and CytoHubba. Finally, NetworkAnalyst was utilized to create the hub genes' targeted microRNAs (miRNAs). ResultsA total of 62 DEGs were discovered, including 60 up-regulated and 2 down-regulated DEGs. The signaling pathways were mainly engaged in cytokine interaction and cytotoxicity, and the DEGs were mostly enriched in immunological and inflammatory responses. IL2RA, CXCL9, IL10RA, CCL3, CCL4, CCL2, STAT1, CD4, CSF1R, and ITGAX were chosen as hub genes based on the results of the protein-protein interaction (PPI) network and CytoHubba. Five miRNAs, including mir-24-3p, mir-223-3p, mir-155-5p, mir-34a-5p, mir-26b-5p, and mir-6499-3p, were suggested as likely important miRNAs in HT. ConclusionsThese hub genes, pathways and miRNAs contribute to a better understanding of the pathophysiology of HT and offer potential treatment options for HT.


2020 ◽  
Author(s):  
Vijayakrishna Kolur ◽  
Basavaraj Vastrad ◽  
Chanabasayya Vastrad ◽  
Iranna Kotturshetti ◽  
Anandkumar Tengli

Abstract BackgroundCoronary artery disease (CAD) is one of the most common disorders in the cardiovascular system. This study aims to explore potential signaling pathways and important biomarkers that drive CAD development. MethodsThe CAD GEO Dataset GSE113079 was featured to screen differentially expressed genes (DEGs). The pathway and Gene Ontology (GO) enrichment analysis of DEGs were analyzed using the ToppGene. We screened hub and target genes from protein-protein interaction (PPI) networks, target gene - miRNA regulatory network and target gene - TF regulatory network, and Cytoscape software. Validations of hub genes were performed to evaluate their potential prognostic and diagnostic value for CAD. Results1,036 DEGs were captured according to screening criteria (525upregulated genes and 511downregulated genes). Pathway and Gene Ontology (GO) enrichment analysis of DEGs revealed that these up and down regulated genes are mainly enriched in thyronamine and iodothyronamine metabolism, cytokine-cytokine receptor interaction, nervous system process, cell cycle and nuclear membrane. Hub genes were validated to find out potential prognostic biomarkers, diagnostic biomarkers and novel therapeutic target for CAD. ConclusionsIn summary, our findings discovered pivotal gene expression signatures and signaling pathways in the progression of CAD. CAPN13, ACTBL2, ERBB3, GATA4, GNB4, NOTCH2, EXOSC10, RNF2, PSMA1 and PRKAA1 might contribute to the progression of CAD, which could have potential as biomarkers or therapeutic targets for CAD.


2020 ◽  
Author(s):  
Xi Pan ◽  
Jian-Hao Liu

Abstract Background Nasopharyngeal carcinoma (NPC) is a heterogeneous carcinoma that the underlying molecular mechanisms involved in the tumor initiation, progression, and migration are largely unclear. The purpose of the present study was to identify key biomarkers and small-molecule drugs for NPC screening, diagnosis, and therapy via gene expression profile analysis. Methods Raw microarray data of NPC were retrieved from the Gene Expression Omnibus (GEO) database and analyzed to screen out the potential differentially expressed genes (DEGs). The key modules associated with histology grade and tumor stage was identified by using weighted correlation network analysis (WGCNA). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of genes in the key module were performed to identify potential mechanisms. Candidate hub genes were obtained, which based on the criteria of module membership (MM) and high connectivity. Then we used receiver operating characteristic (ROC) curve to evaluate the diagnostic value of hub genes. The Connectivity map database was further used to screen out small-molecule drugs of hub genes. Results A total of 430 DEGs were identified based on two GEO datasets. The green gene module was considered as key module for the tumor stage of NPC via WGCNA analysis. The results of functional enrichment analysis revealed that genes in the green module were enriched in regulation of cell cycle, p53 signaling pathway, cell part morphogenesis. Furthermore, four DEGs-related hub genes in the green module were considered as the final hub genes. Then ROC revealed that the final four hub genes presented with high areas under the curve, suggesting these hub genes may be diagnostic biomarkers for NPC. Meanwhile, we screened out several small-molecule drugs that have provided potentially therapeutic goals for NPC. Conclusions Our research identified four potential prognostic biomarkers and several candidate small-molecule drugs for NPC, which may contribute to the new insights for NPC therapy.


2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Ming Chen ◽  
Junkai Zeng ◽  
Yeqing Yang ◽  
Buling Wu

Abstract Background Pulpitis is an inflammatory disease, the grade of which is classified according to the level of inflammation. Traditional methods of evaluating the status of dental pulp tissue in clinical practice have limitations. The rapid and accurate diagnosis of pulpitis is essential for determining the appropriate treatment. By integrating different datasets from the Gene Expression Omnibus (GEO) database, we analysed a merged expression matrix of pulpitis, aiming to identify biological pathways and diagnostic biomarkers of pulpitis. Methods By integrating two datasets (GSE77459 and GSE92681) in the GEO database using the sva and limma packages of R, differentially expressed genes (DEGs) of pulpitis were identified. Then, the DEGs were analysed to identify biological pathways of dental pulp inflammation with Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis and Gene Set Enrichment Analysis (GSEA). Protein–protein interaction (PPI) networks and modules were constructed to identify hub genes with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape. Results A total of 470 DEGs comprising 394 upregulated and 76 downregulated genes were found in pulpitis tissue. GO analysis revealed that the DEGs were enriched in biological processes related to inflammation, and the enriched pathways in the KEGG pathway analysis were cytokine-cytokine receptor interaction, chemokine signalling pathway and NF-κB signalling pathway. The GSEA results provided further functional annotations, including complement system, IL6/JAK/STAT3 signalling pathway and inflammatory response pathways. According to the degrees of nodes in the PPI network, 10 hub genes were identified, and 8 diagnostic biomarker candidates were screened: PTPRC, CD86, CCL2, IL6, TLR8, MMP9, CXCL8 and ICAM1. Conclusions With bioinformatics analysis of merged datasets, biomarker candidates of pulpitis were screened and the findings may be as reference to develop a new method of pulpitis diagnosis.


Sign in / Sign up

Export Citation Format

Share Document