Protein Network Analysis to Prioritize Key Genes and Pathway for Stress-Mediated Neurodegeneration

Background:Oxidative Stress (OS) has been implicated in the pathophysiology of many neurodegenerative diseases. OS can cause cellular damage that results in cell death due to overproduction of reactive oxygen species (ROS) that may play the crucial role in the disease progression. An impaired mechanism in correlation with reduced expression of antioxidant proteins is the very common feature among most of the age-related disorders. Variousin-vitroandin-vivostudies suggest the major contribution of oxidative stress in neurodegeneration. Role of Nrf2 gene is well established as a neuroprotective gene especially in concern with stress-mediated neurodegeneration. Nrf2 is a bZIP transcription factor that forms the heterodimer with small Maf protein and transcription factor AP1 that regulates transcription by binding to ARE which coordinates the transcription of genes involved in phase II detoxification and an antioxidant defense that is used to protect the cell from oxidative stress.Aim:The currentinsilicostudy was attempted to prioritize key genes and pathway in stress-mediated neurodegeneration through network-based analysis.Methods:Protein-protein interaction network was constructed and analyzed using 63 Nrf2 regulating candidate genes obtained from NCBI database based on literature studies usingSTRING 10.0database andCytoscape v 3.6.0software plug-inNetwork Analyzer.Further, the functional enrichment analysis of identified gene was done usingPANTHER GENE ONTOLOGYsoftware and DAVID tool.Results:Based on network topological parameter, TP53, JUN, MYC, NFE2L2, AKT1, PIK3CA & UBC were identified as the key gene in the network. Among them, TP53 gene was obtained as a super hub gene with the highest Betweenness Centrality (BC) and node degree. The functional enrichment analysis was done usingPANTHER GENE ONTOLOGYsoftware and DAVID tool reveals their significant role in neurotrophin signaling pathway, MAPK signaling pathway, cellular response to stress & in the regulation of stress.Conclusion:The network analysis will help in prioritizing genes in the pathway that helps in understanding the underlying mechanism of disease. Thus, further study on these genes and their biological mechanism and pathway may, therefore, provide a potential target for the treatment of stress-mediated neurodegeneration.

Download Full-text

Identification and Validation of Immune-Related Prognostic Genes in the Tumor Microenvironment of Colon Adenocarcinoma

Frontiers in Genetics ◽

10.3389/fgene.2021.778153 ◽

2022 ◽

Vol 12 ◽

Author(s):

Shenghua Pan ◽

Tingting Tang ◽

Yanke Wu ◽

Liang Zhang ◽

Zekai Song ◽

...

Keyword(s):

Gene Expression ◽

Network Analysis ◽

Tumor Microenvironment ◽

Signaling Pathway ◽

Colon Adenocarcinoma ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Kaplan Meier

The tumor microenvironment (TME) has been shown to be involved in angiogenesis, tumor metastasis, and immune response, thereby affecting the treatment and prognosis of patients. This study aims to identify genes that are dysregulated in the TME of patients with colon adenocarcinoma (COAD) and to evaluate their prognostic value based on RNA omics data. We obtained 512 COAD samples from the Cancer Genome Atlas (TCGA) database and 579 COAD patients from the independent dataset (GSE39582) in the Gene Expression Omnibus (GEO) database. The immune/stromal/ESTIMATE score of each patient based on their gene expression was calculated using the ESTIMATE algorithm. Kaplan–Meier survival analysis, Cox regression analysis, gene functional enrichment analysis, and protein–protein interaction (PPI) network analysis were performed. We found that immune and stromal scores were significantly correlated with COAD patients’ overall survival (log rank p < 0.05). By comparing the high immune/stromal score group with the low score group, we identified 688 intersection differentially expressed genes (DEGs) from the TCGA dataset (663 upregulated and 25 downregulated). The functional enrichment analysis of intersection DEGs showed that they were mainly enriched in the immune process, cell migration, cell motility, Toll-like receptor signaling pathway, and PI3K-Akt signaling pathway. The hub genes were revealed by PPI network analysis. Through Kaplan–Meier and Cox analysis, four TME-related genes that were significantly related to the prognosis of COAD patients were verified in GSE39582. In addition, we uncovered the relationship between the four prognostic genes and immune cells in COAD. In conclusion, based on the RNA expression profiles of 1091 COAD patients, we screened four genes that can predict prognosis from the TME, which may serve as candidate prognostic biomarkers for COAD.

