Unique functionality of 22-nt miRNAs in triggering RDR6-dependent siRNA biogenesis from target transcripts in Arabidopsis

Respective contributions of Arabidopsis DCL2 and DCL4 to RNA silencing

Dicer proteins are central to the different mechanisms involving RNA interference. Plants have evolved multiple DICER-LIKE (DCL) copies, thus enabling functional diversification. In Arabidopsis, DCL2 and DCL4 process double-stranded RNA into 22 and 21 nucleotide small interfering (si)RNAs, respectively, and have overlapping functions with regards to virus and transgene silencing. Nonetheless, some studies have reported that dcl2 or dcl4 single mutations are sometimes sufficient to hinder silencing. To better dissect the role of DCL2 and DCL4, we analyzed silencing kinetics and efficiencies using different transgenic systems in single and double mutant backgrounds. The results indicate that DCL2 stimulates transitivity and secondary siRNA production through DCL4 while being sufficient for silencing on its own. Notably, silencing of 35S-driven transgenes functions more efficiently in dcl4 mutants, indicating that DCL4 mostly obscures DCL2 in wild-type plants. Nonetheless, in a dcl4 mutant compromised in phloem-originating silencing, ectopically expressed DCL2 allows restoration of silencing, suggesting that DCL2 is not, or poorly, expressed in phloem. Remarkably, this ectopic DCL2 contribution to phloem-originating silencing is dependent on the activity of RNA-DEPENDENT RNA POLYMERASE6. These results indicate that, despite differences in the silencing activity of their small RNA products, DCL2 and DCL4 mostly act redundantly yet hierarchically when present simultaneously.

Suppress to Survive-Implication of Plant Viruses in PTGS

In higher plants, evolutionarily conserved processes playing an essential role during gene expression rely on small noncoding RNA molecules (sRNA). Within a wide range of sRNA-dependent cellular events, there is posttranscriptional gene silencing, the process that is activated in response to the presence of double-stranded RNAs (dsRNAs) in planta. The sequence-specific mechanism of silencing is based on RNase-mediated trimming of dsRNAs into translationally inactive short molecules. Viruses invading and replicating in host are also a source of dsRNAs and are recognized as such by cellular posttranscriptional silencing machinery leading to degradation of the pathogenic RNA. However, viruses are not totally defenseless. In parallel with evolving plant defense strategies, viruses have managed a wide range of multifunctional proteins that efficiently impede the posttranscriptional gene silencing. These viral counteracting factors are known as suppressors of RNA silencing. The aim of this review is to summarize the role and the mode of action of several functionally characterized RNA silencing suppressors encoded by RNA viruses directly involved in plant-pathogen interactions. Additionally, we point out that the widely diverse functions, structures, and modes of action of viral suppressors can be performed by different proteins, even in related viruses. All those adaptations have been evolved to achieve the same goal: to maximize the rate of viral genetic material replication by interrupting the evolutionary conserved plant defense mechanism of posttranscriptional gene silencing.

The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling

Barley (Hordeum vulgare L.) Mla alleles encode coiled-coil (CC), nucleotide binding, leucine-rich repeat (NB-LRR) receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). How Mla or NB-LRR genes in grass species are regulated at post-transcriptional level is not clear. The microRNA family, miR9863, comprises four members that differentially regulate distinct Mla alleles in barley. We show that miR9863 members guide the cleavage of Mla1 transcripts in barley, and block or reduce the accumulation of MLA1 protein in the heterologous Nicotiana benthamiana expression system. Regulation specificity is determined by variation in a unique single-nucleotide-polymorphism (SNP) in mature miR9863 family members and two SNPs in the Mla miR9863-binding site that separates these alleles into three groups. Further, we demonstrate that 22-nt miR9863s trigger the biogenesis of 21-nt phased siRNAs (phasiRNAs) and together these sRNAs form a feed-forward regulation network for repressing the expression of group I Mla alleles. Overexpression of miR9863 members specifically attenuates MLA1, but not MLA10-triggered disease resistance and cell-death signaling. We propose a key role of the miR9863 family in dampening immune response signaling triggered by a group of MLA immune receptors in barley.

Plants rely on cell-surface and intracellular immune receptors to sense pathogen invasion and to mediate defense responses. However, uncontrolled activation of immune responses is harmful to plant growth and development. Small RNAs have recently been shown to fine-tune the expression of intracellular immune receptors and contribute to the regulation of defense signaling in dicot plants, while similar processes have not been well documented in monocot grain crops, such as barley and wheat. Here, we show that, in barley, some members of the miR9863 family target a subset of Mla alleles that confer race-specific disease resistance to the powdery mildew fungus. These miRNAs act on Mla transcripts by cleavage and translational repression. Production of a type of trans-acting small RNAs, designated as phasiRNAs, enhances the effects of miRNA regulation on Mla targets. We propose that Mla-mediated immune signaling is fine-tuned by the miRNAs at later stage of MLA activation to avoid overloading of immune responses in barley cells.

Genome-wide analysis of leafbladeless1-regulated and phased small RNAs underscores the importance of the TAS3 ta-siRNA pathway to maize development

Maize leafbladeless1 (lbl1) encodes a key component in the trans-acting short-interfering RNA (ta-siRNA) biogenesis pathway. Correlated with a great diversity in ta-siRNAs and the targets they regulate, the phenotypes conditioned by mutants perturbing this small RNA pathway vary extensively across species. Mutations in lbl1 result in severe developmental defects, giving rise to plants with radial, abaxialized leaves. To investigate the basis for this phenotype, we compared the small RNA content between wild-type and lbl1 seedling apices. We show that LBL1 affects the accumulation of small RNAs in all major classes, and reveal unexpected crosstalk between ta-siRNA biogenesis and other small RNA pathways regulating transposons. Interestingly, in contrast to data from other plant species, we found no evidence for the existence of phased siRNAs generated via the one-hit model. Our analysis identified nine TAS loci, all belonging to the conserved TAS3 family. Information from RNA deep sequencing and PARE analyses identified the tasiR-ARFs as the major functional ta-siRNAs in the maize vegetative apex where they regulate expression of AUXIN RESPONSE FACTOR3 (ARF3) homologs. Plants expressing a tasiR-ARF insensitive arf3a transgene recapitulate the phenotype of lbl1, providing direct evidence that deregulation of ARF3 transcription factors underlies the developmental defects of maize ta-siRNA biogenesis mutants. The phenotypes of Arabidopsis and Medicago ta-siRNA mutants, while strikingly different, likewise result from misexpression of the tasiR-ARF target ARF3. Our data indicate that diversity in TAS pathways and their targets cannot fully account for the phenotypic differences conditioned by ta-siRNA biogenesis mutants across plant species. Instead, we propose that divergence in the gene networks downstream of the ARF3 transcription factors or the spatiotemporal pattern during leaf development in which these proteins act constitute key factors underlying the distinct contributions of the ta-siRNA pathway to development in maize, Arabidopsis, and possibly other plant species as well.

Mutations in maize leafbladeless1 (lbl1) that disrupt ta-siRNA biogenesis give rise to plants with thread-like leaves that have lost top/bottom polarity. We used genomic approaches to identify lbl1-dependent small RNAs and their targets to determine the basis for these polarity defects. This revealed substantial diversity in small RNA pathways across plant species and identified unexpected roles for LBL1 in the regulation of repetitive elements within the maize genome. We further show that only ta-siRNA loci belonging to the TAS3 family function in the maize vegetative apex. The TAS3-derived tasiR-ARFs are the main ta-siRNA active in the apex, and misregulation of their ARF3 targets emerges as the basis for the lbl1 leaf polarity defects. Supporting this, we show that plants expressing arf3a transcripts insensitive to tasiR-ARF-directed cleavage recapitulate the phenotypes observed in lbl1. The TAS3 ta-siRNA pathway, including the regulation of ARF3 genes, is conserved throughout land plant evolution, yet the phenotypes of plants defective for ta-siRNA biogenesis are strikingly different. Our data leads us to propose that divergence in the processes regulated by the ARF3 transcription factors or the spatiotemporal pattern during development in which these proteins act, underlies the diverse developmental contributions of this small RNA pathway across plants.

