Spatio-temporal accumulation of microRNAs is highly coordinated in developing plant tissues

FvemiR160-FveARF18A-FveAP1/FveFUL module regulates flowering time in woodland strawberry

Flowering is an indicator of plant transformation from vegetative to reproductive growth. miR160 has been shown to have a significant effect on the growth and development of fruits, leaves, and roots of plants or their stress response to environment, but the participation of miR160 in regulating flowering time in plants is unclear. In this study, we found that two FvemiR160s (FvemiR160a/FvemiR160b) mature sequences in strawberry (Fragaria vesca) were consistent. It was displayed that the miR160 mature sequence is highly conserved in various species, and the miR160 mature sequence formed by the 5' arm of the MIR160 precursor was more conserved. Three FveARFs in woodland strawberry were negatively regulated by FvemiR160a, among which FveARF18A was the most significant. Phylogenetic analysis indicated that FvemiR160 is closely related to apple (Malus domestica), grape (Vitis vinifera), and Arabidopsis thaliana, while FveARF18A is closely related to RcARF18. Subsequently, we demonstrated that FvemiR160a can target cutting FveARF18A to negatively regulate its expression by RLM-5' RACE, cleavage site mutation, and GFP fluorescence assay. Moreover, we observed that FveMIR160a overexpressed plants have advanced flowering, while mFveARF18A overexpressed plants have delayed flowering. We also verified that FveARF18A negatively regulates the expression of FveAP1 and FveFUL by binding their promoters by yeast one-hybrid, LUC, and GUS assay, and FveAP1 and FveFUL transgenic Arabidopsis showed early flowering phenotype. In addition, the expression level of FvemiR160a was decreased obviously while that of FveARF18A was increased obviously by MeJA, GA and IAA. In conclusion, our study reveals the important role of the FvemiR160-FveARF18A-FveAP1/FveFUL module in the flowering process of woodland strawberry and provides a new pathway for studying flowering.

miR160: An Indispensable Regulator in Plant

MicroRNAs (miRNA), recognized as crucial regulators of gene expression at the posttranscriptional level, have been found to be involved in the biological processes of plants. Some miRNAs are up- or down-regulated during plant development, stress response, and secondary metabolism. Over the past few years, it has been proved that miR160 is directly related to the developments of different tissues and organs in multifarious species, as well as plant-environment interactions. This review highlights the recent progress on the contributions of the miR160-ARF module to important traits of plants and the role of miR160-centered gene regulatory network in coordinating growth with endogenous and environmental factors. The manipulation of miR160-guided gene regulation may provide a new method to engineer plants with improved adaptability and yield.

MicroRNA-Mediated Responses to Cadmium Stress in Arabidopsis thaliana

In recent decades, the presence of cadmium (Cd) in the environment has increased significantly due to anthropogenic activities. Cd is taken up from the soil by plant roots for its subsequent translocation to shoots. However, Cd is a non-essential heavy metal and is therefore toxic to plants when it over-accumulates. MicroRNA (miRNA)-directed gene expression regulation is central to the response of a plant to Cd stress. Here, we document the miRNA-directed response of wild-type Arabidopsis thaliana (Arabidopsis) plants and the drb1, drb2 and drb4 mutant lines to Cd stress. Phenotypic and physiological analyses revealed the drb1 mutant to display the highest degree of tolerance to the imposed stress while the drb2 mutant was the most sensitive. RT-qPCR-based molecular profiling of miRNA abundance and miRNA target gene expression revealed DRB1 to be the primary double-stranded RNA binding (DRB) protein required for the production of six of the seven Cd-responsive miRNAs analyzed. However, DRB2, and not DRB1, was determined to be required for miR396 production. RT-qPCR further inferred that transcript cleavage was the RNA silencing mechanism directed by each assessed miRNA to control miRNA target gene expression. Taken together, the results presented here reveal the complexity of the miRNA-directed molecular response of Arabidopsis to Cd stress.

In Situ Detection of Mature miRNAs in Plants Using LNA-Modified DNA Probes

MicroRNAs (miRNAs) play important roles in development in plants, and some miRNAs show developmentally regulated organ- and tissue-specific expression patterns. Therefore, in situ detection of mature miRNAs is important for understanding the functions for both miRNAs and their targets. The construction of promoter-reporter fusions and examination of their in planta expression has been widely used and the results obtained thus far are rather informative; however, in some cases, the length of promoter that contains entire regulatory elements is difficult to determine. In addition, traditional in situ hybridization with the antisense RNA fragment as the probe usually fails to detect miRNAs, because the mature miRNAs are too short (~21-nucleotides) to exhibit stable hybridization signals. In recent years, the Locked nucleic acid (LNA) modified DNA probe has been successfully used in animals and plants to detect small RNAs. Here, we describe a modified protocol using LNA-modified DNA probes to detect mature miRNAs in plant tissues, including the design of LNA probes and detailed steps for the in situ hybridization experiment, using Arabidopsis miR165 as an example.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

Small RNA In Situ Hybridizations on Sections of Arabidopsis Embryos

Small RNAs mediate posttranscriptional gene silencing in plants and animals. This often occurs in specific cell or tissue types and can be necessary for their differentiation. Determining small RNA (sRNA) localization patterns at cellular resolution can therefore provide information on the corresponding gene regulatory processes they are involved in. Recent improvements with in situ hybridization methods have allowed them to be applied to sRNAs. Here we describe an in situ hybridization protocol to detect sRNAs from sections of early staged Arabidopsis thaliana (Arabidopsis) embryos.

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.

miR167 limits anther growth to potentiate anther dehiscence

