Trans-species microRNA loci in the parasitic plant Cuscuta campestris have a U6-like snRNA promoter.

Small regulatory RNAs can move between organisms and regulate gene expression in the recipient. Whether the trans-species small RNAs being exported are distinguished from the normal endogenous small RNAs of the source organism is not known. The parasitic plant Cuscuta campestris (dodder) produces many microRNAs that specifically accumulate at the host-parasite interface, several of which have trans-species activity. We found that induction of C. campestris interface-induced microRNAs is similar regardless of host species and occurs in C. campestris haustoria produced in the absence of any host. The loci-encoding C. campestris interface-induced microRNAs are distinguished by a common cis-regulatory element. This element is identical to a conserved upstream sequence element (USE) used by plant small nuclear RNA loci. The properties of the interface-induced microRNA primary transcripts strongly suggest that they are produced via U6-like transcription by RNA polymerase III. The USE promotes accumulation of interface-induced miRNAs (IIMs) in a heterologous system. This promoter element distinguishes C. campestris IIM loci from other plant small RNAs. Our data suggest that C. campestris IIMs are produced in a manner distinct from canonical miRNAs. All confirmed C. campestris microRNAs with documented trans-species activity are interface-induced and possess these features. We speculate that RNA polymerase III transcription of IIMs may allow these miRNAs to be exported to hosts.

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.

Trans-Species Mobility of RNA Interference between Plants and Associated Organisms.

Trans-species RNA interference (RNAi) occurs naturally when small RNAs (sRNAs) silence genes in species different from their origin. This phenomenon has been observed between plants and various organisms including fungi, animals and other plant species. Understanding the mechanisms used in natural cases of trans-species RNAi, such as sRNA processing and movement, will enable more effective development of crop protection methods using host-induced gene silencing (HIGS). Recent progress has been made in understanding the mechanisms of cell-to-cell and long-distance movement of sRNAs within individual plants. This increased understanding of endogenous plant sRNA movement may be translatable to trans-species sRNA movement. Here, we review diverse cases of natural trans-species RNAi focusing on current theories regarding intercellular and long-distance sRNA movement. We also touch on trans-species sRNA evolution, highlighting its research potential and its role in improving the efficacy of HIGS.

MicroRNAs from the parasitic plant Cuscuta campestris target host messenger RNAs.

Dodders (Cuscuta spp.) are obligate parasitic plants that obtain water and nutrients from the stems of host plants via specialized feeding structures called haustoria. Dodder haustoria facilitate bidirectional movement of viruses, proteins and mRNAs between host and parasite, but the functional effects of these movements are not known. Here we show that Cuscuta campestris haustoria accumulate high levels of many novel microRNAs (miRNAs) while parasitizing Arabidopsis thaliana. Many of these miRNAs are 22 nucleotides in length. Plant miRNAs of this length are uncommon, and are associated with amplification of target silencing through secondary short interfering RNA (siRNA) production. Several A. thaliana mRNAs are targeted by 22-nucleotide C. campestris miRNAs during parasitism, resulting in mRNA cleavage, secondary siRNA production, and decreased mRNA accumulation. Hosts with mutations in two of the loci that encode target mRNAs supported significantly higher growth of C. campestris. The same miRNAs that are expressed and active when C. campestris parasitizes A. thaliana are also expressed and active when it infects Nicotiana benthamiana. Homologues of target mRNAs from many other plant species also contain the predicted target sites for the induced C. campestris miRNAs. These data show that C. campestris miRNAs act as trans-species regulators of host-gene expression, and suggest that they may act as virulence factors during parasitism.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants.

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5'-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Expression and processing of polycistronic artificial microRNAs and trans-acting siRNAs from transiently introduced transgenes in Solanum lycopersicum and Nicotiana benthamiana.

