Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2015 Oct 1;29(19):2037-53.
doi: 10.1101/gad.269415.115. Epub 2015 Sep 24.

A majority of m6A residues are in the last exons, allowing the potential for 3' UTR regulation

Affiliations

A majority of m6A residues are in the last exons, allowing the potential for 3' UTR regulation

Shengdong Ke et al. Genes Dev. .

Abstract

We adapted UV CLIP (cross-linking immunoprecipitation) to accurately locate tens of thousands of m(6)A residues in mammalian mRNA with single-nucleotide resolution. More than 70% of these residues are present in the 3'-most (last) exons, with a very sharp rise (sixfold) within 150-400 nucleotides of the start of the last exon. Two-thirds of last exon m(6)A and >40% of all m(6)A in mRNA are present in 3' untranslated regions (UTRs); contrary to earlier suggestions, there is no preference for location of m(6)A sites around stop codons. Moreover, m(6)A is significantly higher in noncoding last exons than in next-to-last exons harboring stop codons. We found that m(6)A density peaks early in the 3' UTR and that, among transcripts with alternative polyA (APA) usage in both the brain and the liver, brain transcripts preferentially use distal polyA sites, as reported, and also show higher proximal m(6)A density in the last exons. Furthermore, when we reduced m6A methylation by knocking down components of the methylase complex and then examined 661 transcripts with proximal m6A peaks in last exons, we identified a set of 111 transcripts with altered (approximately two-thirds increased proximal) APA use. Taken together, these observations suggest a role of m(6)A modification in regulating proximal alternative polyA choice.

