Criteria for annotation of plant MicroRNAs

MicroRNAs (miRNAs) are approximately 21 nucleotide noncoding RNAs produced by Dicer-catalyzed excision from stem-loop precursors. Many plant miRNAs play critical roles in development, nutrient homeostasis, abiotic stress responses, and pathogen responses via interactions with specific target mRNAs. miRNAs are not the only Dicer-derived small RNAs produced by plants: A substantial amount of the total small RNA abundance and an overwhelming amount of small RNA sequence diversity is contributed by distinct classes of 21- to 24-nucleotide short interfering RNAs. This fact, coupled with the rapidly increasing rate of plant small RNA discovery, demands an increased rigor in miRNA annotations. Herein, we update the specific criteria required for the annotation of plant miRNAs, including experimental and computational data, as well as refinements to standard nomenclature.

Large-Scale Sequencing Reveals 21U-RNAs and Additional MicroRNAs and Endogenous siRNAs in C. elegans

We sequenced approximately 400,000 small RNAs from Caenorhabditis elegans. Another 18 microRNA (miRNA) genes were identified, thereby extending to 112 our tally of confidently identified miRNA genes in C. elegans. Also observed were thousands of endogenous siRNAs generated by RNA-directed RNA polymerases acting preferentially on transcripts associated with spermatogenesis and transposons. In addition, a third class of nematode small RNAs, called 21U-RNAs, was discovered. 21U-RNAs are precisely 21 nucleotides long, begin with a uridine 5'-monophosphate but are diverse in their remaining 20 nucleotides, and appear modified at their 3'-terminal ribose. 21U-RNAs originate from more than 5700 genomic loci dispersed in two broad regions of chromosome IV-primarily between protein-coding genes or within their introns. These loci share a large upstream motif that enables accurate prediction of additional 21U-RNAs. The motif is conserved in other nematodes, presumably because of its importance for producing these diverse, autonomously expressed, small RNAs (dasRNAs).

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome

MicroRNAs (miRNAs) regulate the expression of target mRNAs in plants and animals [1]. Plant miRNA targets have been predicted on the basis of their extensive and often conserved complementarity to the miRNAs [2-4], as well as on miRNA overexpression experiments [5]; many of these target predictions have been confirmed by isolation of the products of miRNA-directed cleavage. Here, we present a transcriptome-wide experimental method, called "degradome sequencing," to directly detect cleaved miRNA targets without relying on predictions or overexpression. The 5' ends of polyadenylated, uncapped mRNAs from Arabidopsis were directly sampled, resulting in an empirical snapshot of the degradome. miRNA-mediated-cleavage products were easily discerned from an extensive background of degraded mRNAs, which collectively covered the majority of the annotated transcriptome. Many previously known Arabidopsis miRNA targets were confirmed, and several novel targets were also discovered. Quantification of cleavage fragments revealed that those derived from TAS transcripts, which are unusual in their production of abundant secondary small interfering RNAs (siRNAs), accumulated to very high levels. A subset of secondary siRNAs are also known to direct cleavage of targets in trans[6]; degradome sequencing revealed many cleaved targets of these trans-acting siRNAs (ta-siRNAs). This empirical method is broadly applicable to the discovery and quantification of cleaved targets of small RNAs without a priori predictions.

Classification and comparison of small RNAs from plants

Regulatory small RNAs, which range in size from 20 to 24 nucleotides, are ubiquitous components of endogenous plant transcriptomes, as well as common responses to exogenous viral infections and introduced double-stranded RNA (dsRNA). Endogenous small RNAs derive from the processing of helical RNA precursors and can be categorized into several groups based on differences in biogenesis and function. A major distinction can be observed between small RNAs derived from single-stranded precursors with a hairpin structure [referred to here as hairpin RNAs (hpRNAs)] and those derived from dsRNA precursors [small interfering RNAs (siRNAs)]. hpRNAs in plants can be divided into two secondary groups: microRNAs and those that are not microRNAs. The currently known siRNAs fall mostly into one of three secondary groups: heterochromatic siRNAs, secondary siRNAs, and natural antisense transcript siRNAs. Tertiary subdivisions can be identified within many of the secondary classifications as well. Comparisons between the different classes of plant small RNAs help to illuminate key goals for future research.

