Dataset Information

Identification of Tumor Microenvironment Scoring Scheme Based on Bioinformatics Analysis of Immune Cell Infiltration Pattern of Ovarian Cancer

ABSTRACT:

Background

Tumor microenvironment (TME) is the crucial mediator of tumor progression, and the TME model based on immune cell infiltration to characterize ovarian cancer is considered to be a promising strategy.

Methods

Sample data of three ovarian cancer cohorts were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The scores of 22 kinds of immune cells were calculated based on CIBERSORT, and the TME clusters (TMECs) of ovarian cancer was determined by ConsensusClusterPlus. Genomic subtype was identified by non-negative matrix factorization (NMF). A TME scoring scheme was constructed using k-means algorithm and principal component analysis (PCA) to quantify the TME infiltration pattern of individuals.

Results

Four TME subtypes of ovarian cancer samples were defined: TMEC1, TMEC2, TMEC3, and TMEC4. There were also significant differences in overall survival (OS) among the four TMEC, and the OS of TMEC3 was the longest. The difference analysis of TMEC3 and the other three TMECs respectively identified the DEGs and took the intersection, and 585 DEGs were obtained. Two genomic subtypes were identified by NMF based on the expression of 585 genes, which were called GeneC1 and GeneC2. The TME scoring scheme constructed by k-means and PCA algorithm was used to calculate the TME score of ovarian cancer in TCGA. High-TME score was significantly correlated with shorter survival time, older age, lower immunoactivated molecules, and immune checkpoint gene expression.

Conclusions

This study highlighted the complexity and diversity of TME infiltration patterns in ovarian cancer and constructed a set of TME scoring scheme to reveal TME infiltration patterns and provided new insights into the landscape of TME.

SUBMITTER: Huo X

PROVIDER: S-EPMC9448528 | biostudies-literature | 2022 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Objective: Because of the modest immunotherapeutic response among ovarian carcinoma (OC) patients, it is significant to evaluate antitumor immune response and develop more effective precision immunotherapeutic regimens. Here, this study aimed to determine diverse immune subtypes of OC. Methods: This study curated the expression profiles of prognostic immunologically relevant genes and conducted consensus clustering analyses for determining immune subtypes among OC patients in TCGA cohort. With Boruta algorithm, characteristic genes were screened for conducting an immune scoring system through principal component analysis algorithm. The single-sample gene set enrichment analysis and ESTIAMTE methods were adopted for quantifying the immune infiltrates and responses to chemotherapeutic agents were estimated with pRRophetic algorithm. Two immunotherapeutic cohorts were used for investigating the efficacy of immune score in predicting therapeutic benefits. Results: Two immune subtypes were conducted among 377 OC patients. Immune subtype 2 was characterized by worse clinical prognosis, more frequent genetic variations and mutations, enhanced immune infiltrates, and increased expression of MHC molecules and programmed cell death protein 1/programmed death ligand 1 (PD-1/PD-L1). In total, 30 prognosis-relevant characteristic immune subtype-derived genes were identified for constructing the immune score of OC patients. High immune score was linked with more dismal prognosis, decreased immune infiltrations, and expression of MHC molecules. High immune score presented favorable sensitivity to doxorubicin and vinorelbine and reduced sensitivity to cisplatin. In addition, immune score possessed the potential in predicting benefits from anti-PD-1/PD-L1 therapy. Conclusion: Collectively, our findings propose two complex and diverse immune subtypes of OC. Quantitative assessment of immune subtypes in individual patients strengthens the understanding of tumor microenvironment features and promotes more effective immunotherapeutic regimens.

Project description:BackgroundIncreasing evdence supports the suggestion that the immune cell infiltration (ICI) patterns play a pivotal role in tumor progression in breast cancer (BRCA). Nonetheless, there has been no comprehensive analysis of the ICI patterns effects on the clinical outcomes and immunotherapy.MethodsMultiomic data for BRCA samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method, and CIBERSORT analysis were used to uncover the landscape of the tumor immune microenvironment (TIME). BRCA subtypes based on the ICI pattern were identified by consensus clustering and principal-component analysis was performed to obtain the ICI scores to quantify the ICI patterns in individual tumors. Their prognostic value was validated by the Kaplan-Meier survival curves. Gene set enrichment analysis (GSEA) was applied for functional annotation. Immunophenoscore (IPS) was employed to explore the immunotherapeutic role of the ICI scores. Finally, the mutation data was analyzed by using the "maftools" R package.ResultsThree different immune infiltration patterns with a distinct prognosis and biological signature were recognized among 1,198 BRCA samples. The characteristics of TIME under these three patterns were highly consistent with three known immune profiles: immune- excluded, immune-desert, and immune-inflamed phenotypes, respectively. The identification of the ICI patterns within individual tumors based on the ICI score, developed under the ICI-related signature genes, contributed into dissecting biological processes, clinical outcome, immune cells infiltration, immunotherapeutic effect, and genetic variation. High ICI score subtype, characterized with a suppression of immunity, suggested an immune-exhausted phenotype. Abundant effective immune cells were discovered in the low ICI score patients, which corresponded to an immune-activated phenotype and might present an immunotherapeutic advantage. Immunophenoscore was implemented as a surrogate of immunotherapeutic outcome, low-ICI scores samples obtained a significantly higher immunophenoscore. Enrichment of the JAK/STAT and VEGF signal pathways were activated in the ICI low-score subgroup. Finally, the synergistic effect between the ICI score and the tumor mutation burden (TMB) was confirmed.ConclusionThis work comprehensively elucidated that the ICI patterns served as an indispensable player in complexity and diversity of TIME. Quantitative identification of the ICI patterns in individual tumor will contribute into mapping the landscape of TIME further optimizing precision immunotherapy.