Targeted gene silencing using small regulatory RNAs is a widely used technique for genetic studies in plants. Artificial microRNAs are one common approach, as they have the advantage of producing just a single functional small RNA, which can be designed for high target specificity and low off-target effects. Simultaneous silencing of multiple targets with artificial microRNAs can be achieved by producing polycistronic microRNA precursors. Alternatively, specialized trans-acting short interfering RNA (tasiRNA) precursors can be designed to produce several specific tasiRNAs at once. Here we tested several artificial microRNA- and tasiRNA-based methods for multiplexed gene silencing in Solanum lycopersicum (tomato) and Nicotiana benthamiana. All analyses used transiently expressed transgenes delivered by infiltration of leaves with Agrobacterium tumefacians. Small RNA sequencing analyses revealed that many previously described approaches resulted in poor small RNA processing. The 5'-most microRNA precursor hairpins on polycistronic artificial microRNA precursors were generally processed more accurately than precursors at the 3'-end. Polycistronic artificial microRNAs where the hairpin precursors were separated by transfer RNAs had the best processing precision. Strikingly, artificial tasiRNA precursors failed to be processed in the expected phased manner in our system. These results highlight the need for further development of multiplexed artificial microRNA and tasiRNA strategies. The importance of small RNA sequencing, as opposed to single-target assays such as RNA blots or real-time polymerase chain reaction, is also discussed.

Mechanisms of resistance and virulence in parasitic plant-host interactions.

Parasitic plants pose a major biotic threat to plant growth and development and lead to losses in crop productivity of billions of USD annually. By comparison with "normal" autotrophic plants, parasitic plants live a heterotrophic lifestyle and rely on water, solutes and to a greater (holoparasitic plants) or lesser extent (hemiparasitic plants) on sugars from other host plants. Most hosts are unable to detect an infestation by plant parasites or unable to fend off these parasitic invaders. However, a few hosts have evolved defense strategies to avoid infestation or protect themselves actively post-attack often leading to full or partial resistance. Here, we review the current state of our understanding of the defense strategies to plant parasitism used by host plants with emphasis on the active molecular resistance mechanisms. Furthermore, we outline the perspectives and the potential of future studies that will be indispensable to develop and breed resistant crops.

Exchange of Small Regulatory RNAs between Plants and Their Pests.

Regulatory small RNAs are well known as antiviral agents, regulators of gene expression, and defenders of genome integrity in plants. Several studies over the last decade have also shown that some small RNAs are exchanged between plants and their pathogens and parasites. Naturally occurring trans-species small RNAs are used by host plants to silence mRNAs in pathogens. These gene-silencing events are thought to be detrimental to the pathogen and beneficial to the host. Conversely, trans-species small RNAs from pathogens and parasites are deployed to silence host mRNAs; these events are thought to be beneficial for the pests. The natural ability of plants to exchange small RNAs with invading eukaryotic organisms can be exploited to provide disease resistance. This review gives an overview of the current state of trans-species small RNA research in plants and discusses several outstanding questions for future research.

Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data.

MicroRNAs (miRNAs) are ∼21-nucleotide-long regulatory RNAs that arise from endonucleolytic processing of hairpin precursors. Many function as essential posttranscriptional regulators of target mRNAs and long noncoding RNAs. Alongside miRNAs, plants also produce large numbers of short interfering RNAs (siRNAs), which are distinguished from miRNAs primarily by their biogenesis (typically processed from long double-stranded RNA instead of single-stranded hairpins) and functions (typically via roles in transcriptional regulation instead of posttranscriptional regulation). Next-generation DNA sequencing methods have yielded extensive data sets of plant small RNAs, resulting in many miRNA annotations. However, it has become clear that many miRNA annotations are questionable. The sheer number of endogenous siRNAs compared with miRNAs has been a major factor in the erroneous annotation of siRNAs as miRNAs. Here, we provide updated criteria for the confident annotation of plant miRNAs, suitable for the era of "big data" from DNA sequencing. The updated criteria emphasize replication and the minimization of false positives, and they require next-generation sequencing of small RNAs. We argue that improved annotation systems are needed for miRNAs and all other classes of plant small RNAs. Finally, to illustrate the complexities of miRNA and siRNA annotation, we review the evolution and functions of miRNAs and siRNAs in plants.