Download Full-text

Identify Key Genes by Weighted Gene Co-Expression Network Analysis for Lung Adenocarcinoma

Nano LIFE ◽

10.1142/s1793984419400026 ◽

2019 ◽

Vol 09 (01n02) ◽

pp. 1940002

Author(s):

Jichen Xu ◽

Xianchun Zong ◽

Qianshu Ren ◽

Hongyu Wang ◽

Lijuan Zhao ◽

...

Keyword(s):

Gene Expression ◽

Network Analysis ◽

Lung Adenocarcinoma ◽

Expression Profiles ◽

Gene Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Oocyte Meiosis ◽

Key Genes

The aim of this paper is to identify key genes in lung adenocarcinoma (LUAD) through weighted gene co-expression network analysis (WGCNA), and to further understand the molecular mechanism of LUAD. 107 gene expression profiles were downloaded from GSE10072 in the GEO database. We performed rigorous processing of the initial gene expression profile data. Subsequently, we used WGCNA to identify disease-driven modules and enforced functional enrichment analysis. The key genes were defined as the most connected genes in the driver module and were validated using the GSE75037 and TCGA database. GSE10072 removed 41 unpaired lung samples and 4 outliers. By analyzing the 62 samples using WGCNA, we obtained 26 modules and identified the brown and magenta modules as the driving modules for the LUAD. We found that the “Cell cycle”, “Oocyte meiosis” and “Progesterone-mediated oocyte maturation” pathways may be related to the occurrence of LUAD. GSE75037 removed 8 outlier and obtained 2909 differentially expressed genes (DEGs), 26 genes (9 genes in the brown module, 17 genes in the magenta module) overlap with key genes in the driver module. The results of the survival analysis suggest that 19 genes were significantly correlated with the patient’s survival time, including KPNA2, FEN1, RRM2, TOP2A, CENPF, MCM4, BIRC5, MELK, MAD2L1, CCNB1, CCNA2, KIF11, CDKN3, NUSAP1, CEP55, AURKA, NEK2, KIF14 and CDCA8, which may be potential biomarkers or therapeutic targets for LUAD. In this study, we provide a theoretical basis for further understanding the biological mechanism of LUAD through bioinformatics analysis of LUAD.

Download Full-text

Identification of mRNA-miRNA-lncRNA Competing Endogenous RNA Network in Adrenocortical Carcinoma Using Bioinformatics Analysis

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

2020 ◽

Author(s):

Senlin Ye ◽

Haohui Wang ◽

Wei Li ◽

Lu Yi

Keyword(s):

Adrenocortical Carcinoma ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Node Degree ◽

Pathway Enrichment Analysis ◽

Competing Endogenous Rna ◽

Hub Genes ◽

Key Genes ◽

Worse Prognosis

