Dynamic transcriptomic m6A decoration: writers, erasers, readers and functions in RNA metabolism

The role of YTH domain containing 2 in epigenetic modification and immune infiltration of pan-cancer

YTH domain containing 2 (YTHDC2) is the largest N6‐Methyladenosine (m⁶A) binding protein of the YTH protein family and the only member containing ATP‐dependent RNA helicase activity. For further analysing its biological role in epigenetic modification, we comprehensively explored YTHDC2 from gene expression, genetic alteration, protein‐protein interaction (PPI) network, immune infiltration, diagnostic value and prognostic value in pan‐cancer, using a series of databases and bioinformatic tools. We found that YTHDC2 with Missense mutation could cause a different prognosis in uterine corpus endometrial carcinoma (UCEC), and its different methylation level could lead to a totally various prognosis in adrenocortical carcinoma (ACC), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), lung squamous cell carcinoma (LUSC) and UCEC. The main molecular mechanisms of YTHDC2 focused on catalytic activity, helicase activity, snRNA binding, spliceosome and mRNA surveillance. Additionally, YTHDC2 was notably correlated with tumour immune infiltration. Moreover, YTHDC2 had a high diagnostic value for seven cancer types and a prognostic value for brain lower grade glioma (LGG), rectum adenocarcinoma (READ) and skin cutaneous melanoma (SKCM). Collectively, YTHDC2 plays a significant role in epigenetic modification and immune infiltration and maybe a potential biomarker for diagnosis and prognosis in certain cancers.

Integrated analysis of m6A mRNA methylation in rats with monocrotaline-induced pulmonary arterial hypertension

Background: N6-methyladenosine (m⁶A) modification is one of the most common chemical modifications of eukaryotic mRNAs, which play an important role in tumors and cardiovascular disease through regulating mRNA stability, splicing and translation. However, the changes of m⁶A mRNA and m⁶A-related enzymes in pulmonary arterial hypertension (PAH) remain largely unexplored.

Methods: MeRIP-seq was used to identify m⁶A methylation in lung tissues from control and MCT-PAH rats. Western blot and immunofluorescence were used to evaluate expression of m⁶A-related enzymes.

Results: Compared with control group, m⁶A methylation was mainly increased in lung tissues from MCT-PAH rats. The up-methylated coding genes in MCT-PAH rats were primarily enriched in processes associated with inflammation, glycolysis, ECM-receptor interaction and PDGF signal pathway, while genes with down-methylation were enriched in processes associated with TGF-β family receptor members. The expression of FTO and ALKBH5 downregulated, METTL3 and YTHDF1 increased and other methylation modification-related proteins was not significantly changed in MCT-PAH rats lung tissues. Immunofluorescence indicated that expression of FTO decreased and YTHDF1 increased in small pulmonary arteries of MCT-PAH rats.

Conclusion: m⁶A levels and the expression of methylation-related enzymes were altered in PAH rats, in which FTO and YTHDF1 may play a crucial role in m⁶A modification.

N6-Methyladenosine RNA modification in cerebrospinal fluid as a novel potential diagnostic biomarker for progressive multiple sclerosis

Background

Progressive multiple sclerosis (PMS) is an uncommon and severe subtype of MS that worsens gradually and leads to irreversible disabilities in young adults. Currently, there are no applicable or reliable biomarkers to distinguish PMS from relapsing–remitting multiple sclerosis (RRMS). Previous studies have demonstrated that dysfunction of N6-methyladenosine (m6A) RNA modification is relevant to many neurological disorders. Thus, the aim of this study was to explore the diagnostic biomarkers for PMS based on m6A regulatory genes in the cerebrospinal fluid (CSF).

Methods

Gene expression matrices were downloaded from the ArrayExpress database. Then, we identified differentially expressed m6A regulatory genes between MS and non-MS patients. MS clusters were identified by consensus clustering analysis. Next, we analyzed the correlation between clusters and clinical characteristics. The random forest (RF) algorithm was applied to select key m6A-related genes. The support vector machine (SVM) was then used to construct a diagnostic gene signature. Receiver operating characteristic (ROC) curves were plotted to evaluate the accuracy of the diagnostic model. In addition, CSF samples from MS and non-MS patients were collected and used for external validation, as evaluated by an m6A RNA Methylation Quantification Kit and by real-time quantitative polymerase chain reaction.

Results

The 13 central m6A RNA methylation regulators were all upregulated in MS patients when compared with non-MS patients. Consensus clustering analysis identified two clusters, both of which were significantly associated with MS subtypes. Next, we divided 61 MS patients into a training set (n = 41) and a test set (n = 20). The RF algorithm identified eight feature genes, and the SVM method was successfully applied to construct a diagnostic model. ROC curves revealed good performance. Finally, the analysis of 11 CSF samples demonstrated that RRMS samples exhibited significantly higher levels of m6A RNA methylation and higher gene expression levels of m6A-related genes than PMS samples.

Conclusions

The dynamic modification of m6A RNA methylation is involved in the progression of MS and could potentially represent a novel CSF biomarker for diagnosing MS and distinguishing PMS from RRMS in the early stages of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02981-5.

METTL3-dependent N6-methyladenosine RNA modification mediates the atherogenic inflammatory cascades in vascular endothelium

Atherosclerosis is characterized by the plaque formation that restricts intraarterial blood flow. The disturbed blood flow with the associated oscillatory stress (OS) at the arterial curvatures and branch points can trigger endothelial activation and is one of the risk factors of atherosclerosis. Many studies reported the mechanotransduction related to OS and atherogenesis; however, the transcriptional and posttranscriptional regulatory mechanisms of atherosclerosis remain unclear. Herein, we investigated the role of N⁶-methyladenosine (m⁶A) RNA methylation in mechanotransduction in endothelial cells (ECs) because of its important role in epitranscriptome regulation. We have identified m⁶A methyltransferase METTL3 as a responsive hub to hemodynamic forces and atherogenic stimuli in ECs. OS led to an up-regulation of METTL3 expression, accompanied by m⁶A RNA hypermethylation, increased NF-κB p65 Ser⁵³⁶ phosphorylation, and enhanced monocyte adhesion. Knockdown of METTL3 abrogated this OS-induced m⁶A RNA hypermethylation and other manifestations, while METTL3 overexpression led to changes resembling the OS effects. RNA-sequencing and m⁶A-enhanced cross-linking and immunoprecipitation (eCLIP) experiments revealed NLRP1 and KLF4 as two hemodynamics-related downstream targets of METTL3-mediated hypermethylation. The METTL3-mediated RNA hypermethylation up-regulated NLRP1 transcript and down-regulated KLF4 transcript through YTHDF1 and YTHDF2 m⁶A reader proteins, respectively. In the in vivo atherosclerosis model, partial ligation of the carotid artery led to plaque formation and up-regulation of METTL3 and NLRP1, with down-regulation of KLF4; knockdown of METTL3 via repetitive shRNA administration prevented the atherogenic process, NLRP3 up-regulation, and KLF4 down-regulation. Collectively, we have demonstrated that METTL3 serves a central role in the atherogenesis induced by OS and disturbed blood flow.

A comprehensive review of m6A/m6Am RNA methyltransferase structures

Gene expression is regulated at many levels including co- or post-transcriptionally, where chemical modifications are added to RNA on riboses and bases. Expression control via RNA modifications has been termed 'epitranscriptomics' to keep with the related 'epigenomics' for DNA modification. One such RNA modification is the N6-methylation found on adenosine (m6A) and 2'-O-methyladenosine (m6Am) in most types of RNA. The N6-methylation can affect the fold, stability, degradation and cellular interaction(s) of the modified RNA, implicating it in processes such as splicing, translation, export and decay. The multiple roles played by this modification explains why m6A misregulation is connected to multiple human cancers. The m6A/m6Am writer enzymes are RNA methyltransferases (MTases). Structures are available for functionally characterized m6A RNA MTases from human (m6A mRNA, m6A snRNA, m6A rRNA and m6Am mRNA MTases), zebrafish (m6Am mRNA MTase) and bacteria (m6A rRNA MTase). For each of these MTases, we describe their overall domain organization, the active site architecture and the substrate binding. We identify areas that remain to be investigated, propose yet unexplored routes for structural characterization of MTase:substrate complexes, and highlight common structural elements that should be described for future m6A/m6Am RNA MTase structures.

The Genetic Changes of Hepatoblastoma

