Novel RNA regulatory mechanisms revealed in the epitranscriptome

Abstract

Methyl-6-adenosine (m⁶A) has been hypothesized to exist since the 1970s,¹ but little has been known about the specific RNAs, or sites within them, that are affected by this RNA modification. Here, we report that recent work has shown RNA modifications like m⁶A, collectively called the “epitranscriptome,” are a pervasive feature of mammalian cells and likely play a role in development and disease. An enrichment of m⁶A near the last CDS of thousands of genes has implicated m⁶A in transcript processing, translational regulation and potentially a mechanism for regulating miRNA maturation. Also, because the sites of m⁶A show strong evolutionary conservation and have been replicated in nearly identical sites between mouse and human, strong evolutionary pressures are likely being maintained for this mark.^2,3 Finally, we note that m⁶A is one of over 100 modifications of RNA that have been reported,⁴ and with the combination of high-throughput, next-generation sequencing (NGS) techniques, immunoprecipitation with appropriate antibodies and splicing-aware peak-finding, the dynamics of the epitranscriptome can now be mapped and characterized to discern their specific cellular roles.

Keywords: :

• RNA modifications
• m⁶A
• N6-methyladenosine
• methyladenosine
• epigenetics
• epigenomics
• epitranscriptome

Background

Dynamic regulation of DNA bases in epigenetics has been established since 1925, and it is known to play a critical role in embryogenesis, imprinting, development and cancer.⁵ Indeed, incomplete reprogramming of epigenetic marks from cells transformed during somatic cell nuclear transfer (SCNT) is a factor in the inability of these cells to be completely pluripotent.⁶ The role of epigenetics is so important to biology that some have referred to these bases as the extra fifth and sixth genetic bases of the genetic code.^7,8 However, RNA modifications, collectively referred to as “epitranscritpomes,” are a much less-studied area. This is partially due to the fact that RNA degrades quickly and that there are no easy methods to find modified RNA bases, whereas de-amination methods, such as bi-sulfite treatment, work very well in DNA.

Early estimates of m⁶A frequency were in the range of 0.1–0.4%,^9-11 but they required significant effort and work. For example, one of the first studies of m⁶A used the adenovirus type 2 (Ad-2) system,¹² in which viral mRNAs are generated from a large common nuclear precursor RNA in the nucleus. The abundance of viral RNAs transcribed from multiple copies of the virons made it possible to investigate the timing of N-6-methylation on adenosine and the flow of m⁶A from nuclear precursor RNA to mature cytoplasmic mRNA by kinetic studies. But, this required 3H-methyl-labeled methionine and ¹⁴C-labeled uridine, which then were pulsed in intervals and RNA was measured by quantitative RNA hybridization. A comparison of the nuclear dwell time of precursor RNA compared with the mRNA levels in the cytoplasm showed that the ratio of 3H to ¹⁴C was significantly higher in cytoplasmic RNA than in nuclear RNA, indicating that m⁶A was added to transcripts before splicing. These preliminary results implicated a role for m⁶A in RNA processing or stabilization, or both, during the biogenesis of mRNA.

Even though early papers suggested m⁶A as a fundamental aspect of RNA biology, little work was done in the field for several decades. This is partly because the known methods to identify m⁶A sites were limited to a pre-ordained focus on a single transcript, required radio-labeled bases and it was thought that RNA modifications were primarily in rRNA¹³ and tRNA.¹⁴ However, in 2011, the FTO gene (associated with fat mass and obesity) was discovered as the enzyme responsible for the oxidative demethylation of m⁶A.¹⁵ The FTO gene is a member of the 2-oxoglutarate (2OG)-dependent oxygenases that can catalyze the de-methylation of RNA and DNA.^16,17 Related to FTO is METTL3, a 70-kDA subunit that forms part of the N6-adenosine methyltransferase, and it is thought that METTL3 is responsible for methylating regular adenosines into m⁶A.¹⁸ Notably, additional METTL3 homologs exist in the human and mouse genomes and could also function as RNA adenosine methyltransferases, or in concert with other RNA-binding proteins. For example, the ALKBH5 gene was recently found to demethylate m⁶A in vitro and in vivo, and knockdown of ALKBH5 not only resulted in higher m⁶A levels but also altered mRNA export levels and reduced fertility in mice.¹⁹ The characterization and discovery of the enzymes responsible for m⁶A sparked a renewed interest in this particular RNA modification.

Global Methods for Epitranscriptomics

The emergence of a new m⁶A-specific antibody in 2011 enabled researchers to perform ChIP-seq-like experiments to identify potential m⁶A sites throughout the transcriptome, also called the “RNA methylome,” which offers a glimpse into the full epitranscriptome.

In 2012, two global methods were independently developed by Meyer et al. (MeRIP-Seq)² and by Dominissini et al. (m⁶A-Seq)³ demonstrating that immunoprecipitation and next-generation sequencing (NGS) can be used to efficiently and robustly localize m⁶A sites in RNA.

