Eukaryotic RNA Research Articles

N6-Methyladenosine RNA Modification Regulates the Differential Muscle Development in Large White and Ningxiang Pigs.

N6-methyladenosine (m6A) is the most common modification in eukaryotic RNAs. Growing research indicates that m6A methylation is crucial for a multitude of biological processes. However, research on the m6A modifications in the regulation of porcine muscle growth is lacking. In this study, we identified differentially expressed genes in the neonatal period of muscle development between Large White (LW) and NingXiang (NX) pigs and further reported m6A methylation patterns via MeRIP-seq. We found that m6A modification regulates muscle cell development, myofibrils, cell cycle, and phosphatase regulator activity during the neonatal phase of muscle development. Interestingly, differentially expressed genes in LW and NX pigs were mainly enriched in pathways involved in protein synthesis. Furthermore, we performed a conjoint analysis of MeRIP-seq and RNA-seq data and identified 27 differentially expressed and m6A-modified genes. Notably, a typical muscle-specific envelope transmembrane protein, WFS1, was differentially regulated by m6A modifications in LW and NX pigs. We further revealed that the m6A modification accelerated the degradation of WFS1 in a YTHDF2-dependent manner. Noteworthy, we identified a single nucleotide polymorphism (C21551T) within the last exon of WFS1 that resulted in variable m6A methylation, contributing to the differing WFS1 expression levels observed in LW and NX pigs. Our study conducted a comprehensive analysis of the m6A modification on NX and LW pigs during the neonatal period of muscle development, and elucidated the mechanism by which m6A regulates the differential expression of WFS1 in the two breeds.

METTL3 facilitates the progression of cervical cancer by m6A modification-mediated up-regulation of NEK2.

N6-methyladenosine (m6A) is the most prevalent modification found in eukaryotic RNA and played a significant role in various cancers. However, the mechanism by which m6A modification influences cervical cancer (CC) tumorigenesis remains unclear. Therefore, we aim to elucidate the role and mechanism of METTL3 in CC progression. In the present study, we observed a significant upregulation of METTL3 in CC tissues and cell lines. Knockdown of METTL3 resulted in reduced growth, migration, and invasion of CC cells, as well as affected apoptosis, while overexpression of METTL3 reversed these effects. Through a combined analysis of meRIP-seq and Ribo-seq data following METTL3 knockdown, NEK2 was identified as a key target of METTL3 in CC cells. Correlation analysis, MeRIP-qPCR, and luciferase reporter assay suggested that METTL3 regulates NEK2 expression through m6A modification. NEK2 synergized with METTL3 to mediate the malignant phenotype of CC cells. The METTL3-NEK2 axis promoted CC progression by activating the Wnt/β-catenin pathway and inhibiting the apoptosis pathway. In conclusion, METTL3 facilitated the malignant progression of CC and contributed to the formation of the METTL3-NEK2 regulatory axis in an m6A-dependent manner, which represented a potential target for CC therapy.

Interactive Structural Analysis of KH3-4 Didomains of IGF2BPs with Preferred RNA Motif Having m6A Through Dynamics Simulation Studies.

m6A modification is the most common internal modification of messenger RNA in eukaryotes, and the disorder of m6A can trigger cancer progression. The GGACU is considered the most frequent consensus sequence of target transcripts which have a GGAC m6A core motif. Newly identified m6A 'readers' insulin-like growth factor 2 mRNA-binding proteins modulate gene expression by binding to the m6A binding sites of target mRNAs, thereby affecting various cancer-related processes. The dynamic impact of the methylation at m6A within the GGAC motif on human IGF2BPs has not been investigated at the structural level. In this study, through in silico analysis, we mapped IGF2BPs binding sites for the GGm6AC RNA core motif of target mRNAs. Subsequent molecular dynamics simulation analysis at 400 ns revealed that only the KH4 domain of IGF2BP1, containing the 503GKGG506 motif and its periphery residues, was involved in the interaction with the GGm6AC backbone. Meanwhile, the methyl group of m6A is accommodated by a shallow hydrophobic cradle formed by hydrophobic residues. Interestingly, in IGF2BP2 and IGF2BP3 complexes, the RNA was observed to shift from the KH4 domain to the KH3 domain in the simulation at 400 ns, indicating a distinct dynamic behavior. This suggests a conformational stabilization upon binding, likely essential for the functional interactions involving the KH3-4 domains. These findings highlight the potential of targeting IGF2BPs' interactions with m6A modifications for the development of novel oncological therapies.