Project description:ObjectiveGastric carcinoma (GC) is the fifth most commonly diagnosed cancer and the third leading cause of cancer-related deaths globally. The tumor microenvironment plays a significant role in the pathogenesis, prognosis, and response to immunotherapy. However, the immune-related molecular mechanisms underlying GC remain elusive. Bioinformatics analysis of the gene expression of GC and paracancerous healthy tissues from the same patient was performed to identify the key genes and signaling pathways, as well as their correlation to the infiltration of the tumor microenvironment (TME) by various immune cells related to GC development.MethodsWe employed GSE19826, a gene expression profile from the Gene Expression Omnibus (GEO), for our analysis. Functional enrichment analysis of Differentially Expressed Genes (DEGs) was conducted using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes database.ResultsCytoscape software facilitated the identification of nine hub DEGs, namely, FN1, COL1A1, COL1A2, THBS2, COL3A1, COL5A1, APOE, SPP1, and BGN. Various network analysis algorithms were applied to determine their high connectivity. Among these hub genes, FN1, COL1A2, THBS2, COL3A1, COL5A1, and BGN were found to be associated with a poor prognosis for GC patients. Subsequent analysis using the TIMER database revealed the infiltration status of the TME concerning the overexpression of these six genes. Specifically, the abovementioned genes demonstrated direct correlations with cancer-associated fibroblasts, M1 and M2 macrophages, myeloid-derived suppressor cells, and activated dendritic cells.ConclusionOur findings suggest that the identified hub genes, particularly BGN, FN1, COL1A2, THBS2, COL3A1, and COL5A1, play crucial roles in GC prognosis and TME cell infiltration. This comprehensive analysis enhances our understanding of the molecular mechanisms underlying GC development and may contribute to the identification of potential therapeutic targets and prognostic markers for GC patients.

Project description:Ovarian cancer (OV) is a deadly gynecological cancer. The tumor immune microenvironment (TIME) plays a pivotal role in OV development. However, the TIME of OV is not fully known. Therefore, we aimed to provide a comprehensive network of the TIME in OV. Gene expression data and clinical information from OV patients were obtained from the Cancer Genome Atlas Program (TCGA) database. Non-negative Matrix Factorization, NMFConsensus, and nearest template prediction algorithms were used to perform molecular clustering. The biological functions of differentially expressed genes (DEGs) were identified using Metascape, gene set enrichment analysis (GSEA), gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and tumor mutation burden were analyzed using Gistic 2.0, R package maftools, and TCGA mutations, respectively. Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data and CIBERSORT were utilized to elucidate the TIME. Moreover, external data from the International Cancer Genome Consortium (ICGC) and ArrayExpress databases were used to validate the signature. All 361 samples from the TCGA OV dataset were classified into Immune Class and non-Immune Class with immune signatures. By comparing the two classes, we identified 740 DEGs that accumulated in immune-related, cancer-related, inflammation-related biological functions and pathways. There were significant differences in the CNVs between the Immune and non-Immune Classes. The Immune Class was further divided into immune-activated and immune-suppressed subtypes. There was no significant difference in the top 20 genes in somatic SNPs among the three groups. In addition, the immune-activated subtype had significantly increased proportions of CD4 memory resting T cells, T cells, M1 macrophages, and M2 macrophages than the other two groups. The qRT-PCR results indicated that the mRNA expression levels of RYR2, FAT3, MDN1 and RYR1 were significantly down-regulated in OV compared with normal tissues. Moreover, the signatures of the TIME were validated using ICGC cohort and the ArrayExpress cohort. Our study clustered the OV patients into an immune-activated subtype, immune-suppressed subtype, and non-Immune Class and provided potential clues for further research on the molecular mechanisms and immunotherapy strategies of OV.