Search

Poly(A) tail is a hallmark of eukaryotic messenger RNA and its length plays an essential role in regulating mRNA metabolism. However, a comprehensive resource for plant poly(A) tail length has yet to be established. Here, we applied a poly(A)-enrichment-free, nanopore-based method to profile full-length RNA with poly(A) tail information in plants. Our atlas contains over 120 million polyadenylated mRNA molecules from seven different tissues of Arabidopsis, as well as the shoot tissue of maize, soybean and rice. In most tissues, the size of plant poly(A) tails shows peaks at approximately 20 and 45 nucleotides, while the poly(A) tails in pollen exhibit a distinct pattern with strong peaks centred at 55 and 80 nucleotides. Moreover, poly(A) tail length is regulated in a gene-specific manner-mRNAs with short half-lives in general have long poly(A) tails, while mRNAs with long half-lives are featured with relatively short poly(A) tails that peak at ~45 nucleotides. Across species, poly(A) tails in the nucleus are almost twice as long as in the cytoplasm. Our comprehensive dataset lays the groundwork for future functional and evolutionary studies on poly(A) tail length regulation in plants.

Elongation, splicing and polyadenylation are fundamental steps of transcription, and studying their coordination requires simultaneous monitoring of these dynamic processes on one transcript. We recently developed a full-length nascent RNA sequencing method in the model plant Arabidopsis that simultaneously detects RNA polymerase II position, splicing status, polyadenylation site and poly(A) tail length at genome-wide scale. This method allows calculation of the kinetics of cotranscriptional splicing and detects polyadenylated transcripts with unspliced introns retained at specific positions posttranscriptionally. Here we describe a detailed protocol for this method called FLEP-seq (full-length elongating and polyadenylated RNA sequencing) that is applicable to plants. Library production requires as little as one nanogram of nascent RNA (after rRNA/tRNA removal), and either Nanopore or PacBio platforms can be used for sequencing. We also provide a complete bioinformatic pipeline from raw data processing to downstream analysis. The minimum time required for FLEP-seq, including RNA extraction and library preparation, is 36 h. The subsequent long-read sequencing and initial data analysis ranges between 31 and 40 h, depending on the sequencing platform.

Phytophthora species are destructive plant pathogens that cause significant crop losses worldwide. To understand plant susceptibility to oomycete pathogens and to explore novel disease resistance strategies, we employed the Arabidopsis thaliana–Phytophthora parasitica model pathosystem and screened for A. thaliana T‐DNA insertion mutant lines resistant to P. parasitica. This led to the identification of the resistant mutant 267‐31, which carries two T‐DNA insertion sites in the promoter region of the ethylene‐responsive factor 19 gene (ERF019). Quantitative reverse transcription PCR (RT‐qPCR) assays showed that the expression of ERF019 was induced during P. parasitica infection in the wild type, which was suppressed in the 267‐31 mutant. Additional erf019 mutants were generated using CRISPR/Cas9 technology and were confirmed to have increased resistance to P. parasitica. In contrast, ERF019 overexpression lines were more susceptible. Transient overexpression assays in Nicotiana benthamiana showed that the nuclear localization of ERF019 is crucial for its susceptible function. RT‐qPCR analyses showed that the expression of marker genes for multiple defence pathways was significantly up‐regulated in the mutant compared with the wild type during infection. Flg22‐induced hydrogen peroxide accumulation and reactive oxygen species burst were impaired in ERF019 overexpression lines, and flg22‐induced MAPK activation was enhanced in erf019 mutants. Moreover, transient overexpression of ERF019 strongly suppressed INF‐triggered cell death in N. benthamiana. These results reveal the importance of ERF019 in mediating plant susceptibility to P. parasitica through suppression of pathogen‐associated molecular pattern‐triggered immunity., ERF019 is induced during Phytophthora parasitica infection and acts as a negative regulator of plant resistance by suppressing PAMP‐triggered immunity.