Hepatoblastoma is the most common malignant liver cancer in childhood. The etiology of hepatoblastoma remains obscure. Hepatoblastoma is closely related to genetic syndromes, hinting that hepatoblastoma is a genetic predisposition disease. However, no precise exposures or genetic events are reported to hepatoblastoma occurrence. During the past decade, significant advances have been made in the understanding of etiology leading to hepatoblastoma, and several important genetic events that appear to be important for the development and progression of this tumor have been identified. Advances in our understanding of the genetic changes that underlie hepatoblastoma may translate into better patient outcomes. Single nucleotide polymorphisms (SNPs) have been generally applied in the research of etiology's exploration, disease treatment, and prognosis assessment. Here, we reviewed and discussed the molecular epidemiology, especially SNPs progresses in hepatoblastoma, to provide references for future studies and promote the study of hepatoblastoma's etiology.

The N6-Methyladenosine-Modified Pseudogene HSPA7 Correlates With the Tumor Microenvironment and Predicts the Response to Immune Checkpoint Therapy in Glioblastoma

Background

Glioblastoma (GBM), one of the most aggressive tumors of the brain, has no effective or sufficient therapies. Identifying robust biomarkers for the response to immune checkpoint blockade (ICB) therapy, a promising treatment option for GBM patients, is urgently needed.

Methods

We comprehensively evaluated lncRNA m⁶A modification patterns in m⁶A-sequencing (m⁶A-seq) data for GBM tissues and systematically investigated the immune and stromal regulators of these m⁶A-regulated lncRNAs. We used the single-sample gene-set enrichment analysis (ssGSEA) algorithm to investigate the difference in enriched tumor microenvironment (TME) infiltrating cells and the functional annotation of HSPA7 in individual GBM samples. Further, we validated that HSPA7 promoted the recruitment of macrophages into GBM TME in vitro, as well as in our GBM tissue section. We also explored its impact on the efficacy of ICB therapy using the patient-derived glioblastoma organoid (GBO) model.

Results

Here, we depicted the first transcriptome-wide m⁶A methylation profile of lncRNAs in GBM, revealing highly distinct lncRNA m⁶A modification patterns compared to those in normal brain tissues. We identified the m⁶A-modified pseudogene HSPA7 as a novel prognostic risk factor in GBM patients, with crucial roles in immunophenotype determination, stromal activation, and carcinogenic pathway activation. We confirmed that HSPA7 promoted macrophage infiltration and SPP1 expression via upregulating the YAP1 and LOX expression of glioblastoma stem cells (GSCs) in vitro and in our clinical GBM tumor samples. We also confirmed that knockdown of HSPA7 might increase the efficiency of anti-PD1 therapy utilizing the GBO model, highlighting its potential as a novel target for immunotherapy.

Conclusions

Our results indicated that HSPA7 could be a novel immunotherapy target for GBM patients.

The functions and prognostic values of m6A RNA methylation regulators in thyroid carcinoma

BACKGROUND

N6-Methyladenosine (m6A) is the most common RNA modification and regulates RNA splicing, translation, translocation, and stability. Aberrant expression of m6A has been reported in various types of human cancers. m6A RNA modification is dynamically and reversibly mediated by different regulators, including methyltransferase, demethylases, and m6A binding proteins. However, the role of m6A RNA methylation regulators in thyroid cancer remains unknown. The aim of this study is to investigate the effect of the 13 main m6A RNA modification regulators in thyroid carcinoma.

METHODS

We obtained clinical data and RNA sequencing data of 13 m6A RNA methylation regulators from The Cancer Genome Atlas (TCGA) THCA database. We performed consensus clustering to identify the clinical relevance of m6A RNA methylation regulators in thyroid carcinoma. Then we used LASSO Cox regression analysis to generate a prognostic signature based on m6A RNA modification regulator expression. Kyoto Encyclopedia of Genes and Genomes, Gene Ontology and Gene Set Enrichment Analyses were performed to explore differential cellular processes and signaling pathways between the two groups based on risk signature.

RESULTS

We found that most of the m6A RNA modification regulators are down-regulated in 450 patients with thyroid carcinoma. We derived a three m6A RNA modification regulator genes-based risk signature (FTO, RBM15 and KIAA1429), that is an independent prognostic biomarker in patients with thyroid carcinoma. Moreover, we found that this risk signature could better predict outcome in male than female. Functional research in vitro demonstrated that the m6A RNA methylation regulators involved in the model acted significant role in the proliferation and migration of thyroid cancer cells.

CONCLUSIONS

Our study revealed the influence of m6A RNA methylation regulators on thyroid carcinoma through biological experiments and three-gene prognostic model.

An introduction to the themed issue on RNA biology in China

The year 2019 marked the fortieth anniversary of the Chinese Society of Biochemistry and Molecular Biology (CSBMB), whose mission is to promote biomolecular research and education in China. The last 40 years have witnessed tremendous growth and achievements in biomolecular research by Chinese scientists and Essays in Biochemistry is delighted to publish this themed issue that focuses on exciting areas within RNA biology, with each review contributed by key experts from China.

LCAT3, a novel m6A-regulated long non-coding RNA, plays an oncogenic role in lung cancer via binding with FUBP1 to activate c-MYC

BACKGROUND

Long non-coding RNAs (lncRNAs) are important epigenetic regulators, which play critical roles in diverse physiological and pathological processes. However, the regulatory mechanism of lncRNAs in lung carcinogenesis remains elusive. Here, we characterized a novel oncogenic lncRNA, designated as Lung Cancer Associated Transcript 3 (LCAT3).

METHODS

We predicted and validated LCAT3 by analyzing RNA-sequencing (RNA-seq) data of lung cancer tissues from TCGA. Methylated RNA immunoprecipitation was performed to assess m6A modification on LCAT3. The LCAT3-FUBP1-MYC axis was assessed by dual-luciferase reporter, RNA immunoprecipitation and Chromatin immunoprecipitation assays. Signaling pathways altered by LCAT3 knockdown were identified using RNA-seq. Furthermore, the mechanism of LCAT3 was investigated using loss-of-function and gain-of-function assays in vivo and in vitro.

RESULTS

LCAT3 was found to be up-regulated in lung adenocarcinomas (LUAD), and its over-expression was associated with the poor prognosis of LUAD patients. LCAT3 upregulation is attributable to N6-methyladenosine (m6A) modification mediated by methyltransferase like 3 (METTL3), leading to LCAT3 stabilization. Biologically, loss-of-function assays revealed that LCAT3 knockdown significantly suppressed lung cancer cell proliferation, migration and invasion in vitro, and inhibited tumor growth and metastasis in vivo. LCAT3 knockdown induced cell cycle arrest at the G1 phase. Mechanistically, LCAT3 recruited Far Upstream Element Binding Protein 1 (FUBP1) to the MYC far-upstream element (FUSE) sequence, thereby activating MYC transcription to promote proliferation, survival, invasion and metastasis of lung cancer cells.

CONCLUSIONS

Taken together, we identified and characterized LCAT3 as a novel oncogenic lncRNA in the lung, and validated the LCAT3-FUBP1-MYC axis as a potential therapeutic target for LUAD.

N6-Methylandenosine-Related lncRNA Signature Is a Novel Biomarkers of Prognosis and Immune Response in Colon Adenocarcinoma Patients

Background

Colon adenocarcinoma (COAD) is the most common type of colon cancer. To date, however, the prognostic values of m6A RNA methylation-related long non-coding RNAs (lncRNAs) in COAD are largely unknown.

Materials and Methods

The m6A-related lncRNAs were identified from The Cancer Genome Atlas (TCGA) data set. Univariate and multivariate Cox regression analyses were performed to explore the prognostic m6A-related lncRNAs. Consistent clustering analysis was performed to classify the COAD patients into different subgroups based on the expression of m6A-related lncRNAs. The potential biological functions as well as differences in the stemness index and tumor immune microenvironment between different subgroups were analyzed. The prognostic m6A-related lncRNAs were used to establish an m6A-related lncRNA risk model to predict prognosis and survival status.

Results

We identified 31 m6A-associated lncRNAs with prognostic values from the TCGA data set. Based on the expression of prognostic m6A-associated lncRNAs, TCGA-COAD patients were classified into three clusters using consistent clustering analysis. There was a low correlation of tumor stemness between the three clusters but a significant correlation with the tumor immune microenvironment as well as the tumor mutational load. Thirty-one prognostic-related m6A-associated lncRNAs were used to construct a risk model, which was further determined by survival analysis, receiver operating characteristic (ROC) curve, and univariate and multifactor Cox analysis. The m6A-related risk model demonstrates good performance in predicting prognosis and survival status. The model-based high-risk group exhibited poorer overall survival (OS) compared with the low-risk group.

Conclusion

