A One Precursor One siRNA Model for Pol IV-Dependent siRNA Biogenesis.

RNA-directed DNA methylation in Arabidopsis thaliana is driven by the plant-specific RNA Polymerase IV (Pol IV). It has been assumed that a Pol IV transcript can give rise to multiple 24-nt small interfering RNAs (siRNAs) that target DNA methylation. Here, we demonstrate that Pol IV-dependent RNAs (P4RNAs) from wild-type Arabidopsis are surprisingly short in length (30 to 40 nt) and mirror 24-nt siRNAs in distribution, abundance, strand bias, and 5'-adenine preference. P4RNAs exhibit transcription start sites similar to Pol II products and are featured with 5'-monophosphates and 3'-misincorporated nucleotides. The 3'-misincorporation preferentially occurs at methylated cytosines on the template DNA strand, suggesting a co-transcriptional feedback to siRNA biogenesis by DNA methylation to reinforce silencing locally. These results highlight an unusual mechanism of Pol IV transcription and suggest a "one precursor, one siRNA" model for the biogenesis of 24-nt siRNAs in Arabidopsis.

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.

Plastid dsRNA transgenes trigger phased small RNA-based gene silencing of nuclear-encoded genes.

Plastid transformation technology has been widely used to express traits of potential commercial importance, though the technology has been limited to traits that function while sequestered in the organelle. Prior research indicates that plastid contents can escape from the organelle, suggesting a possible mechanism for engineering plastid transgenes to function in other cellular locations. To test this hypothesis, we created tobacco (Nicotiana tabacum cv. Petit Havana) plastid transformants that express a fragment of the nuclear-encoded Phytoene desaturase (PDS) gene capable of catalyzing post-transcriptional gene silencing if RNA escapes into the cytoplasm. We found multiple lines of direct evidence that plastid-encoded PDS transgenes affect nuclear PDS gene silencing: knockdown of the nuclear-encoded PDS mRNA and/or its apparent translational inhibition, biogenesis of 21-nucleotide (nt) phased small interfering RNAs (phasiRNAs), and pigment-deficient plants. Furthermore, plastid-expressed dsRNA with no cognate nuclear-encoded pairing partner also produced abundant 21-nt phasiRNAs in the cytoplasm, demonstrating that a nuclear-encoded template is not required for siRNA biogenesis. Our results indicate that RNA escape from plastids to the cytoplasm occurs generally, with functional consequences that include entry into the gene silencing pathway. Furthermore, we uncover a method to produce plastid-encoded traits with functions outside of the organelle and open additional fields of study in plastid development, compartmentalization, and small RNA biogenesis.

The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming.

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

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.

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.

Arabidopsis apoplastic fluid contains sRNA- and circular RNA-protein complexes that are located outside extracellular vesicles.

Previously, we have shown that apoplastic wash fluid (AWF) purified from Arabidopsis leaves contains small RNAs (sRNAs). To investigate whether these sRNAs are encapsulated inside extracellular vesicles (EVs), we treated EVs isolated from Arabidopsis leaves with the protease trypsin and RNase A, which should degrade RNAs located outside EVs but not those located inside. These analyses revealed that apoplastic RNAs are mostly located outside and are associated with proteins. Further analyses of these extracellular RNAs (exRNAs) revealed that they include both sRNAs and long noncoding RNAs (lncRNAs), including circular RNAs (circRNAs). We also found that exRNAs are highly enriched in the posttranscriptional modification N6-methyladenine (m6A). Consistent with this, we identified a putative m6A-binding protein in AWF, GLYCINE-RICH RNA-BINDING PROTEIN 7 (GRP7), as well as the sRNA-binding protein ARGONAUTE2 (AGO2). These two proteins coimmunoprecipitated with lncRNAs, including circRNAs. Mutation of GRP7 or AGO2 caused changes in both the sRNA and lncRNA content of AWF, suggesting that these proteins contribute to the secretion and/or stabilization of exRNAs. We propose that exRNAs located outside of EVs mediate host-induced gene silencing, rather than RNA located inside EVs.

Quality control and evaluation of plant epigenomics data.

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

Coexpression network and trans-activation analyses of maize reproductive phasiRNA loci.

