Dataset Information

A Weighted Gene Co-Expression Network Analysis-Derived Prognostic Model for Predicting Prognosis and Immune Infiltration in Gastric Cancer.

ABSTRACT:

Background

Gastric cancer (GC) is a major public health problem worldwide. In recent decades, the treatment of gastric cancer has improved greatly, but basic research and clinical application of gastric cancer remain challenges due to the high heterogeneity. Here, we provide new insights for identifying prognostic models of GC.

Methods

We obtained the gene expression profiles of GSE62254 containing 300 samples for training. GSE15459 and TCGA-STAD for validation, which contain 200 and 375 samples, respectively. Weighted gene co-expression network analysis (WGCNA) was used to identify gene modules. We performed Lasso regression and Cox regression analyses to identify the most significant five genes to develop a novel prognostic model. And we selected two representative genes within the model for immunohistochemistry staining with 105 GC specimens from our hospital to verify the prediction efficiency. Moreover, we estimated the correlation coefficient between our model and immune infiltration using the CIBERSORT algorithm. The data from GSE15459 and TCGA cohort validated the robustness and predictive accuracy of this prognostic model.

Results

Of the 12 gene modules identified, 1,198 green-yellow module genes were selected for further analysis. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analysis using the Cox proportional hazards regression model. Finally, we constructed a five gene prognostic model: Risk Score = [(-0.7547) * Expression (ARHGAP32)] + [(-0.8272) * Expression (KLF5)] + [1.09 * Expression (MAMLD1)] + [0.5174 * Expression (MATN3)] + [1.66 * Expression (NES)]. The prognosis of samples in the high-risk group was significantly poorer than that of samples in the low-risk group (p = 6.503e-11). The risk model was also regarded as an independent predictor of prognosis (HR, 1.678, p < 0.001). The observed correlation with immune cells suggested that this risk model could potentially predict immune infiltration.

Conclusion

This study identified a potential risk model for prognosis and immune infiltration prediction in GC using WGCNA and Cox regression analysis.

SUBMITTER: Chen Q

PROVIDER: S-EPMC7947930 | biostudies-literature | 2021

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A Weighted Gene Co-Expression Network Analysis-Derived Prognostic Model for Predicting Prognosis and Immune Infiltration in Gastric Cancer.

Chen Qingchuan Q Tan Yuen Y Zhang Chao C Zhang Zhe Z Pan Siwei S An Wen W Xu Huimian H

Frontiers in oncology 20210225

<h4>Background</h4>Gastric cancer (GC) is a major public health problem worldwide. In recent decades, the treatment of gastric cancer has improved greatly, but basic research and clinical application of gastric cancer remain challenges due to the high heterogeneity. Here, we provide new insights for identifying prognostic models of GC.<h4>Methods</h4>We obtained the gene expression profiles of GSE62254 containing 300 samples for training. GSE15459 and TCGA-STAD for validation, which contain 200 ...[more]

PMID: 33718128

Similar Datasets

Project description:The dreadful prognosis of hepatocellular carcinoma (HCC) is primarily due to the low early diagnosis rate, rapid progression, and high recurrence rate. Valuable prognostic biomarkers are urgently needed for HCC. In this study, microarray data were downloaded from GSE14520, GSE22058, International Cancer Genome Consortium (ICGC), and The Cancer Genome Atlas (TCGA). Differentially expressed genes (DEGs) were identified among GSE14520, GSE22058, and ICGC databases. Weighted gene co-expression network analysis (WGCNA) was used to establish gene co-expression modules of DEGs, and genes of key modules were examined to identify hub genes using univariate Cox regression in the ICGC cohort. Expression levels and time-dependent receiver operating characteristic (ROC) and area under the curve (AUC) were determined to estimate the prognostic competence of the hub genes. These hub genes were also validated in the Gene Expression Profiling Interactive Analysis (GEPIA) and TCGA databases. TIMER algorithm and GSCALite database were applied to analyze the association of the hub genes with immunocytotic infiltration and their pathway enrichment. Altogether, 276 DEGs were identified and WGCNA described a unique and significantly DEGs-associated co-expression module containing 148 genes, with 10 hub genes selected by univariate Cox regression in the ICGC cohort (BIRC5, FOXM1, CENPA, KIF4A, DTYMK, PRC1, IGF2BP3, KIF2C, TRIP13, and TPX2). Most of the genes were validated in the GEPIA databases, except IGF2BP3. The results of multivariate Cox regression analysis indicated that the abovementioned hub genes are all independent predictors of HCC. The 10 genes were also confirmed to be associated with immune cell infiltration using the TIMER algorithm. Moreover, four-gene signature was developed, including BIRC5, CENPA, FOXM1, DTYMK. These hub genes and the model demonstrated a strong prognostic capability and are likely to be a therapeutic target for HCC. Moreover, the association of these genes with immune cell infiltration improves our understanding of the occurrence and development of HCC.