In this study, we construct a risk model that consists of 31 m6A-related lncRNAs with independent prognostic values in COAD. Our study shows the critical roles of these 31 m6A-related lncRNAs in the tumor immune microenvironment, indicating the prospect of informing prognostic stratification and the development of immunotherapeutic strategies for COAD patients.

YTHDF2 facilitates UBXN1 mRNA decay by recognizing METTL3-mediated m6A modification to activate NF-κB and promote the malignant progression of glioma

BACKGROUND

The prognosis for diffuse gliomas is very poor and the mechanism underlying their malignant progression remains unclear. Here, we aimed to elucidate the role and mechanism of the RNA N6,2'-O-dimethyladenosine (m6A) reader, YTH N6-methyladenosine RNA binding protein 2 (YTHDF2), in regulating the malignant progression of gliomas.

METHODS

YTHDF2 mRNA levels and functions were assessed using several independent datasets. Western blotting, quantitative polymerase chain reaction, and immunohistochemistry were used to evaluate the expression levels of YTHDF2 and other molecules in human and mouse tumor tissues and cells. Knockdown and overexpression were used to evaluate the effects of YTHDF2, methyltransferase-like 3 (METTL3), and UBX domain protein 1 (UBXN1) on glioma malignancy in cell and orthotopic xenograft models. RNA immunoprecipitation (RIP), methylated RIP, and RNA stability experiments were performed to study the mechanisms underlying the oncogenic role of YTHDF2.

RESULTS

YTHDF2 expression was positively associated with a higher malignant grade and molecular subtype of glioma and poorer prognosis. YTHDF2 promoted the malignant progression of gliomas in both in vitro and in vivo models. Mechanistically, YTHDF2 accelerated UBXN1 mRNA degradation via METTL3-mediated m6A, which, in turn, promoted NF-κB activation. We further revealed that UBXN1 overexpression attenuated the oncogenic effect of YTHDF2 overexpression and was associated with better survival in patients with elevated YTHDF2 expression.

CONCLUSIONS

Our findings confirmed that YTHDF2 promotes the malignant progression of gliomas and revealed important insight into the upstream regulatory mechanism of NF-κB activation via UBXN1 with a primary focus on m6A modification.

Coordination of m6A mRNA methylation and gene transcriptome in rice response to cadmium stress

N⁶-methyladenosine (m⁶A) is the most prevalent internal modification present in the mRNAs of all higher eukaryotes. However, the role of the m⁶A methylomes in rice is still poorly understood. With the development of the MeRIP-seq technique, the in-depth identification of mRNAs with m⁶A modification has become feasible. A study suggested that m⁶A modification is crucial for posttranscriptional regulation related to Cd²⁺-induced malignant transformation, but the association between m⁶A modification in plants and Cd tolerance has not been reported. We investigated the m⁶A methylomes in the roots of a cadmium (Cd)-treated group and compared them with the roots in the control (CK) group by m⁶A sequencing of cv. 9311 and cv. Nipponbare (NIP) plants. The results indicated that Cd leads to an altered modification profile in 3,406 differential m⁶A peaks in cv. 9311 and 2,065 differential m⁶A peaks in cv. NIP. KEGG pathway analysis of the genes with differentially modified m⁶A peaks indicated that the "phenylalanine", "tyrosine and tryptophan biosynthesis", "glycine", "adherens junctions", "glycerophospholipid metabolism" and "threonine metabolism" signalling pathways may be associated with the abnormal root development of cv. 9311 rice due to exposure to Cd. The "arginine", "proline metabolism", "glycerolipid", and "protein processing in endoplasmic reticulum" metabolism pathways were significantly enriched in genes with differentially modified m⁶A peaks in cv. NIP. Unlike that in Arabidopsis, the m⁶A-modified nucleotide position on mRNAs (m⁶A peak) distribution in rice exhibited a preference towards both the stop codon and 3' untranslated regions (3' UTRs). These findings provide a resource for plant RNA epitranscriptomic studies and further increase our knowledge on the function of m⁶A modification in RNA in plants.

Genetic drivers of m6A methylation in human brain, lung, heart and muscle

The most prevalent post-transcriptional mRNA modification, N⁶-methyladenosine (m⁶A), plays diverse RNA-regulatory roles, but its genetic control in human tissues remains uncharted. Here we report 129 transcriptome-wide m⁶A profiles, covering 91 individuals and 4 tissues (brain, lung, muscle and heart) from GTEx/eGTEx. We integrate these with interindividual genetic and expression variation, revealing 8,843 tissue-specific and 469 tissue-shared m⁶A quantitative trait loci (QTLs), which are modestly enriched in, but mostly orthogonal to, expression QTLs. We integrate m⁶A QTLs with disease genetics, identifying 184 GWAS-colocalized m⁶A QTL, including brain m⁶A QTLs underlying neuroticism, depression, schizophrenia and anxiety; lung m⁶A QTLs underlying expiratory flow and asthma; and muscle/heart m⁶A QTLs underlying coronary artery disease. Last, we predict novel m⁶A regulators that show preferential binding in m⁶A QTLs, protein interactions with known m⁶A regulators and expression correlation with the m⁶A levels of their targets. Our results provide important insights and resources for understanding both cis and trans regulation of epitranscriptomic modifications, their interindividual variation and their roles in human disease.

m6 A modification of HSATIII lncRNAs regulates temperature-dependent splicing

Nuclear stress bodies (nSBs) are nuclear membraneless organelles formed around stress-inducible HSATIII architectural long noncoding RNAs (lncRNAs). nSBs repress splicing of hundreds of introns during thermal stress recovery, which are partly regulated by CLK1 kinase phosphorylation of temperature-dependent Ser/Arg-rich splicing factors (SRSFs). Here, we report a distinct mechanism for this splicing repression through protein sequestration by nSBs. Comprehensive identification of RNA-binding proteins revealed HSATIII association with proteins related to N⁶ -methyladenosine (m⁶ A) RNA modification. 11% of the first adenosine in the repetitive HSATIII sequence were m⁶ A-modified. nSBs sequester the m⁶ A writer complex to methylate HSATIII, leading to subsequent sequestration of the nuclear m⁶ A reader, YTHDC1. Sequestration of these factors from the nucleoplasm represses m⁶ A modification of pre-mRNAs, leading to repression of m⁶ A-dependent splicing during stress recovery phase. Thus, nSBs serve as a common platform for regulation of temperature-dependent splicing through dual mechanisms employing two distinct ribonucleoprotein modules with partially m⁶ A-modified architectural lncRNAs.

Mapping single-nucleotide m6A by m6A-REF-seq

The past few years have witnessed rapid progress in the field of RNA modifications. As the most prevailing modification on eukaryotic mRNA, m⁶A is characterized to play a vital role in various cellular activities. However, limitations of the detection method impede functional studies of m⁶A. Here we introduce m⁶A-REF-seq, a powerful and straightforward method to identify m⁶A at single-nucleotide resolution. m⁶A-REF-seq relies on the recognition of RNA endonuclease MazF towards m⁶A at the ACA motif, providing an orthogonal method independent of the m⁶A antibody being adopted by most of current methods. We describe a detailed protocol to perform m⁶A-REF-seq, including NGS library construction and sequencing data analysis. In particular, we describe an optimized assay to validate individual m⁶A sites identified by m⁶A-REF-seq, which can also be applied to detect any candidate m⁶A sites.

An m6A-Related Prognostic Biomarker Associated With the Hepatocellular Carcinoma Immune Microenvironment

Background: N6-methyladenosine (m6A) is related to the progression of multiple cancers. However, the underlying influences of m6A-associated genes on the tumor immune microenvironment in hepatocellular carcinoma (HCC) remain poorly understood. Therefore, we sought to construct a survival prediction model using m6A-associated genes to clarify the molecular and immune characteristics of HCC.

Methods: HCC case data were downloaded from The Cancer Genome Atlas (TCGA). Then, by applying consensus clustering, we identified two distinct HCC clusters. Next, four m6A-related genes were identified to construct a prognostic model, which we validated with Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) datasets. Additionally, the molecular and immune characteristics in different subgroups were analyzed.

Results: m6A RNA methylation regulators were differentially expressed between HCC and normal samples and linked with immune checkpoint expression. Using consensus clustering, we divided HCC samples into two subtypes with distinct clinical features. Cluster 2 was associated with unfavorable prognosis, higher immune checkpoint expression and immune cell infiltration levels. In addition, the immune and carcinogenic signaling pathways were enriched in cluster 2. Furthermore, we constructed a risk model using four m6A-associated genes. Patients with different risk scores had distinct survival times, expression levels of immunotherapy biomarkers, TP53 mutation rates, and sensitivities to chemotherapy and targeted therapy. Similarly, the model exhibited an identical impact on overall survival in the validation cohorts.

