Exon junction complex shapes the m6A epitranscriptome

Phase transition of WTAP regulates m6A modification of interferon-stimulated genes

N6-methyladenosine (m6A) is the most prevalent modification of mRNA which controls diverse physiological processes. Although m6A modification has been reported to regulate type I interferon (IFN) responses by targeting the mRNA of IFN-β and the interferon-stimulated genes (ISGs), the detailed mechanism of how m6A methyltransferase complex (MTC) rapidly responds to conduct the modification on nascent mRNA during IFN-β stimulation remains largely unclear. Here, we demonstrate that WTAP, the adaptor protein of m6A MTC, undergoes dephosphorylation-regulated phase transition from aggregates to liquid-like condensates under IFN-β stimulation, thereby mediating m6A modification of a subset of ISGs to restrict their expression. The phase transition of WTAP promotes the interaction with nucleus-translocated transcription factor STAT1, recruits MTC to the promoter regions of ISGs and directs the co-transcriptional m6A modification on ISG mRNAs. Collectively, our findings reveal a novel regulatory role of WTAP phase transition in manipulating signaling pathways and fine-tuning immune response by orchestrating dynamic m6A modification through the cooperation of transcription factors and MTC. Our findings unveil a novel mechanism by which WTAP phase transition controls immune homeostasis via transcription factor-MTC-driven dynamic m6A modification, thereby proposing a potential therapeutic target for alleviating immune dysregulation.

m6A modification is incorporated into bacterial mRNA without specific functional benefit

N 6-Methyladenosine (m6A), the most abundant modification in eukaryotic messenger RNAs (mRNAs), has also been found at a low level in bacterial mRNAs. However, enzyme(s) that introduce m6A modification on mRNAs in bacteria remain elusive. In this work, we combine deep-sequencing approaches that identify m6A sites with in vitro biochemical studies to identify putative m6A methyltransferases that would modify Escherichia coli mRNAs. We tested four uncharacterized candidates predicted to encode proteins with putative methyltransferase domains, whose deletion decreased the m6A level. However, in vitro analysis with the purified putative methyltransferases revealed that none of them installs m6A on mRNA. Exposure to heat and oxidative stress also changed the m6A level; however, we found no clear correlation between the m6A change and the specific stress. Considering two deep-sequencing approaches with different resolution, we found that m6A methylation on bacterial mRNAs is very low and appears randomly introduced. These results suggest that, in contrast to eukaryotes, the m6A modification in bacterial mRNA lacks a direct enzymatic recognition mechanism and has no clear biological function.

Nanopore direct RNA sequencing of human transcriptomes reveals the complexity of mRNA modifications and crosstalk between regulatory features

The identification and functional characterization of chemical modifications on an mRNA molecule, in particular N6-methyladenosine (m6A) modification, significantly broadened our understanding of RNA function and regulation. While interactions between RNA modifications and other RNA features have been proposed, direct evidence showing correlation is limited. Here, using Oxford Nanopore long-read direct RNA sequencing (dRNA-seq), we simultaneously interrogate the transcriptome and epitranscriptome of a human leukemia cell line to investigate the correlation between m6A modifications, mRNA abundance, mRNA stability, polyadenylation (poly(A)) tail length, and alternative splicing. High-quality dRNA-seq is important for unbiased and large-scale correlative analyses. Global assessments indicated a negative association between poly(A) tail length and mRNA abundance while uncovering pathway-specific responses upon depletion of the m6A-forming enzyme METTL3. Overall, our study presented a rich dRNA-seq data resource that has been validated and can be further exploited to inquire into the complexity of RNA modifications and potential interplays between RNA regulatory elements.

Research progress on m6A and drug resistance in gastrointestinal tumors

Gastrointestinal (GI) tumors represent a significant global health burden and are among the leading causes of cancer-related mortality worldwide. their drug resistance is one of the major challenges in cancer therapy. In recent years, epigenetic modifications, especially N6-methyladenosine (m6A) RNA modifications, have become a hot research topic. m6A modification plays an important role in gene expression and cancer progression by regulating RNA splicing, translation, stability, and degradation, which are regulated by "writers," "erasers" and "readers." In GI tumors, resistance to chemotherapy, targeted therapy, and immunotherapy is closely associated with m6A RNA modification. Therefore, the molecular mechanism of m6A modification and its targeted drug development provide new therapeutic strategies for overcoming drug resistance and therapeutic efficacy in GI tumors. In this review, the biological functions of m6A were explored, the specific resistance mechanisms of m6A in different types of GI tumors were explored, new ideas and targets for future treatment resistance were identified, and the limitations of this field were highlighted.

tRNA modifications tune m6A-dependent mRNA decay

Chemically modified nucleotides in mRNA are critical regulators of gene expression, primarily through interactions with reader proteins that bind to these modifications. Here, we present a mechanism by which the epitranscriptomic mark N6-methyladenosine (m6A) is read by tRNAs during translation. Codons that are modified with m6A are decoded inefficiently by the ribosome, rendering them "non-optimal" and inducing ribosome collisions on cellular transcripts. This couples mRNA translation to decay. 5-Methoxycarbonylmethyl-2-thiouridine (mcm5s2U) in the tRNA anticodon loop counteracts this effect. This unanticipated link between the mRNA and tRNA epitranscriptomes enables the coordinated decay of mRNA regulons, including those encoding oncogenic signaling pathways. In cancer, dysregulation of the m6A and mcm5s2U biogenesis pathways-marked by a shift toward more mcm5s2U-is associated with more aggressive tumors and poor prognosis. Overall, this pan-epitranscriptomic interaction represents a novel mechanism of post-transcriptional gene regulation with implications for human health.

Long-read RNA sequencing reveals allele-specific N 6-methyladenosine modifications

Long-read sequencing technology enables highly accurate detection of allele-specific RNA expression, providing insights into the effects of genetic variation on splicing and RNA abundance. Furthermore, the ability to directly sequence RNA enables the detection of RNA modifications in tandem with ascertaining the allelic origin of each molecule. Here, we leverage these advantages to determine allele-biased patterns of N 6-methyladenosine (m6A) modifications in native mRNA. We used human and mouse cells with known genetic variants to assign the allelic origin of each mRNA molecule combined with a supervised machine learning model to detect read-level m6A modification ratios. Our analyses reveal the importance of sequences adjacent to the DRACH motif in determining m6A deposition, in addition to allelic differences that directly alter the motif. Moreover, we discover allele-specific m6A modification events with no genetic variants in close proximity to the differentially modified nucleotide, demonstrating the unique advantage of using long-reads and surpassing the capabilities of antibody-based short-read approaches. This technological advance will further our understanding of the role of genetics in determining mRNA modifications.

Spatial control of m6A deposition on enhancer and promoter RNAs through co-acetylation of METTL3 and H3K27 on chromatin

