Reversible N6-methyladenosine of RNA: The regulatory mechanisms on gene expression and implications in physiology and pathology.

N6-methyladenosine (m6A) is the most abundant inner RNA modification in eukaryotes. Due to the development of RNA sequencing technology, the distribution pattern of m6A in the transcriptome has been uncovered. Dynamically, the reversible N6-methylation is mediated by two types of proteins, which are classified as "writers" and "erasers". Under the association of specific co-factors, writers show spatiotemporal N6-methyltransferase activity. Mechanically, m6A can be recognized by "reader" proteins or can directly modify RNA conformation, and it widely affects gene expression by mediating RNA stability, translation, splicing and export. m6A is involved in a series of physiology processes. Dysregulation of m6A is gradually defined as the pathogenesis of some diseases, e.g., cancer and cardiovascular disease. Therefore, a good understanding of m6A is essential for molecular biology and pathology research. In this article we systemically present an overview of the functions and mechanisms of identified m6A regulators. The discovered biological and pathological processes affected by m6A are also summarized. We hope that readers with related research interests benefit from our review.

Introduction

Strict regulation of gene expression is indispensable for development and physiological processes. Numerous protein and nucleic acid factors have been identified to affect gene expression at distinct stages of RNA processing, such as transcription, splicing and maturation, exportation and protein translation. In recent years, it has been realized that the most abundant inner chemical RNA modification, N6-methyladenosine (m6A), has novel regulatory roles in gene expression. m6A was first characterized in the 1970s in eukaryotes ranging from yeast to mammals.1, 2, 3, 4 In the following 40 years, although several components of m6A methyltransferase or demethylase have been identified, such as methyltransferase-like protein 3/14 (METTL3-14), WTAP (Wilms tumour 1-associating protein) and fat mass and obesity-associated protein (FTO), the detailed regulatory mechanisms of these proteins have not been described.5, 6, 7, 8, 9 In 2012, two independent groups assayed transcriptomic m6A modification sites by means of an m6A immunoprecipitation with next-generation sequencing (MeRIP-Seq)., According to their data, more than 12,000 N6-methylation peaks exist in over 7000 mRNAs and long non-coding RNAs (lncRNA) of human and mouse. These m6A sites share a consensus motif of “RRACH” (R = A or G; H = A, U, or C) and are enriched in the regions of the 3′ UTR, long internal exons or near stop codons. Moreover, the existence of intronic m6A sites suggests that the m6A modification might happen before pre-mRNA splicing. Most of the m6A sites are conserved between human and mouse. Intriguingly, the level of m6A modification of partial m6A sites is different in normal tissues and tumours. Moreover, silencing the methyltransferase METTL3 not only causes a significant decrease of the total m6A level but also affects gene expression and induces apoptosis of HepG2 cells. These results infer a potential association among m6A, RNA processing and cell fate. Based on the foundational impact provided by the m6A mapping approach, the upstream regulatory and downstream functional mechanisms have been gradually uncovered in the past several years: m6A is reversibly managed and maintained in dynamic homeostasis (Fig. 1). RNAs are N6-methylated by the methyltransferase complex, and the central components of the complex are called “writers”, while the major m6A demethylases are called “erasers”. N6-methylated RNAs will be recognized and will interact with a series of “reader” proteins, which are involved in RNA processing, leading to alteration of gene expression. In addition to polyadenylate-tailed mRNA and lncRNA, m6A has also been detected in microRNAs (miRNAs), circular RNAs (circRNAs) and ribosomal RNAs (rRNA) to carry out novel roles (described below). It has also been found that m6A mediates gene expression, widely affecting physiological and pathological processes. Here, we summarize the latest progress from studies of the m6A regulatory system. The significant physiological and pathological mechanisms in which m6A is involved are also presented.

Figure 1

N6-methylation is reversible. The dynamic homeostasis of m6A that regulated by writers and erasers involves in cellular physiological processes including cell proliferation, differentiation, metabolism and death.

RNAs are reversibly m6A-modified through a series of enzymatic systems

m6A is one of the few reversible inner RNA modifications. Methylation and demethylation of m6A are enzymatically regulated. N6-methyltransferation involves two methyltransferase complexes. The approximately 200-kDa complex, which consists of METTL3 and METTL14, is defined as methyltransferase A (MT-A).5, 6, 7 MT-A performs the majority of the catalytic activity. The other complex, which is approximately 800 kDa in size, is defined as methyltransferase B (MT-B) and might be the accessorial fraction required for cellular localization, RNA targeting and recruitment. N6-demethylation is mainly controlled by FTO and ALKBH5. Although the exact components and structures of the methyltransferase complexes are not fully characterized, the central “writers” or “erasers” and their roles have been identified.

METTL3: the catalytic core

It was not known how RNAs were N6-methylated until MT-A was isolated from the nucleus of HeLa cells in 1997. METTL3 was the only characterized subunit of MT-A at that time. METTL3 belongs to the S-adenosyl-l-methionine (SAM)-dependent methyltransferase (SAM MTase) superfamily, which contains a SAM-binding domain and transfers the methyl group of SAM to the adenosine of the target molecule. METTL3 is highly conserved among eukaryotes from yeast to mammals. Among all the identified subunits of the m6A writer complex, METTL3 performs fundamental N6-methylase catalytic activity. Genomic deletion of METTL3 results in the complete loss of m6A in the transcriptome. Generally, METTL3 is concentrated in the nucleus speckle and catalyses N6-methylation of RNAs, while cytoplasmic METTL3 has also been detected in several cancer cell lines.,, Cytoplasmic METTL3 functions in a methylation-independent manner: it can interact with eukaryotic translation initiation factor 3 subunit H (EIF3H) and the 3′ UTR of target mRNAs, forming a METTL3-EIF3H-mRNA loop to facilitate translation.,

