m6A-LAIC-seq reveals the census and complexity of the m6A epitranscriptome

m6A-LAIC-seq quantifies the levels of N6-methyladenosine in isoforms of the same gene in a cell-type-specific manner. N6-Methyladenosine (m6A) is a widespread, reversible chemical modification of RNA molecules, implicated in many aspects of RNA metabolism. Little quantitative information exists as to either how many transcript copies of particular genes are m6A modified ('m6A levels') or the relationship of m6A modification(s) to alternative RNA isoforms. To deconvolute the m6A epitranscriptome, we developed m6A-level and isoform-characterization sequencing (m6A-LAIC-seq). We found that cells exhibit a broad range of nonstoichiometric m6A levels with cell-type specificity. At the level of isoform characterization, we discovered widespread differences in the use of tandem alternative polyadenylation (APA) sites by methylated and nonmethylated transcript isoforms of individual genes. Strikingly, there is a strong bias for methylated transcripts to be coupled with proximal APA sites, resulting in shortened 3′ untranslated regions, while nonmethylated transcript isoforms tend to use distal APA sites. m6A-LAIC-seq yields a new perspective on transcriptome complexity and links APA usage to m6A modifications.