The anther-enriched phased, small interfering RNAs (phasiRNAs) play vital roles in sustaining male fertility in grass species. Their long non-coding precursors are synthesized by RNA polymerase II and are likely regulated by transcription factors (TFs). A few putative transcriptional regulators of the 21- or 24-nucleotide phasiRNA loci (referred to as 21- or 24-PHAS loci) have been identified in maize (Zea mays), but whether any of the individual TFs or TF combinations suffice to activate any PHAS locus is unclear. Here, we identified the temporal gene coexpression networks (modules) associated with maize anther development, including two modules highly enriched for the 21- or 24-PHAS loci. Comparisons of these coexpression modules and gene sets dysregulated in several reported male sterile TF mutants provided insights into TF timing with regard to phasiRNA biogenesis, including antagonistic roles for OUTER CELL LAYER4 and MALE STERILE23. Trans-activation assays in maize protoplasts of individual TFs using bulk-protoplast RNA-sequencing showed that two of the TFs coexpressed with 21-PHAS loci could activate several 21-nucleotide phasiRNA pathway genes but not transcription of 21-PHAS loci. Screens for combinatorial activities of these TFs and, separately, the recently reported putative transcriptional regulators of 24-PHAS loci using single-cell (protoplast) RNA-sequencing, did not detect reproducible activation of either 21-PHAS or 24-PHAS loci. Collectively, our results suggest that the endogenous transcriptional machineries and/or chromatin states in the anthers are necessary to activate reproductive PHAS loci.

Temperature and light reverse the fertility of rice P/TGMS line ostms19 via reactive oxygen species homeostasis.

P/TGMS (Photo/thermo-sensitive genic male sterile) lines are crucial resources for two-line hybrid rice breeding. Previous studies revealed that slow development is a general mechanism for sterility-fertility conversion of P/TGMS in Arabidopsis. However, the difference in P/TGMS genes between rice and Arabidopsis suggests the presence of a distinct P/TGMS mechanism in rice. In this study, we isolated a novel P/TGMS line, ostms19, which shows sterility under high-temperature conditions and fertility under low-temperature conditions. OsTMS19 encodes a novel pentatricopeptide repeat (PPR) protein essential for pollen formation, in which a point mutation GTA(Val) to GCA(Ala) leads to ostms19 P/TGMS phenotype. It is highly expressed in the tapetum and localized to mitochondria. Under high temperature or long-day photoperiod conditions, excessive ROS accumulation in ostms19 anthers during pollen mitosis disrupts gene expression and intine formation, causing male sterility. Conversely, under low temperature or short-day photoperiod conditions, ROS can be effectively scavenged in anthers, resulting in fertility restoration. This indicates that ROS homeostasis is critical for fertility conversion. This relationship between ROS homeostasis and fertility conversion has also been observed in other tested rice P/TGMS lines. Therefore, we propose that ROS homeostasis is a general mechanism for the sterility-fertility conversion of rice P/TGMS lines.

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.

miRador: a fast and precise tool for the prediction of plant miRNAs.

Plant microRNAs (miRNAs) are short, noncoding RNA molecules that restrict gene expression via posttranscriptional regulation and function in several essential pathways, including development, growth, and stress responses. Accurately identifying miRNAs in populations of small RNA sequencing libraries is a computationally intensive process that has resulted in the misidentification of inaccurately annotated miRNA sequences. In recent years, criteria for miRNA annotation have been refined with the aim to reduce these misannotations. Here, we describe miRador, a miRNA identification tool that utilizes the most up-to-date, community-established criteria for accurate identification of miRNAs in plants. We combined target prediction and Parallel Analysis of RNA Ends (PARE) data to assess the precision of the miRNAs identified by miRador. We compared miRador to other commonly used miRNA prediction tools and found that miRador is at least as precise as other prediction tools while being substantially faster than other tools. miRador should be broadly useful for the plant community to identify and annotate miRNAs in plant genomes.

Correction: ARGONAUTE10 promotes the degradation of miR165/6 through the SDN1 and SDN2 exonucleases in Arabidopsis.

[This corrects the article DOI: 10.1371/journal.pbio.2001272.].

Introduction of barnase/barstar in soybean produces a rescuable male sterility system for hybrid breeding.

Hybrid breeding for increased vigour has been used for over a century to boost agricultural outputs without requiring higher inputs. While this approach has led to some of the most substantial gains in crop productivity, breeding barriers have fundamentally limited soybean (Glycine max) from reaping the benefits of hybrid vigour. Soybean flowers self-pollinate prior to opening and thus are not readily amenable to outcrossing. In this study, we demonstrate that the barnase/barstar male sterility/rescue system can be used in soybean to produce hybrid seeds. By expressing the cytotoxic ribonuclease, barnase, under a tapetum-specific promoter in soybean anthers, we are able to completely block pollen maturation, creating male sterile plants. We show that fertility can be rescued in the F1 generation of these barnase-expressing lines when they are crossed with pollen from plants that express the barnase inhibitor, barstar. Importantly, we found that the successful rescue of male fertility is dependent on the relative dosage of barnase and barstar. When barnase and barstar were expressed under the same tapetum-specific promoter, the F1 offspring remained male sterile. When we expressed barstar under a relatively stronger promoter than barnase, we were able to achieve a successful rescue of male fertility in the F1 generation. This work demonstrates the successful implementation of a biotechnology approach to produce fertile hybrid offspring in soybean.