Genome-wide discovery and analysis of phased small interfering RNAs in Chinese sacred lotus

Phased small interfering RNA (phasiRNA) generating loci (briefly as PHAS) in plants are a novel class of genes that are normally regulated by microRNAs (miRNAs). Similar to miRNAs, phasiRNAs encoded by PHAS play important regulatory roles by targeting protein coding transcripts in plant species. We performed a genome-wide discovery of PHAS loci in Chinese sacred lotus and identified a total of 106 PHAS loci. Of these, 47 loci generate 21 nucleotide (nt) phasiRNAs and 59 loci generate 24 nt phasiRNAs, respectively. We have also identified a new putative TAS3 and a putative TAS4 loci in the lotus genome. Our results show that some of the nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance proteins and MYB transcription factors potentially generate phasiRNAs. Furthermore, our results suggest that some large subunit (LSU) rRNAs can derive putative phasiRNAs, which is potentially resulted from crosstalk between small RNA biogenesis pathways that are employed to process rRNAs and PHAS loci, respectively. Some of the identified phasiRNAs have putative trans-targets with less than 4 mismatches, suggesting that the identified PHAS are involved in many different pathways. Finally, the discovery of 24 nt PHAS in lotus suggests that there are 24 nt PHAS in dicots.

Epigenetic regulation in plants

The study of epigenetics in plants has a long and rich history, from initial descriptions of non-Mendelian gene behaviors to seminal discoveries of chromatin-modifying proteins and RNAs that mediate gene silencing in most eukaryotes, including humans. Genetic screens in the model plant Arabidopsis have been particularly rewarding, identifying more than 130 epigenetic regulators thus far. The diversity of epigenetic pathways in plants is remarkable, presumably contributing to the phenotypic plasticity of plant postembryonic development and the ability to survive and reproduce in unpredictable environments.

Epigenetic regulation in plants

The study of epigenetics in plants has a long and rich history, from initial descriptions of non-Mendelian gene behaviors to seminal discoveries of chromatin-modifying proteins and RNAs that mediate gene silencing in most eukaryotes, including humans. Genetic screens in the model plant Arabidopsis have been particularly rewarding, identifying more than 130 epigenetic regulators thus far. The diversity of epigenetic pathways in plants is remarkable, presumably contributing to the phenotypic plasticity of plant postembryonic development and the ability to survive and reproduce in unpredictable environments.

An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

Identification of novel and conserved microRNAs in Coffea canephora and Coffea arabica

As microRNAs (miRNAs) are important regulators of many biological processes, a series of small RNAomes from plants have been produced in the last decade. However, miRNA data from several groups of plants are still lacking, including some economically important crops. Here microRNAs from Coffea canephora leaves were profiled and 58 unique sequences belonging to 33 families were found, including two novel microRNAs that have never been described before in plants. Some of the microRNA sequences were also identified in Coffea arabica that, together with C. canephora, correspond to the two major sources of coffee production in the world. The targets of almost all miRNAs were also predicted on coffee expressed sequences. This is the first report of novel miRNAs in the genus Coffea, and also the first in the plant order Gentianales. The data obtained establishes the basis for the understanding of the complex miRNA-target network on those two important crops.

Genome-wide discovery and differential regulation of conserved and novel microRNAs in chickpea via deep sequencing

MicroRNAs (miRNAs) are essential components of complex gene regulatory networks that orchestrate plant development. Although several genomic resources have been developed for the legume crop chickpea, miRNAs have not been discovered until now. For genome-wide discovery of miRNAs in chickpea (Cicer arietinum), we sequenced the small RNA content from seven major tissues/organs employing Illumina technology. About 154 million reads were generated, which represented more than 20 million distinct small RNA sequences. We identified a total of 440 conserved miRNAs in chickpea based on sequence similarity with known miRNAs in other plants. In addition, 178 novel miRNAs were identified using a miRDeep pipeline with plant-specific scoring. Some of the conserved and novel miRNAs with significant sequence similarity were grouped into families. The chickpea miRNAs targeted a wide range of mRNAs involved in diverse cellular processes, including transcriptional regulation (transcription factors), protein modification and turnover, signal transduction, and metabolism. Our analysis revealed several miRNAs with differential spatial expression. Many of the chickpea miRNAs were expressed in a tissue-specific manner. The conserved and differential expression of members of the same miRNA family in different tissues was also observed. Some of the same family members were predicted to target different chickpea mRNAs, which suggested the specificity and complexity of miRNA-mediated developmental regulation. This study, for the first time, reveals a comprehensive set of conserved and novel miRNAs along with their expression patterns and putative targets in chickpea, and provides a framework for understanding regulation of developmental processes in legumes.

A stearoyl-acyl carrier protein desaturase, NbSACPD-C, is critical for ovule development in Nicotiana benthamiana

Stearoyl-acyl carrier protein desaturase (SACPD) activity is essential for production of the major unsaturated fatty acids (UFAs) in plant lipids. We report here the characterization of three SACPD genes from Nicotiana benthamiana, NbSACPD-A, -B, and -C. All three genes share high similarity to AtSSI2/FAB2 (Suppressor of Salicylic acid-Insensitivity2/Fatty Acid Biosynthesis2), the primary SACPD isoform in Arabidopsis. Knocking down the expression of individual or combinations of NbSACPDs by an artificial microRNA approach resulted in significantly reduced accumulation of 18C UFAs and elevated levels of 18:0-FA (Fatty acids) in leaves, indicating that all three genes participated in fatty acid desaturation. The triple knockdown (KD) plants displayed severe growth phenotypes, including spontaneous cell death and dwarfing. While no vegetative morphologic abnormality was observed in NbSACPD-A, -B, or -C KD plants, strikingly, NbSACPD-C KD plants produced small fruits with aborted ovules. Reciprocal crosses with wild-type and NbSACPD-C KD plants revealed that knocking down NbSACPD-C expression caused female, but not male, sterility. Furthermore, arrested ovule development and significantly altered lipid composition in ovaries were observed in NbSACPD-C KD plants, consistent with the predominant NbSACPD-C expression in ovules. The ovule development defect was fully complemented by coexpressing an amiRNA-resistant NbSACPD-C variant in the NbSACPD-C KD background, further supporting a specific requirement for NbSACPD-C in female fertility. Our results thus indicated that NbSACPD-C plays a critical role maintaining membrane lipid composition in ovule development for female fertility in N. benthamiana, complementing and extending prior understanding on the well-demonstrated roles of SACPDs in biotic and abiotic stresses.

MicroRNAs suppress NB domain genes in tomato that confer resistance to Fusarium oxysporum

MicroRNAs (miRNAs) suppress the transcriptional and post-transcriptional expression of genes in plants. Several miRNA families target genes encoding nucleotide-binding site–leucine-rich repeat (NB-LRR) plant innate immune receptors. The fungus Fusarium oxysporum f. sp. lycopersici causes vascular wilt disease in tomato. We explored a role for miRNAs in tomato defense against F. oxysporum using comparative miRNA profiling of susceptible (Moneymaker) and resistant (Motelle) tomato cultivars. slmiR482f and slmiR5300 were repressed during infection of Motelle with F. oxysporum. Two predicted mRNA targets each of slmiR482f and slmiR5300 exhibited increased expression in Motelle and the ability of these four targets to be regulated by the miRNAs was confirmed by co-expression in Nicotiana benthamiana. Silencing of the targets in the resistant Motelle cultivar revealed a role in fungal resistance for all four genes. All four targets encode proteins with full or partial nucleotide-binding (NB) domains. One slmiR5300 target corresponds to tm-2, a susceptible allele of the Tomato Mosaic Virus resistance gene, supporting functions in immunity to a fungal pathogen. The observation that none of the targets correspond to I-2, the only known resistance (R) gene for F. oxysporum in tomato, supports roles for additional R genes in the immune response. Taken together, our findings suggest that Moneymaker is highly susceptible because its potential resistance is insufficiently expressed due to the action of miRNAs.