Interaction between the N6-methyladenosine (m6A) methyltransferase METTL3 and METTL14 is critical for METTL3 to deposit m6A on various types of RNAs. It remains to be uncovered whether there is spatial control of m6A deposition on different types of RNAs. Here, through genome-wide CRISPR-Cas9 screening in the A549 cell line, we find that H3K27ac acetylase p300-mediated METTL3 acetylation suppresses the binding of METTL3 on H3K27ac-marked chromatin by inhibiting its interaction with METTL14. Consistently, p300 catalyzing the acetylation of METTL3 specifically occurs on H3K27ac-marked chromatin. Disruptive mutations on METTL3 acetylation sites selectively promote the m6A of chromatin-associated RNAs from p300-bound enhancers and promoters marked by H3K27ac, resulting in transcription inhibition of ferroptosis-inhibition-related genes. In addition, PAK2 promotes METTL3 acetylation by phosphorylating METTL3. Inhibition of PAK2 promotes ferroptosis in a manner that depends on the acetylation of METTL3. Our study reveals a spatial-selective way to specifically regulate the deposition of m6A on enhancer and promoter RNAs.

Identification of Endogenous Sequences Interacting with METTL3/METTL14 RNA Methyltransferase

N6-methyladenosine (m6A) is the most abundant RNA modification in mRNA and regulates various biological processes. The RNA-binding properties of m6A writer proteins play an important role in determining RNA modification sites. METTL3 and METTL14 form the core of the m6A writer complex, with METTL3 as the catalytic methyltransferase and METTL14 as the RNA-binding scaffold. Thus far, the comprehensive RNA-binding properties of METTL3/14 remain unknown. Using RNA-capturing microsphere particles (R-CAMPs), immobilizing RNA fragments derived from endogenous RNAs of human lung carcinoma cells, RNA fragments that interacted with the METTL3/14 methyltransferase domain are isolated. Bioinformatics analysis reveals that the pool of isolated sequences contains significantly more regions with the potential to form RNA G-quadruplexes (rG4s) than the randomly extracted RNA sequence pool and that the (GGA) repeat sequences are most enriched. Circular dichroism spectroscopy, gel mobility shift assays, and methylation experiments demonstrate that METTL3/14 binds to RNA sequences containing (GGA) repeats that form rG4 structure much stronger than RNAs that do not form rG4s. This study shows the potential of RNA G-quadruplex structures as a modulator of epitranscriptomic modification of m6As.

Single-molecule m6A detection empowered by endogenous labeling unveils complexities across RNA isoforms

The landscape of N6-methyadenosine (m6A) on different RNA isoforms is still incompletely understood. Here, in HEK293T cells, we endogenously label the methylated m6A sites on single Oxford Nanopore Technology (ONT) direct RNA sequencing (DRS) reads by APOBEC1-YTH-induced C-to-U mutations 10-100 nt away, obtaining 1,020,237 5-mer single-read m6A signals. We then trained m6Aiso, a deep residual neural network model that accurately identifies and quantifies m6A at single-read resolution. Analyzing m6Aiso-determined m6A on single reads and isoforms uncovers distance-dependent linkages of m6A sites along single molecules. It also uncovers specific methylation of identical m6A sites on intron-retained isoforms, partly due to their differential distances to exon junctions and isoform-specific binding of TARBP2. Moreover, we find that transcription factor SMAD3 promotes m6A deposition on its transcribed RNA isoforms during epithelial-mesenchymal transition, resulting in isoform-specific regulation of m6A on isoforms with alternative promoters. Our study underscores the effectiveness of m6Aiso in elucidating the intricate dynamics and complexities of m6A across RNA isoforms.

YTHDC1 cooperates with the THO complex to prevent RNA-damage-induced DNA breaks

Certain environmental toxins and chemotherapeutics are nucleic acid-damaging agents, causing adducts in DNA and RNA. While most of these adducts occur in RNA, the consequences of RNA damage are largely unexplored. Here, we demonstrate that nuclear RNA damage can result in loss of genome integrity in human cells. Specifically, we show that YTHDC1 regulates alkylation damage responses with the THO complex (THOC). In addition to its established binding to N6-methyladenosine (m6A), YTHDC1 binds to chemically induced N1-methyladenosine (m1A). Without YTHDC1, cells have greater alkylation damage sensitivity and increased DNA breaks, which are rescued by an RNA-specific dealkylase. These RNA-damage-induced DNA breaks (RDIBs) depend on R-loop formation, which is converted to DNA breaks by the XPG nuclease. Strikingly, in the absence of YTHDC1 or THOC, a nuclear RNA m1A methyltransferase is sufficient to induce DNA breaks. Our results provide mechanistic insight into how damaged RNAs can impact genomic integrity.

PYM1 limits non-canonical Exon Junction Complex occupancy in a gene architecture dependent manner to tune mRNA expression

The Exon Junction Complex (EJC) deposited upstream of exon-exon junctions during pre-mRNA splicing in the nucleus remains stably bound to RNA to modulate mRNA fate at multiple post-transcriptional steps until its disassembly during translation. Here, we investigated two EJC disassembly mechanisms in human embryonic kidney 293 (HEK293) cells, one mediated by PYM1, a factor that can bind both the ribosome and the RBM8A/MAGOH heterodimer of the EJC core, and another by the elongating ribosome itself. We find that EJCs lacking PYM1 interaction show no defect in translation-dependent disassembly but is required for translation-independent EJC destabilization. Surprisingly, PYM1 interaction deficient EJCs are enriched on sites away from the canonical EJC binding position including on transcripts without introns or with fewer and longer exons. Acute reduction of PYM1 levels in HEK293 cells results in a modest inhibition of nonsense-mediated mRNA decay and stabilization of mRNAs that localize to endoplasmic reticulum associated TIS-granules and are characterized by fewer and longer exons. We confirmed the previously reported PYM1-flavivirus capsid protein interaction and found that human cells expressing the capsid protein or infected with flaviviruses show similar changes in gene expression as upon PYM1 depletion. Thus, PYM1 acts as an EJC specificity factor that is hijacked by flaviviruses to alter global EJC occupancy and reshape host cell mRNA regulation.

Uncovering the isoform-resolution kinetic landscape of nonsense-mediated mRNA decay with EZbakR

Cellular RNA levels are a product of synthesis and degradation kinetics, which can differ among transcripts of the same gene. An important cause of isoform-specific decay is the nonsense-mediated mRNA decay (NMD) pathway, which degrades transcripts with premature termination codons (PTCs) and other features. Understanding NMD functions requires strategies to quantify isoform kinetics; however, current approaches remain limited. Methods like nucleotide-recoding RNA-seq (NR-seq) enable insights into RNA kinetics, but existing bioinformatic tools do not provide robust, isoform-specific degradation rate constant estimates. We extend the EZbakR-suite by implementing a strategy to infer isoform-level kinetics from short-read NR-seq data. This approach uncovers unexpected variability in NMD efficiency among transcripts with conserved PTC-containing exons and rapid decay of a subset of mRNAs lacking PTCs. Our findings highlight the effects of competition between NMD and other decay pathways, provide mechanistic insights into established NMD efficiency correlates, and identify transcript features promoting efficient decay.