Abstract Background: Adrenocortical carcinoma (ACC) is a rare malignant tumor originating from the adrenal cortex. However, there are no eﬀective therapies to treat patients with ACC. LncRNA participates in a variety of biological processes of cancers. We constructed ceRNA network and identify key competing endogenous RNAs (ceRNAs) in adrenocortical carcinoma (ACC) using bioinformatic processing tools. Methods: Firstly, the differentially expressed genes (DEGs) were identified by analyzing GSE12368 and GSE19750 datasets. SangerBox was used to generate volcano maps. DAVID database was used for functional enrichment analysis. STRING database was used to conduct Protein-protein interaction (PPI) network, and hub genes were identified by Cytoscape plug-in CytoHubba. UCSC database was used to construct hierarchical clustering of hub genes. Upstream miRNAs of mRNAs were predicted by miRTarBase and upstream lncRNAs of miRNA by miRNet. Expression analysis for lncRNAs were performed via GEPIA. Prognostic analysis for genes, miRNAs and lncRNAs were performed via cBioPortal, OncomiR and GEPIA, respectively. Results: In this study, 49 and 276 upregulated and downregulated significant DEGs were identified. KEGG pathway enrichment analysis showed that they were significantly enriched in cancer-associated pathways. According to node degree, the top 10 upregulated genes and downregulated genes were classfied as hub genes. However, only 9 hub genes were defined as key genes because alteration was significantly associated with worse prognosis and all the 9 key genes were upregulated hub genes. Then, 15 miRNAs were predicted to target the 7 out of 9 key genes. But only 4 miRNAs were defined as key miRNAs because alteration significantly influenced prognosis in cancer. 185 lncRNAs were predicted to potentially interaction with the 4 miRNAs. Only 3 lncRNAs(XIST, HOXA11-AS and TMPO-AS1) were up-regulated and only 1 lncRNA (HOXA11-AS ) indicated alteration was significantly associated with worse prognosis in adrenocortical carcinoma. HOXA11-AS were finally identified as key lncRNA. Finally, RRM2-miR-24-3p/let-7a-5p-HOXA11-AS, CDK1/MCM4-miR-24-3P-HOXA11-AS competing endogenous RNA (ceRNA) sub-networks were constructed in adrenocortical carcinoma. Conclution：This study has constructed RRM2-miR-24-3p/let-7a-5p-HOXA11-AS, CDK1/MCM4-miR-24-3p-HOXA11-AS competing endogenous RNA (ceRNA) sub-networks. Our results suggested that these sub-networks might be potential therapeutic targets or prognostic biomarkers in ACC.

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

The Role of Circular RNA CDR1as/ciRS-7 in Regulating Tumor Microenvironment: A Pan-Cancer Analysis

Biomolecules ◽

10.3390/biom9090429 ◽

2019 ◽

Vol 9 (9) ◽

pp. 429 ◽

Cited By ~ 24

Author(s):

Zou ◽

Zheng ◽

Deng ◽

Yang ◽

Xie ◽

...

Keyword(s):

Tumor Microenvironment ◽

Signaling Pathway ◽

Malignant Tumors ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Circular Rna ◽

Gene Set Enrichment Analysis ◽

Functional Enrichment ◽

Receptor Interaction

Circular RNA CDR1as/ciRS-7 functions as an oncogenic regulator in various cancers. However, there has been a lack of systematic and comprehensive analysis to further elucidate its underlying role in cancer. In the current study, we firstly performed a bioinformatics analysis of CDR1as among 868 cancer samples by using RNA-seq datasets of the MiOncoCirc database. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment analysis (GSEA), CIBERSORT, Estimating the Proportion of Immune and Cancer cells (EPIC), and the MAlignant Tumors using Expression data (ESTIMATE) algorithm were applied to investigate the underlying functions and pathways. Functional enrichment analysis suggested that CDR1as has roles associated with angiogenesis, extracellular matrix (ECM) organization, integrin binding, and collagen binding. Moreover, pathway analysis indicated that it may regulate the TGF-β signaling pathway and ECM-receptor interaction. Therefore, we used CIBERSORT, EPIC, and the ESTIMATE algorithm to investigate the association between CDR1as expression and the tumor microenvironment. Our data strongly suggest that CDR1as may play a specific role in immune and stromal cell infiltration in tumor tissue, especially those of CD8+ T cells, activated NK cells, M2 macrophages, cancer-associated fibroblasts (CAFs) and endothelial cells. Generally, systematic and comprehensive analyses of CDR1as were conducted to shed light on its underlying pro-cancerous mechanism. CDR1as regulates the TGF-β signaling pathway and ECM-receptor interaction to serve as a mediator in alteration of the tumor microenvironment.

