AGO-unbound cytosolic pool of mature miRNAs in plant cells reveals a novel regulatory step at AGO1 loading.

RNA interference (RNAi) is mediated by small, 20-24-nt-long, non-coding regulatory (s)RNAs such as micro (mi) and small interfering (si) RNAs via the action of ARGONAUTE (AGO) proteins. High-throughput sequencing of size-separated sRNA pools of plant crude extracts revealed that the majority of the canonical miRNAs were associated with high molecular weight RNA-induced silencing complexes co-migrating with AGO1 (HMW RISC). In contrast, the majority of 24-nt-long siRNAs were found in association with low molecular weight complexes co-migrating with AGO4 (LMW RISC). Intriguingly, we identified a large set of cytoplasmic sRNAs, including mature miRNA sequences, in the low molecular size range corresponding to protein-unbound sRNAs. By comparing the RISC-loaded and protein-unbound pools of miRNAs, we identified miRNAs with highly different loading efficiencies. Expression of selected miRNAs in transient and transgenic systems validated their altered loading abilities implying that this process is controlled by information associated with the diverse miRNA precursors. We also showed that the availability of AGO proteins is a limiting factor determining the loading efficiency of miRNAs. Our data reveal the existence of a regulatory checkpoint determining the RISC-loading efficiencies of various miRNAs by sorting only a subset of the produced miRNAs into the biologically active RISCs.

Suppression of NB-LRR genes by miRNAs promotes nitrogen-fixing nodule development in Medicago truncatula.

Plant genomes contain two major classes of innate immune receptors to recognize different pathogens. The pattern recognition receptors perceive conserved pathogen-associated molecular patterns and the resistance genes with nucleotide-binding (NB) and leucine-rich repeat (LRR) domains recognize specific pathogen effectors. The precise regulation of resistance genes is important since the unregulated expression of NB-LRR genes can inhibit growth and may result in autoimmunity in the absence of pathogen infection. It was shown that a subset of miRNAs could target NB-LRR genes and act as an important regulator of plant immunity in the absence of pathogens. Plants not only interact with pathogens, but they can also establish symbiotic interactions with microbes. Nitrogen-fixing symbiotic interaction and nodule formation of legumes may also require the suppression of host defence to prevent immune responses. We found that upon symbiotic interactions, miRNAs repressing NB-LRR expression are upregulated in the developing nodules of Medicago truncatula. Furthermore, we show that the suppression of the activity of the NB-LRR genes targeted by these miRNAs is important during nodule development. Our results suggest that the downregulation of NB-LRR resistance genes in the developing nodule produces a suitable niche that facilitates bacterial colonization and the development of an N-fixing nodule.

Ambient temperature regulates the expression of a small set of sRNAs influencing plant development through NF-YA2 and YUC2.

Plants substantially alter their developmental programme upon changes in the ambient temperature. The 21-24 nt small RNAs (sRNAs) are important gene expression regulators, which play a major role in development and adaptation. However, little is known about how the different sRNA classes respond to changes in the ambient temperature. We profiled the sRNA populations in four different tissues of Arabidopsis thaliana plants grown at 15°C, 21°C, and 27°C. We found that only a small fraction (0.6%) of the sRNA loci are ambient temperature-controlled. We identified thermoresponsive microRNAs and identified their target genes using degradome libraries. We verified that the target of the thermoregulated miR169, NF-YA2, is also ambient temperature-regulated. NF-YA2, as the component of the conserved transcriptional regulator NF-Y complex, binds the promoter of the flowering time regulator FT and the auxin biosynthesis gene YUC2. Other differentially expressed loci include thermoresponsive phased siRNA loci that target various auxin pathway genes and tRNA fragments. Furthermore, a temperature-dependent 24-nt heterochromatic siRNA locus in the promoter of YUC2 may contribute to the epigenetic regulation of auxin homeostasis. This holistic approach facilitated a better understanding of the role of different sRNA classes in ambient temperature adaptation of plants.

Molecular characterization of a beet ringspot nepovirus isolated from Begonia ricinifolia in Hungary.

A nepovirus was isolated from Begonia ricinifolia showing chlorotic ringspot and line pattern symptoms. The purified virus had spherical particles of ca. 30 nm and contained a single coat protein subunit of ca. 56 kDa. The complete nucleotide sequence of the bipartite viral genome was determined. RNA 1 is 7394 nucleotides long, flanked by 5' and 3' untranslated regions (UTR), and followed by a 3' poly-A tail. It contains a single 6810 nt long open reading frame (ORF), which is translated into a 255 kDa polyprotein composed of 2269 amino acids. The 4684 nt long RNA 2 has a 4053 nt long ORF which encodes a single polyprotein of 1350 amino acids with a molecular weight of 149 kDa. Sequence comparisons revealed that the virus isolated from B. ricinifolia has the highest sequence similarity to beet ringspot virus and should be considered as a strain of BRSV. This is the first report on the occurrence of BRSV in B. ricinifolia and the presence of this virus outside Scotland.