Initiation of RPS2-specified disease resistance in Arabidopsis is coupled to the AvrRpt2-directed elimination of RIN4

Plants have evolved a sophisticated innate immune system to recognize invading pathogens and to induce a set of host defense mechanisms resulting in disease resistance. Pathogen recognition is often mediated by plant disease resistance (R) proteins that respond specifically to one or a few pathogen-derived molecules. This specificity has led to suggestions of a receptor-ligand mode of R protein function. Delivery of the bacterial effector protein AvrRpt2 by Pseudomonas syringae specifically induces disease resistance in Arabidopsis plants expressing the RPS2 R protein. We demonstrate that RPS2 physically interacts with Arabidopsis RIN4 and that AvrRpt2 causes the elimination of RIN4 during activation of the RPS2 pathway. AvrRpt2-mediated RIN4 elimination also occurs in the rps2, ndr1, and Atrar1 mutant backgrounds, demonstrating that this activity can be achieved independent of an RPS2-mediated signaling pathway. Therefore, we suggest that RPS2 initiates signaling based upon perception of RIN4 disappearance rather than direct recognition of AvrRpt2.

CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets

MicroRNAs (miRNAs) are approximately 20- to 22-nt long endogenous RNA sequences that play a critical role in the regulation of gene expression in eukaryotic genomes. Confident identification of miRNA targets is vital to understand their functions. Currently available computational algorithms for miRNA target prediction have diverse degrees of sensitivity and specificity and as a consequence each predicted target generally requires experimental confirmation. miRNAs and other small RNAs that direct endonucleolytic cleavage of target mRNAs produce diagnostic uncapped, polyadenylated mRNA fragments. Degradome sequencing [also known as PARE (parallel analysis of RNA ends) and GMUCT (genome-wide mapping of uncapped transcripts)] samples the 5'-ends of uncapped mRNAs and can be used to discover in vivo miRNA targets independent of computational predictions. Here, we describe a generalizable computational pipeline, CleaveLand, for the detection of cleaved miRNA targets from degradome data. CleaveLand takes as input degradome sequences, small RNAs and an mRNA database and outputs small RNA targets. CleaveLand can thus be applied to degradome data from any species provided a set of mRNA transcripts and a set of query miRNAs or other small RNAs are available.

AVAILABILITY

The code and documentation for CleaveLand is freely available under a GNU license at http://www.bio.psu.edu/people/faculty/Axtell/AxtellLab/Software.html

A two-hit trigger for siRNA biogenesis in plants

In Arabidopsis, microRNA-directed cleavage can define one end of RNAs that then generate phased siRNAs. However, most miRNA-targeted RNAs do not spawn siRNAs, suggesting the existence of additional determinants within those that do. We find that in moss, phased siRNAs arise from regions flanked by dual miR390 cleavage sites. AtTAS3, an siRNA locus important for development and conserved among higher plants, also has dual miR390 complementary sites. Both sites bind miR390 in vitro and are functionally required in Arabidopsis, but cleavage is undetectable at the 5' site--demonstrating that noncleavable sites can be functional in plants. Phased siRNAs also emanate from the bounded regions of every Arabidopsis gene with two known microRNA/siRNA complementary sites, but only rarely from genes with single sites. Therefore, two "hits,"--often, but not always, two cleavage events--constitute a conserved trigger for siRNA biogenesis, a finding with implications for recognition and silencing of aberrant RNA.

Antiquity of microRNAs and their targets in land plants

