Screening Hub Genes of Hepatocellular Carcinoma Based on Public Databases

Tumor recurrence and metastasis often occur in HCC patients after surgery, and the prognosis is not optimistic. Hence, searching effective biomarkers for prognosis of is of great importance. Firstly, HCC-related data was acquired from the TCGA and GEO databases. Based on GEO data, 256 differentially expressed genes (DEGs) were obtained firstly. Subsequently, to clarify function of DEGs, clusterProfiler package was used to conduct functional enrichment analyses on DEGs. Protein-protein interaction (PPI) network analysis screened 20 key genes. The key genes were filtered via GEPIA database, by which 11 hub genes (F9, CYP3A4, ASPM, AURKA, CDC20, CDCA5, NCAP, PRC1, PTTG1, TOP2A, and KIFC1) were screened out. Then, univariate Cox analysis was applied to construct a prognostic model, followed by a prediction performance validation. With the risk score calculated by the model and common clinical features, univariate and multivariate analyses were carried out to assess whether the prognostic model could be used independently for prognostic prediction. In conclusion, the current study screened HCC prognostic gene signature based on public databases.

Download Full-text

Two Similar Signatures for Predicting the Prognosis and Immunotherapy Efficacy of Stomach Adenocarcinoma Patients

Frontiers in Cell and Developmental Biology ◽

10.3389/fcell.2021.704242 ◽

2021 ◽

Vol 9 ◽

Author(s):

Taohua Yue ◽

Shuai Zuo ◽

Jing Zhu ◽

Shihao Guo ◽

Zhihao Huang ◽

...

Keyword(s):

Gene Signature ◽

Functional Enrichment ◽

Venn Diagram ◽

High Morbidity ◽

Barr Virus ◽

Genes Expression ◽

Stomach Adenocarcinoma ◽

Key Genes ◽

Epstein Barr ◽

Cox Analysis

BackgroundGlobally, stomach adenocarcinoma (STAD)’s high morbidity and mortality should arouse our urgent attention. How long can STAD patients survive after surgery and whether novel immunotherapy is effective are questions that our clinicians cannot escape.MethodsVarious R packages, GSEA software, Metascape, STRING, Cytoscape, Venn diagram, TIMER2.0 website, TCGA, and GEO databases were used in our study.ResultsIn the TCGA and GEO, macrophage abundance of STAD tissues was significantly higher than that of adjacent tissues and was an independent prognostic factor, significantly related to the overall survival (OS) of STAD patients. Between the high- and low- macrophage abundance, we conducted differential expression, univariate and multivariate Cox analysis, and obtained 12 candidate genes, and finally constructed a 3-gene signature. Both low macrophage abundance group and group D had higher TMB and PD-L1 expression. Furthermore, top 5 common gene-mutated STAD tissues had lower macrophage abundance. Macrophage abundance and 3 key genes expression were also lower in the Epstein-Barr Virus (EBV) and HM-indel STAD subtypes and significantly correlated with the tumor microenvironment score. The functional enrichment and ssGSEA revealed 2 signatures were similar and closely related to BOQUEST_STEM_CELL_UP, including genes up-regulated in proliferative stromal stem cells. Hsa-miR-335-5p simultaneously regulated 3 key genes and significantly related to the expression of PD-L1, CD8A and PDCD1.Conclusionmacrophage abundance and 3-gene signature could simultaneously predict the OS and immunotherapy efficacy, and both 2 signatures had remarkable similarities. Hsa-miR-335-5p and BOQUEST_STEM_CELL_UP might be novel immunotherapy targets.

Download Full-text

Screen Key Genes Associated with Distraction-Induced Osteogenesis of Stem Cells Using Bioinformatics Methods

International Journal of Molecular Sciences ◽

10.3390/ijms22126505 ◽

2021 ◽

Vol 22 (12) ◽

pp. 6505

Author(s):

Jishizhan Chen ◽

Jia Hua ◽

