PhasiRNAs in Plants: Their Biogenesis, Genic Sources, and Roles in Stress Responses, Development, and Reproduction

Local trafficking and long-distance transport of small RNAs in plants

Canonical small RNAs in plants, including miRNAs and siRNAs, are key triggers of RNA interference and regulate nearly every major biological process in plants. To establish systemic silencing, small RNAs undergo both short-distance intracellular trafficking or intercellular communication and long-distance transport from one organ to another, even across parasites or pathogens. This enables the delivery of effector molecules throughout the plant, promoting the spread of gene silencing. Biologically, the spatiotemporal regulation of small RNAs results in gradient distributions within cells or along the direction of organogenesis. Furthermore, the spreading capacity of small RNAs, generated in somatic or nurse cells, can guide target gene silencing in germlines in plants. In this review, we summarize recent advances in understanding the regulation and the functional roles of local trafficking and long-distance transport of plant small RNAs in developmental polarity, the maintenance of cell identity, and with a particular focus, the mechanisms of small RNA movement and delivery between companion cells and gametes in plants. Additionally, we discuss the methods and challenges of monitoring small RNA transport in vivo through live imaging, as well as the potential applications of small RNA transport and delivery in the development of RNA-based pesticides.

Atypical epigenetic and small RNA control of degenerated transposons and their fragments in clonally reproducing Spirodela polyrhiza

A handful of model plants have provided insight into silencing of transposable elements (TEs) through RNA-directed DNA methylation (RdDM). Guided by 24 nt long small-interfering RNAs (siRNAs), this epigenetic regulation installs DNA methylation and histone modifications like H3K9me2, which can be subsequently maintained independently of siRNAs. However, the genome of the clonally propagating duckweed Spirodela polyrhiza (Lemnaceae) has low levels of DNA methylation, very low expression of RdDM components, and near absence of 24 nt siRNAs. Moreover, some genes encoding RdDM factors, DNA methylation maintenance, and RNA silencing mechanisms are missing from the genome. Here, we investigated the distribution of TEs and their epigenetic marks in the Spirodela genome. Although abundant degenerated TEs have largely lost DNA methylation and H3K9me2 is low, they remain marked by the heterochromatin-associated H3K9me1 and H3K27me1 modifications. In contrast, we find high levels of DNA methylation and H3K9me2 in the relatively few intact TEs, which are source of 24 nt siRNAs, like RdDM-controlled TEs in other angiosperms. The data suggest that, potentially as adaptation to vegetative propagation, RdDM extent, silencing components, and targets are different from other angiosperms, preferentially focused on potentially intact TEs. It also provides evidence for heterochromatin maintenance independently of DNA methylation in flowering plants. These discoveries highlight the diversity of silencing mechanisms that exist in plants and the importance of using disparate model species to discover these mechanisms.

The Emerging Applications of Artificial MicroRNA-Mediated Gene Silencing in Plant Biotechnology

Improving crop yield potential is crucial to meet the increasing demands of a rapidly expanding global population in an ever-changing and challenging environment. Therefore, different technological approaches have been proposed over the last decades to accelerate plant breeding. Among them, artificial microRNAs (amiRNAs) represent an innovative tool with remarkable potential to assist plant improvement. MicroRNAs (miRNAs) are a group of endogenous, small (20-24 nucleotides), non-coding RNA molecules that play a crucial role in gene regulation. They are associated with most biological processes of a plant, including reproduction, development, cell differentiation, biotic and abiotic stress responses, metabolism, and plant architecture. In this context, amiRNAs are synthetic molecules engineered to mimic the structure and function of endogenous miRNAs, allowing for the targeted silencing of specific nucleic acids. The current review explores the diverse applications of amiRNAs in plant biology and agriculture, such as the management of infectious agents and pests, the engineering of plant metabolism, and the enhancement of plant resilience to abiotic stress. Moreover, we address future perspectives on plant amiRNA-based gene silencing strategies, highlighting the need for further research to fully comprehend the potential of this technology and to translate its scope toward the widespread adoption of amiRNA-based strategies for plant breeding.

CsMPDB 1.0: An interactive web application for visualizing and exploring the microRNAs and phasiRNAs of tea plant (Camellia sinensis var. sinensis 'Shuchazao')

The microRNAs and phasiRNAs of plant are small non-coding RNAs with important functions through regulating gene expression at the post-transcriptional level. However, identifying miRNAs, phasiRNAs and their target genes from numerous sequencing raw data requires multiple software and command-line operations, which are time-consuming and labor-intensive for non-model plants. Therefore, we present CsMPDB (miRNAs and phasiRNAs database of Camellia sinensis), an interactive web application with multiple analysis modules developed to visualize and explore miRNA and phasiRNA in tea plants based on 259 sRNA-seq samples and 24 degradome-seq samples in NCBI. The source code for the CsMPDB was written in R/shiny. It is compatible, extendable, and portable to be easily set up on different operating systems, and can be accessed at http://myshiny.cpolar.io/CsMPDB. This application plays an important role in accelerating the functional study of sRNAs in the transcriptional regulation of tea plants, and has important reference value for the development of sRNA databases of other species.

An organ-specific transcriptome atlas of Curcuma wenyujin: MicroRNAs, phasiRNAs, and metabolic pathways

Curcuma wenyujin Y. H. Chen et C. Ling (C. wenyujin) is a medicinal plant widely used for clinical treatments. In this study, integrated omics data analyses enabled us to discover the microRNAs (miRNAs) and the phased small interfering RNAs (phasiRNAs) on a transcriptome-wide scale. A total of 186 mature miRNAs and 23 precursors were reported. Besides, 31 miRNAs of 14 families were organ-specifically expressed, and 13 of these miRNAs could perform organ-specific target regulation. More than 80% of the phasiRNA loci were organ-specifically expressed, especially in tubers. In some cases, phasiRNAs with distinct increments, but with accordant organ-specific expression patterns, were generated from a highly overlapped region, indicating that different machineries might be synchronously engaged in phasiRNA processing. Based on the transcriptome assembly, 28 and 56 tuber-specific genes were identified to be involved in alkaloid and terpenoid metabolisms, respectively. Analysis of the enzyme-coding genes of the β-elemene biosynthetic pathway showed that the downstream genes were tuber-specific, while the upstream genes were not. We assumed that the precursor metabolites produced in the other organs might be delivered to the tubers for the final steps of β-elemene biosynthesis. Summarily, our report provided an organ-specific transcriptome atlas of C. wenyujin, which could serve as the basis for the molecular studies on organ development and secondary metabolisms in this plant.

Small Interfering RNAs as Critical Regulators of Plant Life Process: New Perspectives on Regulating the Transcriptomic Machinery