Fusarium oxysporum is a fungal pathogen that represents a species complex, with members that infect numerous crops. In spite of its importance to agriculture, very little is known about roles of small RNAs in plant immunity against F. oxysporum. In this study, we set up a screen for tomato microRNAs (miRNAs) that correlate with resistance to F. oxysporum f.sp. lycopersici by performing deep sequencing of small RNAs from a resistant and susceptible tomato cultivar. We focused on two miRNAs that are uniquely down-regulated in the resistant cultivar during fungal infection. All predicted targets of these miRNAs encode proteins with NB domains, a motif associated with pathogen resistance in plants. Using a heterologous system, we validated that the miRNAs could regulate expression of four targets. Silencing of the target genes in tomato resulted in decreased immunity to F. oxysporum in the normally resistant cultivar. The finding that none of our targets correspond to I-2, the only known resistance (R) gene for F. oxysporum in tomato, supports roles for additional R genes in the immune response. Our results suggest that the potential resistance of the susceptible cultivar is insufficiently expressed due to the action of miRNAs.

Specific argonautes selectively bind small RNAs derived from potato spindle tuber viroid and attenuate viroid accumulation in vivo

The identification of viroid-derived small RNAs (vd-sRNAs) of 21 to 24 nucleotides (nt) in plants infected by viroids (infectious non-protein-coding RNAs of just 250 to 400 nt) supports their targeting by Dicer-like enzymes, the first host RNA-silencing barrier. However, whether viroids, like RNA viruses, are also targeted by the RNA-induced silencing complex (RISC) remains controversial. At the RISC core is one Argonaute (AGO) protein that, guided by endogenous or viral sRNAs, targets complementary RNAs. To examine whether AGO proteins also load vd-sRNAs, leaves of Nicotiana benthamiana infected by potato spindle tuber viroid (PSTVd) were agroinfiltrated with plasmids expressing epitope-tagged versions of AGO1, AGO2, AGO3, AGO4, AGO5, AGO6, AGO7, AGO9, and AGO10 from Arabidopsis thaliana. Immunoprecipitation analyses of the agroinfiltrated halos revealed that all AGOs except AGO6, AGO7, and AGO10 associated with vd-sRNAs: AGO1, AGO2, and AGO3 preferentially with those of 21 and 22 nt, while AGO4, AGO5, and AGO9 additionally bound those of 24 nt. Deep-sequencing analyses showed that sorting of vd-sRNAs into AGO1, AGO2, AGO4, and AGO5 depended essentially on their 5'-terminal nucleotides, with the profiles of the corresponding AGO-loaded vd-sRNAs adopting specific hot spot distributions along the viroid genome. Furthermore, agroexpression of AGO1, AGO2, AGO4, and AGO5 on PSTVd-infected tissue attenuated the level of the genomic RNAs, suggesting that they, or their precursors, are RISC targeted. In contrast to RNA viruses, PSTVd infection of N. benthamiana did not affect miR168-mediated regulation of the endogenous AGO1, which loaded vd-sRNAs with specificity similar to that of its A. thaliana counterpart. Importance: To contain invaders, particularly RNA viruses, plants have evolved an RNA-silencing mechanism relying on the generation by Dicer-like (DCL) enzymes of virus-derived small RNAs of 21 to 24 nucleotides (nt) that load and guide Argonaute (AGO) proteins to target and repress viral RNA. Viroids, despite their minimal genomes (non-protein-coding RNAs of only 250 to 400 nt), infect and incite disease in plants. The accumulation in these plants of 21- to 24-nt viroid-derived small RNAs (vd-sRNAs) supports the notion that DCLs also target viroids but does not clarify whether vd-sRNAs activate one or more AGOs. Here, we show that in leaves of Nicotiana benthamiana infected by potato spindle tuber viroid, the endogenous AGO1 and distinct AGOs from Arabidopsis thaliana that were overexpressed were associated with vd-sRNAs displaying the same properties (5'-terminal nucleotide and size) previously established for endogenous and viral small RNAs. Overexpression of AGO1, AGO2, AGO4, and AGO5 attenuated viroid accumulation, supporting their role in antiviroid defense.

Endogenous small-noncoding RNAs and their roles in chilling response and stress acclimation in Cassava

BACKGROUND

Small noncoding RNA (sncRNA), including microRNAs (miRNAs) and endogenous small-interfering RNAs (endo-siRNAs) are key gene regulators in eukaryotes, playing critical roles in plant development and stress tolerance. Trans-acting siRNAs (ta-siRNAs), which are secondary siRNAs triggered by miRNAs, and siRNAs from natural antisense transcripts (nat-siRNAs) are two well-studied classes of endo-siRNAs.

RESULTS

In order to understand sncRNAs' roles in plant chilling response and stress acclimation, we performed a comprehensive study of miRNAs and endo-siRNAs in Cassava (Manihot esculenta), a major source of food for the world populations in tropical regions. Combining Next-Generation sequencing and computational and experimental analyses, we profiled and characterized sncRNA species and mRNA genes from the plants that experienced severe and moderate chilling stresses, that underwent further severe chilling stress after chilling acclimation at moderate stress, and that grew under the normal condition. We also included castor bean (Ricinus communis) in our study to understand conservation of sncRNAs. In addition to known miRNAs, we identified 32 (22 and 10) novel miRNAs as well as 47 (26 and 21) putative secondary siRNA-yielding and 8 (7 and 1) nat-siRNA-yielding candidate loci in Cassava and castor bean, respectively. Among the expressed sncRNAs, 114 miRNAs, 12 ta-siRNAs and 2 nat-siRNAs showed significant expression changes under chilling stresses.

CONCLUSION

Systematic and computational analysis of microRNAome and experimental validation collectively showed that miRNAs, ta-siRNAs, and possibly nat-siRNAs play important roles in chilling response and chilling acclimation in Cassava by regulating stress-related pathways, e.g. Auxin signal transduction. The conservation of these sncRNA might shed lights on the role of sncRNA-mediated pathways affected by chilling stress and stress acclimation in Euphorbiaceous plants.

Analysis of phytoplasma-responsive sRNAs provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

The yellow dwarf disease associated with phytoplasmas is one of the most devastating diseases of mulberry and the pathogenesis involved in the disease is poorly understood. To analyze the molecular mechanisms mediating gene expression in mulberry-phytoplasma interaction, the comprehensive sRNA changes of mulberry leaf in response to phytoplasma-infection were examined. A total of 164 conserved miRNAs and 23 novel miRNAs were identified, and 62 conserved miRNAs and 13 novel miRNAs were found to be involved in the response to phytoplasma-infection. Meanwhile, target genes of the responsive miRNAs were identified by sequencing of the degradome library. In addition, the endogenous siRNAs were sequenced, and their expression profiles were characterized. Interestingly, we found that phytoplasma infection induced the accumulation of mul-miR393-5p which was resulted from the increased transcription of MulMIR393A, and mul-miR393-5p most likely initiate the biogenesis of siRNAs from TIR1 transcript. Based on the results, we can conclude that phytoplasma-responsive sRNAs modulate multiple hormone pathways and play crucial roles in the regulation of development and metabolism. These responsive sRNAs may work cooperatively in the response to phytoplasma-infection and be responsible for some symptoms in the infected plants.