Conclusion: The constructed m6A-based signature may be promising as a biomarker for prognostics and to distinguish immune characteristics in HCC.

METTL3 promotes tumour development by decreasing APC expression mediated by APC mRNA N6-methyladenosine-dependent YTHDF binding

The adenomatous polyposis coli (APC) is a frequently mutated tumour suppressor gene in cancers. However, whether APC is regulated at the epitranscriptomic level remains elusive. In this study, we analysed TCGA data and separated 200 paired oesophageal squamous cell carcinoma (ESCC) specimens and their adjacent normal tissues and demonstrated that methyltransferase-like 3 (METTL3) is highly expressed in tumour tissues. m⁶A-RNA immunoprecipitation sequencing revealed that METTL3 upregulates the m⁶A modification of APC, which recruits YTHDF for APC mRNA degradation. Reduced APC expression increases the expression of β-catenin and β-catenin-mediated cyclin D1, c-Myc, and PKM2 expression, thereby leading to enhanced aerobic glycolysis, ESCC cell proliferation, and tumour formation in mice. In addition, downregulated APC expression correlates with upregulated METTL3 expression in human ESCC specimens and poor prognosis in ESCC patients. Our findings reveal a mechanism by which the Wnt/β-catenin pathway is upregulated in ESCC via METTL3/YTHDF-coupled epitranscriptomal downregulation of APC.

The epitranscriptomic regulation of adenomatous polyposis coli (APC) tumour suppressor gene in cancers is unclear. Here the authors show that N6-methyladenosine methylation writer METTL3 downregulates APC by recruiting YTHDF2 for APC mRNA degradation, and that this promotes glycolysis and tumour growth in oesophageal cancers.

The m6A RNA Modification Quantity and mRNA Expression Level of RNA Methylation-Related Genes in Head and Neck Squamous Cell Carcinoma Cell Lines and Patients

RNA methylation at the nitrogen sixth of adenosine (m⁶A, N⁶-methyladenosine) is the most abundant RNA modification which plays a crucial role in all RNA metabolic aspects. Recently, m⁶A modification has been assigned to mediate the biological processes of cancer cells, but their significance in HNSCC development is still poorly described. Thus, the main aim of this study was to globally quantify m⁶A modification by the mass spectrometry approach and determine the mRNA expression level of selected m⁶A RNA methyltransferase (METTL3), demethylase (FTO), and m⁶A readers (YTHDF2, YTHDC2) in 45 HNSCC patients and 4 cell lines (FaDu, Detroit 562, A-253 and SCC-15) using qPCR. In the results, we have not observed differences in the global amount of m⁶A modification and the mRNA level of the selected genes between the cancerous and paired-matched histopathologically unchanged tissues from 45 HNSCC patients. However, we have found a positive correlation between selected RNA methylation machinery genes expression and m⁶A abundance on total RNA and characterized the transcript level of those genes in the HNSCC cell lines. Moreover, the lack of global m⁶A differences between cancerous and histopathologically unchanged tissues suggests that m⁶A alterations in specific RNA sites may specifically influence HNSCC tumorigenesis.

Transcriptome-Wide N6-Methyladenosine (m6A) Profiling of Susceptible and Resistant Wheat Varieties Reveals the Involvement of Variety-Specific m6A Modification Involved in Virus-Host Interaction Pathways

N⁶-methyladenosine (m⁶A) methylation is the most prevalent internal modification of post-transcriptional modifications in mRNA, tRNA, miRNA, and long non-coding RNA in eukaryotes. m⁶A methylation has been proven to be involved in plant resistance to pathogens. However, there are no reports on wheat (Triticum aestivum) m⁶A transcriptome-wide map and its potential biological function in wheat resistance to wheat yellow mosaic virus (WYMV). To the best of our knowledge, this study is the first to determine the transcriptome-wide m⁶A profile of two wheat varieties with different resistances to WYMV. By analyzing m⁶A-sequencing (m⁶A-seq) data, we identified 25,752 common m⁶A peaks and 30,582 common m⁶A genes in two groups [WYMV-infected resistant wheat variety (WRV) and WYMV-infected sensitive wheat variety (WSV)], and all these peaks were mainly enriched in 3′ untranslated regions and stop codons of coding sequences. Gene Ontology analysis of m⁶A-seq and RNA-sequencing data revealed that genes that showed significant changes in both m⁶A and mRNA levels were associated with plant defense responses. Kyoto Encyclopedia of Genes and Genomes analysis revealed that these selected genes were enriched in the plant–pathogen interaction pathway. We further verified these changes in m⁶A and mRNA levels through gene-specific m⁶A real-time quantitative PCR (RT-qPCR) and normal RT-qPCR. This study highlights the role of m⁶A methylation in wheat resistance to WYMV, providing a solid basis for the potential functional role of m⁶A RNA methylation in wheat resistance to infection by RNA viruses.

Identification of an m6A-Related Signature as Biomarker for Hepatocellular Carcinoma Prognosis and Correlates with Sorafenib and Anti-PD-1 Immunotherapy Treatment Response

Background

N6-methyladenosine (m6A) modification plays an essential role in diverse key biological processes and may take part in the development and progression of hepatocellular carcinoma (HCC). Here, we systematically analyzed the expression profiles and prognostic values of 13 widely reported m6A modification-related genes in HCC.

Methods

The mRNA expression of 13 m6A modification-related genes and clinical parameters of HCC patients were downloaded from TCGA, ICGC, GSE109211, and GSE78220. Univariate and LASSO analyses were used to develop risk signature. Time-dependent ROC was performed to assess the predictive accuracy and sensitivity of risk signature.

Results

FTO, YTHDC1, YTHDC2, ALKBH5, KIAA1429, HNRNPC, METTL3, RBM15, YTHDF2, YTHDF1, and WTAP were significantly overexpressed in HCC patients. YTHDF1, HNRNPC, RBM15, METTL3, and YTHDF2 were independent prognostic factors for OS and DFS in HCC patients. Next, a risk signature was also developed and validated with five m6A modification-related genes in TCGA and ICGC HCC cohort. It could effectively stratify HCC patients into high-risk patients with shorter OS and DFS and low-risk patients with longer OS and DFS and showed good predictive efficiency in predicting OS and DFS. Moreover, significantly higher proportions of macrophages M0 cells, neutrophils, and Tregs were found to be enriched in HCC patients with high risk scores, while significantly higher proportions of memory CD4 T cells, gamma delta T cells, and naive B cells were found to be enriched in HCC patients with low scores. Finally, significantly lower risk scores were found at sorafenib treatment responders and anti-PD-1 immunotherapy responders compared to that in nonresponders, and anti-PD-1 immunotherapy-treated patients with lower risk scores had better OS than patients with higher risk scores.

Conclusion

A risk signature developed with the expression of 5 m6A-related genes could improve the prediction of prognosis of HCC and correlated with sorafenib treatment and anti-PD-1 immunotherapy response.

FTO (Fat-Mass and Obesity-Associated Protein) Participates in Hemorrhage-Induced Thalamic Pain by Stabilizing Toll-Like Receptor 4 Expression in Thalamic Neurons

[Figure: see text].

Systematic analysis of molecular characterization and clinical relevance of m6A regulators in digestive system pan-cancers

Digestive system tumors, which mainly include esophagus, stomach, colorectum, liver, pancreas, bile duct, and some other tumors, often have a poor prognosis. N6-methyladenosine (m6A) has critical functions in development and tumorigenesis and may help improve the molecular mechanisms of digestive system tumors. However, current understanding of the reconstitution of m6A in digestive system tumors is far from comprehensive. Herein, this study systematically analyzed multi-layered genomic characteristics and clinical relevance of m6A regulators in 1906 patients involving seven digestive system tumor types. We discovered that m6A regulators showed extensive genetic changes and highly consistent expression regulation. The m6A expression was closely related to the activity of cancer pathways. At the same time, we also identified m6A regulators significantly related to the common cancer pathways of digestive system tumors and specific cancer pathways of digestive tract and digestive glands. These cancer pathways may explain the prognostic differences of patients with digestive tract tumors. In addition, m6A regulators demonstrated strong potential in prognostic stratification and drug development, especially in multiple research cohorts on pancreatic cancer, pointing to a strong prognostic stratification capability of m6A regulators. Finally, a m6A scoring model significantly related to highly active ubiquitin-mediated proteolysis, mismatch repair, cell cycle, ebasal transcription factors was constructed and had a strong prognostic stratification ability in digestive gland tumors. The score showed a significant negative correlation with the tumor immune microenvironment. This study demonstrated that the similarities and difference of the action mechanism m6A regulators in the digestive tract and digestive gland tumor progression could guide potential drug development.