De novo basecalling of RNA modifications at single molecule and nucleotide resolution

RNA modifications influence RNA function and fate, but detecting them in individual molecules remains challenging for most modifications. Here we present a novel methodology to generate training sets and build modification-aware basecalling models. Using this approach, we develop the m6ABasecaller, a basecalling model that predicts m6A modifications from raw nanopore signals. We validate its accuracy in vitro and in vivo, revealing stable m6A modification stoichiometry across isoforms, m6A co-occurrence within RNA molecules, and m6A-dependent effects on poly(A) tails. Finally, we demonstrate that our method generalizes to other RNA and DNA modifications, paving the path towards future efforts detecting other modifications.

N-6-methyladenosine (m6A) promotes the nuclear retention of mRNAs with intact 5' splice site motifs

In humans, misprocessed mRNAs containing intact 5' Splice Site (5'SS) motifs are nuclear retained and targeted for decay by ZFC3H1, a component of the Poly(A) Exosome Targeting complex, and U1-70K, a component of the U1 snRNP. In S. pombe, the ZFC3H1 homolog, Red1, binds to the YTH domain-containing protein Mmi1 and targets certain RNA transcripts to nuclear foci for nuclear retention and decay. Here we show that YTHDC1 and YTHDC2, two YTH domain-containing proteins that bind to N-6-methyladenosine (m6A) modified RNAs, interact with ZFC3H1 and U1-70K, and are required for the nuclear retention of mRNAs with intact 5'SS motifs. Disruption of m6A deposition inhibits both the nuclear retention of these transcripts and their accumulation in YTHDC1-enriched foci that are adjacent to nuclear speckles. Endogenous RNAs with intact 5'SS motifs, such as intronic poly-adenylated transcripts, tend to be m6A-modified at low levels. Thus, the m6A modification acts on a conserved quality control mechanism that targets misprocessed mRNAs for nuclear retention and decay.

Decoupling actin assembly from microtubule disassembly by TBC1D3C-mediated direct GEF-H1 activation

Actin and microtubules are essential cytoskeletal components and coordinate their dynamics through multiple coupling and decoupling mechanisms. However, how actin and microtubule dynamics are decoupled remains incompletely understood. Here, we identified TBC1D3C as a new regulator that can decouple actin filament assembly from microtubule disassembly. We showed that TBC1D3C induces the release of GEF-H1 from microtubules into the cytosol without perturbing microtubule arrays, leading to RhoA activation and actin filament assembly. Mechanistically, we found that TBC1D3C directly binds to GEF-H1, disrupting its interaction with the Tctex-DIC-14-3-3 complex and thereby displacing GEF-H1 from microtubules independently of microtubule disassembly. Super-resolution microscopy and live-cell imaging further confirmed that TBC1D3C triggers GEF-H1 release and actin filament assembly while maintaining microtubule integrity. Therefore, our findings demonstrated that TBC1D3C functions as a direct GEF activator and a novel regulator in decoupling actin assembly from microtubule disassembly, providing new insights into cytoskeletal regulation.

Epitranscriptomic RNA m6A Modification in Cancer Therapy Resistance: Challenges and Unrealized Opportunities

Significant advances in the development of new cancer therapies have given rise to multiple novel therapeutic options in chemotherapy, radiotherapy, immunotherapy, and targeted therapies. Although the development of resistance is often reported along with temporary disease remission, there is often tumor recurrence of an even more aggressive nature. Resistance to currently available anticancer drugs results in poor overall and disease-free survival rates for cancer patients. There are multiple mechanisms through which tumor cells develop resistance to therapeutic agents. To date, efforts to overcome resistance have only achieved limited success. Epitranscriptomics, especially related to m6A RNA modification dysregulation in cancer, is an emerging mechanism for cancer therapy resistance. Here, recent studies regarding the contributions of m6A modification and its regulatory proteins to the development of resistance to different cancer therapies are comprehensively reviewed. The promise and potential limitations of targeting these entities to overcome resistance to various anticancer therapies are also discussed.

m6A sites in the coding region trigger translation-dependent mRNA decay

N6-Methyladenosine (m6A) is the predominant internal RNA modification in eukaryotic messenger RNAs (mRNAs) and plays a crucial role in mRNA stability. Here, using human cells, we reveal that m6A sites in the coding sequence (CDS) trigger CDS-m6A decay (CMD), a pathway that is distinct from previously reported m6A-dependent degradation mechanisms. Importantly, CDS m6A sites act considerably faster and more efficiently than those in the 3' untranslated region, which to date have been considered the main effectors. Mechanistically, CMD depends on translation, whereby m6A deposition in the CDS triggers ribosome pausing and transcript destabilization. The subsequent decay involves the translocation of the CMD target transcripts to processing bodies (P-bodies) and recruitment of the m6A reader protein YT521-B homology domain family protein 2 (YTHDF2). Our findings highlight CMD as a previously unknown pathway, which is particularly important for controlling the expression of developmental regulators and retrogenes.

Single-cell m6A profiling in the mouse brain uncovers cell type-specific RNA methylomes and age-dependent differential methylation

N6-methyladenosine (m6A) is an abundant mRNA modification in the brain that has important roles in neurodevelopment and brain function. However, because of technical limitations, global profiling of m6A sites within the individual cell types that make up the brain has not been possible. Here, we develop a mouse model that enables transcriptome-wide m6A detection in any tissue of interest at single-cell resolution. We use these mice to map m6A across different brain regions and within single cells of the mouse cortex and discover a high degree of shared methylation across brain regions and cell types. However, we also identify a small number of differentially methylated mRNAs in neurons that encode important regulators of neuronal signaling, and we discover that microglia have lower levels of m6A than other cell types. Finally, we perform single-cell m6A mapping in aged mice and identify many transcripts with age-dependent changes in m6A.

METTL3/MYCN cooperation drives neural crest differentiation and provides therapeutic vulnerability in neuroblastoma

Neuroblastoma (NB) is the most common extracranial childhood cancer, caused by the improper differentiation of developing trunk neural crest cells (tNCC) in the sympathetic nervous system. The N6-methyladenosine (m6A) epitranscriptomic modification controls post-transcriptional gene expression but the mechanism by which the m6A methyltransferase complex METTL3/METTL14/WTAP is recruited to specific loci remains to be fully characterized. We explored whether the m6A epitranscriptome could fine-tune gene regulation in migrating/differentiating tNCC. We demonstrate that the m6A modification regulates the expression of HOX genes in tNCC, thereby contributing to their timely differentiation into sympathetic neurons. Furthermore, we show that posterior HOX genes are m6A modified in MYCN-amplified NB with reduced expression. In addition, we provide evidence that sustained overexpression of the MYCN oncogene in tNCC drives METTL3 recruitment to a specific subset of genes including posterior HOX genes creating an undifferentiated state. Moreover, METTL3 depletion/inhibition induces DNA damage and differentiation of MYCN overexpressing cells and increases vulnerability to chemotherapeutic drugs in MYCN-amplified patient-derived xenografts (PDX) in vivo, suggesting METTL3 inhibition could be a potential therapeutic approach for NB.