METTL14: the N6-methyltransferase-inactive collaborator of METTL3

METTL14, a homologue of METTL3, is another component of the MT-A complex. These two SAM MTases form a heterodimer and co-localize in the nucleus speckle. METTL14 was first characterized as a catalytic subunit of the methyltransferase complex that methylates distinct m6A sites from METTL3. Crudely purified METTL14 showed N6-methylase activity in vitro. However, sequence analysis suggested that METTL14 only contains a disrupted SAM-binding domain and has no catalytic activity. The later structural analysis not only revealed the interaction of METTL3 and METTL14 but also revealed the fact that METTL14 lacks a SAM-binding domain. Indeed, silencing of METTL14 largely decreases the level of m6A. In addition, compared with the METTL3-METTL14 heterodimer, the catalytic activity of monomeric METTL3 was rather weaker in vitro., This finding indicates that METTL14 may act as a catalytic enhancer of METTL3. At present, METTL14 is proposed to be the stabilizer, activator and RNA binding scaffold of METTL3.,

WTAP: the crucial partner for localization of the catalytic core

WTAP was initially identified as a protein associated with Wilms' tumour suppressor WT1 by a yeast two-hybrid assay. Because of its ubiquitous expression and nuclear localization, WTAP was also previously considered to be a house-keeping or RNA-processing related protein. Following studies demonstrated that WTAP regulates the alternative splicing of pre-mRNA and the stability of mRNA., WTAP also suppresses the expression of WT1-target genes by interacting with WT1. In 2014, WTAP was found to be involved in m6A regulation. Pull-down assays showed that WTAP binds both METTL3 and METTL14. Knockdown of WTAP causes a significant decrease in the level of m6A. As WTAP lacks N6-methyltransferase activity and shares consensus target m6A sites with METTL3, it is suggested that WTAP helps METTL3-METTL14 localize to nucleus speckles and recruit other N6-methyltransferase components. Silencing of WTAP disrupts the nucleus speckle-localization of METTL3-METTL14.

Other writers with accessory functions

In recent years, other associated writers have gradually been discovered: In 2013, Horiuchi et al identified a series of WTAP-interacting proteins. In addition to METTL3 (described as MTA-70), protein virilizer homologue (VIRMA, also known as KIAA1429), KIAA0853 (also known as Zinc finger CCCH domain-containing protein 13, ZC3H13), RNA-binding protein 15 (RBM15), E3 ubiquitin-protein ligase HAKAI (also known as Casitas B-lineage lymphoma-transforming sequence-like protein 1, Cbll1) and two arginine/serine-rich domain-containing proteins, BCLAF1 and THRAP3 are included in the WTAP-binding protein complex. These WTAP-binding proteins are related with the alternative splicing of pre-mRNA. This study preliminarily screened out possible candidates associated with the N6-methyltransferase complex. Further research indicated that VIRMA, ZC3H13, RBM15 and HAKAI are involved in N6-methylation: VIRMA may be involved in the selective N6-methylation of sites in 3′ UTRs and located near stop codons in mRNAs. ZC3H13 acts as a connector between WTAP and RBM15. ZC3H13 is also essential for the nuclear localization of other writers. Silencing of ZC3H13 leads to cytoplasmic localization of VIRMA and HAKAI. RBM15 interacts with WTAP in a ZC3H13 dependent manner. RBM15 functions as a connecter between catalytic writers and readers. As an RNA binding protein (RBP), RBM15 recruits specific target RNAs to catalytic centres and guides the reader to selectively recognize m6A sites. HAKAI ordinarily functions as an E3 ubiquitin ligase. The detailed roles of HAKAI in N6-methylation are still unclear.

Catalytic N6-methylases with novel targets and organization

With the exception of METTL3, there are also N6-methylases unrelated to MT complexes for specific types or structures of RNAs. METTL16, a homologue of METTL3, can N6-methylate a conserved hairpin of MAT2A mRNA and U6 snRNA with stimulation of SAM. Further exploration revealed that the methylation sites of METTL16 are enriched in introns and are not conserved, which is distinct from the METTL3-modulated pattern. The detailed mechanisms of METTL16 still need to be explored. Zinc finger CCHC domain-containing protein 4 (ZCCHC4) has been reported to be an N6-methylase for 28S rRNA. Loss of ZCCHC4 globally inhibits protein translation. The proteomic association and roles involved in N6-methylation of ZCCHC4 are unclear.

Mammalian m6A erasers

To date, FTO and ALKBH5 are the only two characterized N6-demethylases. Both of these N6-demethylases are members of the alpha-ketoglutarate-dependent dioxygenase AlkB-like superfamily. These AlkB-like dioxygenases have a catalytic domain that is homologous to that in Fe-2OG dioxygenases and catalyses the hydroxylation of the target methyl by using 2-oxoglutarate and oxygen molecular as substrates. The hydroxymethyl is unstable when it is connected to a nitrogen atom, resulting in demethylation.,

It has been demonstrated that FTO plays multiple catalytic roles and has multiple substrates. FTO was initially discovered to catalyse the demethylation of 3-methylthymidine (3mT) in single-stranded DNA and 3-methyluracil (3mU) in single-stranded RNA., In 2011, the N6-demethylase activity of FTO was characterized in vitro. FTO was the first identifiedm 6A eraser. A recent study discovered that FTO can demethylase the N6,2-O-dimethyladenosine (m6Am) in the 5′ cap of mRNAs. The crosstalk among these types of methylation deserves to be explored. Immunostaining indicated that FTO is concentrated in nucleus speckles, suggesting that FTO might demethylate nuclear m6A-modified RNAs. Cytoplasmic FTO and its high N6-demethylase activity were also observed in certain types of cells, such as leukaemia cells. An in vivo study showed that FTO might not mediate N6-demethylation in the early embryonic stage. The m6A level in FTO-knockout mouse embryos was not significantly affected. These data suggest that FTO may spatiotemporally catalyse N6-demethylation.