Project description:BackgroundWhile large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed to screen the hub genes and the immune cell infiltration pattern in ULMS by bioinformatics methods.MethodsFirstly, the GSE67463 dataset, including 25 ULMS tissues and 29 normal myometrium (NL) tissues, was downloaded from the public database. The differentially expressed genes (DEGs) were screened by the 'limma' package and hub modules were identified by weighted gene co-expression network analysis (WGCNA). Subsequently, gene function annotations were performed to investigate the biological role of the genes from the intersection of two groups (hub module and DEGs). The above genes were calculated in the protein-protein interaction (PPI) network to select the hub genes further. The hub genes were validated using external data (GSE764 and GSE68295). In addition, the differential immune cell infiltration between UL and ULMS tissues was investigated using the CIBERSORT algorithm. Finally, we used western blot to preliminarily detect the hub genes in cell lines.ResultsWGCNA analysis revealed a green-yellow module possessed the highest correlation with ULMS, including 1063 genes. A total of 172 DEGs were selected by thresholds set in the 'limma' package. The above two groups of genes were intersected to obtain 72 genes for functional annotation analysis. Interestingly, it indicated that 72 genes were mainly involved in immune processes and the Neddylation pathway. We found a higher infiltration of five types of cells (memory B cells, M0-type macrophages, mast cells activated, M1-type macrophages, and T cells follicular helper) in ULMS tissues than NL tissues, while the infiltration of two types of cells (NK cells activated and mast cells resting) was lower than in NL tissues. In addition, a total of five genes (KDR, CCL21, SELP, DPT, and DCN) were identified as the hub genes. Internal and external validation demonstrated that the five genes were over-expressed in NL tissues compared with USML tissues. Finally, the correlation analysis results indicate that NK cells activated and mast cells activated positively correlated with the hub genes. However, M1-type macrophages had a negative correlation with the hub genes. Moreover, only the DCN may be associated with the Neddylation pathway.ConclusionA series of evidence confirm that the five hub genes and the infiltration of seven types of immune cells are related to USML occurrence. These hub genes may affect the occurrence of USML through immune-related and Neddylation pathways, providing molecular evidence for the treatment of USML in the future.

Project description:Background: Insulin-like growth factor binding protein-7 (IGFBP7) contributes to multiple biological processes in various tumors. However, the role of IGFBP7 in gastric cancer (GC) is still undetermined. The study aims to explore the role of IGFBP7 in GC via an integrated bioinformatics analysis. Methods: IGFBP7 expression levels in GC and its normal gastric tissues were analyzed using multiple databases, including the Tumor Immune Estimation Resource (TIMER), Oncomine, The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, as well as by our clinical gastric specimens. The methylation analysis was conducted with MEXPRESS, UALCAN and Xena online tools. The survival analysis was conducted using the Kaplan-Meier Plotter and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Coexpressed genes of IGFBP7 were selected with the cBioPortal tool and enrichment analysis was conducted with the clusterProfiler package in R software. Gene set enrichment analysis (GSEA) was performed to explore the IGFBP7-related biological processes involved in GC. Correlations between IGFBP7 and immune cell infiltrates were analyzed using the TIMER database. Results: IGFBP7 expression was significantly upregulated in GC and correlated with stage, grade, tumor status and Helicobacter pylori infection. High IGFBP7 expression and low IGFBP7 methylation levels were significantly associated with short survival of patients with GC. Univariate and multivariate analyses revealed that IGFBP7 was an independent risk factor for GC. The coexpressed genes LHFPL6, SEPTIN4, HSPB2, LAYN and GGT5 predicted unfavorable outcomes of GC. Enrichment analysis showed that the coexpressed genes were involved in extracellular matrix (ECM)-related processes. GSEA indicated that IGFBP7 was positively related to ECM and inflammation-related pathways. TIMER analysis indicated that the mRNA level of IGFBP7 was strongly correlated with genes related to various infiltrating immune cells in GC, especially with gene markers of tumor associated macrophages (TAMs). Conclusions: Increased IGFBP7 expression correlates with poor prognosis and immune cell infiltration in GC, which might be a potential biomarker for the diagnosis of GC.

