Prognostic values of m6A RNA methylation regulators in differentiated Thyroid Carcinoma

Comprehensive analysis of m6A RNA methylation regulators in esophageal carcinoma

Background

N6-methyladenosine (m6A) is the most pervasive modification of RNA methylation in eukaryotic cells. m6A modification plays a pivotal role in tumorigenesis and progression in many types of cancers. Until now, the role of m6A modification in esophageal carcinoma (ESCA) has remained obscure. The aim of the study was to construct and validate prognostic signatures based on m6A regulators for ESCA.

Methods

Transcriptomic data, somatic mutations and clinical information were obtained from The Cancer Genome Atlas (TCGA). Copy number variations were obtained from the UCSC (University of California, Santa Cruz) Xena database. We curated 21 m6A regulators and performed consensus clustering analysis to quantify the m6A modification pattern.

Results

Of the 184 patients, 23 (12.5%) were genetically altered in m6A regulators, with the highest frequency of mutations in ZC3H13 and LRPPRC. We constructed a m6A score system to investigate the prognosis of ESCA. The m6A score was closely related to immune cell infiltration in the tumor immune microenvironment. Patients with a high m6A score had an unfavorable prognosis. The combination of tumor mutation burden and m6A score would improve the prognostic value.

Conclusions

Our study established and validated a strong prognostic signature based on m6A regulators. This can be used to accurately predict the prognosis of ESCA.

The emerging roles of N6-methyladenosine RNA modifications in thyroid cancer

Thyroid cancer (TC) is the most predominant malignancy of the endocrine system, with steadily growing occurrence and morbidity worldwide. Although diagnostic and therapeutic methods have been rapidly developed in recent years, the underlying molecular mechanisms in the pathogenesis of TC remain enigmatic. The N6-methyladenosine(m6A) RNA modification is designed to impact RNA metabolism and further gene regulation. This process is intricately regulated by a variety of regulators, such as methylases and demethylases. Aberrant m6A regulators expression is related to the occurrence and development of TC and play an important role in drug resistance. This review comprehensively analyzes the effect of m6A methylation on TC progression and the potential clinical value of m6A regulators as prognostic markers and therapeutic targets in this disease.

YTHDF3as a prognostic predictive biomarker of thyroid cancer and its correlation with immune infiltration

PURPOSE

Thyroid cancer (TC) is one of the most common endocrine malignancies, and its morbidity continues to rise. N6-methyladenosine (m6A) RNA methylation, an epigenetic modification, is an important regulator of gene expression in TC. Therefore, it's worth finding the characteristics and predictive value of the m6A RNA methylation regulators in thyroid cancer (TC).

METHOD

RNA-seq data of TC was downloaded from the Cancer Genome Atlas (TCGA) database to screen out the differential expressed regulators. The absolute contraction selection operator (Lasso) Cox regression was used to construct the risk model of m6A methylation regulators. The predictive value of the risk scoring model was evaluated by Kaplan Meier (K-M) analysis and receiver operating characteristic (ROC) curves. The underlying mechanism of m6A methylation regulators in TC was predicted by gene set enrichment analysis (GSEA). Further validation was performed by using immunohistochemistry (IHC) and q-PCR. The correlation between risk-related gene and immune infiltration was evaluated by Tumour Immune Estimation Resource (TIMER).

RESULTS

IGF2BP2, YTHDF1 and YTHDF3 were screened out as strong independent prognostic factors of TC. Then a risk score model was established to further screen the predictors. Finally, according to the results of overall survival (OS) and clinical characteristics of TC, YTHDF3 was screened out as a potential predictor. Meanwhile, IHC and qPCR confirmed that YTHDF3 was expressed differential in TC. The expression of YTHDF3 was positively associated with the infiltration level of CD4+ T cells and macrophages. It was strongly correlated with a variety of immune markers in TC.

CONCLUSION

We confirmed that YTHDF3 can be used as a potential prognostic biomarker of TC. It not only plays a decisive role in the initiation and development of TC, but also provides a new perspective for understanding the modification of m6A RNA in TC.

Comprehensive analysis of m6A regulators and relationship with tumor microenvironment, immunotherapy strategies in colorectal adenocarcinoma

BACKGROUND

The N6-methyladenosine (m6A) RNA modification is the most prevalent and abundant type found in eukaryotic cells. It plays a crucial role in the initiation and progression of cancers. In this study, we aimed to comprehensively investigate the landscape of m6A regulators and their association with tumor microenvironment (TME), immunotherapeutic strategies in colon adenocarcinoma (COAD).

RESULTS