Small interfering RNAs (siRNAs) are a distinct class of regulatory RNAs in plants and animals. Gene silencing by small interfering RNAs is one of the fundamental mechanisms for regulating gene expression. siRNAs are critical regulators during developmental processes. siRNAs have similar structures and functions to small RNAs but are derived from double-stranded RNA and may be involved in directing DNA methylation of target sequences. siRNAs are a less well-studied class than the miRNA group, and researchers continue to identify new classes of siRNAs that appear at specific developmental stages and in particular tissues, revealing a more complex mode of siRNA action than previously thought. This review characterizes the siRNA classes and their biogenesis process and focuses on presenting their known functions in the regulation of plant development and responses to biotic and abiotic stresses. The review also highlights the exciting potential for future research in this field, proposing methods for detecting plant siRNAs and a bioinformatic pathway for identifying siRNAs and their functions.

Transcriptome analysis and functional identification of transfer RNA-derived fragments in grape leaves exposed to UV-C radiation

Transfer RNA-derived fragments (tRFs) are noncoding small RNAs derived from transfer RNAs (tRNAs) in microorganisms, animals and plants. In plants, tRFs are known to respond to environmental stimuli, including heat, oxidative stress and UV radiation; however, their specific functions in horticultural plants, such as grapevine, remain poorly understood. In this study, we used RNA-seq to identify differentially expressed genes (DEGs) in grape leaves exposed to UV-C radiation. A total of 1329 and 8055 of genes were differentially expression after 1 and 6 h of UV-C treatment, respectively. We identified a large number of secondary metabolism-related genes in the DEGs, including genes involved in stilbene and flavonoid biosynthesis. Noticeably, the stilbene biosynthesis-related gene was induced earlier than the other genes in the phenylalanine metabolic pathway. We also conducted small RNA-seq and identified differentially expressed (DE) miRNAs and their targets. To explore whether the tRFs involved in UV-C response, further analysis of the small RNA-seq data revealed 23 down-regulated and 41 up-regulated DE tRFs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that the target genes of these tRFs are involved in multiple biological processing, including hormone signal transduction and metabolite synthesis. To validate the function of tRFs, tRF39 and tRF45 were selected and overexpressed in tobacco leaves, and the expression levels of their target genes were inhibited. Our study suggests that the tRFs may regulate multiple biological processes in response to UV-C exposure in grapevine. Our findings provide a foundation for further elucidating the regulatory mechanisms of tRFs in horticultural crops.

Identification and characterization of DICER-LIKE genes and their roles in Marchantia polymorpha development and salt stress response

DICER-LIKE (DCL) proteins play a central role in plant small RNA (sRNA) biogenesis. The genome of the early land plant Marchantia polymorpha encodes four DCL proteins: MpDCL1a, MpDCL1b, MpDCL3, and MpDCL4. While MpDCL1a, MpDCL3 and MpDCL4 show high similarities to their orthologs in Physcomitrium patens and Arabidopsis thaliana, MpDCL1b shares only a limited homology with PpDCL1b, but it is very similar, in terms of functional domains, to orthologs in other moss and fern species. We generated Mpdclge mutant lines for all MpDCL genes with the CRISPR/Cas9 system and conducted phenotypic analyses under control, salt stress, and phytohormone treatments to uncover specific MpDCL functions. The mutants displayed severe developmental aberrations, altered responses to salt and phytohormones, and disturbed sexual organ development. By combining mRNA and sRNA analyses, we demonstrate that MpDCLs and their associated sRNAs play pivotal roles in regulating development, abiotic stress tolerance and phytohormone response in M. polymorpha. We identified MpDCL1a in microRNA biogenesis, MpDCL4 in trans-acting small interfering RNA generation, and MpDCL3 in the regulation of pathogen-related genes. Notably, salt sensitivity in M. polymorpha is dependent on MpDCL1b and Mpdcl1bge mutants display enhanced tolerance and reduced miRNA expression in response to salt stress. We propose that M. polymorpha employs specific mechanisms for regulating MpDCL1b associated miRNAs under high salinity conditions, potentially shared with other species harboring MpDCL1b homologs.

DH2-dependent trans-acting siRNAs regulate leaf and lemma development in rice

In most crops, the development of lateral organs such as leaves and floral organs plays important roles in the architecture of plants and grains and then determines the yield. Establishment of polarity in these lateral organs is one of the most critical events for their morphogenesis. However, the molecular mechanisms regarding this in rice remain unclear. Here, we isolated two allelic mutants named degenerated hull 2-1 and degenerated hull 2-2 (dh2-1 and dh2-2) in rice, exhibiting abaxially rolled leaves and rod-shaped lemmas. DH2 encoded the relatively conservative ARGONAUT 7 (AGO7) protein in plants and expressed in the lateral organs including leaf and floral organs. In addition, the overexpression lines of DH2 showed adaxially rolled leaves. Next, it was proved that DH2 was involved in the synthesis of tasiR-ARFs, the expression level of which was decreased sharply in lateral organs of dh2 mutants. Then, it was found that the expression of OsARF2, OsARF3, OsARF14, and OsARF15, the potential targets of tasiR-ARFs, was increased in lateral organs of dh2 mutants. However, it was not expected that the results of in situ hybridization showed that the four ARF genes were not expressed in WT lemma, whereas they were all ectopically expressed in rod-shaped lemma in dh2 mutants. Meanwhile, tasiR-ARFs were expressed in the whole lemma but not in the abaxial side. This means that there was no opposite expression of tasiR-ARFs and ARFs in the adaxial-abaxial of lemma. Therefore, according to our data, we believe that the pathway of OsAGO7-tasiR-ARFs in rice was more likely involved in the development of the whole lemma, not only the abaxial side, by restricting the ectopic expression of OsARFs in the whole lemma, which was different from that in the lateral organs of Arabidopsis.

Epigenetics in the modern era of crop improvements

Epigenetic mechanisms are integral to plant growth, development, and adaptation to environmental stimuli. Over the past two decades, our comprehension of these complex regulatory processes has expanded remarkably, producing a substantial body of knowledge on both locus-specific mechanisms and genome-wide regulatory patterns. Studies initially grounded in the model plant Arabidopsis have been broadened to encompass a diverse array of crop species, revealing the multifaceted roles of epigenetics in physiological and agronomic traits. With recent technological advancements, epigenetic regulations at the single-cell level and at the large-scale population level are emerging as new focuses. This review offers an in-depth synthesis of the diverse epigenetic regulations, detailing the catalytic machinery and regulatory functions. It delves into the intricate interplay among various epigenetic elements and their collective influence on the modulation of crop traits. Furthermore, it examines recent breakthroughs in technologies for epigenetic modifications and their integration into strategies for crop improvement. The review underscores the transformative potential of epigenetic strategies in bolstering crop performance, advocating for the development of efficient tools to fully exploit the agricultural benefits of epigenetic insights.

Small RNA sequencing analysis provides novel insights into microRNA-mediated regulation of defense responses in chickpea against Fusarium wilt infection

MAIN CONCLUSION