In flowering plants, anther dehiscence and pollen release are essential for sexual reproduction. Anthers dehisce after cell wall degradation weakens stomium cell junctions in each anther locule, and desiccation creates mechanical forces that open the locules. Either effect or both together may break stomium cell junctions. The microRNA miR167 negatively regulates ARF6 and ARF8, which encode auxin response transcription factors. Arabidopsis mARF6 or mARF8 plants with mutated miR167 target sites have defective anther dehiscence and ovule development. Null mir167a mutations recapitulated mARF6 and mARF8 anther and ovule phenotypes, indicating that MIR167a is the main miR167 precursor gene that delimits ARF6 and ARF8 expression in these organs. Anthers of mir167a or mARF6/8 plants overexpressed genes encoding cell wall loosening functions associated with cell expansion, and grew larger than wild-type anthers did starting at flower stage 11. Experimental desiccation enabled dehiscence of miR167-deficient anthers, indicating competence to dehisce. Conversely, high humidity conditions delayed anther dehiscence in wild-type flowers. These results support a model in which miR167-mediated anther growth arrest permits anther dehiscence. Without miR167 regulation, excess anther growth delays dehiscence by prolonging desiccation.

Identification of novel miRNAs and their target genes from Populus szechuanica infected with Melampsora larici-populina

Two novel miRNAs were selected from a pre-constructed RNA library of Populus szechuanica infected with the foliar rust fungus Melampsora larici-populina in order to detect the genes regulated as targets of the miRNAs novel_mir_11 and novel_mir_357. The novel miRNAs were identified from P. szechuanica using stem-loop methods and their precursors were able to fold into a complete stem loop structure. The predicted target genes of the novel miRNAs were verified with RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). The full-length sequences of target genes, RPM1 and RPS2/5, in P. szechuanica were obtained through rapid amplification of cDNA ends (RACE) and officially named PsRPM1 and PsRPS2/5. These genes contain nucleotide binding site-leucine-rich repeats (NBS-LRR) domains typical of resistance genes. The expression levels of miRNAs and their target genes in different periods post infection were analysed with quantitative real-time PCR (qRT-PCR). After infection with the foliar rust fungus, the expression levels of the novel miRNAs and their target genes were dynamic. Both novel_mir_11 and novel_mir_357 negatively regulated the expression of their target genes. In this study, the regulatory effects of two novel miRNAs through their target genes were characterized to provide further mechanistic information regarding the interaction between Populus and a foliar rust fungus. Results of this study improve our understanding of the defence response mechanisms of Populus and will stimulate future work to characterize strategies to prevent and control Populus diseases.

Functional Characterization of microRNA171 Family in Tomato

Deeply conserved plant microRNAs (miRNAs) function as pivotal regulators of development. Nevertheless, in the model crop Solanum lycopersicum (tomato) several conserved miRNAs are still poorly annotated and knowledge about their functions is lacking. Here, the tomato miR171 family was functionally analyzed. We found that the tomato genome contains at least 11 SlMIR171 genes that are differentially expressed along tomato development. Downregulation of sly-miR171 in tomato was successfully achieved by transgenic expression of a short tandem target mimic construct (STTM171). Consequently, sly-miR171-targeted mRNAs were upregulated in the silenced plants. Target upregulation was associated with irregular compound leaf development and an increase in the number of axillary branches. A prominent phenotype of STTM171 expressing plants was their male sterility due to a production of a low number of malformed and nonviable pollen. We showed that sly-miR171 was expressed in anthers along microsporogenesis and significantly silenced upon STTM171 expression. Sly-miR171-silenced anthers showed delayed tapetum ontogenesis and reduced callose deposition around the tetrads, both of which together or separately can impair pollen development. Collectively, our results show that sly-miR171 is involved in the regulation of anther development as well as shoot branching and compound leaf morphogenesis.

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.

Sensitive whole mount in situ localization of small RNAs in plants

Small RNAs, such as microRNAs (miRNAs), regulate gene expression and play important roles in many plant processes. Although our knowledge of their biogenesis and mode of action has significantly progressed, we still have comparatively little information about their biological functions. In particular, knowledge about their spatio-temporal expression patterns rely on either indirect detection by use of reporter constructs or labor-intensive direct detection by in situ hybridization on sectioned material. None of the current approaches allows a systematic investigation of small RNA expression patterns. Here, we present a sensitive method for in situ detection of miRNAs and siRNAs in intact plant tissues that utilizes both double-labeled probes and a specific cross-linker. We determined the expression patterns of several small RNAs in diverse plant tissues.

Integrated mRNA and miRNA transcriptome reveal a cross-talk between developing response and hormone signaling for the seed kernels of Siberian apricot

Recently, our transcriptomic analysis has identified some functional genes responsible for oil biosynthesis in developing SASK, yet miRNA-mediated regulation for SASK development and oil accumulation is poorly understood. Here, 3 representative periods of 10, 30 and 60 DAF were selected for sRNA sequencing based on the dynamic patterns of growth tendency and oil content of developing SASK. By miRNA transcriptomic analysis, we characterized 296 known and 44 novel miRNAs in developing SASK, among which 36 known and 6 novel miRNAs respond specifically to developing SASK. Importantly, we performed an integrated analysis of mRNA and miRNA transcriptome as well as qRT-PCR detection to identify some key miRNAs and their targets (miR156-SPL, miR160-ARF18, miR164-NAC1, miR171h-SCL6, miR172-AP2, miR395-AUX22B, miR530-P2C37, miR393h-TIR1/AFB2 and psi-miRn5-SnRK2A) potentially involved in developing response and hormone signaling of SASK. Our results provide new insights into the important regulatory function of cross-talk between development response and hormone signaling for SASK oil accumulation.

RNA silencing movement in plants

Multicellular organisms, like higher plants, need to coordinate their growth and development and to cope with environmental cues. To achieve this, various signal molecules are transported between neighboring cells and distant organs to control the fate of the recipient cells and organs. RNA silencing produces cell non-autonomous signal molecules that can move over short or long distances leading to the sequence specific silencing of a target gene in a well defined area of cells or throughout the entire plant, respectively. The nature of these signal molecules, the route of silencing spread, and the genes involved in their production, movement and reception are discussed in this review. Additionally, a short section on features of silencing spread in animal models is presented at the end of this review.

