A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants.

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

Light controls mesophyll-specific post-transcriptional splicing of photoregulatory genes by AtPRMT5.

Light plays a central role in plant growth and development, providing an energy source and governing various aspects of plant morphology. Previous study showed that many polyadenylated full-length RNA molecules within the nucleus contain unspliced introns (post-transcriptionally spliced introns, PTS introns), which may play a role in rapidly responding to changes in environmental signals. However, the mechanism underlying post-transcriptional regulation during initial light exposure of young, etiolated seedlings remains elusive. In this study, we used FLEP-seq2, a Nanopore-based sequencing technique, to analyze nuclear RNAs in Arabidopsis (Arabidopsis thaliana) seedlings under different light conditions and found numerous light-responsive PTS introns. We also used single-nucleus RNA sequencing (snRNA-seq) to profile transcripts in single nucleus and investigate the distribution of light-responsive PTS introns across distinct cell types. We established that light-induced PTS introns are predominant in mesophyll cells during seedling de-etiolation following exposure of etiolated seedlings to light. We further demonstrated the involvement of the splicing-related factor A. thaliana PROTEIN ARGININE METHYLTRANSFERASE 5 (AtPRMT5), working in concert with the E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), a critical repressor of light signaling pathways. We showed that these two proteins orchestrate light-induced PTS events in mesophyll cells and facilitate chloroplast development, photosynthesis, and morphogenesis in response to ever-changing light conditions. These findings provide crucial insights into the intricate mechanisms underlying plant acclimation to light at the cell-type level.

Genome evolution and initial breeding of the Triticeae grass Leymus chinensis dominating the Eurasian Steppe.

Leymus chinensis, a dominant perennial grass in the Eurasian Steppe, is well known for its remarkable adaptability and forage quality. Hardly any breeding has been done on the grass, limiting its potential in ecological restoration and forage productivity. To enable genetic improvement of the untapped, important species, we obtained a 7.85-Gb high-quality genome of L. chinensis with a particularly long contig N50 (318.49 Mb). Its allotetraploid genome is estimated to originate 5.29 million years ago (MYA) from a cross between the Ns-subgenome relating to Psathyrostachys and the unknown Xm-subgenome. Multiple bursts of transposons during 0.433-1.842 MYA after genome allopolyploidization, which involved predominantly the Tekay and Angela of LTR retrotransposons, contributed to its genome expansion and complexity. With the genome resource available, we successfully developed a genetic transformation system as well as the gene-editing pipeline in L. chinensis. We knocked out the monocot-specific miR528 using CRISPR/Cas9, resulting in the improvement of yield-related traits with increases in the tiller number and growth rate. Our research provides valuable genomic resources for Triticeae evolutionary studies and presents a conceptual framework illustrating the utilization of genomic information and genome editing to accelerate the improvement of wild L. chinensis with features such as polyploidization and self-incompatibility.

Cell type-specific cytonuclear coevolution in three allopolyploid plant species.

Cytonuclear disruption may accompany allopolyploid evolution as a consequence of the merger of different nuclear genomes in a cellular environment having only one set of progenitor organellar genomes. One path to reconcile potential cytonuclear mismatch is biased expression for maternal gene duplicates (homoeologs) encoding proteins that target to plastids and/or mitochondria. Assessment of this transcriptional form of cytonuclear coevolution at the level of individual cells or cell types remains unexplored. Using single-cell (sc-) and single-nucleus (sn-) RNAseq data from eight tissues in three allopolyploid species, we characterized cell type-specific variations of cytonuclear coevolutionary homoeologous expression and demonstrated the temporal dynamics of expression patterns across development stages during cotton fiber development. Our results provide unique insights into transcriptional cytonuclear coevolution in plant allopolyploids at the single-cell level.

Single-nucleus transcriptomes reveal spatiotemporal symbiotic perception and early response in Medicago.