In eukaryotes, genes are transcribed by RNA polymerase-II (Pol-II) and introns are removed by the spliceosome largely cotranscriptionally1–3; analysis using long-read sequencing revealed that splicing occurs immediately after Pol-II passes introns in yeast4,5. Here, we developed a Nanopore-based method to profile chromatin-bound RNA that enables the simultaneous detection of splicing status, Pol-II position and polyadenylation at the genome-wide scale in Arabidopsis. We found that more than half of the introns remain unspliced after Pol-II transcribes 1 kb past the 3′ splice site, which is much slower than the rate of splicing reported in yeast4,5. Many of the full-length chromatin-bound RNA molecules are polyadenylated, yet still contain unspliced introns at specific positions. These introns are nearly absent in the cytoplasm and are resistant to nonsense-mediated decay, suggesting that they are post-transcriptionally spliced before the transcripts are released into the cytoplasm; we therefore termed these introns post-transcriptionally spliced introns (pts introns). Analysis of around 6,500 public RNA-sequencing libraries found that the splicing of pts introns requires the function of splicing-related proteins such as PRMT5 and SKIP, and is also influenced by various environmental signals. The majority of the intron retention events in Arabidopsis are at pts introns, suggesting that chromatin-tethered post-transcriptional splicing is a major contributor to the widespread intron retention that is observed in plants, and could be a mechanism to produce fully spliced functional mRNAs for rapid response. Chromatin-bound RNA profiling in Arabidopsis using a Nanopore-based method reveals slower intron splicing than reported in yeast. Many introns are post-transcriptionally spliced before transcripts enter the cytoplasm, which contributes to intron retention.

In the past decade, genome-wide characterization of the three-dimensional chromatin structure in plants using high-throughput methods has greatly advanced our knowledge in plant genome architecture (Liu and Weigel, 2015; Ouyang et al., 2020). However, due to the limitation of Illumina short-read sequencing, the genome-wide contact map obtained by Hi-C/ChIA-PET is pairwise, and the multi-way interaction can only be inferred from the two-way data. To directly capture multi-way interaction in Arabidopsis, we applied a long-read-based method called Pore-C that directly sequences the DNA fragments joined by proximity-based ligation (Ulahannan et al., 2019).

Background The dynamic process of transcription termination produces transient RNA intermediates that are difficult to distinguish from each other via short-read sequencing methods. Results Here, we use single-molecule nascent RNA sequencing to characterize the various forms of transient RNAs during termination at genome-wide scale in wildtype Arabidopsis and in atxrn3, fpa, and met1 mutants. Our data reveal a wide range of termination windows among genes, ranging from ~ 50 nt to over 1000 nt. We also observe efficient termination before downstream tRNA genes, suggesting that chromatin structure around the promoter region of tRNA genes may block pol II elongation. 5′ Cleaved readthrough transcription in atxrn3 with delayed termination can run into downstream genes to produce normally spliced and polyadenylated mRNAs in the absence of their own transcription initiation. Consistent with previous reports, we also observe long chimeric transcripts with cryptic splicing in fpa mutant; but loss of CG DNA methylation has no obvious impact on termination in the met1 mutant. Conclusions Our method is applicable to establish a comprehensive termination landscape in a broad range of species.

The broad application of single-cell RNA profiling in plants has been hindered by the prerequisite of protoplasting that requires digesting the cell walls from different types of plant tissues. Here, we present a protoplasting-free approach, flsnRNA-seq, for large-scale full-length RNA profiling at a single-nucleus level in plants using isolated nuclei. Combined with 10x Genomics and Nanopore long-read sequencing, we validate the robustness of this approach in Arabidopsis root cells and the developing endosperm. Sequencing results demonstrate that it allows for uncovering alternative splicing and polyadenylation-related RNA isoform information at the single-cell level, which facilitates characterizing cell identities. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-021-02288-0.

The broad application of large-scale single-cell RNA profiling in plants has been restricted by the prerequisite of protoplasting. We recently found that the Arabidopsis nucleus contains abundant polyadenylated mRNAs, many of which are incompletely spliced. To capture the isoform information, we combined 10x Genomics and Nanopore long-read sequencing to develop a protoplasting-free full-length single-nucleus RNA profiling method in plants. Our results demonstrated using Arabidopsis root that nuclear mRNAs faithfully retain cell identity information, and single-molecule full-length RNA sequencing could further improve cell type identification by revealing splicing status and alternative polyadenylation at single-cell level.

In eukaryotes, genes are transcribed by RNA polymerase-II (Pol-II) and introns are removed by the spliceosome largely cotranscriptionally

Application of Next Generating Sequencing (NGS) technology in transcriptome profiling has greatly improved our understanding of transcriptional regulation at genome-wide scale in the last decade, and tens of thousands of RNA-sequencing (RNA-seq) libraries have been produced by the research community. However, accessing such huge amount of RNA-seq data poses a big challenge for groups that lack dedicated bioinformatic personnel or expensive computational resources. Here, we introduce the Arabidopsis RNA-seq database (ARS), a free, web-accessible, and user-friendly to quickly explore expression level of any gene in 20,000+ publicly available Arabidopsis RNA-seq libraries.