RNA processing in auxin and cytokinin pathways

Auxin and cytokinin belong to the 'magnificent seven' plant hormones, having tightly interconnected pathways leading to common as well as opposing effects on plant morphogenesis. Tremendous progress in the past years has yielded a broad understanding of their signalling, metabolism, regulatory pathways, transcriptional networks, and signalling cross-talk. One of the rapidly expanding areas of auxin and cytokinin research concerns their RNA regulatory networks. This review summarizes current knowledge about post-transcriptional gene silencing, the role of non-coding RNAs, the regulation of translation, and alternative splicing of auxin- and cytokinin-related genes. In addition, the role of tRNA-bound cytokinins is also discussed. We highlight the most recent publications dealing with this topic and underline the role of RNA processing in auxin- and cytokinin-mediated growth and development.

Identification and characterization of microRNAs at different flowering developmental stages in moso bamboo (Phyllostachys edulis) by high-throughput sequencing

Researching moso bamboo flowering has been difficult because of its unknown flowering interval and the rarity of florescent samples. To identify microRNAs (miRNAs) and study their expression patterns during the flower developmental process of moso bamboo, small RNAs from non-flowering leaves and four flower developmental periods were sequenced using Illumina technology. In total, 409 known miRNAs and 492 differentially expressed novel miRNAs were identified in moso bamboo. Of the known miRNAs that were differentially expressed between non-flowering and flowering samples, 64 were predicted to have a total of 308 targets. Among the miRNAs, seven known and five novel miRNAs were selected, as were four of their target genes, and their expression profiles were validated using qRT-PCR. The results indicated that the miRNA expression levels were negatively correlated with those of their targets. The research comprehensively revealed that the differentially expressed miRNAs and their targets participated in diverse biological pathways and played significant regulatory roles in moso bamboo flowering. The data provide a significant resource for understanding the molecular mechanisms in moso bamboo flowering and senescence, and serve as the primary foundation for further studies on metabolic regulatory networks that involve miRNAs.

Small RNA deep sequencing reveals the important role of microRNAs in the halophyte Halostachys caspica

MicroRNAs (miRNAs), an extensive class of small regulatory RNAs, play versatile roles in plant growth and development as well as stress responses. However, the regulatory mechanism is unclear on miRNA-mediated response to abiotic stress in plants. Halostachys caspica is a halophytic plant species and a great model for investigating plant response to salinity stress. However, no research has been performed on miRNAs in H. caspica. In this study, we employed deep sequencing to identify both conserved and novel miRNAs from salinity-exposed H. caspica and its untreated control. Among the 13-19 million sequences generated from both treatments, a total of 170 conserved miRNAs, belonging to 151 miRNA families, were identified; among these miRNAs, 31 were significantly up-regulated and 48 were significantly down-regulated by salinity stress. We also identified 102 novel miRNAs from H. caspica; among them, 12 miRNAs were significantly up-regulated and 13 were significantly down-regulated by salinity. qRT-PCR expression analysis validated the deep sequencing results and also demonstrated that miRNAs and their targeted genes were responsive to high salt stress and existed a negative expression correlation between miRNAs and their targets. miRNA-target prediction, GO and KEGG analysis showed that miRNAs were involved in salt stress-related biological pathway, including calcium signalling pathway, MAPK signalling pathway, plant hormone signal transduction and flavonoid biosynthesis, etc. This suggests that miRNAs play an important role in plant salt stress tolerance in H. caspica. This result could be used to improve salt tolerance in crops and woods.

Identification of microRNAs and transcript targets in Camelina sativa by deep sequencing and computational methods

Camelina sativa is an annual oilseed crop that is under intensive development for renewable resources of biofuels and industrial oils. MicroRNAs, or miRNAs, are endogenously encoded small RNAs that play key roles in diverse plant biological processes. Here, we conducted deep sequencing on small RNA libraries prepared from camelina leaves, flower buds and two stages of developing seeds corresponding to initial and peak storage products accumulation. Computational analyses identified 207 known miRNAs belonging to 63 families, as well as 5 novel miRNAs. These miRNAs, especially members of the miRNA families, varied greatly in different tissues and developmental stages. The predicted miRNA target genes are involved in a broad range of physiological functions including lipid metabolism. This report is the first step toward elucidating roles of miRNAs in C. sativa and will provide additional tools to improve this oilseed crop for biofuels and biomaterials.

miRNA expression during prickly pear cactus fruit development

miRNAs are a class of small non-coding RNAs that regulate gene expression. They are involved in the control of many developmental processes, including fruit development. The increasing amount of information on miRNAs, on their expression, abundance, and conservation between various species, provides a new opportunity to study the role of miRNAs in non-model plant species. In this work, we used a combination of Northern blot and tissue print hybridization analysis to identify conserved miRNAs expressed during prickly pear cactus (Opuntia ficus indica) fruit development. Comparative profiling detected the expression of 34 miRNAs, which were clustered in three different groups that were associated with the different phases of fruit development. Variation in the level of miRNA expression was observed. Gradual expression increase of several miRNAs was observed during fruit development, including miR164. miR164 was selected for stem-loop RT-PCR and for a detailed spatial-temporal expression analysis. At early floral stages, miR164 was mainly localized in meristematic tissues, boundaries and fusion zones, while it was more homogenously expressed in fruit tissues. Our results provide the first evidence of miRNA expression in the prickly pear cactus and provide the basis for future research on miRNAs in Opuntia. Moreover, our analyses suggest that miR164 plays different roles during prickly pear cactus fruit development.

microRNA biogenesis, degradation and activity in plants

microRNAs (miRNAs) are important regulators of gene expression. After excised from primary miRNA transcript by dicer-like1 (DCL1, an RNAse III enzyme), miRNAs bind and guide their effector protein named argonaute 1 (AGO1) to silence the expression of target RNAs containing their complementary sequences in plants. miRNA levels and activities are tightly controlled to ensure their functions in various biological processes such as development, metabolism and responses to abiotic and biotic stresses. Studies have identified many factors that involve in miRNA accumulation and activities. Characterization of these factors in turn greatly improves our understanding of the processes related to miRNAs. Here, we review recent progress of mechanisms underlying miRNA expression and functions in plants.

