scholarly journals Comprehensive Analysis of m6A Regulators Characterized by the Immune Cell Infiltration in Head and Neck Squamous Cell Carcinoma to Aid Immunotherapy and Chemotherapy

2021 ◽  
Vol 11 ◽  
Author(s):  
Zhiqiang Yang ◽  
Xiaoping Ming ◽  
Shuo Huang ◽  
Minlan Yang ◽  
Xuhong Zhou ◽  
...  

BackgroundN6-Methyladenosine (m6A), which is a prevalent regulator of mRNA expression, has gathered increasing study interests. Though the role of m6A as being important in many biological processes (such as growth and proliferation of cancers) has been well documented, its potential role in tumor immune microenvironment (TIME) has rarely been analyzed.MethodsWe downloaded RNA expression, single nucleotide polymorphism (SNP), and copy number variation (CNV) data from The Cancer Genome Atlas (TCGA). We then curated 21 m6A regulators and clustered patients into three m6A subtypes and m6A-related gene subtypes and compared them based on overall survival (OS). The combination of CIBERSORT as well as ssGSEA quantified the infiltration levels of immune cells and immune-related functions. The m6A scores were determined by using principal component analysis (PCA) algorithm. Furthermore, we evaluate the correlation of m6A regulators with immune and response to therapy.ResultsThree m6A clusters were identified based on the TCGA-HNSCC cohort, and there were significant associations among them in overall outcomes and caner-related pathways. We found that three m6A clusters were consistent with three phenotypes: immune-inflamed, immune-dessert, and immune-excluded. HNSCC patients were divided into high– and low–m6A score groups based on the cutoff of m6A score. Patients with lower m6A score had better overall survival outcome. Further analysis indicated that patients with higher m6A score presented higher tumor mutation burden (TMB). In addition, patients in low–m6A score subgroup had high chemotherapeutics sensitivity. GEO cohort confirmed patients with low m6A score demonstrated significant overall survival advantages and clinical benefits. Low m6A score carry an increased neoantigen load, eliciting a response to immunotherapy, and its value in predicting survival outcomes of immunotherapy was also confirmed in three anti-PD-1 cohorts.ConclusionsOur study demonstrated that m6A regulators are closely related to TIME and the m6A score was an effective prognostic biomarker and predictive indicator for immunotherapy and chemotherapeutics. Comprehensive evaluation of m6A regulators in tumors will extend our understanding of TIME and effectively guide increasing study investigations on immunotherapy and chemotherapy strategies for HNSCC.

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.


2020 ◽  
Author(s):  
Wei Li ◽  
Haitao Xiao ◽  
Xuewen Xu ◽  
Yange Zhang

Abstract Background: Tumor-infiltrating immune cells were demonstrated to be associated with patient survival and responses of targeted therapy. However, in melanoma, there is no reliable and individualized prognostic signatures based on comprehensive evaluation of the immune profile inferred from bulk tumor transcriptomes. In this study, we aimed to develop immunoscores associated with prognosis and responses of anti-PD1 targeted therapy in stage III-IV melanoma.Methods: The immunoscore and immunoscore-based prognostic nomogram were constructed based on the melanoma cohort from the Cancer Genome Atlas (TCGA), and validated in the population from the Gene Expression Omnibus (GEO). Besides, in another cohort obtained from the GEO database, we developed an immunoscore for predicting responses of anti-PD1 therapy. Twenty-two types of immune cell fraction were estimated using CIBERSORT. The least absolute shrinkage and selection operator (Lasso) regression model was utilized to develop individualized immunoscores.Results: With the Lasso regression, an immunoscore was constructed consisting of nine types of immune cell subtypes. In both of the training (192 cases) and validation (227 cases) cohorts, significant difference was observed between immunoscore-low and immunoscore-high groups in overall survival (OS). Multivariable analysis demonstrated that the immunoscore was an independent prognostic factor (P < 0.001) for OS. The prognostic value of the immunoscore was also confirmed by the ROC curves. Nomogram integrating immunotype and other clinical characteristics also showed good discrimination, calibration and usability in both of the training and validation cohorts. Finally, in another GEO cohort (218 cases), an immunoscore was constructed based on nine immune cell types for predicting anti-PD1 therapy response. Conclusion: The proposed immunoscores represent promising models for estimating overall survival and anti-PD1 treatment response in patients with stage III-IV melanoma.