Leveraging HILIC/ERLIC Separations for Online Nanoscale LC-MS/MS Analysis of Phosphopeptide Isoforms from RNA Polymerase II C-terminal Domain.

The eukaryotic RNA polymerase II (Pol II) multi-protein complex transcribes mRNA and coordinates several steps of co-transcriptional mRNA processing and chromatin modification. The largest Pol II subunit, Rpb1, has a C-terminal domain (CTD) comprising dozens of repeated heptad sequences (Tyr1-Ser2-Pro3-Thr4-Ser5-Pro6-Ser7), each containing five phospho-accepting amino acids. The CTD heptads are dynamically phosphorylated, creating specific patterns correlated with steps of transcription initiation, elongation, and termination. This CTD phosphorylation 'code' choreographs dynamic recruitment of important co-regulatory proteins during gene transcription. Genetic tools were used to engineer protease cleavage sites across the CTD (msCTD), creating tryptic peptides with unique sequences amenable to mass spectrometry analysis. However, phosphorylation isoforms within each msCTD sequence are difficult to resolve by standard reversed phase chromatography typically used for LC-MS/MS applications. Here, we use a panel of synthetic CTD phosphopeptides to explore the potential of hydrophilic interaction and electrostatic repulsion hydrophilic interaction (HILIC and ERLIC) chromatography as alternatives to reversed phase separation for CTD phosphopeptide analysis. Our results demonstrate that ERLIC provides improved performance for separation of singly- and doubly-phosphorylated CTD peptides for sequence analysis by LC-MS/MS. Analysis of native yeast msCTD confirms that phosphorylation on Ser5 and Ser2 represents the major endogenous phosphoisoforms. We expect this methodology will be especially useful in the investigation of pathways where multiple protein phosphorylation events converge in close proximity.

Self-assembly and condensation of intermolecular poly(UG) RNA quadruplexes.

Poly(UG) or 'pUG' dinucleotide repeats are highly abundant sequences in eukaryotic RNAs. In Caenorhabditis elegans, pUGs are added to RNA 3' ends to direct gene silencing within Mutator foci, a germ granule condensate. Here, we show that pUG RNAs efficiently self-assemble into gel condensates through quadruplex (G4) interactions. Short pUG sequences form right-handed intermolecular G4s (pUG G4s), while longer pUGs form left-handed intramolecular G4s (pUG folds). We determined a 1.05 Å crystal structure of an intermolecular pUG G4, which reveals an eight stranded G4 dimer involving 48 nucleotides, 7 different G and U quartet conformations, 7 coordinated potassium ions, 8 sodium ions and a buried water molecule. A comparison of the intermolecular pUG G4 and intramolecular pUG fold structures provides insights into the molecular basis for G4 handedness and illustrates how a simple dinucleotide repeat sequence can form complex structures with diverse topologies.

Establishment of closed 35S ribosomal RNA gene chromatin in stationary Saccharomyces cerevisiae cells.

As a first step in eukaryotic ribosome biogenesis RNA polymerase (Pol) I synthesizes a large ribosomal RNA (rRNA) precursor from multicopy rRNA gene loci. This process is essential for cellular growth and regulated in response to the cell's physiological state. rRNA gene transcription is downregulated upon growth to stationary phase in the yeast Saccharomyces cerevisiae. This reduction correlates with characteristic changes in rRNA gene chromatin structure from a transcriptionally active 'open' state to a non-transcribed 'closed' state. The conserved lysine deacetylase Rpd3 was shown to be required for this chromatin transition. We found that Rpd3 is needed for tight repression of Pol I transcription upon growth to stationary phase as a prerequisite for the establishment of the closed chromatin state. We provide evidence that Rpd3 regulates Pol I transcription by adjusting cellular levels of the Pol I preinitiation complex component core factor (CF). Importantly, our study identifies CF as the complex limiting the number of open rRNA genes in exponentially growing and stationary cells.

The role of structure in regulatory RNA elements.

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e. RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Genetic Variants in METTL16 Affect the Risk of Non-Syndromic Orofacial Clefts.

