Identification of m6A methyltransferase-related genes predicts prognosis and immune infiltrates in head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) accounts for 90% of head and neck malignant tumors. As the early symptoms of HNSCC are not obvious, and it is prone to recurrence and metastasis, making the overall survival (OS) rate of patients very low. Existing studies have shown m6A methylation plays a crucial role in various cancers, but it is rarely studied in HNSCC. This study aimed to explore the expression of m6A methylation-related genes in HNSCC and its correlation with prognosis, and to explore its relationship with immune infiltration. The gene expression data of HNSCC patient tumor samples (tumor =510) and adjacent normal tissue samples (normal =50) were extracted from The Cancer Genome Atlas (TCGA) database, and the expression characteristics of m6A regulatory factors were described. Kaplan-Meier survival analysis was used to analyze the relationship between m6A regulatory factors and OS and disease-specific survival (DSS). Least absolute shrinkage and selection operator (LASSO) regression was used to construct the m6A regulatory factor-HNSCC risk prediction model. In addition, the relationship between m6A methylation-related genes and tumor immune infiltration were discussed. The differential expression of 20 genes were identified by TCGA, and 18 genes (IGF2BP2, IGF2BP1, IGF2BP3, VIRMA, YTHDF1, YTHDF2, YTHDF3, ZC3H13, METTL14, ALKBH5, METTL3, RBMX, WTAP, YTHDC1, FTO, HNRNPC, HNRNPA2B1, and RBM15) were overexpressed in HNSCC. The survival rate of different gene expression levels was different. The high expression of YTHDC1 and YTHDC2 indicated better OS. Furthermore, for DSS, increased expression of YTHDC2 was also correlated with better clinical outcomes (P<0.05). At the same time, we drew a 3-gene risk score model in the TCGA-HNSCC cohort, and the survival curve showed compared with low-risk patients, high-risk patients had significantly worse OS (P<0.05). Gene enrichment analysis showed EPITHELIAL_MESENCHYMAL_TRANSITIO, MTORC1_SIGNALING, MYC_TARGETS_V1, MYC_TARGETS_V2, MYOGENESIS pathways, high TP53 mutations, and suppressive immunity were related to the high-risk group. The low-risk group was related to ALLOGRAFT_REJECTION, COMPLEMENT, IL6_JAK_STAT3_SIGNALING, INTERFERON_ALPHA_RESPONSE, INTERFERON_GAMMA_RESPONSE pathways, low TP53 mutations, and active immunity. The m6A methyltransferase-related genes can predict the prognosis of HNSCC and are related to immune infiltration.

Zhang Y ，Li L ，Ye Z ，Zhang L ，Yao N ，Gai L ... -  《-》

Genomic and transcriptomic alterations in m6A regulatory genes are associated with tumorigenesis and poor prognosis in head and neck squamous cell carcinoma.

Genetic alterations in N6-methyladenosine (m6A) regulatory genes are observed in many cancers. Recent studies have shown that newly identified m6A regulatory gene family (IGF2BPs; IGF2BP1, IGF2BP2, and IGF2BP3) were highly expressed in various types of cancer that stabilize and promote translation of multiple oncogenes, resulting in tumor development, survival and drug resistance. However, the oncogenic roles and prognostic values of IGF2BPs in head and neck squamous cell carcinoma (HNSCC) remain largely unknown. In this study, we examined the m6A regulatory genes alteration, their mRNAs expression and the prognostic values in HNSCC. We also analyzed the interaction network and functional enrichment of m6A regulators. Our results showed that m6A regulatory genes were altered in 41% (205/504) of HNSCC patients, of which IGF2BP2 was amplified in 20% (101/504) of HNSCC patents and positively correlated with its mRNA expression. Importantly, we have validated the expression of IGF2BP2 in HNSCC and normal tissue samples. Interestingly, we also found that the IGF2BP2 was frequently co-amplified with the most common oncogenes in HNSCC patients. In addition, this study found that other m6A regulatory genes such as METTL3, METTL14, WTAP, KIAA1429, ZC3H13, RBM15, ALKBH5, FTO, YTHDF1, YTHDF2, YTHDF3, YTHDC1, IGF2BP1, and IGF2BP3 were significantly upregulated in HNSCC samples. Moreover, patients with high expression of IGF2BP1, IGF2BP2, and IGF2BP3 had poor overall survival (OS) than those with low expression. Therefore, it is evident that IGF2BP family plays a key role in the oncogenesis of HNSCC and might serve as novel prognostic biomarkers and potential therapeutic targets in HNSCC.

Paramasivam A ，George R ，Priyadharsini JV 《American Journal of Cancer Research》

Comprehensive analysis of m6A related gene mutation characteristics and prognosis in colorectal cancer.