Small RNA sequencing analysis in two chickpea genotypes, JG 62 (Fusarium wilt-susceptible) and WR 315 (Fusarium wilt-resistant), under Fusarium wilt stress led to identification of 544 miRNAs which included 406 known and 138 novel miRNAs. A total of 115 miRNAs showed differential expression in both the genotypes across different combinations. A miRNA, Car-miR398 targeted copper chaperone for superoxide dismutase (CCS) that, in turn, regulated superoxide dismutase (SOD) activity during chickpea-Foc interaction. Fusarium wilt (FW) of chickpea (Cicer arietinum L.) caused by Fusarium oxysporum f. sp. ciceris (Foc) is a destructive soil-borne disease that severely reduces the chickpea yield and quality globally. In the present study, we have investigated microRNAs and the microRNA/target gene crosstalk involved in chickpea resistance to FW. The control and stress samples from two genotypes, JG 62 (FW-susceptible) and WR 315 (FW-resistant), collected at 10 days post-inoculation (dpi), were selected for small RNA sequencing. A total of 12 libraries were constructed and sequenced using Illumina HiSeq 2500 platform. The sequencing and in silico analyses revealed the identification of 544 miRNAs which included 406 known and 138 novel miRNAs. A total of 50 miRNAs were physically co-localized with Foc-resistance QTLs present on chromosome 2 (also known as Foc hotspot). A total of 115 miRNAs showed differential expression in both the genotypes across different combinations. Prediction and functional annotation of miRNA targets revealed their role in transcription regulation, disease resistance, defense response, metabolism, etc. Ten miRNAs and their targets were validated using poly(A)-based qRT-PCR in two genotypes grown under lab and field conditions. Many miRNAs and their targets showed genotype-specific expression. The expression profiling also highlighted, both, similar and different expression patterns for the same sets of miRNA and mRNA at different stages of Foc infection. A high correlation in expression patterns of the miRNAs and their targets in lab- and field-grown plant samples was observed. Interestingly, Car-miR398 targeted copper chaperone for superoxide dismutase (CCS) that, in turn, regulated superoxide dismutase (SOD) activity during chickpea-Foc interaction. The cleavage site in targets was mapped for three miRNAs by analyzing publicly available degradome data for chickpea. The study, for the first time, provides novel insights into microRNA-mediated regulation of resistance and susceptibility mechanisms in chickpea against FW and opens up avenues for the development of the wilt-resistant cultivars in chickpea.

RNA Metabolism and the Role of Small RNAs in Regulating Multiple Aspects of RNA Metabolism

RNA metabolism is focused on RNA molecules and encompasses all the crucial processes an RNA molecule may or will undergo throughout its life cycle. It is an essential cellular process that allows all cells to function effectively. The transcriptomic landscape of a cell is shaped by the processes such as RNA biosynthesis, maturation (RNA processing, folding, and modification), intra- and inter-cellular transport, transcriptional and post-transcriptional regulation, modification, catabolic decay, and retrograde signaling, all of which are interconnected and are essential for cellular RNA homeostasis. In eukaryotes, sRNAs, typically 20-31 nucleotides in length, are a class of ncRNAs found to function as nodes in various gene regulatory networks. sRNAs are known to play significant roles in regulating RNA population at the transcriptional, post-transcriptional, and translational levels. Along with sRNAs, such as miRNAs, siRNAs, and piRNAs, new categories of ncRNAs, i.e., lncRNAs and circRNAs, also contribute to RNA metabolism regulation in eukaryotes. In plants, various genetic screens have demonstrated that sRNA biogenesis mutants, as well as RNA metabolism pathway mutants, exhibit similar growth and development defects, misregulated primary and secondary metabolism, as well as impaired stress response. In addition, sRNAs are both the "products" and the "regulators" in broad RNA metabolism networks; gene regulatory networks involving sRNAs form autoregulatory loops that affect the expression of both sRNA and the respective target. This review examines the interconnected aspects of RNA metabolism with sRNA regulatory pathways in plants. It also explores the potential conservation of these pathways across different kingdoms, particularly in plants and animals. Additionally, the review highlights how cellular RNA homeostasis directly impacts adaptive responses to environmental changes as well as different developmental aspects in plants.

Ribosome binding of phasiRNA precursors accelerates the 24-nt phasiRNA burst in meiotic maize anthers

Reproductive phasiRNAs (phased, secondary, small interfering RNAs), produced from numerous PHAS loci, are essential for plant anther development. PHAS transcripts are enriched on endoplasmic reticulum-bound ribosomes in maize (Zea mays), but the impact of ribosome binding on phasiRNA biogenesis remains elusive. Through ribosome profiling of maize anthers at 10 developmental stages, we demonstrated that 24-PHAS transcripts are bound by ribosomes, with patterns corresponding to the timing and abundance of 24-PHAS transcripts. Ribosome binding to 24-PHAS transcripts is conserved among different maize inbred lines, with ribosomes enriched upstream of the miR2275 target sites. We detected short open reading frames (sORFs) in the ribosome-binding regions of some 24-PHAS transcripts and observed a 3-nt periodicity in most sORFs, but mass spectrometry failed to detect peptides corresponding to the sORFs. Deletion of the entire ribosome-binding region of 24PHAS_NO296 locus eliminated ribosome binding and decreased 24-nt phasiRNA production, without affecting 24PHAS_NO296 transcript levels. In contrast, disrupting only the sORFs in 24PHAS_NO296 did not substantially affect the generation of 24-nt phasiRNAs. A newly formed sORF in these mutants may have re-directed ribosome binding to its transcripts. Overall, these findings demonstrate that sORFs facilitate ribosome binding to 24-PHAS transcripts, thereby promoting phasiRNA biogenesis in meiotic anthers.

Comparative RNA profiling identifies stage-specific phasiRNAs and coexpressed Argonaute genes in Bambusoideae and Pooideae species

Phased, small interfering RNAs (PhasiRNAs) play a crucial role in supporting male fertility in grasses. Earlier work in maize (Zea mays) and rice (Oryza sativa)-and subsequently many other plant species-identified premeiotic 21-nucleotide (nt) and meiotic 24-nt phasiRNAs. More recently, a group of premeiotic 24-nt phasiRNAs was discovered in the anthers of 2 Pooideae species, barley (Hordeum vulgare) and bread wheat (Triticum aestivum). Whether premeiotic 24-nt phasiRNAs and other classes of reproductive phasiRNAs are conserved across Pooideae species remains unclear. We conducted comparative RNA profiling of 3 anther stages in 6 Pooideae species and 1 Bambusoideae species. We observed complex temporal accumulation patterns of 21-nt and 24-nt phasiRNAs in Pooideae and Bambusoideae grasses. In Bambusoideae, 21-nt phasiRNAs accumulated during meiosis, whereas 24-nt phasiRNAs were present in both premeiotic and postmeiotic stages. We identified premeiotic 24-nt phasiRNAs in all 7 species examined. These phasiRNAs exhibit distinct biogenesis mechanisms and potential Argonaute effectors compared to meiotic 24-nt phasiRNAs. We show that specific Argonaute genes coexpressed with stage-specific phasiRNAs are conserved across Bambusoideae and Pooideae species. Our degradome analysis identified a set of conserved miRNA target genes across species, while 21-nt phasiRNA targets were species-specific. Cleavage of few targets was observed for 24-nt phasiRNAs. In summary, this study demonstrates that premeiotic 24-nt phasiRNAs are present across Bambusoideae and Pooideae families, and the temporal accumulation of other classes of 21-nt and 24-nt phasiRNA differs between bamboo and Pooideae species. Furthermore, targets of the 3 classes of phasiRNAs may be rapidly evolving or undetectable.

Unlocking the small RNAs: local and systemic modulators for advancing agronomic enhancement

Small regulatory RNAs (sRNAs) are essential regulators of gene expression across a wide range of organisms to precisely modulate gene activity based on sequence-specific recognition. In model plants like Arabidopsis thaliana, extensive research has primarily concentrated on 21 to 24-nucleotide (nt) sRNAs, particularly microRNAs (miRNAs). Recent advancements in cell and tissue isolation techniques, coupled with advanced sequencing technologies, are revealing a diverse array of preciously uncharacterized sRNA species. These include previously novel structural RNA fragments as well as numerous cell- and tissue-specific sRNAs that are active during distinct developmental stages, thereby enhancing our understanding of the precise and dynamic regulatory roles of sRNAs in plant development regulation. Additionally, a notable feature of sRNAs is their capacity for amplification and movement between cells and tissues, which facilitates long-distance communication-an adaptation critical to plants due to their sessile nature. In this review, we will discuss the classification and mechanisms of action of sRNAs, using legumes as a primary example due to their essential engagement for the unique organ establishment of root nodules and long-distance signaling, and further illustrating the potential applications of sRNAs in modern agricultural breeding and environmentally sustainable plant protection strategies.

Functions and applications of RNA interference and small regulatory RNAs

Small regulatory RNAs play a variety of crucial roles in eukaryotes, influencing gene regulation, developmental timing, antiviral defense, and genome integrity via a process termed RNA interference (RNAi). This process involves Argonaute/small RNA (AGO/sRNA) complexes that target transcripts via sequence complementarity and modulate gene expression and epigenetic modifications. RNAi is a highly conserved gene regulatory phenomenon that recognizes self- and non-self nucleic acids, thereby defending against invasive sequences. Since its discovery, RNAi has been widely applied in functional genomic studies and a range of practical applications. In this review, we focus on the current understanding of the biological roles of the RNAi pathway in transposon silencing, fertility, developmental regulation, immunity, stress responses, and acquired transgenerational inheritance. Additionally, we provide an overview of the applications of RNAi technology in biomedical research, agriculture, and therapeutics.

Integrating multiregulatory analysis reveals the negative regulatory function of miR482a in the response of poplar to canker pathogen infection

Canker disease caused by the bacterium Lonsdalea populi is one of the most destructive diseases affecting poplar stems. However, the detailed stress response mechanisms of poplar have not been widely characterized. To explore the diverse regulatory RNA landscape and the function of key regulators in poplar subjected to L. populi stress, we integrated time-course experiment with mock-inoculation (CK) and inoculation (IN) with L. populi at the first, third, and sixth day (IN1, IN3, IN6) on Populus × euramericana cv. '74/76' (107), small RNA-seq, whole transcriptome-wide analysis, degradome analysis and transgenic experiments. A total of 98 differentially expressed (DE) miRNA, 17 974 DEmRNA, and 807 DElncRNA were identified in poplar infected by L. populi, presenting dynamic changes over the infection course. Regulatory networks among RNAs were further constructed. Notably, a network centered on ptc-miR482a in CK-vs-IN3 contained most DEGs. We show that miR482a and miR1448 are located in one transcript as a polycistron. Overexpression of pre-miR482a-miR1448 (OX482-1448) and pre-miR482a (OX482) increased poplar susceptibility to canker pathogen with reduced accumulation of reactive oxygen species, while the suppression of miR482a (STTM482) conferred poplar disease resistance. PHA7 was validated as the target of miR482a with degradome sequencing and tobacco transient co-transformation, its expression being downregulated in OX482-1448 and OX482 lines. Additionally, a series of phasiRNAs were triggered by miR482a targeting PHA7, forming regulatory cascades with more RLP, NBS-LRR, and PK genes, further verifying the defense function of miR482a. These findings provide insights for understanding the roles of ncRNAs and regulatory networks involved in poplar immunity.

The biogenesis, regulation and functions of transitive siRNA in plants

Small RNA (sRNA)-mediated RNA interference (RNAi) is a sequence-specific gene silencing mechanism that modulates gene expression in eukaryotes. As core molecules of RNAi, various sRNAs are encoded in the plant genome or derived from invading RNA molecules, and their biogenesis depends on distinct genetic pathways. Transitive small interfering RNAs (siRNAs), which are sRNAs produced from double-strand RNA (dsRNA) in a process that depends on RNA-dependent RNA polymerases (RDRs), can amplify and spread silencing signals to additional transcripts, thereby enabling a phenomenon termed "transitive RNAi". Members of this class of siRNAs function in various biological processes ranging from development to stress adaptation. In Arabidopsis thaliana, two RDRs participate in the generation of transitive siRNAs, acting cooperatively with various siRNA generation-related factors, such as the RNA-induced silencing complex (RISC) and aberrant RNAs. Transitive siRNAs are produced in diverse subcellular locations and structures under the control of various mechanisms, highlighting the intricacies of their biogenesis and functions. In this review, we discuss recent advances in understanding the molecular events of transitive siRNA biogenesis and its regulation, with a particular focus on factors involved in RDR recruitment. We aim to provide a comprehensive description of the generalized mechanism governing the biogenesis of transitive siRNAs. Additionally, we present an overview of the diverse biological functions of these siRNAs and raise some pressing questions in this area for further investigation.

Asymmetric bulges within hairpin RNA transgenes influence small RNA size, secondary siRNA production and viral defence

Small RNAs (sRNAs) are essential for normal plant development and range in size classes of 21-24 nucleotides. The 22nt small interfering RNAs (siRNAs) and miRNAs are processed by Dicer-like 2 (DCL2) and DCL1 respectively and can initiate secondary siRNA production from the target transcript. 22nt siRNAs are under-represented due to competition between DCL2 and DCL4, while only a small number of 22nt miRNAs exist. Here we produce abundant 22nt siRNAs and other siRNA size classes using long hairpin RNA (hpRNA) transgenes. By introducing asymmetric bulges into the antisense strand of hpRNA, we shifted the dominant siRNA size class from 21nt of the traditional hpRNA to 22, 23 and 24nt of the asymmetric hpRNAs. The asymmetric hpRNAs effectively silenced a β-glucuronidase (GUS) reporter transgene and the endogenous ethylene insensitive-2 (EIN2) and chalcone synthase (CHS) genes. Furthermore, plants containing the asymmetric hpRNA transgenes showed increased amounts of 21nt siRNAs downstream of the hpRNA target site compared to plants with the traditional hpRNA transgenes. This indicates that these asymmetric hpRNAs are more effective at inducing secondary siRNA production to amplify silencing signals. The 22nt asymmetric hpRNA constructs enhanced virus resistance in plants compared to the traditional hpRNA constructs.

Proteolytic control of the RNA silencing machinery

Studies in plants were often pioneering in the field of RNA silencing and revealed a broad range of small RNA (sRNA) categories. When associated with ARGONAUTE (AGO) proteins, sRNAs play important functions in development, genome integrity, stress responses, and antiviral immunity. Today, most of the protein factors required for the biogenesis of sRNA classes, their amplification through the production of double-stranded RNA, and their function in transcriptional and posttranscriptional regulation have been identified. Nevertheless, and despite the importance of RNA silencing, we still know very little about their posttranslational regulation. This is in stark contrast with studies in metazoans, where different modifications such as prolyl hydroxylation, phosphorylation, sumoylation, ubiquitylation, and others have been reported to alter the activity and stability of key factors, such as AGO proteins. Here, we review current knowledge of how key components of the RNA silencing machinery in plants are regulated during development and by microbial hijacking of endogenous proteases.

A telomere-to-telomere cotton genome assembly reveals centromere evolution and a Mutator transposon-linked module regulating embryo development

Assembly of complete genomes can reveal functional genetic elements missing from draft sequences. Here we present the near-complete telomere-to-telomere and contiguous genome of the cotton species Gossypium raimondii. Our assembly identified gaps and misoriented or misassembled regions in previous assemblies and produced 13 centromeres, with 25 chromosomal ends having telomeres. In contrast to satellite-rich Arabidopsis and rice centromeres, cotton centromeres lack phased CENH3 nucleosome positioning patterns and probably evolved by invasion from long terminal repeat retrotransposons. In-depth expression profiling of transposable elements revealed a previously unannotated DNA transposon (MuTC01) that interacts with miR2947 to produce trans-acting small interfering RNAs (siRNAs), one of which targets the newly evolved LEC2 (LEC2b) to produce phased siRNAs. Systematic genome editing experiments revealed that this tripartite module, miR2947-MuTC01-LEC2b, controls the morphogenesis of complex folded embryos characteristic of Gossypium and its close relatives in the cotton tribe. Our study reveals a trans-acting siRNA-based tripartite regulatory pathway for embryo development in higher plants.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

Genome-wide identification of gene families related to miRNA biogenesis in Mangifera indica L. and their possible role during heat stress

Mango is a popular tropical fruit that requires quarantine hot water treatment (QHWT) for postharvest sanitation, which can cause abiotic stress. Plants have various defense mechanisms to cope with stress; miRNAs mainly regulate the expression of these defense responses. Proteins involved in the biogenesis of miRNAs include DICER-like (DCL), ARGONAUTE (AGO), HYPONASTIC LEAVES 1 (HYL1), SERRATE (SE), HUA ENHANCER1 (HEN1), HASTY (HST), and HEAT-SHOCK PROTEIN 90 (HSP90), among others. According to our analysis, the mango genome contains five DCL, thirteen AGO, six HYL, two SE, one HEN1, one HST, and five putative HSP90 genes. Gene structure prediction and domain identification indicate that sequences contain key domains for their respective gene families, including the RNase III domain in DCL and PAZ and PIWI domains for AGOs. In addition, phylogenetic analysis indicates the formation of clades that include the mango sequences and their respective orthologs in other flowering plant species, supporting the idea these are functional orthologs. The analysis of cis-regulatory elements of these genes allowed the identification of MYB, ABRE, GARE, MYC, and MeJA-responsive elements involved in stress responses. Gene expression analysis showed that most genes are induced between 3 to 6 h after QHWT, supporting the early role of miRNAs in stress response. Interestingly, our results suggest that mango rapidly induces the production of miRNAs after heat stress. This research will enable us to investigate further the regulation of gene expression and its effects on commercially cultivated fruits, such as mango, while maintaining sanitary standards.

Enigmatic role of auxin response factors in plant growth and stress tolerance

Abiotic and biotic stresses globally constrain plant growth and impede the optimization of crop productivity. The phytohormone auxin is involved in nearly every aspect of plant development. Auxin acts as a chemical messenger that influences gene expression through a short nuclear pathway, mediated by a family of specific DNA-binding transcription factors known as Auxin Response Factors (ARFs). ARFs thus act as effectors of auxin response and translate chemical signals into the regulation of auxin responsive genes. Since the initial discovery of the first ARF in Arabidopsis, advancements in genetics, biochemistry, genomics, and structural biology have facilitated the development of models elucidating ARF action and their contributions to generating specific auxin responses. Yet, significant gaps persist in our understanding of ARF transcription factors despite these endeavors. Unraveling the functional roles of ARFs in regulating stress response, alongside elucidating their genetic and molecular mechanisms, is still in its nascent phase. Here, we review recent research outcomes on ARFs, detailing their involvement in regulating leaf, flower, and root organogenesis and development, as well as stress responses and their corresponding regulatory mechanisms: including gene expression patterns, functional characterization, transcriptional, post-transcriptional and post- translational regulation across diverse stress conditions. Furthermore, we delineate unresolved questions and forthcoming challenges in ARF research.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function

Reproductive phasiRNAs (phased, small interfering RNAs) are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt (nucleotides) phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at the premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 (Dicer-like 5) for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely i) not triggered by microRNAs, ii) not loaded by AGO18 proteins, and iii) not capable of mediating PHAS precursor cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

Regulatory microRNAs and phasiRNAs of paclitaxel biosynthesis in Taxus chinensis

Paclitaxel (trade name Taxol) is a rare diterpenoid with anticancer activity isolated from Taxus. At present, paclitaxel is mainly produced by the semi-synthetic method using extract of Taxus tissues as raw materials. The studies of regulatory mechanisms in paclitaxel biosynthesis would promote the production of paclitaxel through tissue/cell culture approaches. Here, we systematically identified 990 transcription factors (TFs), 460 microRNAs (miRNAs), and 160 phased small interfering RNAs (phasiRNAs) in Taxus chinensis to explore their interactions and potential roles in regulation of paclitaxel synthesis. The expression levels of enzyme genes in cone and root were higher than those in leaf and bark. Nearly all enzyme genes in the paclitaxel synthesis pathway were significantly up-regulated after jasmonate treatment, except for GGPPS and CoA Ligase. The expression level of enzyme genes located in the latter steps of the synthesis pathway was significantly higher in female barks than in male. Regulatory TFs were inferred through co-expression network analysis, resulting in the identification of TFs from diverse families including MYB and AP2. Genes with ADP binding and copper ion binding functions were overrepresented in targets of miRNA genes. The miRNA targets were mainly enriched with genes in plant hormone signal transduction, mRNA surveillance pathway, cell cycle and DNA replication. Genes in oxidoreductase activity, protein-disulfide reductase activity were enriched in targets of phasiRNAs. Regulatory networks were further constructed including components of enzyme genes, TFs, miRNAs, and phasiRNAs. The hierarchical regulation of paclitaxel production by miRNAs and phasiRNAs indicates a robust regulation at post-transcriptional level. Our study on transcriptional and posttranscriptional regulation of paclitaxel synthesis provides clues for enhancing paclitaxel production using synthetic biology technology.