The differential expression, mutation, CNV frequency and prognostic value of 27 m6A regulators were systematically analyzed in COAD. Patients were classified into two clusters based on m6A regulators through consistent clustering analysis, with cluster A showing significant survival benefits. Most of the m6A regulators were negatively correlated with immune cells, except for WTAP, IGF2BP3, FTO, ALKBH5, which showed a positive correlation. We developed an m6A scoring system to calculate the m6Ascore for each patient. Patients with a high-m6Ascore had a better outcome, with the AUC of 0.775. An independent cohort of 416 COAD patients acquired from GSE38832 database was used to validate the prognosis prediction ability of m6Ascore. Moreover, the m6Ascore was negatively correlated with infiltration of anti-tumor immune cells. Additionally, patients with a high-m6Ascore responded better to anti-PD1 and anti-CTLA4 therapies, and those with MSI-H had a higher m6Ascore. Finally, we investigated the value of m6Ascore in predicting the response of patients to 15 commonly used drugs.

CONCLUSIONS

We comprehensively analyzed m6A regulators in COAD, including RNA expression, CNV changes, mutations and their correlation with TME. Our results showed that the m6A scoring system had significant predictive power for the prognosis of COAD patients, potentially leading to new personalized immunotherapy strategies.

Methyltransferase-like 3 induces the development of cervical cancer by enhancing insulin-like growth factor 2 mRNA-binding proteins 3-mediated apoptotic chromatin condensation inducer 1 mRNA stability

N6-methyladenosine (m6A) plays a critical role in the tumorigenesis of cervical cancer (CC). Here, we aimed to investigate the potential role of methyltransferase-like 3 (METTL3) in CC. Gene expression was determined via real-time quantitative polymerase chain reaction. Cellular functions were detected using colony formation, 5-ethynyl-2'-deoxyuridine (EdU), and Transwell assays. The interactions among METTL3, insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3), and apoptotic chromatin condensation inducer 1 (ACIN1) were confirmed using the MeRIP and RIP assays. An in vivo assay was performed to verify the role of METTL3 in CC development. METTL3 is overexpressed in CC, and therefore, its knockdown inhibits the proliferation and migration of CC cells. Silencing METTL3 inhibits tumor growth in vivo. Moreover, a positive association was observed between METTL3 and ACIN1. METTL3 interacts with IGF2BP3 to promote the mRNA stability of ACIN1, the overexpression of which induces the aggressiveness of CC cells. METTL3 promotes ACIN1 mRNA stability to accelerate CC progression, implying that METTL3 is a promising biomarker in CC.

Targeting IGF2BP2 Promotes Differentiation of Radioiodine Refractory Papillary Thyroid Cancer via Destabilizing RUNX2 mRNA

Simple Summary

Differentiation therapy is one of the most promising treatment approaches for radioiodine refractory papillary thyroid cancer (RR-PTC). In this study, we found that insulin-like growth factor 2 mRNA-binding protein 2 promoted dedifferentiation of PTC via integrating to 3′-untranslated regions of runt-related transcription factor 2, which bound to the promoter region of sodium/iodide symporter, downregulating its expression.

Abstract

N6-methyladenosine (m6A) regulators play an important role in multiple biological and pathological processes of radioiodine refractory papillary thyroid cancer (RR-PTC). However, the function of m6A regulators in differentiation of RR-PTC remains unclear. In this study, online data, clinical samples, and RR-PTC cell lines (K1 and TPC1) were used to identify the m6A regulators that contributed to the differentiation of RR-PTC. Insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) was found to be associated with thyroid-specific genes in online data analyses, and metastatic PTCs with high expression of IGF2BP2 were prone to be ¹³¹I-nonavid in clinical analyses. Furthermore, targeting IGF2BP2 increased ¹²⁵I uptake in RR-PTC cell lines and enhanced the sodium/iodide symporter (NIS) expression. Mechanistically, IGF2BP2 bound to the m6A modification site of runt-related transcription factor 2 (RUNX2) 3′-UTR and enhanced the RUNX2 mRNA stability. Moreover, RUNX2 could bind to the promoter region of NIS to block the differentiation of RR-PTC. Together, these results demonstrated that IGF2BP2 represents a diagnostic marker for RR-PTC, suggesting a novel differentiation therapeutic strategy of targeting IGF2BP2.

The Role of m6A RNA Methylation-Related lncRNAs in the Prognosis and Tumor Immune Microenvironment of Papillary Thyroid Carcinoma