The N6-demethylase activity of ALKBH5 was discovered immediately after FTO in 2013. Modifying the expression of ALKBH5 affects the level of m6A in certain types of cells. In addition, ALKBH5 deficiency in mice is not lethal and only increases the apoptotic rate of testicular cells and leads to defects in spermatogenesis. These data suggest that ALKBH5 selectively demethylates m6A, similar to FTO. The demethylation process mediated by ALKBH5 may be necessary for RNA processing, which regulates testis development and male fertility. On the other hand, unlike FTO, ALKBH5 has specific nuclear localization and N6-demethylase activity., Loss of function analysis showed that ALKBH5 mediates RNA export in HeLa cells and RNA splicing in sperm cells in an m6A-demethylase activity-dependent manner, perhaps due to the interaction between ALKBH5 and splicing factors.

Selective N6-methylaiton in the transcriptome

Although the “RRACH” motif appears frequently in transcriptomes, the distribution of m6A sites is enriched in specific regions, including long internal exons, the 3′-UTR and regions near stop codons of mRNAs. Intronic N6-methylation also exists in pre-mRNA., Moreover, quite a large portion of m6A sites exhibit spatiotemporally specific patterns. Since the catalytic core, METTL3/14, is ubiquitously expressed and has not been demonstrated to independently recognize specific regions, selective N6-methylation is believed to be carried out by some associated components. In this section we would like to discuss the demonstrated and potential mechanisms involved in selective N6-methylation.

VIRMA or RNA polymerase II regulates N6-methylation in the 3′-UTR and near stop codons

Evidence has suggested that some accessary writers function as guides to recruit the N6-methyltransfase complex to regions in the 3′ UTR and near stop codons. VIRMA is the largest characterized component of the N6-methyltransferase complex and is indicated to be the scaffold that bridges the N6-methyltransferase complex and target RNAs. MeRIP-Seq data in HeLa cells indicate that knockdown or depletion of VIRMA significantly disrupts the distribution of m6A sites in the 3′-UTR and near stop codons without altering the total m6A level, revealing that VIRMA plays critical roles in 3′-UTR- and stop codon-specific targeting. However, the detailed mechanism of the selective targeting mediated by VIMRA is unclear. VIRMA-associated cleavage and polyadenylation specificity factor subunit 5 and 6 (CPSF5 and CPSF6), which are the components of the CFIm complex, are able to bind to “UGUA” elements in the 3′-UTR of pre-mRNAs. It has been observed that VIRMA and CPSF5 share a substantial overlap of targets, which implies that the CFIm complex might be involved in selective N6-methylation. Nevertheless, silencing of CPSF5 and CPSF6 only increase the density of N6-methylation in the 3′ UTR. Perhaps other subunits of the CFIm complex mediate this process.

Another possible mechanism is that selective N6-methylation in the 3′ UTR and near stop codons might be co-transcriptionally regulated. RNA polymerase II (Pol II) has been reported to interact with the N6-methyltransfase complex. A Pol II retardation-induced transcription pause increases the recruitment of m6A writers and the level of m6A. It could be suggested that the localization of transcription termination signals near the corresponding region of the 3′ UTR and stop codons might result in the accumulation of m6A sites.

Potential transcriptional and post-transcriptional mechanisms of N6-methylation in long internal exons

Exonic N6-methylation might be correlated with some exon-associated factors. Generally, the corresponding positions in the chromatin of long exons show higher histone methylation levels than other regions of mRNA. It is intriguing that histone methylation of H3K36me3 was reported to be related with N6-methylation. METTL14 can interact with H3k36me3 in chromatin. Knockdown of histone-lysine N-methyltransferase SETD2 decreases the m6A level in mRNA by approximately 40%. Therefore, it was proposed that H3k36me3 guides N6-methylation in the corresponding positions of nascent RNAs, which leads to a high N6-methylation density in long exons. However, the correlation between H3K36me3 sites and N6-methylaiton sites has not been confirmed.

It is known that various RBPs bind specific motifs in pre-mRNAs, mediating RNA processing post-transcriptionally. Long internal exons have a high density of the RBPs binding motif. Before discussing the association between RBPs and long internal exon-specific N6-methylation, we must make clear whether N6-methylation also occurs in the post-transcriptional period. Notably, a study of m6A in circRNAs performed by Zhou et al may provide some clues. They observed that among a total of 1191 detected circRNAs in human embryonic stem cells, 304 circRNAs were generated from genes that encoded a mRNA without the m6A signal. The rest of the 887 circRNAs were generated from genes that encoded N6-methylated mRNAs, and the m6A sites of 392 circRNAs were distributed in different exons from the exons in the cognate N6-methylated mRNAs, displaying somewhat distinct distribution patterns. circRNAs are generated from pre-mRNA via back-splicing post-transcriptionally, which results in the circularization of exons and pairing of flanking introns, forming a lariat structure.41, 42 The distinct m6A distribution pattern in mRNA and circRNAs suggests that RNAs are also N6-methylated post-transcriptionally. Therefore, the N6-methyltranferase-associated RBPs were possibly able to guide the selective N6-methylation. For instance, RBM15, an accessary m6A writer mentioned above, is also an RBP. RBM15 has been demonstrated to be involved in N6-methylation of lncRNA XIST. It could be observed from the CLIP sequencing data that an abundance of RBM15 binding sites were adjacent to N6-methylation sites in exons. Knockdown of RBM15 reduced the level of m6A in mRNA, suggesting that RBM15 may recruit the N6-methyltransfase complex to its targets.