Genome-wide analysis of alternative splicing in Zea mays: landscape and genetic regulation

Alternative splicing enhances transcriptome diversity in all eukaryotes and plays a role in plant tissue identity and stress adaptation. To catalog new maize (Zea mays) transcripts and identify genomic loci that regulate alternative splicing, we analyzed over 90 RNA-seq libraries from maize inbred lines B73 and Mo17, as well as Syn10 doubled haploid lines (progenies from B73 × Mo17). Transcript discovery was augmented with publicly available data from 14 maize tissues, expanding the maize transcriptome by more than 30,000 and increasing the percentage of intron-containing genes that undergo alternative splicing to 40%. These newly identified transcripts greatly increase the diversity of the maize proteome, sometimes coding for entirely different proteins compared with their most similar annotated isoform. In addition to increasing proteome diversity, many genes encoding novel transcripts gained an additional layer of regulation by microRNAs, often in a tissue-specific manner. We also demonstrate that the majority of genotype-specific alternative splicing can be genetically mapped, with cis-acting quantitative trait loci (QTLs) predominating. A large number of trans-acting QTLs were also apparent, with nearly half located in regions not shown to contain genes associated with splicing. Taken together, these results highlight the currently underappreciated role that alternative splicing plays in tissue identity and genotypic variation in maize.

Exploring new models of easiRNA biogenesis

Although silent transposons in plants can be reactivated by stress or during development, their potential deleterious effects are prevented by transposon-derived epigenetically activated small interfering RNAs (easiRNAs). A new study shows how serendipitous interactions between reactivated transposons and endogenous microRNAs might initiate easiRNA biogenesis, establishing an unexpected link between these two classes of silencing small RNAs.

A comprehensive genome-wide study on tissue-specific and abiotic stress-specific miRNAs in Triticum aestivum

Productivity of wheat crop is largely dependent on its growth and development that, in turn, is mainly regulated by environmental conditions, including abiotic stress factors. miRNAs are key regulators of gene expression networks involved in diverse aspects of development and stress responses in plants. Using high-throughput sequencing of eight small RNA libraries prepared from diverse abiotic stresses and tissues, we identified 47 known miRNAs belonging to 20 families, 49 true novel and 1030 candidate novel miRNAs. Digital gene expression analysis revealed that 257 miRNAs exhibited tissue-specific expression and 74 were associated with abiotic stresses. Putative target genes were predicted for miRNAs identified in this study and their grouping into functional categories indicated that the putative targets were involved in diverse biological processes. RLM-RACE of predicted targets of three known miRNAs (miR156, miR160 and miR164) confirmed their mRNA cleavage, thus indicating their regulation at post-transcriptional level by the corresponding miRNAs. Mapping of the sequenced data onto the wheat progenitors and closely related monocots revealed a large number of evolutionary conserved miRNAs. Additional expression profiling of some of these miRNAs in other abiotic stresses underline their involvement in multiple stresses. Our findings provide valuable resource for an improved understanding of the role of miRNAs in stress tolerance as well as plant development.

Computational identification of microRNAs and their targets in perennial Ryegrass (Lolium perenne)

MicroRNAs (miRNAs) are small non-coding RNA molecules of 22 nucleotides in length that have been characterized as regulators of messenger RNA (mRNA) regulating a number of developmental processes in plants and animals by silencing genes using multiple mechanisms. miRNAs have been extensively studied in various plant species; however, few information are available about miRNAs in perennial ryegrass, animal feed, and industrial raw materials. In this study, the 12 potential perennial ryegrass miRNAs were identified for the first time by computational approach. Using the newly identified miRNA sequences, the perennial ryegrass mRNA database was further used for BLAST search and detected 33 potential targets of miRNAs. Prediction of potential miRNA target genes revealed their functions involved in various important plant biological processes. Our result should be useful for further investigation into the biological functions of miRNAs in perennial ryegrass. The selected miRNAs representing four families were verified by RT-PCR experiment, indicating that the prediction method that we used to identify the miRNAs was effective.

The CCAAT box-binding transcription factor NF-YA1 controls rhizobial infection

Symbiosis between legume plants and soil rhizobia culminates in the formation of a novel root organ, the 'nodule', containing bacteria differentiated as facultative nitrogen-fixing organelles. MtNF-YA1 is a Medicago truncatula CCAAT box-binding transcription factor (TF), formerly called HAP2-1, highly expressed in mature nodules and required for nodule meristem function and persistence. Here a role for MtNF-YA1 during early nodule development is demonstrated. Detailed expression analysis based on RNA sequencing, quantitiative real-time PCR (qRT-PCR), as well as promoter-β-glucuronidase (GUS) fusions reveal that MtNF-YA1 is first induced at the onset of symbiotic development during preparation for, and initiation and progression of, symbiotic infection. Moreover, using a new knock-out mutant, Mtnf-ya1-1, it is shown that MtNF-YA1 controls infection thread (IT) progression from initial root infection through colonization of nodule tissues. Extensive confocal and electronic microscopic observations suggest that the bulbous and erratic IT growth phenotypes observed in Mtnf-ya1-1 could be a consequence of the fact that walls of ITs in this mutant are thinner and less coherent than in the wild type. It is proposed that MtNF-YA1 controls rhizobial infection progression by regulating the formation and the wall of ITs.

Independent activity of the homologous small regulatory RNAs AbcR1 and AbcR2 in the legume symbiont Sinorhizobium meliloti