In MeRIP-Seq and m⁶A-Seq, fragmented RNA is immunoprecipitated using an m6A-specific antibody. There are currently two available antibodies, from Synaptic Systems and New England Biolabs (NEB), and both appear to have equal performance.²⁰ The IP RNA fragments are then sequenced, along with the non-IP sample (control), to be analyzed using peak-detection software. One publically available tool, MeRIPPeR (www.sourceforge.net/projects/meripper/), maps reads to windows across the genome and uses Fisher’s exact test to find statistically significant enrichment of reads in the IP vs. the control samples. Windows are then merged into possible peak locations for m⁶A sites. Both data sets found an enrichment of sites near the stop codon and the beginning of the 3′ UTR, and some evidence was reported of an enrichment of m⁶A near transcription start site (TSS).^2,3,20

While IP-based methods have been successful so far, they are limited insofar as they cannot directly reveal the precise adenosine within a peak, and other methods are now being developed to advance the field. In particular, there are two proof-of-principle studies that have shown that m⁶A can also be directly detected on the native molecules. The first uses a single-molecule sequencing approach to observe a HIV-based reverse transcriptase as it makes cDNA.²⁰ As bases are incorporated to form the complementary strand, a kinetic signature of each base context can be observed. It has been shown that the steric change of the single methyl group in m⁶A is large enough to observe a significantly different kinetic signature of base incorporation between m⁶A vs. A. Second, this same principle can be applied to any single-molecule monitoring technology, and direct RNA sequencing for RNA within a protein nanopore has also been reported.²¹ As these technologies improve, global methods for m⁶A detection will likely all become accurate to the exact base, which will help in the analysis and contextualization of these data.

Distribution of m⁶A across the genome

While MeRIP-Seq and m⁶A-Seq have mapped putative m⁶A sites across the transcriptome, their precise biological function and the mechanism of m⁶A site-specificity is not yet known. The enrichment of m⁶A sites near the stop codon and in the 3′ UTR suggest that it may be involved in transcription or translation regulation, but experimental evidence of this interaction has not yet been generated. Similarly, the discovery of the FTO, ALKBH5 and METTL3 genes provide insight into how m⁶A sites are directly manipulated, but these proteins are likely part of a larger network that facilitates dynamic changes in m⁶A sites in response to environmental changes and cellular stimuli. For example, m⁶A levels differ between tissues (as demonstrated by immunoblot experiments) with higher levels in the brain, liver and kidney cells, while lower in the heart and lung.² Also, m⁶A levels have been shown to increase with neuronal development in mouse tissues, which implies that m⁶A sites and levels have a defined spatio-temporal trajectory,² and current evidence shows that m⁶A levels are more important as an animal ages.

All existing data show that m⁶A sites cover the majority of the genome (Fig. 1), but there is strong evidence that the sites are not randomly distributed. The enrichment of the peaks near the stop codon and 5′ end of the 3′ UTR has been replicated in both studies,^2,3 and the discovered sites have exact-matching orthologous regions found in both human and mouse genes. Also, an analysis of the specific sequences that underlie the m⁶A sites show very high PhyloP conservation scores,² which indicates strong evolutionary pressures at these sites. Taken together, these data indicate that m⁶A sites likely play a key role in the RNA biology and cellular functions of these transcripts, and that their functions are constrained over millions of years of evolution.

Figure 1. Sites of m6A are pervasive across the genome. This shows a circos plot of all m6A peaks for humans, with HEK293T cell peaks from Meyer et al. (in red, inner circle) and HepG2 from Dominissini et al. (in blue, inner-most circle). Peaks were found in over 10,000 genes and in all chromosomes (outer sections, varying colors).

Six Putative Cellular Functions of m⁶A

Many new questions about RNA biology have been raised by the discovery and validation of this epitranscriptomic modification. While m⁶A has been localized in thousands of transcripts, its biological function and role are still largely a mystery. Thus far, studies have shown that m⁶A sites can be found in both mature and immature transcripts, and older data indicate that the m⁶A may be placed on transcripts before splicing or other RNA processing steps. It is also notable many peaks in the MeRIP papers^2,3 were found in introns, non poly-adenylated transcripts and non-coding RNAs, providing further evidence that at least some of the m⁶A machinery is working in the early stages of RNA processing.

RNA stability

Accordingly, one likely major role for m⁶A could be in the stabilization of transcripts before they are processed. This hypothesis is supported by work examining methylation patterns in plant rRNA, which showed a stabilizing effect for ribosomal function in response to environmental changes such as temperature.²² Thus, m⁶A could have a similar function in stabilizing mRNA transcripts and ncRNAs. Studies of m⁶A in plant mRNAs have further demonstrated the necessity of this modification and implicated its role in genes involved in transport and stress and stimulus response.²³

Splicing regulator