Conformational regulation of RNA G-quadruplexes in intronic N6-methylation

G-quadruplexes (G4s) are a type of non-canonical secondary structure in DNA or RNA that form by guanine-rich sequences. These sequences are defined as potential G-quadruplex-forming sequences (PQSs) and are characterized by 4 guanine clusters separated by loops. It has been reported that G4s can modulate DNA methylation. Recently, the potential correlation between m6A and G4s also emerged from the epitranscriptomic analysis performed by Burrow's group. Initially, several instances of colocalization of m6A and PQSs in virus were explored. At the time of the writing of this review, they further discovered that m6A in intron splice sites of human pre-mRNAs were frequently localized in PQS regions. These observations suggest that selective N6-methylation might partially result from G4s. Since G4s sequences have been demonstrated to restrict DNA methylation by binding the methyltransferase DNMT1, it is likely that intronic RNA G4s are able to modulate local m6A by recruiting m6A writers in a similar way. Nevertheless, it was observed that non-intronic N6-methylation sites did not show a strong correlation with G4s. Intronic m6A prefers to be added at SAG consensus motifs, which may be indicate that these are the limiting elements involved in the m6A modification. It is possible that some G4s interact with intermediates that can recognize the SAG motifs and act as guides for the m6A writers. However, this hypothesis should be further confirmed by mutational analysis. Additionally, the colocalization of m6A and G4s need to be widely verified in other species.

Transcription factors and histone methylation are associated with transcript-specific N6-methylation

It has been observed that the m6A peaks are always concentrated in partial transcripts from accumulated MeRIP-Seq data, which implies that the N6-methyltransfases are recruited to specific transcripts during transcription or post-transcription periods. Histone methylations and transcription factors are two types of essential transcriptional regulators. In addition to histone methylation, as discussed above, transcription factors are also indicated to be involved in selective N6-methylation and to guide N6-methylation co-transcriptionally. For instance, the SMAD family member 2/3 bound METTL3/14 and co-transcriptionally recruits METTL3/14 to its target transcripts. METTL3 interacts with selective gene promotors by associating with the CCAAT enhancer binding protein zeta. At the post-transcriptional level, RBP is one of the potential regulators of transcript-specific N6-methylation. Other related factors remain to be detected.

Chemical modification affects the activity of N6-methyltransfase

Most m6A sites exhibit spatiotemporally specific patterns, and the level of m6A also dynamically fluctuates during different developmental or biological stages. As m6A is enzymatically regulated, the activity of the N6-methyltransfase is closely related to the diversity of the level of m6A. The catalytic activity of METTL3/14 has been demonstrated to be regulated by a variety of chemical modifications. For instance, small ubiquitin-related modifier 1 (SUMO1) inhibits the N6-methyltransferase activity of METTL3 by SUMOylating the lysine residues of METTL3 without affecting its protein level or its association with other N6-methyltransferase components. This SUMOylation can be decreased by sentrin-specific protease 1. The phosphorylation of residue Ser-399 of METTL14 is crucial for its interaction with METTL3. Mechanically, the phosphorylated Ser-399 in METTL14 forms a salt bridge with ‘Arg-471’ of METTL3. Dephosphorylation of Ser-399 disrupts the interaction of METTL3 and METTL14, resulting in the inactivation of N6-methyltransfase. Proteomic analysis also revealed phosphorylation sites in WTAP. However, it is unclear whether the phosphorylation in WTAP is related to localization and interaction with the N6-methyltransfase complex.51, 52, 53 The mechanisms affecting the level of m6A writers have not been clearly reported. Although the expression of writers is dysregulated in some diseases, the details of this dysregulation are unknown.

The cooperation of a series of factors dynamically modulate N6-methylation and lead to distinct N6-methylation patterns in various types of cells and stages. More underlying factors and mechanisms remain to be explored and verified. In addition, it is unclear whether an m6A eraser and N6-demethylation are involved in the formation of the m6A distribution pattern. Analysis of the associations of erasers might be able to provide answers.

Recognition of m6A and its effects on RNA processing and gene expression

m6A controls the fate of RNAs. Aberrant expression levels of writers or erasers result in the dysregulation of RNA degradation, exportation, translation and more. In most cases, N6-methylation does not affect RNA features directly (Fig. 2). Readers are able to recognize m6A and regulate the stability, exportation and translation of methylated RNAs. The YTH (YT521-B homology) domain containing protein family (YTHDF) is a group of evolutionarily conserved m6A readers. The YTH domain has been suggested to bind RNA in an m6A-dependent manner.54, 55, 56 The family members YTHDF1-3 contain long, low-complexity regions enriched with Pro-Xn-Gly motifs. The association of N6-methylated RNAs promotes the phase separation of YTHDFs. Conversely, YTHDFs lead N6-methylated RNAs to diverse compartments in the cytosol, which then undergo specific processes. YTHDF1 is transported to stress granules by the guidance of the bound N6-methylated mRNAs. In stress granules, YTHDF1 recruits and interacts with EIF3 to form a loop structure composed of mRNAs and translational factors. This association further enhances ribosome loading and translation initiation of target mRNAs. Distinct from YTHDF1, YTHDF2 affects the stability of target RNAs. The interaction of N6-methylated RNAs with YTHDF2 promotes the YTHDF2-RNA complex to localize in processing bodies. The CCR4-NOT deadenylation complex and decapping complex are recruited by YTHDF2 and degenerate N6-methylated RNAs., The third member of the YTHDF family, YTHDF3, has been reported to serve as the functional partner of the other YTHDFs in HeLa cells. It was revealed that YTHDFs share a considerable amount of target transcripts and interact with each other, while YTHDF3 was the earliest member to recognize and bind N6-methylated RNAs. YTHDF3 associates with the translation initiation complex, facilitates the activity of YTHDF1 and subsequently increases the efficiency of YTHDF2, indicating that YTHDF3 seems to be the deliverer of N6-methylated RNAs. Insufficient YTHDF3 decreases the m6A affinity of YTHDF1 and YTHDF2., In addition, Yang et al also reported the novel finding that YTHDF3 can interact with eukaryotic translation initiation factor 4 gamma 2 (EIF4G2) and recruit translation initiation cooperators to the target m6A sites in circular RNAs. This process drives the 5′ cap independent translation of circular RNAs. These functional studies of YTHDF3 not only revealed the translation regulation manner of YTHDF3 but also indicated that the members of the YTHDF family are interrelated counterparts of the m6A metabolic machinery. However, structural analysis might still be the most efficient means to detect the underlying domains in YTHDFs that are crucial for the spatiotemporal distinction of phase separation and target specify.