N6-methyladenosine (m6A) is the most prevalent modification of RNA in eukaryotes which is associated with many cellular processes and diseases. Here, our objective is to explore whether genetic variants in m6A modification genes are associated with the risk of non-syndrome orofacial clefts (NSOCs). The transmission disequilibrium test (TDT) was performed to calculate the association between single nucleotide polymorphisms (SNPs) in m6A modification genes and NSOCs risk in 944 case-parent trios. The function of SNP was predicted by HaploReg, RegulomeDB and histone enrichment data. The expression quantitative trait locus (eQTL) analysis was examined using Genotype-Tissue Expression (GTEx) and eQTLGen. The role of gene in the development of NSOCs was assessed with correlation and enrichment analysis based on gene expression data in mice craniofacial tissue and zebrafish embryo. We identified that rs8078195 (A > C) in METTL16 was suggestively associated with the increased risk of NSOCs (OR = 1.32, p = 1.80E - 03). The region surrounding rs8078195 was subjected to deoxyribonuclease hypersensitivity and enriched with multiple histone modifications. In addition, it had a significant eQTL effect with METTL16 in skin tissue and human peripheral blood, which played an important role in NSOCs development. Bioinformatic analysis indicated that METTL16 contributed to the development of NSOCs probably by regulating cell cycle process. Rs8078195 in METTL16 was associated with the occurrence of NSOCs.

ZC3H14 facilitates backsplicing by binding to exon-intron boundary and 3' UTR.

Circular RNAs (circRNAs) are natural outputs of eukaryotic transcription and RNA processing and have emerged as critical regulators in physiology and diseases. Although multiple cis-elements and trans-factors are reported to modulate the backsplicing of circRNA biogenesis, most of these regulations play roles in flanking introns of circRNAs. Here, using a genome-wide CRISPR knockout screen, we have identified an evolutionarily conserved RNA-binding protein ZC3H14 in regulating circRNA biogenesis. ZC3H14 binds to 3' and 5' exon-intron boundaries and 3' UTRs of cognate mRNAs to promote circRNA biogenesis through dimerization and the association with spliceosome. Yeast knockout of the ZC3H14 ortholog Nab2 has significantly lower levels of circRNAs. Zc3h14-/- mice exhibit disrupted spermatogenesis and reduced testicular circRNA levels. Additionally, expression levels of human ZC3H14 are associated with non-obstructive azoospermia. Our findings reveal a conserved requirement for ZC3H14 in the modulation of backsplicing and link ZC3H14 and circRNA biogenesis to male fertility.

Biochemical characterization of Mycobacterial RNA polymerases.

Tuberculosis is caused by the bacterium Mycobacterium tuberculosis (Mtb). While eukaryotic species employ several specialized RNA polymerases (Pols) to fulfill the RNA synthesis requirements of the cell, bacterial species use a single RNA polymerase (RNAP). To contribute to the foundational understanding of how Mtb and the related non-pathogenic mycobacterial species, Mycobacterium smegmatis (Msm), perform the essential function of RNA synthesis, we performed a series of in vitro transcription experiments to define the unique enzymatic properties of Mtb and Msm RNAPs. In this study, we characterize the mechanism of nucleotide addition used by these bacterial RNAPs with comparisons to previously characterized eukaryotic Pols I, II, and III. We show that Mtb RNAP and Msm RNAP demonstrate similar enzymatic properties and nucleotide addition kinetics to each other but diverge significantly from eukaryotic Pols. We also show that Mtb RNAP and Msm RNAP uniquely bind a nucleotide analog with significantly higher affinity than canonical nucleotides, in contrast to eukaryotic RNA polymerase II. This affinity for analogs may reveal a vulnerability for selective inhibition of the pathogenic bacterial enzyme.IMPORTANCETuberculosis, caused by the bacterium Mycobacterium tuberculosis (Mtb), remains a severe global health threat. The World Health Organization (WHO) has reported that tuberculosis is second only to COVID-19 as the most lethal infection worldwide, with more annual deaths than HIV and AIDS (WHO.int). The first-line treatment for tuberculosis, Rifampin (or Rifampicin), specifically targets the Mtb RNA polymerase. This drug has been used for decades, leading to increased numbers of multi-drug-resistant infections (Stephanie, et al). To effectively treat tuberculosis, there is an urgent need for new therapeutics that selectively target vulnerabilities of the bacteria and not the host. Characterization of the differences between Mtb enzymes and host enzymes is critical to inform these ongoing drug design efforts.