The legume symbiont Sinorhizobium meliloti expresses a plethora of small noncoding RNAs (sRNAs) whose function is mostly unknown. Here, we have functionally characterized two tandemly encoded S. meliloti Rm1021 sRNAs that are similar in sequence and structure. Homologous sRNAs (designated AbcR1 and AbcR2) have been shown to regulate several ABC transporters in the related α-proteobacteria Agrobacterium tumefaciens and Brucella abortus. In Rm1021, AbcR1 and AbcR2 exhibit divergent unlinked regulation and are stabilized by the RNA chaperone Hfq. AbcR1 is transcribed in actively dividing bacteria, either in culture, rhizosphere or within the invasion zone of mature alfalfa nodules. Conversely, AbcR2 expression is induced upon entry into stationary phase and under abiotic stress. Only deletion of AbcR1 resulted into a discrete growth delay in rich medium, but both are dispensable for symbiosis. Periplasmic proteome profiling revealed down-regulation of the branched-chain amino acid binding protein LivK by AbcR1, but not by AbcR2. A double-plasmid reporter assay confirmed the predicted specific targeting of the 5′-untranslated region of the livK mRNA by AbcR1 in vivo. Our findings provide evidences of independent regulatory functions of these sRNAs, probably to fine-tune nutrient uptake in free-living and undifferentiated symbiotic rhizobia.

Coordination of flower maturation by a regulatory circuit of three microRNAs

The development of multicellular organisms relies on interconnected genetic programs that control progression through their life cycle. MicroRNAs (miRNAs) and transcription factors (TFs) play key roles in such regulatory circuits. Here, we describe how three evolutionary conserved miRNA-TF pairs interact to form multiple checkpoints during reproductive development of Arabidopsis thaliana. Genetic, cellular, and physiological experiments show that miR159- and miR319-regulated MYB and TCP transcription factors pattern the expression of miR167 family members and their ARF6/8 targets. Coordinated action of these miRNA-TF pairs is crucial for the execution of consecutive hormone-dependent transitions during flower maturation. Cross-regulation includes both cis- and trans-regulatory interactions between these miRNAs and their targets. Our observations reveal how different miRNA-TF pairs can be organized into modules that coordinate successive steps in the plant life cycle.

Development of multicellular organisms relies on properly timed execution of different genetic programs. An example is provided by developmental progression of flowers, which begins with the initiation of individual organs, followed by differentiation, growth, and finally production of the gametes. This article investigates the contribution of three microRNAs (miRNAs) and the transcription factors (TFs) that are regulated by these miRNAs to this process. Two of the miRNA-TF pairs act early to control in parallel the activity of the third miRNA-TF pair, which in turn modulates hormone programs that drive organ maturation and reproduction. Importantly, the two upstream TFs directly interact to regulate expression of the downstream miRNA. The results described here demonstrate how miRNA-TF pairs can be organized into regulatory circuits, with independent miRNA-TF pairs converging on common downstream genes.

Discovering the first microRNA-targeted drug

MicroRNAs (miRNAs) are important post-transcriptional regulators of nearly every biological process in the cell and play key roles in the pathogenesis of human disease. As a result, there are many drug discovery programs that focus on developing miRNA-based therapeutics. The most advanced of these programs targets the liver-expressed miRNA-122 using the locked nucleic acid (LNA)–modified antisense oligonucleotide miravirsen. Here, we describe the discovery of miravirsen, which is currently in phase 2 clinical trials for treatment of hepatitis C virus (HCV) infection.

Analyzing small and long RNAs in plant development using non-radioactive in situ hybridization

In the past decade, hundreds of non-coding RNAs (small and long RNAs) have been identified as crucial elements in developmental processes and stress response in plants. Among small RNAs, the microRNAs or miRNAs control levels of specific mRNA by inhibiting translation or reducing the stability of their mRNA targets through integration into different ribonucleoproteins (RNP). Spatio-temporal expression of small and long RNAs, using reporter genes or in situ hybridization, is essential to understand their functions. We are interested in understanding the role of various non-coding RNAs (including miRNAs) in the regulation of root and nodule development in legumes, which are agriculturally important crops. Here, we present the protocol we are currently using for detection of small and long RNA in model legume plants and tissues, like nodules and roots. The probe selection, as well as the fixation and permeabilization steps allowing to preserve tissues and cell integrity and to increase accessibility to RNA targets, will be specifically discussed.

Whole-mount in situ detection of microRNAs on Arabidopsis tissues using Zip Nucleic Acid probes

MicroRNAs (miRNAs) affect fundamental processes of development. In plants miRNAs regulate organ development, transition to flowering, and responses to abiotic/biotic stresses. To understand the biological role of miRNAs, in addition to identifying their targeted transcripts, it is necessary to characterize the spatiotemporal regulation of their expression. Many methods have been used to define the set of organ-specific miRNAs by tissue dissection and miRNA profiling but none of them can describe their tissue and cellular distribution at the high resolution provided by in situ hybridization (ISH). This article describes the setup and optimization of a whole-mount ISH protocol to target endogenous miRNAs on intact Arabidopsis seedlings using DIG-labeled Zip Nucleic Acid (ZNA) oligonucleotide probes. Automation of the main steps of the procedure by robotized liquid handling has also been implemented in the protocol for best reproducibility of results, enabling running of ISH experiments at high throughput.

"And yet it moves": cell-to-cell and long-distance signaling by plant microRNAs

MicroRNAs (miRNAs) are key regulators of numerous genes in many eukaryotes. Some plant miRNAs are involved in developmental and physiological processes that require intercellular or inter-organ signaling. Movement of other small RNAs within plants has been established. Recent findings also demonstrate intercellular signaling by miRNAs and strongly support that a subset of these regulatory molecules move from one cell to another or over long distances. Phloem exudates contain diverse miRNAs and at least two of them, involved in responses to nutrient availability, are transmitted through grafts, indicating long-distance movement. Two miRNAs that regulate developmental processes are present in cells outside their domains of expression. Several results strongly support that one of them moves from cell to cell. Research on a mutant affected in plasmodesmata trafficking indicates that these intercellular channels are required for transmission of miRNA activity to adjacent cells. Moreover, ARGONAUTE proteins might be involved in the regulation of miRNA trafficking. Hypothesis on the features and mechanisms that may determine miRNA mobility are presented. Future challenges include identifying other mobile miRNAs; demonstrating that miRNA movement is required for non-cell autonomous action; and characterizing the mechanisms of translocation and genetic pathways that regulate miRNA movement.