RNA modifications: emerging players in the regulation of reproduction and development

The intricate world of RNA modifications, collectively termed the epitranscriptome, covers over 170 identified modifications and impacts RNA metabolism and, consequently, almost all biological processes. In this review, we focus on the regulatory roles and biological functions of a panel of dominant RNA modifications (including m 6A, m 5C, Ψ, ac 4C, m 1A, and m 7G) on three RNA types-mRNA, tRNA, and rRNA-in mammalian development, particularly in the context of reproduction as well as embryonic development. We discuss in detail how those modifications, along with their regulatory proteins, affect RNA processing, structure, localization, stability, and translation efficiency. We also highlight the associations among dysfunctions in RNA modification-related proteins, abnormal modification deposition and various diseases, emphasizing the roles of RNA modifications in critical developmental processes such as stem cell self-renewal and cell fate transition. Elucidating the molecular mechanisms by which RNA modifications influence diverse developmental processes holds promise for developing innovative strategies to manage developmental disorders. Finally, we outline several unexplored areas in the field of RNA modification that warrant further investigation.

The crosstalk between alternative splicing and circular RNA in cancer: pathogenic insights and therapeutic implications

RNA splicing is a fundamental step of gene expression. While constitutive splicing removes introns and joins exons unbiasedly, alternative splicing (AS) selectively determines the assembly of exons and introns to generate RNA variants corresponding to the same transcript. The biogenesis of circular RNAs (circRNAs) is inextricably associated with AS. Back-splicing, the biogenic process of circRNA, is a special form of AS. In cancer, both AS and circRNA deviate from the original track. In the present review, we delve into the intricate interplay between AS and circRNAs in the context of cancer. The relationship between AS and circRNAs is intricate, where AS modulates the biogenesis of circRNAs and circRNAs in return regulate AS events. Beyond that, epigenetic and posttranscriptional modifications concurrently regulate AS and circRNAs. On the basis of this modality, we summarize current knowledge on how splicing factors and other RNA binding proteins regulate circRNA biogenesis, and how circRNAs interact with splicing factors to influence AS events. Specifically, the feedback loop regulation between circRNAs and AS events contributes greatly to oncogenesis and cancer progression. In summary, resolving the crosstalk between AS and circRNA will not only provide better insight into cancer biology but also provoke novel strategies to combat cancer.

The N6-methyladenosine Epitranscriptomic Landscape of Lung Adenocarcinoma

Comprehensive N6-methyladenosine (m6A) epitranscriptomic profiling of primary tumors remains largely uncharted. Here, we profiled the m6A epitranscriptome of 10 nonneoplastic lung tissues and 51 lung adenocarcinoma (LUAD) tumors, integrating the corresponding transcriptomic, proteomic, and extensive clinical annotations. We identified distinct clusters and genes that were exclusively linked to disease progression through m6A modifications. In comparison with nonneoplastic lung tissues, we identified 430 transcripts to be hypo-methylated and 222 to be hyper-methylated in tumors. Among these genes, EML4 emerged as a novel metastatic driver, displaying significant hypermethylation in tumors. m6A modification promoted the translation of EML4, leading to its widespread overexpression in primary tumors. Functionally, EML4 modulated cytoskeleton dynamics by interacting with ARPC1A, enhancing lamellipodia formation, cellular motility, local invasion, and metastasis. Clinically, high EML4 protein abundance correlated with features of metastasis. METTL3 small-molecule inhibitor markedly diminished both EML4 m6A and protein abundance and efficiently suppressed lung metastases in vivo. Significance: Our study reveals a dynamic and functional epitranscriptomic landscape in LUAD, offering a valuable resource for further research in the field. We identified EML4 hypermethylation as a key driver of tumor metastasis, highlighting a novel therapeutic strategy of targeting EML4 to prevent LUAD metastasis.

METTL3-STAT5B interaction facilitates the co-transcriptional m6A modification of mRNA to promote breast tumorigenesis

Enhanced expression of methyltransferase-like 3 (METTL3) promotes the m6A modification of specific mRNAs, contributing to breast tumorigenesis. While the mRNA substrates targeted by METTL3 are well characterized, the factors dictating the selection of these specific mRNA remain elusive. This study aimed to examine the regulatory role of the transcription factor STAT5B in METTL3-induced m6A modification. METTL3 specifically interacts with STAT5B in response to mitogenic stimulation by epidermal growth factor (EGF). Chromatin immunoprecipitation and CRISPR/Cas9 mutagenesis showed that STAT5B recruits METTL3 to gene promoters like CCND1, where METTL3 interacts with RPB1, dependent on CDK9-mediated RPB1 (Ser2) phosphorylation during transcription elongation. Inhibition and depletion of either STAT5B or CDK9 prevented the EGF-induced m6A modification of CCND1. The translation efficiency of CCND1 was increased following m6A modification, thereby increasing cell proliferation. STAT5B facilitated METTL3-induced tumor formation by increasing CCND1 expression in an orthotopic mouse model. In clinical context, a positive correlation was observed between p-STAT5B and METTL3 expression in high-grade breast tumors. This study elucidates a novel mechanism that underlies the specificity of m6A modification in breast cancer cells, thereby underscoring its potential therapeutic value.

Timing is everything: When is m6A deposited?

In this issue of Molecular Cell, Tang et al. suggest that m6A deposition is predominantly post-transcriptional.1 They further propose that nuclear dwell time dictates the post-transcriptional accumulation of m6A. These findings have important implications for m6A biogenesis and function.

Nuclear retention coupled with sequential polyadenylation dictates post-transcriptional m6A modification in the nucleus

N6-methyladenosine (m6A) modification is deemed to be co-transcriptionally installed on pre-mRNAs, thereby influencing various downstream RNA metabolism events. However, the causal relationship between m6A modification and RNA processing is often unclear, resulting in premature or even misleading generalizations on the function of m6A modification. Here, we develop 4sU-coupled m6A-level and isoform-characterization sequencing (4sU-m6A-LAIC-seq) and 4sU-GLORI to quantify the m6A levels for both newly synthesized and steady-state RNAs at transcript and single-base-resolution levels, respectively, which enable dissecting the relationship between m6A modification and alternative RNA polyadenylation. Unexpectedly, our results show that many m6A addition events occur post-transcriptionally, especially on transcripts with high m6A levels. Importantly, we find higher m6A levels on shorter 3' UTR isoforms, which likely result from sequential polyadenylation of longer 3' UTR isoforms with prolonged nuclear dwelling time. Therefore, m6A modification can also take place post-transcriptionally to intimately couple with other key RNA metabolism processes to establish and dynamically regulate epi-transcriptomics in mammalian cells.