Long Non-Coding RNA Epigenetics

Long noncoding RNAs exceeding a length of 200 nucleotides play an important role in ensuring cell functions and proper organism development by interacting with cellular compounds such as miRNA, mRNA, DNA and proteins. However, there is an additional level of lncRNA regulation, called lncRNA epigenetics, in gene expression control. In this review, we describe the most common modified nucleosides found in lncRNA, 6-methyladenosine, 5-methylcytidine, pseudouridine and inosine. The biosynthetic pathways of these nucleosides modified by the writer, eraser and reader enzymes are important to understanding these processes. The characteristics of the individual methylases, pseudouridine synthases and adenine-inosine editing enzymes and the methods of lncRNA epigenetics for the detection of modified nucleosides, as well as the advantages and disadvantages of these methods, are discussed in detail. The final sections are devoted to the role of modifications in the most abundant lncRNAs and their functions in pathogenic processes.

m6A modification regulates lung fibroblast-to-myofibroblast transition through modulating KCNH6 mRNA translation

Idiopathic pulmonary fibrosis (IPF) is a chronic, fatal lung disease characterized by progressive and non-reversible abnormal matrix deposition in lung parenchyma. Myofibroblasts origin mainly from resident fibroblasts via fibroblast-to-myofibroblast transition (FMT) are the dominant collagen-producing cells in pulmonary fibrosis. N⁶-methyladenosine (m⁶A) modification has been implicated in various biological process. However, the role of m⁶A modification in pulmonary fibrosis remains elusive. In this study, we reveal that m⁶A modification is up-regulated in bleomycin induced pulmonary fibrosis mice model, FMT-derived myofibroblasts and idiopathic pulmonary fibrosis patient lung samples. Lowering m⁶A level through silencing METTL3 inhibits FMT process in vitro and vivo. Mechanistically, KCNH6 is involved in m⁶A-regulated FMT process. m⁶A modification regulates the expression of KCNH6 by modulating its translation in a YTHDF1 dependent manner. Together, our study highlights the critical role of m⁶A modification in pulmonary fibrosis. Manipulation of m⁶A modification through targeting METTL3 may become a promising strategy for the treatment of pulmonary fibrosis.

KAT1 triggers YTHDF2-mediated ITGB1 mRNA instability to alleviate the progression of diabetic retinopathy

Diabetic retinopathy (DR) is a major complication of diabetes and a leading cause of blindness and visual impairment. This study focuses on the function of lysine acetyltransferase 1 (KAT1) in the progression of DR and the epigenetic mechanism. A mouse model with DR was induced by streptozotocin (STZ). Abundantly expressed genes in STZ-induced mice were analyzed. KAT1 was found to be significantly downregulated in the retinal tissues of model mice. Retinal microvascular endothelial cells (RMECs) and retinal Müller cells (rMCs) were cultured in high-glucose medium for in vitro studies. Upregulation of KAT1 suppressed inflammation, neovascularization, and vascular leakage in mouse retinal tissues, and it reduced the activity and inflammatory responses in rMCs, as well as the proliferation and metastatic potential of RMECs. KAT1 activated the transcription activity of YTHDF2 through histone acetylation of the promoter, and YTHDF2 triggered the instability of ITGB1 mRNA to induce mRNA degradation in an m6A manner. The activities of rMCs and RMECs were increased by sh-YTHDF2 but suppressed by sh-ITGB1. The FAK/PI3K/AKT signaling pathway was suppressed upon ITGB1 silencing. Collectively, this study demonstrated that KAT1 triggers YTHDF2-mediated ITGB1 mRNA instability to alleviate the progression of DR.

The m6A-related gene signature for predicting the prognosis of breast cancer

N6-methyladenosine (m⁶A) modification has been shown to participate in tumorigenesis and metastasis of human cancers. The present study aimed to investigate the roles of m⁶A RNA methylation regulators in breast cancer. We used LASSO regression to identify m⁶A-related gene signature predicting breast cancer survival with the datasets downloaded from Gene Expression Omnibus and The Cancer Genome Atlas (TCGA). RNA-Seq data of 3409 breast cancer patients from GSE96058 and 1097 from TCGA were used in present study. A 10 m⁶A-related gene signature associated with prognosis was identified from 22 m⁶A RNA methylation regulators. The signature divided patients into low- and high-risk group. High-risk patients had a worse prognosis than the low-risk group. Further analyses indicated that IGF2BP1 may be a key m⁶A RNA methylation regulator in breast cancer. Survival analysis showed that IGF2BP1 is an independent prognostic factor of breast cancer, and higher expression level of IGF2BP1 is associated with shorter overall survival of breast cancer patients. In conclusion, we identified a 10 m⁶A-related gene signature associated with overall survival of breast cancer. IGF2BP1 may be a key m⁶A RNA methylation regulator in breast cancer.

Pan-Cancer Prognostic, Immunity, Stemness, and Anticancer Drug Sensitivity Characterization of N6-Methyladenosine RNA Modification Regulators in Human Cancers

N6-methyladenosine RNA modification plays a significant role in the progression of multiple tumorigenesis. Our study identified the imperative role of m6A regulators in the tumor immune microenvironment, survival, stemness score, and anticancer drug sensitivity of pan-cancer. The Wilcox test was to identify the differential expression between 17 m6A regulators across 33 TCGA cancer types and their normal tissues from UCSC Xena GDC pan-cancer. Survival analysis of m6A-related regulators in 33 TCGA cancer types was identified using the "survival" and "survminer" package. The Spearman correlation test and Pearson correlation test were used to identify the correlation relationship between m6A regulators expression and tumor microenvironment, tumor stem cell score, and drug sensitivity of anticancer drugs. ConsensusPathDB was used for exploring m6A regulators functional enrichment. The 17 (METTL3, WTAP, METTL14, RBM15, RBM15B, VIRMA, HNRNPC, HNRNPA2B1, YTHDC1, ZC3H13, YTHDF1, YTHDC2, YTHDF2, IGF2BP3, IGF2BP1, FTO, and ALKBH5) m6A regulators were differentially expressed in 18 TCGA cancer types and adjacent normal tissues. Correlation analysis indicated that the relationship between the expression of 17 m6A regulators and tumor microenvironment indicated that the higher expression of m6A regulators, the higher the degree of tumor stem cells. The anticancer drug sensitivity analysis indicated that ZC3H13 expression had a positive relationship with anticancer drugs such as selumetinib, dabrafenib, cobimetinib, trametinib, and hypothemycin (p < 0.001). YTHDF2 expression was significantly negatively correlated with the anticancer drug dasatinib (p < 0.001). The pan-cancer immune subtype analysis showed that the 17 m6A regulators were significantly different in immune subtype C1 (wound healing), C3 (inflammatory), C2 (IFN-gamma dominant), C5 (immunological quiet), C4 (lymphocyte depleted), and C6 (TGF-beta dominant) (p < 0.001). Our study provides a comprehensive insight for revealing the significant role of m6A regulators in the tumor immune microenvironment, stemness score, and anticancer drug sensitivity of human cancers.

YTHDF2 mediates LPS-induced osteoclastogenesis and inflammatory response via the NF-κB and MAPK signaling pathways

Aberrant elevation of osteoclast differentiation and function is responsible for disrupting bone homeostasis in various inflammatory bone diseases. YTH domain family 2 (YTHDF2) is a well-known m⁶A-binding protein that plays an essential role in regulating cell differentiation and inflammatory processes by mediating mRNA degradation. However, the regulatory role of YTHDF2 in inflammatory osteoclast differentiation remains unelucidated. Here, we detected the expression of m⁶A-related genes and found that YTHDF2 was upregulated in RANKL-primed osteoclast precursors stimulated with lipopolysaccharide (LPS). Ythdf2 knockdown in RAW264.7 cells and primary bone marrow-derived macrophages (BMMs) enhanced osteoclast formation and bone resorption, which was assessed by TRAP staining assay and pit formation assay. Ythdf2 depletion upregulated osteoclast-related gene expression and proinflammatory cytokine secretion. In contrast, overexpression of Ythdf2 produced the reverse effect. Furthermore, Ythdf2 knockdown enhanced the phosphorylation of IKKα/β, IκBα, ERK, P38 and JNK. NF-κB and MAPK signaling pathway inhibitors effectively abrogated the enhanced expression of Nfact1, c-Fos, IL-1β and TNF-α caused by Ythdf2 knockdown. Mechanistically, the mRNA stability assay revealed that Ythdf2 depletion led to stabilization of Tnfrsf11a, Traf6, Map4k4, Map2k3, Map2k4 and Nfatc1 mRNA. In summary, our findings demonstrated that YTHDF2 has a negative regulatory role in LPS-induced osteoclast differentiation and the inflammatory response via the NF-κB and MAPK signaling pathways.