Molecules ◽  
2021 ◽  
Vol 26 (18) ◽  
pp. 5659
Author(s):  
Rebecca M. Rodriguez ◽  
Mark Menor ◽  
Brenda Y. Hernandez ◽  
Youping Deng ◽  
Vedbar S. Khadka

One in five cancers is attributed to infectious agents, and the extent of the impact on the initiation, progression, and disease outcomes may be underestimated. Infection-associated cancers are commonly attributed to viral, and to a lesser extent, parasitic and bacterial etiologies. There is growing evidence that microbial community variation rather than a single agent can influence cancer development, progression, response to therapy, and outcome. We evaluated microbial sequences from a subset of infection-associated cancers—namely, head and neck squamous cell carcinoma (HNSC), liver hepatocellular carcinoma (LIHC), and stomach adenocarcinoma (STAD) from The Cancer Genome Atlas (TCGA). A total of 470 paired tumor and adjacent normal samples were analyzed. In STAD, concurrent presence of EBV and Selemonas sputigena with a high diversity index were associated with poorer survival (HR: 2.23, 95% CI 1.26–3.94, p = 0.006 and HR: 2.31, 95%CI 1.1–4.9, p = 0.03, respectively). In LIHC, lower microbial diversity was associated with poorer overall survival (HR: 2.57, 95%CI: 1.2, 5.5, p = 0.14). Bacterial within-sample diversity correlates with overall survival in infection-associated cancers in a subset of TCGA cohorts.


2021 ◽  
Author(s):  
Zaoqu Liu ◽  
Long Liu ◽  
Chunguang Guo ◽  
Libo Wang ◽  
Zhaonan Li ◽  
...  

Abstract BackgroundEsophageal adenocarcinoma (EAC) remains a leading cause of cancer-related deaths worldwide, and demonstrates a predominant rising incidence in Western countries. Recently, immunotherapy has dramatically changed the landscape of treatment for many advanced cancers, the benefit in EAC thus far been limited to a small fraction of patients. MethodsUsing somatic mutations data of The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC), we delineated somatic mutation landscape of EAC patients from US and England. Bioinformatics algorithms were utilized to perform function annotation, immune cell infiltration analysis, and immunotherapy response assessment.ResultsWe found that RYR2 was a common frequently mutated gene (FMG) in both cohorts, and patients with RYR2 mutation suggested higher tumor mutation burden (TMB), better prognosis, and superior expression of immune checkpoints. Moreover, RYR2 mutation upregulated the signaling pathways implicated in immune response and enhanced antitumor immunity in EAC. Multiple bioinformatics algorithms for assessing immunotherapy response demonstrated that patients with RYR2 mutation might benefit more from immunotherapy. In order to provide additional reference for antitumor therapy of different RYR2 status, we identified nine latent antitumor drugs associated with RYR2 status in EAC. ConclusionsThis study reveals a novel gene whose mutation could be served as a potential biomarker for prognosis, TMB, and immunotherapy of EAC patients.


Author(s):  
Fan Kou ◽  
Lei Wu ◽  
Ye Zhu ◽  
Baihui Li ◽  
Ziqi Huang ◽  
...  

