Roles of DCL4 and DCL3b in rice phased small RNA biogenesis

Broadening the Nicotiana benthamiana research toolbox through the generation of dicer-like mutants using CRISPR/Cas9 approaches

RNA silencing in plants plays a pivotal role in various biological processes, including development, epigenetic modifications and stress response. Key components of this network are Dicer-like (DCL) proteins. Nicotiana benthamiana encodes four DCLs, each responsible for the generation of distinct small RNA (sRNA) populations, which regulate different functions. However, elucidating the precise role of each DCL has been proven challenging, as overlapping functions exist within DCLs. In our present study, we have successfully generated dcl2, dcl3 and dcl4 homozygous mutants, employing two different CRISPR/Cas9 approaches. The first approach is based on a transgene-mediated delivery of the single-guide RNA (sgRNA), while the second approach employs a viral vector for sgRNA delivery. By utilizing a suite of screening techniques, including polymerase chain reaction (PCR), T7 endonuclease I (T7E1) assay, high-resolution melt analysis (HRMA) and DNA sequencing, we successfully generated dcl2, dcl3 and dcl4 homozygous mutants harboring identical mutations in every allele. To evaluate these dcl mutants, we examined their sRNA profiles and phenotypes. We further have indications that homozygous mutations of a gene do not always lead to the desired loss-of-function, highlighting the importance of mutant evaluation. dcl mutants represent invaluable tools to explore how overlapping silencing pathways are connected to essential plant functions, including development, stress responses and pathogen defense. Additionally, they hold potential for biotechnological applications, such as crop improvement and gene silencing tools. We anticipate that our study will make significant contributions to enhance understanding of the role of DCLs in plants.

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.

Decoding the Dialog Between Plants and Arbuscular Mycorrhizal Fungi: A Molecular Genetic Perspective

Arbuscular mycorrhizal (AM) symbiosis, a mutually beneficial interaction between plant roots and AM fungi, plays a key role in plant growth, nutrient acquisition, and stress tolerance, which make it a major focus for sustainable agricultural strategies. This intricate association involves extensive transcriptional reprogramming in host plant cells during the formation of arbuscules, which are specialized fungal structures for nutrient exchange. The symbiosis is initiated by molecular signaling pathways triggered by fungal chitooligosaccharides and strigolactones released by plant roots, which act as chemoattractants and signaling molecules to promote fungal spore germination, colonization, and arbuscule development. Calcium spiking, mediated by LysM domain receptor kinases, serves as a critical second messenger in coordinating fungal infection and intracellular accommodation. GRAS transcription factors are key components that regulate the transcriptional networks necessary for arbuscule development and maintenance, while small RNAs (sRNAs) from both plant and fungi, contribute to modifications in gene expression, including potential bidirectional sRNA exchange to modulate symbiosis. Understanding the molecular mechanisms related to AM symbiosis may provide valuable insights for implementation of strategies related to enhancing plant productivity and resilience.

Profiling of 24-nt phasiRNAs in Rice Anther by Massively Parallel Sequencing and Stem-Loop Reverse Transcription PCR

Small RNAs are highly abundant and play important roles in plant reproduction. Profiling of small RNAs in reproductive tissues is a critical step in understanding their biology. Here, we describe a protocol for small RNA profiling in rice anthers, with a focus on an abundantly expressed but little-understood reproductive small RNA class named 24-nucleotide phased secondary small interfering RNAs (24-nt phasiRNAs). The protocol details small RNA library preparation steps for high-throughput sequencing, with a subsequent bioinformatic analysis framework. A low-throughput stem-loop reverse transcription PCR protocol for quick semi-quantification of small RNAs is also included as a complementary approach. The collective protocol, although focused on 24-nt phasiRNAs, should be amenable to other small RNA classes in various plant reproductive tissues.

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.

Genome-wide identification, characterization and expression analysis of key gene families in RNA silencing in centipedegrass

BACKGROUND

Argonaute (AGO), Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) are essential components of RNA silencing pathways in plants. These components are crucial for the generation and regulatory functions of small RNAs, especially in plant development and response to environmental stresses. Despite their well-characterized functions in other plant species, there is limited information about these genes and their stress responses in centipedegrass (Eremochloa ophiuroides), a key turfgrass species.

RESULTS

Using genome-wide analysis we identified 20 AGO, 6 DCL, and 10 RDR members in centipedegrass and provided a comprehensive overview of their characteristics. We performed the chromosomal location, gene duplication, syntenic analysis, conserve motif, gene structure, and cis-acting elements analysis. And conducted phylogenetic analyses to clarify the evolutionary relationships among the EoAGO, EoDCL, and EoRDR gene families. Three-dimensional modeling prediction of EoAGO, EoDCL, and EoRDR proteins supported the phylogenetic classification. Furthermore, we examined the expression patterns of these genes in different tissues (spike, stem, leaf, root, and flower) and under different stress conditions (cold, salt, drought, aluminum, and herbicide) using RT-qPCR. The results revealed that most of EoAGO, EoDCL, and EoRDR genes were upregulated in response to multiple abiotic stresses, while some exhibited unique responses, suggesting potential specialized regulatory functions.

CONCLUSION

In this study, we performed a comprehensive genome‑wide identification, and phylogenetic and expression pattern analyses of the EoAGO, EoDCL and EoRDR gene families. Our analysis provides a foundation for future research on the RNA silence elements of turfgrass, and affords scientific basis and insights for clarifying the expression patterns of EoAGO, EoDCL and EoRDR genes under adversity stress. Further functional validation and molecular breeding of these genes can be carried out for enhancing the stress resistance of centipedegrass.

Genome-wide identification of phasiRNAs in Arabidopsis thaliana, and insights into biogenesis, temperature sensitivity, and organ specificity

The knowledge of biogenesis and target regulation of the phased small interfering RNAs (phasiRNAs) needs continuous update, since the phasiRNA loci are dynamically evolved in plants. Here, hundreds of phasiRNA loci of Arabidopsis thaliana were identified in distinct tissues and under different temperature. In flowers, most of the 24-nt loci are RNA-dependent RNA polymerase 2 (RDR2)-dependent, while the 21-nt loci are RDR6-dependent. Among the RDR-dependent loci, a significant portion is Dicer-like 1-dependent, indicating the involvement of microRNAs in their expression. Besides, two TAS candidates were discovered. Some interesting features of the phasiRNA loci were observed, such as the strong strand bias of phasiRNA generation, and the capacity of one locus for producing phasiRNAs by different increments. Both organ specificity and temperature sensitivity were observed for phasiRNA expression. In leaves, the TAS genes are highly activated under low temperature. Several trans-acting siRNA-target pairs are also temperature-sensitive. In many cases, the phasiRNA expression patterns correlate well with those of the processing signals. Analysis of the rRNA-depleted degradome uncovered several phasiRNA loci to be RNA polymerase II-independent. Our results should advance the understanding on phasiRNA biogenesis and regulation in plants.

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.

Switching action modes of miR408-5p mediates auxin signaling in rice

MicroRNAs (miRNAs) play fundamental roles in many developmental and physiological processes in eukaryotes. MiRNAs in plants generally regulate their targets via either mRNA cleavage or translation repression; however, which approach plays a major role and whether these two function modes can shift remains elusive. Here, we identify a miRNA, miR408-5p that regulates AUXIN/INDOLE ACETIC ACID 30 (IAA30), a critical repressor in the auxin pathway via switching action modes in rice. We find that miR408-5p usually inhibits IAA30 protein translation, but in a high auxin environment, it promotes the decay of IAA30 mRNA when it is overproduced. We further demonstrate that IDEAL PLANT ARCHITECTURE1 (IPA1), an SPL transcription factor regulated by miR156, mediates leaf inclination through association with miR408-5p precursor promoter. We finally show that the miR156-IPA1-miR408-5p-IAA30 module could be controlled by miR393, which silences auxin receptors. Together, our results define an alternative auxin transduction signaling pathway in rice that involves the switching of function modes by miR408-5p, which contributes to a better understanding of the action machinery as well as the cooperative network of miRNAs in plants.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

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.

Dicer-like 3a mediates intergenerational resistance against root-knot nematodes in rice via hormone responses

In a continuously changing and challenging environment, passing down the memory of encountered stress factors to offspring could provide an evolutionary advantage. In this study, we demonstrate the existence of "intergenerational acquired resistance" in the progeny of rice (Oryza sativa) plants attacked by the belowground parasitic nematode Meloidogyne graminicola. Transcriptome analyses revealed that genes involved in defense pathways are generally downregulated in progeny of nematode-infected plants under uninfected conditions but show a stronger induction upon nematode infection. This phenomenon was termed "spring loading" and depends on initial downregulation by the 24-nucleotide (nt) siRNA biogenesis gene dicer-like 3a (dcl3a) involved in the RNA-directed DNA methylation pathway. Knockdown of dcl3a led to increased nematode susceptibility and abolished intergenerational acquired resistance, as well as jasmonic acid/ethylene spring loading in the offspring of infected plants. The importance of ethylene signaling in intergenerational resistance was confirmed by experiments on a knockdown line of ethylene insensitive 2 (ein2b), which lacks intergenerational acquired resistance. Taken together, these data indicate a role for DCL3a in regulating plant defense pathways during both within-generation and intergenerational resistance against nematodes in rice.

Non-coding RNAs-mediated environmental surveillance determines male fertility in plants

Plants are continuously exposed to environmental stresses leading to significant yield losses. With the changing climatic conditions, the intensity and duration of these stresses are expected to increase, posing a severe threat to crop productivity worldwide. Male gametogenesis is one of the most sensitive developmental stages. Exposure to environmental stresses during this stage leads to male sterility and yield loss. Elucidating the underlying molecular mechanism of environment-affected male sterility is essential to address this challenge. High-throughput RNA sequencing studies, loss-of-function phenotypes of sRNA biogenesis genes and functional genomics studies with non-coding RNAs have started to unveil the roles of small RNAs, long non-coding RNAs and the complex regulatory interactions between them in regulating male fertility under different growth regimes. Here, we discuss the current understanding of the non-coding RNA-mediated environmental stress surveillance and regulation of male fertility in plants. The candidate ncRNAs emerging from these studies can be leveraged to generate environment-sensitive male sterile lines for hybrid breeding or mitigate the impact of climate change on male fertility, as the situation demands.

Genome-wide identification of DCL, AGO, and RDR gene families in wheat (Triticum aestivum L.) and their expression analysis in response to heat stress