African swine fever virus RNA polymerase subunits C315R and H359L inhibition host translation by activating the PKR-eIF2a pathway and suppression inflammatory responses.

ASFV C315R is homologous to the transcription factor TFIIB of large unclassified DNA viruses, and H359L is identical to the subunit 3 (RPB3) of eukaryotic RNA polymerase II. The C315R and H359L may play an important role in ASFV replication and transcription. Here, we evaluated the biological function of the C315R and H359L genes during virus replication in vitro and during infection in pigs. Results showed that C315R and H359L are highly conserved among ASFV genotype II strains; quantitative PCR (qPCR) and western blotting analyses revealed that C315R and H359L are early transcribed genes prior to viral DNA replication, but their protein expression is delayed. The immunofluorescence and western blotting analysis revealed that both proteins localized in the cell cytoplasm and nucleus at 24 h post infection, however, pH359L was mainly detected in the cell cytoplasm. Furthermore, overexpression of pH359L in MA104 cells significantly increased viral titer, RNA transcription levels, and viral protein expression levels, while overexpression of pC315R slightly enhanced ASFV replication. In contrast, siRNA targeting ASFV-H359L or C315R reduced replication efficiency in porcine macrophage culture compared to the parent ASFV-CN/SC/2019, demonstrating that C315R and H359L genes are necessary for ASFV replication. Finally, the functional role of C315R or H359L on PKR and eIF2α phosphorylation status and SG formation, as well as cytokine production were evaluated. These studies demonstrated that C315R and H359L are involved in virus replication processes in swine and play important roles in ASFV replication.

Regulatory effect of N6-methyladenosine on tumor angiogenesis.

Previous studies have demonstrated that genetic alterations governing epigenetic processes frequently drive tumor development and that modifications in RNA may contribute to these alterations. In the 1970s, researchers discovered that N6-methyladenosine (m6A) is the most prevalent form of RNA modification in advanced eukaryotic messenger RNA (mRNA) and noncoding RNA (ncRNA). This modification is involved in nearly all stages of the RNA life cycle. M6A modification is regulated by enzymes known as m6A methyltransferases (writers) and demethylases (erasers). Numerous studies have indicated that m6A modification can impact cancer progression by regulating cancer-related biological functions. Tumor angiogenesis, an important and unregulated process, plays a pivotal role in tumor initiation, growth, and metastasis. The interaction between m6A and ncRNAs is widely recognized as a significant factor in proliferation and angiogenesis. Therefore, this article provides a comprehensive review of the regulatory mechanisms underlying m6A RNA modifications and ncRNAs in tumor angiogenesis, as well as the latest advancements in molecular targeted therapy. The aim of this study is to offer novel insights for clinical tumor therapy.

Asgard archaea defense systems and their roles in the origin of eukaryotic immunity

Dozens of new antiviral systems have been recently characterized in bacteria. Some of these systems are present in eukaryotes and appear to have originated in prokaryotes, but little is known about these defense mechanisms in archaea. Here, we explore the diversity and distribution of defense systems in archaea and identify 2610 complete systems in Asgardarchaeota, a group of archaea related to eukaryotes. The Asgard defense systems comprise 89 unique systems, including argonaute, NLR, Mokosh, viperin, Lassamu, and CBASS. Asgard viperin and argonaute proteins have structural homology to eukaryotic proteins, and phylogenetic analyses suggest that eukaryotic viperin proteins were derived from Asgard viperins. We show that Asgard viperins display anti-phage activity when heterologously expressed in bacteria. Eukaryotic and bacterial argonaute proteins appear to have originated in Asgardarchaeota, and Asgard argonaute proteins have argonaute-PIWI domains, key components of eukaryotic RNA interference systems. Our results support that Asgardarchaeota played important roles in the origin of antiviral defense systems in eukaryotes.

Research Progress on the Role of M6A in Regulating Economic Traits in Livestock.

Reversible regulation of N6-methyladenosine (m6A) methylation of eukaryotic RNA via methyltransferases is an important epigenetic event affecting RNA metabolism. As such, m6A methylation plays crucial roles in regulating animal growth, development, reproduction, and disease progression. Herein, we review the latest research advancements in m6A methylation modifications and discuss regulatory aspects in the context of growth, development, and reproductive traits of livestock. New insights are highlighted and perspectives for the study of m6A methylation modifications in shaping economically important traits are discussed.

Selective recognition of PTRE1 transcripts mediated by protein–protein interaction between the m6A reader ECT2 and PTRE1

N6-methyladenosine (m6A) is a prevalent internal post-transcriptional modification in eukaryotic RNAs executed by m6A-binding proteins known as “readers.” Our previous research demonstrated that the Arabidopsis m6A reader ECT2 positively regulates transcript levels of the proteasome regulator PTRE1 and several 20S proteasome subunits, thereby enhancing 26S proteasome activity. However, mechanism underlying the selective recognition of m6A targets by readers, such as ECT2, remains elusive. In this study, we further demonstrate that ECT2 physically interacts with PTRE1 and several 20S proteasome subunits. This interaction, which occurs on the ribosome, involves the N terminus of PTRE1, suggesting that ECT2 might bind to the nascent PTRE1 polypeptide. Deleting ECT2’s protein interaction domain impairs its mRNA-binding ability, whereas mutations in the m6A-RNA-binding site do not affect protein–protein interactions. Moreover, introducing a novel protein-binding domain into ECT2 increases transcript levels of proteins interacting with this domain. Our findings indicate that interaction with the PTRE1 protein enhances ECT2’s binding to PTRE1 m6A mRNAs during translation, thereby regulating PTRE1 mRNA levels.

Transcriptome-wide m6A methylation profile reveals its potential role underlying drought response in wheat (Triticum aestivum L.).

This study revealed the transcriptome-wide m6A methylation profile under drought stress and found that TaETC9 might regulate drought tolerance through mediating RNA methylation in wheat. Drought is one of the most destructive environmental constraints limiting crop growth and development. N6-methyladenosine (m6A) is a prevalent and important post-transcriptional modification in various eukaryotic RNA molecules, playing the crucial role in regulating drought response in plants. However, the significance of m6A in wheat (Triticum aestivum L.), particularly its involvment in drought response, remains underexplored. In this study, we investigated the transcriptome-wide m6A profile under drought stress using parallel m6A immunoprecipitation sequencing (MeRIP-seq). Totally, 4221 m6A peaks in 3733 m6A-modified genes were obtained, of which 373 methylated peaks exhibited differential expression between the control (CK) and drought-stressed treatments. These m6A loci were significantly enriched in proximity to stop codons and within the 3'-untranslated region. Integration of MeRIP-seq and RNA-seq revealed a positive correlation between m6A methylation and mRNA abundance and the genes displaying both differential methylation and expression were obtained. Finally, qRT-PCR analyses were further performed and the results found that the m6A-binding protein (TaETC9) showed significant up-regulation, while the m6A demethylase (TaALKBH10B) was significantly down-regulated under drought stress, contributing to increased m6A levels. Furthermore, the loss-of-function mutant of TaECT9 displayed significantly higher drought sensitivity compared to the wild type, highlighting its role in regulating drought tolerance. This study reported the first wheat m6A profile associated with drought stress, laying the groundwork for unraveling the potential role of RNA methylation in drought responses and enhancing stress tolerance in wheat through epigenetic approaches.

Functions and mechanisms of RNA m6A regulators in breast cancer (Review).

Breast cancer (BC) is a major malignant tumor in females and the incidence rate of BC has increased worldwide in recent years. N6‑methyladenosine (m6A) is a methylation modification that occurs extensively in eukaryotic RNA. The abnormal expression of m6A and related regulatory proteins can activate or inhibit certain signal pathways or oncogenes, thus affecting the proliferation, metastasis and prognosis of BC. Numerous studies have shown that m6A regulator disorder exists in BC, and this disorder can be reversed. Therefore, m6A is predicted as a potential therapeutic target for BC. However, the molecular mechanism of m6A RNA methylation regulating the occurrence and development of BC has not been comprehensively elucidated. In this review article, the functions of various m6A regulators and the specific mechanisms of certain regulators of the progress of BC were summarized. Furthermore, the dual role of RNA methylation in tumor progression was discussed, concluding that RNA methylation can not only lead to tumorigenesis but at times give rise to inhibition of tumor formation. In addition, further comprehensive analysis on mechanisms of m6A regulators in BC is conducive to screening effective potential targets and formulating targeted treatment strategies, which will provide new methods for the prevention and treatment of BC.

A nuclear RNA degradation code for eukaryotic transcriptome surveillance.