Modulation of m6A RNA modification by DAP3 in cancer cells

N6-methyladenosine (m6A) RNA methylation is a prevalent RNA modification that significantly impacts RNA metabolism and cancer development. Maintaining the global m6A levels in cancer cells relies on RNA accessibility to methyltransferases and the availability of the methyl donor S-adenosylmethionine (SAM). Here, we reveal that death associated protein 3 (DAP3) plays a crucial role in preserving m6A levels through two distinct mechanisms. First, although DAP3 is not a component of the m6A writer complex, it directly binds to m6A target regions, thereby facilitating METTL3 binding. Second, DAP3 promotes MAT2A's last intron splicing, increasing MAT2A protein, cellular SAM, and m6A levels. Silencing DAP3 hinders tumorigenesis, which can be rescued by MAT2A overexpression. This evidence suggests DAP3's role in tumorigenesis, partly through m6A regulation. Our findings unveil DAP3's complex role as an RNA-binding protein and tumor promoter, impacting RNA processing, splicing, and m6A modification in cancer transcriptomes.

Long-read RNA sequencing reveals allele-specific N6-methyladenosine modifications

Long-read sequencing technology enables highly accurate detection of allele-specific RNA expression, providing insights into the effects of genetic variation on splicing and RNA abundance. Furthermore, the ability to directly sequence RNA promises the detection of RNA modifications in tandem with ascertaining the allelic origin of each molecule. Here, we leverage these advantages to determine allele-biased patterns of N6-methyladenosine (m6A) modifications in native mRNA. We utilized human and mouse cells with known genetic variants to assign allelic origin of each mRNA molecule combined with a supervised machine learning model to detect read-level m6A modification ratios. Our analyses revealed the importance of sequences adjacent to the DRACH-motif in determining m6A deposition, in addition to allelic differences that directly alter the motif. Moreover, we discovered allele-specific m6A modification (ASM) events with no genetic variants in close proximity to the differentially modified nucleotide, demonstrating the unique advantage of using long reads and surpassing the capabilities of antibody-based short-read approaches. This technological advancement promises to advance our understanding of the role of genetics in determining mRNA modifications.

Metabolic regulation of mRNA splicing

Alternative mRNA splicing enables the diversification of the proteome from a static genome and confers plasticity and adaptiveness on cells. Although this is often explored in development, where hard-wired programs drive the differentiation and specialization, alternative mRNA splicing also offers a way for cells to react to sudden changes in outside stimuli such as small-molecule metabolites. Fluctuations in metabolite levels and availability in particular convey crucial information to which cells react and adapt. We summarize and highlight findings surrounding the metabolic regulation of mRNA splicing. We discuss the principles underlying the biochemistry and biophysical properties of mRNA splicing, and propose how these could intersect with metabolite levels. Further, we present examples in which metabolites directly influence RNA-binding proteins and splicing factors. We also discuss the interplay between alternative mRNA splicing and metabolite-responsive signaling pathways. We hope to inspire future research to obtain a holistic picture of alternative mRNA splicing in response to metabolic cues.

Comparative analyses suggest a link between mRNA splicing, stability, and RNA covalent modifications in flowering plants

BACKGROUND

In recent years, covalent modifications on RNA nucleotides have emerged as pivotal moieties influencing the structure, function, and regulatory processes of RNA Polymerase II transcripts such as mRNAs and lncRNAs. However, our understanding of their biological roles and whether these roles are conserved across eukaryotes remains limited.

RESULTS

In this study, we leveraged standard polyadenylation-enriched RNA-sequencing data to identify and characterize RNA modifications that introduce base-pairing errors into cDNA reads. Our investigation incorporated data from three Poaceae (Zea mays, Sorghum bicolor, and Setaria italica), as well as publicly available data from a range of stress and genetic contexts in Sorghum and Arabidopsis thaliana. We uncovered a strong enrichment of RNA covalent modifications (RCMs) deposited on a conserved core set of nuclear mRNAs involved in photosynthesis and translation across these species. However, the cohort of modified transcripts changed based on environmental context and developmental program, a pattern that was also conserved across flowering plants. We determined that RCMs can partly explain accession-level differences in drought tolerance in Sorghum, with stress-associated genes receiving a higher level of RCMs in a drought tolerant accession. To address function, we determined that RCMs are significantly enriched near exon junctions within coding regions, suggesting an association with splicing. Intriguingly, we found that these base-pair disrupting RCMs are associated with stable mRNAs, are highly correlated with protein abundance, and thus likely associated with facilitating translation.

CONCLUSIONS

Our data point to a conserved role for RCMs in mRNA stability and translation across the flowering plant lineage.

Base-resolution m5C profiling across the mammalian transcriptome by bisulfite-free enzyme-assisted chemical labeling approach

5-methylcytosine (m5C) is a prevalent RNA modification crucial for gene expression regulation. However, accurate and sensitive m5C sites identification remains challenging due to severe RNA degradation and reduced sequence complexity during bisulfite sequencing (BS-seq). Here, we report m5C-TAC-seq, a bisulfite-free approach combining TET-assisted m5C-to-f5C oxidation with selective chemical labeling, therefore enabling direct base-resolution m5C detection through pre-enrichment and C-to-T transitions at m5C sites. With m5C-TAC-seq, we comprehensively profiled the m5C methylomes in human and mouse cells, identifying a substantially larger number of confident m5C sites. Through perturbing potential m5C methyltransferases, we deciphered the responsible enzymes for most m5C sites, including the characterization of NSUN5's involvement in mRNA m5C deposition. Additionally, we characterized m5C dynamics during mESC differentiation. Notably, the mild reaction conditions and preservation of nucleotide composition in m5C-TAC-seq allow m5C detection in chromatin-associated RNAs. The accurate and robust m5C-TAC-seq will advance research into m5C methylation functional investigation.

Decoding the specificity of m6A RNA methylation and its implication in cancer therapy

N6-methyladenosine (m6A) is the most abundant endogenous modification in eukaryotic RNAs. It plays important roles in various biological processes and diseases, including cancers. More and more studies have revealed that the deposition of m6A is specifically regulated in a context-dependent manner. Here, we review the diverse mechanisms that determine the topology of m6A along RNAs and the cell-type-specific m6A methylomes. The exon junction complex (EJC) as well as histone modifications play important roles in determining the topological distribution of m6A along nascent RNAs, while the transcription factors and RNA-binding proteins, which usually bind specific DNAs and RNAs in a cell-type-specific manner, largely account for the cell-type-specific m6A methylomes. Due to the lack of specificity of m6A writers and readers, there are still challenges to target the core m6A machinery for cancer therapies. Therefore, understanding the mechanisms underlying the specificity of m6A modifications in cancers would be important for future cancer therapies through m6A intervention.