Key components of the RNA interference (RNAi) pathway include the Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. While these components have been studied in various plant species, their functional validation in wheat remains unexplored particularly under heat stress. In this study, a comprehensive genome-wide analysis to identify, and characterize DCL, AGO, and RDR genes in wheat and their expression patterns was carried out. Using phylogenetic analysis with orthologous genes from Arabidopsis and rice, we identified a total of 82 AGO, 31 DCL, and 31 RDR genes distributed across the 21 chromosomes of wheat. To understand the regulatory network, a network analysis of miRNAs that target RNA-silencing genes was performed. Our analysis revealed that 13 miRNAs target AGO genes, 8 miRNAs target DCL genes, and 10 miRNAs target RDR genes at different sites, respectively. Additionally, promoter analysis of the RNA-silencing genes was done and identified the presence of 132 cis-elements responsive to stress and phytohormones. To examine their expression patterns, we performed RNA-seq analysis in the flag leaf samples of wheat exposed to both normal and heat stress conditions. To understand the regulation of RNA silencing, we experimentally analysed the transcriptional changes in response to gradient heat stress treatments. Our results showed constitutive expression of the AGO1, AGO9, and DCL2 gene families, indicating their importance in the overall biological processes of wheat. Notably, RDR1, known to be involved in small interfering RNA (siRNA) biogenesis, exhibited higher expression levels in wheat leaf tissues. These findings suggest that these genes may play a role in responses to stress in wheat, highlighting their significance in adapting to environmental challenges. Overall, our study provides additional knowledge to understand the mechanisms underlying heat stress responses and emphasizes the essential roles of these gene families in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01362-0.

Noncoding RNAs and their roles in regulating the agronomic traits of crops

Molecular breeding is one of the most effective methods for improving the performance of crops. Understanding the genome features of crops, especially the physiological functions of individual genes, is of great importance to molecular breeding. Evidence has shown that genomes of both animals and plants transcribe numerous non-coding RNAs, which are involved in almost every aspect of development. In crops, an increasing number of studies have proven that non-coding RNAs are new genetic resources for regulating crop traits. In this review, we summarize the current knowledge of non-coding RNAs, which are potential crop trait regulators, and focus on the functions of long non-coding RNAs (lncRNAs) in determining crop grain yield, phased small-interfering RNAs (phasiRNAs) in regulating fertility, small interfering RNAs (siRNAs) and microRNAs (miRNAs) in facilitating plant immune response and disease resistance, and miRNAs mediating nutrient and metal stress. Finally, we also discuss the next-generation method for ncRNA application in crop domestication and breeding.

Transcriptional Profiling Reveals the Wheat Defences against Fusarium Head Blight Disease Regulated by a NAC Transcription Factor

The wheat NAC transcription factor TaNACL-D1 enhances resistance to the economically devastating Fusarium head blight (FHB) disease. The objective of this study was to decipher the alterations in gene expression, pathways and biological processes that led to enhanced resistance as a result of the constitutive expression of TaNACL-D1 in wheat. Transcriptomic analysis was used to determine the genes and processes enhanced in wheat due to TaNACL-D1 overexpression, both in the presence and absence of the causal agent of FHB, Fusarium graminearum (0- and 1-day post-treatment). The overexpression of TaNACL-D1 resulted in more pronounced transcriptional reprogramming as a response to fungal infection, leading to the enhanced expression of genes involved in detoxification, immune responses, secondary metabolism, hormone biosynthesis, and signalling. The regulation and response to JA and ABA were differentially regulated between the OE and the WT. Furthermore, the results suggest that the OE may more efficiently: (i) regulate the oxidative burst; (ii) modulate cell death; and (iii) induce both the phenylpropanoid pathway and lignin synthesis. Thus, this study provides insights into the mode of action and downstream target pathways for this novel NAC transcription factor, further validating its potential as a gene to enhance FHB resistance in wheat.

Biogenesis, Mode of Action and the Interactions of Plant Non-Coding RNAs

The central dogma of genetics, which outlines the flow of genetic information from DNA to RNA to protein, has long been the guiding principle in molecular biology. In fact, more than three-quarters of the RNAs produced by transcription of the plant genome are not translated into proteins, and these RNAs directly serve as non-coding RNAs in the regulation of plant life activities at the molecular level. The breakthroughs in high-throughput transcriptome sequencing technology and the establishment and improvement of non-coding RNA experiments have now led to the discovery and confirmation of the biogenesis, mechanisms, and synergistic effects of non-coding RNAs. These non-coding RNAs are now predicted to play important roles in the regulation of gene expression and responses to stress and evolution. In this review, we focus on the synthesis, and mechanisms of non-coding RNAs, and we discuss their impact on gene regulation in plants.

Spatial distribution of three ARGONAUTEs regulates the anther phasiRNA pathway

Argonaute protein (AGO) in association with small RNAs is the core machinery of RNA silencing, an essential mechanism for precise development and defense against pathogens in many organisms. Here, we identified two AGOs in rice anthers, AGO1b and AGO1d, that interact with phased small interfering RNAs (phasiRNAs) derived from numerous long non-coding RNAs. Moreover, 3D-immunoimaging and mutant analysis indicated that rice AGO1b and AGO1d cell type-specifically regulate anther development by acting as mobile carriers of these phasiRNAs from the somatic cell layers to the germ cells in anthers. Our study also highlights a new mode of reproductive RNA silencing via the specific nuclear and cytoplasmic localization of three AGOs, AGO1b, AGO1d, and MEL1, in rice pollen mother cells.

Phylogenetic analyses of seven protein families refine the evolution of small RNA pathways in green plants