Wenhui Song

Keyword(s):

Stem Cells ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Negative Control ◽

Gene Set Enrichment Analysis ◽

Venn Diagram ◽

Hub Genes ◽

Protein Protein Interaction ◽

Two Samples ◽

Key Genes

Applying mesenchymal stem cells (MSCs), together with the distraction osteogenesis (DO) process, displayed enhanced bone quality and shorter treatment periods. The DO guides the differentiation of MSCs by providing mechanical clues. However, the underlying key genes and pathways are largely unknown. The aim of this study was to screen and identify hub genes involved in distraction-induced osteogenesis of MSCs and potential molecular mechanisms. Material and Methods: The datasets were downloaded from the ArrayExpress database. Three samples of negative control and two samples subjected to 5% cyclic sinusoidal distraction at 0.25 Hz for 6 h were selected for screening differentially expressed genes (DEGs) and then analysed via bioinformatics methods. The Gene Ontology (GO) terms and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment were investigated. The protein–protein interaction (PPI) network was visualised through the Cytoscape software. Gene set enrichment analysis (GSEA) was conducted to verify the enrichment of a self-defined osteogenic gene sets collection and identify osteogenic hub genes. Results: Three hub genes (IL6, MMP2, and EP300) that were highly associated with distraction-induced osteogenesis of MSCs were identified via the Venn diagram. These hub genes could provide a new understanding of distraction-induced osteogenic differentiation of MSCs and serve as potential gene targets for optimising DO via targeted therapies.

Download Full-text

HEMGN and SLC2A1 might be potential diagnostic biomarkers of steroid-induced osteonecrosis of femoral head: study based on WGCNA and DEGs screening

BMC Musculoskeletal Disorders ◽

10.1186/s12891-021-03958-7 ◽

2021 ◽

Vol 22 (1) ◽

Author(s):

Zhixin Wu ◽

Yinxian Wen ◽

Guanlan Fan ◽

Hangyuan He ◽

Siqi Zhou ◽

...

Keyword(s):

Femoral Head ◽

Roc Curve ◽

Enrichment Analysis ◽

Diagnostic Value ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Overlapping Genes ◽

Hub Genes ◽

Diagnostic Biomarkers ◽

Key Genes

Abstract Background Steroid-induced osteonecrosis of the femoral head (SONFH) is a chronic and crippling bone disease. This study aims to reveal novel diagnostic biomarkers of SONFH. Methods The GSE123568 dataset based on peripheral blood samples from 10 healthy individuals and 30 SONFH patients was used for weighted gene co-expression network analysis (WGCNA) and differentially expressed genes (DEGs) screening. The genes in the module related to SONFH and the DEGs were extracted for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Genes with |gene significance| > 0.7 and |module membership| > 0.8 were selected as hub genes in modules. The DEGs with the degree of connectivity ≥5 were chosen as hub genes in DEGs. Subsequently, the overlapping genes of hub genes in modules and hub genes in DEGs were selected as key genes for SONFH. And then, the key genes were verified in another dataset, and the diagnostic value of key genes was evaluated by receiver operating characteristic (ROC) curve. Results Nine gene co-expression modules were constructed via WGCNA. The brown module with 1258 genes was most significantly correlated with SONFH and was identified as the key module for SONFH. The results of functional enrichment analysis showed that the genes in the key module were mainly enriched in the inflammatory response, apoptotic process and osteoclast differentiation. A total of 91 genes were identified as hub genes in the key module. Besides, 145 DEGs were identified by DEGs screening and 26 genes were identified as hub genes of DEGs. Overlapping genes of hub genes in the key module and hub genes in DEGs, including RHAG, RNF14, HEMGN, and SLC2A1, were further selected as key genes for SONFH. The diagnostic value of these key genes for SONFH was confirmed by ROC curve. The validation results of these key genes in GSE26316 dataset showed that only HEMGN and SLC2A1 were downregulated in the SONFH group, suggesting that they were more likely to be diagnostic biomarkers of SOFNH than RHAG and RNF14. Conclusions Our study identified that two key genes, HEMGN and SLC2A1, might be potential diagnostic biomarkers of SONFH.