Project description:Background: Increasing evidence supports that competing endogenous RNAs (ceRNAs) and tumor immune infiltration act as pivotal players in tumor progression of hepatocellular carcinoma (HCC). Nonetheless, comprehensive analysis focusing on ceRNAs and immune infiltration in HCC is lacking. Methods: RNA and miRNA sequencing information, corresponding clinical annotation, and mutation data of HCC downloaded from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) project were employed to identify significant differentially expressed mRNAs (DEMs), miRNAs (DEMis), and lncRNAs (DELs) to establish a ceRNA regulatory network. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene ontology (GO) enrichment pathways were analyzed to functionally annotate these DEMs. A multigene-based risk signature was developed utilizing least absolute shrinkage and selection operator method (LASSO) algorithm. Moreover, survival analysis and receiver operating characteristic (ROC) analysis were applied for prognostic value validation. Seven algorithms (TIMER, XCELL, MCPcounter, QUANTISEQ, CIBERSORT, EPIC, and CIBERSORT-ABS) were utilized to characterize tumor immune microenvironment (TIME). Finally, the mutation data were analyzed by employing "maftools" package. Results: In total, 136 DELs, 128 DEMis, and 2,028 DEMs were recognized in HCC. A specific lncRNA-miRNA-mRNA network consisting of 3 lncRNAs, 12 miRNAs, and 21 mRNAs was established. A ceRNA-based prognostic signature was established to classify samples into two risk subgroups, which presented excellent prognostic performance. In additional, prognostic risk-clinical nomogram was delineated to assess risk of individual sample quantitatively. Besides, risk score was significantly associated with contexture of TIME and immunotherapeutic targets. Finally, potential interaction between risk score with tumor mutation burden (TMB) was revealed. Conclusion: In this work, comprehensive analyses of ceRNAs coexpression network will facilitate prognostic prediction, delineate complexity of TIME, and contribute insight into precision therapy for HCC.

Project description:Our previous studies shown that syndecan-1 (SDC1) may be a novel class of biomarkers for the diagnosis and treatment of glioma, but its specific roles and the in-depth molecular mechanism remain elusive. Here, we used Estimation of STromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) algorithms and single-sample Gene Set Enrichment Analysis (ssGSEA) algorithms to evaluate the immune score of tumor samples and quantify the relative infiltration of immune cells in the tumor microenvironment (TME), respectively, in different data sets obtained from the Chinese Glioma Genome Atlas and The Cancer Gene Atlas. Next, we calculate the correlation of the immune score and immune cells with SDC1, respectively. To identify the specific process regulated by SDC1, the differentially expressed genes (DEGs) analysis between the high and low expression of SDC1 of glioma samples were used to discover the hub genes through Weighted Gene Coexpression Network Analysis (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed cardinal biological processes and pathways involved in genes and tumor grade correlation and survival analysis verified its significance in glioma. The results show that SDC1 is associated with the immune infiltration of glioma in the TME, especially activated CD4+T cells and CD8+T cells. The three data sets filter 8,887 DEGs, the genes in the blue modules were selected as hub genes in WGCNA. GO and KEGG analysis found eight genes in the blue modules involved in antigen processing and presentation in T cells in glioma. Kaplan-Meier estimator and log-rank test statistic determined that the introduced genes are associated with poor prognosis in glioma. Protein-protein network interaction analysis showed that SDC1 may regulate antigen processing and presentation through CTSL or CD4 in glioma. Finally, this study provided insights and clues for the next research direction of SDC1 and identified the key pathways and genes that might participate in the immune escape of glioma. These results might provide a new insight on the study of immune infiltration of glioma in the future.

Dataset Information

A Weighted Gene Co-Expression Network Analysis-Derived Prognostic Model for Predicting Prognosis and Immune Infiltration in Gastric Cancer.

Background

Methods

Results

Conclusion

Publications

A Weighted Gene Co-Expression Network Analysis-Derived Prognostic Model for Predicting Prognosis and Immune Infiltration in Gastric Cancer.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets