scholarly journals Identification of hub genes and construction of an mRNA-miRNA-lncRNA network of gastric carcinoma using integrated bioinformatics analysis

PLoS ONE ◽  
2021 ◽  
Vol 16 (12) ◽  
pp. e0261728
Author(s):  
Gang Wei ◽  
Youhong Dong ◽  
Zhongshi He ◽  
Hu Qiu ◽  
Yong Wu ◽  
...  

Background Gastric carcinoma (GC) is one of the most common cancer globally. Despite its worldwide decline in incidence and mortality over the past decades, gastric cancer still has a poor prognosis. However, the key regulators driving this process and their exact mechanisms have not been thoroughly studied. This study aimed to identify hub genes to improve the prognostic prediction of GC and construct a messenger RNA-microRNA-long non-coding RNA(mRNA-miRNA-lncRNA) regulatory network. Methods The GSE66229 dataset, from the Gene Expression Omnibus (GEO) database, and The Cancer Genome Atlas (TCGA) database were used for the bioinformatic analysis. Differential gene expression analysis methods and Weighted Gene Co-expression Network Analysis (WGCNA) were used to identify a common set of differentially co-expressed genes in GC. The genes were validated using samples from TCGA database and further validation using the online tools GEPIA database and Kaplan-Meier(KM) plotter database. Gene set enrichment analysis(GSEA) was used to identify hub genes related to signaling pathways in GC. The RNAInter database and Cytoscape software were used to construct an mRNA-miRNA-lncRNA network. Results A total of 12 genes were identified as the common set of differentially co-expressed genes in GC. After verification of these genes, 3 hub genes, namely CTHRC1, FNDC1, and INHBA, were found to be upregulated in tumor and associated with poor GC patient survival. In addition, an mRNA-miRNA-lncRNA regulatory network was established, which included 12 lncRNAs, 5 miRNAs, and the 3 hub genes. Conclusions In summary, the identification of these hub genes and the establishment of the mRNA-miRNA-lncRNA regulatory network provide new insights into the underlying mechanisms of gastric carcinogenesis. In addition, the identified hub genes, CTHRC1, FNDC1, and INHBA, may serve as novel prognostic biomarkers and therapeutic targets.

2020 ◽  
Vol 2020 ◽  
pp. 1-12 ◽  
Author(s):  
Kaisheng Liu ◽  
Minshan Lai ◽  
Shaoxiang Wang ◽  
Kai Zheng ◽  
Shouxia Xie ◽  
...  

Colon cancer is the third most common cancer, with a high incidence and mortality. Construction of a specific and sensitive prediction model for prognosis is urgently needed. In this study, profiles of patients with colon cancer with clinical and gene expression data were downloaded from Gene Expression Omnibus and The Cancer Genome Atlas (TCGA). CXC chemokines in patients with colon cancer were investigated by differential expression gene analysis, overall survival analysis, receiver operating characteristic analysis, gene set enrichment analysis (GSEA), and weighted gene coexpression network analysis. CXCL1, CXCL2, CXCL3, and CXCL11 were upregulated in patients with colon cancer and significantly correlated with prognosis. The area under curve (AUC) of the multigene forecast model of CXCL1, CXCL11, CXCL2, and CXCL3 was 0.705 in the GSE41258 dataset and 0.624 in TCGA. The prediction model was constructed using the risk score of the multigene model and three clinicopathological risk factors and exhibited 92.6% and 91.8% accuracy in predicting 3-year and 5-year overall survival of patients with colon cancer, respectively. In addition, by GSEA, expression of CXCL1, CXCL11, CXCL2, and CXCL3 was correlated with several signaling pathways, including NOD-like receptor, oxidative phosphorylation, mTORC1, interferon-gamma response, and IL6/JAK/STAT3 pathways. Patients with colon cancer will benefit from this prediction model for prognosis, and this will pave the way to improve the survival rate and optimize treatment for colon cancer.


2021 ◽  
Vol 10 ◽  
Author(s):  
Ji’an Yang ◽  
Qian Yang

Glioblastoma multiforme is the most common primary intracranial malignancy, but its etiology and pathogenesis are still unclear. With the deepening of human genome research, the research of glioma subtype screening based on core molecules has become more in-depth. In the present study, we screened out differentially expressed genes (DEGs) through reanalyzing the glioblastoma multiforme (GBM) datasets GSE90598 from the Gene Expression Omnibus (GEO), the GBM dataset TCGA-GBM and the low-grade glioma (LGG) dataset TCGA-LGG from the Cancer Genome Atlas (TCGA). A total of 150 intersecting DEGs were found, of which 48 were upregulated and 102 were downregulated. These DEGs from GSE90598 dataset were enriched using the overrepresentation method, and multiple enriched gene ontology (GO) function terms were significantly correlated with neural cell signal transduction. DEGs between GBM and LGG were analyzed by gene set enrichment analysis (GSEA), and the significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways involved in synapse signaling and oxytocin signaling pathways. Then, a protein-protein interaction (PPI) network was constructed to assess the interaction of proteins encoded by the DEGs. MCODE identified 2 modules from the PPI network. The 11 genes with the highest degrees in module 1 were designated as core molecules, namely, GABRD, KCNC1, KCNA1, SYT1, CACNG3, OPALIN, CD163, HPCAL4, ANK3, KIF5A, and MS4A6A, which were mainly enriched in ionic signaling-related pathways. Survival analysis of the GSE83300 dataset verified the significant relationship between expression levels of the 11 core genes and survival. Finally, the core molecules of GBM and the DrugBank database were assessed by a hypergeometric test to identify 10 drugs included tetrachlorodecaoxide related to cancer and neuropsychiatric diseases. Further studies are required to explore these core genes for their potentiality in diagnosis, prognosis, and targeted therapy and explain the relationship among ionic signaling-related pathways, neuropsychiatric diseases and neurological tumors.


2021 ◽  
Author(s):  
Sheng Fang ◽  
Xiao Fang ◽  
Xin Xu ◽  
Lin Zhong ◽  
An-quan Wang ◽  
...  

Abstract Relevance Rheumatoid arthritis (RA) is a systemic autoimmune disease with an aggressive, chronic synovial inflammation as the main pathological change. However, the specific etiology, pathogenesis, and related biomarkers in diagnosis and treatment are still not fully elucidated. This study attempts to provide new perspectives and insights into RA at the genetic, molecular, and cellular levels through the tenet of personalized medicine. Methods Gene expression profiles of four individual knee synovial tissues were downloaded from a comprehensive gene expression database, R language was used to screen for significantly differentially expressed genes (DEGs), Gene Ontology Enrichment Analysis, Kyoto Gene Encyclopedia, and Gene Set Enrichment Analysis were performed to analyze the biological functions and signaling pathways of these DEGs, STRING online database was used to establish protein-protein interaction networks, Cytoscape software to obtain ten hub genes, Goplot to get six inflammatory immune-related hub genes, and CIBERSORT algorithm to impute immune infiltration. Results Molecular pathways that play important roles in RA were obtained: Toll-like receptors, AMPK, MAPK, TNF, FoxO, TGF-beta, PI3K and NF-κB pathways, Ten hub genes: Ccr1, Ccr2, Ccr5, Ccr7, Cxcl5, Cxcl6, Cxcl13, Ccl13, Adcy2, and Pnoc. among which Adcy2 and Pnoc have not been reported in RA studies, suggesting that they may be worthy targets for further study. It was also found that among the synoviocytes in RA, the proportions of plasma cells, CD8 T cells, follicular helper T cells, monocytes, γ delta T cells, and M0 macrophages were higher, while the proportions of CD4 memory resting T cells, regulatory T cells (Tregs), activated NK cells, resting dendritic cells, M1 macrophages, eosinophils, activated mast cells, resting mast cells were lower in proportion, and each cell played an important role in RA. Conclusions This study may help understand the key genes, molecular pathways, the role of inflammatory immune infiltrating cells in RA’s pathogenesis and provide new targets and ideas for the diagnosis and personalized treatment of RA.


2021 ◽  
Author(s):  
Chao Zhang ◽  
Feng Xu ◽  
Fang Fang