Figure 2

Schematic of N6-methylation/demethylation and m6A-regulated RNA processing. RNAs are N6-methylated by the writers and reversely N6-demethlyated by the erasers. The processes of methylation/demethylation may be guided by specific transcriptional factors, RBPs and histone modifications. N6-methylated RNAs are further recognized and processed by a serious of readers and associated factors.

YTH domain-containing protein 1 (YTHDC1) is another cluster of readers that are localized in the nucleus. YTHDC1 is functionally related to RNA maturation. When pre-mRNA is transcribed, YTHDC1 promotes exon-inclusion by recruiting serine/arginine-rich splicing factor 3 (SRSF3) to the target pre-mRNA and displacing the interacted SRSF10. Furthermore, the interaction of YTHDC1, SRSF3 and nuclear RNA export factor 1 (NXF1) promotes nuclear export of N6-methylated mRNAs. Chen et al demonstrated that YTHDC1 facilitates the export of the N6-methylated circular RNA circNSUN2. This research direction was significant for completing our understanding of the splicing and export mechanisms of circular RNA. Further explorations of the associated splicing and export factors were expected. YTHDC1 was also found to interact with the N6-methylated nuclear lncRNA XIST and to enhance XIST-mediated X chromosome silencing. In addition to XIST, other N6-methylated nuclear lncRNAs, such as NEAT1 and MALAT1, bind to YTHDC1. However, the effect of this binding is unknown. Recently, YTHDC1 was detected to be involved in m6A-mediated transcription regulation. Nuclear chromosome-associated regulatory RNAs (carRNAs), including promoter-associated RNA (paRNA), enhancer RNA (eRNA) and RNA transcribed from transposable elements (repeats RNA), can recruit CBP/EP300 and YY1 to promote chromatin opening and transcription. YTHDC1 interacts with N6-methylated carRNAs and degrades them via association with nuclear exosome targeting complex (NEXT), resulting in selective inhibition of transcription.

YTHDC2 is the largest YTH protein, and it contains multiple helicase domains and two ankyrin repeats. Correspondingly, YTHDC2 plays a translation regulatory role, similar to 3′–5′ RNA helicases. YTHDC2 has also been indicated to affect mRNA stability by associating with XRN1 exoribonuclease. However, the effect of YTHDC2 is relatively weak because of its limitations in m6A binding. YTHDC2 has been shown to have a distinct m6A target pattern and lower affinity for N6-methylated RNAs than other YTH readers.,, Moreover, YTHDC2-deficiency only causes a spermatogenesis defect., A CLIP study suggested that YTHDC2 does not interact with N6-methylated RNAs in HEK293 cells., It is possible that YTHDC2 mediates m6A in a novel and tissue-specific way.

m6A are believed to be related to the maturation of miRNAs. m6A in pri-miRNAs improves DGCR8 binding and the level of mature miRNAs. The promotion induced by m6A occurs in a global and non-cell-type specific manner. Since DGCR8 cannot specifically recognize m6A, ubiquitously expressed nuclear m6A readers may act as cooperators for the recruitment of N6-methylated pri-miRNA, such as YTHDC1.

Therefore, we can summarize that m6A functions by methylation, recognition and processing of RNA in general. Following writer-regulated methylation, readers further specifically deal with N6-methylated RNAs. The distinct structures of readers determine their localization, target points and associated cofactors. In addition to RNA processing, the associated cofactors might also reversibly affect readers by selectively binding m6A. As a result, the fate of RNAs is organically regulated.

In a few cases, an RNA structure modified by m6A would function directly by altering the interactions of RNAs and proteins. Several members of the heterogeneous nuclear ribonucleoprotein (HNRNP) family, including HNRNPC, HNRNPA2B1 and HNRNPG, which are involved in RNA processing, are affected by m6A in this way.74, 75, 76 For example, HNRNPC binds the single-stranded poly-U motif in pre-mRNA and regulates pre-mRNA processing. Liu et al reported that m6A destabilized the U-A pairing that could expose the single-stranded U-enriched region in pre-mRNA from UA duplexes and facilitate the binding of HNRNPC.

Insulin-like growth factor 2 mRNA-binding proteins 1-3 (IGF2BP 1-3) have also been reported to bind and stabilize N6-methylated mRNAs. However, it seems that IGF2BPs cannot specifically recognize m6A. It is possible that IGF2BPs indirectly interact with N6-methylated mRNAs depending on other cofactors.

Biological functions and dysregulation of m6A in diseases

Abundant m6A controls gene expression throughout the lifespan of an organism. It has been well proven that the highly regulated stages of m6A modification, recognition and processing are indispensable for diverse developmental and physiological activities. Conversely, once writers, erasers or readers become disorganized, the induced unbalanced gene expression causes disease or other biological defects (Table 1). In this part, we summarize representative studies in biology and pathology.

