scholarly journals Characterization of m6A Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration in Ovarian Cancer

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 ◽  
Author(s):  
Wancheng Zhao ◽  
Lili Yin

Abstract Background: Hypoxia-related genes have been reported to play important roles in a variety of cancers. However, their roles in ovarian cancer (OC) have remained unknown. The aim of our research was to explore the significance of hypoxia-related genes in OC patients.Methods: In this study, 15 hypoxia-related genes were screened from The Cancer Genome Atlas (TCGA) database to group the ovarian cancer patients using the consensus clustering method. Principal component analysis (PCA) was performed to calculate the hypoxia score for each patient to quantify the hypoxic status. Results: The OC patients from TCGA-OV dataset were divided into two distinct hypoxia statuses (cluster.A and cluster.B) based on the expression level of the 15 hypoxia-related genes. Most hypoxia-related genes were expressed more highly in the cluster.A group than in the cluster.B group. We also found that patients in the cluster.A group exhibited higher expression of immune checkpoint-related genes, epithelial-mesenchymal transition-related genes, and immune activation-related genes, as well as elevated immune infiltrates. PCA algorithm indicated that patients in the cluster.A group had higher hypoxia scores than that in in the cluster.B group.Conclusions: In summary, our research elucidated the vital role of hypoxia-related genes in immune infiltrates of OC. Our investigation of hypoxic status may be able to improve the efficacy of immunotherapy for OC.


2020 ◽  
Vol 13 (1) ◽  
Author(s):  
Feng-Juan Han ◽  
Jia Li ◽  
Ying Shen ◽  
Ying Guo ◽  
Yi-Chao Liu ◽  
...  

Abstract Objective Ovarian cancer (OC) has been regarded as the most malignant gynecological neoplasm and often confers grave outcomes owing to the frequent metastasis and high recurrence. A previous study has demonstrated that miR-1271-5p is implicated in OC progression, however, the possible mechanism of it remains unknown. The purpose of this investigation was to explore how miR-1271-5p regulates the progression of OC. Methods Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases were employed to analyze the differentially expressed miRNAs or genes as well as their corresponding prognostic values. miR-1271-5p expression in OC cells was examined by qRT-PCR. Cell counting kit 8 (CCK-8), colony formation, and transwell tests were conducted to evaluate the proliferation, migration and invasion potentials. Bioinformatics prediction and luciferase activity analysis were utilized to predict and verify the target gene of miR-1271-5p. Western blot assay was carried out to measure protein expression. Results miR-1271-5p was significantly decreased in OC and its down-regulation was associated with the grave outcome of OC patients. Upregulation of miR-1271-5p inhibited cell viability, but miR-1271-5p knockdown promoted the proliferation of OC cells. TIAM1 was a direct target gene of miR-1271-5p and expressed in OC tissues at higher level. High expression of TIAM1 induced the poorer prognosis of patients with OC. Further functional analyses showed that the suppressive role of miR-1271-5p on OC cell malignant behaviors was overturned by the upregulation of TIAM1. The protein levels of Cyclin D1, HES1, NOTCH and NUMB were remarkably changed due to the abnormal expression of miR-1271-5p and TIAM1. Conclusion To sum up, miR-1271-5p inhibits proliferation, invasion and migration of OC cells by directly repressing TIAM1 to inactivate the Notch signaling pathway, which provides an alternative therapeutic candidate for the advancement of OC treatment.


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.


2020 ◽  
Vol 8 (1) ◽  
pp. e000613
Author(s):  
Nicholas Bevins ◽  
Shulei Sun ◽  
Zied Gaieb ◽  
John A Thorson ◽  
Sarah S Murray

BackgroundTumor mutation burden (TMB) is a biomarker frequently reported by clinical laboratories, which is derived by quantifying of the number of single nucleotide or indel variants (mutations) identified by next-generation sequencing of tumors. TMB values can inform prognosis or predict the response of a patient’s tumor to immune checkpoint inhibitor therapy. Methods for the calculation of TMB are not standardized between laboratories, with significant variables being the gene content of the panels sequenced and the inclusion or exclusion of synonymous variants in the calculations. The impact of these methodological differences has not been investigated and the concordance of reported TMB values between laboratories is unknown.MethodsSequence variant lists from more than 9000 tumors of various types were downloaded from The Cancer Genome Atlas. Variant lists were filtered to include only appropriate variant types (ie, non-synonymous only or synonymous and non-synonymous variants) within the genes found in five commonly used targeted solid tumor gene panels as well as an in-house gene panel. Calculated TMB was paired with corresponding overall survival (OS) data of each patient.ResultsRegression analysis indicates high concordance of TMB as derived from the examined panels. TMB derived from panels was consistently and significantly lower than that derived from a whole exome. TMB, as derived from whole exome or the examined panels, showed a significant correlation with OS in the examined data.ConclusionsTMB derived from the examined gene panels was analytically equivalent between panels, but not between panels and whole-exome sequencing. Correlation between TMB and OS is significant if TMB method-specific cut-offs are used. These results suggest that TMB values, as derived from the gene panels examined, are analytically and prognostically equivalent.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Rong Tang ◽  
Xiaomeng Liu ◽  
Wei Wang ◽  
Jie Hua ◽  
Jin Xu ◽  
...  