Premeiotic 24-nt phasiRNAs are present in the Zea genus and unique in biogenesis mechanism and molecular function

Reproductive phasiRNAs are broadly present in angiosperms and play crucial roles in sustaining male fertility. While the premeiotic 21-nt phasiRNAs and meiotic 24-nt phasiRNA pathways have been extensively studied in maize (Zea mays) and rice (Oryza sativa), a third putative category of reproductive phasiRNAs-named premeiotic 24-nt phasiRNAs-have recently been reported in barley (Hordeum vulgare) and wheat (Triticum aestivum). To determine whether premeiotic 24-nt phasiRNAs are also present in maize and related species and begin to characterize their biogenesis and function, we performed a comparative transcriptome and degradome analysis of premeiotic and meiotic anthers from five maize inbred lines and three teosinte species/subspecies. Our data indicate that a substantial subset of the 24-nt phasiRNA loci in maize and teosinte are already highly expressed at premeiotic phase. The premeiotic 24-nt phasiRNAs are similar to meiotic 24-nt phasiRNAs in genomic origin and dependence on DCL5 for biogenesis, however, premeiotic 24-nt phasiRNAs are unique in that they are likely (i) not triggered by microRNAs, (ii) not loaded by AGO18 proteins, and (iii) not capable of mediating cis-cleavage. In addition, we also observed a group of premeiotic 24-nt phasiRNAs in rice using previously published data. Together, our results indicate that the premeiotic 24-nt phasiRNAs constitute a unique class of reproductive phasiRNAs and are present more broadly in the grass family (Poaceae) than previously known.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

MicroRNA482/2118 is lineage-specifically involved in gibberellin signalling via the regulation of GID1 expression by targeting noncoding PHAS genes and subsequently instigated phasiRNAs

MicroRNA482/2118 (miR482/2118) is a 22-nt miRNA superfamily, with conserved functions in disease resistance and plant development. It usually instigates the production of phased small interfering RNAs (phasiRNAs) from its targets to expand or reinforce its silencing effect. Using a new high-quality reference genome sequence and comprehensive small RNA profiling, we characterized a newly evolved regulatory pathway of miR482/2118 in litchi. In this pathway, miR482/2118 cleaved a novel noncoding trans-acting gene (LcTASL1) and triggered phasiRNAs to regulate the expression of gibberellin (GA) receptor gene GIBBERELLIN INSENSITIVE DWARF1 (GID1) in trans; another trans-acting gene LcTASL2, targeted by LcTASL1-derived phasiRNAs, produced phasiRNAs as well to target LcGID1 to reinforce the silencing effect of LcTASL1. We found this miR482/2118-TASL-GID1 pathway was likely involved in fruit development, especially the seed development in litchi. In vivo construction of the miR482a-TASL-GID1 pathway in Arabidopsis could lead to defects in flower and silique development, analogous to the phenotype of gid1 mutants. Finally, we found that a GA-responsive transcription factor, LcGAMYB33, could regulate LcMIR482/2118 as a feedback mechanism of the sRNA-silencing pathway. Our results deciphered a lineage-specifically evolved regulatory module of miR482/2118, demonstrating the high dynamics of miR482/2118 function in plants.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

The unique dual targeting of AGO1 by two types of PRMT enzymes promotes phasiRNA loading in Arabidopsis thaliana

Arginine/R methylation (R-met) of proteins is a widespread post-translational modification (PTM), deposited by a family of protein arginine/R methyl transferase enzymes (PRMT). Regulations by R-met are involved in key biological processes deeply studied in metazoan. Among those, post-transcriptional gene silencing (PTGS) can be regulated by R-met in animals and in plants. It mainly contributes to safeguard processes as protection of genome integrity in germlines through the regulation of piRNA pathway in metazoan, or response to bacterial infection through the control of AGO2 in plants. So far, only PRMT5 has been identified as the AGO/PIWI R-met writer in higher eukaryotes. We uncovered that AGO1, the main PTGS effector regulating plant development, contains unique R-met features among the AGO/PIWI superfamily, and outstanding in eukaryotes. Indeed, AGO1 contains both symmetric (sDMA) and asymmetric (aDMA) R-dimethylations and is dually targeted by PRMT5 and by another type I PRMT in Arabidopsis thaliana. We showed also that loss of sDMA didn't compromise AtAGO1 subcellular trafficking in planta. Interestingly, we underscored that AtPRMT5 specifically promotes the loading of phasiRNA in AtAGO1. All our observations bring to consider this dual regulation of AtAGO1 in plant development and response to environment, and pinpoint the complexity of AGO1 post-translational regulation.

High-quality assembled and annotated genomes of Nicotiana tabacum and Nicotiana benthamiana reveal chromosome evolution and changes in defense arsenals

Nicotiana tabacum and Nicotiana benthamiana are widely used models in plant biology research. However, genomic studies of these species have lagged. Here we report the chromosome-level reference genome assemblies for N. benthamiana and N. tabacum with an estimated 99.5% and 99.8% completeness, respectively. Sensitive transcription start and termination site sequencing methods were developed and used for accurate gene annotation in N. tabacum. Comparative analyses revealed evidence for the parental origins and chromosome structural changes, leading to hybrid genome formation of each species. Interestingly, the antiviral silencing genes RDR1, RDR6, DCL2, DCL3, and AGO2 were lost from one or both subgenomes in N. benthamiana, while both homeologs were kept in N. tabacum. Furthermore, the N. benthamiana genome encodes fewer immune receptors and signaling components than that of N. tabacum. These findings uncover possible reasons underlying the hypersusceptible nature of N. benthamiana. We developed the user-friendly Nicomics (http://lifenglab.hzau.edu.cn/Nicomics/) web server to facilitate better use of Nicotiana genomic resources as well as gene structure and expression analyses.

The Emerging Role of Non-Coding RNAs (ncRNAs) in Plant Growth, Development, and Stress Response Signaling

Plant species utilize a variety of regulatory mechanisms to ensure sustainable productivity. Within this intricate framework, numerous non-coding RNAs (ncRNAs) play a crucial regulatory role in plant biology, surpassing the essential functions of RNA molecules as messengers, ribosomal, and transfer RNAs. ncRNAs represent an emerging class of regulators, operating directly in the form of small interfering RNAs (siRNAs), microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). These ncRNAs exert control at various levels, including transcription, post-transcription, translation, and epigenetic. Furthermore, they interact with each other, contributing to a variety of biological processes and mechanisms associated with stress resilience. This review primarily concentrates on the recent advancements in plant ncRNAs, delineating their functions in growth and development across various organs such as root, leaf, seed/endosperm, and seed nutrient development. Additionally, this review broadens its scope by examining the role of ncRNAs in response to environmental stresses such as drought, salt, flood, heat, and cold in plants. This compilation offers updated information and insights to guide the characterization of the potential functions of ncRNAs in plant growth, development, and stress resilience in future research.