m6 A modification of lncRNA PCAT6 promotes bone metastasis in prostate cancer through IGF2BP2-mediated IGF1R mRNA stabilization

BACKGROUND

Bone metastasis is the leading cause of tumor-related death in prostate cancer (PCa) patients. Long noncoding RNAs (lncRNAs) have been well documented to be involved in the progression of multiple cancers. Nevertheless, the role of lncRNAs in PCa bone metastasis remains largely unclear.

METHODS

The expression of prostate cancer-associated transcripts was analyzed in published datasets and further verified in clinical samples and cell lines by RT-qPCR and in situ hybridization assays. Colony formation assay, MTT assay, cell cycle analysis, EdU assay, Transwell migration and invasion assays, wound healing assay, and in vivo experiments were carried out to investigate the function of prostate cancer-associated transcript 6 (PCAT6) in bone metastasis and tumor growth of PCa. Bioinformatic analysis, RNA pull-down, and RIP assays were conducted to identify the proteins binding to PCAT6 and the potential targets of PCAT6. The therapeutic potential of targeting PCAT6 by antisense oligonucleotides (ASO) was further explored in vivo.

RESULTS

PCAT6 was upregulated in PCa tissues with bone metastasis and increased PCAT6 expression predicted poor prognosis in PCa patients. Functional experiments found that PCAT6 knockdown significantly inhibited PCa cell invasion, migration, and proliferation in vitro, as well as bone metastasis and tumor growth in vivo. Mechanistically, METTL3-mediated m6 A modification contributed to PCAT6 upregulation in an IGF2BP2-dependent manner. Furthermore, PCAT6 upregulated IGF1R expression by enhancing IGF1R mRNA stability through the PCAT6/IGF2BP2/IGF1R RNA-protein three-dimensional complex. Importantly, PCAT6 inhibition by ASO in vivo showed therapeutic potential against bone metastasis in PCa. Finally, the clinical correlation of METTL3, IGF2BP2, IGF1R, and PCAT6 was further demonstrated in PCa tissues and cells.

CONCLUSIONS

Our study uncovers a novel molecular mechanism by which the m6 A-induced PCAT6/IGF2BP2/IGF1R axis promotes PCa bone metastasis and tumor growth, suggesting that PCAT6 may serve as a promising prognostic marker and therapeutic target against bone-metastatic PCa.

Importance of N6-methyladenosine RNA modification in lung cancer (Review)

The N6-methyladenosine (m6A) modification is the most common mRNA modification in eukaryotes and exerts biological functions by affecting RNA metabolism. The m6A modification is installed by m6A methyltransferases, removed by demethylases and recognized by m6A-binding proteins. The interaction between these three elements maintains the dynamic equilibrium of m6A in cells. Accumulating evidence indicates that m6A RNA methylation has a significant impact on RNA metabolism and is involved in the pathogenesis of cancer. Lung cancer is the leading cause of cancer-related deaths worldwide. The treatment options for lung cancer have developed considerably over the past few years; however, the survival rate of patients with lung cancer still remains very low. Although diagnostic methods and targeted therapies have been rapidly developed in recent years, the underlying mechanism and importance of m6A RNA methylation in the pathogenesis of lung cancer remains ambiguous. The current review summarized the biological functions of m6A modification and considers the potential roles of m6A regulators in the occurrence and development of lung cancer.

The Role of N6 -Methyladenosine Modified Circular RNA in Pathophysiological Processes

Circular RNA (circRNA) is a type of covalently closed and endogenous non-coding RNA (ncRNA) with tissue- and cell-specific expression patterns generated by a non-canonical splicing event. Previous reports have indicated that circRNAs exert their functions in different ways, thereby participating in various pathophysiological processes. N6 -methyladenosine (m6A) methylation occurs in the N6-position, which is the most abundant and conserved internal transcriptional modification in eukaryotes, including mRNA and ncRNAs. Accumulating evidences confirm that m6A modification also exists in the circRNA and greatly affects the biological functions of circRNA. Their dysregulated expression can be a cause of various pathophysiological processes, such as spermatogenesis, myoblast differentiation, cancer, cardiovascular disease, mental illness and so on. Understanding the role of m6A-modified circRNAs in pathophysiological processes may contribute to better understanding the physiological mechanisms and develop new biomarkers. This review summarizes the regulatory mechanism of m6A modification on circRNA metabolism and the role of m6A-modified circRNAs in pathophysiological processes. This article may pave the way for a better understanding of the role of epigenetically modified circRNAs in pathophysiological process.

Comprehensive Profiling Reveals Distinct Microenvironment and Metabolism Characterization of Lung Adenocarcinoma

Lung adenocarcinoma has entered into an era of immunotherapy with the development of immune checkpoint inhibitors (ICIs). The identification of immune subtype is crucial to prolonging survival in patients. The tumor microenvironment (TME) and metabolism have a profound impact on prognosis and therapy. The majority of previous studies focused on only one aspect, while both of them are essential to the understanding of tumorigenesis and development. We hypothesized that lung adenocarcinoma can be stratified into immune subgroups with alterations in the TME infiltration. We aimed to explore the "TME-Metabolism-Risk" patterns in each subtypes and the mechanism behind. Glycolysis and cholesterol were selected for the analysis of metabolic states based on the first half of the study. Bioinformatic analysis was performed to investigate the transcriptomic and clinical data integrated by three lung adenocarcinoma cohorts (GSE30219, GSE31210, GSE37745, N = 415). The results were validated in an independent cohort (GSE50081, N = 127). In total, 415 lung adenocarcinoma samples were integrated and analyzed. Four major immune subtypes were indentified using bioinformatic analysis. Subtype NC1, characterized by a high level of glycolysis, with extremely low microenvironment cell infiltration. Subtype NC2, characterized by the "Silence" and "Cholesterol biosynthesis Predominant" metabolic states, with a middle degree infiltration of microenvironment cell. Subtype NC3, characterized by the lack of "Cholesterol biosynthesis Predominant" metabolic state, with abundant microenvironment cell infiltration. Subtype NC4, characterized by "Mixed" metabolic state, with a relatively low microenvironment cell infiltration. Least absolute shrinkage and selection operator (LASSO) regression and multivariate analyses were performed to calculate the risk of each sample, and we attempted to find out the potential immune escape mechanism in different subtypes. The result revealed that the lack of immune cells infiltration might contribute to the immune escape in subtypes NC1 and NC4. NC3 was characterized by the high expression of immune checkpoint molecules and fibroblasts. NC2 had defects in activation of innate immune cells. There existed an obviously survival advantage in subtype NC2. Gene set enrichment analysis (GSEA) and Gene Ontology analysis indicated that the PI3K-AKT-mTOR, TGF-β, MYC-related pathways might be correlated with this phenomenon. In addition, some differentially expressed genes (DEGs) were indentified in subtype NC3, which might be potential targets for survival phenotype transformation.

Rational Design of Novel Anticancer Small-Molecule RNA m6A Demethylase ALKBH5 Inhibitors

The RNA 6-N-methyladenosine (m6A) demethylase ALKBH5 has been shown to be oncogenic in several cancer types, including leukemia and glioblastoma. We present here the target-tailored development and first evaluation of the antiproliferative effects of new ALKBH5 inhibitors. Two compounds, 2-[(1-hydroxy-2-oxo-2-phenylethyl)sulfanyl]acetic acid (3) and 4-{[(furan-2-yl)methyl]amino}-1,2-diazinane-3,6-dione (6), with IC50 values of 0.84 μM and 1.79 μM, respectively, were identified in high-throughput virtual screening of the library of 144 000 preselected compounds and subsequent verification of hits in an m6A antibody-based enzyme-linked immunosorbent assay (ELISA) enzyme inhibition assay. The effect of these compounds on the proliferation of selected target cancer cell lines was then measured. In the case of three leukemia cell lines (HL-60, CCRF-CEM, and K562) the cell proliferation was suppressed at low micromolar concentrations of inhibitors, with IC50 ranging from 1.38 to 16.5 μM. However, the effect was low or negligible in the case of another leukemia cell line, Jurkat, and the glioblastoma cell line A-172. These results demonstrate the potential of ALKBH5 inhibition as a cancer-cell-type-selective antiproliferative strategy.

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression and Kaplan–Meier (K–M) analysis were performed to identify m6A methyltransferase-related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes (BPs) and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.

