scholarly journals A Six-lncRNA Signature for Immunophenotype Prediction of Glioblastoma Multiforme

2021 ◽  
Vol 11 ◽  
Author(s):  
Ming Gao ◽  
Xinzhuang Wang ◽  
Dayong Han ◽  
Enzhou Lu ◽  
Jian Zhang ◽  
...  

Glioblastoma multiforme (GBM) is the most aggressive primary tumor of the central nervous system. As biomedicine advances, the researcher has found the development of GBM is closely related to immunity. In this study, we evaluated the GBM tumor immunoreactivity and defined the Immune-High (IH) and Immune-Low (IL) immunophenotypes using transcriptome data from 144 tumors profiled by The Cancer Genome Atlas (TCGA) project based on the single-sample gene set enrichment analysis (ssGSEA) of five immune expression signatures (IFN-γ response, macrophages, lymphocyte infiltration, TGF-β response, and wound healing). Next, we identified six immunophenotype-related long non-coding RNA biomarkers (im-lncRNAs, USP30-AS1, HCP5, PSMB8-AS1, AL133264.2, LINC01684, and LINC01506) by employing a machine learning computational framework combining minimum redundancy maximum relevance algorithm (mRMR) and random forest model. Moreover, the expression level of identified im-lncRNAs was converted into an im-lncScore using the normalized principal component analysis. The im-lncScore showed a promising performance for distinguishing the GBM immunophenotypes with an area under the curve (AUC) of 0.928. Furthermore, the im-lncRNAs were also closely associated with the levels of tumor immune cell infiltration in GBM. In summary, the im-lncRNA signature had important clinical implications for tumor immunophenotyping and guiding immunotherapy in glioblastoma patients in future.

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yangming Hou ◽  
Yingjuan Xu ◽  
Dequan Wu

AbstractThe infiltration degree of immune and stromal cells has been shown clinically significant in tumor microenvironment (TME). However, the utility of stromal and immune components in Gastric cancer (GC) has not been investigated in detail. In the present study, ESTIMATE and CIBERSORT algorithms were applied to calculate the immune/stromal scores and the proportion of tumor-infiltrating immune cell (TIC) in GC cohort, including 415 cases from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were screened by Cox proportional hazard regression analysis and protein–protein interaction (PPI) network construction. Then ADAMTS12 was regarded as one of the most predictive factors. Further analysis showed that ADAMTS12 expression was significantly higher in tumor samples and correlated with poor prognosis. Gene Set Enrichment Analysis (GSEA) indicated that in high ADAMTS12 expression group gene sets were mainly enriched in cancer and immune-related activities. In the low ADAMTS12 expression group, the genes were enriched in the oxidative phosphorylation pathway. CIBERSORT analysis for the proportion of TICs revealed that ADAMTS12 expression was positively correlated with Macrophages M0/M1/M2 and negatively correlated with T cells follicular helper. Therefore, ADAMTS12 might be a tumor promoter and responsible for TME status and tumor energy metabolic conversion.


2020 ◽  
Vol 22 (Supplement_2) ◽  
pp. ii28-ii28
Author(s):  
Alvaro Alvarado ◽  
Kaleab Tessema ◽  
Kunal Patel ◽  
Riki Kawaguchi ◽  
Richard Everson ◽  
...  

Abstract Despite efforts to gain a deeper understanding of its molecular architecture, glioblastoma (GBM) remains uniformly fatal. While genome-based molecular subtyping has revealed that GBMs may be parsed into several molecularly distinct categories, this insight has yielded little progress towards extending patient survival. In particular, the great phenotypic heterogeneity of GBM – both inter and intratumorally – has hindered therapeutic efforts. To this end, we interrogated tumor samples using a pathway-based approach to resolve tumoral heterogeneity. Gene set enrichment analysis (GSEA) was applied to gene expression data and used to provide an overview of each sample that can be compared to other samples by generating sample clusters based on overall patterns of enrichment. The Cancer Genome Atlas (TCGA) samples were clustered using the canonical and oncogenic signatures and in both cases the clustering was distinct from the molecular subtype previously reported and clusters were informative of patient survival. We also analyzed single cell RNA sequencing datasets and uniformly found two clusters of cells enriched for cell cycle regulation and survival pathways. We have validated our approach by generating gene lists from common elements found in the top contributing genesets for a particular cluster and testing the top targets in appropriate gliomasphere patient-derived lines. Samples enriched for cell cycle related genesets showed a decrease in sphere formation capacity when E2F1, out top target, was silenced and when treated with fulvestrant and calcitriol, which were identified as potential drugs targeting this genelist. Conversely, no changes were observed in samples not enriched for this gene list. Finally, we interrogated spatial heterogeneity and found higher enrichment of the proliferative signature in contrast enhancing compared with non-enhancing regions. Our studies relate inter- and intratumoral heterogeneity to critical cellular pathways dysregulated in GBM, with the ultimate goal of establishing a pipeline for patient- and tumor-specific precision medicine.