Non-coding RNAs and leaf senescence: Small molecules with important roles

Non-coding RNAs (ncRNAs) are a special class of functional RNA molecules that are not translated into proteins. ncRNAs have emerged as pivotal regulators of diverse developmental processes in plants. Recent investigations have revealed the association of ncRNAs with the regulation of leaf senescence, a complex and tightly regulated developmental process. However, a comprehensive review of the involvement of ncRNAs in the regulation of leaf senescence is still lacking. This manuscript aims to summarize the molecular mechanisms underlying ncRNAs-mediated leaf senescence and the potential applications of ncRNAs to manipulate the onset and progression of leaf senescence. Various classes of ncRNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), are discussed in terms of their regulatory mechanisms in leaf senescence. Furthermore, we explore the interactions between ncRNA and the key regulators of senescence, including transcription factors as well as core components in phytohormone signaling pathways. We also discuss the possible challenges and approaches related to ncRNA-mediated leaf senescence. This review contributes to a further understanding of the intricate regulatory network involving ncRNAs in leaf senescence.

The plant siRNA landscape

Whereas micro (mi)RNAs are considered the clean, noble side of the small RNA world, small interfering (si)RNAs are often seen as a noisy set of molecules whose barbarian acronyms reflect a large diversity of often elusive origins and functions. Twenty-five years after their discovery in plants, however, new classes of siRNAs are still being identified, sometimes in discrete tissues or at particular developmental stages, making the plant siRNA world substantially more complex and subtle than originally anticipated. Focusing primarily on the model Arabidopsis, we review here the plant siRNA landscape, including transposable elements (TE)-derived siRNAs, a vast array of non-TE-derived endogenous siRNAs, as well as exogenous siRNAs produced in response to invading nucleic acids such as viruses or transgenes. We primarily emphasize the extraordinary sophistication and diversity of their biogenesis and, secondarily, the variety of their known or presumed functions, including via non-cell autonomous activities, in the sporophyte, gametophyte, and shortly after fertilization.

Identification of Tomato microRNAs in Late Response to Trichoderma atroviride

The tomato (Solanum lycopersicum) is an important crop worldwide and is considered a model plant to study stress responses. Small RNAs (sRNAs), 21-24 nucleotides in length, are recognized as a conserved mechanism for regulating gene expression in eukaryotes. Plant endogenous sRNAs, such as microRNA (miRNA), have been involved in disease resistance. High-throughput RNA sequencing was used to analyze the miRNA profile of the aerial part of 30-day-old tomato plants after the application of the fungus Trichoderma atroviride to the seeds at the transcriptional memory state. Compared to control plants, ten differentially expressed (DE) miRNAs were identified in those inoculated with Trichoderma, five upregulated and five downregulated, of which seven were known (miR166a, miR398-3p, miR408, miR5300, miR6024, miR6027-5p, and miR9471b-3p), and three were putatively novel (novel miR257, novel miR275, and novel miR1767). miRNA expression levels were assessed using real-time quantitative PCR analysis. A plant sRNA target analysis of the DE miRNAs predicted 945 potential target genes, most of them being downregulated (84%). The analysis of KEGG metabolic pathways showed that most of the targets harbored functions associated with plant-pathogen interaction, membrane trafficking, and protein kinases. Expression changes of tomato miRNAs caused by Trichoderma are linked to plant defense responses and appear to have long-lasting effects.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato

BACKGROUND

Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV).

RESULTS

We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes.

CONCLUSION

We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

The trans-kingdom communication of noncoding RNAs in plant-environment interactions

As conserved regulatory agents, noncoding RNAs (ncRNAs) have an important impact on many aspects of plant life, including growth, development, and environmental response. Noncoding RNAs can travel through not only plasmodesma and phloem but also intercellular barriers to regulate distinct processes. Increasing evidence shows that the intercellular trans-kingdom transmission of ncRNAs is able to modulate many important interactions between plants and other organisms, such as plant response to pathogen attack, the symbiosis between legume plants and rhizobia and the interactions with parasitic plants. In these interactions, plant ncRNAs are believed to be sorted into extracellular vesicles (EVs) or other nonvesicular vehicles to pass through cell barriers and trigger trans-kingdom RNA interference (RNAi) in recipient cells from different species. There is evidence that the features of extracellular RNAs and associated RNA-binding proteins (RBPs) play a role in defining the RNAs to retain in cell or secrete outside cells. Despite the few reports about RNA secretion pathway in plants, the export of extracellular ncRNAs is orchestrated by a series of pathways in plants. The identification and functional analysis of mobile small RNAs (sRNAs) are attracting increasing attention in recent years. In this review, we discuss recent advances in our understanding of the function, sorting, transport, and regulation of plant extracellular ncRNAs.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

Mechanisms that regulate the production of secondary siRNAs in plants

Many organisms produce secondary small interfering RNAs (siRNAs) that are triggered by primary small RNAs to regulate various biological processes. Plants have evolved several types of secondary siRNA biogenesis pathways that play important roles in development, stress responses and defense against viruses and transposons. The critical step of these pathways is the production of double-stranded RNAs by RNA-dependent RNA polymerases. This step is normally tightly regulated, but when its control is released, secondary siRNA production is initiated. In this article, we will review the recent advances in secondary siRNA production triggered by microRNAs encoded in the genome and siRNAs derived from invasive nucleic acids. In particular, we will focus on the factors, events, and RNA/DNA elements that promote or inhibit the early steps of secondary siRNA biogenesis.

The phased short-interfering RNA siRD29(-) regulates GIBBERELLIN 3-OXIDASE 3 during stolon-to-tuber transitions in potato

Phased short-interfering RNAs (phasiRNAs) fine tune various stages of growth, development, and stress responses in plants. Potato (Solanum tuberosum) tuberization is a complex process, wherein a belowground modified stem (stolon) passes through developmental stages like swollen stolon and minituber before it matures to a potato. Previously, we identified several phasiRNA-producing loci (PHAS) from stolon-to-tuber transition stages. However, whether phasiRNAs mediate tuber development remains unknown. Here, we show that a gene encoding NB-ARC DOMAIN-CONTAINING DISEASE RESISTANCE PROTEIN (StRGA4; a PHAS locus) is targeted by Stu-microRNA482c to generate phasiRNAs. Interestingly, we observed that one of the phasiRNAs, referred as short-interfering RNA D29(-), i.e. siRD29(-), targets the gibberellin (GA) biosynthesis gene GIBBERELLIN 3-OXIDASE 3 (StGA3ox3). Since regulation of bioactive GA levels in stolons controls tuber development, we hypothesized that a gene regulatory module, Stu-miR482c-StRGA4-siRD29(-)-StGA3ox3, could govern tuber development. Through transient expression assays and small RNA sequencing, generation of siRD29(-) and its phase was confirmed in planta. Notably, the expression of StGA3ox3 was higher in swollen stolon compared to stolon, whereas siRD29(-) showed a negative association with StGA3ox3 expression. Antisense (AS) lines of StGA3ox3 produced more tubers compared to wild type. As expected, StRGA4 overexpression (OE) lines had high levels of siRD29(-) and mimicked the phenotypes of StGA3ox3-AS lines, indicating the functionality of this module in potato. In vitro tuberization assays (with or without a GA inhibitor) using StGA3ox3 antisense lines and overexpression lines of StGA3ox3 or StRGA4 revealed that StGA3ox3 controls the tuber stalk development. Taken together, our findings suggest that a phasiRNA, siRD29(-), mediates the regulation of StGA3ox3 during stolon-to-tuber transitions in potato.