Colorectal cancer is considered as the second most common cancer worldwide. Studies have shown that m6A RNA methylation abnormalities play an important role in the pathogenesis of many human diseases, including cancer. The current study was designed to characterize the mutation of m6A related genes and explore their prognostic role in colorectal cancer. RNA-seq data and somatic mutation data of TCGA-COAD and TCGA-READ were downloaded from UCSC xena for comprehensive analysis. M6A related genes were selected from previous literatures, including "Writer" protein (METTL3, METTL5, METTL14, METTL16, ZC3H13, RBM15, WTAP, KIAA1429), "Reader" protein YTHDF1, YTHDF2, YTHDF3, YTHDC1, YTHDC2, HNRNPC, IGF2BP1, IGF2BP2, IGF2BP3), and "Eraser" protein (FTO, ALKBH5). Kaplan-Meier diagrams were used to explore the correlation between m6A-related genes and colorectal cancer prognosis. The correlations between m6A-related genes and clinical parameters and immune-related indicators were explored by Spearman correlation analysis. And finally, the expression patterns of five key genes (RBMX, FMR1, IGF2BP1, LRPPRC and YTHDC2) were detected by qPCR in CRC specimens. In CRC, the expressions of m6A-related genes were significantly different between CRC and normal control except METTL14, YTHDF2, YTHDF3. Some of CRC patients (178 in 536) have a m6A-related genes mutation. ZC3H13 has highest mutation frequency of all m6A-related genes. M6A-related genes mainly enrich in regulation of mRNA metabolic process pathway. Patients with high expressions of FMR1, LRPPRC, METTL14, RBMX, YTHDC2, YTHDF2, YTHDF3 have poor prognosis in CRC. There was a significant correlation between the FMR1, LRPPRC, RBMX, YTHDC2, IGF2BP1 expression and the clinical characteristics of CRC. In addition, these genes are significantly associated with immune-related indicators. According to the expression patterns of FMR1, LRPPRC, RBMX, YTHDC2, and IGF2BP1, patients with CRC were clustered into two groups, and their survival was significantly different. By evaluating the tumor microenvironment in two clusters using ssGSEA, expressions of immune checkpoints and GSVA enrichment analysis, we observed that the immune and stem cell index of two cluster were much different. The qPCR results showed that RBMX expression was markedly elevated in cancerous tissues than in the normal colonic tissues. Our study identified novel prognostic markers associated with immune of CRC cancer patients. Moreover, the potential mechanisms of prognostic markers in regulating the etiology of CRC cancer were investigated. These findings enrich our understanding of the relationships between m6a related genes and CRC, and may provide novel ideas in the therapy of CRC patients.

Jiang T ，Xing L ，Zhao L ，Ye Z ，Yu D ，Lin S ... -  《BMC Medical Genomics》

Analysis of Genetic Alteration Signatures and Prognostic Values of m6A Regulatory Genes in Head and Neck Squamous Cell Carcinoma.

Genetic alteration involving N6-methyladenosine (m6A) regulatory genes is a frequent characteristic of multiple tumors. Nevertheless, little is known regarding their genetic alteration signatures and prognostic values in head and neck squamous cell carcinoma (HNSCC). In this study, RNA sequence profiles and copy number variation (CNV) data of 506 HNSCC patients were downloaded from The Cancer Genome Atlas (TCGA) database. Correlation analysis involving alteration of m6A regulatory genes, clinicopathological characteristics, and patient survival was performed using R language. The results suggest that alteration of m6A regulatory genes was correlated with clinical staging. Patients with high expression of ALKBH5, FTO, METTL14, WTAP, YTHDC1, YTHDF1, and YTHDF2 had poor overall survival (OS) than those with low expression. Univariate and multivariate Cox regression analyses showed that ALKBH5 and YTHDC2 were independent risk factors for OS. However, patients with high YTHDC2 expression had better OS. Moreover, according to machine learning results, YTHDC2 was found to be the most important gene among the 10 m6A regulators. Additionally, high expression of YTHDC2 was correlated with activation of apoptosis and ubiquitin-mediated proteolysis. Here, we identified alterations to m6A regulatory genes in HNSCC for the first time and found that seven m6A regulators were associated with poor prognosis, especially ALKBH5, whereas YTHDC2 was associated with better prognosis. These m6A-related regulators could act as novel prognostic biomarkers for HNSCC. Our findings provide clues for understanding RNA epigenetic modifications in HNSCC.

Zhou X ，Han J ，Zhen X ，Liu Y ，Cui Z ，Yue Z ，Ding L ，Xu S ... -  《Frontiers in Oncology》

m6A Regulators Is Differently Expressed and Correlated With Immune Response of Esophageal Cancer.

N6 methyladenosine (m6A) RNA methylation regulators play an important role in the development of tumors. However, their function in esophageal cancer (EC) has not been fully elucidated. Here, we analyzed the gene expression data of 24 major m6A RNA methylation regulators from 775 patients with EC from TCGA dataset. The present study showed the aberrations of m6A regulators in genome were correlated to prognosis in human ECs. Meanwhile, 17 m6A regulators showed increased expression in EC samples, including YTHDC1, IGF2BP2, FTO, METTL14, YTHDF3, RBM15, WTAP, HNRNPA2B1, HNRNPC, ALKBH5, YTHDF2, METTL16, IGF2BP3, VIRMA, RBM15B, YTHDF1, KIAA1429, HAKAI, and ZC3H13. Among them, we found HNRNPC, YTHDC2, WTAP, VIRMA, IGF2BP3, and HNRNPA2B1 were significantly correlated to worse outcomes and advanced stage in EC. Furthermore, we showed levels of m6A regulators is correlated with the expression of Immuno-regulators (Immunoinhibitors, Immunostimulators, and MHC molecules) and immune infiltration levels in EC. Bioinformatics further confirm m6A regulators were involved in regulating RNA splicing, RNA stability, and cell proliferation. Our study showed m6A regulators are promising targets and biomarkers for cancer immunotherapy in EC.

Zhao H ，Xu Y ，Xie Y ，Zhang L ，Gao M ，Li S ，Wang F ... -  《Frontiers in Cell and Developmental Biology》