Download Full-text

Prediction of Prognostic Biomarkers and Construction of an Autophagy Prognostic Model for Colorectal Cancer Using Bioinformatics

Technology in Cancer Research & Treatment ◽

10.1177/1533033820984177 ◽

2020 ◽

Vol 19 ◽

pp. 153303382098417

Author(s):

Ting-ting Liu ◽

Shu-min Liu

Keyword(s):

Colorectal Cancer ◽

Regression Analysis ◽

Prognostic Model ◽

Cox Regression ◽

Small Gtpase ◽

Functional Enrichment ◽

Clinical Correlation ◽

Cox Regression Analysis ◽

Prognostic Analysis ◽

Key Genes

Objective: The incidence of colorectal cancer is increasing every year, and autophagy may be related closely to the pathogenesis of colorectal cancer. Autophagy is a natural catabolic mechanism that allows the degradation of cellular components in eukaryotic cells. However, autophagy plays a dual role in tumorigenesis. It not only promotes normal cell survival and tumor growth but also induces cell death and suppresses tumors survival. In addition, the pathogenesis of various conditions, including inflammation, neurodegenerative diseases, or tumors, is associated with abnormal autophagy. The present work aimed to examine the significance of autophagy-related genes (ARGs) in prognosis prediction, to construct an autophagy prognostic model, and to identify independent prognostic factors for colorectal cancer (CRC). Methods: This study discovered a total of 36 ARGs in CRC cases using The Cancer Genome Atlas (TCGA) and Human Autophagy-dedicated (HADd) databases along with functional enrichment analysis. Then, an autophagy prognostic model was constructed using univariate Cox regression analysis, and the key prognostic genes were screened. Finally, independent prognostic markers were determined through independent prognostic analysis and clinical correlation analysis of key genes. Results: Of the 36 differentially expressed ARGs, 13 were related to prognosis, as determined by univariate Cox regression analysis. A total of 6 key genes were obtained by a multivariate Cox regression analysis. Independent prognostic values were shown by 3 genes, namely, microtubule-associated protein 1 light chain 3 (MAP1LC3C), small GTPase superfamily and Rab family (RAB7A), and WD-repeat domain phosphoinositide-interacting protein 2 (WIPI2) by independent prognostic analysis and clinical correlation. Conclusions: In this study, molecular bioinformatics technology was employed to determine and construct a prognostic model of autophagy for colon cancer patients, which revealed 3 autophagy-related features, namely, MAP1LC3C, WIPI2, and RAB7A.

Download Full-text

Nomogram Models Based on Gene Expression in Prediction of Breast Cancer Bone Metastasis

10.21203/rs.3.rs-1109324/v1 ◽

2021 ◽

Author(s):

Teng-di Fan ◽

Di-kai Bei ◽

Song-wei Li

Keyword(s):

Breast Cancer ◽

Gene Expression ◽

Bone Metastasis ◽

Gene Signature ◽

Functional Enrichment ◽

Training Set ◽

Hub Genes ◽

Kaplan Meier ◽

Patient Prognosis ◽

Kaplan Meier Plotter