Several protein families participate in the biogenesis and function of small RNAs (sRNAs) in plants. Those with primary roles include Dicer-like (DCL), RNA-dependent RNA polymerase (RDR), and Argonaute (AGO) proteins. Protein families such as double-stranded RNA-binding (DRB), SERRATE (SE), and SUPPRESSION OF SILENCING 3 (SGS3) act as partners of DCL or RDR proteins. Here, we present curated annotations and phylogenetic analyses of seven sRNA pathway protein families performed on 196 species in the Viridiplantae (aka green plants) lineage. Our results suggest that the RDR3 proteins emerged earlier than RDR1/2/6. RDR6 is found in filamentous green algae and all land plants, suggesting that the evolution of RDR6 proteins coincides with the evolution of phased small interfering RNAs (siRNAs). We traced the origin of the 24-nt reproductive phased siRNA-associated DCL5 protein back to the American sweet flag (Acorus americanus), the earliest diverged, extant monocot species. Our analyses of AGOs identified multiple duplication events of AGO genes that were lost, retained, or further duplicated in subgroups, indicating that the evolution of AGOs is complex in monocots. The results also refine the evolution of several clades of AGO proteins, such as AGO4, AGO6, AGO17, and AGO18. Analyses of nuclear localization signal sequences and catalytic triads of AGO proteins shed light on the regulatory roles of diverse AGOs. Collectively, this work generates a curated and evolutionarily coherent annotation for gene families involved in plant sRNA biogenesis/function and provides insights into the evolution of major sRNA pathways.

Small RNA-mediated DNA methylation during plant reproduction

Reproductive tissues are a rich source of small RNAs, including several classes of short interfering (si)RNAs that are restricted to this stage of development. In addition to RNA polymerase IV-dependent 24-nt siRNAs that trigger canonical RNA-directed DNA methylation, abundant reproductive-specific siRNAs are produced from companion cells adjacent to the developing germ line or zygote and may move intercellularly before inducing methylation. In some cases, these siRNAs are produced via non-canonical biosynthesis mechanisms or from sequences with little similarity to transposons. While the precise role of these siRNAs and the methylation they trigger is unclear, they have been implicated in specifying a single megaspore mother cell, silencing transposons in the male germ line, mediating parental dosage conflict to ensure proper endosperm development, hypermethylation of mature embryos, and trans-chromosomal methylation in hybrids. In this review, we summarize the current knowledge of reproductive siRNAs, including their biosynthesis, transport, and function.

SEMI-ROLLED LEAF 10 stabilizes catalase isozyme B to regulate leaf morphology and thermotolerance in rice (Oryza sativa L.)

Plant architecture and stress tolerance play important roles in rice breeding. Specific leaf morphologies and ideal plant architecture can effectively improve both abiotic stress resistance and rice grain yield. However, the mechanism by which plants simultaneously regulate leaf morphogenesis and stress resistance remains elusive. Here, we report that SRL10, which encodes a double-stranded RNA-binding protein, regulates leaf morphology and thermotolerance in rice through alteration of microRNA biogenesis. The srl10 mutant had a semi-rolled leaf phenotype and elevated sensitivity to high temperature. SRL10 directly interacted with catalase isozyme B (CATB), and the two proteins mutually increased one other's stability to enhance hydrogen peroxide (H₂ O₂ ) scavenging, thereby contributing to thermotolerance. The natural Hap3 (AGC) type of SRL10 allele was found to be present in the majority of aus rice accessions, and was identified as a thermotolerant allele under high temperature stress in both the field and the growth chamber. Moreover, the seed-setting rate was 3.19 times higher and grain yield per plant was 1.68 times higher in near-isogenic line (NIL) carrying Hap3 allele compared to plants carrying Hap1 allele under heat stress. Collectively, these results reveal a new locus of interest and define a novel SRL10-CATB based regulatory mechanism for developing cultivars with high temperature tolerance and stable yield. Furthermore, our findings provide a theoretical basis for simultaneous breeding for plant architecture and stress resistance.

CRISPR/Cas9-targeted mutagenesis of TaDCL4, TaDCL5 and TaRDR6 induces male sterility in common wheat

Phased, small interfering RNAs (phasiRNAs) are important for plant anther development, especially for male sterility. PhasiRNA biogenesis is dependent on genes like RNA polymerase 6 (RDR6), DICER-LIKE 4 (DCL4), or DCL5 to produce 21- or 24 nucleotide (nt) double-strand small RNAs. Here, we generated mutants of DCL4, DCL5 and RDR6 using CRISPR/Cas9 system and studied their effects on plant reproductive development and phasiRNA production in wheat. We found that RDR6 mutation caused sever consequence throughout plant development starting from seed germination and the dcl4 mutants grew weaker with thorough male sterility, while dcl5 plants developed normally but exhibited male sterility. Correspondingly, DCL4 and DCL5, respectively, specified 21- and 24-nt phasiRNA biogenesis, while RDR6 contributed to both. Also, the three key genes evolved differently in wheat, with TaDCL5-A/B becoming non-functioning and TaRDR6-A being lost after polyploidization. Furthermore, we found that PHAS genes (phasiRNA precursors) identified via phasiRNAs diverged rapidly among sub-genomes of polyploid wheat. Despite no similarity being found among phasiRNAs of grasses, their targets were enriched for similar biological functions. In light of the important roles of phasiRNA pathways in gametophyte development, genetic dissection of the function of key genes may help generate male sterile lines suitable for hybrid wheat breeding.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

An efficient deep learning based predictor for identifying miRNA-triggered phasiRNA loci in plant

Phasic small interfering RNAs are plant secondary small interference RNAs that typically generated by the convergence of miRNAs and polyadenylated mRNAs. A growing number of studies have shown that miRNA-initiated phasiRNA plays crucial roles in regulating plant growth and stress responses. Experimental verification of miRNA-initiated phasiRNA loci may take considerable time, energy and labor. Therefore, computational methods capable of processing high throughput data have been proposed one by one. In this work, we proposed a predictor (DIGITAL) for identifying miRNA-initiated phasiRNAs in plant, which combined a multi-scale residual network with a bi-directional long-short term memory network. The negative dataset was constructed based on positive data, through replacing 60% of nucleotides randomly in each positive sample. Our predictor achieved the accuracy of 98.48% and 94.02% respectively on two independent test datasets with different sequence length. These independent testing results indicate the effectiveness of our model. Furthermore, DIGITAL is of robustness and generalization ability, and thus can be easily extended and applied for miRNA target recognition of other species. We provide the source code of DIGITAL, which is freely available at https://github.com/yuanyuanbu/DIGITAL.