In situ detection of mature miRNAs in plants using LNA-modified DNA probes

MicroRNAs (miRNAs) play important roles in development in plants, and some miRNAs show developmentally regulated organ- and tissue-specific expression patterns. Therefore, in situ detection of mature miRNAs is important for understanding the functions of both miRNAs and their targets. The construction of promoter-reporter fusions and examination of their in planta expression have been widely used and the results obtained thus far are rather informative; however, in some cases, the length of promoter that contains the entire regulatory elements is difficult to determine. In addition, traditional in situ hybridization with the antisense RNA fragment as the probe usually fails to detect miRNAs because the mature miRNAs are too short (~21 nt) to exhibit stable hybridization signals. In recent years, the locked nucleic acid (LNA)-modified DNA probe has been successfully used in animals and plants to detect small RNAs. Here, we describe a modified protocol using LNA-modified DNA probes to detect mature miRNAs in plant -tissues, including the design of LNA probes and detailed steps for the in situ hybridization experiment, using Arabidopsis miR165 as an example.

Four AUXIN RESPONSE FACTOR genes downregulated by microRNA167 are associated with growth and development in Oryza sativa

MicroRNA167 (miR167), as a conserved miRNA, has been implicated in auxin signalling by regulating the expression of certain auxin response factor (ARF) genes to determine the plant developmental process. Among the 10 MIR167 genes of rice, the precursor structures derived from MIR167a, MIR167b and MIR167c produce miR167 with high efficiency. To explore the biological function of miR167 in rice, four of its predicted target genes, OsARF6, OsARF12, OsARF17 and OsARF25, were identified in vivo. Although the expression levels of miR167 and its target OsARFs did not show an obvious negative correlation, the enhanced miR167 level in transgenic rice overexpressing miR167 resulted in a substantial decrease in mRNA levels of the four OsARF genes. Moreover, the transgenic rice plants were small in stature with remarkably reduced tiller number. These results suggest that miR167 is important for the appropriate expression of at least four OsARFs, which mediate the auxin response, to contribute to the normal growth and development of rice.

Expression and tissue-specific localization of nitrate-responsive miRNAs in roots of maize seedlings

Nitrogen availability seriously affects crop productivity and environment. The knowledge of post-transcriptional regulation of plant response to nutrients is important to improve nitrogen use efficiency of crop. This research was aimed at understanding the role of miRNAs in the molecular control of plant response to nitrate. The expression profiles of six mature miRNAs were deeply studied by quantitative real time polymerase chain reaction and in situ hybridization (ISH). To this aim, a novel optimized protocol was set up for the use of digoxygenin-labelled Zip Nucleic Acid-modified oligonucleotides as probes for ISH. Significant differences in miRNAs' transcripts accumulation were evidenced between nitrate-supplied and nitrate-depleted roots. Real-time PCR analyses and in situ detection of miRNA confirmed the array data and allowed us to evidence distinct miRNAs spatio-temporal expression patterns in maize roots. Our results suggest that a prolonged nitrate depletion may induce post-transcriptionally the expression of target genes by repressing the transcription of specific miRNAs. In particular, the repression of the transcription of miR528a/b, miR528a*/b*, miR169i/j/k, miR169i*/j*/k*, miR166j/k/n and miR408/b upon nitrate shortage could represent a crucial step integrating nitrate signals into developmental changes in maize roots.

Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

Identification and characterization of microRNAs in Phaseolus vulgaris by high-throughput sequencing

Background

MicroRNAs (miRNAs) are endogenously encoded small RNAs that post-transcriptionally regulate gene expression. MiRNAs play essential roles in almost all plant biological processes. Currently, few miRNAs have been identified in the model food legume Phaseolus vulgaris (common bean). Recent advances in next generation sequencing technologies have allowed the identification of conserved and novel miRNAs in many plant species. Here, we used Illumina's sequencing by synthesis (SBS) technology to identify and characterize the miRNA population of Phaseolus vulgaris.

Results

Small RNA libraries were generated from roots, flowers, leaves, and seedlings of P. vulgaris. Based on similarity to previously reported plant miRNAs,114 miRNAs belonging to 33 conserved miRNA families were identified. Stem-loop precursors and target gene sequences for several conserved common bean miRNAs were determined from publicly available databases. Less conserved miRNA families and species-specific common bean miRNA isoforms were also characterized. Moreover, novel miRNAs based on the small RNAs were found and their potential precursors were predicted. In addition, new target candidates for novel and conserved miRNAs were proposed. Finally, we studied organ-specific miRNA family expression levels through miRNA read frequencies.

Conclusions

This work represents the first massive-scale RNA sequencing study performed in Phaseolus vulgaris to identify and characterize its miRNA population. It significantly increases the number of miRNAs, precursors, and targets identified in this agronomically important species. The miRNA expression analysis provides a foundation for understanding common bean miRNA organ-specific expression patterns. The present study offers an expanded picture of P. vulgaris miRNAs in relation to those of other legumes.

Kiwifruit floral gene APETALA2 is alternatively spliced and accumulates in aberrant indeterminate flowers in the absence of miR172