Organizing your space: The potential for integrating spatial transcriptomics and 3D imaging data in plants.

Plant cells communicate information for the regulation of development and responses to external stresses. A key form of this communication is transcriptional regulation, accomplished via complex gene networks operating both locally and systemically. To fully understand how genes are regulated across plant tissues and organs, high resolution, multi-dimensional spatial transcriptional data must be acquired and placed within a cellular and organismal context. Spatial transcriptomics (ST) typically provides a two-dimensional spatial analysis of gene expression of tissue sections that can be stacked to render three-dimensional data. For example, X-ray and light-sheet microscopy provide sub-micron scale volumetric imaging of cellular morphology of tissues, organs, or potentially entire organisms. Linking these technologies could substantially advance transcriptomics in plant biology and other fields. Here, we review advances in ST and 3D microscopy approaches and describe how these technologies could be combined to provide high resolution, spatially organized plant tissue transcript mapping.

The evolutionary history of small RNAs in Solanaceae.

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

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

Phased secondary small interfering RNAs (phasiRNAs) constitute a major category of small RNAs in plants, but most of their functions are still poorly defined. Some phasiRNAs, known as trans-acting siRNAs, are known to target complementary mRNAs for degradation and to function in development. However, the targets or biological roles of other phasiRNAs remain speculative. New insights into phasiRNA biogenesis, their conservation, and their variation across the flowering plants continue to emerge due to the increased availability of plant genomic sequences, deeper and more sophisticated sequencing approaches, and improvements in computational biology and biochemical/molecular/genetic analyses. In this review, we survey recent progress in phasiRNA biology, with a particular focus on two classes associated with male reproduction: 21-nucleotide (accumulate early in anther ontogeny) and 24-nucloetide (produced in somatic cells during meiosis) phasiRNAs. We describe phasiRNA biogenesis, function, and evolution and define the unanswered questions that represent topics for future research.

Sex Determination by Two Y-Linked Genes in Garden Asparagus.

The origin and early evolution of sex chromosomes have been hypothesized to involve the linkage of factors with antagonistic effects on male and female function. Garden asparagus (Asparagus officinalis) is an ideal species to investigate this hypothesis, as the X and Y chromosomes are cytologically homomorphic and evolved from an ancestral autosome pair in association with a shift from hermaphroditism to dioecy. Mutagenesis screens paired with single-molecule fluorescence in situ hybridization directly implicate Y-specific genes that respectively suppress female (pistil) development and are necessary for male (anther) development. Comparison of contiguous X and Y chromosome assemblies shows that hemizygosity underlies the loss of recombination between the genes suppressing female organogenesis (SUPPRESSOR OF FEMALE FUNCTION) and promoting male function (TAPETAL DEVELOPMENT AND FUNCTION1 [aspTDF1]). We also experimentally demonstrate the function of aspTDF1. These findings provide direct evidence that sex chromosomes can function through linkage of two sex determination genes.

Developmental Genetics of Corolla Tube Formation: Role of the tasiRNA-ARF Pathway and a Conceptual Model.

Over 80,000 angiosperm species produce flowers with petals fused into a corolla tube. The corolla tube contributes to the tremendous diversity of flower morphology and plays a critical role in plant reproduction, yet it remains one of the least understood plant structures from a developmental genetics perspective. Through mutant analyses and transgenic experiments, we show that the tasiRNA-ARF pathway is required for corolla tube formation in the monkeyflower species Mimulus lewisii Loss-of-function mutations in the M. lewisii orthologs of ARGONAUTE7 and SUPPRESSOR OF GENE SILENCING3 cause a dramatic decrease in abundance of TAS3-derived small RNAs and a moderate upregulation of AUXIN RESPONSE FACTOR3 (ARF3) and ARF4, which lead to inhibition of lateral expansion of the bases of petal primordia and complete arrest of the upward growth of the interprimordial regions, resulting in unfused corollas. Using the DR5 auxin-responsive promoter, we discovered that auxin signaling is continuous along the petal primordium base and the interprimordial region during the critical stage of corolla tube formation in the wild type, similar to the spatial pattern of MlARF4 expression. Auxin response is much weaker and more restricted in the mutant. Furthermore, exogenous application of a polar auxin transport inhibitor to wild-type floral apices disrupted petal fusion. Together, these results suggest a new conceptual model highlighting the central role of auxin-directed synchronized growth of the petal primordium base and the interprimordial region in corolla tube formation.

Soybean DICER-LIKE2 Regulates Seed Coat Color via Production of Primary 22-Nucleotide Small Interfering RNAs from Long Inverted Repeats.

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