Identification of Nicotiana benthamiana microRNAs and their targets using high throughput sequencing and degradome analysis.

BACKGROUND: Nicotiana benthamiana is a widely used model plant species for research on plant-pathogen interactions as well as other areas of plant science. It can be easily transformed or agroinfiltrated, therefore it is commonly used in studies requiring protein localization, interaction, or plant-based systems for protein expression and purification. To discover and characterize the miRNAs and their cleaved target mRNAs in N. benthamiana, we sequenced small RNA transcriptomes and degradomes of two N. benthamiana accessions and validated them by Northern blots. RESULTS: We used a comprehensive molecular approach to detect and to experimentally validate N. benthamiana miRNAs and their target mRNAs from various tissues. We identified 40 conserved miRNA families and 18 novel microRNA candidates and validated their target mRNAs with a genomic scale approach. The accumulation of thirteen novel miRNAs was confirmed by Northern blot analysis. The conserved and novel miRNA targets were found to be involved in various biological processes including transcription, RNA binding, DNA modification, signal transduction, stress response and metabolic process. Among the novel miRNA targets we found the mRNA of REPRESSOR OF SILENCING (ROS1). Regulation of ROS1 by a miRNA provides a new regulatory layer to reinforce transcriptional gene silencing by a post-transcriptional repression of ROS1 activity. CONCLUSIONS: The identified conserved and novel miRNAs along with their target mRNAs also provides a tissue specific atlas of known and new miRNA expression and their cleaved target mRNAs of N. benthamiana. Thus this study will serve as a valuable resource to the plant research community that will be beneficial well into the future.

RNA interference-mediated intrinsic antiviral immunity in plants.

RNA interference (RNAi) is an evolutionarily conserved, sequence-specific gene-inactivation system that plays an essential role in many biological processes, such as genome defense against mobile DNA elements or regulation of factors involved in plant and animal development. In higher plants and invertebrates, it also functions as a powerful antiviral mechanism. To overcome antiviral RNAi, viruses have evolved suppressor proteins which counteract host RNAi-based antiviral processes and target one or more key points in the RNAi machinery. Here, we review recent progress in our understanding of the mechanism and function of antiviral RNAi in plants and on the viral responses through the expression of silencing suppressor proteins. As a counter-attack RNAi may also regulate innate immunity in plants and contribute to a novel layer of defense against pathogen attack. We also discuss emerging evidence that viruses use RNAi to manipulate host gene expression to modify the cellular environment for the benefit of invading viruses.

PAREsnip: a tool for rapid genome-wide discovery of small RNA/target interactions evidenced through degradome sequencing.

Small RNAs (sRNAs) are a class of short (20-25 nt) non-coding RNAs that play important regulatory roles in gene expression. An essential first step in understanding their function is to confidently identify sRNA targets. In plants, several classes of sRNAs such as microRNAs (miRNAs) and trans-acting small interfering RNAs have been shown to bind with near-perfect complementarity to their messenger RNA (mRNA) targets, generally leading to cleavage of the mRNA. Recently, a high-throughput technique known as Parallel Analysis of RNA Ends (PARE) has made it possible to sequence mRNA cleavage products on a large-scale. Computational methods now exist to use these data to find targets of conserved and newly identified miRNAs. Due to speed limitations such methods rely on the user knowing which sRNA sequences are likely to target a transcript. By limiting the search to a tiny subset of sRNAs it is likely that many other sRNA/mRNA interactions will be missed. Here, we describe a new software tool called PAREsnip that allows users to search for potential targets of all sRNAs obtained from high-throughput sequencing experiments. By searching for targets of a complete 'sRNAome' we can facilitate large-scale identification of sRNA targets, allowing us to discover regulatory interaction networks.

Distribution of Small RNAs Along Transposable Elements in Vitis vinifera During Somatic Embryogenesis.

Metaviridae is a family of reverse-transcribing viruses, closely related to retroviruses; they exist within their host's DNA as transposable elements. Transposable element study requires the use of specialized tools, in part because of their repetitive nature. By combining data from transcript RNA-Seq, small RNA-Seq, and parallel analysis of RNA ends-Seq from grapevine somatic embryos, we set up a bioinformatics flowchart that could be able to assemble and identify transposable elements.

Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns.

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.

Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions.