Download Full-text

Identification of PKP 2/3 as potential biomarkers of ovarian cancer based on bioinformatics and experiments

Cancer Cell International ◽

10.1186/s12935-020-01602-3 ◽

2020 ◽

Vol 20 (1) ◽

Author(s):

Lingling Gao ◽

Xiao Li ◽

Qian Guo ◽

Xin Nie ◽

Yingying Hao ◽

...

Keyword(s):

Ovarian Cancer ◽

Signaling Pathway ◽

Poor Prognosis ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Biological Behavior ◽

Cell Junction ◽

Immune Infiltration ◽

Potential Biomarker

Abstract Background Plakophilins (PKPs) are widely involved in gene transcription, translation, and signal transduction, playing a crucial role in tumorigenesis and progression. However, the function and potential mechanism of PKP1/2/3 in ovarian cancer (OC) remains unclear. It’s of great value to explore the expression and prognostic values of PKP1/2/3 and their potential mechanisms, immune infiltration in OC. Methods The expression levels, prognostic values and genetic variations of PKP1/2/3 in OC were explored by various bioinformatics tools and databases, and PKP2/3 were selected for further analyzing their regulation network and immune infiltration. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) enrichment were also conducted. Finally, the expression and prognosis of PKP2 were validated by immunohistochemistry. Results The expression level and prognosis of PKP1 showed little significance in ovarian cancer, and the expression of PKP2/3 mRNA and protein were upregulated in OC, showing significant correlations with poor prognosis of OC. Functional enrichment analysis showed that PKP2/3 and their correlated genes were significantly enriched in adaptive immune response, cytokine receptor activity, organization of cell–cell junction and extracellular matrix; KEGG analysis showed that PKP2/3 and their significantly correlated genes were involved in signaling pathways including cytokine-mediated signaling pathway, receptor signaling pathway and pathways in cancer. Moreover, PKP2/3 were correlated with lymphocytes and immunomodulators. We confirmed that high expression of PKP2 was significantly associated with advanced stage, poor differentiation and poor prognosis of OC patients. Conclusion Members of plakophilins family showed various degrees of abnormal expressions and prognostic values in ovarian cancer. PKP2/3 played crucial roles in tumorigenesis, aggressiveness, malignant biological behavior and immune infiltration of OC, and can be regarded as potential biomarker for early diagnosis and prognosis evaluation in OC.

Download Full-text

Identification of nine key genes by bioinformatics analysis for predicting poor prognosis in smoking-induced lung adenocarcinoma

Lung Cancer Management ◽

10.2217/lmt-2020-0009 ◽

2020 ◽

Vol 9 (2) ◽

pp. LMT30

Author(s):

Chuanli Ren ◽

Weixiu Sun ◽

Xu Lian ◽

Chongxu Han

Keyword(s):

Lung Adenocarcinoma ◽

Differentially Expressed Genes ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Kaplan Meier ◽

Key Genes ◽

Core Genes

Aim: To screen and identify key genes related to the development of smoking-induced lung adenocarcinoma (LUAD). Materials & methods: We obtained data from the GEO chip dataset GSE31210. The differentially expressed genes were screened by GEO2R. The protein interaction network of differentially expressed genes was constructed by STRING and Cytoscape. Finally, core genes were screened. The overall survival time of patients with the core genes was analyzed by Kaplan–Meier method. Gene ontology and Kyoto encyclopedia of genes and genomes bioaccumulation was calculated by DAVID. Results: Functional enrichment analysis indicated that nine key genes were actively involved in the biological process of smoking-related LUAD. Conclusion: 23 core genes and nine key genes among them were correlated with adverse prognosis of LUAD induced by smoking.

Download Full-text

Differential co-expression network analysis for identification of functional modules in HIV infection and tubercolosis-HIV co-infection

Journal of Biological Research - Bollettino della Società Italiana di Biologia Sperimentale ◽

10.4081/jbr.2016.5950 ◽

2016 ◽