AbstractSomatic copy number alterations (SCNA), which are widespread in cancer, can predict the efficacy of immune checkpoint inhibitors in non-small-cell lung cancer (NSCLC). However, the usefulness of SCNA for predicting the survival of patients treated with cytokine-induced killer (CIK) cells or chemotherapy (CT) is unknown. This study aimed to explore the correlation between SCNA and clinical outcome in NSCLC patients treated with CIK + CT or CT alone. We performed whole-exome sequencing on 45 NSCLC patients treated with CIK + CT, as well as 305 NSCLC patients treated with CT alone, from The Cancer Genome Atlas, which showed SCNA had a superiority in predicting the progression-free survival (PFS) over tumor mutation burden (TMB) and SCNA + TMB in NSCLC patients treated with CIK + CT, especially in lung adenocarcinoma, while SCNA could not predict the efficacy of CT alone. Additionally, we investigated the association between SCNA and immune cell infiltration by RNA sequencing and immunohistochemistry. The results revealed that SCNA was negatively associated with the expression of dendritic cells. Collectively, this study revealed a negative correlation between SCNA and response to CIK + CT and showed that SCNA is a predictive indicator in LUAD patients treated with CIK + CT.


Author(s):  
Yihong Luo ◽  
Xiang Sun ◽  
Jian Xiong

Introduction: Studies have demonstrated the epigenetic regulation of immune responses in various cancers. However, little is known about the RNA N6-methyladenosine (m6A) modification patterns of the microenvironment (TME) cell infiltration in ovarian cancer (OC).Methods: We evaluated the correlation between m6A modification patterns and TME cell infiltration based on 459 OC samples from the Cancer Genome Atlas and Gene-Expression Omnibus database. We constructed an m6Ascore system to quantify m6A modification patterns using principal component analysis.Results: Based on unsupervised clustering, three m6A modification patterns were identified. Gene set variation analysis showed that the antigen presentation signal pathway, the NOTCH signaling pathway, and the metabolism-related pathway differed significantly across m6A modificaiton patterns. The m6Ascore is closely correlated with TME cell infiltration. OC patients with lower m6Ascores had worse outcomes. There was better risk stratification with combined m6Ascore and tumor mutation burden. The responses to immune checkpoint inhibitor treatment significantly differed between high and low m6Ascore groups.Conclusion: M6A modification plays an essential role in TME cell infiltration in OC. Evaluating the m6A modification patterns in OC patients could enhance our understanding of TME infiltration characterization and guide immunotherapy strategies.


Epigenomics ◽  
2020 ◽  
Author(s):  
Qijie Zhao ◽  
Jinan Guo ◽  
Yueshui Zhao ◽  
Jing Shen ◽  
Parham Jabbarzadeh Kaboli ◽  
...  

Background: PD-L1 and PD-L2 are ligands of PD-1. Their overexpression has been reported in different cancers. However, the underlying mechanism of PD-L1 and PD-L2 dysregulation and their related signaling pathways are still unclear in gastrointestinal cancers. Materials & methods: The expression of PD-L1 and PD-L2 were studied in The Cancer Genome Atlas and Genotype-Tissue Expression databases. The gene and protein alteration of PD-L1 and PD-L2 were analyzed in cBioportal. The direct transcription factor regulating PD-L1/ PD-L2 was determined with ChIP-seq data. The association of PD-L1/PD-L2 expression with clinicopathological parameters, survival, immune infiltration and tumor mutation burden were investigated with data from The Cancer Genome Atlas. Potential targets and pathways of PD-L1 and PD-L2 were determined by protein enrichment, WebGestalt and gene ontology. Results: Comprehensive analysis revealed that PD-L1 and PD-L2 were significantly upregulated in most types of gastrointestinal cancers and their expressions were positively correlated. SP1 was a key transcription factor regulating the expression of PD-L1. Conclusion: Higher PD-L1 or PD-L2 expression was significantly associated with poor overall survival, higher tumor mutation burden and more immune and stromal cell populations. Finally, HIF-1, ERBB and mTOR signaling pathways were most significantly affected by PD-L1 and PD-L2 dysregulation. Altogether, this study provided comprehensive analysis of the dysregulation of PD-L1 and PD-L2, its underlying mechanism and downstream pathways, which add to the knowledge of manipulating PD-L1/PD-L2 for cancer immunotherapy.