MicroRNAs negatively regulate the accumulation of mRNAs therefore when they are expressed in the same cells their expression profiles show an inverse correlation. We previously described one positively correlated miRNA/target pair, but it is not known how widespread this phenomenon is. Here, we investigated the correlation between the expression profiles of differentially expressed miRNAs and their targets during tomato fruit development using deep sequencing, Northern blot and RT-qPCR. We found an equal number of positively and negatively correlated miRNA/target pairs indicating that positive correlation is more frequent than previously thought. We also found that the correlation between microRNA and target expression profiles can vary between mRNAs belonging to the same gene family and even for the same target mRNA at different developmental stages. Since microRNAs always negatively regulate their targets, the high number of positively correlated microRNA/target pairs suggests that mutual exclusion could be as widespread as temporal regulation. The change of correlation during development suggests that the type of regulatory circuit directed by a microRNA can change over time and can be different for individual gene family members. Our results also highlight potential problems for expression profiling-based microRNA target identification/validation.

Structural and functional analysis of viral siRNAs.

A large amount of short interfering RNA (vsiRNA) is generated from plant viruses during infection, but the function, structure and biogenesis of these is not understood. We profiled vsiRNAs using two different high-throughput sequencing platforms and also developed a hybridisation based array approach. The profiles obtained through the Solexa platform and by hybridisation were very similar to each other but different from the 454 profile. Both deep sequencing techniques revealed a strong bias in vsiRNAs for the positive strand of the virus and identified regions on the viral genome that produced vsiRNA in much higher abundance than other regions. The hybridisation approach also showed that the position of highly abundant vsiRNAs was the same in different plant species and in the absence of RDR6. We used the Terminator 5'-Phosphate-Dependent Exonuclease to study the 5' end of vsiRNAs and showed that a perfect control duplex was not digested by the enzyme without denaturation and that the efficiency of the Terminator was strongly affected by the concentration of the substrate. We found that most vsiRNAs have 5' monophosphates, which was also confirmed by profiling short RNA libraries following either direct ligation of adapters to the 5' end of short RNAs or after replacing any potential 5' ends with monophosphates. The Terminator experiments also showed that vsiRNAs were not perfect duplexes. Using a sensor construct we also found that regions from the viral genome that were complementary to non-abundant vsiRNAs were targeted in planta just as efficiently as regions recognised by abundant vsiRNAs. Different high-throughput sequencing techniques have different reproducible sequence bias and generate different profiles of short RNAs. The Terminator exonuclease does not process double stranded RNA, and because short RNAs can quickly re-anneal at high concentration, this assay can be misleading if the substrate is not denatured and not analysed in a dilution series. The sequence profiles and Terminator digests suggest that CymRSV siRNAs are produced from the structured positive strand rather than from perfect double stranded RNA or by RNA dependent RNA polymerase.

Ecotype-specific blockage of tasiARF production by two different RNA viruses in Arabidopsis.

Arabidopsis thaliana is one of the most studied model organisms of plant biology with hundreds of geographical variants called ecotypes. One might expect that this enormous genetic variety could result in differential response to pathogens. Indeed, we observed previously that the Bur ecotype develops much more severe symptoms (upward curling leaves and wavy leaf margins) upon infection with two positive-strand RNA viruses of different families (turnip vein-clearing virus, TVCV, and turnip mosaic virus, TuMV). To find the genes potentially responsible for the ecotype-specific response, we performed a differential expression analysis of the mRNA and sRNA pools of TVCV and TuMV-infected Bur and Col plants along with the corresponding mock controls. We focused on the genes and sRNAs that showed an induced or reduced expression selectively in the Bur virus samples in both virus series. We found that the two ecotypes respond to the viral infection differently, yet both viruses selectively block the production of the TAS3-derived small RNA specimen called tasiARF only in the virus-infected Bur plants. The tasiARF normally forms a gradient through the adaxial and abaxial parts of the leaf (being more abundant in the adaxial part) and post-transcriptionally regulates ARF4, a major leaf polarity determinant in plants. The lack of tasiARF-mediated silencing could lead to an ectopically expressed ARF4 in the adaxial part of the leaf where the misregulation of auxin-dependent signaling would result in an irregular growth of the leaf blade manifesting as upward curling leaf and wavy leaf margin. QTL mapping using Recombinant Inbred Lines (RILs) suggests that the observed symptoms are the result of a multigenic interaction that allows the symptoms to develop only in the Bur ecotype. The particular nature of genetic differences leading to the ecotype-specific symptoms remains obscure and needs further study.

Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis.