Vol 89 (2) ◽

Cited By ~ 1

Author(s):

Mohit Jha ◽

Anvita Gupta ◽

Sudha Singh ◽

Khushhali Menaria Pandey

Keyword(s):

Network Analysis ◽

Hiv Infection ◽

Enrichment Analysis ◽

Gene Interaction ◽

Gene Clusters ◽

Functional Enrichment Analysis ◽

Gene Expression Omnibus ◽

Functional Enrichment ◽

Specific Gene ◽

Diagnostic Strategies

Co-infection with tuberculosis (TB) is the preeminent cause of demise in human immunodeficiency virus (HIV) infected individuals. However, diagnosis of TB, particularly in the presence of an HIV co-infection, can be limiting owing to the high inaccuracy associated with conventional diagnostic strategies. Here we determine dysregulated pathways in TB-HIV co-infection and HIV infection utilizing coexpression networks. Primarily, we utilized preservation statistics to identify gene modules that exhibit a weak conservation of network topology within HIV infected and TB-HIV co-infected networks. Raw data was downloaded from Gene Expression Omnibus (GSE50834) and duly pre-processed. Co-expression networks for each condition (HIV infected and TB-HIV co-infected) were constructed independently. Preservation of HIV infected network edges was evaluated with respect to TB-HIV co-infected and vice versa using weighted correlation network analysis. Two out of the 22 modules were identified as exhibiting weak preservation in both conditions. Functional enrichment analysis identified that weakly preserved modules were pertinent to the condition under study. For instance, weakly preserved TBHIV co-infected module T1 enriched for genes associated with mitochondrion exhibited the highest fraction of gene interaction pairs exclusive to TB-HIV co-infection. Concisely, we illustrated the application of using preservation statistics to detect modules functionally linked with dysregulated pathways in disease, as exemplified by the mitochondrion module T1. Our analyses discovered gene clusters that are non-randomly linked with the disease. Highly specific gene pairs pointed to interactions between known markers of disease and favoured identification of possible markers that are likely to be associated with the disease.

Download Full-text

Identification of novel genes and pathways in carotid atheroma using integrated bioinformatic methods

Scientific Reports ◽

10.1038/srep18764 ◽

2016 ◽

Vol 6 (1) ◽

Cited By ~ 10

Author(s):

Wenqing Nai ◽

Diane Threapleton ◽

Jingbo Lu ◽

Kewei Zhang ◽

Hongyuan Wu ◽

...

Keyword(s):

Signaling Pathway ◽

Transforming Growth Factor ◽

Interaction Network ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Transforming Growth Factor Β ◽

Functional Enrichment ◽

Differentially Expressed ◽

Pathway Enrichment Analysis ◽

Atheromatous Plaques

Abstract Atherosclerosis is the primary cause of cardiovascular events and its molecular mechanism urgently needs to be clarified. In our study, atheromatous plaques (ATH) and macroscopically intact tissue (MIT) sampled from 32 patients were compared and an integrated series of bioinformatic microarray analyses were used to identify altered genes and pathways. Our work showed 816 genes were differentially expressed between ATH and MIT, including 443 that were up-regulated and 373 that were down-regulated in ATH tissues. GO functional-enrichment analysis for differentially expressed genes (DEGs) indicated that genes related to the “immune response” and “muscle contraction” were altered in ATHs. KEGG pathway-enrichment analysis showed that up-regulated DEGs were significantly enriched in the “FcεRI-mediated signaling pathway”, while down-regulated genes were significantly enriched in the “transforming growth factor-β signaling pathway”. Protein-protein interaction network and module analysis demonstrated that VAV1, SYK, LYN and PTPN6 may play critical roles in the network. Additionally, similar observations were seen in a validation study where SYK, LYN and PTPN6 were markedly elevated in ATH. All in all, identification of these genes and pathways not only provides new insights into the pathogenesis of atherosclerosis, but may also aid in the development of prognostic and therapeutic biomarkers for advanced atheroma.

Download Full-text