Assembly and function of small RNA – Argonaute protein complexes

Small RNAs such as microRNAs (miRNAs), short interfering RNAs (siRNAs) or Piwi-interacting RNAs (piRNAs) are important regulators of gene expression in various organisms. Small RNAs bind to a member of the Argonaute protein family and are incorporated into larger structures that mediate diverse gene silencing events. The loading of Argonaute proteins with small RNAs is aided by a number of auxiliary factors as well as ATP hydrolysis. This review will focus on the mechanisms of Argonaute loading in different organisms. Furthermore, we highlight the versatile functions of small RNA-Argonaute protein complexes in organisms from all three kingdoms of life.

Whole-genome discovery of miRNAs and their targets in wheat (Triticum aestivum L.)

BACKGROUND

MicroRNAs (miRNAs) are small, non-coding RNAs playing essential roles in plant growth, development, and stress responses. Sequencing of small RNAs is a starting point for understanding their number, diversity, expression and possible roles in plants.

RESULTS

In this study, we conducted a genome-wide survey of wheat miRNAs from 11 tissues, characterizing a total of 323 novel miRNAs belonging to 276 families in wheat. A miRNA conservation analysis identified 191 wheat-specific miRNAs, 2 monocot-specific miRNAs, and 30 wheat-specific variants from 9 highly conserved miRNA families. To understand possible roles of wheat miRNAs, we determined 524 potential targets for 124 miRNA families through degradome sequencing, and cleavage of a subset of them was validated via 5' RACE. Based on the genome-wide identification and characterization of miRNAs and their associated target genes, we further identified 64 miRNAs preferentially expressing in developing or germinating grains, which could play important roles in grain development.

CONCLUSION

We discovered 323 wheat novel miRNAs and 524 target genes for 124 miRNA families in a genome-wide level, and our data will serve as a foundation for future research into the functional roles of miRNAs in wheat.

New generation of artificial MicroRNA and synthetic trans-acting small interfering RNA vectors for efficient gene silencing in Arabidopsis

Artificial microRNAs (amiRNAs) and synthetic trans-acting small interfering RNAs (syn-tasiRNAs) are used for small RNA-based, specific gene silencing or knockdown in plants. Current methods to generate amiRNA or syn-tasiRNA constructs are not well adapted for cost-effective, large-scale production or for multiplexing to specifically suppress multiple targets. Here, we describe simple, fast, and cost-effective methods with high-throughput capability to generate amiRNA and multiplexed syn-tasiRNA constructs for efficient gene silencing in Arabidopsis (Arabidopsis thaliana) and other plant species. amiRNA or syn-tasiRNA inserts resulting from the annealing of two overlapping and partially complementary oligonucleotides are ligated directionally into a zero background BsaI/ccdB-based expression vector. BsaI/ccdB vectors for amiRNA or syn-tasiRNA cloning and expression contain a modified version of Arabidopsis MIR390a or TAS1c precursors, respectively, in which a fragment of the endogenous sequence was substituted by a ccdB cassette flanked by two BsaI sites. Several amiRNA and syn-tasiRNA sequences designed to target one or more endogenous genes were validated in transgenic plants that (1) exhibited the expected phenotypes predicted by loss of target gene function, (2) accumulated high levels of accurately processed amiRNAs or syn-tasiRNAs, and (3) had reduced levels of the corresponding target RNAs.

Endogenous small RNA clusters in plants

In plants, small RNAs (sRNAs) usually refer to non-coding RNAs (ncRNAs) with lengths of 20–24 nucleotides. sRNAs are involved in the regulation of many essential processes related to plant development and environmental responses. sRNAs in plants are mainly grouped into microRNAs (miRNAs) and small interfering RNAs (siRNAs), and the latter can be further classified into trans-acting siRNAs (ta-siRNAs), repeat-associated siRNAs (ra-siRNAs), natural anti-sense siRNAs (nat-siRNAs), etc. Many sRNAs exhibit a clustered distribution pattern in the genome. Here, we summarize the features and functions of cluster-distributed sRNAs, aimed to not only provide a thorough picture of sRNA clusters (SRCs) in plants, but also shed light on the identification of new classes of functional sRNAs.

Differential sRNA regulation in leaves and roots of sugarcane under water depletion

Plants have developed multiple regulatory mechanisms to respond and adapt to stress. Drought stress is one of the major constraints to agricultural productivity worldwide and recent reports have highlighted the importance of plant sRNA in the response and adaptation to water availability. In order to increase our understanding of the roles of sRNA in response to water depletion, cultivars of sugarcane were submitted to treatment of ceasing drip irrigation for 24 hours. Deep sequencing analysis was carried out to identify the sRNA regulated in leaves and roots of sugarcane cultivars with different drought sensitivities. The pool of sRNA selected allowed the analysis of different sRNA classes (miRNA and siRNA). Twenty-eight and 36 families of conserved miRNA were identified in leaf and root libraries, respectively. Dynamic regulation of miRNA was observed and the expression profiles of eight miRNA were verified in leaf samples from three biological replicates by stem-loop qRT-PCR assay using the cultivars: SP90-1638--sensitive cultivar--and SP83-2847 and SP83-5073--tolerant cultivars. Altered miRNA regulation was correlated with changes in mRNA levels of specific targets. Two leaf libraries from individual sugarcane cultivars with contrasting drought-tolerance properties were also analyzed. An enrichment of 22-nt sRNA species was observed in leaf libraries. 22-nt miRNA triggered siRNA production by cleavage of their targets in response to water depletion. A number of genes of the sRNA biogenesis pathway were down-regulated in tolerant genotypes and up-regulated in sensitive in response to water depletion treatment. Our analysis contributes to increase the knowledge on the roles of sRNA in sugarcane submitted to water depletion.

miRNAs trigger widespread epigenetically activated siRNAs from transposons in Arabidopsis

In plants, post-transcriptional gene silencing (PTGS) is mediated by DICER-LIKE1 (DCL1)-dependent miRNAs, that also trigger 21-nt secondary siRNA via RNA DEPENDENT RNA POLYMERASE6 (RDR6), DCL4, and ARGONAUTE1 (AGO1)^1–3, while transcriptional gene silencing (TGS) of transposons is mediated by 24-nt heterochromatic (het)siRNA RDR2, DCL3 and AGO4⁴. Transposons can also give rise to abundant 21-nt “epigenetically activated” small interfering RNAs (easiRNAs) in DECREASE IN DNA METHYLATION1 (ddm1) and DNA METHYLTRANSFERASE1 (met1) mutants, as well as in the vegetative nucleus of pollen grains⁵, and in dedifferentiated plant cell cultures⁶. Here we show that easiRNAs resemble secondary siRNAs, in that thousands of transposon transcripts are specifically targeted by more than fifty miRNAs for cleavage and processing by RDR6. Loss of RDR6, DCL4 or DCL1 in a ddm1 background results in loss of 21-nt easiRNA, and severe infertility, but 24-nt hetsiRNA are partially restored, supporting an antagonistic relationship between PTGS and TGS. Thus miRNA-directed easiRNA biogenesis is a latent mechanism that specifically targets transposon transcripts, but only when they are epigenetically reactivated during reprogramming of the germline. This ancient recognition mechanism may have been retained both by transposons to evade long-term heterochromatic silencing, and by their hosts for genome defence.

MicroRNA-target binding structures mimic microRNA duplex structures in humans