2021 ◽  
Vol 49 (1) ◽  
pp. 030006052098153
Author(s):  
Qing Bi ◽  
Yang Liu ◽  
Tao Yuan ◽  
Huizhen Wang ◽  
Bin Li ◽  
...  

Objective The role of tumor-infiltrating lymphocytes (TILs) has not yet been characterized in sarcomas. The aim of this bioinformatics study was to explore the effect of TILs on sarcoma survival and genome alterations. Methods Whole-exome sequencing, transcriptome sequencing, and survival data of sarcoma were obtained from The Cancer Genome Atlas. Immune infiltration scores were calculated using the Tumor Immune Estimation Resource. Potential associations between abundance of infiltrating TILs and survival or genome alterations were examined. Results Levels of CD4+ T cell infiltration were associated with overall survival of patients with pan-sarcomas, and higher CD4+ T cell infiltration levels were associated with better survival. Somatic copy number alterations, rather than mutations, were found to correlate with CD4+ T cell infiltration levels. Conclusions This data mining study indicated that CD4+ T cell infiltration levels predicted from RNA sequencing could predict sarcoma prognosis, and higher levels of CD4+ T cells infiltration indicated a better chance of survival.


Cancers ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 662
Author(s):  
Mario Mischkulnig ◽  
Barbara Kiesel ◽  
Daniela Lötsch ◽  
Thomas Roetzer ◽  
Martin Borkovec ◽  
...  

Diffusely infiltrating gliomas are characterized by a variable clinical course, and thus novel prognostic biomarkers are needed. The heme biosynthesis cycle constitutes a fundamental metabolic pathway and might play a crucial role in glioma biology. The aim of this study was thus to investigate the role of the heme biosynthesis mRNA expression signature on prognosis in a large glioma patient cohort. Glioma patients with available sequencing data on heme biosynthesis expression were retrieved from The Cancer Genome Atlas (TCGA). In each patient, the heme biosynthesis mRNA expression signature was calculated and categorized into low, medium, and high expression subgroups. Differences in progression-free and overall survival between these subgroups were investigated including a multivariate analysis correcting for WHO grade, tumor subtype, and patient age and sex. In a total of 693 patients, progression-free and overall survival showed a strictly monotonical decrease with increasing mRNA expression signature subgroups. In detail, median overall survival was 134.2 months in the low, 79.9 months in the intermediate, and 16.5 months in the high mRNA expression signature subgroups, respectively. The impact of mRNA expression signature on progression-free and overall survival was independent of the other analyzed prognostic factors. Our data indicate that the heme biosynthesis mRNA expression signature might serve as an additional novel prognostic marker in patients with diffusely infiltrating gliomas to optimize postoperative management.


2021 ◽  
Author(s):  
Xiao-Cheng Wang ◽  
Ya Liu ◽  
Fei-Wu Long ◽  
Liang-Ren Liu ◽  
Chuan-Wen Fan

Background: The relationship between long noncoding RNAs (lncRNAs) and the mRNA stemness index (mRNAsi) in colorectal cancer (CRC) is still unclear. Materials & methods: The mRNAsi, mRNAsi-related lncRNAs and their clinical significance were analyzed by bioinformatic approaches in The Cancer Genome Atlas (TCGA)-COREAD dataset. Results: mRNAsi was negatively related to pathological features but positively related to overall survival and recurrence-free survival in CRC. A five mRNAsi-related lncRNAs prognostic signature was further developed and showed independent prognostic factors related to overall survival in CRC patients, due to the five mRNAsi-related lncRNAs involved in several pathways of the cancer stem cells and malignant cancer cell phenotypes. Conclusion: The present study highlights the potential roles of mRNAsi-related lncRNAs as alternative prognostic markers.


Sign in / Sign up

Export Citation Format

Share Document