scholarly journals Investigation of Genetic Determinants of Glioma Immune Phenotype by Integrative Immunogenomic Scale Analysis

2021 ◽  
Vol 12 ◽  
Author(s):  
Binghao Zhao ◽  
Yuekun Wang ◽  
Yaning Wang ◽  
Congxin Dai ◽  
Yu Wang ◽  
...  

The immunosuppressive mechanisms of the surrounding microenvironment and distinct immunogenomic features in glioblastoma (GBM) have not been elucidated to date. To fill this gap, useful data were extracted from The Cancer Genome Atlas (TCGA), the Chinese Glioma Genome Atlas (CGGA), GSE16011, GSE43378, GSE23806, and GSE12907. With the ssGSEA method and the ESTIMATE and CIBERSORT algorithms, four microenvironmental signatures were used to identify glioma microenvironment genes, and the samples were reasonably classified into three immune phenotypes. The molecular and clinical features of these phenotypes were characterized via key gene set expression, tumor mutation burden, fraction of immune cell infiltration, and functional enrichment. Exhausted CD8+ T cell (GET) signature construction with the predictive response to commonly used antitumor drugs and peritumoral edema assisted in further characterizing the immune phenotype features. A total of 2,466 glioma samples with gene expression profiles were enrolled. Tumor purity, ESTIMATE, and immune and stromal scores served as the 4 microenvironment signatures used to classify gliomas into immune-high, immune-middle and immune-low groups, which had distinct immune heterogeneity and clinicopathological characteristics. The immune-H phenotype had higher expression of four immune signatures; however, most checkpoint molecules exhibited poor survival. Enriched pathways among the subtypes were related to immunity. The GET score was similar among the three phenotypes, while immune-L was more sensitive to bortezomib, cisplatin, docetaxel, lapatinib, and rapamycin prescriptions and displayed mild peritumor edema. The three novel immune phenotypes with distinct immunogenetic features could have utility for understanding glioma microenvironment regulation and determining prognosis. These results contribute to classifying glioma subtypes, remodeling the immunosuppressive microenvironment and informing novel cancer immunotherapy in the era of precision immuno-oncology.

2020 ◽  
Vol 10 ◽  
Author(s):  
Quanwei Zhou ◽  
Xuejun Yan ◽  
Weidong Liu ◽  
Wen Yin ◽  
Hongjuan Xu ◽  
...  

Diffuse glioma is one of the most prevalent malignancies of the brain, with high heterogeneity of tumor-infiltrating immune cells. However, immune-associated subtypes of diffuse glioma have not been determined, nor has the effect of different immune-associated subtypes on disease prognosis and immune infiltration of diffuse glioma patients. We retrieved the expression profiles of immune-related genes from The Cancer Genome Atlas (TCGA) (n = 672) and GSE16011 (n = 268) cohorts and used them to identify subtypes of diffuse glioma via Consensus Cluster Plus analysis. We used the limma, clusterProfiler, ESTIMATE, and survival packages of R for differential analysis, functional enrichment, immune and stromal score evaluation respectively in three subtypes, and performed log-rank tests in immune subtypes of diffuse glioma. The immune-associated features of diffuse glioma in the two cohorts were characterized via bioinformatic analyses of the mRNA expression data of immune-related genes. Three subtypes (C1–3) of diffuse glioma were identified from TCGA data, and were verified using the GSE16011 cohort. We then evaluated their immune characteristics and clinical features. Our mRNA profiling analyses indicated that the different subtypes of diffuse glioma presented differential expression profile of specific genes and signal pathways in the TCGA cohort. Patients with subtype C1, who were mostly diagnosed with grade IV glioma, had poorer outcomes than patients with subtype C2 or C3. Subtype C1 was characterized by a higher degree of immune cell infiltration as estimated by GSVA, and more frequent wildtype IDH1. By contrast, subtype C3 included more grade II and IDH1-mutated glioma, and was associated with more infiltration of CD4+T cells. Most subtype C2 had the features between subtypes C1 and C3. Meanwhile, immune checkpoints and their ligand molecules, including PD1/(PD-L1/PDL2), CTLA4/(CD80/CD86), and B7H3/TLT2, were significantly upregulated in subtype C1 and downregulated in subtype C3. In addition, patients with subtype C1 exhibited more frequent gene mutations. Univariate and multivariate Cox regression analyses revealed that diffuse glioma subtype was an effective, independent, and better prognostic factor. Therefore, we established a novel immune-related classification of diffuse glioma, which provides potential immunotherapy targets for diffuse glioma.


Cancers ◽  
2018 ◽  
Vol 10 (12) ◽  
pp. 494 ◽  
Author(s):  
Kyuryung Kim ◽  
Sora Jeon ◽  
Tae-Min Kim ◽  
Chan Jung