Another role for m⁶A could serve to regulate alternative splicing. For example, the mRNA adenosine methylate (MTA) homolog in plants is proposed to interact with FIP37, which is implicated in RNA splicing.²⁴ Moreover, m⁶A has previously been shown to be involved in splicing, particularly in alternatively spliced (AS) exons, implicating m⁶A in the regulation of AS events.³ Interestingly, previous work has shown that AS exons evolve faster than constitutively spliced (CS) exons, and they enable an “evolutionary playground” for genes to modify their minor transcript isoforms while keeping the major isoform present in cells.²⁵ However, since m⁶A sites show overall strong conservation,² it is likely that m⁶A conservation may be contingent upon its genomic location, such that m⁶A site conservation could be relaxed in fast-evolving genes with complex splicing patterns. More tissues, and more time points, will be needed to examine the distribution of different spliced forms of m⁶A-enriched transcripts.

miRNA Attenuator

Given the proximity to the stop codon and the 5′ end of the 3′UTR, m⁶A may also interact with miRNA sites in the 3′UTR.² It is notable that 67% of UTRs with an m⁶A peak have a predicted miRNA-binding site, yet only 30% of genes in the genome have predicted miRNA-binding sites,² representing a 2.2-fold enrichment. Moreover, preliminary data have shown an inverse spatial relationship of miRNA-binding sites and m⁶A sites,² even though the m⁶A-containing transcripts have higher expression of their cognate miRNA targets. Interestingly, some reports have shown that methylated RNA can inhibit miRNA maturation, giving evidence for a direct connection of the methylation of miRNA and miRNA function,²⁶ whereby methylation serves as a negative regulator of the ability of the miRNA to function on the transcript. However, it is not yet known if the methylation of the mRNA is related to the methylation of the miRNA, or if the same enzymes are responsible for controlling methylation on both small and large RNAs. To explain these nascent data, further work is needed to (1) characterize the myriad RNA binding proteins that are predicted and their targets, (2) discern which of the 107 known RNA modifications are present on a target,²⁰ including m⁶A and (3) examine the molecular dynamics and binding rate changes associated with these ribonucleotide modifications. Depending on these data and where modifications occur in transcripts, both antagonistic and synergistic models of epitranscriptomic control of miRNA function are possible.

RNA editing prevention

Another potential regulatory role for m⁶A is in the post-transcriptional modification of A→I editing, which results from the deamination of the N⁶ in adenosine by the ADAR family of adenosine deaminases. Intriguingly, m⁶A is known to block the activity of the ADAR enzymes that perform A→I conversion in mRNA.²⁷ Inosine is treated as guanosine in translation and when forming RNA secondary structures,²⁸ so m⁶A may serve as a mechanism to maintain adenosines in transcripts and hinder RNA editing. However, the challenge is both in finding existing A→I editing sites and seeing a dynamic change induced by an introduction of an m⁶A site. So far, there is no m⁶A peak that overlaps with any known RNA editing site, which gives some support for this hypothesis (Fig. 2).

Figure 2. Effects of m6A on RNA function. We hypothesize that the sites of m6A (blue) will prevent the N-6 de-amination that occurs in RNA editing (red), where adenosine (A) gets converted into inosine (I) then guanosine (G). These sites may then have many roles in the function of RNA, from splicing to translation changes (blue).

Disease accelerators

Variants of the FTO allele have been associated with obesity and increased body fat,^29,30 where as other mutations in the FTO gene have resulted in lower frontal lobe volumes and are also associated with Alzheimer disease.^31-33 Rhein (cassic acid) was recently discovered to be an inhibitor for FTO, binding to its active site in vitro and preventing demethylation of m⁶A, which may help to study FTO and m⁶A levels and their association with diseases.³⁴ Also, m⁶A levels have been observed with substantial variation in different cancer cell lines² and affect TP53 function,¹⁹ meaning that they may play a role in cancer as well. These works potentially implicate m⁶A in a variety of diseases, but work remains to establish whether m⁶A is the causative mechanism that results in the higher risk of disease or if it is a secondary effect simply correlated with disease states.

Epigenome control

Lastly, one additional role for modifications of RNA such as m⁶A could be to serve as a regulator of epigenetic reprogramming and chromatin modification. Specifically, recent work on Alkbh5 (RNA demethylation) has shown that apopotosis- and spermatogenesis-related genes are differentially expressed in Alkbh5-deficient testicular cells when compared with wild-type testicular cells.¹⁹ Moreover, several key epigenetic regulatory genes are downregulated in the Alkbh5-deficient cells, including DNA-methyltransferases (Dmnt1), RNA polymerase II elongation factor (Ell3) and Bromodomain-containing protein Brd4. Given the importance of epigenetic reprogramming in development, and with the additional evidence that sperm motility and testicular function is impaired without proper RNA methylation control,¹⁹ it is possible that an “epigenome-epitranscriptome” regulatory network may exist in development that is then perturbed in diseases as well.

Conclusion

In summary, despite its known molecular dynamics and the knowledge thousands of its targets, the functions of m⁶A are still being discovered. Tremendous strides have recently been made toward clarifying key questions related to m⁶A biology, but much work remains to detail the real impact of m⁶A in other tissues and organisms. However, for the first time, the regulation, abundance, dynamics and localization of m⁶A can now be tracked. These epitranscriptomic data will allow us to test key hypotheses related to the function of m⁶A, providing the foundation for understanding the role of this critical RNA modification and reveal new aspects of RNA biology.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.