Mobile ARGONAUTE 1d binds 22-nt miRNAs to generate phasiRNAs important for low-temperature male fertility in rice

Phased small interfering RNAs (phasiRNAs) are abundantly expressed in anthers and linked to environment-related male fertility in grasses, yet how they function under different environmental conditions remains unclear. Here, we identified a rice (Oryza sativa) low temperature-induced Argonaute (AGO) protein, OsAGO1d, that is responsible for generating phasiRNAs and preserving male fertility at low temperature. Loss of OsAGO1d function causes low-temperature male sterility associated with delayed programmed cell death of tapetal cells during anther development. OsAGO1d binds miR2118 and miR2275 family members and triggers phasiRNA biogenesis; it also binds 21-nt phasiRNAs with a 5' terminal U. In total, phasiRNAs from 972 loci are OsAGO1d-dependent. OsAGO1d protein moves from anther wall cells into meiocytes, where it loads miR2275 to produce 24-nt phasiRNAs. Together, our results show that OsAGO1d acts as a mobile signal to fine-tune phasiRNA production and this function is important for male fertility at low temperature.

The coordinated regulation mechanism of rice plant architecture and its tolerance to stress

Rice plant architecture and stress tolerance have historically been primary concerns for rice breeders. The "Green Revolution" and super-rice breeding practices have demonstrated that ideal plant architecture can effectively improve both stress tolerance and yield. The synergistic selection and breeding of rice varieties with ideal architecture and stress tolerance can increase and stabilize yield. While rice plant plant architecture and stress tolerance are separately regulated by complicated genetic networks, the molecular mechanisms underlying their relationships and synergism have not yet been explored. In this paper, we review the regulatory mechanism between plant architecture, stress tolerance, and biological defense at the different level to provide a theoretical basis for the genetic network of the synergistic regulation and improvement of multiple traits.

DNA methylation dynamics during germline development

DNA methylation plays essential homeostatic functions in eukaryotic genomes. In animals, DNA methylation is also developmentally regulated and, in turn, regulates development. In the past two decades, huge research effort has endorsed the understanding that DNA methylation plays a similar role in plant development, especially during sexual reproduction. The power of whole-genome sequencing and cell isolation techniques, as well as bioinformatics tools, have enabled recent studies to reveal dynamic changes in DNA methylation during germline development. Furthermore, the combination of these technological advances with genetics, developmental biology and cell biology tools has revealed functional methylation reprogramming events that control gene and transposon activities in flowering plant germlines. In this review, we discuss the major advances in our knowledge of DNA methylation dynamics during male and female germline development in flowering plants.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.

Degradome Assisted Plant MicroRNA Prediction Under Alternative Annotation Criteria

Current microRNA (miRNA) prediction methods are generally based on annotation criteria that tend to miss potential functional miRNAs. Recently, new miRNA annotation criteria have been proposed that could lead to improvements in miRNA prediction methods in plants. Here, we investigate the effect of the new criteria on miRNA prediction in Arabidopsis thaliana and present a new degradome assisted functional miRNA prediction approach. We investigated the effect by applying the new criteria, and a more permissive criteria on miRNA prediction using existing miRNA prediction tools. We also developed an approach to miRNA prediction that is assisted by the functional information extracted from the analysis of degradome sequencing. We demonstrate the improved performance of degradome assisted miRNA prediction compared to unassisted prediction and evaluate the approach using miRNA differential expression analysis. We observe how the miRNA predictions fit under the different criteria and show a potential novel miRNA that has been missed within Arabidopsis thaliana. Additionally, we introduce a freely available software 'PAREfirst' that employs the degradome assisted approach. The study shows that some miRNAs could be missed due to the stringency of the former annotation criteria, and combining a degradome assisted approach with more permissive miRNA criteria can expand confident miRNA predictions.

Epigenetic Dynamics and Regulation of Plant Male Reproduction

Flowering plant male germlines develop within anthers and undergo epigenetic reprogramming with dynamic changes in DNA methylation, chromatin modifications, and small RNAs. Profiling the epigenetic status using different technologies has substantially accumulated information on specific types of cells at different stages of male reproduction. Many epigenetically related genes involved in plant gametophyte development have been identified, and the mutation of these genes often leads to male sterility. Here, we review the recent progress on dynamic epigenetic changes during pollen mother cell differentiation, microsporogenesis, microgametogenesis, and tapetal cell development. The reported epigenetic variations between male fertile and sterile lines are summarized. We also summarize the epigenetic regulation-associated male sterility genes and discuss how epigenetic mechanisms in plant male reproduction can be further revealed.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

3D Imaging and In Situ Hybridization for Uncovering the Functions of MicroRNA in Rice Anther

Small RNAs specifically expressed in reproductive tissues are key regulators of germline development in eukaryotes. Rice microRNA2118 (miR2118), which is enriched during reproduction in grasses, is a trigger to produce phased small interfering RNAs (phasiRNAs). These phasiRNAs demonstrate the temporal regulation with premeiotic phasiRNAs and meiotic phasiRNAs in anther development. Furthermore, the site-specific regulation via miR2118 and phasiRNAs is of importance in soma and germ development in anthers. Accordingly, histological imaging methods are essential tools for understanding spatiotemporal regulation during reproduction and elucidating the reproductive roles of miRNAs and phasiRNAs. We successfully developed a method to visualize the three-dimensional (3D) structure of entire rice anthers, which can also be used for distinguishing the internal structure of the anthers in other plants. Here, we describe the detailed methods of in situ hybridization for miR2118 localization and the visualization of the 3D structure of entire anthers of rice.