METTL3 mediates bone marrow mesenchymal stem cell adipogenesis to promote chemoresistance in acute myeloid leukaemia

Adipogenesis of bone marrow mesenchymal stem cells (MSCs) promotes chemoresistance of acute myeloid leukaemia (AML) cells. MSCs from AML patients (AML‐MSCs) display enhanced adipogenesis compared with bone marrow MSCs from healthy donors. However, the precise molecular mechanism by which adipogenesis of MSCs from AML marrow differs from normal counterparts remains obscure. We found that METTL3 significantly inhibits MSC adipogenesis. Here, we aimed to identify the molecular mechanism linking METTL3 and MSC adipogenesis. Analysis of m⁶A epigenetic changes in MSCs determined via RIP‐qPCR and MeRIP‐qPCR indicated that METTL3 affects AKT protein expression in MSCs by mediating m⁶A modification of AKT1‐mRNA. Downregulated METTL3 expression in AML‐MSCs induced an increase in AKT protein, resulting in enhanced MSC adipogenesis, thereby contributing to chemoresistance in AML cells. Therefore, targeting AKT regulation by mRNA modification in MSC adipogenesis might provide a novel therapeutic strategy to overcome AML chemoresistance.

mTORC1 stimulates cell growth through SAM synthesis and m6A mRNA-dependent control of protein synthesis

The mechanistic target of rapamycin complex 1 (mTORC1) regulates metabolism and cell growth in response to nutrient, growth, and oncogenic signals. We found that mTORC1 stimulates the synthesis of the major methyl donor, S-adenosylmethionine (SAM), through the control of methionine adenosyltransferase 2 alpha (MAT2A) expression. The transcription factor c-MYC, downstream of mTORC1, directly binds to intron 1 of MAT2A and promotes its expression. Furthermore, mTORC1 increases the protein abundance of Wilms' tumor 1-associating protein (WTAP), the positive regulatory subunit of the human N6-methyladenosine (m6A) RNA methyltransferase complex. Through the control of MAT2A and WTAP levels, mTORC1 signaling stimulates m6A RNA modification to promote protein synthesis and cell growth. A decline in intracellular SAM levels upon MAT2A inhibition decreases m6A RNA modification, protein synthesis rate, and tumor growth. Thus, mTORC1 adjusts m6A RNA modification through the control of SAM and WTAP levels to prime the translation machinery for anabolic cell growth.

Regulatory Role of N6-methyladenosine (m6A) Modification in Osteosarcoma

Osteosarcoma is the most common primary bone malignancy, typically occurring in childhood or adolescence. Unfortunately, the clinical outcomes of patients with osteosarcoma are usually poor because of the aggressive nature of this disease and few treatment advances in the past four decades. N6-methyladenosine (m⁶A) is one of the most extensive forms of RNA modification in eukaryotes found both in coding and non-coding RNAs. Accumulating evidence suggests that m⁶A-related factors are dysregulated in multiple osteosarcoma processes. In this review, we highlight m⁶A modification implicated in osteosarcoma, describing its pathophysiological role and molecular mechanism, as well as future research trends and potential clinical application in osteosarcoma.

N6-methyladenosine (m6A) RNA methylation signature as a predictor of stomach adenocarcinoma outcomes and its association with immune checkpoint molecules

OBJECTIVE

Although N⁶-methyladenosine (m⁶A) RNA methylation is the most common mRNA modification process, few studies have examined the role of m⁶A in stomach adenocarcinomas (STADs).

METHODS