2021 ◽  
Vol 20 ◽  
pp. 153303382199208
Author(s):  
Shufang Wang ◽  
Xinlong Huo

Background: Estrogen-related receptor alpha (ESRRA) was reported to play an important role in multiple biological processes of neoplastic diseases. The roles of ESRRA in endometrial cancer have not been fully investigated yet. Methods: Expression data and clinicopathological data of patients with uteri corpus endometrial carcinoma (UCEC) were obtained from The Cancer Genome Atlas (TCGA). Comprehensive bioinformatics analysis was performed, including receiver operating characteristics (ROC) curve analysis, Kaplan-Meier survival analysis, gene ontology (GO) enrichment analysis, and Gene Set Enrichment Analysis (GSEA). Immunohistochemistry was used to detect the protein expression level of ESRRA and CCK-8 assay was performed to evaluate the effect of ESRRA on the proliferation ability. Results: A total of 552 UCEC tissues and 35 normal tissues were obtained from the TCGA database. The mRNA and protein expression level of ESRRA was highly elevated in UCEC compared with normal tissues, and was closely associated with poor prognosis. ROC analysis indicated a very high diagnostic value of ESRRA for patients with UCEC. GO and GSEA functional analysis showed that ESRRA might be mainly involved in cellular metabolism processes, in turn, tumorigenesis and progression of UCEC. Knockdown of ESRRA inhibited the proliferation of UCEC cells in vitro. Further immune cell infiltration demonstrated that ESRRA enhanced the infiltration level of neutrophil cell and reduced that of T cell (CD4+ naïve), NK cell, and cancer associated fibroblast (CAF). The alteration of immune microenvironment will greatly help in developing immune checkpoint therapy for UCEC. Conclusions: Our study comprehensively analyzed the expression level, clinical value, and possible mechanisms of action of ESRRA in UCEC. These findings showed that ESRRA might be a potential diagnostic and therapeutic target.


2021 ◽  
Vol 11 ◽  
Author(s):  
Qiming Wang ◽  
Yan Cai ◽  
Xuewen Fu ◽  
Liang Chen

In recent years, the incidence and the mortality rate of cervical cancer have been gradually increasing, becoming one of the major causes of cancer-related death in women. In particular, patients with advanced and recurrent cervical cancers present a very poor prognosis. In addition, the vast majority of cervical cancer cases are caused by human papillomavirus (HPV) infection, of which HPV16 infection is the main cause and squamous cell carcinoma is the main presenting type. In this study, we performed screening of differentially expressed genes (DEGs) based on The Cancer Genome Atlas (TCGA) database and GSE6791, constructed a protein–protein interaction (PPI) network to screen 34 hub genes, filtered to the remaining 10 genes using the CytoHubba plug-in, and used survival analysis to determine that RPS27A was most associated with the prognosis of cervical cancer patients and has prognostic and predictive value for cervical cancer. The most significant biological functions and pathways of RPS27A enrichment were subsequently investigated with gene set enrichment analysis (GSEA), and integration of TCGA and GTEx database analyses revealed that RPS27A was significantly expressed in most cancer types. In this study, our analysis revealed that RPS27A can be used as a prognostic biomarker for HPV16 cervical cancer and has biological significance for the growth of cervical cancer cells.


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.