Abstract Background: Sepsis-associated acute lung injury (ALI) is a potentially lethal complication associated with a poor prognosis and high mortality worldwide, especially in the outbreak of COVID-19. However, the fundamental mechanisms of this complication were still not fully elucidated. Thus, we conducted this study to identify hub genes and biological pathways of sepsis-associated ALI, mainly focus on two pathways of LPS and HMGB1. Methods: Gene expression profile GSE3037 were downloaded from Gene Expression Omnibus (GEO) database, including 8 patients with sepsis-induced acute lung injury, with 8 unstimulated blood neutrophils, 8 LPS- induced neutrophils and 8 HMGB1-induced neutrophils. Differentially expressed genes (DEGs) identifications, Gene Ontology (GO) function analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, Gene Set Enrichment Analysis (GSEA) and protein-protein interaction (PPI) network constructions were performed to obtain hub genes and relevant biological pathways.Results: We identified 534 and 317 DEGs for LPS- and HMGB1-induced ALI, respectively. The biological pathways involved in LPS- and HMGB1-induced ALI were also identified accordingly. By PPI network analysis, we found that ten hub genes for LPS-induced ALI (CXCL8, TNF, IL6, IL1B, ICAM1, CXCL1, CXCL2, IL1A, IL1RN and CXCL3) and another ten hub genes for HMGB1-induced ALI (CCL20, CXCL2, CXCL1, CCL4, CXCL3, CXCL9, CCL21, CXCR6, KNG1 and SST). Furthermore, by combining analysis, the results revealed that genes of TNF, CCL20, IL1B, NFKBIA, CCL4, PTGS2, TNFAIP3, CXCL2, CXCL1 and CXCL3 were potential biomarkers for sepsis-associated ALI. Conclusions: Our study revealed that ten hub genes associated with sepsis-induced ALI were TNF, CCL20, IL1B, NFKBIA, CCL4, PTGS2, TNFAIP3, CXCL2, CXCL1 and CXCL3, which may serve as genetic biomarkers and be further verified in prospective experimental trials.


2020 ◽  
Author(s):  
Zhenhua Yin ◽  
Dejun Wu ◽  
Jianping Shi ◽  
Xiyi Wei ◽  
Nuyun Jin ◽  
...  

Abstract Background: Extensive research has revealed that genes play a pivotal role in tumor development and growth. However, the underlying involvement of gene expression in gastric carcinoma (GC) remains to be investigated further.Methods: In this study, we identified overlapping differentially expressed genes (DEGs) by comparing tumor tissue with adjacent normal tissue using the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) database.Results: Our analysis identified 79 up-regulated and ten down-regulated genes. Functional enrichment analysis and prognosis analysis were conducted on the identified genes, and the fatty aldehyde dehydrogenase (FALDH) gene, ALDH3A2, was chosen for more detailed analysis. We performed Gene Set Enrichment Analysis (GSEA) and immunocorrelation analysis (infiltration, copy number alterations, and checkpoints) to elucidate the mechanisms of action of ALDH3A2 in depth. The immunohistochemical (IHC) result based on 140 paraffin-embedded human GC samples indicated that ALDH3A2 was over-expressed in low-grade GC cases and the OS of patients with low expression of ALDH3A2 was significantly shorter than those with high ALDH3A2 expression. In vitro results indicated that the expression of ALDH3A2 was negatively correlated with PDCD1, PDCD1LG2, and CTLA-4.Conclusion: We conclude that ALDH3A2 might be useful as a potential reference value for the relief and immunotherapy of GC, and also as an independent predictive marker for the prognosis of GC.


2021 ◽  
Author(s):  
Zimeng Wei ◽  
Min Zhao ◽  
Linnan Zang