Traditionally, researchers match a microRNA guide strand to mRNA sequences using sequence comparisons to predict its potential target genes. However, many of the predictions can be false positives due to limitations in sequence comparison alone. In this work, we consider the association of two related RNA structures that share a common guide strand: the microRNA duplex and the microRNA-target binding structure. We have analyzed thousands of such structure pairs and found many of them share high structural similarity. Therefore, we conclude that when predicting microRNA target genes, considering just the microRNA guide strand matches to gene sequences may not be sufficient--the microRNA duplex structure formed by the guide strand and its companion passenger strand must also be considered. We have developed software to translate RNA binding structure into encoded representations, and we have also created novel automatic comparison methods utilizing such encoded representations to determine RNA structure similarity. Our software and methods can be utilized in the other RNA secondary structure comparisons as well.

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.

Global and local perturbation of the tomato microRNA pathway by a trans-activated DICER-LIKE 1 mutant

DICER-like 1 (DCL1) is a major player in microRNA (miRNA) biogenesis and accordingly, its few known loss-of-function mutants are either lethal or display arrested development. Consequently, generation of dcl1 mutants by reverse genetics and functional analysis of DCL1 in late-developing organs are challenging. Here, these challenges were resolved through the unique use of trans-activated RNA interference. Global, as well as organ-specific tomato DCL1 (SlDCL1) silencing was induced by crossing the generated responder line (OP:SlDCL1IR) with the appropriate driver line. Constitutive trans-activation knocked down SlDCL1 levels by ~95%, resulting in severe abnormalities including post-germination growth arrest accompanied by decreased miRNA and 21-nucleotide small RNA levels, but prominently elevated levels of 22-nucleotide small RNAs. The increase in the 22-nucleotide small RNAs was correlated with specific up-regulation of SlDCL2b and SlDCL2d, which are probably involved in their biogenesis. Leaf- and flower-specific OP:SlDCL1IR trans-activation inhibited blade outgrowth, induced premature bud senescence and produced pale petals, respectively, emphasizing the importance of SlDCL1-dependent small RNAs in these processes. Together, these results establish OP:SlDCL1IR as an efficient tool for analysing processes regulated by SlDCL1-mediated gene regulation in tomato.

Biogenesis of trans-acting siRNAs, endogenous secondary siRNAs in plants

Trans-acting small interfering RNAs (tasiRNAs) are plant-specific endogenous siRNAs that control non-identical mRNAs via cleavage. The production of tasiRNAs is triggered by cleavage of capped and polyadenylated primary TAS transcripts (pri-TASs) by specific miRNAs. Following miRNA-directed cleavage, either 5' or 3' cleavage fragments are converted into double-stranded RNAs (dsRNAs) by RNA-DEPENDENT RNA POLYMERASE 6. The dsRNAs are processed to tasiRNAs by DICER-LIKE 4 in a phasing manner. There are two forms of pri-TASs; One has a single miRNA target site that is targeted by 22-nucleotide microRNAs, and the other has two miR390 target sites. Secondary siRNAs that are important for the amplification of RNA silencing are defined as siRNAs whose production is initiated by the cleavage of primary small RNA-containing RNA-induced silencing complexes. Thus, tasiRNA production is a model system of secondary siRNA production in plants. This review focuses on the production of tasiRNAs that are endogenous secondary siRNAs.

The Arabidopsis miR472-RDR6 silencing pathway modulates PAMP- and effector-triggered immunity through the post-transcriptional control of disease resistance genes

RNA-DEPENDENT RNA POLYMERASE 6 (RDR6) is a key RNA silencing factor initially characterized in transgene silencing and virus resistance. This enzyme also contributes to the biosynthesis of endogenous short interfering RNAs (siRNAs) from non-coding RNAs, transposable elements and protein-coding transcripts. One class of protein-coding transcripts that have recently emerged as major sources of RDR6-dependent siRNAs are nucleotide-binding leucine-rich repeat (NB-LRR) proteins, a family of immune-receptors that perceive specific pathogen effector proteins and mount Effector-Triggered Immunity (ETI). Nevertheless, the dynamic post-transcriptional control of NB-LRR transcripts during the plant immune response and the functional relevance of NB-LRRs in signaling events triggered by Pathogen-Associated Molecular Patterns (PAMPs) remain elusive. Here, we show that PTI is constitutive and sensitized in the Arabidopsis rdr6 loss-of-function mutant, implicating RDR6 as a novel negative regulator of PTI. Accordingly, rdr6 mutant exhibits enhanced basal resistance towards a virulent Pseudomonas syringae strain. We further provide evidence that dozens of CC-NB-LRRs (CNLs), including the functionally characterized RPS5 gene, are post-transcriptionally controlled by RDR6 both constitutively and during PTI. These CNL transcripts are also regulated by the Arabidopsis microRNA miR472 and knock-down of this miRNA recapitulates the PTI and basal resistance phenotypes observed in the rdr6 mutant background. Furthermore, both miR472 and rdr6 mutants were more resistant to Pto DC3000 expressing AvrPphB, a bacterial effector recognized by the disease resistance protein RPS5, whereas transgenic plants overexpressing miR472 were more susceptible to this bacterial strain. Finally, we show that the enhanced basal and RPS5-mediated resistance phenotypes observed in the rdr6 mutant are dependent on the proper chaperoning of NB-LRR proteins, and might therefore be due to the enhanced accumulation of CNL proteins whose cognate mRNAs are no longer controlled by RDR6-dependent siRNAs. Altogether, this study supports a model whereby the miR472- and RDR6-mediated silencing pathway represents a key regulatory checkpoint modulating both PTI and ETI responses through the post-transcriptional control of disease resistance genes.

Virus resistance relies in some plant-viral interactions on the RNA-DEPENDANT RNA POLYMERASE 6 (RDR6), a major actor of RNA silencing that acts at the post-transcriptional level. Here, we demonstrate that RDR6 also plays a role in basal defense and race-specific resistance. RDR6 and the microRNA miR472, which targets the mRNAs of disease resistance genes of coiled-coil nucleotide-binding leucine-rich-repeats family (e.g. RPS5), act in cooperation to control post-transcriptionally these immune receptors. Induction of these resistance genes is primed in rdr6- and miR472-elicited mutants and this effect is associated with an enhanced basal and race-specific immunity in these backgrounds.

Novel insights from non-conserved microRNAs in plants

Plant microRNAs (miRNAs), a class of small non-coding regulatory RNAs, are canonically 20-24 nucleotides in length and bind to complementary target RNA sequences, guiding target attenuation via mRNA degradation or translation inhibition. Of the annotated miRNA families, evolutionarily conserved families have been well known to extensively regulate analogous targets and play critical roles in plant development and adaptation to adverse environments. By contrast, majority of these families that are merely present in a specific lineage or in a few closely related species have not been well functionally explored until recently. The fast-growing progresses being made in the actions of non-conserved miRNAs nowadays in diverse plant species may represent a highly promising research field in future. This review thereby summarizes the emerging advances in our understanding of the biogenesis, associated effectors, modes to targets, and biological functions of plant non-conserved miRNAs. In addition, it outlines the regulatory units recently discovered between conserved miRNAs and their alternative targets.

The Arabidopsis miR472-RDR6 silencing pathway modulates PAMP- and effector-triggered immunity through the post-transcriptional control of disease resistance genes

