The variables on RNA molecules: concert or cacophony? Answers in long-read sequencing

RNA methylation to form N6-methyladenosine (m6A) in mRNA accounts for the most abundant mRNA internal modification and has emerged as a widespread regulatory mechanism that controls gene expression in diverse physiological processes. Transcriptome-wide m6A mapping has revealed the distribution and pattern of m6A in cellular RNAs, referred to as the epitranscriptome. These maps have revealed the specific mRNAs that are regulated by m6A, providing mechanistic links connecting m6A to cellular differentiation, cancer progression and other processes. The effects of m6A on mRNA are mediated by an expanding list of m6A readers and m6A writer-complex components, as well as potential erasers that currently have unclear relevance to m6A prevalence in the transcriptome. Here we review new and emerging methods to characterize and quantify the epitranscriptome, and we discuss new concepts - in some cases, controversies - regarding our understanding of the mechanisms and functions of m6A readers, writers and erasers.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion base pair (bp) sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million bp of sequence containing 1,956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies. Twenty years after the initial drafts, a truly complete sequence of a human genome reveals what has been missing.

Through alternative processing of pre-mRNAs, individual mammalian genes often produce multiple mRNA and protein isoforms that may have related, distinct or even opposing functions. Here we report an in-depth analysis of 15 diverse human tissue and cell line transcriptomes based on deep sequencing of cDNA fragments, yielding a digital inventory of gene and mRNA isoform expression. Analysis of mappings of sequence reads to exon-exon junctions indicated that 92-94% of human genes undergo alternative splicing (AS), ∼86% with a minor isoform frequency of 15% or more. Differences in isoform-specific read densities indicated that a majority of AS and of alternative cleavage and polyadenylation (APA) events vary between tissues, while variation between individuals was ∼2- to 3-fold less common. Extreme or ‘switch-like’ regulation of splicing between tissues was associated with increased sequence conservation in regulatory regions and with generation of full-length open reading frames. Patterns of AS and APA were strongly correlated across tissues, suggesting coordinated regulation of these processes, and sequence conservation of a subset of known regulatory motifs in both alternative introns and 3′ UTRs suggested common involvement of specific factors in tissue-level regulation of both splicing and polyadenylation.