Analysis of RDR1/RDR2/RDR6-independent small RNAs in Arabidopsis thaliana improves MIRNA annotations and reveals unexplained types of short interfering RNA loci.

Plant small RNAs (sRNAs) modulate key physiological mechanisms through post-transcriptional and transcriptional silencing of gene expression. Small RNAs fall into two major categories: those are reliant on RNA-dependent RNA polymerases (RDRs) for biogenesis and those that are not. Known RDR1/2/6-dependent sRNAs include phased and repeat-associated short interfering RNAs, while known RDR1/2/6-independent sRNAs are primarily microRNAs (miRNA) and other hairpin-derived sRNAs. In this study we produced and analyzed sRNA-seq libraries from rdr1/rdr2/rdr6 triple mutant plants. We found 58 previously annotated miRNA loci that were reliant on RDR1, -2, or -6 function, casting doubt on their classification. We also found 38 RDR1/2/6-independent sRNA loci that are not MIRNAs or otherwise hairpin-derived, and did not fit into other known paradigms for sRNA biogenesis. These 38 sRNA-producing loci have as-yet-undescribed biogenesis mechanisms, and are frequently located in the vicinity of protein-coding genes. Altogether, our analysis suggests that these 38 loci represent one or more undescribed types of sRNA in Arabidopsis thaliana.

AGO4 is specifically required for heterochromatic siRNA accumulation at Pol V-dependent loci in Arabidopsis thaliana.

In plants, 24 nucleotide long heterochromatic siRNAs (het-siRNAs) transcriptionally regulate gene expression by RNA-directed DNA methylation (RdDM). The biogenesis of most het-siRNAs depends on the plant-specific RNA polymerase IV (Pol IV), and ARGONAUTE4 (AGO4) is a major het-siRNA effector protein. Through genome-wide analysis of sRNA-seq data sets, we found that AGO4 is required for the accumulation of a small subset of het-siRNAs. The accumulation of AGO4-dependent het-siRNAs also requires several factors known to participate in the effector portion of the RdDM pathway, including RNA POLYMERASE V (POL V), DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) and SAWADEE HOMEODOMAIN HOMOLOGUE 1 (SHH1). Like many AGO proteins, AGO4 is an endonuclease that can 'slice' RNAs. We found that a slicing-defective AGO4 was unable to fully recover AGO4-dependent het-siRNA accumulation from ago4 mutant plants. Collectively, our data suggest that AGO4-dependent siRNAs are secondary siRNAs dependent on the prior activity of the RdDM pathway at certain loci.

Comprehensive Annotation of Physcomitrella patens Small RNA Loci Reveals That the Heterochromatic Short Interfering RNA Pathway Is Largely Conserved in Land Plants.

Many plant small RNAs are sequence-specific negative regulators of target mRNAs and/or chromatin. In angiosperms, the two most abundant endogenous small RNA populations are usually 21-nucleotide microRNAs (miRNAs) and 24-nucleotide heterochromatic short interfering RNAs (siRNAs). Heterochromatic siRNAs are derived from repetitive regions and reinforce DNA methylation at targeted loci. The existence and extent of heterochromatic siRNAs in other land plant lineages has been unclear. Using small RNA-sequencing (RNA-seq) of the moss Physcomitrella patens, we identified 1090 loci that produce mostly 23- to 24-nucleotide siRNAs. These loci are mostly in intergenic regions with dense DNA methylation. Accumulation of siRNAs from these loci depends upon P. patens homologs of DICER-LIKE3 (DCL3), RNA-DEPENDENT RNA POLYMERASE2, and the largest subunit of DNA-DEPENDENT RNA POLYMERASE IV, with the largest subunit of a Pol V homolog contributing to expression at a smaller subset of the loci. A MINIMAL DICER-LIKE (mDCL) gene, which lacks the N-terminal helicase domain typical of DCL proteins, is specifically required for 23-nucleotide siRNA accumulation. We conclude that heterochromatic siRNAs, and their biogenesis pathways, are largely identical between angiosperms and P. patens, with the notable exception of the P. patens-specific use of mDCL to produce 23-nucleotide siRNAs.

Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAs.