In Arabidopsis, the identity of perianth and reproductive organs are specified by antagonistic action of two floral homeotic genes, APETALA2 (AP2) and AGAMOUS (AG). AP2 is also negatively regulated by an evolutionary conserved interaction with a microRNA, miR172, and has additional roles in general plant development. A kiwifruit gene with high levels of homology to AP2 and AP2-like genes from other plant species was identified. The transcript was abundant in the kiwifruit flower, particularly petal, suggesting a role in floral organ identity. Splice variants were identified, all containing both AP2 domains, including a variant that potentially produces a shorter transcript without the miRNA172 targeting site. Increased AP2 transcript accumulation was detected in the aberrant flowers of the mutant ‘Pukekohe dwarf’ with multiple perianth whorls and extended petaloid features. In contrast to normal kiwifruit flowers, the aberrant flowers failed to accumulate miR172 in the developing whorls, although accumulation was detected at the base of the flower. An additional role during dormancy in kiwifruit was proposed based on AP2 transcript accumulation in axillary buds before and after budbreak.

Give It AGO: The Search for miRNA-Argonaute Sorting Signals in Arabidopsis thaliana Indicates a Relevance of Sequence Positions Other than the 5'-Position Alone

The specific recognition of miRNAs by Argonaute (AGO) proteins, the effector proteins of the RNA-induced silencing complex, constitutes the final step of the biogenesis of miRNAs and is crucial for their target interaction. In the genome of Arabidopsis thaliana (Ath), 10 different AGO proteins are encoded and the sorting decision, which miRNA associates with which AGO protein, was reported to depend exclusively on the identity of the 5'-sequence position of mature miRNAs. Hence, with only four different bases possible, a 5'-position-only sorting signal would not suffice to specifically target all 10 different AGOs individually or would suggest redundant AGO action. Alternatively, other and as of yet unidentified sorting signals may exist. We analyzed a dataset comprising 117 Ath-miRNAs with clear sorting preference to either AGO1, AGO2, or AGO5 as identified in co-immunoprecipitation experiments combined with sequencing. While mutual information analysis did not identify any other single position but the 5'-nucleotide to be informative for the sorting at sufficient statistical significance, significantly better than random classification results using Random Forests nonetheless suggest that additional positions and combinations thereof also carry information with regard to the AGO sorting. Positions 2, 6, 9, and 13 appear to be of particular importance. Furthermore, uracil bases at defined positions appear to be important for the sorting to AGO2 and AGO5, in particular. No predictive value was associated with miRNA length or base pair binding pattern in the miRNA:miRNA* duplex. From inspecting available AGO gene expression data in Arabidopsis, we conclude that the temporal and spatial expression profile may also contribute to the fine-tuning of miRNA sorting and function.

High-resolution experimental and computational profiling of tissue-specific known and novel miRNAs in Arabidopsis

Small non-coding RNAs (ncRNAs) are key regulators of plant development through modulation of the processing, stability, and translation of larger RNAs. We present small RNA data sets comprising more than 200 million aligned Illumina sequence reads covering all major cell types of the root as well as four distinct developmental zones. MicroRNAs (miRNAs) constitute a class of small ncRNAs that are particularly important for development. Of the 243 known miRNAs, 133 were found to be expressed in the root, and most showed tissue- or zone-specific expression patterns. We identified 66 new high-confidence miRNAs using a computational pipeline, PIPmiR, specifically developed for the identification of plant miRNAs. PIPmiR uses a probabilistic model that combines RNA structure and expression information to identify miRNAs with high precision. Knockdown of three of the newly identified miRNAs results in altered root growth phenotypes, confirming that novel miRNAs predicted by PIPmiR have functional relevance.

Detection of microRNAs in plants by in situ hybridisation

MicroRNAs (miRNAs) are short, about 21 nucleotides in length, non-coding, regulatory RNA molecules representing a new layer in post-transcriptional gene expression regulation. Spatial and temporal analysis of miRNA accumulation by in situ analyses is the prerequisite of understanding the precise biological functions of miRNAs. Since miRNAs are very short molecules, their in situ analysis is technically demanding. Our method is based on the usage of highly sensitive LNA-modified oligonucleotide probes. LNA modification significantly enhances the sensitivity and specificity of miRNA detecting probes and provides relatively easy in situ miRNA detection. Here, we describe a protocol for this challenging technique step by step, in order to help every user to achieve success.

Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa

Background

Pollen development in flowering plants requires strict control of the gene expression program and genetic information stability by mechanisms possibly including the miRNA pathway. However, our understanding of the miRNA pathway in pollen development remains limited, and the dynamic profile of miRNAs in developing pollen is unknown.

Results

Using next-generation sequencing technology, we pyrosequenced small RNA populations from rice uninucleate microspores to tricellular pollen and control sporophytic tissues at the genome-wide level. We identified 292 known miRNAs, including members of all 20 families conserved in plants, and 75 novel miRNAs. Of the 292 known miRNAs, 202 were expressed, with 103 enriched, in developing pollen. More than half of these novel miRNAs displayed pollen-or stage-specific expression. Furthermore, analyzing the 367 miRNAs and their predicted targets indicated that correlation in expression profiles of pollen-enriched known miRNAs and their targets significantly differs from that of sporophyte-enriched known miRNAs and their targets in some functional terms, while novel miRNAs appeared to negatively regulate their targets. Importantly, gene ontology abundance analysis demonstrated chromatin assembly and disassembly was important in the targets of bicellular pollen-expressed novel miRNAs. Principal component analysis revealed pollen of all three stages was discriminated from sporophytes, largely because of the novel and non-conserved known miRNAs.

Conclusions

Our study, for the first time, revealed the differences in composition and expression profiles of miRNAs between developing pollen and sporophytes, with novel and non-conserved known miRNAs the main contributors. Our results suggest the important roles of the miRNA pathway in pollen development.

Insights from ANA-grade angiosperms into the early evolution of CUP-SHAPED COTYLEDON genes

BACKGROUND AND AIMS

The closely related NAC family genes NO APICAL MERISTEM (NAM) and CUP-SHAPED COTYLEDON3 (CUC3) regulate the formation of boundaries within and between plant organs. NAM is post-transcriptionally regulated by miR164, whereas CUC3 is not. To gain insight into the evolution of NAM and CUC3 in the angiosperms, we analysed orthologous genes in early-diverging ANA-grade angiosperms and gymnosperms.

METHODS