Abstract Background High tumor mutation burden (TMB) has gradually become a sensitive biomarker for predicting the response to immunotherapy in many cancers, including lung, bladder and head and neck cancers. However, whether high TMB predicts the response to immunotherapy and prognosis in pancreatic ductal adenocarcinoma (PDAC) remained obscure. Hence, it is significant to investigate the role of genes related to TMB (TRGs) in PDAC. Methods The transcriptome and mutation data of PDAC was downloaded from The Cancer Genome Atlas-Pancreatic Adenocarcinoma (TCGA). Five independent external datasets of PDAC were chosen to validate parts of our results. qRT-PCR and immunohistochemical staining were also performed to promote the reliability of this study. Results The median overall survival (OS) was significantly increased in TMB_low group compared with the counterpart with higher TMB score after tumor purity adjusted (P = 0.03). 718 differentially expressed TRGs were identified and functionally enriched in some oncogenic pathways. 67 TRGs were associated with OS in PDAC. A prognostic model for the OS was constructed and showed a high predictive accuracy (AUC = 0.849). We also found TMB score was associated with multiple immune components and signatures in tumor microenvironment. In addition, we identified a PDAC subgroup featured with TMBlowMicrosatellite instabilityhigh (MSIhigh) was associated with prolonged OS and a key molecule, ANKRD55, potentially mediating the survival benefits. Conclusion This study analyzed the biological function, prognosis value, implications for mutation landscape and potential influence on immune microenvironment of TRGs in PDAC, which contributed to get aware of the role of TMB in PDAC. Future studies are expected to investigate how these TRGs regulate the initiation, development or repression of PDAC.


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.


2021 ◽  
Vol 5 (1) ◽  
Author(s):  
Mingdong Lu ◽  
Bin Zhao ◽  
Mengshan Liu ◽  
Le Wu ◽  
Yingying Li ◽  
...  

AbstractHistone methyltransferase SETD2 plays a critical role in maintaining genomic integrity and stability. Here, we investigated the characteristics of SETD2 somatic mutation in the cancer genome atlas pan-cancer cohort. Our data revealed that, compared with SETD2 nonmutant patients, SETD2 mutant patients had higher tumor mutation burden and microsatellite instability. In addition, the transcriptions of most genes related to immune activities were upregulated in patients with SETD2 mutant tumors. Further examination of cancer patients treated with immune checkpoint inhibitors suggested SETD2 mutation was associated with favorable clinical outcomes. These results have implication for the personalization of cancer immunotherapy.


2021 ◽  
Author(s):  
Wancheng Zhao ◽  
Lili Yin

Abstract Background: Hypoxia-related genes have been reported to play important roles in a variety of cancers. However, their roles in ovarian cancer (OC) have remained unknown. Our aim was to explore the significance of hypoxia-related genes in OC patients.Methods: In this study, 15 hypoxia-related genes were screened from The Cancer Genome Atlas (TCGA) database to group the ovarian cancer patients using the consensus clustering method. Principal component analysis (PCA) was used to calculate the hypoxia score to quantify the hypoxic status. Results: The OC patients were divided into two hypoxia statuses based on the expression level of the 15 hypoxia-related genes. Most hypoxia-related genes were more highly expressed in the cluster.A group. We also found that patients in the cluster.A group exhibited higher expression level of immune checkpoint-related genes, epithelial-mesenchymal transition-related genes, and immune activation-related genes, as well as elevated immune infiltrates. PCA algorithm indicated that the cluster.A group had higher hypoxia scores.Conclusions: In summary, our research elucidated the vital role of hypoxia-related genes in immune infiltrates of OC. Our investigation of hypoxic status may be able to improve the efficacy of immunotherapy for OC.


Sign in / Sign up

Export Citation Format

Share Document