Beyond Loading: Functions of Plant ARGONAUTE Proteins

ARGONAUTE (AGO) proteins are key components of the RNA-induced silencing complex (RISC) that mediates gene silencing in eukaryotes. Small-RNA (sRNA) cargoes are selectively loaded into different members of the AGO protein family and then target complementary sequences to in-duce transcriptional repression, mRNA cleavage, or translation inhibition. Previous reviews have mainly focused on the traditional roles of AGOs in specific biological processes or on the molecular mechanisms of sRNA sorting. In this review, we summarize the biological significance of canonical sRNA loading, including the balance among distinct sRNA pathways, cross-regulation of different RISC activities during plant development and defense, and, especially, the emerging roles of AGOs in sRNA movement. We also discuss recent advances in novel non-canonical functions of plant AGOs. Perspectives for future functional studies of this evolutionarily conserved eukaryotic protein family will facilitate a more comprehensive understanding of the multi-faceted AGO proteins.

Role of miRNAs in sucrose stress response, reactive oxygen species, and anthocyanin biosynthesis in Arabidopsis thaliana

In plants, sucrose is the main transported disaccharide that is the primary product of photosynthesis and controls a multitude of aspects of the plant life cycle including structure, growth, development, and stress response. Sucrose is a signaling molecule facilitating various stress adaptations by crosstalk with other hormones, but the molecular mechanisms are not well understood. Accumulation of high sucrose concentrations is a hallmark of many abiotic and biotic stresses, resulting in the accumulation of reactive oxygen species and secondary metabolite anthocyanins that have antioxidant properties. Previous studies have shown that several MYeloBlastosis family/MYB transcription factors are positive and negative regulators of sucrose-induced anthocyanin accumulation and subject to microRNA (miRNA)-mediated post-transcriptional silencing, consistent with the notion that miRNAs may be "nodes" in crosstalk signaling by virtue of their sequence-guided targeting of different homologous family members. In this study, we endeavored to uncover by deep sequencing small RNA and mRNA transcriptomes the effects of exogenous high sucrose stress on miRNA abundances and their validated target transcripts in Arabidopsis. We focused on genotype-by-treatment effects of high sucrose stress in Production of Anthocyanin Pigment 1-Dominant/pap1-D, an activation-tagged dominant allele of MYB75 transcription factor, a positive effector of secondary metabolite anthocyanin pathway. In the process, we discovered links to reactive oxygen species signaling through miR158/161/173-targeted Pentatrico Peptide Repeat genes and two novel non-canonical targets of high sucrose-induced miR408 and miR398b*(star), relevant to carbon metabolic fluxes: Flavonoid 3'-Hydroxlase (F3'H), an important enzyme in determining the B-ring hydroxylation pattern of flavonoids, and ORANGE a post-translational regulator of Phytoene Synthase expression, respectively. Taken together, our results contribute to understanding the molecular mechanisms of carbon flux shifts from primary to secondary metabolites in response to high sugar stress.

PhasiHunter: a robust phased siRNA regulatory cascade mining tool based on multiple reference sequences

SUMMARY

In recent years, phased small interfering RNA has been found to play crucial roles in many biological processes in plants. However, efficiently predicting phasiRNA regulatory cascades with computational methods is still challenging. Here, we introduce PhasiHunter, a phasiRNA regulatory network prediction tool that has several distinctive features compared to existing tools: (i) PhasiHunter employs two major phasiRNA prediction algorithms, namely phase score and hypergeometric distribution-based methods, to ensure the integrity and accuracy of prediction; (ii) PhasiHunter can identify phasiRNAs and their regulatory networks based on multiple reference sequences and the predicted results can be automatically integrated; (iii) PhasiHunter can efficiently identify the phasiRNAs generated through alternative splicing events; and (iv) the excellent data structure and parallel computing architecture allow PhasiHunter to predict phasiRNAs and their regulatory pathways with high efficiency.

AVAILABILITY AND IMPLEMENTATION

PhasiHunter is an open-source tool that is available at https://github.com/HuangLab-CBI/PhasiHunter.

Multiomics analysis reveals a link between Brassica-specific miR1885 and rapeseed tolerance to low temperature

Brassica crops include various edible vegetable and plant oil crops, and their production is limited by low temperature beyond their tolerant capability. The key regulators of low-temperature resistance in Brassica remain largely unexplored. To identify posttranscriptional regulators of plant response to low temperature, we performed small RNA profiling, and found that 16 known miRNAs responded to cold treatment in Brassica rapa. The cold response of seven of those miRNAs were further confirmed by qRT-PCR and/or northern blot analyses. In parallel, a genome-wide association study of 220 accessions of Brassica napus identified four candidate MIRNA genes, all of which were cold-responsive, at the loci associated with low-temperature resistance. Specifically, these large-scale data analyses revealed a link between miR1885 and the plant response to low temperature in both B. rapa and B. napus. Using 5' rapid amplification of cDNA ends approach, we validated that miR1885 can cleave its putative target gene transcripts, Bn.TIR.A09 and Bn.TNL.A03, in B. napus. Furthermore, overexpression of miR1885 in Semiwinter type B. napus decreased the mRNA abundance of Bn.TIR.A09 and Bn.TNL.A03 and resulted in increased sensitivity to low temperature. Knocking down of miR1885 in Spring type B. napus led to increased mRNA abundance of its targets and improved rapeseed tolerance to low temperature. Together, our results suggested that the loci of miR1885 and its targets could be potential candidates for the molecular breeding of low temperature-tolerant Spring type Brassica crops.

Evidence of a novel silencing effect on transgenes in the Arabidopsis thaliana sperm cell

We encountered unexpected transgene silencing in Arabidopsis thaliana sperm cells; transgenes encoding proteins with no specific intracellular localization (cytoplasmic proteins) were silenced transcriptionally or posttranscriptionally. The mRNA of cytoplasmic protein transgenes tagged with a fluorescent protein gene was significantly reduced, resulting in undetectable fluorescent protein signals in the sperm cell. Silencing of the cytoplasmic protein transgenes in the sperm cell did not affect the expression of either its endogenous homologous genes or cotransformed transgenes encoding a protein with targeted intracellular localization. This transgene silencing in the sperm cell persisted in mutants of the major gene silencing machinery including DNA methylation. The incomprehensible, yet real, transgene silencing phenotypes occurring in the sperm cell could mislead the interpretation of experimental results in plant reproduction, and this Commentary calls attention to that risk and highlights details of this novel cytoplasmic protein transgene silencing.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.