Keywords: alternative polyadenylation; last exon; m6A-CLIP/IP; microRNA.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
Single-nucleotide resolution achieved by m6A-CLIP/immunoprecipitation (IP). (A) m6A-CLIP/IP illustration. (B) The m6A PS (defined in the text) maps m6A sites at its 0 position. “Enrichment of RRACU or RAC” is the fold enrichment of RRACU (red) or RAC (blue) motif density at that position compared with the background motif density (horizontal gray dot line). The background motif density was calculated as the motif density of RRACU or RAC at regions 400 nt upstream of or downstream from PSs where the density value was flattened. The vertical black dotted line represents the PS 0 position. The motif logo with its P-value represents the de novo motif. (C) Substitution of cross-linking-induced mutation sites (CIMSs) maps m6A sites at its −1 position. (D) m6A-induced truncation sites (MITSs) map m6A sites at its +1 position. (E) Summary of m6A sites identified by different approaches. “All RACs” is the total number of RACs within m6A peak regions. (F) Precisely mapped m6A sites are conserved in vertebrates regardless of the mapping approaches. Box plot bars are the first quartile to third quartile; the black line in the middle is the median value. (**) P < 10−6; (*) P < 0.01, Wilcoxon rank sum test. (G) m6A-CLIP/IP mapped a precise m6A site in TUG1 RNA by two independent means: at −2 of a deletion CIMS and also at +1 of a MITS. “m6A-IP enrich” is m6A-IP reads normalized to its input. “m6A-CLIP/IP” is by m6A-CLIP/IP (the A marked in red). SCARLET (site-specific cleavage and radioactive labeling followed by ligation-assisted extraction and thin-layer chromatography) is by Liu et al. (2013) at the underlined motif. (H) m6A-CLIP/IP mapped a precise m6A site in BSG mRNA that is located at the +1 position of MITSs. (I) m6A-CLIP/IP mapped precise m6A sites in TPT1 mRNA. From left to right, these sites were located at −1 of substitution CIMSs, +1 of MITSs, and 0 of PSs and were co-mapped by −1 of substitution CIMSs and +1 of MITSs. See also Supplemental Figures 1 and 2.
Figure 2.
Figure 2.
m6A is enriched when entering last exons but not at stop codons. (A) Entering the last exon, m6A density increased sharply (bottom panel) in contrast to its lagging increase when approaching the stop codon (top panel). Data are from mouse brains. “m6A peak density” was calculated as the number of m6A peak regions in a 10-nt interval divided by the total number of mRNAs that contained this position. (B) m6A is enriched in the last exons but not around the stop codons. The top panel is the m6A peak distribution in the region around the stop codons; the two bottom panels are the distribution of the m6A peak regions and stop codons around the last exon start. mRNAs were grouped according to their stop codon locations to the last exon start. (C) m6A is enriched in the last exons but not around the stop codon when the stop codon is not in the last exon. (Top and bottom panels) mRNA with noncoding last exons anchored at the last exon start or the stop codon. (D) Most exonic m6As locate in last exon (mouse brains). The three pie graphs show the relative proportions of m6A peaks and the RAC and RRACU motifs in the last exon and other exons, with 100% representing all m6A peaks and RAC and RRACU sequences in mRNA. See also Supplemental Figure 3 for human and previously published data.
Figure 3.
Figure 3.
Short last exons have a lower m6A density (mouse brains). The m6A density in the last exons were compared between longer and shorter last exons (long: >850 nt, 5916 genes; short: ≤850 nt, 3624 genes). All of the genes considered are multiexon-coding genes with RPKM (reads per kilobase per million mapped reads) ≥ 1 in the mouse brain. The results showing more proximal m6A in transcripts using distal polyA sites are normalized to transcript abundance: “m6A peak density” was calculated as the number of m6A peak regions in a 10-nt interval divided by the total number of mRNAs that contained this position (the same definition is used throughout the text). The error bar indicates the standard error of the mean at each position. See also Supplemental Figure 4 for human data and RAC motif controls.
Figure 4.
Figure 4.
Higher m6A levels in the last exons correlates with more distal polyA site usage. (A) m6A is more abundant near proximal sites that are used less frequently (black line; left panel); when distal sites are used more, the m6A level is higher until the polyA site is reached (black line; right panel). The error bar indicates the standard error of the mean at each position. The shaded area highlights the major area of difference between two groups. (**) P < 10−35, Fisher's exact test. “Used more“ means ≥60% usage of all polyA sites in the last exon, and “used less” means ≤40% usage. Other cutoffs, including 50% or 70%, produce the same findings. The expression of mRNAs in regions between proximal and distal polyA sites is adequate for m6A detection (RPKM ≥ 1). Mouse brain data are presented here; see also Supplemental Figure 5A for human data. (B) Distal APA sites in the last exons are used more often in brains for coexpressed mRNAs. APA site pairs with statistically significant differences in usage are highlighted in orange (higher usage of distal sites in the brain) or blue (higher usage of distal sites in the liver). FDR = 20%, Fisher's exact test. The same conclusion is true for FDR = 5% (Supplemental Fig. 5B). (C) m6A peak regions in the last exons have a higher enrichment value in the brain. “m6A peak enrich. value” is the m6A-IP reads normalized to the input reads for each m6A peak region. m6A peak regions with statistically significant differences in “m6A peak enrich. value” are highlighted in orange (higher in the brain) or blue (higher in the liver). FDR = 5%, Fisher's exact test. (D) Positional plot of m6A peaks around APA sites in the last exons. “m6A peak density (brain–liver)” for each position was calculated as the number of m6A peak regions that are significantly higher in the brain (orange points in C) in a 10-nt interval divided by the total number of mRNAs that contained this position. (Orange lines) More distal usage; (gray lines) no change in the brain. The expression of mRNAs in regions between the proximal and distal polyA sites is adequate for m6A detection in both tissues (RPKM ≥ 1). The error bar indicates the standard error of the mean. The shaded area highlights the major area of difference between two groups. (**) P < 10−20, Fisher's exact test.
Figure 5.
Figure 5.
The effects of global limitation of m6A on alternative polyadenylation. (A) Triple knockdown of METTL3, METTL14, and WTAP reduced m6A of polyA+ RNAs to 26% in human A549 cells. The three different controls are biological replicates of control siRNA. The four different knockdown samples are as follows: KD1a and KD1b are biological replicates of one siRNA set, KD2a and KD2b are biological replicates of the other siRNA set, and the two siRNA sets generate essentially the same result. (**) P < 10−7, Student's t-test. (B) m6A may contribute to the inhibition of proximal polyA sites in the last exons. (C) The global limitation of m6A changed the choice of APA sites in the last exons. APA site pairs with statistically significant differences in usage are highlighted in orange (higher usage of proximal sites in knockdown) or blue (higher usage of proximal sites in the control). FDR = 20%, Fisher's exact test (the same conclusion is true for FDR = 5%, as shown in Supplemental Fig. S6). The mRNAs using two APA sites in the last exon and with m6A present around the proximal polyA sites are summarized. See also Supplemental Figures S6 and S7.
Figure 6.
Figure 6.
Overlap of m6A and Ago-binding sites in 3′ UTRs of mouse brains. (% m6A sites in overlap) The percentage of m6A sites (±200-nt flank) overlapped with the RNA-binding protein-binding site in mouse brains (mapped by published HITS-CLIP data) (see the Materials and Methods for reference details). The control was the randomly picked RAC motif (±200-nt flank) in the same 3′ UTR as m6A. Ten independent controls were performed to estimate the statistics. The error bar indicates the standard error of the mean. (***) P < 10−50; (**) P < 10−10, one sample t-test.

References

    1. Almada AE, Wu X, Kriz AJ, Burge CB, Sharp PA. 2013. Promoter directionality is controlled by U1 snRNP and polyadenylation signals. Nature 499: 360–363. - PMC - PubMed
    1. Batista PJ, Molinie B, Wang J, Qu K, Zhang J, Li L, Bouley DM, Lujan E, Haddad B, Daneshvar K, et al. 2014. m6A RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell 15: 707–719. - PMC - PubMed
    1. Bernhardt D, Darnell JEJ. 1969. tRNA synthesis in HeLa cells: a precursor to tRNA and the effects of methionine starvation on tRNA synthesis. J Mol Biol 42: 43–56. - PubMed
    1. Bokar JA, Shambaugh ME, Polayes D, Matera AG, Rottman FM. 1997. Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferase. RNA 3: 1233–1247. - PMC - PubMed
    1. Bringmann P, Luhrmann R. 1987. Antibodies specific for N6-methyladenosine react with intact snRNPs U2 and U4/U6. FEBS Lett 213: 309–315. - PubMed

Publication types

Associated data