RNA-DEPENDENT RNA POLYMERASE 6 (RDR6) is a key RNA silencing factor initially characterized in transgene silencing and virus resistance. This enzyme also contributes to the biosynthesis of endogenous short interfering RNAs (siRNAs) from non-coding RNAs, transposable elements and protein-coding transcripts. One class of protein-coding transcripts that have recently emerged as major sources of RDR6-dependent siRNAs are nucleotide-binding leucine-rich repeat (NB-LRR) proteins, a family of immune-receptors that perceive specific pathogen effector proteins and mount Effector-Triggered Immunity (ETI). Nevertheless, the dynamic post-transcriptional control of NB-LRR transcripts during the plant immune response and the functional relevance of NB-LRRs in signaling events triggered by Pathogen-Associated Molecular Patterns (PAMPs) remain elusive. Here, we show that PTI is constitutive and sensitized in the Arabidopsis rdr6 loss-of-function mutant, implicating RDR6 as a novel negative regulator of PTI. Accordingly, rdr6 mutant exhibits enhanced basal resistance towards a virulent Pseudomonas syringae strain. We further provide evidence that dozens of CC-NB-LRRs (CNLs), including the functionally characterized RPS5 gene, are post-transcriptionally controlled by RDR6 both constitutively and during PTI. These CNL transcripts are also regulated by the Arabidopsis microRNA miR472 and knock-down of this miRNA recapitulates the PTI and basal resistance phenotypes observed in the rdr6 mutant background. Furthermore, both miR472 and rdr6 mutants were more resistant to Pto DC3000 expressing AvrPphB, a bacterial effector recognized by the disease resistance protein RPS5, whereas transgenic plants overexpressing miR472 were more susceptible to this bacterial strain. Finally, we show that the enhanced basal and RPS5-mediated resistance phenotypes observed in the rdr6 mutant are dependent on the proper chaperoning of NB-LRR proteins, and might therefore be due to the enhanced accumulation of CNL proteins whose cognate mRNAs are no longer controlled by RDR6-dependent siRNAs. Altogether, this study supports a model whereby the miR472- and RDR6-mediated silencing pathway represents a key regulatory checkpoint modulating both PTI and ETI responses through the post-transcriptional control of disease resistance genes.

The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis

In eukaryotic RNA silencing, RNase-III classes of enzymes in the Dicer family process double-stranded RNA of cellular or exogenous origin into small-RNA (sRNA) molecules. sRNAs are then loaded into effector proteins known as ARGONAUTEs (AGOs), which, as part of RNA-induced silencing complexes, target complementary RNA or DNA for silencing. Plants have evolved a large variety of pathways over the Dicer-AGO consortium, which most likely underpins part of their phenotypic plasticity. Dicer-like proteins produce all known classes of plant silencing sRNAs, which are invariably stabilized via 2'-O-methylation mediated by HUA ENHANCER 1 (HEN1), potentially amplified by the action of several RNA-dependent RNA polymerases, and function through a variety of AGO proteins. Here, we review the known characteristics and biochemical properties of the core silencing factors found in the model plant Arabidopsis thaliana. We also describe how interactions between these core factors and more specialized proteins allow the production of a plethora of silencing sRNAs involved in a large array of biological functions. We emphasize in particular the biogenesis and activities of silencing sRNAs of endogenous origin.

miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton

In this study, we characterized the miR482 family in cotton using existing small RNA datasets and the recently released draft genome sequence of Gossypium raimondii, a diploid cotton species whose progenitor is the putative contributor of the Dt (representing the D genome of tetraploid) genome of the cultivated tetraploid cotton species G. hirsutum and G. barbadense. Of the three ghr-miR482 members reported in G. hirsutum, ghr-miR482a has no homolog in G. raimondii, ghr-miR482b and ghr-miR482c each has a single homolog in G. raimondii. Gra-miR482d has five homologous loci (gra-miR482d, f-i) in G. raimondii and also exists in G. hirsutum (ghr-miR482d). A variant, miR482.2 that is a homolog of miR2118 in other species, is produced from several GHR-MIR482 loci in G. hirsutum. Approximately 12% of the G. raimondii NBS-LRR genes were predicted targets of various members of the gra-miR482 family. Based on the rationale that the regulatory relationship between miR482 and NBS-LRR genes will be conserved in G. raimondii and G. hirsutum, we investigated this relationship using G. hirsutum miR482 and G. raimondii NBS-LRR genes, which are not currently available in G. hirsutum. Ghr-miR482/miR482.2-mediated cleavage was confirmed for three of the four NBS-LRR genes analysed. As in tomato, miR482-mediated cleavage of NBS-LRR genes triggered production of phased secondary small RNAs in cotton. In seedlings of the susceptible cultivar Sicot71 (G. hirsutum) infected with the fungal pathogen Verticillium dahliae, the expression levels of ghr-miR482b/miR482b.2, ghr-miR482c and ghr-miR482d.2 were down-regulated, and several NBS-LRR targets of ghr-miR482c and ghr-miR482d were up-regulated. These results imply that, like tomato plants infected with viruses or bacteria, cotton plants are able to induce expression of NBS-LRR defence genes by suppression of the miRNA-mediated gene silencing pathway upon fungal pathogen attack.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

BACKGROUND

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

RESULTS

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

CONCLUSIONS

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

RNA silencing suppression by plant pathogens: defence, counter-defence and counter-counter-defence

RNA silencing is a central regulator of gene expression in most eukaryotes and acts both at the transcriptional level through DNA methylation and at the post-transcriptional level through direct mRNA interference mediated by small RNAs. In plants and invertebrates, the same pathways also function directly in host defence against viruses by targeting viral RNA for degradation. Successful viruses have consequently evolved diverse mechanisms to avoid silencing, most notably through the expression of viral suppressors of RNA silencing. RNA silencing suppressors have also been recently identified in plant pathogenic bacteria and oomycetes, suggesting that disruption of host silencing is a general virulence strategy across several kingdoms of plant pathogens. There is also increasing evidence that plants have evolved specific defences against RNA-silencing suppression by pathogens, providing yet another illustration of the never-ending molecular arms race between plant pathogens and their hosts.

Identification of miRNAs and their target genes in peach (Prunus persica L.) using high-throughput sequencing and degradome analysis

MicroRNAs play critical roles in various biological and metabolic processes. The function of miRNAs has been widely studied in model plants such as Arabidopsis and rice. However, the number of identified miRNAs and related miRNA targets in peach (Prunus persica) is limited. To understand further the relationship between miRNAs and their target genes during tissue development in peach, a small RNA library and three degradome libraries were constructed from three tissues for deep sequencing. We identified 117 conserved miRNAs and 186 novel miRNA candidates in peach by deep sequencing and 19 conserved miRNAs and 13 novel miRNAs were further evaluated for their expression by RT-qPCR. The number of gene targets that were identified for 26 conserved miRNA families and 38 novel miRNA candidates, were 172 and 87, respectively. Some of the identified miRNA targets were abundantly represented as conserved miRNA targets in plant. However, some of them were first identified and showed important roles in peach development. Our study provides information concerning the regulatory network of miRNAs in peach and advances our understanding of miRNA functions during tissue development.

A 22-nt artificial microRNA mediates widespread RNA silencing in Arabidopsis

It is known that 22-nucleotide (nt) microRNAs (miRNAs) derived from asymmetric duplexes trigger phased small-interfering RNA (phasiRNA) production from complementary targets. Here we investigate the efficacy of 22-nt artificial miRNA (amiRNA)-mediated RNA silencing relative to conventional hairpin RNA (hpRNA) and 21-nt amiRNA-mediated RNA silencing. CHALCONE SYNTHASE (CHS) was selected as a target in Arabidopsis thaliana due to the obvious and non-lethal loss of anthocyanin accumulation upon widespread RNA silencing. Over-expression of CHS in the pap1-D background facilitated visual detection of both local and systemic RNA silencing. RNA silencing was initiated in leaf tissues from hpRNA and amiRNA plant expression vectors under the control of an Arabidopsis RuBisCo small subunit 1A promoter (SSU). In this system, hpRNA expression triggered CHS silencing in most leaf tissues but not in roots or seed coats. Similarly, 21-nt amiRNA expression from symmetric miRNA/miRNA* duplexes triggered CHS silencing in all leaf tissues but not in roots or seed coats. However, 22-nt amiRNA expression from an asymmetric duplex triggered CHS silencing in all tissues, including roots and seed coats, in the majority of plant lines. This widespread CHS silencing required RNA-DEPENDENT RNA POLYMERASE6-mediated accumulation of phasiRNAs from the endogenous CHS transcript. These results demonstrate the efficacy of asymmetric 22-nt amiRNA-directed RNA silencing and associated phasiRNA production and activity, in mediating widespread RNA silencing of an endogenous target gene. Asymmetric 22-nt amiRNA-directed RNA silencing requires little modification of existing amiRNA technology and is expected to be effective in suppressing other genes and/or members of gene families.