Emerging evidence has indicated that N6-methylandenosine (m6A) RNA methylation plays a critical role in cancer development. However, the function of m6A RNA methylation-related long noncoding RNAs (m6A-lncRNAs) in papillary thyroid carcinoma (PTC) has never been reported. This study aimed to investigate the role of m6A-lncRNAs in the prognosis and tumor microenvironment (TME) of PTC. Three subgroups (clusters 1, 2, and 3) were identified by consensus clustering of 19 prognosis-related m6A-lncRNA regulators, of which cluster 1 is preferentially related to unfavorable prognosis, lower immune scores, and distinct immune infiltrate level. A risk-score model was established based on 8 prognosis-related m6A-lncRNAs. Patients with a high-risk score showed a worse prognosis, and the ROC indicated a reliable prediction performance for patients with PTC (AUC = 0.802). As expected, the immune scores, the infiltration levels of immune cells, and ESTIMATE scores in the low-risk subgroups were notably higher (p < 0.001) when compared with those in high-risk subgroups. Furthermore, GSEA analysis revealed that tumor associated pathways, hallmarks, and biological processes were remarkably enriched in the high-risk subgroup. Further analysis indicated that the risk score and age were independent prognostic factors for PTC. An integrated nomogram was constructed that accurately predicted the survival status (AUC = 0.963). Moreover, a lncRNA-miRNA-mRNA regulated network was established based on seven prognosis-related m6A-lncRNAs. In addition, 30 clinical samples and different PTC cells were validated. This is the first study to reveal that m6A-lncRNAs plays a vital role in the prognosis and TME of PTC. To a certain degree, m6A-lncRNAs can be considered as new, promising prognostic biomarkers and treatment targets.

Analysis of the prognostic significance and potential mechanisms of lncRNAs associated with m6A methylation in papillary thyroid carcinoma

BACKGROUND

m6A methylation-related long non-coding RNAs (lncRNAs) play a significant role in the progression of various tumors and can be used as prognostic markers. However, whether m6A-related lncRNAs also play the same function as prognostic markers in papillary thyroid carcinoma (PTC) remains unclear.

METHODS

Consensus cluster analysis was performed to divide PTC samples obtained from The Cancer Genome Atlas database into two clusters according to the expression of m6A-related lncRNAs. Then, the least absolute shrinkage and selection operator (LASSO) regression analysis was performed to create and verify a prognostic model. Furthermore, the relationship among risk scores, clusters, programmed death-ligand 1 (PD-L1), tumor microenvironment (TME), clinicopathological characteristics, immune infiltration, immune checkpoint, and tumor mutation burden (TMB) was analyzed. In addition, a nomogram was created, and subsequently, the drug sensitivity of lncRNAs in the prognostic model was analyzed. Finally, the relationship between these lncRNAs and prognosis in pan-cancer was investigated.

RESULTS

The prognosis, RAS, BRAF, M, and TME were found to be different in two clusters. The prognostic model included three lncRNAs: PSMG3-AS1, BHLHE40-AS1, and AC016747.3. The risk score was associated with clusters, PD-L1, tumor microenvironment, clinicopathological characteristics, immune cell infiltration, immune checkpoint, and TMB, and thus, risk score was confirmed as useful prognostic indicator. Differentially expressed lncRNAs are involved in many malignancies and can be identified as cancer prognostic makers.

CONCLUSION

According to our research, we can regard m6A-related lncRNAs involved in the procession of PTC as a biomarker of progression-free survival for PTC patients, and pan-cancer.

m6A regulator-mediated RNA methylation modification patterns and immune microenvironment infiltration characterization in severe asthma

N6‐methyladenosine (m6A) modification is one of the most prevalent RNA modification forms of eukaryotic mRNA and is an important post‐transcriptional mechanism for regulating genes. However, the role of m6A modification in the regulation of severe asthma has never been reported. Thus, we aimed to investigate the m6A regulator‐mediated RNA methylation modification patterns and immune microenvironment infiltration characterization in severe asthma. In this study, 87 healthy controls and 344 severe asthma cases from the U‐BIOPRED (Unbiased Biomarkers for the Prediction of Respiratory Disease Outcomes) programme were used to systematically evaluate the m6A modification patterns mediated by 27 m6A regulators and to investigate the effects of m6A modification on immune microenvironment characteristics. We found that 16 m6A regulators were abnormal and identified two key m6A regulators (YTHDF3 and YTHDC1) and three m6A modification patterns. The study of infiltration characteristics of immune microenvironment found that pattern 2 had more infiltrating immune cells and more active immune response. Besides, it was found that the eosinophils which are very important for severe asthma were affected by YTHDF3 and EIF3B. We also verified key m6A regulators with merip‐seq and found that they were mainly distributed in exons and enriched in 3′UTR. In conclusion, our findings suggested that m6A modification plays a key role in severe asthma, and may be able to guide the future strategy of immunotherapy.