Establishing legume-rhizobial symbiosis requires precise coordination of complex responses in a time- and cell type-specific manner. Encountering Rhizobium, rapid changes of gene expression levels in host plants occur in the first few hours, which prepare the plants to turn off defence and form a symbiotic relationship with the microbes. Here, we applied single-nucleus RNA sequencing to characterize the roots of Medicago truncatula at 30 min, 6 h and 24 h after nod factor treatment. We found drastic global gene expression reprogramming at 30 min in the epidermis and cortex and most of these changes were restored at 6 h. Moreover, plant defence response genes are activated at 30 min and subsequently suppressed at 6 h in non-meristem cells. Only in the cortical cells but not in other cell types, we found the flavonoid synthase genes required to recruit rhizobia are highly expressed 30 min after inoculation with nod factors. A gene module enriched for symbiotic nitrogen fixation genes showed that MtFER (MtFERONIA) and LYK3 (LysM domain receptor-like kinase 3) share similar responses to symbiotic signals. We further found that MtFER can be phosphorylated by LYK3 and it participates in rhizobial symbiosis. Our results expand our understanding of dynamic spatiotemporal symbiotic responses at the single-cell level.

Single-molecule targeted accessibility and methylation sequencing of centromeres, telomeres and rDNAs in Arabidopsis.

The short read-length of next-generation sequencing makes it challenging to characterize highly repetitive regions (HRRs) such as centromeres, telomeres and ribosomal DNAs. Based on recent strategies that combined long-read sequencing and exogenous enzymatic labelling of open chromatin, we developed single-molecule targeted accessibility and methylation sequencing (STAM-seq) in plants by further integrating nanopore adaptive sampling to investigate the HRRs in wild-type Arabidopsis and DNA methylation mutants that are defective in CG- or non-CG methylation. We found that CEN180 repeats show higher chromatin accessibility and lower DNA methylation on their forward strand, individual rDNA units show a negative correlation between their DNA methylation and accessibility, and both accessibility and CHH methylation levels are lower at telomere compared to adjacent subtelomeric region. Moreover, DNA methylation-deficient mutants showed increased chromatin accessibility at HRRs, consistent with the role of DNA methylation in maintaining heterochromatic status in plants. STAM-seq can be applied to study accessibility and methylation of repetitive sequences across diverse plant species.

Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule maturation.

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 zones of nodules, we found that uninfected cells specialize into functionally distinct subgroups 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.

Plant Intron-Splicing Efficiency Database (PISE): exploring splicing of ∼1,650,000 introns in Arabidopsis, maize, rice, and soybean from ∼57,000 public RNA-seq libraries.

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

Interplay of phytohormones and epigenetic regulation: A recipe for plant development and plasticity.

Both phytohormone signaling and epigenetic mechanisms have long been known to play crucial roles in plant development and plasticity in response to ambient stimuli. Indeed, diverse signaling pathways mediated by phytohormones and epigenetic processes integrate multiple upstream signals to regulate various plant traits. Emerging evidence indicates that phytohormones and epigenetic processes interact at multiple levels. In this review, we summarize the current knowledge of the interplay between phytohormones and epigenetic processes from the perspective of phytohormone biology. We also review chemical regulators used in epigenetic studies and propose strategies for developing novel regulators using multidisciplinary approaches.

MCIBox: a toolkit for single-molecule multi-way chromatin interaction visualization and micro-domains identification.

The emerging ligation-free three-dimensional (3D) genome mapping technologies can identify multiplex chromatin interactions with single-molecule precision. These technologies not only offer new insight into high-dimensional chromatin organization and gene regulation, but also introduce new challenges in data visualization and analysis. To overcome these challenges, we developed MCIBox, a toolkit for multi-way chromatin interaction (MCI) analysis, including a visualization tool and a platform for identifying micro-domains with clustered single-molecule chromatin complexes. MCIBox is based on various clustering algorithms integrated with dimensionality reduction methods that can display multiplex chromatin interactions at single-molecule level, allowing users to explore chromatin extrusion patterns and super-enhancers regulation modes in transcription, and to identify single-molecule chromatin complexes that are clustered into micro-domains. Furthermore, MCIBox incorporates a two-dimensional kernel density estimation algorithm to identify micro-domains boundaries automatically. These micro-domains were stratified with distinctive signatures of transcription activity and contained different cell-cycle-associated genes. Taken together, MCIBox represents an invaluable tool for the study of multiple chromatin interactions and inaugurates a previously unappreciated view of 3D genome structure.

Genome-wide characterization of nascent RNA processing in plants.

Following transcription initiation, RNA polymerase II (Pol II) elongates through the genic region and terminates after the polyadenylation signal. This process is accompanied by splicing, 3' cleavage, and polyadenylation, to eventually form a mature mRNA. Recent advances in short-read and long-read high-throughput sequencing methods have shed light on the global landscape of these co-transcriptional events at nucleotide resolution. In this mini review, we summarize recent developments in genome-wide approaches that broadened our understanding of nascent RNA processing in plants.