Roles of N6-methyladenosine writers, readers and erasers in the mammalian germline

N6-methyladenosine (m6A) is the most abundant internal modification of mRNAs in eukaryotes. Numerous studies have shown that m6A plays key roles in many biological and pathophysiological processes, including fertility. The factors involved in m6A-dependent mRNA regulation include writers, which deposit the m6A mark, erasers, which remove it, and readers, which bind to m6A-modified transcripts and mediate the regulation of mRNA fate. Many of these proteins are highly expressed in the germ cells of mammals, and some have been linked to fertility disorders in human patients. In this review, we summarise recent findings on the important roles played by proteins involved in m6A biology in mammalian gametogenesis and fertility. Continued study of the m6A pathway in the mammalian germline will shed further light on the importance of epitranscriptomics in reproduction and may lead to effective treatment of human fertility disorders.

Genetic regulation of m6A RNA methylation and its contribution in human complex diseases

N6-methyladenosine (m6A) has been established as the most prevalent chemical modification in message RNA (mRNA), playing an essential role in determining the fate of RNA molecules. Dysregulation of m6A has been revealed to lead to abnormal physiological conditions and cause various types of human diseases. Recent studies have delineated the genetic regulatory maps for m6A methylation by mapping the quantitative trait loci of m6A (m6A-QTLs), thereby building up the regulatory circuits linking genetic variants, m6A, and human complex traits. Here, we review the recent discoveries concerning the genetic regulatory maps of m6A, describing the methodological and technical details of m6A-QTL identification, and introducing the key findings of the cis- and trans-acting drivers of m6A. We further delve into the tissue- and ethnicity-specificity of m6A-QTL, the association with other molecular phenotypes in light of genetic regulation, the regulators underlying m6A genetics, and importantly, the functional roles of m6A in mediating human complex diseases. Lastly, we discuss potential research avenues that can accelerate the translation of m6A genetics studies toward the development of therapies for human genetic diseases.

Insight into the regulatory mechanism of m6A modification: From MAFLD to hepatocellular carcinoma

In recent years, there has been a significant increase in the incidence of metabolic-associated fatty liver disease (MAFLD), which has been attributed to the increasing prevalence of type 2 diabetes mellitus (T2DM) and obesity. MAFLD affects more than one-third of adults worldwide, making it the most prevalent liver disease globally. Moreover, MAFLD is considered a significant risk factor for hepatocellular carcinoma (HCC), with MAFLD-related HCC cases increasing. Approximately 1 in 6 HCC patients are believed to have MAFLD, and nearly 40 % of these HCC patients do not progress to cirrhosis, indicating direct transformation from MAFLD to HCC. N6-methyladenosine (m6A) is commonly distributed in eukaryotic mRNA and plays a crucial role in normal development and disease progression, particularly in tumors. Numerous studies have highlighted the close association between abnormal m6A modification and cellular metabolic alterations, underscoring its importance in the onset and progression of MAFLD. However, the specific impact of m6A modification on the progression of MAFLD to HCC remains unclear. Can targeting m6A effectively halt the progression of MAFLD-related HCC? In this review, we investigated the pivotal role of abnormal m6A modification in the transition from MAFLD to HCC, explored the potential of m6A modification as a therapeutic target for MAFLD-related HCC, and proposed possible directions for future investigations.

Interplay between N6-adenosine RNA methylation and mRNA splicing

N6-methyladenosine (m6A) is the most abundant modification to mRNAs. Loss-of-function studies of main m6A regulators have indicated the role of m6A in pre-mRNA splicing. Recent studies have reported the role of splicing in preventing m6A deposition. Understanding the interplay between m6A and mRNA splicing holds the potential to clarify the significance of these fundamental molecular mechanisms in cell development and function, thereby shedding light on their involvement in the pathogenesis of myriad diseases.

N6-methyladenosine (m6A) RNA modification in chronic myeloid leukemia: unveiling a novel therapeutic target

N6-methyladenosine (m6A), the most prevalent internal mRNA modification, plays a critical role in physiological processes by regulating gene expression through modulation of mRNA metabolism at multiple stages. In recent years, m6A has garnered significant attention for a deeper understanding of the initiation, progression, and drug resistance of various cancers, including hematological malignancies. Dysregulation of m6A has been implicated in both cancer promotion and suppression. m6A methylation is a complex regulatory process involving methyltransferases (writers), demethylases (erasers), and proteins that recognize specific m6A modifications (readers). This intricate interplay presents challenges for precisely modulating m6A levels, either globally or at specific sites. This review specifically focuses on the role of m6A in chronic myeloid leukemia (CML), a blood cancer characterized by the BCR-ABL1 fusion. We emphasize its impact on leukemia cell survival and drug resistance mechanisms. Notably, inhibitors targeting m6A regulators show promise in preclinical models, suggesting a potential therapeutic avenue for CML. Integrating our understanding of m6A biology with current treatment strategies may lead to more effective therapies, especially for patients with advanced-stage or resistant CML.

Quantitative profiling of m6A at single base resolution across the life cycle of rice and Arabidopsis

N6-methyladenosine (m6A) plays critical roles in regulating mRNA metabolism. However, comprehensive m6A methylomes in different plant tissues with single-base precision have yet to be reported. Here, we present transcriptome-wide m6A maps at single-base resolution in different tissues of rice and Arabidopsis using m6A-SAC-seq. Our analysis uncovers a total of 205,691 m6A sites distributed across 22,574 genes in rice, and 188,282 m6A sites across 19,984 genes in Arabidopsis. The evolutionarily conserved m6A sites in rice and Arabidopsis ortholog gene pairs are involved in controlling tissue development, photosynthesis and stress response. We observe an overall mRNA stabilization effect by 3' UTR m6A sites in certain plant tissues. Like in mammals, a positive correlation between the m6A level and the length of internal exons is also observed in plant mRNA, except for the last exon. Our data suggest an active m6A deposition process occurring near the stop codon in plant mRNA. In addition, the MTA-installed plant mRNA m6A sites correlate with both translation promotion and translation suppression, depicting a more complicated regulatory picture. Our results therefore provide in-depth resources for relating single-base resolution m6A sites with functions in plants and uncover a suppression-activation model controlling m6A biogenesis across species.

m6A regulates heterochromatin in mammalian embryonic stem cells

