Mol. Cells 2022; 45(7): 435-443
Published online June 24, 2022
https://doi.org/10.14348/molcells.2022.0017
© The Korean Society for Molecular and Cellular Biology
Correspondence to : ginalee@uci.edu
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/.
In response to environmental changes, signaling pathways rewire gene expression programs through transcription factors. Epigenetic modification of the transcribed RNA can be another layer of gene expression regulation. N6-adenosine methylation (m6A) is one of the most common modifications on mRNA. It is a reversible chemical mark catalyzed by the enzymes that deposit and remove methyl groups. m6A recruits effector proteins that determine the fate of mRNAs through changes in splicing, cellular localization, stability, and translation efficiency. Emerging evidence shows that key signal transduction pathways including TGFβ (transforming growth factor-β), ERK (extracellular signal-regulated kinase), and mTORC1 (mechanistic target of rapamycin complex 1) regulate downstream gene expression through m6A processing. Conversely, m6A can modulate the activity of signal transduction networks via m6A modification of signaling pathway genes or by acting as a ligand for receptors. In this review, we discuss the current understanding of the crosstalk between m6A and signaling pathways and its implication for biological systems.
Keywords ERK, mTOR, N6-methyladenosine, RNA modification, signaling, TGFβ
Indeed, cells tightly regulate m6A modification using specialized enzymes, m6A writers and erasers. The m6A-modified RNAs then recruit m6A-binding proteins (readers) that guide these RNAs for RNA biogenesis processes such as pre-mRNA splicing, nuclear export, stabilization, degradation, and translation. Aberrant m6A modifications by overactivation or suppression of these enzymes lead to human diseases such as cancer, diabetes, and neurological disorders. There are comprehensive review papers about the molecular functions of m6A enzymes (Meyer and Jaffrey, 2017; Wiener and Schwartz, 2021; Zaccara et al., 2019) and their pathophysiological functions (Barbieri and Kouzarides, 2020; He and He, 2021; Huang et al., 2020; Kasowitz et al., 2018). In this review, we will focus on how the signal transduction pathways, which play key roles in diverse physiological and pathological conditions, coordinate cellular processes through m6A. Given that m6A also controls signaling pathways through RNA modification or acting as a ligand, understanding the crosstalk between signal transduction networks and m6A RNA processing will provide us insights into the complex biological systems.
The m6A methyltransferase (writer) consists of the enzymatic m6A-methyltransferase like (METTL) complex (MAC) and the scaffolding MAC-associated complex (MACOM) (Fig. 1). METTL3 is a catalytic core of MAC, which methylates target mRNAs on the adenosine of DRACH sequence (Bokar et al., 1994; Dominissini et al., 2012; Meyer et al., 2012). METTL14 acts as a scaffolding protein of MAC by recognizing the substrate RNA and interacting with Wilms’ tumor 1-associating protein (WTAP) of MACOM (Bujnicki et al., 2002; Liu et al., 2014; Ping et al., 2014). MACOM consists of several adaptor proteins including WTAP, VIRMA (vir-like m6A methyltransferase associated), RBM15 (RNA-binding motif protein 15), HAKAI, and ZC3H13 (zinc finger CCCH domain-containing protein 13) (Knuckles et al., 2018; Patil et al., 2016; Růžička et al., 2017; Śledź and Jinek, 2016; Wang et al., 2021; Yue et al., 2018). The MACOM complex does not have catalytic activity, but it is necessary for efficient m6A methylation by recruiting RNA substrates and stabilizing the MAC-MACOM complex in the nucleus and nuclear speckles.
Similar to DNA and histone methylations, m6A RNA methylation is a reversible process regulated by the demethylase enzymes (erasers): Fat mass and obesity-associated protein (FTO) and AlkB homolog 5 (ALKBH5) (Jia et al., 2011; Zheng et al., 2013) (Fig. 1). They belong to AlkB homolog iron(II) and αKG-dependent dioxygenases, which include nine proteins with different substrate preferences toward single-stranded (ss) or double-stranded (ds) DNA and RNA substrates (Guengerich, 2015). Interestingly, although FTO and ALKBH5 exhibit comparable catalytic activities for m6A demethylation on ssRNA, the reaction steps are quite different. While ALKBH5 directly converts m6A to adenosine, FTO produces two intermediates
The m6A binding proteins that determine the fates of m6A-methylated mRNA are classified as readers, which include several proteins such as YT521-B homology (YTH) and insulin-like growth factor-2 mRNA-binding protein (IGF2BP) family proteins (Fig. 1). As a common function, YTHDF1/2/3 promote the degradation of m6A-containing mRNAs (Lee et al., 2020; Patil et al., 2018; Wang et al., 2014; Zaccara and Jaffrey, 2020). YTHDF1 and YTHDF3 facilitate protein translation of m6A-methylated mRNAs by promoting ribosome assembly (Li et al., 2017; Shi et al., 2017; Wang et al., 2015). YTHDF2 undergoes liquid-liquid phase separation with mRNAs containing multiple m6A residues (Ries et al., 2019; Wang et al., 2014). Similar to YTHDF1 and YTHDF3, YTHDC2 induces degradation of m6A-modified mRNAs while enhancing their translation efficiency (Hsu et al., 2017; Mao et al., 2019; Tanabe et al., 2016; Wojtas et al., 2017; Zhou et al., 2021). The last YTH family protein, YTHDC1, facilitates pre-mRNA splicing and nuclear export of m6A-modified mRNAs (Kasowitz et al., 2018; Roundtree et al., 2017a; 2017b; Xiao et al., 2016; Xu et al., 2014). In contrast to the YTH family proteins, IGF2BP family proteins increase both the stability and translation efficiency of m6A-modified mRNAs, maximizing the expression of m6A-modified genes (Huang et al., 2018).
TGFβ signaling pathway plays an essential role in cell fate decisions, including pluripotency maintenance, differentiation, senescence, apoptosis, and tumorigenesis (Derynck and Zhang, 2003). TGFβ family proteins are ligands for TGFβ receptors (TGFβR), which includes TGFβ, nodal, activin, bone morphogenetic protein (BMP), and growth differentiation factor (GDF) (Derynck and Zhang, 2003; Zhang et al., 2017). The activated TGFβR phosphorylates downstream signaling proteins, the receptor-regulated SMADs (R-SMADs). Each TGFβ family ligand activates distinct R-SMADs. For example, BMP and GDF promote the phosphorylation of SMAD1, SMAD5, and SMAD8 (Hata and Chen, 2016). On the other hand, TGFβ, nodal, and activin promote SMAD2 and SMAD3 phosphorylation (Hata and Chen, 2016). The phosphorylated R-SMADs form a heterodimer such as SMAD2-SMAD3 (SMAD2/3), which subsequently binds with a common binding partner SMAD4 (co-SMAD). The SMAD complex then moves into the nucleus and associates with transcription factors and chromatin remodeling proteins to induce transcription of target genes (Derynck and Zhang, 2003; Hata and Chen, 2016) (Fig. 2).
In addition to their well-established role in transcription, a novel function of TGFβ-SMAD signaling has been revealed in m6A modification (Fig. 2). TGFβ ligands activin and nodal maintain embryonic stem cell (ESC) stemness through SMAD2/3-dependent pluripotency gene expression (James et al., 2005). On the other hand, considering the rapid differentiation of ESCs upon TGFβ withdrawal, ESCs would have also developed some mechanisms for such a rapid, efficient transition. From interactome analysis of SMAD2/3, Bertero et al. (2018) found that the phosphorylated SMAD2/3 interacts with the m6A writer complex, METTL3, METTL14, and WTAP. m6A methylation of SMAD2/3 target genes,
Interestingly, SMAD2/3 do not directly control the activity of m6A writer complex (Bertero et al., 2018). Instead, SMAD2/3 induce m6A methylation of its target genes by recruiting the m6A writer complex to the active transcription sites (Fig. 2). In another study, Huang et al. (2019) showed that trimethylation of histone H3 at Lys36 (H3K36me3) recruits m6A writer complex to the active transcription elongation sites. It may seem odd that transcription factors and elongation markers recruit m6A enzymes to label newly transcribed mRNAs with m6A for degradation. However, this priming system would be most efficient when a timely cell fate transition is required. For example, in early development of zebrafish embryos, maternal mRNAs are marked with m6A and degraded by the YTHDF2 reader protein during the maternal-to-zygotic transition (Zhao et al., 2017).
TGFβ signaling also triggers m6A modification of target genes for epithelial-mesenchymal transition (EMT) of cancer cells (Fig. 2). TGFβ treatment increases m6A levels in diverse cancer cells, including cervical, liver, breast, and lung cancers. Particularly, TGFβ induces methylation of
ERK is a member of the mitogen-activated protein kinases (MAPKs) family (Roberts and Der, 2007). MAPK pathway is a phosphorylation cascade composed of GTPase-activated kinase (MAPKKK) that phosphorylates and activates an intermediate kinase (MAPKK), which in turn phosphorylates and activates the effector kinase (MAPK) (Lavoie et al., 2020). In the ERK MAPK pathway, the epithermal growth factor (EGF) binds to the receptor tyrosine kinase, EGF receptor (EGFR), to activate RAS GTPase (Boriack-Sjodin et al., 1998). Then the GTP-loaded RAS promotes the kinase activity of RAF (MAPKKK), which is followed by MEK (MAPKK) activation (Lavoie and Therrien, 2015). Finally, MEK activates ERK (MAPK) that controls a wide range of cellular processes by phosphorylating downstream target proteins (Wee and Wang, 2017). Because of its key role in regulating cell proliferation, survival, and differentiation, ERK signaling is one of the frequently activated signaling pathways in human cancers (Davies et al., 2002; Li et al., 2018).
The activity of RNA processing enzymes such as RNA polymerase, splicing enzymes, and translation factors are often regulated by phosphorylation (Nosella and Forman-Kay, 2021; Thapar, 2015). In a recent paper, Sun et al. (2020) provided direct evidence for the phosphorylation-dependent regulation of m6A processing enzymes. To find a new regulator of m6A modification, the authors performed a CRISPR knockout screen using a GFP reporter system that contains m6A modification site on a circular RNA sequence. Once the GFP RNA is assembled by back splicing of the circular RNA (Yang et al., 2017), m6A methylation on that RNA drives translation and expression of GFP. From the screen, several genes in the MAPK signaling pathway were identified (Sun et al., 2020). Mechanistically, it turns out that the effector protein of MAPK signaling pathway, ERK, phosphorylates m6A writer proteins METTL3 and WTAP. The ERK-dependent phosphorylation of METTL3 and WTAP strengthened their interaction; however, ERK did not influence the binding between METTL3 and METTL14 (Sun et al., 2020). The association between METTL3 and METTL14 is strong and not affected by other phosphorylation events either, including a serine residue of METTL14 that forms a salt bridge with METTL3 (Schöller et al., 2018; Wang et al., 2016).
Even though the binding of METTL3-METTL14 was not regulated by ERK-dependent phosphorylation, ERK controlled the activity of MAC complex through METTL3 stabilization (Sun et al., 2020) (Fig. 3). While the non-phosphorylated METTL3 is degraded by ubiquitination, the phosphorylated METTL3 recruits ubiquitin-specific protease 5 (USP5) which removes ubiquitin from METTL3. From the m6A-GFP reporter CRISPR screen, Sun et al. (2020) also found several E3 ubiquitin ligases that decrease m6A levels. Knockdown of ubiquitin ligase candidates
In addition to m6A writers, ERK-dependent regulation of the m6A reader protein is reported (Fig. 3). Fang et al. (2021) found that ERK phosphorylation status correlates with YTHDF2 expression level in the glioblastoma tissues. Upon EGF stimulation or EGFR overexpression, ERK phosphorylates YTHDF2 to induce stabilization of YTHDF2. The stabilized YTHDF2 degrades m6A-modified
Under stress conditions, such as heat shock, m6A modification of mRNA is globally increased (Meyer et al., 2015; Ries et al., 2019; Zhou et al., 2015). Yu et al. (2021) found that m6A rapidly accumulates upon reactive oxygen species (ROS) stress, within five minutes. Such rapid response indicates that ROS may directly influence the activity of m6A enzymes. Surprisingly, ERK signaling plays a crucial role in this stress response by promoting sumoylation of m6A eraser, ALKBH5 (Fig. 3). In hematopoietic stem and progenitor cells, ROS-induced ERK and another MAPK, JNK (c-Jun N-terminal kinase), phosphorylate ALKBH5. ALKBH5 phosphorylation promotes its interaction with UBC9, a SUMO E2 conjugating enzyme, and disassociates ALKBH5 from the desumoylase SENP1. Consequently, the increased sumoylation on ALKBH5 prevents its binding with the substrate mRNAs and leads to the elevation of m6A levels on mRNAs (Yu et al., 2021). Among the m6A readers, the mRNA stabilizing reader, IGF2BP (Huang et al., 2018), plays a major role in ROS-ERK-ALKBH5-dependent gene expression regulation. In response to ROS, IGF2BP1/2/3 stabilize mRNA of
mTORC1 is a serine/threonine kinase that promotes anabolic process including synthesis of proteins, nucleic acids, and lipids (Saxton and Sabatini, 2017; Shimobayashi and Hall, 2014). In response to extracellular stimuli such as growth factors and nutrients, PI3K-Akt initiate the signaling cascade that activates mTORC1. Akt inhibits tuberous sclerosis complex (TSC) 1/2, a GTPase activating protein that inhibits mTORC1-activating small GTPase Rheb (Saxton and Sabatini, 2017; Shimobayashi and Hall, 2014). To turn on the gene expression program, mTORC1 activates RNA processes from transcription, splicing, to translation through its downstream proteins including ribosomal protein S6 kinase (S6K), SRPK (serine/arginine-rich protein specific kinase), eIF4B (eukaryotic translation initiation factor 4B), and 4E-BP (eIF4E-binding protein) (Lee et al., 2017; Ma and Blenis, 2009).
mTORC1’s role in promoting m6A RNA modification has also been elucidated. Cho et al. (2021) and Villa et al. (2021) found that mTORC1 increases expression of the MACOM component, WTAP (Fig. 4). mTORC1/S6K enhances WTAP translation through eIF4A/4B-dependent unwinding of secondary structure in WTAP’s 5’-untranslated region. The elevation of WTAP expression enhances m6A methyltransferase activity (Cho et al., 2021; Villa et al., 2021). In another paper, Tang et al. (2021) found that mTORC1-mediated activation of the chaperonin protein, chaperonin containing tailless complex polypeptide 1 (CCT) complex, stabilizes the MAC component proteins, METTL3 and METTL14 (Fig. 4). In addition to regulating enzyme expression, mTORC1 promotes m6A modification by increasing S-adenosylmethionine (SAM) level (Villa et al., 2021) (Fig. 4). SAM is a methyl donor for m6A, whose increase can stimulate the activity of methyltransferase enzymes (Bokar et al., 1997; Kim and Lee, 2021; Tuck, 1992). mTORC1 induces SAM synthesis through cMyc-dependent expression of MAT2A (methionine adenosyl transferase 2A) (Villa et al., 2021). SAM has been shown to induce condensation of METTL3 through liquid-liquid phase separation and promotes the association of MAC with WTAP at nuclear speckles (Han et al., 2022). Therefore, induction of SAM levels could be another way of increasing the activity of m6A methyltransferase complex by mTORC1. Further studies are necessary to elucidate how mTORC1-dependent expression and localization changes of MAC (METTL3 and METTL14) and MACOM (WTAP) proteins coordinate m6A processing in physiological and disease conditions.
Surprisingly, mTORC1-induced m6A modification induced degradation of mRNAs, which seems to be the opposite of canonical mTORC1 function in promoting macromolecule synthesis. However, those mRNAs methylated by mTORC1 include cell growth-suppressing genes such as cMyc suppressor and autophagy machinery (Cho et al., 2021; Tang et al., 2021). For example, mTORC1 induces methylation of
In the signal transduction cascade, feedback inhibition of upstream signaling by the downstream components plays crucial role in preventing overactivation of signal transduction cascade (Mendoza et al., 2011). Such negative feedback loops are often hijacked by cancers to promote cell proliferation and survival. From a genome-wide sequencing study of endometrial cancers, Liu et al. (2018) found that
In addition to the gene expression regulation, an unexpected function of m6A has been uncovered as an extracellular signaling molecule. Considering the existence of G-protein coupled receptors (GPCR) that are activated by nucleotide ligands such as adenosine receptor (AR) (Borea et al., 2015), Ogawa et al. (2021) performed a screen to identify a new nucleotide ligand for ARs. From the screen, several adenosine derivatives including 1-methyladenosine (m1A), m6A, and m6Am activated adenosine A3 receptor (A3R), with m6A being the most potent activator. In fact, m6A activated A3R approximately 10-fold higher than adenosine, with EC50 (half maximal effective concentration) of 10 nM in contrast to that of adenosine being 100 nM. The ligand binding domain of A3R has hydrophobic amino acids including valine, leucine, and isoleucine, which could form van der Waals interactions with the methyl group on m6A. In contrast, adenosine would be less stable in the ligand binding pocket due to the lack of those intermolecular interactions. Indeed, when the ligand binding pocket of other ARs was mutated to contain those hydrophobic amino acids, they were also activated by m6A (Ogawa et al., 2021).
When m6A was treated onto A3R-expressing cells, the AR downstream signals such as ERK and intracellular calcium transport were induced, which was abolished by AR antagonist. Upon cytotoxic stresses such as ROS that activate AR signaling, m6A was produced in cells by lysosomal degradation of mRNA and rRNA to initiate the m6A-AR-ERK signaling pathway (Fig. 5B) (Ogawa et al., 2021). Although Ogawa et al. (2021) tested the activity of only single nucleoside m6A molecules, it is possible that m6A-containing oligonucleotides can also act as a ligand for GPCRs. The RNase T2 ribonucleases in the lysosome generate both mono- and oligo-nucleotides (Fujiwara et al., 2017) which can be secreted outside of the cells. It will be exciting to dissect the potential roles of m6A-containing nucleic acids as receptor-binding signaling molecules in various developmental and disease processes.
The regulation of m6A modification by multiple signaling pathways demonstrates cells’ abilities to dynamically determine their fate by rewiring gene expression via post-transcriptional gene modifications beyond the gene transcription level. Although research has begun to identify the effects of individual signaling pathways on m6A processing, there remains open questions regarding the potential for crosstalk between interwoven signaling pathways. For example, phosphorylation of METTL3 by ERK signaling induces m6A methylation of pluripotency genes such as
We thank members of the Lee laboratory for helpful discussions, especially Laurence Seabrook, Yujin Chun, and Sunhee Jung, for their feedback on the manuscript. We apologize to authors whose work was not cited due to space limitations. This work was supported by the Department of Defense grant TS200022.
K.-H.J., C.R.H., and G.L. researched data for the review. C.R.H. made substantial contributions to the m6A processing enzyme section and figures. K.-H.J. made substantial contributions to the rest of the review. K.-H.J., C.R.H., and G.L. wrote, reviewed, and edited the manuscript.
The authors have no potential conflicts of interest to disclose.
Mol. Cells 2022; 45(7): 435-443
Published online July 31, 2022 https://doi.org/10.14348/molcells.2022.0017
Copyright © The Korean Society for Molecular and Cellular Biology.
Ki-Hong Jang1 , Chloe R. Heras1,2
, and Gina Lee1,*
1Department of Microbiology and Molecular Genetics, Chao Family Comprehensive Cancer Center, School of Medicine, University of California Irvine, Irvine, CA 92617, USA, 2School of Biological Sciences, University of California Irvine, Irvine, CA 92697, USA
Correspondence to:ginalee@uci.edu
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/.
In response to environmental changes, signaling pathways rewire gene expression programs through transcription factors. Epigenetic modification of the transcribed RNA can be another layer of gene expression regulation. N6-adenosine methylation (m6A) is one of the most common modifications on mRNA. It is a reversible chemical mark catalyzed by the enzymes that deposit and remove methyl groups. m6A recruits effector proteins that determine the fate of mRNAs through changes in splicing, cellular localization, stability, and translation efficiency. Emerging evidence shows that key signal transduction pathways including TGFβ (transforming growth factor-β), ERK (extracellular signal-regulated kinase), and mTORC1 (mechanistic target of rapamycin complex 1) regulate downstream gene expression through m6A processing. Conversely, m6A can modulate the activity of signal transduction networks via m6A modification of signaling pathway genes or by acting as a ligand for receptors. In this review, we discuss the current understanding of the crosstalk between m6A and signaling pathways and its implication for biological systems.
Keywords: ERK, mTOR, N6-methyladenosine, RNA modification, signaling, TGFβ
Indeed, cells tightly regulate m6A modification using specialized enzymes, m6A writers and erasers. The m6A-modified RNAs then recruit m6A-binding proteins (readers) that guide these RNAs for RNA biogenesis processes such as pre-mRNA splicing, nuclear export, stabilization, degradation, and translation. Aberrant m6A modifications by overactivation or suppression of these enzymes lead to human diseases such as cancer, diabetes, and neurological disorders. There are comprehensive review papers about the molecular functions of m6A enzymes (Meyer and Jaffrey, 2017; Wiener and Schwartz, 2021; Zaccara et al., 2019) and their pathophysiological functions (Barbieri and Kouzarides, 2020; He and He, 2021; Huang et al., 2020; Kasowitz et al., 2018). In this review, we will focus on how the signal transduction pathways, which play key roles in diverse physiological and pathological conditions, coordinate cellular processes through m6A. Given that m6A also controls signaling pathways through RNA modification or acting as a ligand, understanding the crosstalk between signal transduction networks and m6A RNA processing will provide us insights into the complex biological systems.
The m6A methyltransferase (writer) consists of the enzymatic m6A-methyltransferase like (METTL) complex (MAC) and the scaffolding MAC-associated complex (MACOM) (Fig. 1). METTL3 is a catalytic core of MAC, which methylates target mRNAs on the adenosine of DRACH sequence (Bokar et al., 1994; Dominissini et al., 2012; Meyer et al., 2012). METTL14 acts as a scaffolding protein of MAC by recognizing the substrate RNA and interacting with Wilms’ tumor 1-associating protein (WTAP) of MACOM (Bujnicki et al., 2002; Liu et al., 2014; Ping et al., 2014). MACOM consists of several adaptor proteins including WTAP, VIRMA (vir-like m6A methyltransferase associated), RBM15 (RNA-binding motif protein 15), HAKAI, and ZC3H13 (zinc finger CCCH domain-containing protein 13) (Knuckles et al., 2018; Patil et al., 2016; Růžička et al., 2017; Śledź and Jinek, 2016; Wang et al., 2021; Yue et al., 2018). The MACOM complex does not have catalytic activity, but it is necessary for efficient m6A methylation by recruiting RNA substrates and stabilizing the MAC-MACOM complex in the nucleus and nuclear speckles.
Similar to DNA and histone methylations, m6A RNA methylation is a reversible process regulated by the demethylase enzymes (erasers): Fat mass and obesity-associated protein (FTO) and AlkB homolog 5 (ALKBH5) (Jia et al., 2011; Zheng et al., 2013) (Fig. 1). They belong to AlkB homolog iron(II) and αKG-dependent dioxygenases, which include nine proteins with different substrate preferences toward single-stranded (ss) or double-stranded (ds) DNA and RNA substrates (Guengerich, 2015). Interestingly, although FTO and ALKBH5 exhibit comparable catalytic activities for m6A demethylation on ssRNA, the reaction steps are quite different. While ALKBH5 directly converts m6A to adenosine, FTO produces two intermediates
The m6A binding proteins that determine the fates of m6A-methylated mRNA are classified as readers, which include several proteins such as YT521-B homology (YTH) and insulin-like growth factor-2 mRNA-binding protein (IGF2BP) family proteins (Fig. 1). As a common function, YTHDF1/2/3 promote the degradation of m6A-containing mRNAs (Lee et al., 2020; Patil et al., 2018; Wang et al., 2014; Zaccara and Jaffrey, 2020). YTHDF1 and YTHDF3 facilitate protein translation of m6A-methylated mRNAs by promoting ribosome assembly (Li et al., 2017; Shi et al., 2017; Wang et al., 2015). YTHDF2 undergoes liquid-liquid phase separation with mRNAs containing multiple m6A residues (Ries et al., 2019; Wang et al., 2014). Similar to YTHDF1 and YTHDF3, YTHDC2 induces degradation of m6A-modified mRNAs while enhancing their translation efficiency (Hsu et al., 2017; Mao et al., 2019; Tanabe et al., 2016; Wojtas et al., 2017; Zhou et al., 2021). The last YTH family protein, YTHDC1, facilitates pre-mRNA splicing and nuclear export of m6A-modified mRNAs (Kasowitz et al., 2018; Roundtree et al., 2017a; 2017b; Xiao et al., 2016; Xu et al., 2014). In contrast to the YTH family proteins, IGF2BP family proteins increase both the stability and translation efficiency of m6A-modified mRNAs, maximizing the expression of m6A-modified genes (Huang et al., 2018).
TGFβ signaling pathway plays an essential role in cell fate decisions, including pluripotency maintenance, differentiation, senescence, apoptosis, and tumorigenesis (Derynck and Zhang, 2003). TGFβ family proteins are ligands for TGFβ receptors (TGFβR), which includes TGFβ, nodal, activin, bone morphogenetic protein (BMP), and growth differentiation factor (GDF) (Derynck and Zhang, 2003; Zhang et al., 2017). The activated TGFβR phosphorylates downstream signaling proteins, the receptor-regulated SMADs (R-SMADs). Each TGFβ family ligand activates distinct R-SMADs. For example, BMP and GDF promote the phosphorylation of SMAD1, SMAD5, and SMAD8 (Hata and Chen, 2016). On the other hand, TGFβ, nodal, and activin promote SMAD2 and SMAD3 phosphorylation (Hata and Chen, 2016). The phosphorylated R-SMADs form a heterodimer such as SMAD2-SMAD3 (SMAD2/3), which subsequently binds with a common binding partner SMAD4 (co-SMAD). The SMAD complex then moves into the nucleus and associates with transcription factors and chromatin remodeling proteins to induce transcription of target genes (Derynck and Zhang, 2003; Hata and Chen, 2016) (Fig. 2).
In addition to their well-established role in transcription, a novel function of TGFβ-SMAD signaling has been revealed in m6A modification (Fig. 2). TGFβ ligands activin and nodal maintain embryonic stem cell (ESC) stemness through SMAD2/3-dependent pluripotency gene expression (James et al., 2005). On the other hand, considering the rapid differentiation of ESCs upon TGFβ withdrawal, ESCs would have also developed some mechanisms for such a rapid, efficient transition. From interactome analysis of SMAD2/3, Bertero et al. (2018) found that the phosphorylated SMAD2/3 interacts with the m6A writer complex, METTL3, METTL14, and WTAP. m6A methylation of SMAD2/3 target genes,
Interestingly, SMAD2/3 do not directly control the activity of m6A writer complex (Bertero et al., 2018). Instead, SMAD2/3 induce m6A methylation of its target genes by recruiting the m6A writer complex to the active transcription sites (Fig. 2). In another study, Huang et al. (2019) showed that trimethylation of histone H3 at Lys36 (H3K36me3) recruits m6A writer complex to the active transcription elongation sites. It may seem odd that transcription factors and elongation markers recruit m6A enzymes to label newly transcribed mRNAs with m6A for degradation. However, this priming system would be most efficient when a timely cell fate transition is required. For example, in early development of zebrafish embryos, maternal mRNAs are marked with m6A and degraded by the YTHDF2 reader protein during the maternal-to-zygotic transition (Zhao et al., 2017).
TGFβ signaling also triggers m6A modification of target genes for epithelial-mesenchymal transition (EMT) of cancer cells (Fig. 2). TGFβ treatment increases m6A levels in diverse cancer cells, including cervical, liver, breast, and lung cancers. Particularly, TGFβ induces methylation of
ERK is a member of the mitogen-activated protein kinases (MAPKs) family (Roberts and Der, 2007). MAPK pathway is a phosphorylation cascade composed of GTPase-activated kinase (MAPKKK) that phosphorylates and activates an intermediate kinase (MAPKK), which in turn phosphorylates and activates the effector kinase (MAPK) (Lavoie et al., 2020). In the ERK MAPK pathway, the epithermal growth factor (EGF) binds to the receptor tyrosine kinase, EGF receptor (EGFR), to activate RAS GTPase (Boriack-Sjodin et al., 1998). Then the GTP-loaded RAS promotes the kinase activity of RAF (MAPKKK), which is followed by MEK (MAPKK) activation (Lavoie and Therrien, 2015). Finally, MEK activates ERK (MAPK) that controls a wide range of cellular processes by phosphorylating downstream target proteins (Wee and Wang, 2017). Because of its key role in regulating cell proliferation, survival, and differentiation, ERK signaling is one of the frequently activated signaling pathways in human cancers (Davies et al., 2002; Li et al., 2018).
The activity of RNA processing enzymes such as RNA polymerase, splicing enzymes, and translation factors are often regulated by phosphorylation (Nosella and Forman-Kay, 2021; Thapar, 2015). In a recent paper, Sun et al. (2020) provided direct evidence for the phosphorylation-dependent regulation of m6A processing enzymes. To find a new regulator of m6A modification, the authors performed a CRISPR knockout screen using a GFP reporter system that contains m6A modification site on a circular RNA sequence. Once the GFP RNA is assembled by back splicing of the circular RNA (Yang et al., 2017), m6A methylation on that RNA drives translation and expression of GFP. From the screen, several genes in the MAPK signaling pathway were identified (Sun et al., 2020). Mechanistically, it turns out that the effector protein of MAPK signaling pathway, ERK, phosphorylates m6A writer proteins METTL3 and WTAP. The ERK-dependent phosphorylation of METTL3 and WTAP strengthened their interaction; however, ERK did not influence the binding between METTL3 and METTL14 (Sun et al., 2020). The association between METTL3 and METTL14 is strong and not affected by other phosphorylation events either, including a serine residue of METTL14 that forms a salt bridge with METTL3 (Schöller et al., 2018; Wang et al., 2016).
Even though the binding of METTL3-METTL14 was not regulated by ERK-dependent phosphorylation, ERK controlled the activity of MAC complex through METTL3 stabilization (Sun et al., 2020) (Fig. 3). While the non-phosphorylated METTL3 is degraded by ubiquitination, the phosphorylated METTL3 recruits ubiquitin-specific protease 5 (USP5) which removes ubiquitin from METTL3. From the m6A-GFP reporter CRISPR screen, Sun et al. (2020) also found several E3 ubiquitin ligases that decrease m6A levels. Knockdown of ubiquitin ligase candidates
In addition to m6A writers, ERK-dependent regulation of the m6A reader protein is reported (Fig. 3). Fang et al. (2021) found that ERK phosphorylation status correlates with YTHDF2 expression level in the glioblastoma tissues. Upon EGF stimulation or EGFR overexpression, ERK phosphorylates YTHDF2 to induce stabilization of YTHDF2. The stabilized YTHDF2 degrades m6A-modified
Under stress conditions, such as heat shock, m6A modification of mRNA is globally increased (Meyer et al., 2015; Ries et al., 2019; Zhou et al., 2015). Yu et al. (2021) found that m6A rapidly accumulates upon reactive oxygen species (ROS) stress, within five minutes. Such rapid response indicates that ROS may directly influence the activity of m6A enzymes. Surprisingly, ERK signaling plays a crucial role in this stress response by promoting sumoylation of m6A eraser, ALKBH5 (Fig. 3). In hematopoietic stem and progenitor cells, ROS-induced ERK and another MAPK, JNK (c-Jun N-terminal kinase), phosphorylate ALKBH5. ALKBH5 phosphorylation promotes its interaction with UBC9, a SUMO E2 conjugating enzyme, and disassociates ALKBH5 from the desumoylase SENP1. Consequently, the increased sumoylation on ALKBH5 prevents its binding with the substrate mRNAs and leads to the elevation of m6A levels on mRNAs (Yu et al., 2021). Among the m6A readers, the mRNA stabilizing reader, IGF2BP (Huang et al., 2018), plays a major role in ROS-ERK-ALKBH5-dependent gene expression regulation. In response to ROS, IGF2BP1/2/3 stabilize mRNA of
mTORC1 is a serine/threonine kinase that promotes anabolic process including synthesis of proteins, nucleic acids, and lipids (Saxton and Sabatini, 2017; Shimobayashi and Hall, 2014). In response to extracellular stimuli such as growth factors and nutrients, PI3K-Akt initiate the signaling cascade that activates mTORC1. Akt inhibits tuberous sclerosis complex (TSC) 1/2, a GTPase activating protein that inhibits mTORC1-activating small GTPase Rheb (Saxton and Sabatini, 2017; Shimobayashi and Hall, 2014). To turn on the gene expression program, mTORC1 activates RNA processes from transcription, splicing, to translation through its downstream proteins including ribosomal protein S6 kinase (S6K), SRPK (serine/arginine-rich protein specific kinase), eIF4B (eukaryotic translation initiation factor 4B), and 4E-BP (eIF4E-binding protein) (Lee et al., 2017; Ma and Blenis, 2009).
mTORC1’s role in promoting m6A RNA modification has also been elucidated. Cho et al. (2021) and Villa et al. (2021) found that mTORC1 increases expression of the MACOM component, WTAP (Fig. 4). mTORC1/S6K enhances WTAP translation through eIF4A/4B-dependent unwinding of secondary structure in WTAP’s 5’-untranslated region. The elevation of WTAP expression enhances m6A methyltransferase activity (Cho et al., 2021; Villa et al., 2021). In another paper, Tang et al. (2021) found that mTORC1-mediated activation of the chaperonin protein, chaperonin containing tailless complex polypeptide 1 (CCT) complex, stabilizes the MAC component proteins, METTL3 and METTL14 (Fig. 4). In addition to regulating enzyme expression, mTORC1 promotes m6A modification by increasing S-adenosylmethionine (SAM) level (Villa et al., 2021) (Fig. 4). SAM is a methyl donor for m6A, whose increase can stimulate the activity of methyltransferase enzymes (Bokar et al., 1997; Kim and Lee, 2021; Tuck, 1992). mTORC1 induces SAM synthesis through cMyc-dependent expression of MAT2A (methionine adenosyl transferase 2A) (Villa et al., 2021). SAM has been shown to induce condensation of METTL3 through liquid-liquid phase separation and promotes the association of MAC with WTAP at nuclear speckles (Han et al., 2022). Therefore, induction of SAM levels could be another way of increasing the activity of m6A methyltransferase complex by mTORC1. Further studies are necessary to elucidate how mTORC1-dependent expression and localization changes of MAC (METTL3 and METTL14) and MACOM (WTAP) proteins coordinate m6A processing in physiological and disease conditions.
Surprisingly, mTORC1-induced m6A modification induced degradation of mRNAs, which seems to be the opposite of canonical mTORC1 function in promoting macromolecule synthesis. However, those mRNAs methylated by mTORC1 include cell growth-suppressing genes such as cMyc suppressor and autophagy machinery (Cho et al., 2021; Tang et al., 2021). For example, mTORC1 induces methylation of
In the signal transduction cascade, feedback inhibition of upstream signaling by the downstream components plays crucial role in preventing overactivation of signal transduction cascade (Mendoza et al., 2011). Such negative feedback loops are often hijacked by cancers to promote cell proliferation and survival. From a genome-wide sequencing study of endometrial cancers, Liu et al. (2018) found that
In addition to the gene expression regulation, an unexpected function of m6A has been uncovered as an extracellular signaling molecule. Considering the existence of G-protein coupled receptors (GPCR) that are activated by nucleotide ligands such as adenosine receptor (AR) (Borea et al., 2015), Ogawa et al. (2021) performed a screen to identify a new nucleotide ligand for ARs. From the screen, several adenosine derivatives including 1-methyladenosine (m1A), m6A, and m6Am activated adenosine A3 receptor (A3R), with m6A being the most potent activator. In fact, m6A activated A3R approximately 10-fold higher than adenosine, with EC50 (half maximal effective concentration) of 10 nM in contrast to that of adenosine being 100 nM. The ligand binding domain of A3R has hydrophobic amino acids including valine, leucine, and isoleucine, which could form van der Waals interactions with the methyl group on m6A. In contrast, adenosine would be less stable in the ligand binding pocket due to the lack of those intermolecular interactions. Indeed, when the ligand binding pocket of other ARs was mutated to contain those hydrophobic amino acids, they were also activated by m6A (Ogawa et al., 2021).
When m6A was treated onto A3R-expressing cells, the AR downstream signals such as ERK and intracellular calcium transport were induced, which was abolished by AR antagonist. Upon cytotoxic stresses such as ROS that activate AR signaling, m6A was produced in cells by lysosomal degradation of mRNA and rRNA to initiate the m6A-AR-ERK signaling pathway (Fig. 5B) (Ogawa et al., 2021). Although Ogawa et al. (2021) tested the activity of only single nucleoside m6A molecules, it is possible that m6A-containing oligonucleotides can also act as a ligand for GPCRs. The RNase T2 ribonucleases in the lysosome generate both mono- and oligo-nucleotides (Fujiwara et al., 2017) which can be secreted outside of the cells. It will be exciting to dissect the potential roles of m6A-containing nucleic acids as receptor-binding signaling molecules in various developmental and disease processes.
The regulation of m6A modification by multiple signaling pathways demonstrates cells’ abilities to dynamically determine their fate by rewiring gene expression via post-transcriptional gene modifications beyond the gene transcription level. Although research has begun to identify the effects of individual signaling pathways on m6A processing, there remains open questions regarding the potential for crosstalk between interwoven signaling pathways. For example, phosphorylation of METTL3 by ERK signaling induces m6A methylation of pluripotency genes such as
We thank members of the Lee laboratory for helpful discussions, especially Laurence Seabrook, Yujin Chun, and Sunhee Jung, for their feedback on the manuscript. We apologize to authors whose work was not cited due to space limitations. This work was supported by the Department of Defense grant TS200022.
K.-H.J., C.R.H., and G.L. researched data for the review. C.R.H. made substantial contributions to the m6A processing enzyme section and figures. K.-H.J. made substantial contributions to the rest of the review. K.-H.J., C.R.H., and G.L. wrote, reviewed, and edited the manuscript.
The authors have no potential conflicts of interest to disclose.
Yun Jeong Lee, Eun-Young Lee, Bo Hee Choi, Hyonchol Jang, Jae-Kyung Myung, and Hye Jin You
Mol. Cells 2020; 43(5): 459-468 https://doi.org/10.14348/molcells.2020.0041Quy Thi Cam Nguyen, Sun-ji Lee, Seo-wha Choi, Yeon-ju Na, Mi-ran Song, Quyen Thi Ngoc Hoang, Seo Young Sim, Min-Sik Kim, Jeong-Il Kim, Moon-Soo Soh, and Soo Young Kim
Mol. Cells 2019; 42(9): 646-660 https://doi.org/10.14348/molcells.2019.0173Jung-Ae Kim
Mol. Cells 2018; 41(4): 264-270 https://doi.org/10.14348/molcells.2018.0076