Evolution of MIR159/319 microRNA genes and their post-transcriptional regulatory link to siRNA pathways
© Li et al; licensee BioMed Central Ltd. 2011
Received: 13 August 2010
Accepted: 12 May 2011
Published: 12 May 2011
MicroRNAs (miRNAs) are prevalent and important endogenous gene regulators in eukaryotes. MiR159 and miR319 are highly conserved miRNAs essential for plant development and fertility. Despite high similarity in conservation pattern and mature miRNA sequences, miR159 and miR319 have distinct expression patterns, targets and functions. In addition, both MIR319 and MIR159 precursors produce multiple miRNAs in a phased loop-to-base manner. Thus, MIR159 and MIR319 appear to be related in origin and considerably diverged. However the phylogeny of MIR159 and MIR319 genes and why such unusual style of miRNA production has been conserved during evolution is not well understood.
We reconstructed the phylogeny of MIR159/319 genes and analyzed their mature miRNA expression. The inferred phylogeny suggests that the MIR159/319 genes may have formed at least ten extant early-branching clades through gene duplication and loss. A series of duplications occurred in the common ancestor of seed plants leading to the original split of flowering plant MIR159 and MIR319. The results also indicate that the expression of MIR159/319 is regulated at post-transcriptional level to switch on the expression of alternative miRNAs during development in a highly spatio-temporal specific manner, and to selectively respond to the disruption of defensive siRNA pathways. Such intra-stem-loop regulation appears diverged across the early-branching clades of MIR159/319 genes.
Our results support that the MIR159 and MIR319 genes evolve from a common ancestor, which is likely to be a phased stem-loop small RNA. Through duplication and loss of genes this miRNA gene family formed clades specific to moss, lycopods, gymnosperms and angiosperms including the two major clades of flowering plants containing the founding members of MIR319 and MIR159 genes in A.thaliana. Our analyses also suggest that some MIR159/319 have evolved into unusual miRNA genes that are regulated at post-transcriptional level to express multiple mature products with variable proportions under different circumstances. Moreover, our analyses reveal conserved regulatory link of MIR159/319 genes to siRNA pathway through post-transcriptional regulation.
MicroRNA is a class of approximately 21-nucleotide (nt) small RNA that regulate endogenous gene expression in eukaryotes ranging from single cellular green alga to mammals and flowering plants. Many biological processes including development, senescence, metabolism and stress responses are regulated by miRNAs [1, 2]. In plants, miRNAs are distinguished from other small RNAs in that they are excised precisely from specific positions on stem-loop precursors by DICER-LIKE 1 (DCL1) with the assistance of HYL1 and SERRATE [2, 3]. The mature miRNAs are 2'-o-methylated at 3' end by HEN1 and exported into cytoplasm. After being loaded into RNA Induced Silencing Complex (RISC) the miRNAs anneal to the complementary sites on target mRNAs to impose translational repression or/and cleavage of target mRNAs .
MiR159 and miR319 are highly conserved miRNAs that play important roles in plant growth, morphogenesis and reproduction . In Arabidopsis, the 21-nt mature miR159 and miR319 share 17 identical nucleotides. However, they have distinct target genes due to sequence specificity and different expression patterns . MiR159 restricts the expression of some MYB transcription factors, while miR319 targets a subset of TCP transcription factor genes [5–8]. Expression of miR159 is abundant and widespread over the whole plant, while miR319 is expressed at much lower level and confined to specific tissues and developmental stages . MiR159 can not induce mRNA cleavage of the miR319-targeted TCP transcription factors due to sequence specificity. Although the miR319 can also mediate the cleavage of MYB33 and MYB65 mRNAs, miR319 does not normally make significant contribution to the regulation of MYB because of its low and specialized expression . The interplay of miR159 and its target MYB is involved in the regulation of vegetative growth, flowering time, anther development, seed shape and germination [8, 10, 11]. In contrast to miR159, miR319 and corresponding targets regulate embryonic patterning, jasmonate synthesis, leaf morphogenesis and senescence [6, 12, 13]. A recent study also showed the regulatory role of miR319 in the development of petal and stamen .
A distinguished feature of the MIR159 and MIR319 genes conserved from moss to flowering plant is that their stem-loop precursors usually have elongated stem structure. A loop-proximal segment on the MIR159/319 stem-loop precursor outside of the miRNA and miRNA* is also conserved, albeit to a much weaker extent [6, 15, 16]. Recent studies show that the MIR159 and MIR319 precursors are processed from loop to base to liberate three phased miRNA duplexes [15, 17, 18]. The miR319 or miR159 can not be efficiently excised without correct processing of the loop-proximal miRNA duplex . According to this scenario, when miR159 or miR319 duplexes are produced the same amount of miRNAs derived from the other two duplexes must be generated. However, only the mature miR159 or miR319 are abundantly expressed while the alternative miRNAs are expressed at low levels, indicating that only miR159 or miR319 are incorporated into RISC and stabilized. Conservation of such an uncommon pattern in biogenesis during the long time evolution of land plants is inexplicable by the known function of miR159 and miR319. The underlying significance of such unusual style of maturation for the MIR159/319 genes remains unknown.
With similarities in sequence, conservation pattern and biogenesis, miR159 and miR319 might originate from a common ancestor. Conversely, differences in expression patterns, target genes and functions indicate that miR159 and miR319 are evolutionarily distinct groups. Although clear evidences are absent to support the common origin of miR159 and miR319 , they are categorized into one miRNA gene family in the miRbase and some other studies [7, 19]. Hypothesis has been proposed that the miR159 might evolve from miR319 because miR159 seems more specialized in the spectrum of targets . However, the phylogeny of the two hypothesized subfunctionalized groups remains unknown.
In order to study the evolution of MIR159/319 family in land plants, we collected and analyzed a broad set of MIR159/319 stem-loops from a wide range of plant species. Then, we reconstructed the phylogenetic tree of MIR159/319 from a structural alignment using RNA models for paired and unpaired nucleotides, thus revealing the evolutionary history of this miRNA gene family. Our analysis also suggests that post-transcriptional regulations of MIR159/319 expression have evolved during plant evolution, through which MIR159/319 miRNAs respond to defensive siRNA deficiency. This provides explanations for the conserved three-duplex dicing of the MIR159/319 precursors as a framework for an unknown product-tunable mechanism controlling output from their stem-loops. Taken together, our analysis revealed the evolutionary history of MIR159/319 genes, conservation and variability of post-transcriptional regulation across their extant clades during phylogeny.
Homology search for MIR159/319 genes
Classification of MIR159/319 genes according to conservation pattern
Classification and distribution of the five type MIR159/319 genes
Number of genes
Number of species
The remaining types are with degenerated ACRs (type 3 and type 4) or even lacking these two regions (type 2 and type 5). Type 2 stem-loops appear as in a shortened form without the elongated part probably due to deletion (Table 1 and Figure 2). Sequences of the loop-proximal regions are conserved in type 3 stem-loops (Table 1 and Figure 2). However the secondary structure is disrupted. The elongated stem is maintained in type 4 stem-loops, but sequence similarity in the region corresponding to ACR of type 1 is too weak to be identifiable (Table 1 and Figure 2). Type 5 stem-loops do not have identifiable ACR sequences and the elongated stem is interrupted by multi-loops (Table 1 and Figure 2). Conservation of the type 1 MIR159/319 genes across land plants indicates the origin of MIR159 and MIR319 from a long stem-loop. The type 2 to 5 MIR159/319 genes seems to be variants evolving from type 1. During evolution, the loop-proximal part on their precursor stem-loops may degenerate in sequence and/or structural conservation to different extent. Since the processing of miR159/319 genes depends on correct processing of the loop-proximal duplex , these MIR159/319 genes may lose their activity due to unsuccessful processing. For example, the processing of ath-MIR159c, which is a type 4 gene, is far less efficient than that of its paralogues . Likewise, mutations outside of the mature miR319 may lead to failure in the processing of MIR319 in B.oleracea . Therefore, some if not all type 2 to 5 MIR159/319 genes might be "pseudo-miRNA" genes unable to produce mature miRNAs.
Phylogeny of the MIR159/319 genes
During evolution the manner of base substitutions in structured RNA is different from that in protein-coding sequences. The selections for mutations in structured RNA mainly arise from base-pairings to maintain the RNA structures while in protein-coding sequences selection pressure disadvantages mutations in triplet codes that disrupt protein functions through alteration of amino-acids. Based on the structural alignment and consensus RNA secondary structure of the type 1 MIR159/319 genes, we reconstructed the phylogenetic tree using a combination of GTR and doublet models in a Bayesian approach . A well-supported consensus tree was obtained (Additional file 3). The three A.thaliana MIR319 paralogues are far more distantly related to each other than that of the ath-MIR159a and ath-MIR159b in the tree, consistent with the inference from the miRNA-carrying segmental duplications in Arabidopsis . The observation that ath-MIR159a and ath-MIR159b have very similar expression patterns and are functionally redundant while the three MIR319 paralogues have distinct expression patterns also supports the inferred relative distances among the paralogues of MIR159 and MIR319 in A.thaliana [10, 14]. In addition, the relationships of ath-MIR319b with its sense and antisense duplicates were correctly resolved in the tree (Additional file 3). These observations provide additional supports for the reliability of the inferred phylogenetic tree.
Duplication and loss of genes are deduced from the tree using the criteria that duplications increase the number of paralogues in one clade, while gene losses were inferred by the absence of an ancestral gene in descendant clades. Because the split of orthologous genes must predate or accompany the split of related species that carry them, the time of the duplications relative to the emergence of new plant clades can be deduced. Six duplications can be deduced from the condensed tree showing early-branching clades, one in the common ancestor of vascular plants and five in the common ancestor of seed plants (Figure 4). The tree shows that the split of the MIR319 and MIR159 clades in flowering plants and gymnosperms originated from a series of duplications occurring in the common ancestor of seed plants (Figure 4). This suggests that the MIR319 and MIR159 genes classified by the founding members identified in A.thaliana originated from a common ancestor, which may have emerged after the split of moss and lycopod MIR159/319 clades. The two angiosperm MIR319 clades branched earlier than the angiosperm MIR159 clades. Subsequent duplications led to the formation of gymnosperm clades and angiosperm MIR159 clades (Figure 4). Six possible losses in the major clades of MIR159/319 can be deduced from the condensed tree, one in lycopod and five in gymnosperms and angiosperms (Figure 4). Because our analyses do not include all the species in a given clade and entire genomes are only available for a small number of species, the deduced gene losses are incomprehensive and indicate only the possible inclination of gene evolution in the species examined. However, deduced losses are of high confidence in the species for which genome-wide scan of MIR159/319 gene have been performed (see methods). Of the six possible losses, four are supported by the absence of ancestral MIR159/319 genes in one or more species where the entire genomes were examined (Figure 4). The other two losses are deduced to have occurred in gymnosperms whose whole genome sequences are currently not available (Figure 4). Therefore, these two losses and loss of gymnosperm MIR159/319 in 'early branching eudicot MIR319' need to be verified in the future using the whole genomes of gymnosperm species (Figure 4). Considering the current data, these gene losses may have occurred in an interleaved manner that descendant species mutually lost the equivalent of extant genes in other descendant species after gene duplications in their common ancestor. It appears that clade-specific groups of MIR159/319 genes may have formed through the dynamics of gene duplications in ancestors and gene losses in descendants as opposed to the divergence of one ancestral gene into orthologues. The duplications and subsequent losses that occurred before the split of monocots and eudicots are more frequent in the MIR319 clade than the MIR159, indicating that the two clades might have evolved in different ways (Figure 4 and Additional file 3).
Shift in the proportion of mature MIR159/319 products during ontogeny
We also observed the switching on of the MIR159 or MIR319 ACR miRNA expression in monocots. In rice, analyses using high throughput small RNA sequencing data deposited under [GEO:GSE14462] suggest that osa-MIR159a mainly produces ACR3 miRNAs in the flag leaves (Figure 5G), which is an organ important for grain filling and protection of the spike from being eroded by pathogens. Sequencing data from [GEO:GSE23217] suggest that ACR3 miRNA is the major product of osa-MIR319 in rice seedlings  (Figure 5H). In maize, a shift in the proportion of mature MIR319 products was observed in the tissues of shoots and roots using the data from [GEO:GSE15286]  (Figure 5I). Such a transition is not observed for maize MIR159 (Additional file 5H). Another two examples of non-canonical expression of monocot MIR159/319 come from CSPSR data  (Figure 5J and 5K). In Hordeum vulgare plants, expression level of MIR159 ACR3 miRNA is higher than miR159 in leaves but not in inflorescence tissues, indicating shift of mature product proportions in the two tissues (Figure 5J). In flower and young panicle of Sorghum bicolor, ACR3 miRNAs are more abundant than miR319 (Figure 5K). These results indicate that both in eudicots and monocots, the major products of MIR159 or MIR319 genes shift to ACR miRNAs in a highly spatio-temporal specific manner during plant development. This is likely to occur through post-transcriptional mechanisms, because the proportion of the multiple mature products can not be changed solely via transcriptional regulations. In some cases, shift in strand selection also occurred. When sequencing tags can not be unambiguously assigned to a unique MIR159/319 stem-loop, such as maize MIR319 genes in Figure 5I, related paralogues are analyzed collectively. In this case, the observed shift in mature product proportion could possibly be attributed to down- or up-regulation at transcriptional level of different paralogues with fixed mature product proportions. However, there must be a post-transcriptional mechanism to distinguish different paralogues so that they produce multiple mature miRNAs in different proportions.
For the gymnosperms and earlier-branching clades, lycopod and moss, we also analyzed the expressions of the MIR159/319 genes using small RNA sequencing data from CSPSR and [GEO:GSE5103, GEO:GSE7320, GEO:GSE12468] [15, 29, 34]. In Cycas rumphii leaves, the three phased miRNAs at 3' arm, ACR3, sp4 and miR159 are expressed at significant levels, among which ACR3 is the most abundant (Figure 5L). In the lycopod S.moellendorffii, the predominant products of both smo-MIR159 and smo-MIR319 are ACR5 miRNAs, consistent with the distance in the inferred phylogenetic tree (Figure 4, Figure 5M and Additional file 5J). The proportion of mature products from the moss MIR319 differs from all other clades in that the ACR miRNAs are rarely expressed while the sp2 and sp4 miRNAs are expressed at significantly higher level (Figure 5N and Additional file 5K). Because the sp2 miRNA of ppt-MIR319d is abundantly expressed, this miRNA is annotated as ppt-miR319d*.2 [15, 35, 36]. This is consistent with our analysis. In 10-day-old moss protonemata the miR319* are the most abundant products (Figure 5N). The distinct proportion of mature products from the MIR159/319 genes in different clades suggests that the post-transcriptional control of MIR159/319 gene output may have diverged during evolution accounting for the variability of mature miRNA expression profiles of this family.
Selective up-regulation of mature MIR159/319 miRNAs in siRNA-deficient mutants
From the evolutionary perspective, miRNAs are generally classified into highly conserved ancient miRNAs and weakly conserved species- or clade-specific miRNAs . The phylogeny of highly conserved miRNA genes is largely unknown including the ancient miRNA gene family -- MIR319 and MIR159, which play important roles in plant development . The miR319 and miR159 were identified in independent studies using different methods [6, 24]. Homologues identified in following studies are classified into two groups according to sequence similarity to the miR319 and miR159 founding members . Although miR159 and miR319 are seemingly related in evolution, because of their similarity in mature miRNA sequence, secondary structure, conservation pattern and biogenesis, their origin is still unclear [5, 17]. Evolution of MIR319 genes in closely related Brassicaceae species has been studied revealing that the stem-loop and an upstream element of MIR319 are recalcitrant to fast mutations while sequences of other parts are highly variable . We aligned a large number of land plant MIR159/319 stem-loops and reconstructed the phylogenetic tree. Our results support a common origin of MIR159 and MIR319 from two aspects. First, another duplex outside miR159/319 are highly conserved in most MIR159 and MIR319 stem-loops across land plants (Figure 2, Figure 3, Table 1 and Additional file 2). Moreover, partitions defined by the conserved regions are consistent with the phasing of mature MIR159/319 miRNAs from moss to flowering plants [17, 18] (Figure 3 and Additional file 4). These observations suggest that MIR159 and MIR319 originated from a common phased stem-loop RNA similar to those discovered in the green alga Chlamydomonas reinhardtii and rice recently [42–44]. Second, the inferred phylogenetic tree shows that the MIR159 and MIR319 founding members in A.thaliana evolve from two early-branching clades specific to flowering plants, which derived from the common ancestor of seed plant MIR159/319 (Figure 4 and Additional file 3). The tree also shows that the angiosperm MIR319 major clade branched earlier than the angiosperm MIR159 clade (Figure 4). This is consistent with the finding that miR159 is more specialized in target spectrum than the miR319 in A.thaliana .
The inferred phylogenetic tree reveals the evolutionary history of extant MIR159/319 branches, which might be dated back to the first colonization of land plants about 450 Ma (million years ago). It can be inferred that extensive duplications of MIR159/319 genes might have occurred in the common ancestor of gymnosperms and angiosperms about 370 Ma, in the period of late Devonian to early Carboniferous. These duplications initialized the split of MIR159 and MIR319 in flowering plants. It has been shown that major body-plan innovations during animal phylogeny are concomitant with miRNA repertoire expansions . Hypothesis has been proposed that expansion of miRNAs may contribute to the improvement of body-plan complexity [45, 46]. Whether similar scenario is also the case in plants is unknown. From the phylogenetic tree, it can be inferred that an expansion of MIR159/319 family genes might be commitment with the evolution of seed-bearing plants (Figure 4). Recent studies have shown that both miR159 and miR319 play important roles in reproductive growth of flowering plants [8, 14], indicating their functional relatedness to the evolution of seed bearing in plants. It has been shown in both plants and animals that miRNA repertoire encoded in the genome is in the dynamics of genesis and loss of miRNA genes during evolution [39, 47, 48]. Some new born miRNA genes seem to be targetless and evolutionarily transient, likely to serve as resources for the selection of functional miRNAs [39, 47]. Similarly, from the phylogenetic tree of MIR159/319 family, we found this ancient miRNA family underwent dynamics of duplication and loss, through which clade or species-specific miRNA gene subgroups have formed. Specifically, during evolution clade-specific subgroups formed in the manner that specific members of the expanded gene family maintained in a given descendant clade were lost in sister clades. Moreover, duplicated gene members maintained in one clade of plants, e.g. miR159 and miR319 in flowering plants, tend to diverge into subfunctionalized groups improving the complexity of miRNA regulation. Similarly, the three members of Arabidopsis miR319, which arose from more recent duplications, also diverged into genes with distinct expression patterns and possibly different functions .
The ACR miRNAs are much more variable than the miR159 or miR319 (Figure 3), indicating different selection pressure imposed on the two miRNAs on the same stem-loop. Moreover, a small number of MIR159/319 genes (type 2 to 5) lost conservation of the loop-proximal part corresponding to the ACRs of type 1. Since the correct processing of the loop-proximal miRNA duplex is the prerequisite for the dicing of miR159 or miR319 , the ACR-lacking MIR159/319 might become pseudo genes unable to produce mature miRNAs like the ath-MIR159c, which can not be efficiently processed when transformed into plants . Alternatively, the expression of non-type-1 MIR159/319 genes might be highly spatio-temporal specific or inducible under specific conditions. Since the non-type 1 MIR159/319 genes are rare and we did not find exact evidence that any species only keep non-type 1 MIR159/319 in their genomes, the non-type 1 MIR159/319 genes might derive from type 1.
With regulatory roles in gene expression, expression of miRNAs is also regulated delicately. Many miRNAs are expressed in specific population of cells dedicated to certain functions during development [1, 2]. In animals, growing evidences suggest that miRNA expression is not only regulated at transcriptional level but also at post-transcriptional level. The processing of specific pri-miRNAs by the microprocessor is regulated via diverse mechanisms in various biological processes . Recently, selective stabilization of miRNAs by 3'-terminal adenylation has been reported revealing post-dicing regulation of miRNA abundance by synthesis and degradation homeostasis . However, plant genomes do not possess genes encoding protein factors equivalent to Drosha in animals, which is the RNase III enzyme catalyzing pri- to pre-miRNA cleavage . Question arises whether post-transcriptional mechanisms have also evolved in plants to control the expression of specific miRNAs. Variegated dependencies on miRNA pathway components of miRNA accumulations indicate that there might be miRNA-specific regulations at post-transcriptional steps in plant miRNA biogenesis and stability. For example, it has been reported that null or hypomorphic alleles of DCL1, AGO1, SERRATE, DRB1, ABH1 and CBP20 have different impact on the accumulation of distinct miRNA subsets [51–54]. A family of exoribonucleases responsible for the degradation of plant miRNAs is identified, indicating that mature miRNAs can also be regulated through degradation in plants . We found spatio-temporal specific shifts in the proportion of miRNA products within MIR159/319 stem-loop precursors during development (Figure 5, Figure 6 and Additional file 5). Expression analysis also showed selective up-regulation of mature MIR159/319 miRNAs in siRNA-deficient mutants of moss and flowering plants (Figure 7 and Additional file 6). These observations provide evidences for the evolution of miRNA-specific post-transcriptional regulatory mechanism in plants.
Across the major clades of MIR159/319 genes (Figure 4), the post-transcriptional regulation of miRNA production seems to be diverged. First, general proportions of MIR159/319 mature miRNAs are distinct for the early-branching clades. Moss MIR159/319 genes have pronounced expression in sp2 and sp4 miRNAs. Lycopod MIR159/319 genes produce ACR miRNAs in higher levels than miR159/319 in normal state. The gymnosperm Cycas rumphii MIR159 produce three phased miRNAs at comparable levels from 3' arm of the stem-loop precursor in leaves (Figure 5L). In most cases, angiosperm MIR319 and MIR159 generate miR319 or miR159 as major products (Figure 4, Figure 5 and Additional file 4). Second, genes from angiosperm major MIR319 clade and major MIR159 clade can have different proportions of mature miRNAs in specific organs of the same plant (Figure 4, Figure 5 and Additional file 5). Third, genes from different early-branching MIR159/319 clades have different responsiveness to the siRNA deficiency in the mature miRNA expressions (Figure 4 and Figure 7). In addition, paralogues in one early-branching MIR159/319 clade also show differences in tissue-specific proportion of mature miRNAs and responsiveness to siRNA deficiency. For example, the proportion of ath-MIR319a ACR3 miRNA is much higher in 5-week-old inflorescence than that of ath-MIR319b (Figure 5A). The paralogues of MIR319 also showed distinct expressions in A.lyrata (Figure 5B and Additional file 5C). In siRNA deficiency mutants, ACR3 miRNA of ath-MIR319a is up-regulated to much higher level than that of ath-MIR319b (Figure 7A and B). Such differences in their expressions are consistent with the inferred phylogenetic distance and different tissue specific expression patterns of ath-MIR319a and ath-MIR319b reported recently  (Additional file 3). These observations suggest that the post-transcriptional regulation of MIR159/319 may have specialized during evolution and can serve as another feature that could distinguish their subgroups.
In addition to the small RNAs dedicated to endogenous gene regulation in normal developmental and physiological processes, e.g. miRNA, tasiRNAs, various siRNA pathways have evolved in plants and animals to harness transposons, antagonize infection of pathogens and respond to biotic and abiotic stresses [56, 57]. Epigenetic inactivation of transposons and exogenous DNAs mediated by approximately 24-nt siRNAs require the action of DCL3 and RDR2 for their biogenesis and AGO4 in effecter complex [56, 57]. DCL2, DCL4 and RDR6 are involved in the production of 21~22-nt siRNAs from virus and A.tumefaciens T-DNA, which mediate post-transcriptional silencing of pathogen-derived genes [56, 57]. Selective up-regulation of MIR159/319 ACR miRNAs in rdr2 and dcl2 dcl3 dcl4 triple mutants indicates that intrinsic link between miRNA and the defensive siRNA pathways has evolved in plants (Figure 7). Mutation of the gene encoding RDR2 can cause global changes in small RNA transcriptome and gene expressions including reduction of 24-nt heterochromatic siRNA, enrichment of miRNA and tasiRNAs, reactivation of specific type of transposons and altered expression of genes in chromatin modification pathways [41, 58, 59]. Some of the miRNAs are up-regulated much higher above the overall enrichment of miRNAs in rdr2 mutants, such as the miR163 in A.thaliana [41, 58]. The up-regulations of MIR159/319 miRNAs in rdr2 or dcl234 triple mutants are very similar to these cases. However, the up-regulations of MIR159/319 miRNAs are of strong selectivity upon the small RNA products from the same stem-loop suggesting that it is unlikely to be a non-specific process (Figure 7). It was speculated that the different accumulation levels of miRNAs in wild type and rdr2 mutants are caused by a secondary level of control by a siRNA-mediated pathway, which depends on RDR2/mop1 . Likewise, the ACR miRNA accumulation of MIR319 genes might be downstream of RDR2-dependent siRNA regulation. Since only ACR miRNAs are induced to high levels while increment of miR319 is similar to the overall enrichment of miRNAs in rdr2 or dcl234 triple mutants, a post-transcriptional mechanism may distinguish the two miRNA duplexes on the same stem-loop and link the siRNA pathway. Because the selective up-regulation of MIR159/319 genes is observed in both monocot and eudicots, the regulatory link between MIR159/319 genes and siRNA pathway may have an ancient origin. It has been found in Physcomitrella patens that silencing of miRNA targets via DNA methylation occurs in PpDCLb mutants and ABA treated wild type plants . This suggests that silencing mechanisms of the heterochromatic siRNA pathway can also be exploited by the miRNA pathway. Our findings revealed another layer of interactions between the two small RNA pathways that production of miRNA can be affected by defects in heterochromatic siRNA pathway.
This uncharacterized mechanism controlling mature product proportion of MIR159/319 genes is distinguished from those previously found in animals in that multiple mature products from a single miRNA stem-loop precursor can be selectively expressed in different situations. PARE data suggest that speed for dicing at the steps during phased processing might be regulated, provided that detection of more intermediate products can serve as the signal for higher dicing speed for the corresponding step of slicing. Accumulations of miR159, miR168 and miR165 are insensitive to decreased DCL1 activity unlike other miRNAs such as miR173, indicative of alternative dicer activity for the processing of these miRNAs . Therefore, it is possible that only processing by alternative factors yields abundant ACR miRNAs but less miR159/319 when alternative activity overrides the canonical. Similar scenarios are possible for other factors involved in miRNA biogenesis and activity, dysfunction of which affect miR159 in a mild way. For example, the ABH1 is required for the processing of some pri-miRNAs but not MIR159 [52, 53]. Alternatively, since the theoretical ratio of ACR miRNA to miR159 or miR319 should be one, the variable proportion of MIR159/319 mature products could be resulted from regulation of mature miRNA stability. Probably, which of the MIR159/319 mature miRNAs are incorporated into the RISC and stabilized is tunable. Since the ACR3 miRNA has been shown to associate with AGO4 [17, 61], incorporation of ACR3 miRNA into AGO4 might be tremendously increased in the rdr2 or dcl234 mutants, of which the AGO4 associated siRNAs are depleted [41, 58, 59]. Increased AGO1 up-regulates miR159 accumulation indicating excessive miR159 might be subjected to active degradation . Likewise, it is possible that the ACR miRNAs might be directly degraded in normal state but stabilized under specific circumstances. This kind of regulation has been reported for a class of 21-nt endogenous siRNAs derived from mRNAs, which accumulated only in the absence of EIN5, a 5'-3' exoribonuclease . Possibly, degradation of processing intermediates could be another way through which maturation of miRNA is regulated. Some miRNA processing intermediates are up-regulated in exosome-deficient mutants, including the loop fragment of MIR159 . Since MIR159/319 precursors are processed in a phased manner, differential exosome-mediated degradation of the intermediates might result in variable product proportions. A combined regulation at multiple steps of miRNA metabolic cycle can not be ruled out. In addition to the unexpected up-regulation of ACR miRNAs, we also observed some cases of higher expression of miR319* than miR319, namely a shift in strand selection (Figure 5A and 5B). The change in strand bias of miR159 and miR159* have been observed in drb1 mutant plants, which direct strand selection of miRNA duplexes . Possibly, the spatial-temporal specific selection of strands in miRNA duplexes might be resulted from incongruent expressions of the miRNA and DRB1.
The switching on of the ACR miRNAs in specific tissues and their responsiveness to siRNA deficiency suggests that these miRNAs might be functional. There is only one experimentally verified example that osa-MIR159a ACR3 miRNA induces site-specific cleavage of an mRNA encoding a GT-2-like transcription factor in rice . However, it can not be excluded that the MIR159/319 ACR miRNAs may exert their function over their targets through other mechanisms such as translational repression or RNA mediated DNA methylation [2, 60]. In light of the intra-stem-loop regulation, the phased loop-to-base processing of miR159/319 stem-loop precursors may have evolved as framework for the tunable mature miRNA production.
In this study we showed that MIR159/319 genes with two conserved duplexes in their stem-loops distribute throughout land plants and constitute the most plentiful type classified by conservation pattern. Together with phylogenetic reconstruction and mature miRNA expression analysis, evidences support that the MIR159 and MIR319 have a common origin. Evolution history of this gene family was revealed that the two clades containing the MIR159 and MIR319 founding members identified in A.thaliana are specific to flowering plants and originated from duplications occurring before the split of gymnosperm and angiosperms. Interestingly, post-transcriptional regulation of multiple mature miRNAs from some genes of MIR159/319 family has evolved, through which proportion of mature MIR159/319 miRNAs can be changed in a spatio-temporal specific manner or in response to siRNA pathway deficiency. These findings provide explanation for the conserved phased loop-to-base processing of MIR159/319 genes as framework for post-transcriptionally regulated expressions of their mature miRNAs. Taken together, our study provides insights on the evolution of MIR159/319 genes and an unexplored feature of them that the expressions of multiple forms of mature products are post-transcriptionally regulated.
Homology search of miRNA genes
Known MIR159/319 sequences were collected from the miRbase Version 11.0 , homologues in EST were identified by Jones-Rhoades and Bartel , and the Brassicaceae MIR319 sequences were identified by Warthmann et al. . In one line, SSEARCH searches were performed against the TIGR Plant Transcript Assembly (http://plantta.tigr.org) and the 11 plant genomes listed in Additional File 7 using ath-miR159a or ath-miR319a as queries . In the other line, all of the mature miRNA sequences of the MIR159/319 family available in miRbase Version 11.0 were used to BLAST search against the nr, gss, htgs and est databases in the NCBI GenBank . The flanking sequences (350-nt or as long as possible) on either side were retrieved for each resulting High-scoring Sequence Pairs (HSP) from the searches. The RNA secondary structures were predicted using RNAfold . The miRNA-like stem-loop structures were selected using modified miRcheck rules . Briefly, 1) the paired bases for miR or miR* (miR: miR159/319 homologous sequence; miR*: antisense sequence in the stem-loop structure corresponding to miR) must be located within one arm, either the 5' or the 3'; 2) there must be no more than 6 unpaired bases in the miR; 3) no more than 3 bulges in the miR; 4) no more than 3-nt differences between the lengths of miR and miR*; 5) no more than 3 asymmetrically unpaired bases in the miR; 6) no more than 3 contiguous unpaired bases in the miR; and 7) the distance between miR and miR* must be no less than 5 and no more than 300 nucleotides. The filtrates were inspected manually to exclude the loose and unstable structures. Finally, if the miR sequence has no more than 3 mismatches or indels compared to a MIR159/319 mature miRNA in the miRbase, the stem-loop was considered a homologue. If the miR sequence contains 4 mismatches or indels, the stem-loop was considered a candidate homologue. Homologous sequences contain more mismatches were not considered. The naming of the new homologues or candidates, 159 or 319, followed the name of the MIR159/319 gene in the miRbase that had the smallest edit distance in the mature miRNA sequence. The abbreviations for the species names are: afo: Aquilegia formosa; ace: Arabidopsis cebennensis; aha: Arabidopsis halleri; aly: Arabidopsis lyrata subsp. Petraea; ath: Arabidopsis thaliana; ahy: Arachis hypogaea; bvl: Barbarea vulgaris; bvu: Beta vulgaris; bst: Boechera stricta; bna: Brassica napus; bol: Brassica oleracea; bra: Brassica rapa subsp. Pekinensis; cbp: Capsella bursa-pastoris; cru: Capsella rubella; cso: Centaurea solstitialis; cch: Cheiranthus cheiri/Erysimum cheiri; ccl: Citrus clementina; csi: Citrus sinensis; csu: Citrus sunki; cca: Coffea canephora; cor: Conringia orientalis; cme: Cucumis melo subsp. Melo; crp: Cycas rumphii; dso: Descurainia sophia; dka: Diospyros kaki; ees: Euphorbia esula; far: Festuca arundinacea; fve: Fragaria vesca; gma: Glycine max; htu: Helianthus tuberosus; hvu: Hordeum vulgare; ini: Ipomoea nil; lsl: Lactuca saligna; lst: Lactuca sativa; lse: Lactuca serriola; ltu: Liriodendron tulipifera; lja: Lotus japonicus; mma: Malcolmia maritima; mdo: Malus × domestica; mes: Manihot esculenta; mtr: Medicago truncatula; Nof: Nasturtium officinale; nta: Nicotiana tabacum; oal: Oryza alta; oco: Oryza coarctata; ogl: Oryza glaberrima; omi: Oryza minuta; oni: Oryza nivara; opu: Oryza punctata; ori: Oryza ridleyi; oru: Oryza rufipogon; osa: Oryza sativa; osg: Oryza sativa × Oryza glaberrima; pvi: Panicum virgatum; pgl: Pennisetum glaucum; pvu: Phaseolus vulgaris; ppt: Physcomitrella patens; pga: Picea glauca; psi: Picea sitchensis; pta: Pinus taeda; ptf: Poncirus trifoliata; pde: Populus deltoides; pni: Populus nigra; ptr: Populus tremula; ptt: Populus tremula × Populus tremuloides; ptc: Populus trichocarpa; rra: Raphanus raphanistrum subsp. Raphanistrum; rco: Ricinus communis; rin: Rorippa indica; sof: Saccharum officinarum; smo: Selaginella moellendorffii; svi: Sibara virginica; sly: Solanum lycopersicum; sbi: Sorghum bicolor; tar: Thlaspi arvense; tae: Triticum aestivum; vun: Vigna unguiculata; vvi: Vitis vinifera; zma: Zea mays. This three-letter code shown in all lower case is used for sequences from miRbase , all upper case letters are used for sequences from Jones-Rhoades and Bartel , and others are written with only the first letter capitalized.
The taxonomic tree was plotted using NCBI taxonomy (http://www.ncbi.nlm.nih.gov/Taxonomy/CommonTree/wwwcmt.cgi) and drawn by NJplot .
Structural alignment and phylogenetic inference
All of the stem-loop sequences of the MIR159 family were aligned using the software T-coffee version 6.06 . Sequences that can't be aligned correctly in the conserved regions outside of miR159/319 and miR159/319* were excluded in subsequent analysis. The remaining sequences were aligned again. Secondary structures were produced for each sequence in the alignment using RNAfold . The software 4SALE was then used to correct the alignment manually by considering the agreement in secondary structure . After removing the columns that contain more than a 75% gap, the consensus sequences and structures were produced using ConStruct, RNAalifold, alidot and SARSE [70–73]. Structure logos were produced using slogo . The partitions of the well-aligned stem-loops were divided empirically by assigning approximately 21-nt for the well-conserved regions. The phylogenetic tree was inferred using MrBayes version 3.1 with the GTR model for unpaired regions and the Doublet model for paired regions with 4 discrete gamma distribution categories . Two chains were run for 20 million generations. Trees were sampled every 100 generations. The consensus tree was summarized using a 50 majority rule from 300,000 trees after removing the burn-in. Only one sequence was retained for identical sequences in the alignment for phylogenetic inference. These representative taxa were replaced with nodes of the original redundant taxa, for which posterior probabilities were set to 100 and branch lengths were set to 0.001 for viewing. The phylogenetic tree was edited using Dendroscope and MEGA4 [75, 76]. Gene duplications were deduced from the tree using the criteria that duplications increase the number of paralogues in one clade, while gene losses were inferred by the absence of an ancestral gene in descendant clades.
Re-analysis of existing small RNA sequencing data
Small RNA sequences were mapped on the stem-loop sequences using a customized perl script. Only perfectly matched small RNAs were considered. The small RNA sequencing resources are listed in Additional File 8. Sequencing abundances were normalized to the average sequencing scale for all samples in a given database unless the values had already been normalized using the original data generators. The expression value for a specific partition in the alignment was calculated by summing all the normalized sequencing frequencies of unique small RNAs that had 80% or more of their nucleotides overlap with that partition. The small RNAs that could not be assigned to any partition are categorized as "leak" to denote that these mature miRNAs were inconsistent with the phased processing. The mutant-to-wild type ratios were calculated using the maximum possible values, which included reads that also matched the closely related paralogues when an individual miRNA gene was studied. For [GEO:GSE12468] of P. patens , the mutant-to-wild type ratio was calculated by using the mean of two replicates for each of the wild type and mutant values.
PARE (Parallel Analysis of RNA Ends) data analysis
One part of PARE sequence tags was downloaded from the AtPARE database (http://mpss.udel.edu/at_pare). The other part was obtained from [GEO:GSE11007] . Sequencing abundances for all samples were normalized to TP10M (transcript per 10 million reads). The tags that represent uncapped 5' RNA ends were mapped to the Arabidopsis thaliana MIR159 and MIR319 genes, stem-loops or an approximately 1 kb sequence with the stem-loop at the centre. The pri-MIR319a is from [GenBank:AY922324] . The frequency of the cleavage between two nucleotides was calculated by summing the sequencing abundances of all the sequence tags whose 5' first nucleotide situated exactly at the second nucleotide. Total cleavages were calculated by summing the values of all the samples for a cleavage site. The maximum possible frequencies that included tags that also matched closely related homologues were used in Figure 6. RNA secondary structures were drawn using XRNA (http://rna.ucsc.edu/rnacenter/xrna/xrna.html).
List of Abbreviations
alternative conserved region
RNA induced silencing complex
We thank authors of the small RNA sequencing and PARE data used in this study including Drs. Blake Meyers and Pamela Green in University of Delaware, James Carrington in Oregon State University, Lila Vodkin in University of Illinois, Hai Peng in Jianghan University, Tian Li in tsinghua university, Xingwang Deng in Yale University, David Bartel in White Head Institute, Michael Axtell in Pennsylvania State University. We also thank suggestions of the two anonymous reviewers. This work was supported by the grants from the National Key Basic Research and Development Program (No. 2011CB811304), National Basic Research Program of China (No. 2011CBA01105) and the Chinese Academy of Sciences (KSCX2-YW-R-096).
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.View ArticlePubMedGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.View ArticlePubMedGoogle Scholar
- Palatnik JF, Wollmann H, Schommer C, Schwab R, Boisbouvier J, Rodriguez R, Warthmann N, Allen E, Dezulian T, Huson D, Carrington JC, Weigel D: Sequence and expression differences underlie functional specialization of arabidopsis microRNAs miR159 and miR319. Dev Cell. 2007, 13 (1): 115-125. 10.1016/j.devcel.2007.04.012.View ArticlePubMedGoogle Scholar
- Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425 (6955): 257-263. 10.1038/nature01958.View ArticlePubMedGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.View ArticlePubMedGoogle Scholar
- Millar AA, Gubler F: The Arabidopsis GAMYB-like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell. 2005, 17 (3): 705-721. 10.1105/tpc.104.027920.View ArticlePubMedPubMed CentralGoogle Scholar
- Axtell MJ, Bartel DP: Antiquity of microRNAs and their targets in land plants. Plant Cell. 2005, 17 (6): 1658-1673. 10.1105/tpc.105.032185.View ArticlePubMedPubMed CentralGoogle Scholar
- Allen RS, Li J, Stahle MI, Dubroue A, Gubler F, Millar AA: Genetic analysis reveals functional redundancy and the major target genes of the Arabidopsis miR159 family. Proc Natl Acad Sci USA. 2007, 104 (41): 16371-16376. 10.1073/pnas.0707653104.View ArticlePubMedPubMed CentralGoogle Scholar
- Reyes JL, Chua NH: ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 2007, 49 (4): 592-606. 10.1111/j.1365-313X.2006.02980.x.View ArticlePubMedGoogle Scholar
- Ori N, Cohen AR, Etzioni A, Brand A, Yanai O, Shleizer S, Menda N, Amsellem Z, Efroni I, Pekker I, Alvarez JP, Blum E, Zamir D, Eshed Y: Regulation of LANCEOLATE by miR319 is required for compound-leaf development in tomato. Nat Genet. 2007, 39 (6): 787-791. 10.1038/ng2036.View ArticlePubMedGoogle Scholar
- Schommer C, Palatnik JF, Aggarwal P, Chetelat A, Cubas P, Farmer EE, Nath U, Weigel D: Control of jasmonate biosynthesis and senescence by miR319 targets. PLoS Biol. 2008, 6 (9): e230-10.1371/journal.pbio.0060230.View ArticlePubMedPubMed CentralGoogle Scholar
- Nag A, King S, Jack T: miR319a targeting of TCP4 is critical for petal growth and development in Arabidopsis. Proc Natl Acad Sci USA. 2009, 106 (52): 22534-22539. 10.1073/pnas.0908718106.View ArticlePubMedPubMed CentralGoogle Scholar
- Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants. Plant Cell. 2007, 19 (6): 1750-1769. 10.1105/tpc.107.051706.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Y, Li W, Jin YX: Computational identification of novel family members of microRNA genes in Arabidopsis thaliana and Oryza sativa. Acta Biochim Biophys Sin (Shanghai). 2005, 37 (2): 75-87. 10.1093/abbs/37.2.75.View ArticleGoogle Scholar
- Bologna NG, Mateos JL, Bresso EG, Palatnik JF: A loop-to-base processing mechanism underlies the biogenesis of plant microRNAs miR319 and miR159. EMBO J. 2009, 28 (23): 3646-3656. 10.1038/emboj.2009.292.View ArticlePubMedPubMed CentralGoogle Scholar
- Addo-Quaye C, Snyder JA, Park YB, Li YF, Sunkar R, Axtell MJ: Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome. RNA. 2009, 15 (12): 2112-2121. 10.1261/rna.1774909.View ArticlePubMedPubMed CentralGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, D140-144. 34 Database
- Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species. Bmc Plant Biology. 2008, 8:Google Scholar
- Warthmann N, Das S, Lanz C, Weigel D: Comparative analysis of the MIR319a microRNA locus in Arabidopsis and related Brassicaceae. Mol Biol Evol. 2008, 25 (5): 892-902. 10.1093/molbev/msn029.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Maher C, Stein L, Ware D: Evolution of Arabidopsis microRNA families through duplication events. Genome Res. 2006, 16 (4): 510-519. 10.1101/gr.4680506.View ArticlePubMedPubMed CentralGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes Dev. 2002, 16 (13): 1616-1626. 10.1101/gad.1004402.View ArticlePubMedPubMed CentralGoogle Scholar
- German MA, Pillay M, Jeong DH, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008, 26 (8): 941-946. 10.1038/nbt1417.View ArticlePubMedGoogle Scholar
- Ma Z, Coruh C, Axtell MJ: Arabidopsis lyrata small RNAs: transient MIRNA and small interfering RNA loci within the Arabidopsis genus. Plant Cell. 2010, 22 (4): 1090-1103. 10.1105/tpc.110.073882.View ArticlePubMedPubMed CentralGoogle Scholar
- Fahlgren N, Jogdeo S, Kasschau KD, Sullivan CM, Chapman EJ, Laubinger S, Smith LM, Dasenko M, Givan SA, Weigel D, Carrington JC: MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana. Plant Cell. 2010, 22 (4): 1074-1089. 10.1105/tpc.110.073999.View ArticlePubMedPubMed CentralGoogle Scholar
- Tuteja JH, Zabala G, Varala K, Hudson M, Vodkin LO: Endogenous, tissue-specific short interfering RNAs silence the chalcone synthase gene family in glycine max seed coats. Plant Cell. 2009, 21 (10): 3063-3077. 10.1105/tpc.109.069856.View ArticlePubMedPubMed CentralGoogle Scholar
- Comparative Sequencing of Plant Small RNAs. [http://smallrna.udel.edu/index.php]
- Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ: Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008, 18 (10): 758-762. 10.1016/j.cub.2008.04.042.View ArticlePubMedPubMed CentralGoogle Scholar
- Chekanova JA, Gregory BD, Reverdatto SV, Chen H, Kumar R, Hooker T, Yazaki J, Li P, Skiba N, Peng Q, Alonso J, Brukhin V, Grossniklaus U, Ecker JR, Belostotsky DA: Genome-wide high-resolution mapping of exosome substrates reveals hidden features in the Arabidopsis transcriptome. Cell. 2007, 131 (7): 1340-1353. 10.1016/j.cell.2007.10.056.View ArticlePubMedGoogle Scholar
- Li T, Li H, Zhang YX, Liu JY: Identification and analysis of seven H2O2-responsive miRNAs and 32 new miRNAs in the seedlings of rice (Oryza sativa L. ssp. indica). Nucleic Acids Res. 2011, 39 (7): 2821-2833. 10.1093/nar/gkq1047.View ArticlePubMedGoogle Scholar
- Wang X, Elling AA, Li X, Li N, Peng Z, He G, Sun H, Qi Y, Liu XS, Deng XW: Genome-wide and organ-specific landscapes of epigenetic modifications and their relationships to mRNA and small RNA transcriptomes in maize. Plant Cell. 2009, 21 (4): 1053-1069. 10.1105/tpc.109.065714.View ArticlePubMedPubMed CentralGoogle Scholar
- Cho SH, Addo-Quaye C, Coruh C, Arif MA, Ma Z, Frank W, Axtell MJ: Physcomitrella patens DCL3 is required for 22-24 nt siRNA accumulation, suppression of retrotransposon-derived transcripts, and normal development. PLoS Genet. 2008, 4 (12): e1000314-10.1371/journal.pgen.1000314.View ArticlePubMedPubMed CentralGoogle Scholar
- Arazi T, Talmor-Neiman M, Stav R, Riese M, Huijser P, Baulcombe DC: Cloning and characterization of micro-RNAs from moss. Plant J. 2005, 43 (6): 837-848. 10.1111/j.1365-313X.2005.02499.x.View ArticlePubMedGoogle Scholar
- Fattash I, Voss B, Reski R, Hess WR, Frank W: Evidence for the rapid expansion of microRNA-mediated regulation in early land plant evolution. BMC Plant Biol. 2007, 7: 13-10.1186/1471-2229-7-13.View ArticlePubMedPubMed CentralGoogle Scholar
- Xie Z, Johansen LK, Gustafson AM, Kasschau KD, Lellis AD, Zilberman D, Jacobsen SE, Carrington JC: Genetic and functional diversification of small RNA pathways in plants. PLoS Biol. 2004, 2 (5): E104-10.1371/journal.pbio.0020104.View ArticlePubMedPubMed CentralGoogle Scholar
- Kasschau KD, Fahlgren N, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Carrington JC: Genome-wide profiling and analysis of Arabidopsis siRNAs. PLoS Biol. 2007, 5 (3): e57-10.1371/journal.pbio.0050057.View ArticlePubMedPubMed CentralGoogle Scholar
- Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One. 2007, 2 (2): e219-10.1371/journal.pone.0000219.View ArticlePubMedPubMed CentralGoogle Scholar
- Gasciolli V, Mallory AC, Bartel DP, Vaucheret H: Partially redundant functions of Arabidopsis DICER-like enzymes and a role for DCL4 in producing trans-acting siRNAs. Curr Biol. 2005, 15 (16): 1494-1500. 10.1016/j.cub.2005.07.024.View ArticlePubMedGoogle Scholar
- Nobuta K, Lu C, Shrivastava R, Pillay M, De Paoli E, Accerbi M, Arteaga-Vazquez M, Sidorenko L, Jeong DH, Yen Y, Green PJ, Chandler VL, Meyers BC: Distinct size distribution of endogeneous siRNAs in maize: Evidence from deep sequencing in the mop1-1 mutant. Proc Natl Acad Sci USA. 2008, 105 (39): 14958-14963. 10.1073/pnas.0808066105.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465. 10.1101/gr.075572.107.View ArticlePubMedPubMed CentralGoogle Scholar
- Heisel SE, Zhang Y, Allen E, Guo L, Reynolds TL, Yang X, Kovalic D, Roberts JK: Characterization of unique small RNA populations from rice grain. PLoS One. 2008, 3 (8): e2871-10.1371/journal.pone.0002871.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao T, Li G, Mi S, Li S, Hannon GJ, Wang XJ, Qi Y: A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev. 2007, 21 (10): 1190-1203. 10.1101/gad.1543507.View ArticlePubMedPubMed CentralGoogle Scholar
- Hertel J, Lindemeyer M, Missal K, Fried C, Tanzer A, Flamm C, Hofacker IL, Stadler PF: The expansion of the metazoan microRNA repertoire. BMC Genomics. 2006, 7: 25-10.1186/1471-2164-7-25.View ArticlePubMedPubMed CentralGoogle Scholar
- Niwa R, Slack FJ: The evolution of animal microRNA function. Curr Opin Genet Dev. 2007, 17 (2): 145-150. 10.1016/j.gde.2007.02.004.View ArticlePubMedGoogle Scholar
- Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu J, Shen Y, Wu Q, Kumar S, He B, Shi S, Carthew RW, Wang SM, Wu CI: The birth and death of microRNA genes in Drosophila. Nat Genet. 2008, 40 (3): 351-355. 10.1038/ng.73.View ArticlePubMedGoogle Scholar
- Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nature Cell Biology. 2009, 11 (3): 228-234. 10.1038/ncb0309-228.View ArticlePubMedGoogle Scholar
- Katoh T, Sakaguchi Y, Miyauchi K, Suzuki T, Kashiwabara S, Baba T: Selective stabilization of mammalian microRNAs by 3' adenylation mediated by the cytoplasmic poly(A) polymerase GLD-2. Genes Dev. 2009, 23 (4): 433-438. 10.1101/gad.1761509.View ArticlePubMedPubMed CentralGoogle Scholar
- Vaucheret H, Mallory AC, Bartel DP: AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell. 2006, 22 (1): 129-136. 10.1016/j.molcel.2006.03.011.View ArticlePubMedPubMed CentralGoogle Scholar
- Laubinger S, Sachsenberg T, Zeller G, Busch W, Lohmann JU, Ratsch G, Weigel D: Dual roles of the nuclear cap-binding complex and SERRATE in pre-mRNA splicing and microRNA processing in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2008, 105 (25): 8795-8800. 10.1073/pnas.0802493105.View ArticlePubMedPubMed CentralGoogle Scholar
- Gregory BD, O'Malley RC, Lister R, Urich MA, Tonti-Filippini J, Chen H, Millar AH, Ecker JR: A link between RNA metabolism and silencing affecting Arabidopsis development. Dev Cell. 2008, 14 (6): 854-866. 10.1016/j.devcel.2008.04.005.View ArticlePubMedGoogle Scholar
- Eamens AL, Smith NA, Curtin SJ, Wang MB, Waterhouse PM: The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes. RNA. 2009, 15 (12): 2219-2235. 10.1261/rna.1646909.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramachandran V, Chen X: Degradation of microRNAs by a family of exoribonucleases in Arabidopsis. Science. 2008, 321 (5895): 1490-1492. 10.1126/science.1163728.View ArticlePubMedPubMed CentralGoogle Scholar
- Chapman EJ, Carrington JC: Specialization and evolution of endogenous small RNA pathways. Nat Rev Genet. 2007, 8 (11): 884-896. 10.1038/nrg2179.View ArticlePubMedGoogle Scholar
- Ruiz-Ferrer V, Voinnet O: Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol. 2009, 60: 485-510. 10.1146/annurev.arplant.043008.092111.View ArticlePubMedGoogle Scholar
- Lu C, Kulkarni K, Souret FF, MuthuValliappan R, Tej SS, Poethig RS, Henderson IR, Jacobsen SE, Wang W, Green PJ, Meyers BC: MicroRNAs and other small RNAs enriched in the Arabidopsis RNA-dependent RNA polymerase-2 mutant. Genome Res. 2006, 16 (10): 1276-1288. 10.1101/gr.5530106.View ArticlePubMedPubMed CentralGoogle Scholar
- Jia Y, Lisch DR, Ohtsu K, Scanlon MJ, Nettleton D, Schnable PS: Loss of RNA-dependent RNA polymerase 2 (RDR2) function causes widespread and unexpected changes in the expression of transposons, genes, and 24-nt small RNAs. PLoS Genet. 2009, 5 (11): e1000737-10.1371/journal.pgen.1000737.View ArticlePubMedPubMed CentralGoogle Scholar
- Khraiwesh B, Arif MA, Seumel GI, Ossowski S, Weigel D, Reski R, Frank W: Transcriptional control of gene expression by microRNAs. Cell. 2010, 140 (1): 111-122. 10.1016/j.cell.2009.12.023.View ArticlePubMedGoogle Scholar
- Qi Y, He X, Wang XJ, Kohany O, Jurka J, Hannon GJ: Distinct catalytic and non-catalytic roles of ARGONAUTE4 in RNA-directed DNA methylation. Nature. 2006, 443 (7114): 1008-1012. 10.1038/nature05198.View ArticlePubMedGoogle Scholar
- Lacombe S, Nagasaki H, Santi C, Duval D, Piegu B, Bangratz M, Breitler JC, Guiderdoni E, Brugidou C, Hirsch J, Cao X, Brice C, Panaud O, Karlowski WM, Sato Y, Echeverria M: Identification of precursor transcripts for 6 novel miRNAs expands the diversity on the genomic organisation and expression of miRNA genes in rice. BMC Plant Biol. 2008, 8: 123-10.1186/1471-2229-8-123.View ArticlePubMedPubMed CentralGoogle Scholar
- Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci USA. 1988, 85 (8): 2444-2448. 10.1073/pnas.85.8.2444.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Schuster P: Fast Folding and Comparison of Rna Secondary Structures. Monatshefte Fur Chemie. 1994, 125 (2): 167-188. 10.1007/BF00818163.View ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Research. 2006, 34: D140-D144. 10.1093/nar/gkj112.View ArticlePubMedGoogle Scholar
- Perriere G, Gouy M: WWW-query: an on-line retrieval system for biological sequence banks. Biochimie. 1996, 78 (5): 364-369. 10.1016/0300-9084(96)84768-7.View ArticlePubMedGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302 (1): 205-217. 10.1006/jmbi.2000.4042.View ArticlePubMedGoogle Scholar
- Seibel PN, Muller T, Dandekar T, Schultz J, Wolf M: 4SALE--a tool for synchronous RNA sequence and secondary structure alignment and editing. BMC Bioinformatics. 2006, 7: 498-10.1186/1471-2105-7-498.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilm A, Linnenbrink K, Steger G: ConStruct: Improved construction of RNA consensus structures. BMC Bioinformatics. 2008, 9: 219-10.1186/1471-2105-9-219.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernhart SH, Hofacker IL, Will S, Gruber AR, Stadler PF: RNAalifold: improved consensus structure prediction for RNA alignments. BMC Bioinformatics. 2008, 9: 474-10.1186/1471-2105-9-474.View ArticlePubMedPubMed CentralGoogle Scholar
- Hofacker IL, Fekete M, Flamm C, Huynen MA, Rauscher S, Stolorz PE, Stadler PF: Automatic detection of conserved RNA structure elements in complete RNA virus genomes. Nucleic Acids Res. 1998, 26 (16): 3825-3836. 10.1093/nar/26.16.3825.View ArticlePubMedPubMed CentralGoogle Scholar
- Andersen ES, Lind-Thomsen A, Knudsen B, Kristensen SE, Havgaard JH, Torarinsson E, Larsen N, Zwieb C, Sestoft P, Kjems J, Gorodkin J: Semiautomated improvement of RNA alignments. RNA. 2007, 13 (11): 1850-1859. 10.1261/rna.215407.View ArticlePubMedPubMed CentralGoogle Scholar
- Gorodkin J, Heyer LJ, Brunak S, Stormo GD: Displaying the information contents of structural RNA alignments: the structure logos. Computer Applications in the Biosciences. 1997, 13 (6): 583-586.PubMedGoogle Scholar
- Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.