Search

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.

In plants, 22-nucleotide small RNAs trigger the production of secondary small interfering RNAs (siRNAs) and enhance silencing. DICER-LIKE2 (DCL2)-dependent 22-nucleotide siRNAs are rare in Arabidopsis (Arabidopsis thaliana) and are thought to function mainly during viral infection; by contrast, these siRNAs are abundant in many crops such as soybean (Glycine max) and maize (Zea mays). Here, we studied soybean 22-nucleotide siRNAs by applying CRISPR-Cas9 to simultaneously knock out the two copies of soybean DCL2, GmDCL2a and GmDCL2b, in the Tianlong1 cultivar. Small RNA sequencing revealed that most 22-nucleotide siRNAs are derived from long inverted repeats (LIRs) and disappeared in the Gmdcl2a/2b double mutant. De novo assembly of a Tianlong1 reference genome and transcriptome profiling identified an intronic LIR formed by the chalcone synthase (CHS) genes CHS1 and CHS3 This LIR is the source of primary 22-nucleotide siRNAs that target other CHS genes and trigger the production of secondary 21-nucleotide siRNAs. Disruption of this process in Gmdcl2a/2b mutants substantially increased CHS mRNA levels in the seed coat, thus changing the coat color from yellow to brown. Our results demonstrated that endogenous LIR-derived transcripts in soybean are predominantly processed by GmDCL2 into 22-nucleotide siRNAs and uncovered a role for DCL2 in regulating natural traits.

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.

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.

// Na Li 1, 2 , Yunhuan Yan 1, 2 , Angke Zhang 1, 2 , Jiming Gao 1, 2 , Chong Zhang 1, 2 , Xue Wang 1, 2 , Gaopeng Hou 1, 2 , Gaiping Zhang 1, 3 , Jinbu Jia 4, 5 , En-Min Zhou 1, 2 , Shuqi Xiao 1, 2 1 College of Veterinary Medicine, Northwest A&F University, Yangling, Shaanxi 712100, China 2 Experimental Station of Veterinary Pharmacology and Veterinary Biotechnology, Ministry of Agriculture, Yangling, Shaanxi 712100, China 3 College of Animal Science and Veterinary Medicine, Henan Agricultural University, Zhengzhou, Henan 450002, China 4 College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China 5 State Key Laboratory of Crop Stress Biology for Arid Areas, College of Plant Protection, Northwest A&F University, Yangling, Shaanxi 712100, China Correspondence to: Shuqi Xiao, email: xiaoshuqi@nwsuaf.edu.cn En-Min Zhou, email: zhouem@nwsuaf.edu.cn Keywords: microRNA, PRRSV, viral small RNAs (vsRNAs), viral replication, autoregulation Received: July 05, 2016 Accepted: October 03, 2016 Published: October 17, 2016 ABSTRACT Many viruses encode microRNAs (miRNAs) that are small non-coding single-stranded RNAs which play critical roles in virus-host interactions. Porcine reproductive and respiratory syndrome virus (PRRSV) is one of the most economically impactful viruses in the swine industry. The present study sought to determine whether PRRSV encodes miRNAs that could regulate PRRSV replication. Four viral small RNAs (vsRNAs) were mapped to the stem-loop structures in the ORF1a, ORF1b and GP2a regions of the PRRSV genome by bioinformatics prediction and experimental verification. Of these, the structures with the lowest minimum free energy (MFE) values predicted for PRRSV-vsRNA1 corresponded to typical stem-loop, hairpin structures. Inhibition of PRRSV-vsRNA1 function led to significant increases in viral replication. Transfection with PRRSV-vsRNA1 mimics significantly inhibited PRRSV replication in primary porcine alveolar macrophages (PAMs). The time-dependent increase in the abundance of PRRSV-vsRNA1 mirrored the gradual upregulation of PRRSV RNA expression. Knockdown of proteins associated with cellular miRNA biogenesis demonstrated that Drosha and Argonaute (Ago2) are involved in PRRSV-vsRNA1 biogenesis. Moreover, PRRSV-vsRNA1 bound specifically to the nonstructural protein 2 (NSP2)-coding sequence of PRRSV genome RNA. Collectively, the results reveal that PRRSV encodes a functional PRRSV-vsRNA1 which auto-regulates PRRSV replication by directly targeting and suppressing viral NSP2 gene expression. These findings not only provide new insights into the mechanism of the pathogenesis of PRRSV, but also explore a potential avenue for controlling PRRSV infection using viral small RNAs.

Small RNAs (sRNAs) are important non-coding RNA regulators, playing key roles in developmental regulation, transposon suppression, environmental response, host–pathogen interaction and other diverse biological processes. However, their roles in oomycetes are poorly understood. Here, we performed sRNA sequencing and RNA sequencing of Phytophthora parasitica at stages of vegetative growth and infection of Arabidopsis roots to examine diversity and function of sRNAs in P. parasitica, a model hemibiotrophic oomycete plant pathogen. Our results indicate that there are two distinct types of sRNA-generating loci in P. parasitica genome, giving rise to clusters of 25–26 nt and 21 nt sRNAs, respectively, with no significant strand-biases. The 25–26 nt sRNA loci lie predominantly in gene-sparse and repeat-rich regions, and overlap with over 7000 endogenous gene loci. These overlapped genes are typically P. parasitica species-specific, with no homologies to the sister species P. infestans. They include approximately 40% RXLR effector genes, 50% CRN effector genes and some elicitor genes. The transcripts of most of these genes could not be detected at both the vegetative mycelium and infection stages as revealed by RNA sequencing, indicating that the 25–26 nt sRNAs are associated with efficient silencing of these genes. The 21 nt sRNA loci typically overlap with the exon regions of highly expressed genes, suggesting that the biogenesis of the 21 nt sRNAs may be dependent on the level of gene transcription and that these sRNAs do not mediate efficient silencing of homologous genes. Analyses of the published P. infestans sRNA and mRNA sequencing data consistently show that the 25–26 nt sRNAs, but not the 21 nt sRNAs, may mediate efficient gene silencing in Phytophthora.