As the most well-studied modification in mRNA, m6A has been shown to regulate multiple biological processes, including RNA degradation, processing, and translation. Recent studies showed that m6A modification is enriched in chromatin-associated RNAs and nascent RNAs, suggesting m6A might play regulatory roles in chromatin contexts. Indeed, in the past several years, a number of studies have clarified how m6A and its modulators regulate different types of chromatin states. Specifically, in the past 2-3 years, several studies discovered the roles of m6A and/or its modulators in regulating constitutive and facultative heterochromatin, shedding interesting lights on RNA-dependent heterochromatin formation in mammalian cells. This review will summarize and discuss the mechanisms underlying m6A's regulation in different types of heterochromatin, with a specific emphasis on the regulation in mammalian embryonic stem cells, which exhibit distinct features of multiple heterochromatin marks.

N6-Methyladenosine-Modified KREMEN2 Promotes Tumorigenesis and Malignant Progression of High-Grade Serous Ovarian Cancer

High-grade serous ovarian cancer (HGSOC) remains the most lethal female cancer by far. Herein, clinical HGSOC samples had higher N6-methyladenosine (m6A) modification than normal ovarian tissue, and its dysregulation had been reported to drive aberrant transcription and translation programs. However, Kringle-containing transmembrane protein 2 (KREMEN2) and its m6A modification have not been fully elucidated in HGSOC. In this study, the data from the high-throughput messenger RNA (mRNA) sequencing of clinical samples were processed using the weighted correlation network analysis and functional enrichment analysis. Results revealed that KREMEN2 was a driver gene in the tumorigenesis of HGSOC and a potential target of m6A demethylase fat-mass and obesity-associated protein (FTO). KREMEN2 and FTO levels were upregulated and downregulated, respectively, and correlation analysis showed a significant negative correlation in HGSOC samples. Importantly, upregulated KREMEN2 was remarkably associated with lymph node metastasis, distant metastasis, peritoneal metastasis, and high International Federation of Gynecology and Obstetrics stage (Ⅲ/Ⅳ), independent of the age of patients. KREMEN2 promoted the growth of HGSOC in vitro and in vivo, which was dependent on FTO. The methylated RNA immunoprecipitation qPCR and RNA immunoprecipitation assays were performed to verify the m6A level and sites of KREMEN2. FTO overexpression significantly decreased m6A modification in the 3' and 5' untranslated regions of KREMEN2 mRNA and downregulated its expression. In addition, we found that FTO-mediated m6A modification of KREMEN2 mRNA was recognized and stabilized by the m6A reader IGF2BP1 rather than by IGF2BP2 or IGF2BP3. This study highlights the m6A modification of KREMEN2 and extends the importance of RNA epigenetics in HGSOC.

Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

DDX21 mediates co-transcriptional RNA m6A modification to promote transcription termination and genome stability

N6-methyladenosine (m6A) is a crucial RNA modification that regulates diverse biological processes in human cells, but its co-transcriptional deposition and functions remain poorly understood. Here, we identified the RNA helicase DDX21 with a previously unrecognized role in directing m6A modification on nascent RNA for co-transcriptional regulation. DDX21 interacts with METTL3 for co-recruitment to chromatin through its recognition of R-loops, which can be formed co-transcriptionally as nascent transcripts hybridize onto the template DNA strand. Moreover, DDX21's helicase activity is needed for METTL3-mediated m6A deposition onto nascent RNA following recruitment. At transcription termination regions, this nexus of actions promotes XRN2-mediated termination of RNAPII transcription. Disruption of any of these steps, including the loss of DDX21, METTL3, or their enzymatic activities, leads to defective termination that can induce DNA damage. Therefore, we propose that the R-loop-DDX21-METTL3 nexus forges the missing link for co-transcriptional modification of m6A, coordinating transcription termination and genome stability.

Single-molecule epitranscriptomic analysis of full-length HIV-1 RNAs reveals functional roles of site-specific m6As

Although the significance of chemical modifications on RNA is acknowledged, the evolutionary benefits and specific roles in human immunodeficiency virus (HIV-1) replication remain elusive. Most studies have provided only population-averaged values of modifications for fragmented RNAs at low resolution and have relied on indirect analyses of phenotypic effects by perturbing host effectors. Here we analysed chemical modifications on HIV-1 RNAs at the full-length, single RNA level and nucleotide resolution using direct RNA sequencing methods. Our data reveal an unexpectedly simple HIV-1 modification landscape, highlighting three predominant N6-methyladenosine (m6A) modifications near the 3' end. More densely installed in spliced viral messenger RNAs than in genomic RNAs, these m6As play a crucial role in maintaining normal levels of HIV-1 RNA splicing and translation. HIV-1 generates diverse RNA subspecies with distinct m6A ensembles, and maintaining multiple of these m6As on its RNAs provides additional stability and resilience to HIV-1 replication, suggesting an unexplored viral RNA-level evolutionary strategy.

A unique circ_0067716/EIF4A3 double-negative feedback loop impacts malignant transformation of human bronchial epithelial cells induced by benzo(a)pyrene

Benzo(a)pyrene as a pervasive environmental contaminant is characterized by its substantial genotoxicity, and epidemiological investigations have established a correlation between benzo(a)pyrene exposure and the susceptibility to human lung cancer. Notably, much research has focused on the link between epigenetic alterations and lung cancer induced by chemicals, although circRNAs are also emerging as relevant contributors to the carcinogenic process of benzo(a)pyrene. In this study, we identified circ_0067716 as being significantly upregulated in response to stress injury and downregulated during malignant transformation induced by benzo(a)pyrene-7,8-diol-9,10-epoxide (BPDE) in human bronchial epithelial cells. The observed differential expression of circ_0067716 in cells treated with BPDE for varying durations suggests a strong correlation between this circRNA and BPDE exposure. The tissue samples of lung cancer patients also suggest that a lower circ_0067716 expression is associated with BPDE-DNA adduct levels. Remarkably, we demonstrate that EIF4A3, located in the nucleus, interacts with the flanking sequences of circ_0067716 and inhibits its biogenesis. Conversely, circ_0067716 is capable of sequestering EIF4A3 in the cytoplasm, thereby preventing its translocation into the nucleus. EIF4A3 and circ_0067716 can form a double-negative feedback loop that could be affected by BPDE. During the initial phase of BPDE exposure, the expression of circ_0067716 was increased in response to stress injury, resulting in cell apoptosis through the involvement of miR-324-5p/DRAM1/BAX axis. Subsequently, as cellular adaptation progressed, long-term induction due to BPDE exposure led to an elevated EIF4A3 and a reduced circ_0067716 expression, which facilitated the proliferation of cells by stabilizing the PI3K/AKT pathway. Thus, our current study describes the effects of circ_0067716 on the genotoxicity and carcinogenesis induced by benzo(a)pyrene and puts forwards to the possible regulatory mechanism on the occurrence of smoking-related lung cancer, providing a unique insight based on epigenetics.

Determinant of m6A regional preference by transcriptional dynamics