MicroRNAs (miRNAs) affect the morphology of flowering plants by the posttranscriptional regulation of genes involved in critical developmental events. Understanding the spatial and temporal dynamics of miRNA activity during development is therefore central for understanding miRNA functions. We describe a microarray suitable for detection of plant miRNAs. Profiling of Arabidopsis thaliana miRNAs during normal development extends previous expression analyses, highlighting differential expression of miRNA families within specific organs and tissue types. Comparison of our miRNA expression data with existing mRNA microarray data provided a global intersection of plant miRNA and mRNA expression profiles and revealed that tissues in which a given miRNA is highly expressed are unlikely to also show high expression of the corresponding targets. Expression profiling was also used in a phylogenetic survey to test the depth of plant miRNA conservation. Of the 23 families of miRNAs tested, expression of 11 was detected in a gymnosperm and eight in a fern, directly demonstrating that many plant miRNAs have remained essentially unchanged since before the emergence of flowering plants. We also describe an empirical strategy for detecting miRNA target genes from unsequenced transcriptomes and show that targets in nonflowering plants as deeply branching as ferns and mosses are homologous to the targets in Arabidopsis. Therefore, several individual miRNA regulatory circuits have ancient origins and have remained intact throughout the evolution and diversification of plants.

Common functions for diverse small RNAs of land plants

Endogenous small RNAs, including microRNAs (miRNAs) and short interfering RNAs (siRNAs), are critical components of plant gene regulation. Some abundant miRNAs involved in developmental control are conserved between anciently diverged plants, while many other less-abundant miRNAs appear to have recently emerged in the Arabidopsis thaliana lineage. Using large-scale sequencing of small RNAs, we extended the known diversity of miRNAs in basal plants to include 88 confidently annotated miRNA families in the moss Physcomitrella patens and 44 in the lycopod Selaginella moellendorffii. Cleavage of 29 targets directed by 14 distinct P. patens miRNA families and a trans-acting siRNA (ta-siRNA) was experimentally confirmed. Despite a core set of 12 miRNA families also expressed in angiosperms, weakly expressed and apparently lineage-specific miRNAs accounted for the majority of miRNA diversity in both species. Nevertheless, the molecular functions of several of these lineage-specific small RNAs matched those of angiosperms, despite dissimilarities in the small RNA sequences themselves, including small RNAs that mediated negative feedback regulation of the miRNA pathway and miR390-dependent ta-siRNAs that guided the cleavage of AUXIN RESPONSE FACTOR mRNAs. Diverse, lineage-specific, small RNAs can therefore perform common biological functions in plants.

Vive la différence: biogenesis and evolution of microRNAs in plants and animals

MicroRNAs are pervasive in both plants and animals, but many aspects of their biogenesis, function and evolution differ. We reveal how these differences contribute to characteristic features of microRNA evolution in the two kingdoms.

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.

Evolution of plant microRNAs and their targets

MicroRNAs (miRNAs) are a specialized class of small silencing RNAs that regulate gene expression. They have a limited phylogenetic distribution among eukaryotes, suggestive of at least two independent origins from an ancestral small RNA-producing pathway. A set of 21 abundantly expressed miRNAs are clearly conserved among the angiosperms; many of these function to regulate transcription factors involved in developmental control. Recent experiments have uncovered a much larger set of weakly expressed, less conserved miRNAs in plants, and this group has provided insights into the origins of miRNAs and their targets. These data have provided a coherent set of hypotheses explaining the birth, selection and death of miRNAs in land plants.

Transcriptome-wide identification of microRNA targets in rice

MicroRNA (miRNA)-guided target RNA expression is vital for a wide variety of biological processes in eukaryotes. Currently, miRBase (version 13) lists 142 and 353 miRNAs from Arabidopsis and rice (Oryza sativa), respectively. The integration of miRNAs in diverse biological networks relies upon the confirmation of their RNA targets. In contrast with the well-characterized miRNA targets that are cleaved in Arabidopsis, only a few such targets have been confirmed in rice. To identify small RNA targets in rice, we applied the 'degradome sequencing' approach, which globally identifies the remnants of small RNA-directed target cleavage by sequencing the 5' ends of uncapped RNAs. One hundred and sixty targets of 53 miRNA families (24 conserved and 29 rice-specific) and five targets of TAS3-small interfering RNAs (siRNAs) were identified. Surprisingly, an additional conserved target for miR398, which has not been reported so far, has been validated. Besides conserved homologous transcripts, 23 non-conserved genes for nine conserved miRNAs and 56 genes for 29 rice-specific miRNAs were also identified as targets. Besides miRNA targets, the rice degradome contained fragments derived from MIRNA precursors. A closer inspection of these fragments revealed a unique pattern distinct from siRNA-producing loci. This attribute can serve as one of the ancillary criteria for separating miRNAs from siRNAs in plants.

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

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