Table 1

m6A-involved biological and pathological processes.

Physiological and pathological processes	Critical m6A regulators	Representative target genes/transcripts	Pathogeny	References
Embryonic development	METTL3, YTHDF2, YTHDC1	carRNAs	–	12,18,69,79
Haematopoietic stem transition	METTL3, YTHDF2	notch1a	–	80
Neural stem cell maintenance	METTL3/14	pax6, neurogenin-2, cbp, p300	–	81,82
Myelination of axons	METTL14	neurofascin 155	–	83
Male fertility	ALKBH5, YTHDC2	Transcriptomic research	–	37,72,85
Female fertility	YTHDF2, YTHDC2	smc3, cep76	–	71,72,86
Dopaminergic midbrain circuitry	FTO	drd3, kcnj6, grin1	–	89
Learning and memory	YTHDF1	Transcriptomic research	–	90
T cell homeostasis	METTL3	socs1, socs3, cish	–	91
Anti-virus	ALKBH5	OGDH	–	92
Hepatocellular carcinoma	METTL3, YTHDF2	SOC2	Upregulation of METTL3	93
Endometrial cancer	METTL3/14, YTHDF1/2	PHLPP2, PRR5, PRR5L, mTOR	METTL3 deficiency, METTL14 mutation	94
Glioblastoma	ALKBH5	FOXM1	Upregulation of ALKBH5	95
Acute myeloid leukaemia	METTL14, FTO	MYB/MYC, ABS2, RARA	Upregulation of METTL14 and FTO	36,96
Ischaemic cardiac injury	METTL3	TFEB	Upregulation of METTL3	97
Cardiac hypertrophy	METTL3	Transcriptomic research	Upregulation of METTL3	98
Heart failure	ALKBH5	serca2a	Downregulation of ALKBH5	99

m6A is essential for selective gene expression in development and cell differentiation

Cell differentiation is a milestone of development. m6A was initially considered to be a significant epigenetic regulator of embryonic stem cell (ESC) differentiation. By means of MeRIP-Seq, high levels of m6A were detected in the transcriptomes of mammalian ESCs that were especially enriched in the transcripts of genes critical for pluripotency. Decreasing the level of m6A via knockout of METTL3 led to the abolished differentiation of ESCs and embryonic lethality in mice. Pluripotency genes such as NANOG and oct4 showed excessive levels in the status of METTL3 deficiency.,, Later, m6A was indicated to be the signal for YTHDF2 to degrade specific mRNAs in zebrafish embryo development. This result corresponds with the previously observed abnormally sustained mRNAs induced by METTL3-KO. In addition to its YTHDF2-dependent “turn-off” role of gene expression, m6A fulfilled YTHDC1-dependent transcriptional inhibition in the earlier stage of embryonic development (see the details of the mechanism in the previous section). These data demonstrate that m6A controls selective gene expression in embryonic development at both the transcriptional and post-transcriptional levels.

In addition, m6A is involved in secondary cell differentiation. For instance, m6A mediates the transition from haemogenic endothelial cells to haematopoietic stem and progenitor cells in zebrafish. During haematopoietic development, notch1a mRNAs are N6-methylated by METTL3 and then decayed by YTHDF2. The inhibition of notch signalling facilitates the generation of haematopoietic stem and progenitor cells. m6A also maintains the self-renewal and prevented premature differentiation of neural stem cells. During neurogenesis, the mRNAs of genes related to neuronal differentiation and development, such as pax6 and neurogenin-2 and histone H3K27 acetyltransferases, including cbp and p300, are N6-methylated. The expression of differentiation-related genes is inhibited by m6A-directed mRNA decay and a decrease of H3K27-ac., At the subsequent stage of the myelination of axons, m6A mediates the splicing of global RNAs in the oligodendrocyte-lineage, affecting the expression of genes that are essential for myelination, such as neurofascin 155. Knock-out of METTL14 retards the specialization of oligodendrocyte precursor cells to mature oligodendrocytes, indicating that m6A acts as a “switch” for gene expression during development and cell differentiation.

m6A regulates germ cell maturity and fertility

In 2013, Zheng et al discovered that ALKBH5 modulates spermatogenesis. Among eight different organs, including the heart, brain, gonadal fat pads, liver, kidney, spleen, lung and testis, the testis showed the highest expression level of ALKBH5 mRNAs. ALKBH5-KO mice were able to survive to adulthood. However, male mice were infertile and had a significantly undersized testis, malformed spermatozoa, maturation arrest and apoptotic spermatocytes. Although differential expression of genes related to spermatogenesis and apoptosis was detected in ALKBH5-deficency testis, the detailed mechanisms of ALKBH5 in spermatogenesis were unclear at that time. In regards to the transcriptome, mRNA shortening at the period of spermiogenesis is necessary for mRNA metabolism efficiency. Tang et al put forward the hypothesis that ALKBH5 is involved in this process. They observed that m6A was enriched in the long 3′UTR of mRNAs during spermiogenesis. Via demethylation of ALKBH5, these mRNAs with a long 3′UTR avoid decay and are spliced to short forms.

Although m6A mediates female fertility, the regulatory mechanism is distinct from that found in spermatogenesis. During the period of MII oocyte production, YTHDF2 plays the critical role of eliminating N6-methylated maternal RNAs, which is an important process for oocyte maturation and early zygotic development. Loss of YTHDF2 would result in female-specific infertility.

The m6A reader YTHDC2 is a special mediator of spermatogenesis and oogenesis., YTHDC2 mRNAs are highly expressed in the testis and increased in sperm cell meiosis. Knockout of YTHDC2 in mice impairs testis and ovary development. Mechanically, YTHDC2 improves the translation of its target mRNAs, including spermatogenesis-related smc3 and centriole reduplication-related Cep76, during meiosis.,