2019 ◽  
Author(s):  
rui kong ◽  
Nan Wang ◽  
Wei Han ◽  
Yuejuan Zheng ◽  
Jie Lu

Abstract Background: In recent years, long non-coding RNAs (lncRNAs) are emerging as crucial regulators in the immunological process of liver hepatocellular carcinoma (LIHC). Increasing studies have found that some lncRNAs could be used as a diagnostic or therapeutic target for clinical management, but little research has investigated the role of immune-related lncRNA in tumor prognosis. In this study, we aimed to develop an immune lncRNA signature for the precise diagnosis and prognosis of liver hepatocellular carcinoma. Methods: Gene expression profiles of LIHC samples obtained from TCGA were screened for immune-related genes using two reference gene sets. The optimal immune-related lncRNA signature was built via correlational analysis, univariate and multivariate cox analysis. Then the Kaplan-Meier plot, ROC curve, clinical analysis, gene set enrichment analysis, and principal component analysis were carried out to evaluate the capability of immune lncRNA signature as a prognostic indicator. Results: Six long non-coding RNA MSC−AS1, AC009005.1, AL117336.3, AL031985.3, AL365203.2, AC099850.3 were identified via correlation analysis and cox regression analysis considering their interactions with immune genes. Next, tumor samples were separated into two risk groups by the signature with different clinical outcomes. Stratification analysis showed the prognostic ability of this signature acted as an independent factor. The AUC value of ROC curve was 0.779. The Kaplan-Meier method was used in survival analysis and results showed a statistical difference between the two risk groups. The predictive performance of this signature was validated by principal component analysis (PCA). Data from gene set enrichment analysis (GSEA) further unveiled several potential biological processes of these biomarkers may involve in. Conclusion: In summary, the study demonstrated the potential role of the six-lncRNA signature served as an independent prognostic factor for LIHC patients.


2021 ◽  
Author(s):  
kai wang ◽  
Jun xing Feng ◽  
Zhi ling Zheng ◽  
Ying ze Chai ◽  
Hui jun Yu ◽  
...  

Abstract Background: Transient receptor potential cation channel subfamily V member 4 (TRPV4) has been reported to regulate tumor progression in many tumor types. However, its association with the tumor immune microenvironment remains unclear.Methods: TRPV4 expression was assessed using data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database. The clinical features and prognostic roles of TRPV4 were assessed using TCGA cohort. Gene set enrichment analysis (GSEA) of TRPV4 was conducted using the R package clusterProfiler. We analyzed the association between TRPV4 and immune cell infiltration scores of TCGA samples downloaded from published articles and the TIMER2 database.Results: TRPV4 was highly expressed and associated with worse overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI) in colon adenocarcinoma (COAD) and ovarian cancer. Furthermore, TRPV4 expression was closely associated with immune regulation-related pathways. Moreover, tumor-associated macrophage (TAM) infiltration levels were positively correlated with TRPV4 expression in TCGA pan-cancer samples. Immunosuppressive genes such as PD-L1, PD-1, CTLA4, LAG3, TIGIT, TGFB1, and TGFBR1 were positively correlated with TRPV4 expression in most tumors.Conclusions: Our results suggest that TRPV4 is an oncogene and a prognostic marker in COAD and ovarian cancer. High TRPV4 expression is associated with tumor immunosuppressive status and may contribute to TAM infiltration based on TCGA data from pan-cancer samples.


2021 ◽  
Vol 19 (1) ◽  
pp. 169-190
Author(s):  
Peiyuan Li ◽  
◽  
Gangjie Qiao ◽  
Jian Lu ◽  
Wenbin Ji ◽  
...  