In this retrospective study, we analyzed 293 STAD samples from The Cancer Genome Atlas with complete clinicopathological feature profiles. The m⁶A methylation risk signature was derived from LASSO-Cox regression analyses with 15 m⁶A regulators. Statistical analysis was performed and figures were prepared using R software (https://www.R-project.org/).

RESULTS

The m⁶A signature was established as follows: risk score = FTO × 0.127 + YTHDF1 × 0.004 + KIAA1429 × 0.044 + YTHDC2 × 0.112 - RBM15 × 0.135 - ALKBH5 × 0.019 - YTHDF2 × 0.028, which was confirmed as an independent prognostic indicator to predict overall survival of patients with STAD. Risk scores and tumor grades were closely associated. Cell cycle, p53 signaling pathways, DNA mismatch repair, and RNA degradation were enriched in the low-risk subgroup. This subgroup showed significantly higher expression of immune checkpoint molecules including PD-1 (programmed death 1), PD-L1 (programmed death-ligand 1), and CTLA-4 (cytotoxic T-lymphocyte-associated antigen 4), suggesting that the signature may be a useful immunotherapy predictor.

CONCLUSIONS

We established an m⁶A methylation signature as an independent prognostic tool to predict overall survival, which may also be useful as an immunotherapy predictor.

WNT/β-catenin-suppressed FTO expression increases m6A of c-Myc mRNA to promote tumor cell glycolysis and tumorigenesis

FTO removes the N6-methyladenosine (m⁶A) modification from genes and plays a critical role in cancer development. However, the mechanisms underlying the regulation of FTO and its subsequent impact on the regulation of the epitranscriptome remain to be further elucidated. Here, we demonstrate that FTO expression is downregulated and inversely correlated with poor survival of lung adenocarcinoma patients. Mechanistically, Wnt signaling induces the binding of EZH2 to β-catenin. This protein complex binds to the LEF/TCF-binding elements at the promoter region of FTO, where EZH2 enhances H3K27me3 and inhibits FTO expression. Downregulated FTO expression substantially enhances the m⁶A levels in the mRNAs of a large number of genes in critical pathways, particularly metabolic pathway genes, such as MYC. Enhanced m⁶A levels on MYC mRNA recruit YTHDF1 binding, which promotes MYC mRNA translation and a subsequent increase in glycolysis and proliferation of tumor cells and tumorigenesis. Our findings uncovered a critical mechanism of epitranscriptome regulation by Wnt/β-catenin-mediated FTO downregulation and underscored the role of m⁶A modifications of MYC mRNA in regulating tumor cell glycolysis and growth.

METTL3 regulates skeletal muscle specific miRNAs at both transcriptional and post-transcriptional levels

METTL3 increasing the mature miRNA levels via N6-Methyladenosine (m6A) modification of primary miRNA (pri-miRNA) transcripts has emerged as an important post-transcriptional regulation of miRNA biogenesis. Our previous studies and others have showed that muscle specific miRNAs are essential for skeletal muscle differentiation. Whether these miRNAs are also regulated by METTL3 is still unclear. Here, we found that m6A motifs were present around most of these miRNAs, which were indeed m6A modified as confirmed by m6A-modified RNA immunoprecipitation (m6A RIP). However, we surprisingly found that these muscle specific miRNAs were repressed instead of increased by METTL3 in C2C12 in vitro differentiation and mouse skeletal muscle regeneration after injury in vivo model. To elucidate the underlined mechanism, we performed reporter assays in 293T cells and validated METTL3 increasing these miRNAs at post-transcriptional level as expected. Furthermore, in myogenic C2C12 cells, we found that METTL3 not only repressed the expression of myogenic transcription factors (TFs) which can enhance the muscle specific miRNAs, but also increased the expression of epigenetic regulators which can repress these miRNAs. Thus, METTL3 could repress the muscle specific miRNAs at transcriptional level indirectly. Taken together, our results demonstrated that skeletal muscle specific miRNAs were repressed by METTL3 and such repression is likely synthesized transcriptional and post-transcriptional regulations.

N6 -methyladenosine Steers RNA Metabolism and Regulation in Cancer

As one of the most studied ribonucleic acid (RNA) modifications in eukaryotes, N⁶ -methyladenosine (m⁶ A) has been shown to play a predominant role in controlling gene expression and influence physiological and pathological processes such as oncogenesis and tumor progression. Writer and eraser proteins, acting opposite to deposit and remove m⁶ A epigenetic marks, respectively, shape the cellular m⁶ A landscape, while reader proteins preferentially recognize m⁶ A modifications and mediate fate decision of the methylated RNAs, including RNA synthesis, splicing, exportation, translation, and stability. Therefore, RNA metabolism in cells is greatly influenced by these three classes of m⁶ A regulators. Aberrant expression of m⁶ A regulators has been widely reported in various types of cancer, leading to cancer initiation, progression, and drug resistance. The close links between m⁶ A and cancer shed light on the potential use of m⁶ A methylation and its regulators as prognostic biomarkers and drug targets for cancer therapy. Given the notable effects of m⁶ A in reversing chemoresistance and enhancing immune therapy, it is a promising target for combined therapy. Herein, we summarize the recent discoveries on m⁶ A and its regulators, emphasizing their influences on RNA metabolism, their dysregulation and impacts in diverse malignancies, and discuss the clinical implications of m⁶ A modification in cancer.

RNA modifications in cardiovascular diseases, the potential therapeutic targets

More than one hundred RNA modifications decorate the chemical and topological properties of these ribose nucleotides, thereby executing their biological functions through post-transcriptional regulation. In cardiovascular diseases, a wide range of RNA modifications including m⁶A (N6-adenosine methylation), m⁵C (5-methylcytidin), Nm (2'-O-ribose-methylation), Ψ (pseudouridine), m⁷G (N7-methylguanosine), and m¹A (N1-adenosine methylation) have been found in tRNA, rRNA, mRNA and other noncoding RNA, which can function as a novel mechanism in metabolic syndrome, heart failure, coronary heart disease, and hypertension. In this review, we will summarize the current understanding of the regulatory roles and significance of several types of RNA modifications in CVDs (cardiovascular diseases) and the interplay between RNA modifications and noncoding RNA, epigenetics. Finally, we will focus on the potential therapeutic strategies by using RNA modifications.

Systematic expression analysis of m 6 A RNA methyltransferases in clear cell renal cell carcinoma

Objectives

To investigate the regulation of the N‐6‐methyladenosine (m⁶A) methyltransferases METTL3, METTL14, WTAP, KIAA1429, and METTL4, referred to as “m⁶A writers,” in clear cell renal cell carcinoma (ccRCC), and other RCC subtypes in respect of the potential prognostic value.

Patients and methods

Tissue samples were collected within the framework of the Biobank at the Center for Integrated Oncology Bonn. The expression of the methyltransferases was systematically determined in clear cell renal carcinoma (ccRCC) on the RNA (real‐time PCR) and protein level (immunohistochemistry). Additionally, protein expression of the m⁶A writers was further investigated in papillary RCC, chromophobe RCC, sarcomatoid RCC, oncocytoma, and normal renal tissue (immunohistochemistry).

Results

The expression of all m⁶A‐methyltransferases was significantly downregulated in ccRCC compared to benign renal tissue. Low m⁶A‐methyltransferase levels were correlated with higher histological grade, advanced pT‐stage, pN‐stage, and metastatic disease. Reduced m⁶A‐methyltransferase expression was associated with shorter overall survival.

Conclusion

In conclusion, m⁶A‐methyltransferases are dysregulated in ccRCC and might act as tumor suppressor genes, which could be of particular importance for future diagnostic and therapeutic options.

Loss of m6A demethylase ALKBH5 promotes post-ischemic angiogenesis via post-transcriptional stabilization of WNT5A

BACKGROUND

Post-ischemic angiogenesis is critical for blood flow recovery and ischemic tissue repair. N6-methyladenosine (m6A) plays essential roles in numerous biological processes. However, the impact and connected mechanism of m6A on post-ischemic angiogenesis are not fully understood.

METHODS

AlkB homolog 5 (ALKBH5) was screened out among several methyltransferases and demethylases involved in dynamic m6A regulation. Cardiac microvascular endothelial cells (CMECs) angiogenesis and WNT family member 5A (WNT5A) stability were analyzed upon ALKBH5 overexpression with adenovirus or knockdown with small interfering RNAs in vitro. The blood flow recovery, capillary, and small artery densities were evaluated in adeno-associated virus (AAV)-ALKBH5 overexpression or ALKBH5 knockout (KO) mice in a hind-limb ischemia model. The same experiments were conducted to explore the translational value of transient silencing of ALKBH5 with adenovirus.

RESULTS

ALKBH5 was significantly upregulated in hypoxic CMECs and led to a global decrease of m6A level. ALKBH5 overexpression further reduced m6A level in normoxic and hypoxic CMECs, impaired proliferation, migration, and tube formation only in hypoxic CMECs. Conversely, ALKBH5 knockdown preserved m6A levels and promoted angiogenic phenotypes in hypoxic but not in normoxic CMECs. Mechanistically, ALKBH5 regulated WNT5A expression through post-transcriptional mRNA modulation in an m6A-dependent manner, which decreased its stability and subsequently impeded angiogenesis in hypoxic CMECs. Furthermore, ALKBH5 overexpression hindered blood flow recovery and reduced CD31 and alpha-smooth muscle actin expression in hind-limb ischemia mice. As expected, ALKBH5-KO mice exhibited improved blood flow recovery, increased capillary, and small artery densities after hind-limb ischemia, and similar beneficial effects were observed in mice with transient adenoviral ALKBH5 gene silencing.

CONCLUSION

We demonstrate that ALKBH5 is a negative regulator of post-ischemic angiogenesis via post-transcriptional modulation and destabilization of WNT5A mRNA in an m6A-dependent manner. Targeting ALKBH5 may be a potential therapeutic option for ischemic diseases, including peripheral artery disease.

Epitranscriptomic regulation of insecticide resistance

N 6-methyladenosine (m6A) is the most prevalent messenger RNA modification in eukaryotes and an important posttranscriptional regulator of gene expression. However, the biological roles of m6A in most insects remain largely unknown. Here, we show that m6A regulates a cytochrome P450 gene (CYP4C64) in the global whitefly pest, Bemisia tabaci, leading to insecticide resistance. Investigation of the regulation of CYP4C64, which confers resistance to the insecticide thiamethoxam, revealed a mutation in the 5' untranslated region of this gene in resistant B. tabaci strains that introduces a predicted m6A site. We provide several lines of evidence that mRNA methylation of the adenine at this position, in combination with modified expression of m6A writers, acts to increase expression of CYP4C64 and resistance. Collectively, these results provide an example of the epitranscriptomic regulation of the xenobiotic response in insects and implicate the m6A regulatory axis in the development of insecticide resistance.

MMP24 Contributes to Neuropathic Pain in an FTO-Dependent Manner in the Spinal Cord Neurons

Nerve injury-induced gene expression change in the spinal cord is critical for neuropathic pain genesis. RNA N⁶-methyladenosine (m⁶A) modification represents an additional layer of gene regulation. We showed that spinal nerve ligation (SNL) upregulated the expression of matrix metallopeptidase 24 (MMP24) protein, but not Mmp24 mRNA, in the spinal cord neurons. Blocking the SNL-induced upregulation of spinal MMP24 attenuated local neuron sensitization, neuropathic pain development and maintenance. Conversely, mimicking MMP24 increase promoted the spinal ERK activation and produced evoked nociceptive hypersensitivity. Methylated RNA Immunoprecipitation Sequencing (MeRIP-seq) and RNA Immunoprecipitation (RIP) assay indicated the decreased m⁶A enrichment in the Mmp24 mRNA under neuropathic pain condition. Moreover, fat-mass and obesity-associated protein (FTO) was colocalized with MMP24 in spinal neurons and shown increased binding to the Mmp24 mRNA in the spinal cord after SNL. Overexpression or suppression of FTO correlates with promotion or inhibition of MMP24 expression in cultured spinal cord neurons. In conclusion, SNL promoted the m⁶A eraser FTO binding to the Mmp24 mRNA, which subsequently facilitated the translation of MMP24 in the spinal cord, and ultimately contributed to neuropathic pain genesis.

METTL3 enhances the stability of MALAT1 with the assistance of HuR via m6A modification and activates NF-κB to promote the malignant progression of IDH-wildtype glioma

Understanding the role of N6-methyladenosine (m⁶A) in tumorigenesis and stem cell maintenance is an emerging field in glioma research. However, it is necessary to study the function of m⁶A in IDH-mutation and IDH-wildtype gliomas separately. Here, we aimed to elucidate the role and mechanism of the m⁶A writer METTL3 in regulating the malignant progression of IDH-wildtype gliomas. We demonstrated that METTL3 expression is positively associated with a higher malignant grade and poorer prognosis of IDH-wildtype gliomas but not IDH-mutant gliomas. METTL3 could also promote the malignant progression of gliomas in both in vitro and in vivo models. Mechanistically, METTL3 upregulated MALAT1 expression by enhancing its stability via m⁶A modification. We further revealed that HuR was essential for METTL3-mediated MALAT1 stabilization, and upregulated MALAT1 subsequently activated NF-κB. Taken together, our findings confirmed that METTL3 promoted the malignant progression of IDH-wildtype gliomas and revealed important insight into the upstream regulatory mechanism of MALAT1 and NF-κB with a primary focus on m6A modification.