Understanding genome organization and gene regulation requires insight into RNA transcription, processing and modification. We adapted nanopore direct RNA sequencing to examine RNA from a wild-type accession of the model plant Arabidopsis thaliana and a mutant defective in mRNA methylation (m6A). Here we show that m6A can be mapped in full-length mRNAs transcriptome-wide and reveal the combinatorial diversity of cap-associated transcription start sites, splicing events, poly(A) site choice and poly(A) tail length. Loss of m6A from 3’ untranslated regions is associated with decreased relative transcript abundance and defective RNA 3′ end formation. A functional consequence of disrupted m6A is a lengthening of the circadian period. We conclude that nanopore direct RNA sequencing can reveal the complexity of mRNA processing and modification in full-length single molecule reads. These findings can refine Arabidopsis genome annotation. Further, applying this approach to less well-studied species could transform our understanding of what their genomes encode.

N6-methyladenosine (m6A) is the most abundant internal modification of eukaryotic mRNA. This modification has previously been shown to alter the export kinetics for mRNAs though the molecular details surrounding this phenomenon remain poorly understood. Recruitment of the TREX mRNA export complex to mRNA is driven by transcription, 5′ capping and pre-mRNA splicing. Here we identify a fourth mechanism in human cells driving the association of TREX with mRNA involving the m6A methylase complex. We show that the m6A complex recruits TREX to m6A modified mRNAs and this process is essential for their efficient export. TREX also stimulates recruitment of the m6A reader protein YTHDC1 to the mRNA and the m6A complex influences the interaction of TREX with YTHDC1. Together our studies reveal a key role for TREX in the export of m6A modified mRNAs.

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

Genes involved in disease resistance are some of the fastest evolving and most diverse components of genomes. Large numbers of nucleotide-binding, leucine-rich repeat (NLR) genes are found in plant genomes and are required for disease resistance. However, NLRs can trigger autoimmunity, disrupt beneficial microbiota or reduce fitness. It is therefore crucial to understand how NLRs are controlled. Here, we show that the RNA-binding protein FPA mediates widespread premature cleavage and polyadenylation of NLR transcripts, thereby controlling their functional expression and impacting immunity. Using long-read Nanopore direct RNA sequencing, we resolved the complexity of NLR transcript processing and gene annotation. Our results uncover a co-transcriptional layer of NLR control with implications for understanding the regulatory and evolutionary dynamics of NLRs in the immune responses of plants.

Alternative splicing of messenger RNAs is associated with the evolution of developmentally complex eukaryotes. Splicing is mediated by the spliceosome, and docking of the pre-mRNA 5' splice site into the spliceosome active site depends upon pairing with the conserved ACAGA sequence of U6 snRNA. In some species, including humans, the central adenosine of the ACAGA box is modified by N6 methylation, but the role of this m6A modification is poorly understood. Here we show that m6A modified U6 snRNA determines the accuracy and efficiency of splicing. We reveal that the conserved methyltransferase, FIO1, is required for Arabidopsis U6 snRNA m6A modification. Arabidopsis fio1 mutants show disrupted patterns of splicing that can be explained by the sequence composition of 5' splice sites and cooperative roles for U5 and U6 snRNA in splice site selection. U6 snRNA m6A influences 3' splice site usage. We generalise these findings to reveal two major classes of 5' splice site in diverse eukaryotes, which display anti-correlated interaction potential with U5 snRNA loop 1 and the U6 snRNA ACAGA box. We conclude that U6 snRNA m6A modification contributes to the selection of degenerate 5' splice sites crucial to alternative splicing.

Yanocomp is a tool for predicting the positions and stoichiometries of RNA modifications in Nanopore direct RNA sequencing data. It uses general mixture models to identify differentially modified sites between two conditions, with good support for replicates. Yanocomp models across adjacent kmers and uses a uniform component to account for outliers, improving the accuracy of single molecule predictions. Consequently, Yanocomp can be used to measure modification stoichiometry, and correlate modifications with other RNA processing events. Yanocomp is available under an MIT license at www.github.com/bartongroup/yanocomp.

The nutrient-rich tubers of the greater yam Dioscorea alata L. provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an "orphan crop." Here we address this resource gap by presenting a highly-contiguous chromosome-scale genome assembly of greater yam combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient lineage-specific genome duplication, followed by extensive genome-wide reorganization. Using our new genomic tools we find quantitative trait loci for susceptibility to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

21Understanding genome organization and gene regulation requires insight into RNA transcription, 22 processing and modification. We adapted nanopore direct RNA sequencing to examine RNA from a 23 wild-type accession of the model plant Arabidopsis thaliana and a mutant defective in mRNA 24 methylation (m 6 A). Here we show that m 6 A can be mapped in full-length mRNAs transcriptome-wide 25 and reveal the combinatorial diversity of cap-associated transcription start sites, splicing events, 26 poly(A) site choice and poly(A) tail length. Loss of m 6 A from 3' untranslated regions is associated 27 with decreased relative transcript abundance and defective RNA 3′ end formation. A functional 28 consequence of disrupted m 6 A is a lengthening of the circadian period. We conclude that nanopore 29 direct RNA sequencing can reveal the complexity of mRNA processing and modification in full-length 30 single molecule reads. These findings can refine Arabidopsis genome annotation. Further, applying 31 this approach to less well-studied species could transform our understanding of what their genomes 32 encode. 33 34 misidentification of 3′ ends through internal priming 3 , spurious antisense and splicing events 46 produced by RT template switching 4,5 , and the inability to detect all base modifications in the 47 copying process 6 . The fragmentation of RNA prior to short-read sequencing makes it difficult to 48 interpret the combination of authentic RNA processing events and remains an unsolved problem 7 . 49We investigated whether long-read direct RNA sequencing (DRS) with nanopores 8 could 50 reveal the complexity of Arabidopsis mRNA processing and modifications. In nanopore DRS, the 51 protein pore (nanopore) sits in a membrane through which an electrical current is passed, and intact 52 RNA is fed through the nanopore by a motor protein 8 . Each RNA sequence within the nanopore 53 (5 bases) can be identified by the magnitude of signal it produces. Arabidopsis is a pathfinder model 54 in plant biology, and its genome annotation strongly influences the annotation and our 55 understanding of what other plant genomes encode. We applied nanopore DRS and Illumina RNAseq 56 to wild-type Arabidopsis (Col-0) and mutants defective in m 6 A 9 and exosome-mediated RNA decay 10 . 57We reveal m 6 A and combinations of RNA processing events (alternative patterns of 5′ capped 58 transcription start sites, splicing, 3′ polyadenylation and poly(A) tail length) in full-length Arabidopsis 59 mRNAs transcriptome-wide. 60 61 Results 62Nanopore DRS detects long, complex mRNAs and short, structured non-coding RNAs 63We purified poly (A)+ RNA from four biological replicates of 14-day-old Arabidopsis Col-0 seedlings. 64We incorporated synthetic External RNA Controls Consortium (ERCC) RNA Spike-In mixes into all 65 replicates 11,12 and carried out nanopore DRS. Illumina RNAseq was performed in parallel on similar 66 material. Using Guppy base-calling (Oxford Nanopore Technologies) and minimap2 alignment 67 software 13 , we identified around 1 mi...