<abstract> <p>Plasmacytoma variant translocation 1 (PVT1) is involved in multiple signaling pathways and plays an important regulatory role in a variety of malignant tumors. However, its role in the prognosis and immune invasion of bladder urothelial carcinoma (BLCA) remains unclear. This study investigated the expression of PVT1 in tumor tissue and its relationship with immune invasion, and determined its prognostic role in patients with BLCA. Patients were identified from the cancer genome atlas (TCGA). The enrichment pathway and function of PVT1 were explained by gene ontology (GO) term analysis, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA), and the degree of immune cell infiltration was quantified. Kaplan–Meier analysis and Cox regression were used to analyze the correlation between PVT1 and survival rate. PVT1-high BLCA patients had a lower 10-year disease-specific survival (DSS P &lt; 0.05) and overall survival (OS P &lt; 0.05). Multivariate Cox regression analysis showed that PVT1 (high vs. low) (P = 0.004) was an independent prognostic factor. A nomogram was used to predict the effect of PVT1 on the prognosis. PVT1 plays an important role in the progression and prognosis of BLCA and can be used as a medium biomarker to predict survival after cystectomy.</p> </abstract>


2021 ◽  
Vol 11 ◽  
Author(s):  
Yuan Nie ◽  
Mei-chun Jiang ◽  
Cong Liu ◽  
Qi Liu ◽  
Xuan Zhu

BackgroundsTumor microenvironment (TME) plays a crucial role in the initiation and progression of Hepatocellular Carcinoma (HCC), especially immune infiltrates. However, there is still a challenge in understanding the modulation of the immune and stromal components in TME, especially TME related genes.MethodsThe proportion of tumor-infiltrating immune cells (TICs) and the immune and stromal scores in 374 HCC patients from The Cancer Genome Atlas (TCGA) database were determined using CIBERSORT and ESTIMATE computational methods. The final screened genes were confirmed by the PPI network and univariate Cox regression of the differentially expressed genes based on different immune or stromal scores. The correlation between the expression levels of the final gene interactions and the clinical characteristics was based on TCGA database and local hospital data. Gene set enrichment analysis (GSEA) and the effect of CXCL5 expression on TICs were conducted.ResultsThere were correlations between the expression of CXCL5 and survival of HCC patients and TMN classification both in TCGA database and local hospital data. The immune-related activities were enriched in the high-expression group; however, the metabolic pathways were enriched in the low-expression group. The result of CIBERSORT analyzing had indicated that CXCL5 expression were correlated with the proportion of NK cells activated, macrophages M0, Mast cells resting, Neutrophils.ConclusionsCXCL5 was a potential prognostic marker for HCC and provides clues regarding immune infiltrates, which offers extra insight for therapeutics of HCC, however, more independent cohorts and functional experiments of CXCL5 are warranted.


2021 ◽  
Vol 2021 ◽  
pp. 1-17
Author(s):  
Fei Liu ◽  
Xiaopeng Yu ◽  
Guijin He

Background. We analyzed the n6-methyladenosine (m6A) modification patterns of immune cells infiltrating the tumor microenvironment of breast cancer (BC) to provide a new perspective for the early diagnosis and treatment of BC. Methods. Based on 23 m6A regulatory factors, we identified m6A-related gene characteristics and m6A modification patterns in BC through unsupervised cluster analysis. To examine the differences in biological processes among various m6A modification modes, we performed genomic variation analysis. We then quantified the relative infiltration levels of different immune cell subpopulations in the tumor microenvironment of BC using the CIBERSORT algorithm and single-sample gene set enrichment analysis. Univariate Cox analysis was used to screen for m6A characteristic genes related to prognosis. Finally, we evaluated the m6A modification pattern of patients with a single BC by constructing the m6Ascore based on principal component analysis. Results. We identified three different m6A modification patterns in 2128 BC samples. A higher abundance of the immune infiltration of the m6Acluster C was indicated by the results of CIBERSORT and the single-sample gene set enrichment analysis. Based on the m6A characteristic genes obtained through screening, the m6Ascore was determined. The BC patients were segregated into m6Ascore groups of low and high categories, which revealed significant survival benefits among patients with low m6Ascores. Additionally, the high-m6Ascore group had a higher mutation frequency and was associated with low PD-L1 expression, and the m6Ascore and tumor mutation burden showed a positive correlation. In addition, treatment effects were better in patients in the high-m6Ascore group. Conclusions. In case of a single patient with BC, the immune cell infiltration characteristics of the tumor microenvironment and the m6A methylation modification pattern could be evaluated using the m6Ascore. Our results provide a foundation for improving personalized immunotherapy of BC.


Sign in / Sign up

Export Citation Format

Share Document