In plants, microRNAs (miRNAs) comprise one of three classes of small RNAs regulating gene expression at the post-transcriptional level. Many plant miRNAs are conserved, and play a role in development, abiotic stress responses or pathogen responses. However, some miRNAs have only been found in certain species. Here, we use deep-sequencing, computational and molecular methods to identify, profile, and describe conserved and non-conserved miRNAs in four grapevine (Vitis vinifera) tissues. A total of 24 conserved miRNA families were identified in all four tissues, and 26 known but non-conserved miRNAs were also found. In addition to known miRNAs, we also found 21 new grapevine-specific miRNAs together with their star strands. We have also shown that almost all of them originated from single genes. Furthermore, 21 other plausible miRNA candidates have been described. We have found that many known and new miRNAs showed tissue-specific expression. Finally, 112 target mRNAs of known and 44 target mRNAs of new grapevine-specific miRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs.

Genome-Wide Identification of RNA Silencing-Related Genes and Their Expressional Analysis in Response to Heat Stress in Barley (Hordeum vulgare L.).

Barley (Hordeum vulgare L.) is an economically important crop cultivated in temperate climates all over the world. Adverse environmental factors negatively affect its survival and productivity. RNA silencing is a conserved pathway involved in the regulation of growth, development and stress responses. The key components of RNA silencing are the Dicer-like proteins (DCLs), Argonautes (AGOs) and RNA-dependent RNA polymerases (RDRs). Despite its economic importance, there is no available comprehensive report on barley RNA silencing machinery and its regulation. In this study, we in silico identified five DCL (HvDCL), eleven AGO (HvAGO) and seven RDR (HvRDR) genes in the barley genome. Genomic localization, phylogenetic analysis, domain organization and functional/catalytic motif identification were also performed. To understand the regulation of RNA silencing, we experimentally analysed the transcriptional changes in response to moderate, persistent or gradient heat stress treatments: transcriptional accumulation of siRNA- but not miRNA-based silencing factor was consistently detected. These results suggest that RNA silencing is dynamically regulated and may be involved in the coordination of development and environmental adaptation in barley. In summary, our work provides information about barley RNA silencing components and will be a ground for the selection of candidate factors and in-depth functional/mechanistic analyses.

Transcriptome reprogramming in the shoot apical meristem of CymRSV-infected Nicotiana benthamiana plants associates with viral exclusion and the lack of recovery.

In some plant-virus interactions plants show a sign of healing from virus infection, a phenomenon called symptom recovery. It is assumed that the meristem exclusion of the virus is essential to this process. The discovery of RNA silencing provided a possible mechanism to explain meristem exclusion and recovery. Here we show evidence that silencing is not the reason for meristem exclusion in Nicotiana benthamiana plants infected with Cymbidium ringspot virus (CymRSV). Transcriptome analysis followed by in situ hybridization shed light on the changes in gene expression in the shoot apical meristem (SAM) on virus infection. We observed the down-regulation of meristem-specific genes, including WUSCHEL (WUS). However, WUS was not down-regulated in the SAM of plants infected with meristem-invading viruses such as turnip vein-clearing virus (TVCV) and cucumber mosaic virus (CMV). Moreover, there is no connection between loss of meristem function and fast shoot necrosis since TVCV necrotized the shoot while CMV did not. Our findings suggest that the observed transcriptional changes on virus infection in the shoot are key factors in tip necrosis and symptom recovery. We observed a lack of GLYCERALDEHYDE 3-PHOSPHATE DEHYDROGENASE (GAPDH) expression in tissues around the meristem, which likely stops virus replication and spread into the meristem.

Molecular characterization and In Vitro synthesis of infectious RNA of a Turnip vein-clearing virus isolated from Alliaria petiolata in Hungary.

