Search

Intron retention is the most common alternative splicing event in plants and plays a crucial role in the responses of plants to environmental signals. Despite a large number of RNA-seq libraries from different treatments and genetic mutants stored in public domains, a resource for querying the intron-splicing ratio of individual intron is still required. Here, we established the first-ever large-scale splicing efficiency database in any organism. Our database includes over 57,000 plant public RNA-seq libraries, comprising 25,283 from Arabidopsis, 17,789 from maize, 10,710 from rice, and 3,974 from soybean, and covers a total of 1.6 million introns in these four species. In addition, we manually curated and annotated all the mutant- and treatment-related libraries as well as their matched controls included in our library collection, and added graphics to display intron-splicing efficiency across various tissues, developmental stages, and stress-related conditions. The result is a large collection of 3,313 treatment conditions and 3,594 genetic mutants for discovering differentially regulated splicing efficiency. Our online database can be accessed at https://plantintron.com/ .

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.

Legumes form symbiosis with rhizobium leading to the development of nitrogen-fixing nodules. By integrating single-nucleus and spatial transcriptomics, we established a cell atlas of soybean nodules and roots. In central infected zone of nodule, we found that uninfected cells specialize into functionally distinct sub-groups during nodule development and revealed a transitional subtype of infected cells with enriched nodulation-related genes. Overall, our results provide a single-cell perspective for understanding rhizobium-legume symbiosis.

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.

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.

Introductory paragraphPoly(A) tail is a hallmark of eukaryotic mRNA and its length plays an essential role in regulating mRNA metabolism1–3. 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 method4,5 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 nt, presumably the sizes protected by one and two poly(A) binding proteins (PABP), respectively2,6, while the poly(A) tails in pollen exhibit a distinct pattern with strong peaks centered at 55 and 80 nt. 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 nt, suggesting that protection of poly(A) tail in this size by PABP is essential for mRNA stability. Across species, poly(A) tails in the nucleus are almost twice as long as in the cytoplasm, implying a conserved rapid shortening process of poly(A) tail occurs before the mRNA is stabilized in cytoplasm. Our comprehensive dataset lays the groundwork for future functional and evolutionary studies on poly(A) tail length regulation in plants.

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.

RXLR effectors encoded by Phytophthora species play a central role in pathogen-plant interactions. An understanding of the biological functions of RXLR effectors is conducive to the illumination of the pathogenic mechanisms and the development of disease control strategies. However, the virulence function of Phytophthora parasitica RXLR effectors is poorly understood. Here, we describe the identification of a P. parasitica RXLR effector gene, PPTG00121 (PpE4), which is highly transcribed during the early stages of infection. Live cell imaging of P. parasitica transformants expressing a full-length PpE4 (E4FL)-mCherry protein indicated that PpE4 is secreted and accumulates around haustoria during plant infection. Silencing of PpE4 in P. parasitica resulted in significantly reduced virulence on Nicotiana benthamiana. Transient expression of PpE4 in N. benthamiana in turn restored the pathogenicity of the PpE4-silenced lines. Furthermore, the expression of PpE4 in both N. benthamiana and Arabidopsis thaliana consistently enhanced plant susceptibility to P. parasitica. These results indicate that PpE4 contributes to pathogen infection. Finally, heterologous expression experiments showed that PpE4 triggers non-specific cell death in a variety of plants, including tobacco, tomato, potato and A. thaliana. Virus-induced gene silencing assays revealed that PpE4-induced cell death is dependent on HSP90, NPK and SGT1, suggesting that PpE4 is recognized by the plant immune system. In conclusion, PpE4 is an important virulence RXLR effector of P. parasitica and recognized by a wide range of host plants.

Additional file 3. Review history.

Additional file 1: Fig. S1. Schematic diagram of snuupy bioinformatic pipeline. Fig. S2. The sorted nuclei were observed under a microscopy with DAPI staining. Fig. S3. Identification of clusters by a marker-gene-based method. Fig. S4. UMAP visualization of the representative cell-type marker genes for each of the 14 cell clusters. Fig. S5. UMAP visualization showing the abundances of representative marker genes in two subcell types of endodermis. Fig. S6. Scheme for deriving the splicing and APA matrices from Nanopore data. Fig. S7. The gene expression matrices of endosperm generated from the two different libraries were similar to each other. (PPTX 4884 kb)

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.

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.

// 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.

The interaction between Arabidopsis thaliana and the oomycete pathogen Phytophthora parasitica emerges as a model for exploring the molecular basis and evolution of recognition and host defense. Phenotypic variation and genetic analysis is essential to dissect the underlying mechanisms in plant-oomycete interaction. In this study, the reaction phenotypes of 28 A. thaliana accessions to P. parasitica strain Pp016 were examined using detached leaf infection assay. The results showed the presence of four distinct groups based on host response and disease development. Of all the accessions examined, Zurich (Zu-1) is highly resistant to P. parasitica. Microscopic characterization showed that rapid and severe hypersensitive response (HR) at the primary infection epidermal cells is associated with disease resistance. Furthermore, Zu-1 is resistant to a set of 20 diverse P. parasitica strains, which were collected from different host plants and exhibited differential specificities on a set of tobacco cultivars. However, Zu-1 is susceptible to P. parasitica when the root is inoculated, suggesting differential expression of associated resistance genes in the root and foliar tissues. Genetic analysis by crossing Zu-1 and the susceptible accession Landsberg (Ler) showed that the resistance in Zu-1 to P. parasitica is semi-dominant, as shown by infection assays of F1 progenies, and is likely conferred by a single locus, defined as RPPA1Zu-1 (for Resistance to Phytophthora parasitica 1), as shown by analysis of F2 segregating populations. By employing specific-locus amplified fragment sequencing (SLAF-seq) strategy to identify molecular markers potentially linked to the locus, the strongest associated region was determined to be located between 7.1-11.2 Mb in chromosome IV. The future cloning of RPPA1Zu-1 locus will facilitate improved understanding of plant broad-spectrum disease resistance to oomycete pathogens.