Papillary thyroid carcinoma (PTC) represents a heterogeneous disease with diverse clinical outcomes highlighting a need to identify robust biomarkers with clinical relevance. We applied non-negative matrix factorization-based deconvolution to publicly available gene expression profiles of thyroid cancers in the Cancer Genome Atlas (TCGA) consortium. Among three metagene signatures identified, two signatures were enriched in canonical BRAF-like and RAS-like thyroid cancers with up-regulation of genes involved in oxidative phosphorylation and cell adhesions, respectively. The third metagene signature representing up-regulation of immune-related genes further segregated BRAF-like and RAS-like PTCs into their respective subgroups of immunoreactive (IR) and immunodeficient (ID), respectively. BRAF-IR PTCs showed enrichment of tumor infiltrating immune cells, tall cell variant PTC, and shorter recurrence-free survival compared to BRAF-ID PTCs. RAS-IR and RAS-ID PTC subtypes included majority of normal thyroid tissues and follicular variant PTC, respectively. Immunopathological features of PTC subtypes such as immune cell fraction, repertoire of T cell receptors, cytolytic activity, and expression level of immune checkpoints such as and PD-L1 and CTLA-4 were consistently observed in two different cohorts. Taken together, an immune-related metagene signature can classify PTCs into four molecular subtypes, featuring the distinct histologic type, genetic and transcriptional alterations, and potential clinical significance.


2020 ◽  
Vol 18 (1) ◽  
Author(s):  
Qingyu Liang ◽  
Gefei Guan ◽  
Xue Li ◽  
Chunmi Wei ◽  
Jianqi Wu ◽  
...  

Abstract Background Molecular classification has laid the framework for exploring glioma biology and treatment strategies. Pro-neural to mesenchymal transition (PMT) of glioma is known to be associated with aggressive phenotypes, unfavorable prognosis, and treatment resistance. Recent studies have highlighted that long non-coding RNAs (lncRNAs) are key mediators in cancer mesenchymal transition. However, the relationship between lncRNAs and PMT in glioma has not been systematically investigated. Methods Gene expression profiles from The Cancer Genome Atlas (TCGA), the Chinese Glioma Genome Atlas (CGGA), GSE16011, and Rembrandt with available clinical and genomic information were used for analyses. Bioinformatics methods such as weighted gene co-expression network analysis (WGCNA), gene set enrichment analysis (GSEA), Cox analysis, and least absolute shrinkage and selection operator (LASSO) analysis were performed. Results According to PMT scores, we confirmed that PMT status was positively associated with risky behaviors and poor prognosis in glioma. The 149 PMT-related lncRNAs were identified by WGCNA analysis, among which 10 (LINC01057, TP73-AS1, AP000695.4, LINC01503, CRNDE, OSMR-AS1, SNHG18, AC145343.2, RP11-25K21.6, RP11-38L15.2) with significant prognostic value were further screened to construct a PMT-related lncRNA risk signature, which could divide cases into two groups with distinct prognoses. Multivariate Cox regression analyses indicated that the signature was an independent prognostic factor for high-grade glioma. High-risk cases were more likely to be classified as the mesenchymal subtype, which confers enhanced immunosuppressive status by recruiting macrophages, neutrophils, and regulatory T cells. Moreover, six lncRNAs of the signature could act as competing endogenous RNAs to promote PMT in glioblastoma. Conclusions We profiled PMT status in glioma and established a PMT-related 10-lncRNA signature for glioma that could independently predict glioma survival and trigger PMT, which enhanced immunosuppression.


Cells ◽  
2021 ◽  
Vol 10 (12) ◽  
pp. 3484
Author(s):  
Jisun Lim ◽  
Yeon Bi Han ◽  
Soo Young Park ◽  
Soyeon Ahn ◽  
Hyojin Kim ◽  
...  

Many studies support a stepwise continuum of morphologic changes between atypical adenomatous hyperplasia (AAH) and lung adenocarcinoma (ADC). Here we characterized gene expression patterns and the association of differentially expressed genes and immune tumor microenvironment behaviors in AAH to ADC during ADC development. Tumor tissues from nine patients with ADC and synchronous multiple ground glass nodules/lesions (GGN/Ls) were analyzed using RNA sequencing. Using clustering, we identified genes differentially and sequentially expressed in AAH and ADC compared to normal tissues. Functional enrichment analysis using gene ontology terms was performed, and the fraction of immune cell types was estimated. We identified up-regulated genes (ACSL5 and SERINC2) with a stepwise change of expression from AAH to ADC and validated those expressions by quantitative PCR and immunohistochemistry. The immune cell profiles revealed increased B cell activities and decreased natural killer cell activities in AAH and ADC. A stepwise change of differential expression during ADC development revealed potential effects on immune function in synchronous precursors and in tumor lesions in patients with lung cancer.


Diagnostics ◽  
2021 ◽  
Vol 11 (4) ◽  
pp. 726
Author(s):  
Hoang Dang Khoa Ta ◽  
Wan-Chun Tang ◽  
Nam Nhut Phan ◽  
Gangga Anuraga ◽  
Sz-Ying Hou ◽  
...  