A tobamovirus was isolated from leaves of Alliaria petiolata plants, showing vein-clearing, interveinal chlorosis, and moderate deformation. Host range experiments revealed a high similarity of isolate ApH both to ribgrass mosaic viruses and turnip vein-clearing viruses. The complete nucleotide sequence of the viral genome was determined. The genomic RNA is composed of 6312 nucleotides and contains four open reading frames (ORF). ORF1 is 3324 nt-long and encodes a polypeptide of about 125.3 kDa. The ORF1 encoded putative replication protein contains an Alphavirus-like methyltransferase domain. ORF2 is 4806 nt-long and encodes a polypeptide of about 182 kDa. The ORF2 encoded putative replication protein contains an RNA-dependent RNA polymerase, catalytic domain. ORF3 encodes the putative cell-to-cell movement protein with a molecular weight of 30.1 kDa. ORF4 overlaps with ORF3 and encodes the coat protein with a size of 17.5 kDa. Sequence comparisons revealed that the ApH isolate has the highest similarity to turnip vein-clearing viruses and should be considered an isolate of Turnip vein-clearing virus (TVCV). This is the first report on the occurrence of TVCV in Hungary. In vitro transcripts prepared from the full-length cDNA clone of TVCV-ApH were highly infectious and induced typical symptoms characteristic to the original isolate of the virus. Since infectious clones of TVCV-ApH and crTMV (another isolate of TVCV) markedly differed in respect to recovery phenotype in Arabidopsis thaliana, it is feasible to carry out gene exchange or mutational studies to determine viral factors responsible for the symptom recovery phenotype.

High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families.

BACKGROUND: High-throughput sequencing technology is capable to identify novel short RNAs in plant species. We used Solexa sequencing to find new microRNAs in one of the model legume species, barrel medic (Medicago truncatula). RESULTS: 3,948,871 reads were obtained from two separate short RNA libraries generated from total RNA extracted from M. truncatula leaves, representing 1,563,959 distinct sequences. 2,168,937 reads were mapped to the available M. truncatula genome corresponding to 619,175 distinct sequences. 174,504 reads representing 25 conserved miRNA families showed perfect matches to known miRNAs. We also identified 26 novel miRNA candidates that were potentially generated from 32 loci. Nine of these loci produced eight distinct sequences, for which the miRNA* sequences were also sequenced. These sequences were not described in other plant species and accumulation of these eight novel miRNAs was confirmed by Northern blot analysis. Potential target genes were predicted for most conserved and novel miRNAs. CONCLUSION: Deep sequencing of short RNAs from M. truncatula leaves identified eight new miRNAs indicating that specific miRNAs exist in legume species.

Expansion of Capsicum annuum fruit is linked to dynamic tissue-specific differential expression of miRNA and siRNA profiles.

Small regulatory RNAs, such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) have emerged as important transcriptional and post-transcriptional regulators controlling a wide variety of physiological processes including fruit development. Data are, however, limited for their potential roles in developmental processes determining economically important traits of crops. The current study aimed to discover and characterize differentially expressed miRNAs and siRNAs in sweet pepper (Capsicum annuum) during fruit expansion. High-throughput sequencing was employed to determine the small regulatory RNA expression profiles in various fruit tissues, such as placenta, seed, and flesh at 28 and 40 days after anthesis. Comparative differential expression analyses of conserved, already described and our newly predicted pepper-specific miRNAs revealed that fruit expansion is accompanied by an increasing level of miRNA-mediated regulation of gene expression. Accordingly, ARGONAUTE1 protein, the primary executor of miRNA-mediated regulation, continuously accumulated to an extremely high level in the flesh. We also identified numerous pepper-specific, heterochromatin-associated 24-nt siRNAs (hetsiRNAs) which were extremely abundant in the seeds, as well as 21-nt and 24-nt phased siRNAs (phasiRNAs) that were expressed mainly in the placenta and the seeds. This work provides comprehensive tissue-specific miRNA and siRNA expression landscape for a developing pepper fruit. We identified several novel, abundantly expressing tissue- and pepper-specific small regulatory RNA species. Our data show that fruit expansion is associated with extensive changes in sRNA abundance, raising the possibility that manipulation of sRNA pathways may be employed to improve the quality and quantity of the pepper fruit.

Correction: Expansion of Capsicum annum fruit is linked to dynamic tissue-specific differential expression of miRNA and siRNA profiles.

[This corrects the article DOI: 10.1371/journal.pone.0200207.].

Identification of ARGONAUTE/Small RNA Cleavage Sites by Degradome Sequencing.

The method described here enables the high-throughput identification of cleaved mRNA targets of ARGONAUTE/small RNA complexes. The protocol is based on a modified 5'-rapid amplification of cDNA ends combined with deep sequencing of 3' cleavage products of mRNAs. Poly(A) RNA is purified from the tissue of interest which is followed by a 5'-RNA adapter ligation. The ligated products are then reverse transcribed, amplified, and digested with MmeI. After gel separation, a 3' double-stranded DNA adapter is ligated to the fragments, which are then amplified and index labeled for the high-throughput sequencing of pooled degradome libraries. Sequencing datasets from pooled libraries can be analyzed with different bioinformatic approaches.