Abstract Background Lung adenocarcinoma (LUAD) is the main histological subtype of lung cancer. However, the molecular mechanism underlying LUAD is not yet clearly defined, but elucidating this process in detail would be of great significance for clinical diagnosis and treatment. Methods Gene expression profiles were retrieved from Gene Expression Omnibus database (GEO), and the common differentially expressed genes (DEGs) were identified by online GEO2R analysis tool. Subsequently, the enrichment analysis of function and signaling pathways of DEGs in LUAD were performed by gene ontology (GO) and The Kyoto Encyclopedia of Genes and Genomics (KEGG) analysis. The protein-protein interaction (PPI) networks of the DEGs were established through the Search Tool for the Retrieval of Interacting Genes (STRING) database and hub genes were screened by plug-in CytoHubba in Cytoscape. Afterwards, we detected the expression of hub genes in LUAD and other cancers via GEPIA, Oncomine and HPA databases. Finally, Kaplan-Meier plotter were performed to analyze the prognosis efficacy of hub genes. Results 74 up-regulated and 238 down-regulated DEGs were identified. As for the up-regulated DEGs, KEGG analysis results revealed they were mainly enrolled in protein digestion and absorption. However, the down-regulated DEGs were primarily enriched in cell adhesion molecules. Subsequently, 9 hub genes: KIAA0101, CDCA7, TOP2A, CDC20, ASPM, TPX2, CENPF, UBE2T and ECT2, were identified and showed higher expression in both LUAD and other cancers. Finally, all these hub genes were found significantly related to the prognosis of LUAD (p < 0.05). Conclusions Our results screened out the hub genes and pathways that were related to the development and prognosis of LUAD, which could provide new insight for the future molecularly targeted therapy and prognosis evaluation of LUAD.


2021 ◽  
Author(s):  
Mi Jiang ◽  
Jia Li ◽  
Zhi Song

Abstract Background: Epilepsy is a complicated neurological disorder with almost 30% refractory. Recent years, several studies showed that epilepsy is associated with its comorbidities by shared molecular mechanisms. However, the association of epilepsy and digestive comorbidities are still unclear. In this study, we aim to explore the association between inflammatory bowel disease (IBD) and epilepsy, and to find promising therapeutic targets for refractory epilepsy. Methods: Two gene expression profiles (GSE134697 and GSE59071) were selected from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified by GEO2R and the DESeq2 package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of intersection DEGs and Gene Set Enrichment Analysis (GSEA) were conducted by clusterProfiler package. The protein-protein interactions (PPI) network was established by using STRING and visualized by Cytoscape. Genes in the most significant module identified by MCODE plug-in were considered as candidate hub genes. Validation of hub genes were performed by using the GSE143272 dataset. Results: Cytokine-cytokine receptor interaction pathway is identified as a communal pathway between IBD and epilepsy. CXCL8, CXCR4 and ITGAX were identified as the hub genes. Conclusions: The identification of the communal pathway and hub genes in this study contributes to a potential novel therapeutic target in refractory epilepsy.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12605
Author(s):  
Tongtong Zhang ◽  
Suyang Yu ◽  
Shipeng Zhao

Background Gastric cancer (GC) is the most prevalent malignancy among the digestive system tumors. Increasing evidence has revealed that lower mRNA expression of ANXA9 is associated with a poor prognosis in colorectal cancer. However, the role of ANXA9 in GC remains largely unknown. Material and Methods The Gene Expression Profiling Interactive Analysis (GEPIA) and Human Protein Atlas databases were used to investigate the expression of ANXA9 in GC, which was then validated in the four Gene Expression Omnibus (GEO) datasets. The diagnostic value of ANXA9 for GC patients was demonstrated using a receiver operating characteristic (ROC) curve. The correlation between ANXA9 expression and clinicopathological parameters was analyzed in The Cancer Genome Atlas (TCGA) and UALCAN databases. The Kaplan-Meier (K-M) survival curve was used to elucidate the relationship between ANXA9 expression and the survival time of GC patients. We then performed a gene set enrichment analysis (GSEA) to explore the biological functions of ANXA9. The relationship of ANXA9 expression and cancer immune infiltrates was analyzed using the Tumor Immune Estimation Resource (TIMER). In addition, the potential mechanism of ANXA9 in GC was investigated by analyzing its related genes. Results ANXA9 was significantly up-regulated in GC tissues and showed obvious diagnostic value. The expression of ANXA9 was related to the age, gender, grade, TP53 mutation, and histological subtype of GC patients. We also found that ANXA9 expression was associated with immune-related biological function. ANXA9 expression was also correlated with the infiltration level of CD8+ T cells, neutrophils, and dendritic cells in GC. Additionally, copy number variation (VNV) of ANXA9 occurred in GC patients. Function enrichment analyses revealed that ANXA9 plays a role in the GC progression by interacting with its related genes. Conclusions Our results provide strong evidence of ANXA9 expression as a prognostic indicator related to immune responses in GC.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Jingyi Chen ◽  
Yuxuan Song ◽  
Mei Li ◽  
Yu Zhang ◽  
Tingru Lin ◽  
...  