Regulation of FLOWERING LOCUS T by a microRNA in Brachypodium distachyon

The highly conserved florigen gene FLOWERING LOCUS T (FT) functions at the core of the flowering pathways. Extensive studies have examined the transcriptional regulation of FT; however, other layers of FT regulation remain unclear. Here, we identified miR5200 a Pooideae-specific microRNA that is expressed in leaves and targets Brachypodium distachyon FT orthologs for mRNA cleavage. miR5200 was abundantly expressed in plants grown under short-day (SD) conditions but was dramatically repressed in plants transferred to long-day (LD) conditions. We also found that the epigenetic chromatin status, specifically the levels of histone methylation marks, at miR5200 precursor loci changed in response to daylength. Moreover, artificial interruption of miR5200 activity by target mimicry in B. distachyon altered flowering time in SD but not in LD conditions, suggesting that miR5200 functions in photoperiod-mediated flowering time regulation. Together, these findings illustrate a posttranscriptional regulation mechanism of FT and provide insights into understanding of the multiple concerted pathways for flowering time control in plants.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

The transition from primary siRNAs to amplified secondary siRNAs that regulate chalcone synthase during development of Glycine max seed coats

The I locus is a 27-kb inverted repeat cluster of chalcone synthase genes CHS1-3-4 that mediates siRNA down-regulation of CHS7 and CHS8 target mRNAs during seed development leading to yellow seed coats lacking anthocyanin pigments. Here, we report small RNA sequencing of ten stages of seed development from a few days post fertilization through maturity, revealing the amplification from primary to secondary short interfering RNAs (siRNAs) occurring during development. The young seed populations had a higher proportion of siRNAs representing the CHS1-3-4 gene family members, consistent with this region as the origin of the primary siRNAs. More intriguingly, the very young seed had a higher proportion of 22-nt CHS siRNAs than did the mid-maturation seed. We infer that the primary CHS siRNAs increase during development to levels sufficient to trigger amplification of secondary CHS siRNAs from the CHS7/8 target mRNAs, enabling the total levels of 21-nt CHS siRNAs to rise dramatically. Further, we demonstrate that the soybean system exhibits tissue-specific CHS siRNA production because primary CHS siRNA levels are not sufficient to trigger secondary amplification in tissues other than the seed coat.

Identification of miRNAs and miRNA-mediated regulatory pathways in Carica papaya

Plant microRNAs (miRNAs) post-transcriptionally regulate target gene expression to modulate growth and development and biotic and abiotic stress responses. By analyzing small RNA deep sequencing data in combination with the genome sequence, we identified 75 conserved miRNAs and 11 novel miRNAs. Their target genes were also predicted. For most conserved miRNAs, the miRNA-target pairs were conserved across plant species. In addition to these conserved miRNA-target pairs, we also identified some papaya-specific miRNA-target regulatory pathways. Both miR168 and miR530 target the Argonaute 1 gene, indicating a second autoregulatory mechanism for miRNA regulation. A non-conserved miRNA was mapped within an intron of Dicer-like 1 (DCL1), suggesting a conserved homeostatic autoregulatory mechanism for DCL1 expression. A 21-nt miRNA triggers secondary siRNA production from its target genes, nucleotide-binding site leucine-rich repeat protein genes. Certain phased-miRNAs were processed from their conserved miRNA precursors, indicating a putative miRNA evolution mechanism. In addition, we identified a Carica papaya-specific miRNA that targets an ethylene receptor gene, implying its function in the ethylene signaling pathway. This work will also advance our understanding of miRNA functions and evolution in plants.

Mechanistic insights into small RNA recognition and modification by the HEN1 methyltransferase

The HEN1 methyltransferase from Arabidopsis thaliana modifies the 3'-terminal nucleotides of small regulatory RNAs. Although it is one of the best characterized members of the 2'-O-methyltransferase family, many aspects of its interactions with the cofactor and substrate RNA remained unresolved. To better understand the substrate interactions and contributions of individual steps during HEN1 catalysis, we studied the binding and methylation kinetics of the enzyme using a series of unmethylated, hemimethylated and doubly methylated miRNA and siRNA substrates. The present study shows that HEN1 specifically binds double-stranded unmethylated or hemimethylated miR173/miR173* substrates with a subnanomolar affinity in a cofactor-dependent manner. Kinetic studies under single turnover and pre-steady state conditions in combination with isotope partitioning analysis showed that the binary HEN1-miRNA/miRNA* complex is catalytically competent; however, successive methylation of the two strands in a RNA duplex occurs in a non-processive (distributive) manner. We also find that the observed moderate methylation strand preference is largely exerted at the RNA-binding step and is fairly independent of the nature of the 3'-terminal nucleobase, but shows some dependency on proximal nucleotide mispairs. The results of the present study thus provide novel insights into the mechanism of RNA recognition and modification by a representative small RNA 2'-O-methyltransferase.

Molecular cloning and characterisation of SlAGO family in tomato

BACKGROUND

AGO (Argonaute) protein participates in plant developmental processes and virus defense as a core element of transcriptional regulator or/and post-transcriptional regulator in RNA induced silencing complex (RISC), which is guided by small RNAs to repress target genes expression. Previously, it was revealed that 15 putative AGO genes in tomato genome.

RESULTS

In present study, out of 15 detected SlAGO genes, only SlAGO4C and SlAGO15 couldn't be detected in roots, stems, leaves, buds, flowers and fruit of tomato by 30 cycles of PCR. SlAGO7 could be detected in early stage of fruit (-2 dpa, 0 dpa and 4 dpa), but it was significantly down-regulated in fruit collected on the 6 days post anthesis. Moreover, SlAGO5 could only be detected in reproductive tissues and SlAGO4D was specifically detected in fruit. According to blast result with miRNA database, three SlAGO genes harbored complementary sequences to miR168 (SlAGO1A and SlAGO1B) or miR403 (SlAGO2A). 5' RACE (Rapid amplification of cDNA ends) mapping was used to detect the 3' cleavage products of SlAGO mRNAs. In addition, subcellular localization of SlAGO proteins was detected. Our results showed that most SlAGO proteins localized to nucleus and cytoplasm. Importantly, nuclear membrane localization of AGO proteins was observed. Furthermore, mutated miR168 complementary site of SlAGO1A resulted in expanded localization of SlAGO1A, indicating that miR168 regulated localization of SlAGO1A.

CONCLUSIONS

Our results contribute to demonstration of potential roles of these newly isolated AGO family in tomato developmental processes and proved the conserved relationships between AGO genes and miRNAs in tomato, which might play important roles in tomato development and virus defense.

Interaction of small RNA-8105 and the intron of IbMYB1 RNA regulates IbMYB1 family genes through secondary siRNAs and DNA methylation after wounding