Plant small RNAs are subject to various modifications. Previous reports revealed widespread 3' modifications (truncations and non-templated tailing) of plant miRNAs when the 2'-O-methyltransferase HEN1 is absent. However, non-templated nucleotides in plant heterochromatic siRNAs have not been deeply studied, especially in wild-type plants. We systematically studied non-templated nucleotide patterns in plant small RNAs by analyzing small RNA sequencing libraries from Arabidopsis, tomato, Medicago, rice, maize and Physcomitrella Elevated rates of non-templated nucleotides were observed at the 3' ends of both miRNAs and endogenous siRNAs from wild-type specimens of all species. 'Off-sized' small RNAs, such as 25 and 23 nt siRNAs arising from loci dominated by 24 nt siRNAs, often had very high rates of 3'-non-templated nucleotides. The same pattern was observed in all species that we studied. Further analysis of 24 nt siRNA clusters in Arabidopsis revealed distinct patterns of 3'-non-templated nucleotides of 23 nt siRNAs arising from heterochromatic siRNA loci. This pattern of non-templated 3' nucleotides on 23 nt siRNAs is not affected by loss of known small RNA 3'-end modifying enzymes, and may result from modifications added to longer heterochromatic siRNA precursors.

Genome-Wide Characterization of Maize Small RNA Loci and Their Regulation in the required to maintain repression6-1 (rmr6-1) Mutant and Long-Term Abiotic Stresses.

Endogenous small RNAs (sRNAs) contribute to gene regulation and genome homeostasis, but their activities and functions are incompletely known. The maize genome has a high number of transposable elements (TEs; almost 85%), some of which spawn abundant sRNAs. We performed sRNA and total RNA sequencing from control and abiotically stressed B73 wild-type plants and rmr6-1 mutants. RMR6 encodes the largest subunit of the RNA polymerase IV complex and is responsible for accumulation of most 24-nucleotide (nt) small interfering RNAs (siRNAs). We identified novel MIRNA loci and verified miR399 target conservation in maize. RMR6-dependent 23-24 nt siRNA loci were specifically enriched in the upstream region of the most highly expressed genes. Most genes misregulated in rmr6-1 did not show a significant correlation with loss of flanking siRNAs, but we identified one gene supporting existing models of direct gene regulation by TE-derived siRNAs. Long-term drought correlated with changes of miRNA and sRNA accumulation, in particular inducing down-regulation of a set of sRNA loci in the wild-typeleaf.

Analysis of complementarity requirements for plant microRNA targeting using a Nicotiana benthamiana quantitative transient assay.

MicroRNAs (miRNAs) guide RNA-induced silencing complexes to target RNAs based on miRNA-target complementarity. Using a dual-luciferase based sensor system in Nicotiana benthamiana, we quantitatively assessed the relationship between miRNA-target complementarity and silencing efficacy measured at both the RNA and protein levels, using several conserved miRNAs and their known target sites from Arabidopsis thaliana. We found that naturally occurring sites have variable efficacies attributable to their complementarity patterns. We also observed that sites with a few mismatches to the miRNA 3' regions, which are common in plants, are often equally effective and sometimes more effective than perfectly matched sites. By contrast, mismatches to the miRNA 5' regions strongly reduce or eliminate repression efficacy but are nonetheless present in several natural sites, suggesting that in some cases, suboptimal miRNA efficacies are either tolerated or perhaps selected for. Central mismatches fully abolished repression efficacy in our system, but such sites then became effective miRNA target mimics. Complementarity patterns that are functional in animals (seed sites, 3'-supplementary sites, and centered sites) did not reliably confer repression, regardless of context (3'-untranslated region or open reading frame) or measurement type (RNA or protein levels). Overall, these data provide a robust and empirical foundation for understanding, predicting, and designing functional miRNA target sites in plants.