Breast cancer (BRCA) is one of the most complex diseases and involves several biological processes. Members of the L-antigen (LAGE) family participate in the development of various cancers, but their expressions and prognostic values in breast cancer remain to be clarified. High-throughput methods for exploring disease progression mechanisms might play a pivotal role in the improvement of novel therapeutics. Therefore, gene expression profiles and clinical data of LAGE family members were acquired from the cBioportal database, followed by verification using the Oncomine and The Cancer Genome Atlas (TCGA) databases. In addition, the Kaplan-Meier method was applied to explore correlations between expressions of LAGE family members and prognoses of breast cancer patients. MetaCore, GlueGo, and GluePedia were used to comprehensively study the transcript expression signatures of LAGEs and their co-expressed genes together with LAGE-related signal transduction pathways in BRCA. The result indicated that higher LAGE3 messenger (m)RNA expressions were observed in BRCA tissues than in normal tissues, and they were also associated with the stage of BRCA patients. Kaplan-Meier plots showed that overexpression of LAGE1, LAGE2A, LAGE2B, and LAGE3 were highly correlated to poor survival in most types of breast cancer. Significant associations of LAGE family genes were correlated with the cell cycle, focal adhesion, and extracellular matrix (ECM) receptor interactions as indicated by functional enrichment analyses. Collectively, LAGE family members’ gene expression levels were related to adverse clinicopathological factors and prognoses of BRCA patients; therefore, LAGEs have the potential to serve as prognosticators of BRCA patients.


2021 ◽  
Vol 11 ◽  
Author(s):  
Qinglin Feng ◽  
Cheng Qian ◽  
Shibing Fan

Accumulating evidence suggests that hypoxia microenvironment and long non-coding lncRNAs (lncRNAs) exert critical roles in tumor development. Herein, we aim to develop a hypoxia-related lncRNA (HRL) model to predict the survival outcomes of patient with lower-grade glioma (LGG). The RNA-sequencing data of 505 LGG samples were acquired from The Cancer Genome Atlas (TCGA). Using consensus clustering based on the expression of hypoxia-related mRNAs, these samples were divided into three subsets that exhibit distinct hypoxia content, clinicopathologic features, and survival status. The differentially expressed lncRNAs across the subgroups were documented as candidate HRLs. With LASSO regression analysis, eight informative lncRNAs were selected for constructing the prognostic HRL model. This signature had a good performance in predicting LGG patients’ overall survival in the TCGA cohort, and similar results could be achieved in two validation cohorts from the Chinese Glioma Genome Atlas. The HRL model also showed correlations with important clinicopathologic characteristics such as patients’ age, tumor grade, IDH mutation, 1p/19q codeletion, MGMT methylation, and tumor progression risk. Functional enrichment analysis indicated that the HLR signature was mainly involved in regulation of inflammatory response, complement, hypoxia, Kras signaling, and apical junction. More importantly, the signature was related to immune cell infiltration, estimated immune score, tumor mutation burden, neoantigen load, and expressions of immune checkpoints and immunosuppressive cytokines. Finally, a nomogram was developed by integrating the HRL signature and clinicopathologic features, with a concordance index of 0.852 to estimate the survival probability of LGG patients. In conclusion, our study established an effective HRL model for prognosis assessment of LGG patients, which may provide insights for future research and facilitate the designing of individualized treatment.


PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e9038 ◽  
Author(s):  
Yixin Tian ◽  
Yiquan Ke ◽  
Yanxia Ma

Glioma is one of the most fatal tumors in central nervous system. Previous studies gradually revealed the association between tumor microenvironment and the prognosis of gliomas patients. However, the correlation between tumor-infiltrating immune cell and stromal signatures are unknown. In our study, we obtained gliomas samples from the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA). The landscape of tumor infiltrating immune cell subtypes in gliomas was calculated by CIBERSORT. As a result, we found high infiltration of macrophages was correlated with poor outcome (P < 0.05). Then functional enrichment analysis of high/low macrophage-infiltrating groups was performed by GSEA. The results showed three gene sets includes 102 core genes about angiogenesis were detected in high macrophage-infiltrating group. Next, we constructed PPI network and analyzed prognostic value of 102 core genes. We found that five stromal signatures indicated poor prognosis which including HSPG2, FOXF1, KDR, COL3A1, SRPX2 (P < 0.05). Five stromal signatures were adopted to construct a classifier. The classifier showed powerful predictive ability (AUC = 0.748). Patients with a high risk score showed poor survival. Finally, we validated this classifier in TCGA and the result was consistent with CGGA. Our investigation of tumor microenvironment in gliomas may stimulate the new strategy in immunotherapy. Five stromal signature correlated with poor prognosis also provide a strong predator of gliomas patient outcome.


Sign in / Sign up

Export Citation Format

Share Document