We obtained NAM- and CUC3-like sequences from diverse angiosperms and gymnosperms by a combination of reverse transcriptase PCR, cDNA library screening and database searching, and then investigated their phylogenetic relationships by performing maximum-likelihood reconstructions. We also studied the spatial expression patterns of NAM, CUC3 and MIR164 orthologues in female reproductive tissues of Amborella trichopoda, the probable sister to all other flowering plants.

KEY RESULTS

Separate NAM and CUC3 orthologues were found in early-diverging angiosperms, but not in gymnosperms, which contained putative orthologues of the entire NAM + CUC3 clade that possessed sites of regulation by miR164. Multiple paralogues of NAM or CUC3 genes were noted in certain taxa, including Brassicaceae. Expression of NAM, CUC3 and MIR164 orthologues from Am. trichopoda was found to co-localize in ovules at the developmental boundary between the chalaza and nucellus.

CONCLUSIONS

The NAM and CUC3 lineages were generated by duplication, and CUC3 was subsequently lost regulation by miR164, prior to the last common ancestor of the extant angiosperms. However, the paralogous NAM clade genes CUC1 and CUC2 were generated by a more recent duplication, near the base of Brassicaceae. The function of NAM and CUC3 in defining a developmental boundary in the ovule appears to have been conserved since the last common ancestor of the flowering plants, as does the post-transcriptional regulation in ovule tissues of NAM by miR164.

Differential effects of Cucumber mosaic virus satellite RNAs in the perturbation of microRNA-regulated gene expression in tomato

Viral infections generally cause disease symptoms by interfering with the microRNA (miRNA)-mediated regulation of gene expression of host plants. In tomato leaves, the accumulation levels of eleven miRNAs and ten target mRNAs were enhanced by different degrees upon Cucumber mosaic virus (CMV)-Fny and Tomato aspermy virus (TAV)-Bj infections. The ability of CMV-Fny to interfere with miRNA pathway was dramatically suppressed in the addition of the benign satellite (sat) RNA variant (satYn12), but was slightly affected when CMV-Fny was co-inoculated with the aggressive satRNA variant (satT1). In plants harboring the infection of CMV-FnyΔ2b (a CMV-Fny 2b-deletion mutant), the unaltered miRNAs and target mRNAs levels compared with mock inoculated plants indicated that 2b ORF was essential for perturbation of miRNA metabolism. When the amounts of viral open reading frames (ORFs) in these infections were quantified, we found satYn12 caused a higher reduction of CMV-Fny accumulation levels than satT1. These results indicate the complex mechanism by which satRNAs participate in CMV-tomato interaction, and suggest that the severity of disease symptoms positively correlates to some extent with the perturbation of miRNA pathway in tomato.

Misregulation of AUXIN RESPONSE FACTOR 8 underlies the developmental abnormalities caused by three distinct viral silencing suppressors in Arabidopsis

In Arabidopsis, micro (mi)RNAs and trans-acting (ta-si)RNAs synthesized directly or indirectly through the DICER-LIKE-1 (DCL1) ribonuclease have roles in patterning and hormonal responses, while DCL2,3,4-dependent small-interfering (si)RNAs are mainly involved in silencing of transposable elements and antiviral defense. Viral suppressors of RNA silencing (VSRs) produced by phytoviruses to counter plant defense may perturb plant developmental programs because of the collision of their inhibitory effects with the regulatory action of endogenous miRNAs and ta-siRNAs. This could explain the similar developmental aberrations displayed by Arabidopsis miRNA/ta-siRNA pathway mutants, including dcl1, and by some VSR-expressing plants. Nonetheless, the molecular bases for these morphological aberrations have remained mysterious, and their contribution to viral disease symptoms/virulence unexplored. The extent of VSR inhibitory actions to other types of endogenous small RNAs remains also unclear. Here, we present an in-depth analysis of transgenic Arabidopsis expressing constitutively HcPro, P19 and P15, three unrelated VSRs. We show that VSR expression has comparable, yet modest effects on known miRNA and ta-siRNA target RNA levels, similar to those observed using an hypomorphic dcl1 mutation. However, by combining results of transcriptome studies with deep-sequencing data from immuno-precipitated small RNAs, additional, novel endogenous targets of miRNA and ta-siRNA were identified, unraveling an unsuspected complexity in the origin and scope-of-action of these molecules. Other stringent analyses pinpointed misregulation of the miR167 target AUXIN RESPONSE FACTOR 8 (ARF8) as a major cause for the developmental aberrations exhibited by VSR transgenic plants, but also for the phenotypes induced during normal viral infection caused by the HcPro-encoding Turnip mosaic virus (TuMV). Neither RNA silencing, its suppression by VSRs, nor the virulence/accumulation of TuMV was altered by mutations in ARF8. These findings have important implications for our understanding of viral disease symptoms and small RNA-directed regulation of plant growth/development.

In the plant and animal RNA silencing pathways, small RNA molecules known as micro (mi)RNA and short-interfering (si)RNAs have key roles in development and antiviral defense, respectively. In turn, viruses counteract this defense by deploying specific virulence factors, referred to as Viral Suppressors of RNA silencing (VSRs), which target distinct steps of the host silencing machinery. In the model plant species Arabidopsis thaliana, transgenic expression of distinct VSRs often incurs a set of strikingly recurrent developmental anomalies that resemble those triggered by viral infections. While these defects have been assumed to result from a general interference of VSRs with silencing-based mechanisms controlling cellular growth, their exact molecular basis has remained largely elusive. Here, we address this issue by demonstrating that misregulation of a single transcript encoding the AUXIN RESPONSE FACTOR 8, a target of miR167, underlies most, if not all, of the defects caused by VSR expression, both in transgenic and in an authentic infection context. Our study also highlights the value of VSRs as generic tools for the discovery or validation of endogenous RNA silencing targets. These results also have implications for our understanding of small RNA-based regulations in plants, and shed light on the possible origin of some of the symptoms elicited by viral diseases.