Diversity,expression and mRNA targeting abilities of Argonaute-targeting miRNAs among selected vascular plants

Background

Micro (mi)RNAs are important regulators of plant development. Across plant lineages, Dicer-like 1 (DCL1) proteins process long ds-like structures to produce micro (mi) RNA duplexes in a stepwise manner. These miRNAs are incorporated into Argonaute (AGO) proteins and influence expression of RNAs that have sequence complementarity with miRNAs. Expression levels of AGOs are greatly regulated by plants in order to minimize unwarranted perturbations using miRNAs to target mRNAs coding for AGOs. AGOs may also have high promoter specificity-sometimes expression of AGO can be limited to just a few cells in a plant. Viral pathogens utilize various means to counter antiviral roles of AGOs including hijacking the host encoded miRNAs to target AGOs. Two host encoded miRNAs namely miR168 and miR403 that target AGOs have been described in the model plant Arabidopsis and such a mechanism is thought to be well conserved across plants because AGO sequences are well conserved.

Results

We show that the interaction between AGO mRNAs and miRNAs is species-specific due to the diversity in sequences of two miRNAs that target AGOs, sequence diversity among corresponding target regions in AGO mRNAs and variable expression levels of these miRNAs among vascular plants. We used miRNA sequences from 68 plant species representing 31 plant families for this analysis. Sequences of miR168 and miR403 are not conserved among plant lineages, but surprisingly they differ drastically in their sequence diversity and expression levels even among closely related plants. Variation in miR168 expression among plants correlates well with secondary structures/length of loop sequences of their precursors.

Conclusions

Our data indicates a complex AGO targeting interaction among plant lineages due to miRNA sequence diversity and sequences of miRNA targeting regions among AGO mRNAs, thus leading to the assumption that the perturbations by viruses that use host miRNAs to target antiviral AGOs can only be species-specific. We also show that rapid evolution and likely loss of expression of miR168 isoforms in tobacco is related to the insertion of MITE-like transposons between miRNA and miRNA* sequences, a possible mechanism showing how miRNAs are lost in few plant lineages even though other close relatives have abundantly expressing miRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1049) contains supplementary material, which is available to authorized users.     

Sequence Fingerprints of MicroRNA Conservation

It is known that the conservation of protein-coding genes is associated with their sequences both various species, such as animals and plants. However, the association between microRNA (miRNA) conservation and their sequences in various species remains unexplored. Here we report the association of miRNA conservation with its sequence features, such as base content and cleavage sites, suggesting that miRNA sequences contain the fingerprints for miRNA conservation. More interestingly, different species show different and even opposite patterns between miRNA conservation and sequence features. For example, mammalian miRNAs show a positive/negative correlation between conservation and AU/GC content, whereas plant miRNAs show a negative/positive correlation between conservation and AU/GC content. Further analysis puts forward the hypothesis that the introns of protein-coding genes may be a main driving force for the origin and evolution of mammalian miRNAs. At the 5′ end, conserved miRNAs have a preference for base U, while less-conserved miRNAs have a preference for a non-U base in mammals. This difference does not exist in insects and plants, in which both conserved miRNAs and less-conserved miRNAs have a preference for base U at the 5′ end. We further revealed that the non-U preference at the 5′ end of less-conserved mammalian miRNAs is associated with miRNA function diversity, which may have evolved from the pressure of a highly sophisticated environmental stimulus the mammals encountered during evolution. These results indicated that miRNA sequences contain the fingerprints for conservation, and these fingerprints vary according to species. More importantly, the results suggest that although species share common mechanisms by which miRNAs originate and evolve, mammals may develop a novel mechanism for miRNA origin and evolution. In addition, the fingerprint found in this study can be predictor of miRNA conservation, and the findings are helpful in achieving a clearer understanding of miRNA function and evolution.     

Genome-Wide Identification and Analysis of MicroRNAs Involved in Witches’-Broom Phytoplasma Response in Ziziphus jujuba

MicroRNAs (miRNAs) play an important role in responding to biotic and abiotic stresses in plants. Jujube witches’-broom a phytoplasma disease of Ziziphus jujuba is prevalent in China and is a serious problem to the industry. However, the molecular mechanism of the disease is poorly understood. In this study, genome-wide identification and analysis of microRNAs in response to witches’-broom was performed. A total of 85 conserved miRNA unique sequences belonging to 32 miRNA families and 24 novel miRNA unique sequences, including their complementary miRNA* strands were identified from small RNA libraries derived from a uninfected and witches’-broom infected Z. jujuba plant. Differentially expressed miRNAs associated with Jujube witches’-broom disease were investigated between the two libraries, and 12 up-regulated miRNAs and 10 down- regulated miRNAs identified with more than 2 fold changes. Additionally, 40 target genes of 85 conserved miRNAs and 49 target genes of 24 novel miRNAs were predicted and their putative functions assigned. Using the modified 5’-RACE method, we confirmed that SPL and MYB were cleaved by miR156 and miR159, respectively. Our results provide insight into the molecular mechanisms of witches’-broom disease in Z. jujuba.     

Computational identification of microRNAs in peach expressed sequence tags and validation of their precise sequences by miR-RACE

Twenty-two potential miRNAs from seven miRNA families were first predicted from more than 80,857 EST sequences of peach (Prunus persica). Using two specific 5′ and 3′ miRNA RACE (miR-RACE) PCR reactions and sequence-directed cloning, we accurately determined the precise sequences, especially both ends, of eight candidate miRNAs. The sequencing results demonstrated that the ppe-miRNAs were conserved to those that were predicted computationally except ppe-miR171b. We validated the existence of two members (ppe-miR171a and miR171b) of the miR171 family in peach that belonged to different precursors. qRT-PCR was further employed in analyzing expression of the eight miRNAs in peach leaves, flowers, and fruits at different developing stages, where some of the miRNAs showed tissue-specific expression.     

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.