CircHACE1 functions as a competitive endogenous RNA to curb differentiated thyroid cancer progression by upregulating Tfcp2L1 through adsorbing miR-346

Circular RNAs (circRNAs) are correlated with the occurrence and progression of differentiated thyroid cancer (THCA). However, the regulatory mechanism of circRNAs in differentiated THCA is unclear. In the present study, we analyzed the circRNA microarray dataset (GSE93522) of thyroid tumors and discovered that circRNA HACE1 (circHACE1) was downregulated in differentiated THCA. We detected circHACE1 expression by quantitative real-time polymerase chain reaction (qRT-PCR). Gain-of-function experiments were performed to analyze the biological function of circHACE1 in differentiated THCA cells in vitro. The regulatory mechanism of circHACE1 in differentiated THCA was explored through bioinformatics analysis, dual-luciferase reporter, RIP (RNA immunoprecipitation), and/or RNA pull-down assays. The biological function of circHACE1 in THCA was confirmed by xenograft assay. We verified that circHACE1 was downregulated in differentiated THCA. Also, differentiated THCA patients with low circHACE1 expression were associated with TNM grade, lymphoid node metastasis, tumor size, and poor prognosis. CircHACE1 overexpression decreased xenograft tumor growth in vivo and induced cell cycle arrest, apoptosis, impeded proliferation, migration, and invasion in differentiated THCA cells in vitro. CircHACE1 could function as a microRNA (miR)-346 sponge and regulated Tfcp2L1 (transcription factor CP2 like 1) expression. MiR-346 overexpression offset circHACE1 elevation-mediated effects on malignant behaviors of differentiated THCA cells. Furthermore, Tfcp2L1 silencing counteracted the suppressive impact of miR-346 inhibitor on the malignancy of differentiated THCA cells. In conclusion, circHACE1 adsorbed miR-346 and elevated Tfcp2L1 expression, thus curbing cell malignancy in differentiated THCA, manifesting that circHACE1 might be a target for differentiated THCA treatment.

Identification of Critical m6A RNA Methylation Regulators with Prognostic Value in Lower-Grade Glioma

Increasing evidences have revealed that N⁶-methyladenosine (m⁶A) RNA methylation regulators participate in the tumorigenesis and development of multiple tumors. So far, there has been little comprehension about the effects of m⁶A RNA methylation regulators on lower-grade gliomas (LGG). Here, we systematically investigated the expression profiles and prognostic significance of 36 m⁶A RNA methylation regulators in LGG patients from the TCGA and CGGA databases. Most of the m⁶A RNA methylation regulators are differentially expressed in LGG tissues as compared with normal brain tissues and glioblastoma (GBM) tissues. The consensus clustering for these m⁶A RNA methylation regulators identified three clusters. Patients in cluster 3 exhibited worse prognosis. In addition, we constructed an m⁶A-related prognostic signature, which exhibited excellent performance in prognostic stratification of LGG patients according to the results of the Kaplan-Meier curves, ROC curves, and univariate and multivariate Cox regression analyses. In addition, a significant correlation was observed between the m⁶A-related prognostic signature and the immune landscape of the LGG microenvironment. The high-risk group exhibited higher immune scores, stromal scores, and ESTIMATE scores but lower tumor purity and lower abundance of activated NK cells. Moreover, the expression level of immune checkpoints was positively correlated with the risk score. To conclude, the current research systematically demonstrated the prognostic roles of m⁶A RNA methylation regulators in LGG.

The Biological Function, Mechanism, and Clinical Significance of m6A RNA Modifications in Head and Neck Carcinoma: A Systematic Review

Head and neck squamous cell carcinoma (HNSCC) is one of the most common cancers, yet the molecular mechanisms underlying its onset and development have not yet been fully elucidated. Indeed, an in-depth understanding of the potential molecular mechanisms underlying HNSCC oncogenesis may aid the development of better treatment strategies. Recent epigenetic studies have revealed that the m6A RNA modification plays important roles in HNSCC. In this review, we summarize the role of m6A modification in various types of HNSCC, including thyroid, nasopharyngeal, hypopharyngeal squamous cell, and oral carcinoma. In addition, we discuss the regulatory roles of m6A in immune cells within the tumor microenvironment, as well as the potential molecular mechanisms. Finally, we review the development of potential targets for treating cancer based on the regulatory functions of m6A, with an aim to improving targeted therapies for HNSCC. Together, this review highlights the important roles that m6A modification plays in RNA synthesis, transport, and translation, and demonstrates that the regulation of m6A-related proteins can indirectly affect mRNA and ncRNA function, thus providing a novel strategy for reengineering intrinsic cell activity and developing simpler interventions to treat HNSCC.