Integrative Analysis of the DICER-like (DCL) Genes From Peach (Prunus persica): A Critical Role in Response to Drought Stress

DICER-likes (DCLs) proteins are the core component for non-coding RNA (ncRNA) biogenesis, playing essential roles in some biological processes. The DCL family has been characterized in model plants, such as Arabidopsis, rice, and poplar. However, the evolutionary aspect and the expression mechanism under drought stress were scarce and have never been reported and characterized in one of the most important worldwide cultivated fruit trees, peach (Prunus persica). Eight DCLs genes in the Prunus persica genome were detected, in addition to 51 DCLs in the other seven Rosaceae genomes. The phylogenetic analysis with Arabidopsis thaliana and RTL1 gene as outgroups suggested that DCL members are divided into four clades: DCL1, DCL2, DCL3, and DCL4 with several gene gain/loss events of DCL gene copies through the evolutionary tract of the Rosacea family. The number of homologous DCL copies within each clade, along with the chromosomal location indicated gene duplication event of the DCL2 gene occurred once for the subfamily Amygdaloideae and twice for Pyrus communis and Prunus dulics and trice for the P. persica on Chromosome number 7 genes. Another duplication event was found for the DCL3 gene that occurred once for all the eight Rosaceae species with no match in A. thaliana. The DCL genetic similarity and activity was evaluated using BLASTp and previously published RNA-seq data among different tissues and over different time points of peach trees exposed to drought conditions. Finally, the expression pattern of PrupeDCLs in response to drought stress was identified, and two of these members, Prupe.7G047900 and Prupe.6G363600, were found as main candidate genes for response to drought stress. Our data presented here provide useful information for a better understanding of the molecular evolution of DCL genes in Rosaceae genomes, and the function of DCLs in P. persica.

24-nt phasiRNAs move from tapetal to meiotic cells in maize anthers

In maize, 24-nt phased, secondary small interfering RNAs (phasiRNAs) are abundant in meiotic stage anthers, but their distribution and functions are not precisely known. Using laser capture microdissection, we analyzed tapetal cells, meiocytes and other somatic cells at several stages of anther development to establish the timing of 24-PHAS precursor transcripts and the 24-nt phasiRNA products. By integrating RNA and small RNA profiling plus single-molecule and small RNA FISH (smFISH or sRNA-FISH) spatial detection, we demonstrate that the tapetum is the primary site of 24-PHAS precursor and Dcl5 transcripts and the resulting 24-nt phasiRNAs. Interestingly, 24-nt phasiRNAs accumulate in all cell types, with the highest levels in meiocytes, followed by tapetum. Our data support the conclusion that 24-nt phasiRNAs are mobile from tapetum to meiocytes and to other somatic cells. We discuss possible roles for 24-nt phasiRNAs in anther cell types.

Control of OsARF3a by OsKANADI1 contributes to lemma development in rice

In rice (Oryza sativa), the lemma and palea protect the internal organs of the floret，provide nutrients for seed development, and determine grain size. We previously revealed that a trans-acting small interfering RNA targeting AUXIN RESPONSE FACTORS (tasiR-ARF) regulates lemma polarity establishment via post-transcriptional repression of AUXIN RESPONSE FACTORS (ARFs) in rice. TasiR-ARF formation requires RNA-DEPENDENT RNA POLYMERASE 6 (RDR6). However, the underlying molecular mechanism of the tasiR-ARF-ARF regulon in lemma development remains unclear. Here, by genetic screening for suppressors of the thermosensitive mutant osrdr6-1, we identified three suppressors, huifu 1 (hf1), hf9, and hf17. Mapping-by-sequencing revealed that HF1 encodes a MYB transcription factor belonging to the KANADI1 family. The hf1 mutation partially rescued the osrdr6-1 lemma defect but not the defect in tasiR-ARF levels. DNA affinity purification sequencing analysis identified 17 725 OsKANADI1-associated sites, most of which contain the SPBP-box binding motif (RGAATAWW) and are located in the promoter, protein-coding, intron, and intergenic regions. Moreover, we found that OsKANADI1 could directly bind to the intron of OsARF3a in vitro and in vivo and promote OsARF3a expression at the transcriptional level. In addition, hf9 and hf17 are intragenic suppressors containing mutations in OsRDR6 that partially rescue tasiR-ARF levels by restoring OsRDR6 protein levels. Collectively, our results demonstrate that OsKANADI1 and tasiR-ARFs synergistically maintain the proper expression of OsARF3a and thus contribute to rice lemma development.

Small regulatory RNAs in rice epigenetic regulation

Plant small RNAs (sRNAs) are short non-coding RNAs that are implicated in various regulatory processes involving post-transcriptional gene silencing and epigenetic gene regulation. In epigenetic regulation, sRNAs are primarily involved in RNA-directed DNA methylation (RdDM) pathways. sRNAs in the RdDM pathways play a role not only in the suppression of transposable element (TE) activity but also in gene expression regulation. Although the major components of the RdDM pathways have been well studied in Arabidopsis, recent studies have revealed that the RdDM pathways in rice have important biological functions in stress response and developmental processes. In this review, we summarize and discuss recent literature on sRNA-mediated epigenetic regulation in rice. First, we describe the RdDM mechanisms in plants. We then introduce recent discoveries on the biological roles of rice genes involved in the RdDM pathway and TE-derived sRNAs working at specific genomic loci for epigenetic control in rice.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