An atlas of plant full-length RNA reveals tissue-specific and monocots-dicots conserved regulation of poly(A) tail length.

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.

Induced Mutation in GmCOP1b Enhances the Performance of Soybean under Dense Planting Conditions.

CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) is the key photomorphogenic inhibitor that has been extensively studied in higher plants. Nevertheless, its role has not been documented in the economically important soybean. Here we investigated the functions of two COP1 homologous genes, GmCOP1a and GmCOP1b, by analyzing Gmcop1a and Gmcop1b mutants with indels using CRISPR in soybean. We revealed that, although both genes are required for skotomorphogenesis in the dark, the GmCOP1b gene seems to play a more prominent role than GmCOP1a in promoting stem elongation under normal light conditions. Consistently, the bZIP transcriptional factors STF1/2, which repress stem elongation in soybean, accumulated to the highest level in the Gmcop1a1b double mutant, followed by the Gmcop1b and Gmcop1a mutants. Furthermore, the Gmcop1b mutants showed reduced shade response and enhanced performance under high-density conditions in field trials. Taken together, this study provides essential genetic resources for elucidating functional mechanisms of GmCOP1 and breeding of high yield soybean cultivars for future sustainable agriculture.

A SYBR Gold-based Label-free in vitro Dicing Assay.

In Arabidopsis, DICER-LIKE PROTEIN 3 (DCL3) cuts the substrate pre-siRNA into a product siRNA duplex, encompassing one 23-nt strand and one 24-nt strand. To monitor the separation of the siRNA duplex with only 1-nt difference, we developed this protocol to evaluate the in vitro dicing activity of DCL3. The method can be applied for measuring the lengths of single-stranded RNA separated through denaturing urea polyacrylamide gel electrophoresis (urea PAGE), which are visualized by a label-free fluorescence SYBR Gold, and quantified in a multi-function imager. This label-free method is easy to conduct, has low cost, and lacks the hazard of the traditional radio-labeled method. This method can also be adapted to the other Dicers and small RNAs.

Safeguard DCL2-Dependent 22-nt siRNA generation by DCL1.

MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) are crucial for plant growth and development via mediating post-transcriptional gene silencing. In wild-type Arabidopsis, DICER-LIKE 2 (DCL2)-dependent 22-nt siRNAs are rare, whereas DCL1 and DCL4-dependent 21-nt miRNAs and siRNAs are highly abundant. DCL4 naturally inhibits DCL2 in producing abundant 22-nt siRNAs from endogenous transcripts, but whether DCL1 suppresses endogenous 22-nt siRNA production and the extent of repression are still unknown. Here, we report that DCL1 and DCL2 cleaved both miRNA precursors and coding transcript-derived double-stranded RNAs. In a dcl1 dcl4 double mutant, massive 22-nt siRNAs were produced from endogenous protein-coding genes (genic siRNAs). Compared with wild-type, the 22-nt genic siRNAs derived from the Nitrate Reductase 1 (NIA1), NIA2, DIACYLGLYCEROL ACYLTRANSFERASES 3 (DGAT3), SUPPRESSOR OF MAX2 1-LIKE 5 (SMXL5), and SMXL4 in dcl1 dcl4 increased up to 95%. Our analysis further indicated that the 22-nt genic siRNAs in dcl1 dcl4 were mainly loaded into ARGONAUTE 1 (AGO1) or AGO2. Thus, our results demonstrated that both DCL1 and DCL4 safeguard post-transcriptional gene silencing, preventing the production of DCL2-dependent 22-nt genic siRNAs from disrupting plant growth and development.

Regulation of DNA Methylation During Plant Endosperm Development.

The endosperm is a vital storage tissue in plant seeds. It provides nutrients to the embryos or the seedlings during seed development and germination. Although the genetic information in the endosperm cannot be passed directly to the next generation, its inherited epigenetic marks affect gene expression and its development and, consequently, embryo and seed growth. DNA methylation is a major form of epigenetic modification that can be investigated to understand the epigenome changes during reproductive development. Therefore, it is of great significance to explore the effects of endosperm DNA methylation on crop yield and traits. In this review, we discuss the changes in DNA methylation and the resulting imprinted gene expression levels during plant endosperm development, as well as their effects on seed development.