Small RNAs (sRNAs) play important roles in plants under stress conditions. However, limited research has been performed on the sRNAs involved in plant wound responses. In the present study, a novel wounding-induced sRNA, sRNA8105, was identified in sweet potato (Ipomoea batatas cv. Tainung 57) using microarray analysis. It was found that expression of sRNA8105 increased after mechanical wounding. Furthermore, Dicer-like 1 (DCL1) is required for the sRNA8105 precursor (pre-sRNA8105) to generate 22 and 24 nt mature sRNA8105. sRNA8105 targeted the first intron of IbMYB1 (MYB domain protein 1) before RNA splicing, and mediated RNA cleavage and DNA methylation of IbMYB1. The interaction between sRNA8105 and IbMYB1 was confirmed by cleavage site mapping, agro-infiltration analyses, and use of a transgenic sweet potato over-expressing pre-sRNA8105 gene. Induction of IbMYB1-siRNA was observed in the wild-type upon wounding and in transgenic sweet potato over-expressing pre-sRNA8105 gene without wounding, resulting in decreased expression of the whole IbMYB1 gene family, i.e. IbMYB1 and the IbMYB2 genes, and thus directing metabolic flux toward biosynthesis of lignin in the phenylpropanoid pathway. In conclusion, sRNA8105 induced by wounding binds to the first intron of IbMYB1 RNA to methylate IbMYB1, cleave IbMYB1 RNA, and trigger production of secondary siRNAs, further repressing the expression of the IbMYB1 family genes and regulating the phenylpropanoid pathway.

RNA silencing and its suppression: novel insights from in planta analyses

Plants employ multiple layers of innate immunity to fight pathogens. For both RNA and DNA viruses, RNA silencing plays a critical role in plant resistance. To escape this antiviral silencing-based immune response, viruses have evolved various counterdefense strategies, the most widespread being production of viral suppressors of RNA silencing (VSRs) that target various stages of the silencing mechanisms. Recent findings from in planta analyses have provided new insights into the mode of action of VSRs and revealed that plants react to the perturbation of the silencing pathways brought by viral infection by deploying a battery of counter-counterdefense measures. As well as discussing which experimental approaches have been most effective in delivering clear and unambiguous results, this review provides a detailed account of the surprising variety of offensive and defensive strategies set forth by both viruses and hosts in their struggle for survival.

Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks

Plant genomes are the source of large numbers of small RNAs, generated via a variety of genetically separable pathways. Several of these pathways converge in the production of phased, secondary, small interfering RNAs (phasiRNAs), originally designated as trans-acting small interfering RNAs or tasiRNAs. PhasiRNA biogenesis requires the involvement of microRNAs as well as the cellular machinery for the production of siRNAs. PhasiRNAs in Arabidopsis thaliana have been well described for their ability to function in trans to suppress target transcript levels. Plant genomic data from an expanding set of species have demonstrated that Arabidopsis is relatively sparing in its use of phasiRNAs, while other genomes contain hundreds or even thousands of phasiRNA-generating loci. In the dicots, targets of those phasiRNAs include several large or conserved families of genes, such as those encoding disease resistance proteins or transcription factors. Suppression of nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes by small RNAs is particularly unusual because of a high level of redundancy. In this review, we discuss plant phasiRNAs and the possible mechanistic significance of phasiRNA-based regulation of the NB-LRRs.

Biogenesis, turnover, and mode of action of plant microRNAs

MicroRNAs (miRNAs) are small RNAs that control gene expression through silencing of target mRNAs. Mature miRNAs are processed from primary miRNA transcripts by the endonuclease activity of the DICER-LIKE1 (DCL1) protein complex. Mechanisms exist that allow the DCL1 complex to precisely excise the miRNA from its precursor. Our understanding of miRNA biogenesis, particularly its intersection with transcription and other aspects of RNA metabolism such as splicing, is still evolving. Mature miRNAs are incorporated into an ARGONAUTE (AGO) effector complex competent for target gene silencing but are also subjected to turnover through a degradation mechanism that is beginning to be understood. The mechanisms of miRNA target silencing in plants are no longer limited to AGO-catalyzed slicing, and the contribution of translational inhibition is increasingly appreciated. Here, we review the mechanisms underlying the biogenesis, turnover, and activities of plant miRNAs.

Plant microRNAs display differential 3' truncation and tailing modifications that are ARGONAUTE1 dependent and conserved across species

Plant small RNAs are 3' methylated by the methyltransferase HUA1 ENHANCER1 (HEN1). In plant hen1 mutants, 3' modifications of small RNAs, including oligo-uridylation (tailing), are associated with accelerated degradation of microRNAs (miRNAs). By sequencing small RNAs of the wild type and hen1 mutants from Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays), we found 3' truncation prior to tailing is widespread in these mutants. Moreover, the patterns of miRNA truncation and tailing differ substantially among miRNA families but are conserved across species. The same patterns are also observable in wild-type libraries from a broad range of species, only at lower abundances. ARGONAUTE (AGO1), even with defective slicer activity, can bind these truncated and tailed variants of miRNAs. An ago1 mutation in hen1 suppressed such 3' modifications, indicating that they occur while miRNAs are in association with AGO1, either during or after RNA-induced silencing complex assembly. Our results showed AGO1-bound miRNAs are actively 3' truncated and tailed, possibly reflecting the activity of cofactors acting in conserved patterns in miRNA degradation.

Comprehensive investigation of microRNAs enhanced by analysis of sequence variants, expression patterns, ARGONAUTE loading, and target cleavage

MicroRNAs (miRNAs) are a class of small RNAs that typically function by guiding the cleavage of target messenger RNAs. They have been shown to play major roles in a variety of plant processes, including development, and responses to pathogens and environmental stresses. To identify new miRNAs and regulation in Arabidopsis (Arabidopsis thaliana), 27 small RNA libraries were constructed and sequenced from various tissues, stresses, and small RNA biogenesis mutants, resulting in 95 million genome-matched sequences. The use of rdr2 to enrich the miRNA population greatly enhanced this analysis and led to the discovery of new miRNAs arising from both known and new precursors, increasing the total number of Arabidopsis miRNAs by about 10%. Parallel Analysis of RNA Ends data provide evidence that the majority guide target cleavage. Many libraries represented novel stress/tissue conditions, such as submergence-stressed flowers, which enabled the identification of new stress regulation of both miRNAs and their targets, all of which were validated in wild-type plants. By combining small RNA expression analysis with ARGONAUTE immunoprecipitation data and global target cleavage data from Parallel Analysis of RNA Ends, a much more complete picture of Arabidopsis miRNAs was obtained. In particular, the discovery of ARGONAUTE loading and target cleavage biases gave important insights into tissue-specific expression patterns, pathogen responses, and the role of sequence variation among closely related miRNA family members that would not be evident without this combinatorial approach.

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.

MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA Biogenesis in Eudicots

Trans-acting small interfering RNAs (tasiRNAs) are a major class of small RNAs performing essential biological functions in plants. The first reported tasiRNA pathway, that of miR173-TAS1/2, produces tasiRNAs regulating a set of pentatricopeptide repeat (PPR) genes and has been characterized only in Arabidopsis thaliana to date. Here, we demonstrate that the microRNA (miRNA)-trans-acting small interfering RNA gene (TAS)-pentatricopeptide repeat-containing gene (PPR)-small interfering RNA pathway is a highly dynamic and widespread feature of eudicots. Nine eudicot plants, representing six different plant families, have evolved similar tasiRNA pathways to initiate phased small interfering RNA (phasiRNA) production from PPR genes. The PPR phasiRNA production is triggered by different 22-nucleotide miRNAs, including miR7122, miR1509, and fve-PPRtri1/2, and through distinct mechanistic strategies exploiting miRNA direct targeting or indirect targeting through TAS-like genes (TASL), one-hit or two-hit, or even two layers of tasiRNA-TASL interactions. Intriguingly, although those miRNA triggers display high sequence divergence caused by the occurrence of frequent point mutations and splicing shifts, their corresponding MIRNA genes show pronounced identity to the Arabidopsis MIR173, implying a common origin of this group of miRNAs (super-miR7122). Further analyses reveal that super-miR7122 may have evolved from a newly defined miR4376 superfamily, which probably originated from the widely conserved miR390. The elucidation of this evolutionary path expands our understanding of the course of miRNA evolution, especially for relatively conserved miRNA families.