Silencing of Dicer-like protein 2a restores the resistance phenotype in the rice mutant, sxi4 (suppressor of Xa21-mediated immunity 4)

The rice immune receptor XA21 confers resistance to Xanthomonas oryzae pv. oryzae (Xoo), and upon recognition of the RaxX21-sY peptide produced by Xoo, XA21 activates the plant immune response. Here we screened 21 000 mutant plants expressing XA21 to identify components involved in this response, and reported here the identification of a rice mutant, sxi4, which is susceptible to Xoo. The sxi4 mutant carries a 32-kb translocation from chromosome 3 onto chromosome 7 and displays an elevated level of DCL2a transcript, encoding a Dicer-like protein. Silencing of DCL2a in the sxi4 genetic background restores resistance to Xoo. RaxX21-sY peptide-treated leaves of sxi4 retain the hallmarks of XA21-mediated immune response. However, WRKY45-1, a known negative regulator of rice resistance to Xoo, is induced in the sxi4 mutant in response to RaxX21-sY peptide treatment. A CRISPR knockout of a short interfering RNA (TE-siRNA815) in the intron of WRKY45-1 restores the resistance phenotype in sxi4. These results suggest a model where DCL2a accumulation negatively regulates XA21-mediated immunity by altering the processing of TE-siRNA815.

Functional specialization of monocot DCL3 and DCL5 proteins through the evolution of the PAZ domain

Monocot DICER-LIKE3 (DCL3) and DCL5 produce distinct 24-nt small interfering RNAs (siRNAs), heterochromatic siRNAs (hc-siRNAs) and phased secondary siRNAs (phasiRNAs), respectively. The former small RNAs are linked to silencing of transposable elements and heterochromatic repeats, and the latter to reproductive processes. It is assumed that these DCLs evolved from an ancient ‘eudicot-type’ DCL3 ancestor, which may have produced both types of siRNAs. However, how functional differentiation was achieved after gene duplication remains elusive. Here, we find that monocot DCL3 and DCL5 exhibit biochemically distinct preferences for 5′ phosphates and 3′ overhangs, consistent with the structural properties of their in vivo double-stranded RNA substrates. Importantly, these distinct substrate specificities are determined by the PAZ domains of DCL3 and DCL5, which have accumulated mutations during the course of evolution. These data explain the mechanism by which these DCLs cleave their cognate substrates from a fixed end, ensuring the production of functional siRNAs. Our study also indicates how plants have diversified and optimized RNA silencing mechanisms during evolution.

A cascade of bHLH-regulated pathways programs maize anther development

The spatiotemporal development of somatic tissues of the anther lobe is necessary for successful fertile pollen production. This process is mediated by many transcription factors acting through complex, multi-layered networks. Here, our analysis of functional knockout mutants of interacting basic helix-loop-helix genes Ms23, Ms32, basic helix-loop-helix 122 (bHLH122), and bHLH51 in maize (Zea mays) established that male fertility requires all four genes, expressed sequentially in the tapetum (TP). Not only do they regulate each other, but also they encode proteins that form heterodimers that act collaboratively to guide many cellular processes at specific developmental stages. MS23 is confirmed to be the master factor, as the ms23 mutant showed the earliest developmental defect, cytologically visible in the TP, with the most drastic alterations in premeiotic gene expression observed in ms23 anthers. Notably, the male-sterile ms23, ms32, and bhlh122-1 mutants lack 24-nt phased secondary small interfering RNAs (phasiRNAs) and the precursor transcripts from the corresponding 24-PHAS loci, while the bhlh51-1 mutant has wild-type levels of both precursors and small RNA products. Multiple lines of evidence suggest that 24-nt phasiRNA biogenesis primarily occurs downstream of MS23 and MS32, both of which directly activate Dcl5 and are required for most 24-PHAS transcription, with bHLH122 playing a distinct role in 24-PHAS transcription.

MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development

MicroRNAs (miRNAs) are a class of 21-24 nucleotides (nt) noncoding small RNAs ubiquitously distributed across the plant kingdom. miR482/2118, one of the conserved miRNA superfamilies originating from gymnosperms, has divergent main functions in core-angiosperms. It mainly regulates NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT (NBS-LRR) genes in eudicots, functioning as an essential component in plant disease resistance; in contrast, it predominantly targets numerous long noncoding RNAs (lncRNAs) in monocot grasses, which are vital for plant reproduction. Usually, miR482/2118 is 22-nt in length, which can trigger the production of phased small interfering RNAs (phasiRNAs) after directed cleavage. PhasiRNAs instigated from target genes of miR482/2118 enhance their roles in corresponding biological processes by cis-regulation on cognate genes and expands their function to other pathways via trans activity on different genes. This review summarizes the origin, biogenesis, conservation, and evolutionary characteristics of the miR482/2118 superfamily and delineates its diverse functions in disease resistance, plant development, stress responses, etc.

Mechanism for the genomic and functional evolution of the MIR2118 family in the grass lineage