Abstract Background Competing endogenous RNA (ceRNA) represents a class of RNAs (e.g., long noncoding RNAs [lncRNAs]) with microRNA (miRNA) binding sites, which can competitively bind miRNA and inhibit its regulation of target genes. Increasing evidence has underscored the involvement of dysregulated ceRNA networks in the occurrence and progression of colorectal cancer (CRC). The purpose of this study was to construct a ceRNA network related to the prognosis of CRC and further explore the potential mechanisms that affect this prognosis. Methods RNA-Seq and miRNA-Seq data from The Cancer Genome Atlas (TCGA) were used to identify differentially expressed lncRNAs (DElncRNAs), microRNAs (DEmiRNAs), and mRNAs (DEmRNAs), and a prognosis-related ceRNA network was constructed based on DElncRNA survival analysis. Subsequently, pathway enrichment, Pearson correlation, and Gene Set Enrichment Analysis (GSEA) were performed to determine the function of the genes in the ceRNA network. Gene Expression Profiling Interactive Analysis (GEPIA) and immunohistochemistry (IHC) were also used to validate differential gene expression. Finally, the correlation between lncRNA and immune cell infiltration in the tumor microenvironment was evaluated based on the CIBERSORT algorithm. Results A prognostic ceRNA network was constructed with eleven key survival-related DElncRNAs (MIR4435-2HG, NKILA, AFAP1-AS1, ELFN1-AS1, AC005520.2, AC245884.8, AL354836.1, AL355987.4, AL591845.1, LINC02038, and AC104823.1), 54 DEmiRNAs, and 308 DEmRNAs. The MIR4435-2HG- and ELFN1-AS1-associated ceRNA subnetworks affected and regulated the expression of the COL5A2, LOX, OSBPL3, PLAU, VCAN, SRM, and E2F1 target genes and were found to be related to prognosis and tumor-infiltrating immune cell types. Conclusions MIR4435-2HG and ELFN1-AS1 are associated with prognosis and tumor-infiltrating immune cell types and could represent potential prognostic biomarkers or therapeutic targets in colorectal carcinoma.


2020 ◽  
Author(s):  
Zhenhua Yin ◽  
Dejun Wu ◽  
Jianping Shi ◽  
Xiyi Wei ◽  
Nuyun Jin ◽  
...  

Abstract Background: Extensive research has revealed that genes play a pivotal role in tumor development and growth. However, the underlying involvement of gene expression in gastric carcinoma (GC) remains to be investigated further. Methods: In this study, we identified overlapping differentially expressed genes (DEGs) by comparing tumor tissue with adjacent normal tissue using the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) database. Results: Our analysis identified 79 up-regulated and ten down-regulated genes. Functional enrichment analysis and prognosis analysis were conducted on the identified genes, and the fatty aldehyde dehydrogenase (FALDH) gene, ALDH3A2, was chosen for more detailed analysis. We performed Gene Set Enrichment Analysis (GSEA) and immunocorrelation analysis (infiltration, copy number alterations, and checkpoints) to elucidate the mechanisms of action of ALDH3A2 in depth. The immunohistochemical (IHC) result based on 140 paraffin-embedded human GC samples indicated that ALDH3A2 was over-expressed in low-grade GC cases and the OS of patients with low expression of ALDH3A2 was significantly shorter than those with high ALDH3A2 expression. In vitro results indicated that the expression of ALDH3A2 was negatively correlated with PDCD1, PDCD1LG2, and CTLA-4. Conclusion: We conclude that ALDH3A2 might be useful as a potential reference value for the relief and immunotherapy of GC, and also as an independent predictive marker for the prognosis of GC.


Sign in / Sign up

Export Citation Format

Share Document