Abstract Objective: To design a weighted co-expression network and build gene expression signature-based nomogram (GESBN) models for predicting the likelihood of bone metastasis in breast cancer (BC) patients. Methods: Dataset GSE124647 was used as a training set, and GSE14020 was taken as a validation set. In the training cohort, limma package in R was adopted to obtain differentially expressed genes (DEGs) between BC non-bone metastasis and bone metastasis patients, which were used for functional enrichment analysis. After weighted co-expression network analysis (WGCNA), univariate Cox regression and Kaplan-Meier plotter analyses were performed to screen potential prognosis-related genes. Then, GESBN models were constructed and evaluated. Further, the expression levels of genes in the models were explored in the training set, which was validated in GSE14020. Finally, the prognostic value of hub genes in BC was explored. Results: A total of 1858 DEGs were obtained. WGCNA result showed that the blue module was most significantly related to bone metastasis and prognosis. After survival analyses, GAJ1, SLC24A3, ITGBL1, and SLC44A1 were subjected to construct a GESBN model for overall survival. While GJA1, IGFBP6, MDFI, ITGFBI, ANXA2, and SLC24A3 were subjected to build a GESBN model for progression-free survival. Kaplan-Meier plotter and receiver operating characteristic analyses presented the reliable prediction ability of the models. Besides, GJA1, IGFBP6, ITGBL1, SLC44A1, and TGFBI expressions were significantly different between the two groups in GSE124647 and GSE14020. The hub genes had a significant impact on patient prognosis. Conclusion: Both the four-gene signature and six-gene signature could accurately predict patient prognosis, which may provide novel treatment insights for BC bone metastasis.

Download Full-text

Identification of key genes and biological processes contributing to colitis associated dysplasia in ulcerative colitis

PeerJ ◽

10.7717/peerj.11321 ◽

2021 ◽

Vol 9 ◽

pp. e11321

Author(s):

Di Zhang ◽

Pengguang Yan ◽

Taotao Han ◽

Xiaoyun Cheng ◽

Jingnan Li

Keyword(s):

Ulcerative Colitis ◽

Mitochondrial Dysfunction ◽

Single Gene ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Cell Junction ◽

Biological Processes ◽

Hub Genes ◽

Key Genes

Background Ulcerative colitis-associated colorectal cancer (UC-CRC) is a life-threatening complication of ulcerative colitis (UC). The mechanisms underlying UC-CRC remain to be elucidated. The purpose of this study was to explore the key genes and biological processes contributing to colitis-associated dysplasia (CAD) or carcinogenesis in UC via database mining, thus offering opportunities for early prediction and intervention of UC-CRC. Methods Microarray datasets (GSE47908 and GSE87466) were downloaded from Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) between groups of GSE47908 were identified using the “limma” R package. Weighted gene co-expression network analysis (WGCNA) based on DEGs between the CAD and control groups was conducted subsequently. Functional enrichment analysis was performed, and hub genes of selected modules were identified using the “clusterProfiler” R package. Single-gene gene set enrichment analysis (GSEA) was conducted to predict significant biological processes and pathways associated with the specified gene. Results Six functional modules were identified based on 4929 DEGs. Green and blue modules were selected because of their consistent correlation with UC and CAD, and the highest correlation coefficient with the progress of UC-associated carcinogenesis. Functional enrichment analysis revealed that genes of these two modules were significantly enriched in biological processes, including mitochondrial dysfunction, cell-cell junction, and immune responses. However, GSEA based on differential expression analysis between sporadic colorectal cancer (CRC) and normal controls from The Cancer Genome Atlas (TCGA) indicated that mitochondrial dysfunction may not be the major carcinogenic mechanism underlying sporadic CRC. Thirteen hub genes (SLC25A3, ACO2, AIFM1, ATP5A1, DLD, TFE3, UQCRC1, ADIPOR2, SLC35D1, TOR1AIP1, PRR5L, ATOX1, and DTX3) were identified. Their expression trends were validated in UC patients of GSE87466, and their potential carcinogenic effects in UC were supported by their known functions and other relevant studies reported in the literature. Single-gene GSEA indicated that biological processes and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to angiogenesis and immune response were positively correlated with the upregulation of TFE3, whereas those related to mitochondrial function and energy metabolism were negatively correlated with the upregulation of TFE3. Conclusions Using WGCNA, this study found two gene modules that were significantly correlated with CAD, of which 13 hub genes were identified as the potential key genes. The critical biological processes in which the genes of these two modules were significantly enriched include mitochondrial dysfunction, cell-cell junction, and immune responses. TFE3, a transcription factor related to mitochondrial function and cancers, may play a central role in UC-associated carcinogenesis.