A non-canonical plant microRNA target site.

Plant microRNAs (miRNAs) typically form near-perfect duplexes with their targets and mediate mRNA cleavage. Here, we describe an unconventional miRNA target of miR398 in Arabidopsis, an mRNA encoding the blue copper-binding protein (BCBP). BCBP mRNA carries an miR398 complementary site in its 5'-untranslated region (UTR) with a bulge of six nucleotides opposite to the 5' region of the miRNA. Despite the disruption of a target site region thought to be especially critical for function, BCBP mRNAs are cleaved by ARGONAUTE1 between nucleotides 10th and 11th, opposite to the miRNA, like conventional plant target sites. Levels of BCBP mRNAs are inversely correlated to levels of miR398 in mutants lacking the miRNA, or transgenic plants overexpressing it. Introducing two mutations that disrupt the miRNA complementarity around the cleavage site renders the target cleavage-resistant. The BCBP site functions outside of the context of the BCBP mRNA and does not depend on 5'-UTR location. Reducing the bulge does not interfere with miR398-mediated regulation and completely removing it increases the efficiency of the slicing. Analysis of degradome data and target predictions revealed that the miR398-BCBP interaction seems to be rather unique. Nevertheless, our results imply that functional target sites with non-perfect pairings in the 5' region of an ancient conserved miRNA exist in plants.

ShortStack: comprehensive annotation and quantification of small RNA genes.

Small RNA sequencing allows genome-wide discovery, categorization, and quantification of genes producing regulatory small RNAs. Many tools have been described for annotation and quantification of microRNA loci (MIRNAs) from small RNA-seq data. However, in many organisms and tissue types, MIRNA genes comprise only a small fraction of all small RNA-producing genes. ShortStack is a stand-alone application that analyzes reference-aligned small RNA-seq data and performs comprehensive de novo annotation and quantification of the inferred small RNA genes. ShortStack's output reports multiple parameters of direct relevance to small RNA gene annotation, including RNA size distributions, repetitiveness, strandedness, hairpin-association, MIRNA annotation, and phasing. In this study, ShortStack is demonstrated to perform accurate annotations and useful descriptions of diverse small RNA genes from four plants (Arabidopsis, tomato, rice, and maize) and three animals (Drosophila, mice, and humans). ShortStack efficiently processes very large small RNA-seq data sets using modest computational resources, and its performance compares favorably to previously described tools. Annotation of MIRNA loci by ShortStack is highly specific in both plants and animals. ShortStack is freely available under a GNU General Public License.

miR156 and miR390 regulate tasiRNA accumulation and developmental timing in Physcomitrella patens.

microRNA156 (miR156) affects developmental timing in flowering plants. miR156 and its target relationships with members of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) gene family appear universally conserved in land plants, but the specific functions of miR156 outside of flowering plants are unknown. We find that miR156 promotes a developmental change from young filamentous protonemata to leafy gametophores in the moss Physcomitrella patens, opposite to its role as an inhibitor of development in flowering plants. P. patens miR156 also influences accumulation of trans-acting small interfering RNAs (tasiRNAs) dependent upon a second ancient microRNA, miR390. Both miR156 and miR390 directly target a single major tasiRNA primary transcript. Inhibition of miR156 function causes increased miR390-triggered tasiRNA accumulation and decreased accumulation of tasiRNA targets. Overexpression of miR390 also caused a slower formation of gametophores, elevated miR390-triggered tasiRNA accumulation, and reduced level of tasiRNA targets. We conclude that a gene regulatory network controlled by miR156, miR390, and their targets controls developmental change in P. patens. The broad outlines and regulatory logic of this network are conserved in flowering plants, albeit with some modifications. Partially conserved small RNA networks thus influence developmental timing in plants with radically different body plans.