The MIR2118 family has undergone tremendous expansion in the grass lineage, in which the miRNA targets numerous noncoding PHAS loci to produce 21-nt phased small interfering RNAs (phasiRNAs) involved in male fertility. However, the evolutionary trajectory of the grass MIR2118 genes and the functions of phasiRNAs have not yet been fully elucidated. We conducted comparative genomic, molecular evolution, expression and parallel analysis of RNA ends (PARE) analyses of MIR2118 and the miR2118-mediated regulatory pathway in grasses, focusing on Oryza sativa. In total, 617 MIR2118 and eight MIR1859 novel members were identified. Phylogenetic analyses showed that grass MIR2118 genes form a distinct clade from the MIR482/2118 genes of nongrass species. We reconstructed hypothetical evolutionary histories of the grass MIR2118 clusters and its MIR1859 variants, and examined the polycistronic composition and the differential expression of the osa-MIR2118 clusters. PARE data showed that osa-miR2118 might also direct the cleavage of some protein-coding gene transcripts. Importantly, we found that PARE analysis is inherently prone to false-positive target predictions when a large number of small RNAs, such as phasiRNAs, are analysed. Our results revealed the evolution and diversification of the MIR2118 family, and provide new insights into the functions of phasiRNAs in the grasses.

RNA silencing: From discovery and elucidation to application and perspectives

RNA silencing (or RNA interference, RNAi) is a conserved mechanism for regulating gene expression in eukaryotes. The discovery of natural trans-kingdom RNAi indicated that small RNAs act as signaling molecules and enable communication between organisms in different kingdoms. The phenomenon and potential mechanisms of trans-kingdom RNAi are among the most exciting research topics. To better understand trans-kingdom RNAi, we review the history of the discovery and elucidation of RNAi mechanisms. Based on canonical RNAi mechanisms, we summarize the major points of divergence around RNAi pathways in the main eukaryotes' kingdoms, including plants, animals, and fungi. We review the representative incidents associated with the mechanisms and applications of trans-kingdom RNAi in crop protection, and discuss the critical factors that should be considered to develop successful trans-kingdom RNAi-based crop protection strategies.

The rice EP3 and OsFBK1 E3 ligases alter plant architecture and flower development, and affect transcript accumulation of microRNA pathway genes and their targets

ERECTA PANICLE 3 (EP3) and ORYZA SATIVA F-BOX KELCH 1 (OsFBK1) proteins share 57% and 54% sequence identity with the Arabidopsis F-box protein HAWAIIAN SKIRT (HWS). Previously we showed that EP3 is a functional orthologue of HWS. Here we demonstrate that OsFBK1 is another functional orthologue of HWS and show the complexity of interaction between EP3 and OsFBK1 genes at different developmental stages of the plant. qRT-PCR expression analyses and studies of EP3-GFP and OsFBK1-RFP promoter reporter lines demonstrate that although EP3 and OsFBK1 expression can be detected in the same tissues some cells exclusively express EP3 or OsFBK1 whilst others co-express both genes. Loss, reduction or gain-of-function lines for EP3 and OsFBK1, show that EP3 and OsFBK1 affect plant architecture, organ size, floral organ number and size, floral morphology, pollen viability, grain size and weight. We have identified the putative orthologue genes of the rice microRNA pathway for ORYZA SATIVA DAWDLE (OsDDL) and ORYZA SATIVA SERRATE (OsSE), and demonstrated that EP3 and OsFBK1 affect their transcript levels as well as those of CROWN ROOT DEFECT 1/ORYZA SATIVA Exportin-5 HASTY (CRD1/OsHST), ORYZA SATIVA DICER-LIKE 1 (OsDCL) and ORYZA SATIVA WEAVY LEAF1 (OsWAF1). We show that EP3 affects OsPri-MIR164, OsNAM1 and OsNAC1 transcript levels. OsNAC1 transcripts are modified by OsFBK1, suggesting two independent regulatory pathways, one via EP3 and OsMIR164 and the other via OsFBK1. Our data propose that EP3 and OsFBK1 conjointly play similar roles in rice to how HWS does in Arabidopsis.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

Coordinated Action of RTBV and RTSV Proteins Suppress Host RNA Silencing Machinery

RNA silencing is as an adaptive immune response in plants that limits the accumulation or spread of invading viruses. Successful virus infection entails countering the RNA silencing machinery for efficient replication and systemic spread in the host. The viruses encode proteins with the ability to suppress or block the host silencing mechanism, resulting in severe pathogenic symptoms and diseases. Tungro is a viral disease caused by a complex of two viruses and it provides an excellent system to understand the host and virus interactions during infection. It is known that Rice tungro bacilliform virus (RTBV) is the major determinant of the disease while Rice tungro spherical virus (RTSV) accentuates the symptoms. This study brings to focus the important role of RTBV ORF-IV in disease manifestation, by acting as both the victim and silencer of the RNA silencing pathway. The ORF-IV is a weak suppressor of the S-PTGS or stable silencing, but its suppression activity is augmented in the presence of specific RTSV proteins. Among these, RTBV ORF-IV and RTSV CP3 proteins interact with each other. This interaction may lead to the suppression of localized silencing as well as the spread of silencing in the host plants. The findings present a probable mechanistic glimpse of the requirement of the two viruses in enhancing tungro disease.

Spatiotemporal regulation and roles of reproductive phasiRNAs in plants

Since co-suppression was discovered as a pioneer silencing phenomenon of RNA interference (RNAi) in petunia in 1990, many types of small RNAs have been identified in the RNAi pathway among various eukaryotes. In plants, a large number of 21- or 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs) are produced via processing of long RNA precursors by Dicer-like proteins. However, the roles of phasiRNAs remain largely unknown. The development of imaging technology and RNA profiling has clarified the spatiotemporal regulation of phasiRNAs, and subsequently the different functions of 21-nt trans-acting phasiRNAs and 24-nt cis-regulatory phasiRNAs during male organ development. This review focuses on the biogenesis, diversification, spatiotemporal expression pattern and function of phasiRNAs in plants.