Recent advances in functional annotation and prediction of the epitranscriptome

RNA modifications, in particular N⁶-methyladenosine (m⁶A), participate in every stages of RNA metabolism and play diverse roles in essential biological processes and disease pathogenesis. Thanks to the advances in sequencing technology, tens of thousands of RNA modification sites can be identified in a typical high-throughput experiment; however, it remains a major challenge to decipher the functional relevance of these sites, such as, affecting alternative splicing, regulation circuit in essential biological processes or association to diseases. As the focus of RNA epigenetics gradually shifts from site discovery to functional studies, we review here recent progress in functional annotation and prediction of RNA modification sites from a bioinformatics perspective. The review covers naïve annotation with associated biological events, e.g., single nucleotide polymorphism (SNP), RNA binding protein (RBP) and alternative splicing, prediction of key sites and their regulatory functions, inference of disease association, and mining the diagnosis and prognosis value of RNA modification regulators. We further discussed the limitations of existing approaches and some future perspectives.

Identification of the EMT-Related Genes Signature for Predicting Occurrence and Progression in Thyroid Cancer

BACKGROUND

The detection rate of thyroid cancer (TC) has been continuously improved due to the development of detection technology. Epithelial-mesenchymal transition (EMT) is thought to be closely related to the malignant progression of tumors. However, the relationship between EMT-related genes (ERGs) characteristics and the diagnosis and prognosis of TC patients has not been studied.

METHODS

Four datasets from Gene Expression Omnibus (GEO) were used to perform transcriptomic profile analysis. The overlapping differentially expressed ERGs (DEERGs) were analyzed using the R package "limma". Then, the hub genes, which had a higher degree, were identified by the protein-protein interaction (PPI) network. Gene expression analysis between the TC and normal data, the disease-free survival (DFS) analysis of TC patients from The Cancer Genome Atlas Thyroid Cancer (TCGA-THCA) cohort, function analysis, and immunohistochemistry (IHC) were performed to verify the importance of the hub genes. Finally, a prognostic risk scoring was constructed to predict DFS in patients with the selected genes.

RESULTS

A total of 43 DEERGs were identified and 10 DEERGs were considered hub ERGs, which had a high degree of connectivity in the PPI network. Then, the differential expressions of FN1, ITGA2, and KIT between TC and normal tissues were verified in the TCGA-THCA cohort and their protein expressions were also verified by IHC. DFS analysis indicated upregulations of FN1 expression (P<0.01) and ITGA2 expression (P<0.01) and downregulation of KIT expression (P=0.01) increased risks of decreased DFS for TCGA-THCA patients. Besides, by building a prognostic risk scoring model, we found that the DFS of TCGA-THCA patients was significantly worse in high-risk groups.

CONCLUSION

In summary, these hub ERGs were potential biomarkers for diagnosis and prognosis of TC, which can provide a basis for further exploring the efficacy of EMT in patients with TC.

M6A Demethylase FTO Plays a Tumor Suppressor Role in Thyroid Cancer

The purpose of this study was to investigate the expression and clinical significance of N6-methyladenosine demethylase FTO in thyroid cancer. Bioinformatic analysis showed that FTO expression was downregulated in thyroid cancer tissues and correlated with lymph node metastasis in thyroid cancer patients. We conducted experimental verification by collecting Asian samples. The results of quantitative reverse transcription-PCR showed that the mRNA expression of FTO in the blood of 30 thyroid cancer patients was lower than that of the control population. At the same time, we found that FTO expression was negative in tissues of 16/56 (28.57%) thyroid cancer cases and 4/40 (10.00%) nontumor thyroid cases through the immunohistochemical method, indicating a lower FTO expression in thyroid cancer tissues than nontumor thyroid tissues (p < 0.05). In addition, the protein expression of FTO was significantly related to the tumor grade and lymph node metastasis in thyroid cancer patients (p < 0.05), but not to other clinicopathological features. Multivariate logistic regression analysis showed that FTO expression was an independent risk factor for tumor grade. Survival analysis showed no significant difference in the disease-free survival time of thyroid cancer patients between high expression and low expression groups of FTO. Furthermore, bioinformatic analysis found that promoter DNA methylation and copy number variation might cause downregulated FTO and then affect TP53 pathways in thyroid cancer. We found that FTO expression was downregulated in thyroid cancer tissues and related to the progression of thyroid cancer, suggesting a tumor suppressor role of FTO in thyroid cancer.