m6A affects neural activities

The mammalian nervous system is an accurately organized control system that regulates behaviour, metabolism, and other functions. Various neural actions are further regulated by complex molecular mechanisms that are still under exploration. As the most abundant RNA modification, m6A also exists in the transcriptome of neural cells. The detailed roles of m6A in neural activity are being gradually revealed. Defects of the m6A eraser FTO was observed to be associated with abnormal phenotypes, such as attention deficit disorder, growth retardation, aberrant energy expenditure and locomotor activity., Hess et al revealed that FTO mediates the conduction of dopaminergic signalling, which controls a serious of neural activities. Demethylase-inactivation of FTO impairs D2-like receptor-regulated neural and behavioural reactions under the stimulation of cocaine. MeRIP-Seq indicated that knockout of FTO increases m6A levels in mRNAs in the mouse midbrain and striatum. The mRNA levels of genes related to dopaminergic signalling, including drd3, kcnj6 and grin1, were upregulated; but conversely their protein levels were decreased. Shi et al showed that m6A affects learning and memory in a YTHDF1-dependent manner. By neuronal stimulation, YTHDF1 facilitates protein translation of target m6A mRNAs to promote synaptic transmission in the mouse hippocampus. Deficiency of YTHDF1 impairs learning and memory.

m6A involves in immunity

m6A has been indicated to modulate gene expression during the maturation of T cells. JAK1/3-STAT5 signalling is critical for T cell differentiation and proliferation, which is inhibited by SOCS family proteins. Through stimulation by IL-7, the mRNAs of SOCS family genes, including socs1, socs3 and cish, were N6-methylated and degraded. JAK1/3-STAT5 signalling was disinhibited and activated. METTL3-deficiency disrupts native T cell homeostasis. An m6A-dependent antiviral mechanism was also revealed. The ALKBH5 protein is inactivated by altered methylation in mouse primary peritoneal macrophages when suffering from infection of vesicular stomatitis virus. Inhibition of ALKBH5 increases the N6-methylation level of α-ketoglutarate dehydrogenase (OGDH) mRNAs. OGDH catalyses the synthesis of itaconate, which is required for viral replication. m6A destabilizes OGDH mRNAs and decreases the level of the OGDH protein. Viral replication is, hence, retarded by itaconate deficiency.

Aberrant m6A is an oncogenic event

The molecular essence of cancer is characterized by the induction of oncogenes and inhibition of tumour suppressor genes. Accumulated studies have demonstrated that aberrant m6A regulation is one of the pathogeneses of cancers. In some situations, abnormal N6-methylase or N6-demethylase activities dysregulate cancer-related gene expression. For instance, METTL3 is significantly increased in human hepatocellular carcinoma. Increased METTL3 results in high N6-methylation and YTHDF2-controlled degradation of tumour suppressor SOCS2 mRNA. Downregulation of SOCS2 promotes the proliferation and invasion of tumour cells. Mutation of METTL14 and METTL3-insufficiency have been detected in endometrial cancer. The decrease of N6-methylation activates tumorigenic AKT signalling through two ways: First, dephosphorylation of the S473 residue of AKT(S473) by phosphatase PHLPP2 is decreased. Low levels of m6A on PHLPP2 mRNA inhibits the promotion of translation managed by YTHDF1. On the contrary, the phosphorylation of AKT by serine/threonine-protein kinase complex mTORC2 is increased. YTHDF2 targets the N6-methylated mRNAs of mTORC2 components, including PRR5, PRR5L and mTOR. The decrease of N6-methylase activity inhibits their degradation. Therefore, phosphorylation of AKT is facilitated and tumorigenic signalling is activated. In glioblastoma stem-like cells, ALKBH5 is highly expressed. The transcription factor FOXM1 is beneficial for the proliferation of glioblastoma stem-like cell and tumorigenicity. FOXM1 mRNA is targeted and demethylated by ALKBH5. The interaction of ALKBH5 and FOXM1 mRNA relies on lncRNA FOXM1-AS. N6-demethylated FOXM1 mRNA recruits HuR to increase its stability and protein expression. Intriguingly, upregulation of METTL14 and FTO have both been detected in acute myeloid leukaemia (AML). Although the global m6A level is increased in the transcriptome of AML cells, METTL14 and FTO have both been shown to be oncogenic factors., Weng et al revealed that METTL14 increases the N6-methylation, stability and translation efficiency of MYB and MYC mRNAs. MYB/MYC signalling inhibits myeloid differentiation and promotes the proliferation of AML cells. Li et al demonstrated that FTO demethylates the mRNAs of ASB2 and RARA, which regulate myeloid differentiation. This process causes transcripts to be unstable and inhibits protein translation of ASB2 and RARA. Thus, FTO aggravates AML progress. Both METTL14 and FTO are negatively regulated by myeloid differentiation agents. SPI1 has been found to be a transcription repressor of METTL14. The mRNA level of FTO is decreased by All-trans retinoic acid. However, their cooperating m6A readers have not been identified. Data indicate that YTHDFs are not responsible for the mRNA stability or the translation of the targets of FTO and METTL14 in AML.

Excessive m6A is correlated with heart disease