Download Full-text

Weighted gene co-expression network analysis to define pivotal modules and genes in diabetic heart failure

Bioscience Reports ◽

10.1042/bsr20200507 ◽

2020 ◽

Vol 40 (7) ◽

Author(s):

Weiwei Liang ◽

FangFang Sun

Keyword(s):

Heart Failure ◽

Kegg Pathway ◽

Diabetic Heart ◽

Ppi Network ◽

Hub Genes ◽

Protein Protein Interaction ◽

Light Yellow ◽

Go Enrichment ◽

Key Genes ◽

Functional Analyses

Abstract This research was carried out to reveal specific hub genes involved in diabetic heart failure, as well as remarkable pathways that hub genes locate. The GSE26887 dataset from the GEO website was downloaded. The gene co-expression network was generated and central modules were analyzed to identify key genes using the WGCNA method. Functional analyses were conducted on genes of the clinical interest modules via Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene ontology (GO) enrichment, associated with protein–protein interaction (PPI) network construction in a sequence. Centrality parameters of the PPI network were determined using the CentiScape plugin in Cytoscape. Key genes, defined as genes in the ≥95% percentile of the degree distribution of significantly perturbed networks, were identified. Twenty gene co-expression modules were detected by WGCNA analysis. The module marked in light yellow exhibited the most significant association with diabetes (P=0.08). Genes involved in this module were primarily located in immune response, plasma membrane and receptor binding, as shown by the GO analysis. These genes were primarily assembled in endocytosis and phagosomes for KEGG pathway enrichment. Three key genes, STK39, HLA-DPB1 and RAB5C, which may be key genes for diabetic heart failure, were identified. To our knowledge, our study is the first to have constructed the co-expression network involved in diabetic heart failure using the WGCNA method. The results of the present study have provided better understanding the molecular mechanism of diabetic heart failure.

Download Full-text

Development and Validation of a Novel 11-Gene Prognostic Model for Serous Ovarian Carcinomas Based on Lipid Metabolism Expression Profile

International Journal of Molecular Sciences ◽

10.3390/ijms21239169 ◽

2020 ◽

Vol 21 (23) ◽

pp. 9169

Author(s):

Mingjun Zheng ◽

Heather Mullikin ◽

Anna Hester ◽

Bastian Czogalla ◽

Helene Heidegger ◽

...

Keyword(s):

Ovarian Cancer ◽

Lipid Metabolism ◽

Prognostic Model ◽

Molecular Subtypes ◽

Gene Signature ◽

Functional Enrichment ◽

Serous Ovarian Cancer ◽

Ovarian Carcinomas ◽

Comparable Amount ◽

Lipid Metabolism Genes

(1) Background: Biomarkers might play a significant role in predicting the clinical outcomes of patients with ovarian cancer. By analyzing lipid metabolism genes, future perspectives may be uncovered; (2) Methods: RNA-seq data for serous ovarian cancer were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. The non-negative matrix factorization package in programming language R was used to classify molecular subtypes of lipid metabolism genes and the limma package in R was performed for functional enrichment analysis. Through lasso regression, we constructed a multi-gene prognosis model; (3) Results: Two molecular subtypes were obtained and an 11-gene signature was constructed (PI3, RGS, ADORA3, CH25H, CCDC80, PTGER3, MATK, KLRB1, CCL19, CXCL9 and CXCL10). Our prognostic model shows a good independent prognostic ability in ovarian cancer. In a nomogram, the predictive efficiency was notably superior to that of traditional clinical features. Related to known models in ovarian cancer with a comparable amount of genes, ours has the highest concordance index; (4) Conclusions: We propose an 11-gene signature prognosis prediction model based on lipid metabolism genes in serous ovarian cancer.