MicroRNAs from the parasitic plant Cuscuta campestris target host messenger RNAs

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

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.

Genetic and molecular evidence that the Pseudomonas syringae type III effector protein AvrRpt2 is a cysteine protease

Upon delivery to the plant cell during infection, the Pseudomonas syringae effector protein AvrRpt2 undergoes proteolytic processing, enhances pathogen virulence and causes the elimination of the Arabidopsis RIN4 protein. A structure-prediction method was employed in order to investigate possible biochemical functions of AvrRpt2. Results of a secondary structure prediction algorithm suggest that the functional C-terminal portion of AvrRpt2 is a cysteine protease. Mutation of predicted catalytic residues within this portion of AvrRpt2 abolished in planta processing, elimination of Arabidopsis RIN4, and the ability to trigger an RPS2-specific resistance response. These data indicate that AvrRpt2 is most likely a sequence divergent cysteine protease whose activity is required for elimination of RIN4 during infection.

The helicase domain of the TMV replicase proteins induces the N-mediated defence response in tobacco

Tobacco mosaic virus (TMV) induces the hypersensitive response (HR) in tobacco plants containing the N gene. This defence response is characterized by cell death at the site of virus infection and inhibition of viral replication and movement. A previous study indicated that a portion of the TMV replicase containing a putative helicase domain is involved in HR induction. Here, this observation is confirmed and extended by showing that non-viral expression of a 50 kDa TMV helicase fragment (p50) is sufficient to induce the N-mediated HR in tobacco. Like the HR elicited by TMV infection, transgenic expression of p50 induces a temperature-sensitive defence response. We demonstrate that recombinant p50 protein has ATPase activity, as suggested by the presence of conserved sequence motifs found in ATPase/helicase enzymes. A point mutation that alters one of these motifs abolishes ATPase activity in vitro but does not affect HR induction. These results suggest that features of the TMV helicase domain, independent of its enzymatic activity, are recognized by N-containing tobacco to induce TMV resistance.

RNA secondary structural determinants of miRNA precursor processing in Arabidopsis

MicroRNAs (miRNAs) are excised from hairpin structures within primary miRNAs (pri-miRNAs). Most animal pri-miRNAs are processed by two cleavages, the first at a loop-distal site approximately 11 nucleotides (nt) from the end of the hairpin and the second approximately 22 nt beyond the first. To identify RNA structural determinants of miRNA processing in plants, we analyzed the functional consequences of changing the secondary structure of the lower (loop-distal), middle (miRNA:miRNA(*)), and upper (loop-proximal) stems of the hairpin in two different pri-miRNAs. Closing bulges immediately below the loop-distal cleavage sites increased the accumulation of accurately cleaved precursor miRNAs but decreased the abundance of the mature miRNAs. A pri-miRNA variant with an unpaired lower stem was not processed, and variants with a perfectly paired middle or upper stem were processed normally. Bioinformatic analysis of pri-miRNA structures, together with physical mapping of initial cleavage sites and in vitro processing of pri-miRNA, reveals that the first, loop-distal cleavage is often at a distance of approximately 15 nt from an unpaired region. Hence, a common determinant of the rate and location of the initial pri-miRNA cleavage is an imperfectly base-paired duplex of approximately 15 nt between the miRNA:miRNA(*) duplex and either a less structured region of the lower stem or its end.

NPK1, an MEKK1-like mitogen-activated protein kinase kinase kinase, regulates innate immunity and development in plants

Mitogen-activated protein kinase (MAPK) cascades are rapidly activated upon plant recognition of invading pathogens. Here, we describe the use of virus-induced gene silencing (VIGS) to study the role of candidate plant MAP kinase kinase kinase (MAPKKK) homologs of human MEKK1 in pathogen-resistance pathways. We demonstrate that silencing expression of a tobacco MAPKKK, Nicotiana Protein Kinase 1 (NPK1), interferes with the function of the disease-resistance genes N, Bs2, and Rx, but does not affect Pto- and Cf4-mediated resistance. Further, NPK1-silenced plants also exhibit reduced cell size, defective cytokinesis, and an overall dwarf phenotype. Our results provide evidence that NPK1 functions in the regulation of N-, Bs2-, and Rx-mediated resistance responses and may play a role in one or more MAPK cascades, regulating multiple cellular processes.

Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome

Expression profiling of the 5' ends of uncapped mRNAs ("degradome" sequencing) can be used to empirically catalog microRNA (miRNA) targets, to probe patterns of miRNA hairpin processing, to examine mRNA decay, and to analyze accumulation of endogenous short interfering RNA (siRNA) precursors. We sequenced and analyzed the degradome of the moss Physcomitrella patens, an important model system for functional genomic analyses in plant evolution. A total of 52 target mRNAs of 27 different Physcomitrella miRNA families were identified. Many targets of both more conserved and less conserved miRNA families encoded putative regulatory proteins. Remnants of MIRNA hairpin processing also populated the degradome data and indicated an unusual "loop-first" mode of precise processing for the MIR319 gene family. Precise loop-first processing was confirmed for native Physcomitrella, rice, and Arabidopsis MIR319 hairpins, as well as an Arabidopsis artificial MIRNA (aMIRNA) based upon a MIR319 backbone. MIR319 is thus a conserved exception to the general rule of loop-last processing of MIRNA hairpins. Loop-first MIR319 processing may contribute to the high efficacy of a widely used MIR319-based strategy for aMIRNA production in plants.

Arabidopsis lyrata small RNAs: transient MIRNA and small interfering RNA loci within the Arabidopsis genus

Twenty-one-nucleotide microRNAs (miRNAs) and 24-nucleotide Pol IV-dependent small interfering RNAs (p4-siRNAs) are the most abundant types of small RNAs in angiosperms. Some miRNAs are well conserved among different plant lineages, whereas others are less conserved, and it is not clear whether less-conserved miRNAs have the same functionality as the well conserved ones. p4-siRNAs are broadly produced in the Arabidopsis genome, sometimes from active hot spot loci, but it is unknown whether individual p4-siRNA hot spots are retained as hot spots between plant species. In this study, we compare small RNAs in two closely related species (Arabidopsis thaliana and Arabidopsis lyrata) and find that less-conserved miRNAs have high rates of divergence in MIRNA hairpin structures, mature miRNA sequences, and target-complementary sites in the other species. The fidelity of miRNA biogenesis from many less-conserved MIRNA hairpins frequently deteriorates in the sister species relative to the species of first discovery. We also observe that p4-siRNA occupied loci have a slight tendency to be retained as p4-siRNA loci between species, but the most active A. lyrata p4-siRNA hot spots are generally not syntenic to the most active p4-siRNA hot spots of A. thaliana. Altogether, our findings indicate that many MIRNAs and most p4-siRNA hot spots are rapidly changing and evolutionarily transient within the Arabidopsis genus.

Bidirectional processing of pri-miRNAs with branched terminal loops by Arabidopsis Dicer-like1

miRNAs originate from primary transcripts (pri-miRNAs) with characteristic stem-loop structures. Accurate processing of pri-miRNAs is required for functional miRNAs. Here, using pri-miR166 family as a paradigm, we report the decisive role of pri-miRNA terminal loops in miRNA biogenesis. We found that multi-branched terminal loops in pri-miR166s substantially suppressed miR166 expression in vivo. Unlike canonical processing of pri-miRNAs, terminal-loop-branched (TLBed) pri-miRNAs can be processed by Dicer-like1 (DCL1) complexes bi-directionally: from base to loop and from loop to base, resulting in productive and abortive processing of miRNAs, respectively. In either case, DCL1 complexes canonically cut pri-miRNAs at a distance of 16-17 base pairs (bp) from a reference single-stranded loop region. DCL1 also adjusts processing sites toward an internal loop through its helicase domain. Thus, these results provide new insight into the poorly understood processing mechanism of pri-miRNAs with complicated secondary structures.

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

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