Heart disease is the most lethal disease worldwide. The pathogenesis of heart disease still has not been fully understood. From the 2010s, studies of the post-translational regulation mechanisms of heart disease began to emerge. Certain cardiac non-coding RNAs have been identified and considered to be therapeutic targets. In 2018, the post-transcriptional roles of m6A in a serious cardiac pathology were revealed.97, 98, 99 It was discovered that elevated N6-methylase activity aggravates ischaemic cardiac injury. When cardiomyocytes or the myocardium was challenged with oxidative stress, such as hypoxia/reoxygenation (H/R) and ischaemia/reperfusion (I/R), the mRNA and protein levels of METTL3 were significantly induced. The increased METTL3 targeted pre-mRNAs of the transcription factor EB (TFEB), which is critical for the transcription of lysosomal genes and the regulation of autophagy. Then, the N6-methylated pre-mRNAs of TFEB were bound and degraded by HNRNPD, resulting in the inhibition of cardiomyocyte autophagy and increasing apoptosis. This m6A-related degradation of TFEB pre-mRNA was negatively regulated by ALKBH5. Reversely, TFEB could promote the transcription of ALKBH and suppress METTL3 at the post-transcriptional level, forming a feedback loop. Mature mammalian cardiomyocytes are terminally differentiated cells with negligible proliferation. Cardiomyocytes that die from I/R are replaced by fibroblasts, which lack the capacity to contract. Meanwhile cardiac hypertrophy is triggered to adapt the haemodynamic change. Pressure overload and energy deficiency of cardiomyocytes eventually lead to cell death and heart failure. During myocardial remodelling, m6A upregulation has been indicated to play important roles. Dorn et al observed that m6A levels on polyA+ RNAs were elevated in FBS-induced hypertrophic neonatal rat cardiomyocytes. Overexpression of METTL3 in neonatal cardiomyocytes increased the cell surface area. On the contrary, silencing of METTL3 resulted in a smaller cell size. In vivo experiments showed that both a cardiac-specific transgene (cTG) and knockout (cKO) of METTL3 did not significantly affect heart development or cardiac function in juvenile mice. However, 8-month-old METTL3 cTG mice exhibited cardiac hypertrophy, while 8-month-old METTL3 cKO mice suffered heart failure. These results suggest that METTL3 is involved in the stress adaption of heart. However, the details of its mechanisms are unclear. In addition, m6A mapping data and expression profiles of m6A regulators during the stages of mouse embryonic heart development were not shown in this research. These are crucial for rigorously revealing the roles of m6A in heart development. Mathiyalagan et al found that the transcriptomic m6A level was remarkably elevated in ischemia-induced mammalian heart failure. Decreased FTO mRNA and protein levels were detected throughout the process of heart failure. MeRIP-Seq indicated that the mRNAs of some Ca2+ homeostasis and contraction-regulation genes, such as serca2a, were hypermethylated, leading to degradation of mRNA and decrease in translation. Decreased expression of serca2a resulted in dysregulation of myocardial remodelling. Replenishment of FTO could demethylate serca2a mRNA and rescue the failing heart. In vitro experiments showed that hypoxia could inhibit the expression of FTO mRNA and protein. In addition, METTL14 has been observed to have elevated protein levels without an mRNA increase in heart failure, whether METTL14 involved in the process is unknown.

Perspectives

Since m6A was globally characterized in transcriptomic by RIP-Seq in 2012,, numerous studies on m6A have been performed. It is exciting that the mechanisms of methylation/demethylation, recognition and RNA processing have basically been revealed. Increasing pathological and physiological events have been demonstrated to be regulated by m6A. These research achievements have enriched medical theory and deepened our understanding of disease progression.

According to present studies, we can conclude that N6-methylation and demethylation act in a tissue- and developmental stage-specific pattern. For instance, N6-methylation is generally predominant in embryonic development,,,, and deficiency of m6A erasers does not lead significant embryonic defects., However, in some tissues or types of cells, such as the brain, testis, ovary and heart, N6-demethylation is physiologically indispensable.,,,,, To spatiotemporally understand N6-methylation and demethylation, several uncertain points about m6A still remain. First, the details of the associated factors of m6A writers/erasers are needed. Although the writers and erasers possess N6-methylation/demethylation activities, they have no potency to distinguish specific targeting motifs among thousands of transcripts. The association of certain transcription factors, histone modifications and RBPs have been detected, which suggests these proteins might be involved in target selection.,,,,, The details of the proteomic interactions and analysis of the target sites of the writers/erasers in different types of cells and development stages are necessary. Second, information regarding the expression regulation and modification of the writers/erasers is needed. Dynamic equilibrium between N6-methylation and demethylation is critical, as aberrant levels of writers or erasers would cause various defects or disease. It was discovered that writers and erasers are regulated at the mRNA and protein levels.93, 94, 95, 96, 97, However, most of the upstream modulators remain to be explored. In addition, distinct subcellular localization of writers and erasers exist in certain types of cells.,, Cytoplasmic localization might result in some novel functions. Perhaps protein modifications affect the subcellular localization of writers and erasers. Next, the mechanisms of differential recognition of N6-methylated transcripts need to be ascertained. Commonly, readers determine the fate of N6-methylated transcripts. Structure distinctions, cellular localization and external stress may affect the selective targeting of readers.,56, 57, 58

Aberrant m6A levels in some types of cancers might be applied to diagnosis. Combining immunohistochemistry with histological analysis of the level of m6A would facilitate more accurate diagnoses. In addition, the activity of erasers has been found to be regulated by various reagents. For example, the activity of FTO is activated by ascorbate and inhibited by N-oxalylglycine, fumarate, succinate, fluorescein derivatives, meclofenamic acid and R-2-hydroxyglutarate.,100, 101, 102 Both FTO and ALKBH5 need Fe2+ to act as cofactor., R-2-hydroxyglutarate has been especially indicated to inhibit leukaemia by affecting FTO/m6A/MYC/CEBPA signalling. The therapeutic potency of the other activators and inhibitors remains to be verified.

In conclusion, developing technology and deepening research have led us to realize the significance of m6A. Dynamic and integral mapping of m6A sites and the proteomic interaction involving each m6A regulator is expected to uncover the panorama of m6A.

Conflict of interest

The authors have declared that no competing interest exists.