Download Full-text

Co-expression network analysis identifies specific hub genes in association with developmental neuronal remodeling in Drosophila melanogaster

10.7287/peerj.preprints.27586v1 ◽

2019 ◽

Author(s):

Yunze Liu ◽

Xiaojie Sun ◽

Aijun Qu

Keyword(s):

Drosophila Melanogaster ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Hub Genes ◽

Protein Protein Interaction ◽

Neuronal Remodeling ◽

Gene Modules

As an evolutionarily conserved mechanism, developmental neuronal remodeling is needed for the proper wiring of the nervous system and is critical for understanding the neurodevelopment mechanisms. Previous studies have shown that during metamorphosis lots of Drosophila melanogaster mushroom body neurons experience stereotypic remodeling. However, the related regulators and downstream executors of pathways are yet unclear, especially studies of transcriptional gene co-expression analysis of nervous systems remain insufficient. In this study, we develop a weighted gene co-expression network (WGCNA) to classify gene modules associated with neuronal remodeling. Moreover, functional and pathway enrichment analysis with protein-protein network construction is applied to detect high informative hub genes in the targeted gene module. Thus, we select a total of five hub genes that play critical roles in neuronal remodeling and identify them with functional enrichment analysis and protein-protein interaction network. Overall, this study provides insight into the underlying molecular mechanism of developmental neuronal remodeling in Drosophila melanogaster.

Download Full-text

Identification of Key Biomarkers and Immune Infiltration in Sciatic Nerve of Diabetic Neuropathy BKS-db/db Mice by Bioinformatics Analysis

Frontiers in Pharmacology ◽

10.3389/fphar.2021.682005 ◽

2021 ◽

Vol 12 ◽

Author(s):

Yixuan Lin ◽

Fanjing Wang ◽

Lianzhi Cheng ◽

Zhaohui Fang ◽

Guoming Shen

Keyword(s):

Diabetic Neuropathy ◽

Signaling Pathway ◽

Functional Enrichment ◽

Hub Genes ◽

Data Set ◽

Immune Infiltration ◽

P53 Signaling ◽

Ppi Networks ◽

Key Genes ◽

P53 Signaling Pathway

Diabetic neuropathy (DN) is one of the chronic complications of diabetes which can cause severe harm to patients. In order to determine the key genes and pathways related to the pathogenesis of DN, we downloaded the microarray data set GSE27382 from Gene Expression Omnibus (GEO) and adopted bioinformatics methods for comprehensive analysis, including functional enrichment, construction of PPI networks, central genes screening, TFs-target interaction analysis, and evaluation of immune infiltration characteristics. Finally, we examined quantitative real- time PCR (qPCR) to validate the expression of hub genes. A total of 318 differentially expressed genes (DEGs) were identified, among which 125 upregulated DEGs were enriched in the mitotic nuclear division, extracellular region, immunoglobulin receptor binding, and p53 signaling pathway, while 193 downregulated DEGs were enriched in ion transport, membrane, synapse, sodium channel activity, and retrograde endocannabinoid signaling. GSEA plots showed that condensed nuclear chromosome kinetochore were the most significant enriched gene set positively correlated with the DN group. Importantly, we identified five central genes (Birc5, Bub1, Cdk1, Ccnb2, and Ccnb1), and KEGG pathway analysis showed that the five hub genes were focused on progesterone-mediated oocyte maturation, cell cycle, and p53 signaling pathway. The proportion of immune cells from DN tissue and normal group showed significant individual differences. In DN samples, T cells CD4 memory resting and dendritic cells resting accounted for a higher proportion, and macrophage M2 accounted for a lower proportion. In addition, all five central genes showed consistent correlation with immune cell infiltration levels. qPCR showed the same expression trend of five central genes as in our analysis. Our research identified key genes related to differential genes and immune infiltration related to the pathogenesis of DN and provided new diagnostic and potential therapeutic targets for DN.

Download Full-text