N6-Methyladenosine (m6A) is the most abundant chemical modification occurring on eukaryotic mRNAs, and has been reported to be involved in almost all stages of mRNA metabolism. The distribution of m6A sites is notably asymmetric along mRNAs, with a strong preference toward the 3' terminus of the transcript. How m6A regional preference is shaped remains incompletely understood. In this study, by performing m6A-seq on chromatin-associated RNAs, we found that m6A regional preference arises during transcription. Nucleosome occupancy is remarkedly increased in the region downstream of m6A sites, suggesting an intricate interplay between m6A methylation and nucleosome-mediated transcriptional dynamics. Notably, we found a remarkable slowdown of Pol-II movement around m6A sites. In addition, inhibiting Pol-II movement increases nearby m6A methylation levels. By analyzing massively parallel assays for m6A, we found that RNA secondary structures inhibit m6A methylation. Remarkably, the m6A sites associated with Pol-II pausing tend to be embedded within RNA secondary structures. These results suggest that Pol-II pausing could affect the accessibility of m6A motifs to the methyltransferase complex and subsequent m6A methylation by mediating RNA secondary structure. Overall, our study reveals a crucial role of transcriptional dynamics in the formation of m6A regional preference.

Studying m6A in the brain: a perspective on current methods, challenges, and future directions

A major mechanism of post-transcriptional RNA regulation in cells is the addition of chemical modifications to RNA nucleosides, which contributes to nearly every aspect of the RNA life cycle. N6-methyladenosine (m6A) is a highly prevalent modification in cellular mRNAs and non-coding RNAs, and it plays important roles in the control of gene expression and cellular function. Within the brain, proper regulation of m6A is critical for neurodevelopment, learning and memory, and the response to injury, and m6A dysregulation has been implicated in a variety of neurological disorders. Thus, understanding m6A and how it is regulated in the brain is important for uncovering its roles in brain function and potentially identifying novel therapeutic pathways for human disease. Much of our knowledge of m6A has been driven by technical advances in the ability to map and quantify m6A sites. Here, we review current technologies for characterizing m6A and highlight emerging methods. We discuss the advantages and limitations of current tools as well as major challenges going forward, and we provide our perspective on how continued developments in this area can propel our understanding of m6A in the brain and its role in brain disease.

Challenges to mapping and defining m6A function in viral RNA

Viral RNA molecules contain multiple layers of regulatory information. This includes features beyond the primary sequence, such as RNA structures and RNA modifications, including N6-methyladenosine (m6A). Many recent studies have identified the presence and location of m6A in viral RNA and have found diverse regulatory roles for this modification during viral infection. However, to date, viral m6A mapping strategies have limitations that prevent a complete understanding of the function of m6A on individual viral RNA molecules. While m6A sites have been profiled on bulk RNA from many viruses, the resulting m6A maps of viral RNAs described to date present a composite picture of m6A across viral RNA molecules in the infected cell. Thus, for most viruses, it is unknown if unique viral m6A profiles exist throughout infection, nor if they regulate specific viral life cycle stages. Here, we describe several challenges to defining the function of m6A in viral RNA molecules and provide a framework for future studies to help in the understanding of how m6A regulates viral infection.

The exon junction complex is required for DMD gene splicing fidelity and myogenic differentiation

Deposition of the exon junction complex (EJC) upstream of exon-exon junctions helps maintain transcriptome integrity by preventing spurious re-splicing events in already spliced mRNAs. Here we investigate the importance of EJC for the correct splicing of the 2.2-megabase-long human DMD pre-mRNA, which encodes dystrophin, an essential protein involved in cytoskeletal organization and cell signaling. Using targeted RNA-seq, we show that knock-down of the eIF4A3 and Y14 core components of EJC in a human muscle cell line causes an accumulation of mis-splicing events clustered towards the 3' end of the DMD transcript (Dp427m). This deregulation is conserved in the short Dp71 isoform expressed ubiquitously except in adult skeletal muscle and is rescued with wild-type eIF4A3 and Y14 proteins but not with an EJC assembly-defective mutant eIF4A3. MLN51 protein and EJC-associated ASAP/PSAP complexes independently modulate the inclusion of the regulated exons 71 and 78. Our data confirm the protective role of EJC in maintaining splicing fidelity, which in the DMD gene is necessary to preserve the function of the critical C-terminal protein-protein interaction domain of dystrophin present in all tissue-specific isoforms. Given the role of the EJC in maintaining the integrity of dystrophin, we asked whether the EJC could also be involved in the regulation of a mechanism as complex as skeletal muscle differentiation. We found that eIF4A3 knockdown impairs myogenic differentiation by blocking myotube formation. Collectively, our data provide new insights into the functional roles of EJC in human skeletal muscle.

YTHDC1 cooperates with the THO complex to prevent RNA damage-induced DNA breaks

Certain environmental toxins are nucleic acid damaging agents, as are many chemotherapeutics used for cancer therapy. These agents induce various adducts in DNA as well as RNA. Indeed, most of the nucleic acid adducts (>90%) formed due to these chemicals, such as alkylating agents, occur in RNA 1 . However, compared to the well-studied mechanisms for DNA alkylation repair, the biological consequences of RNA damage are largely unexplored. Here, we demonstrate that RNA damage can directly result in loss of genome integrity. Specifically, we show that a human YTH domain-containing protein, YTHDC1, regulates alkylation damage responses in association with the THO complex (THOC) 2 . In addition to its established binding to N 6-methyladenosine (m6A)-containing RNAs, YTHDC1 binds to N 1-methyladenosine (m1A)-containing RNAs upon alkylation. In the absence of YTHDC1, alkylation damage results in increased alkylation damage sensitivity and DNA breaks. Such phenotypes are fully attributable to RNA damage, since an RNA-specific dealkylase can rescue these phenotypes. These R NA d amage-induced DNA b reaks (RDIBs) depend on R-loop formation, which in turn are processed by factors involved in transcription-coupled nucleotide excision repair. Strikingly, in the absence of YTHDC1 or THOC, an RNA m1A methyltransferase targeted to the nucleus is sufficient to induce DNA breaks. Our results uncover a unique role for YTHDC1-THOC in base damage responses by preventing RDIBs, providing definitive evidence for how damaged RNAs can impact genomic integrity.

mRNA accessibility within mRNPs as a determinant of gene expression

Gene expression is a complex process requiring many control mechanisms to achieve a desired phenotype. DNA accessibility within chromatin is well established as an important determinant of gene expression. By contrast, while mRNA also associates with a complement of proteins, the exact nature of messenger ribonucleoprotein (mRNP) packaging and its functional relevance is not as clear. Recent reports indicate that exon junction complex (EJC)-mediated mRNP packaging renders exon junction-proximal regions inaccessible for m6A methylation, and that EJCs reside within the inaccessible interior of globular transcription and export (TREX) complex-associated nuclear mRNPs. We propose that 'mRNA accessibility' within mRNPs is an important determinant of gene expression that may modulate the specificity of a broad array of regulatory processes including but not limited to m6A methylation.

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.