The RNA exosome plays critical roles in eukaryotic RNA degradation, but it remains unclear how the exosome specifically recognizes its targets. The PAXT connection is an adaptor that recruits the exosome to polyadenylated RNAs in the nucleus, especially transcripts polyadenylated at intronic poly(A) sites. Here we show that PAXT-mediated RNA degradation is induced by the combination of a 5' splice site and a poly(A) junction, but not by either sequence alone. These sequences are bound by U1 snRNP and cleavage/polyadenylation factors, which in turn cooperatively recruit PAXT. As the 5' splice site-poly(A) junction combination is typically not found on correctly processed full-length RNAs, we propose that it functions as a "nuclear RNA degradation code" (NRDC). Importantly, disease-associated single nucleotide polymorphisms that create novel 5' splice sites in 3' untranslated regions can induce aberrant mRNA degradation via the NRDC mechanism. Together our study identified the first NRDC, revealed its recognition mechanism, and characterized its role in human diseases.

The Beak of Eukaryotic Ribosomes: Life, Work and Miracles.

Ribosomes are not totally globular machines. Instead, they comprise prominent structural protrusions and a myriad of tentacle-like projections, which are frequently made up of ribosomal RNA expansion segments and N- or C-terminal extensions of ribosomal proteins. This is more evident in higher eukaryotic ribosomes. One of the most characteristic protrusions, present in small ribosomal subunits in all three domains of life, is the so-called beak, which is relevant for the function and regulation of the ribosome's activities. During evolution, the beak has transitioned from an all ribosomal RNA structure (helix h33 in 16S rRNA) in bacteria, to an arrangement formed by three ribosomal proteins, eS10, eS12 and eS31, and a smaller h33 ribosomal RNA in eukaryotes. In this review, we describe the different structural and functional properties of the eukaryotic beak. We discuss the state-of-the-art concerning its composition and functional significance, including other processes apparently not related to translation, and the dynamics of its assembly in yeast and human cells. Moreover, we outline the current view about the relevance of the beak's components in human diseases, especially in ribosomopathies and cancer.

P-127 The RNA m6A landscape of mouse oocytes and preimplantation embryos

Abstract Study question Despite the significance of N6-methyladenosine (m6A) in gene regulation, the requirement for large amounts of RNA has hindered m6A profiling in mammalian early embryos. Summary answer We apply picoMeRIP-seq to map m6A in mouse oocytes and preimplantation embryos. We define the landscape of m6A during the maternal-to-zygotic transition. What is known already N6-methyladenosine (m6A), the most prevalent internal modification in eukaryotic messenger RNA (mRNA), plays key regulatory roles in many biological processes (for example, RNA stability, splicing, transport and translation) and is involved in a variety of physiological processes (for example, cell differentiation and reprogramming, embryonic development and stress responses). Study design, size, duration We recently developed an m6A mapping assay for low RNA input, low-input methyl RNA immunoprecipitation and sequencing (picoMeRIP-seq), which enables mapping of the m6A methylome using a low number of cells. In this article, we applied picoMeRIP-seq to multiple developmental stages of mouse oocytes and early embryos to profile their m6A landscapes. Participants/materials, setting, methods We use mouse as study materials.Using picoMeRIP-seq, we generated transcriptome-wide m6A maps for six stages (GV, MII, zygote, two-cell, eight-cell and blastocyst), with two biological replicates per stage and 49–110 oocytes/embryos per replicate. Main results and the role of chance We define the landscape of m6A during the maternal-to-zygotic transition, including stage-specifically expressed transcription factors essential for cell fate determination. Both the maternally inherited transcripts to be degraded post fertilization and the zygotically activated genes during zygotic genome activation are widely marked by m6A. In contrast to m6A-marked zygotic ally-activated genes, m6A-marked maternally inherited transcripts have a higher tendency to be targeted by microRNAs. Moreover, RNAs derived from retrotransposons, such as MTA that is maternally expressed and MERVL that is transcriptionally activated at the two-cell stage, are largely marked by m6A. Our results provide a foundation for future studies exploring the regulatory roles of m6A in mammalian early embryonic development. Limitations, reasons for caution The effect of m6A on the turnover of retrotransposon-derived transcripts abundant in oocytes and early embryos needs to be further addressed in future studies. Wider implications of the findings Our results provide a foundation for future studies exploring the regulatory roles of m6A in mammalian early embryonic development. Trial registration number not applicable