Comprehensive investigation of microRNAs enhanced by analysis of sequence variants, expression patterns, ARGONAUTE loading, and target cleavage

Somatic DNA demethylation generates tissue-specific methylation states and impacts flowering time

Cytosine methylation is a reversible epigenetic modification of DNA. In plants, removal of cytosine methylation is accomplished by the four members of the DEMETER (DME) family of 5-methylcytosine DNA glycosylases, named DME, DEMETER-LIKE2 (DML2), DML3, and REPRESSOR OF SILENCING1 (ROS1) in Arabidopsis thaliana. Demethylation by DME is critical for seed development, preventing experiments to determine the function of the entire gene family in somatic tissues by mutant analysis. Here, we bypassed the reproductive defects of dme mutants to create somatic quadruple homozygous mutants of the entire DME family. dme; ros1; dml2; and dml3 (drdd) leaves exhibit hypermethylated regions compared with wild-type leaves and rdd triple mutants, indicating functional redundancy among all four demethylases. Targets of demethylation include regions co-targeted by RNA-directed DNA methylation and, surprisingly, CG gene body methylation, indicating dynamic methylation at these less-understood sites. Additionally, many tissue-specific methylation differences are absent in drdd, suggesting a role for active demethylation in generating divergent epigenetic states across wild-type tissues. Furthermore, drdd plants display an early flowering phenotype, which involves 5'-hypermethylation and transcriptional down-regulation of FLOWERING LOCUS C. Active DNA demethylation is therefore required for proper methylation across somatic tissues and defines the epigenetic landscape of intergenic and coding regions.

MicroRNA Profiling During Mulberry (Morus atropurpurea Roxb) Fruit Development and Regulatory Pathway of miR477 for Anthocyanin Accumulation

To understand the mechanism of small non-coding RNAs (miRNA)-mediated development and ripening of mulberry fruits, three small RNA libraries from mulberry fruits at different development stages were constructed, and 159 conserved miRNAs as well as 86 novel miRNAs were successfully identified. Among the miRNAs identified, there were 90 miRNAs which showed differential expression patterns at different stages of fruit development and ripening. The target genes of these differential expressed (DE) miRNAs were involved in growth and development, transcription and regulation of transcription, metabolic processes, and etc. Interestingly, it was found that the expression level of mul-miR477 was increased with fruit ripening, and it can target the antisense lncRNA (Mul-ABCB19AS) of the ATP binding cassette (ABC) transporter B 19 gene (Mul-ABCB19). Our results showed that mul-miR477 can repress the expression of Mul-ABCB19AS and increase the expression of Mul-ABCB19, and it acted as a positive regulator participating anthocyanin accumulation through the regulatory network of mul-miR477-Mul-ABCB19AS-Mul-ABCB19.

Regulatory non-coding RNAs: a new frontier in regulation of plant biology

Beyond the most crucial roles of RNA molecules as a messenger, ribosomal, and transfer RNAs, the regulatory role of many non-coding RNAs (ncRNAs) in plant biology has been recognized. ncRNAs act as riboregulators by recognizing specific nucleic acid targets through homologous sequence interactions to regulate plant growth, development, and stress responses. Regulatory ncRNAs, ranging from small to long ncRNAs (lncRNAs), exert their control over a vast array of biological processes. Based on the mode of biogenesis and their function, ncRNAs evolved into different forms that include microRNAs (miRNAs), small interfering RNAs (siRNAs), miRNA variants (isomiRs), lncRNAs, circular RNAs (circRNAs), and derived ncRNAs. This article explains the different classes of ncRNAs and their role in plant development and stress responses. Furthermore, the applications of regulatory ncRNAs in crop improvement, targeting agriculturally important traits, have been discussed.

Barley Seeds miRNome Stability during Long-Term Storage and Aging

Seed aging is a complex biological process that has been attracting scientists’ attention for many years. High-throughput small RNA sequencing was applied to examine microRNAs contribution in barley seeds senescence. Unique samples of seeds that, despite having the same genetic makeup, differed in viability after over 45 years of storage in a dry state were investigated. In total, 61 known and 81 novel miRNA were identified in dry seeds. The highest level of expression was found in four conserved miRNA families, i.e., miR159, miR156, miR166, and miR168. However, the most astonishing result was the lack of significant differences in the level of almost all miRNAs in seed samples with significantly different viability. This result reveals that miRNAs in dry seeds are extremely stable. This is also the first identified RNA fraction that is not deteriorating along with the loss of seed viability. Moreover, the novel miRNA hvu-new41, with higher expression in seeds with the lowest viability as detected by RT-qPCR, has the potential to become an indicator of the decreasing viability of seeds during storage in a dry state.

Novel approaches on identification of conserved miRNAs for broad-spectrum Potyvirus control measures

Potyviridae comprises more than 200 ssRNA viruses, many of which have a broad host range and geographical distributions. Potyvirids (members of Potyviridae) infect several economically important plants such as saffron, cardamom, cucumber, pepper, potato, tomato, yam, etc. Cumulatively, potyvirids cause a substantial economic loss. The major bottleneck in developing an efficient antiviral strategy is that viruses quickly evade host immunity owing to their higher mutation and recombination rates. Due to this reason, the emergence of newer and improved broad-spectrum approaches to combat viral infections is essential. The use of microRNA's (miRNA) to circumvent viral infection against animal viruses has been successfully employed. Fewer studies reported the development of efficient miRNA-based antivirus resistant strategies against plant viruses and none focused on multiple virus resistance. We focused on potyviruses since studies are limited and identification of conserved miRNAs among various host plants would be an initiative to design broad-spectrum antivirus strategies. In this study, we predicted evolutionarily conserved miRNAs by BLAST searching of reported miRNAs from 15 plants against the GSS and EST sequences of banana. A total of nine miRNAs were predicted and screened in nine diverse potyvirids' hosts (Banana, Tomato, Green gram, Jasmine, Chilli, Coriander, Onion, Rose and Colocasia) belonging to eight different orders (Zingiberales, Solanales, Fabales, Lamiales, Apiales, Asperagales, Rosales and Alismatales). Results suggested that miR168 and miR162 are conserved among all the selected plants. This comprehensive study laid the foundations to design broad-spectrum antivirus resistance using miRNAs. To conclude miR168 and miR162 are conserved among many plants and play a crucial role in evading virus infection which could be used as a potential candidate for developing antiviral strategies against potyvirid infections.

Genome-wide analysis of brassinosteroid responsive small RNAs in Arabidopsis thaliana

BACKGROUND

Brassinosteroids (BRs) are a class of phytohormones with important roles in regulating physiological and developmental processes. Small RNAs, including small interfering RNAs and microRNAs (miRNAs), are non-protein coding RNAs that regulate gene expression at the transcriptional and post-transcriptional levels. However, the roles of small RNAs in BR response have not been studied well.

OBJECTIVE

In this study, we aimed to identify BR-responsive small RNA clusters and miRNAs in Arabidopsis. In addition, the effect of BR-responsive small RNAs on their transcripts and target genes were examined.

METHODS

Small RNA libraries were constructed from control and epibrassinolide-treated seedlings expressing wild-type BRI1-Flag protein under its native promoter in the bri1-5 mutant. After sequencing the small RNA libraries, differentially expressed small RNA clusters were identified by examining the expression levels of small RNAs in 100-nt bins of the Arabidopsis genome. To identify the BR-responsive miRNAs, the expression levels of all the annotated mature miRNAs, registered in miRBase, were analyzed. Previously published RNA-seq data were utilized to monitor the BR-responsive expression patterns of differentially expressed small RNA clusters and miRNA target genes.

RESULTS

In results, 38 BR-responsive small RNA clusters, including 30 down-regulated and eight up-regulated clusters, were identified. These differentially expressed small RNA clusters were from miRNA loci, transposons, protein-coding genes, pseudogenes and others. Of these, a transgene, BRI1, accumulates small RNAs, which are not found in the wild type. Small RNAs in this transgene are up-regulated by BRs while BRI1 mRNA is down-regulated by BRs. By analyzing the expression patterns of mature miRNAs, we have identified BR-repressed miR398a-5p and BR-induced miR156g. Although miR398a-5p is down-regulated by BRs, its predicted targets were not responsive to BRs. However, SPL3, a target of BR-inducible miR156g, is down-regulated by BRs.

CONCLUSION

BR-responsive small RNAs and miRNAs identified in this study will provide an insight into the role of small RNAs in BR responses in plants. Especially, we suggest that miR156g/SPL3 module might play a role in BR-mediated growth and development in Arabidopsis.

Diff isomiRs: Large-scale detection of differential isomiRs for understanding non-coding regulated stress omics in plants

Plants have an amazing ability to cope with wide variety of stresses by regulating the expression of genes and thus by altering the physiological status. In the past few years, canonical microRNA variants (isomiRs) have been shown to play pivotal roles by acting as regulators of the transcriptional machinery. In the present research, we present Diff isomiRs, a web-based exploratory repository of differential isomiRs across 16 sequenced plant species representing a total of 433 datasets across 21 different stresses and 158 experimental states. Diff isomiRs provides the high-throughput detection of differential isomiRs using mapping-based and model-based differential analysis revealing a total of 16,157 and 2,028 differential isomiRs, respectively. Easy-to-use and web-based exploration of differential isomiRs provides several features such as browsing of the differential isomiRs according to stress or species, as well as association of the differential isomiRs to targets and plant endogenous target mimics (PeTMs). Diff isomiRs also provides the relationship between the canonical miRNAs, isomiRs and the miRNA-target interactions. This is the first web-based large-scale repository for browsing differential isomiRs and will facilitate better understanding of the regulatory role of the isomiRs with respect to the canonical microRNAs. Diff isomiRs can be accessed at: www.mcr.org.in/diffisomirs.

Genome-Wide Identification and Transcriptional Expression of the PAL Gene Family in Common Walnut (Juglans Regia L.)

Juglans regia L. is an economically important crop cultivated worldwide for its high quality and quantity of wood and nuts. Phenylalanine ammonia-lyase (PAL) is the first enzyme in the phenylpropanoid pathway that plays a critical role in plant growth, development, and adaptation, but there have been few reports of the PAL gene family in common walnut. Here, we report a genome-wide study of J. regiaPAL genes and analyze their phylogeny, duplication, microRNA, and transcriptional expression. A total of 12 PAL genes were identified in the common walnut and clustered into two subfamilies based on phylogenetic analysis. These common walnut PALs are distributed on eight different pseudo-chromosomes. Seven of the 12 PALs (JrPAL2-3, JrPAL4-2, JrPAL2-1, JrPAL4-1, JrPAL8, JrPAL9, and JrPAL6) were specific found in J. regia, and JrPAL3, JrPAL5, JrPAL1-2, JrPAL7, and JrPAL2-2 were found to be closely associated with the woody plant Populus trichocarpa. Additionally, the expression patterns of JrPAL3, JrPAL7, JrPAL9, and JrPAL2-1 showed that they had high expression in female and male flowers. The miRNA ath-miR830-5p regulates two genes, JrPAL5 and JrPAL1, such that they have low expression in the male and female flowers of the common walnut. Our research provides useful information for further research into the function of PAL genes in common walnut and Juglans.

MicroR408 regulates defense response upon wounding in sweet potato

MiRNAs play diverse roles in plant development and defense responses by binding to their mRNA targets based on sequence complementarity. Here, we investigated a wound-related miR408 and its target genes in sweet potato (Ipomoea batatas) by small RNA deep sequencing and transcriptome analysis. The expression patterns of miR408 and the miR408 precursor were significantly repressed by wounding and jasmonate (JA). In contrast, expression of the putative target genes IbKCS (3-ketoacyl-CoA synthase 4), IbPCL (plantacyanin), and IbGAUT (galacturonosyltransferase 7-like) of miR408 was increased following wounding, whereas only IbKCS was increased after JA treatment. Target cleavage site mapping and Agrobacterium-mediated transient assay demonstrated that IbKCS, IbPCL, and IbGAUT were the targets of miR408. The expression of miR408 target genes was repressed in transgenic sweet potatoes overexpressing miR408. These data indicated a relationship between miR408 and its target genes. Notably, miR408-overexpressing plants showed a semi-dwarf phenotype and attenuated resistance to insect feeding, while transgenic plants overexpressing IbKCS exhibited more insect resistance than plants overexpressing only the empty vector. Collectively, sweet potato reduces the abundance of miR408 upon wounding to elevate the expression of IbKCS, IbPCL, and IbGAUT. The expression of IbKCS enhances the defense system against herbivore wounding.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

MiRNA-seq-based profiles of miRNAs in mulberry phloem sap provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

A wide range of miRNAs have been identified as phloem-mobile molecules that play important roles in coordinating plant development and physiology. Phytoplasmas are associated with hundreds of plant diseases, and the pathogenesis involved in the interactions between phytoplasmas and plants is still poorly understood. To analyse the molecular mechanisms of phytoplasma pathogenicity, the miRNAs profiles in mulberry phloem saps were examined in response to phytoplasma infection. A total of 86 conserved miRNAs and 19 novel miRNAs were identified, and 30 conserved miRNAs and 13 novel miRNAs were differentially expressed upon infection with phytoplasmas. The target genes of the differentially expressed miRNAs are involved in diverse signalling pathways showing the complex interactions between mulberry and phytoplasma. Interestingly, we found that mul-miR482a-5p was up-regulated in the infected phloem saps, and grafting experiments showed that it can be transported from scions to rootstock. Based on the results, the complexity and roles of the miRNAs in phloem sap and the potential molecular mechanisms of their changes were discussed. It is likely that the phytoplasma-responsive miRNAs in the phloem sap modulate multiple pathways and work cooperatively in response to phytoplasma infection, and their expression changes may be responsible for some symptoms in the infected plants.

Bioinformatics resources for deciphering the biogenesis and action pathways of plant small RNAs

The next-generation sequencing (NGS) technology has revolutionized our previous understanding of the plant genomes, relying on its innate advantages, such as high throughput and deep sequencing depth. In addition to the protein-coding gene loci, massive transcription signals have been detected within intergenic or intragenic regions. Most of these signals belong to non-coding ones, considering their weak protein-coding potential. Generally, these transcripts could be divided into long non-coding RNAs and small non-coding RNAs (sRNAs) based on their sequence length. In addition to the well-known microRNAs (miRNAs), many plant endogenous sRNAs were collectively referred to as small interfering RNAs. However, an increasing number of unclassified sRNA species are being discovered by NGS. The high heterogeneity of these novel sRNAs greatly hampered the mechanistic studies, especially on the clear description of their biogenesis and action pathways. Fortunately, public databases, bioinformatics softwares and NGS datasets are increasingly available for plant sRNA research. Here, by summarizing these valuable resources, we proposed a general workflow to decipher the RDR (RNA-dependent RNA polymerase)- and DCL (Dicer-like)-dependent biogenesis pathways, and the Argonaute-mediated action modes (such as target cleavages and chromatin modifications) for specific sRNA species in plants. Taken together, we hope that by summarizing a list of the public resources, this work will facilitate the plant biologists to perform classification and functional characterization of the interesting sRNA species.

Identification of miRNA-mediated drought responsive multi-tiered regulatory network in drought tolerant rice, Nagina 22

Comparative characterization of microRNA-mediated stress regulatory networks in contrasting rice cultivars is critical to decipher plant stress response. Consequently, a multi-level comparative analysis, using sRNA sequencing, degradome analysis, enzymatic and metabolite assays and metal ion analysis, in drought tolerant and sensitive rice cultivars was conducted. The study identified a group of miRNAs "Cultivar-specific drought responsive" (CSDR)-miRNAs (osa-miR159f, osa-miR1871, osa-miR398b, osa-miR408-3p, osa-miR2878-5p, osa-miR528-5p and osa-miR397a) that were up-regulated in the flag-leaves of tolerant cultivar, Nagina 22 (N22) and Vandana, but down-regulated in the sensitive cultivar, Pusa Basmati 1 (PB1) and IR64, during drought. Interestingly, CSDR-miRNAs target several copper-protein coding transcripts like plantacyanins, laccases and Copper/Zinc superoxide dismutases (Cu/Zn SODs) and are themselves found to be similarly induced under simulated copper-starvation in both N22 and PB1. Transcription factor OsSPL9, implicated in Cu-homeostasis also interacted with osa-miR408-3p and osa-miR528-5p promoters. Further, N22 flag leaves showed lower SOD activity, accumulated ROS and had a higher stomata closure. Interestingly, compared to PB1, internal Cu levels significantly decreased in the N22 flag-leaves, during drought. Thus, the study identifies the unique drought mediated dynamism and interplay of Cu and ROS homeostasis, in the flag leaves of drought tolerant rice, wherein CSDR-miRNAs play a pivotal role.

Loss of function of Arabidopsis microRNA-machinery genes impairs fertility, and has effects on homologous recombination and meiotic chromatin dynamics

MicroRNAs (miRNAs) are ~22-nt single-stranded noncoding RNAs with regulatory roles in a wide range of cellular functions by repressing eukaryotic gene expression at a post-transcriptional level. Here, we analyzed the effects on meiosis and fertility of hypomorphic or null alleles of the HYL1, HEN1, DCL1, HST and AGO1 genes, which encode miRNA-machinery components in Arabidopsis. Reduced pollen and megaspore mother cell number and fertility were shown by the mutants analyzed. These mutants also exhibited a relaxed chromatin conformation in male meiocytes at the first meiotic division, and increased chiasma frequency, which is likely to be due to increased levels of mRNAs from key genes involved in homologous recombination. The hen1-13 mutant was found to be hypersensitive to gamma irradiation, which mainly causes double-strand breaks susceptible to be repaired by homologous recombination. Our findings uncover a role for miRNA-machinery components in Arabidopsis meiosis, as well as in the repression of key genes required for homologous recombination. These genes seem to be indirect miRNA targets.

Mutations in Argonaute5 Illuminate Epistatic Interactions of the K1 and I Loci Leading to Saddle Seed Color Patterns in Glycine max

The soybean (Glycine max) seed coat has distinctive, genetically programmed patterns of pigmentation, and the recessive k1 mutation can epistatically overcome the dominant I and ii alleles, which inhibit seed color by producing small interfering RNAs (siRNAs) targeting chalcone synthase (CHS) mRNAs. Small RNA sequencing of dissected regions of immature seed coats demonstrated that CHS siRNA levels cause the patterns produced by the ii and ik alleles of the I locus, which restrict pigment to the hilum or saddle region of the seed coat, respectively. To identify the K1 locus, we compared RNA-seq data from dissected regions of two Clark isolines having similar saddle phenotypes mediated by CHS siRNAs but different genotypes (homozygous ik K1 versus homozygous ii k1). By examining differentially expressed genes, mapping information, and genome resequencing, we identified a 129-bp deletion in Glyma.11G190900 encoding Argonaute5 (AGO5), a member of the Argonaute family. Amplicon sequencing of several independent saddle pattern mutants from different genetic backgrounds revealed independent lesions affecting AGO5, thus establishing Glyma.11G190900 as the K1 locus. Nonfunctional AGO5 from k1 alleles leads to altered distributions of CHS siRNAs, thus explaining how the k1 mutation reverses the phenotype of the seed coat regions from yellow to pigmented, even in the presence of the normally dominant I or ii alleles.

Characterization and Function of MicroRNA∗s in Plants

MicroRNAs, a group of non-coding RNA molecules, play essential roles in a wide range of cellular processes in different molecules, cells, and organisms. In plants, microRNAs are a class of 20- to 24-nucleotides endogenous small RNAs that repress gene expression. The microRNA guide strand (miRNA) and its complementary strand (miRNA^∗) both originate from the miRNA/miRNA^∗ duplex. Generally, the guide strands act as post-transcriptional regulators that suppress gene expression by cleaving their target mRNA transcripts, whereas the complementary strands were thought to be degraded as 'passenger strands.' However, the complementary strand has been confirmed to possess significant biological functionality in recent reports. In this review, we summarized the binding characteristics of the miRNA^∗ strands with ARGONAUTE proteins, their tissue-specific accumulations and their biological functions, illustrating the essential roles of miRNA^∗s in biological processes and therefore providing directions for further exploration.

MicroRNA858 Is a Potential Regulator of Phenylpropanoid Pathway and Plant Development

MicroRNAs (miRNAs) are endogenous, noncoding small RNAs that function as critical regulators of gene expression. In plants, miRNAs have shown their potential as regulators of growth, development, signal transduction, and stress tolerance. Although the miRNA-mediated regulation of several processes is known, the involvement of miRNAs in regulating secondary plant product biosynthesis is poorly understood. In this study, we functionally characterized Arabidopsis (Arabidopsis thaliana) miR858a, which putatively targets R2R3-MYB transcription factors involved in flavonoid biosynthesis. Overexpression of miR858a in Arabidopsis led to the down-regulation of several MYB transcription factors regulating flavonoid biosynthesis. In contrast to the robust growth and early flowering of miR858OX plants, reduction of plant growth and delayed flowering were observed in Arabidopsis transgenic lines expressing an artificial miRNA target mimic (MIM858). Genome-wide expression analysis using transgenic lines suggested that miR858a targets a number of regulatory factors that modulate the expression of downstream genes involved in plant development and hormonal and stress responses. Furthermore, higher expression of MYBs in MIM858 lines leads to redirection of the metabolic flux towards the synthesis of flavonoids at the cost of lignin synthesis. Altogether, our study has established the potential role of light-regulated miR858a in flavonoid biosynthesis and plant growth and development.

Functional and evolutionary analyses of the miR156 and miR529 families in land plants

BACKGROUND

MicroRNAs (miRNAs) are important regulatory elements of gene expression. Similarly to coding genes, miRNA genes follow a birth and death pattern of evolution likely reflecting functional relevance and divergence. For instance, miRNA529 is evolutionarily related to miRNA156 (a highly conserved miRNA in land plants), but it is lost in Arabidopsis thaliana. Interestingly, both miRNAs target sequences overlap in some members of the SQUAMOSA promoter-binding protein like (SPL) family, raising important questions regarding the diversification of the miR156/miR529-associated regulatory network in land plants.

RESULTS

In this study, through phylogenic reconstruction of miR156/529 target sequences from several taxonomic groups, we have found that specific eudicot SPLs, despite miRNA529 loss, retained the corresponding target site. Detailed molecular evolutionary analyses of miR156/miR529-target sequence showed that loss of miR529 in core eudicots, such as Arabidopsis, is correlated with a more relaxed selection of the miRNA529 specific target element, while miRNA156-specific target sequence is under stronger selection, indicating that these two target sites might be under distinct evolutionary constraints. Importantly, over-expression in Arabidopsis of MIR529 precursor from a monocot, but not from a basal eudicot, demonstrates specific miR529 regulation of AtSPL9 and AtSPL15 genes, which contain conserved responsive elements for both miR156 and miR529.

CONCLUSIONS

Our results suggest loss of functionality of MIR529 genes in the evolutionary history of eudicots and show that the miR529-responsive element present in some eudicot SPLs is still functional. Our data support the notion that particular miRNA156 family members might have compensated for the loss of miR529 regulation in eudicot species, which concomitantly may have favored diversification of eudicot SPLs.

Efficient and Inducible Use of Artificial MicroRNAs in Marchantia polymorpha

We describe the efficient use of artificial microRNAs (amiRs) in Marchantia polymorpha using both endogenous and heterologous primary microRNA (pri-miR) hairpin backbones. Targeting of two transcription factor genes, MpARF1 and MpRR-B, mediating different hormonal responses, demonstrated that amiRs can create specific and reproducible physiological and morphological defects, facilitating interpretation of gene function. A third amiR was designed to target a gene encoding a component of the Polycomb recessive complex 2, MpE(z), and constitutive expression of this amiR results in sporeling lethality. Adaptation of an estrogen-inducible system allowed analysis of the phenotypic effects of induction of this amiR during other stages of the life cycle. We discuss the advantages and challenges of the use of amiRs as a tool for reverse genetic analysis in M. polymorpha.

MicroRNA miR396 Regulates the Switch between Stem Cells and Transit-Amplifying Cells in Arabidopsis Roots

To ensure an adequate organ mass, the daughters of stem cells progress through a transit-amplifying phase displaying rapid cell division cycles before differentiating. Here, we show that Arabidopsis thaliana microRNA miR396 regulates the transition of root stem cells into transit-amplifying cells by interacting with GROWTH-REGULATING FACTORs (GRFs). The GRFs are expressed in transit-amplifying cells but are excluded from the stem cells through inhibition by miR396. Inactivation of the GRFs increases the meristem size and induces periclinal formative divisions in transit-amplifying cells. The GRFs repress PLETHORA (PLT) genes, regulating their spatial expression gradient. Conversely, PLT activates MIR396 in the stem cells to repress the GRFs. We identified a pathway regulated by GRF transcription factors that represses stem cell-promoting genes in actively proliferating cells, which is essential for the progression of the cell cycle and the orientation of the cell division plane. If unchecked, the expression of the GRFs in the stem cell niche suppresses formative cell divisions and distorts the organization of the quiescent center. We propose that the interactions identified here between miR396 and GRF and PLT transcription factors are necessary to establish the boundary between the stem cell niche and the transit-amplifying region.

MicroRNA nomenclature and the need for a revised naming prescription

A central environment and interface for microRNA (miRNA) registry and repository and a general standardized framework for their systematic annotation was established over a decade ago. However, the numbers of experimentally and computationally identified miRNAs are swiftly accumulating, and new aspects of miRNA-mediated gene regulation are being revealed. Currently, it is of great significance that the annotation framework should be redefined to include newly discovered miRNA species such as the variants of mature miRNAs (isomiRNAs), and organellar miRNAs: cipomiRNAs and mitomiRNAs. It is also of great importance that key terminology referring to the novelty, evolutionary history or biogenesis of miRNAs, as well as the confidence of their identification are standardized in the literature and disseminated in a central miRNA registry. Here, we review the status of miRNA nomenclature, curation and critical points of need for a revision of miRNA nomenclature and terminology.

The Complexity of Posttranscriptional Small RNA Regulatory Networks Revealed by In Silico Analysis of Gossypium arboreum L. Leaf, Flower and Boll Small Regulatory RNAs

MicroRNAs (miRNAs) and secondary small interfering RNAs (principally phased siRNAs or trans-acting siRNAs) are two distinct subfamilies of small RNAs (sRNAs) that are emerging as key regulators of posttranscriptional gene expression in plants. Both miRNAs and secondary-siRNAs (sec-siRNAs) are processed from longer RNA precursors by DICER-LIKE proteins (DCLs). Gossypium arboreum L., also known as tree cotton or Asian cotton, is a diploid, possibly ancestral relative of tetraploid Gossypium hirsutum L., the predominant type of commercially grown cotton worldwide known as upland cotton. To understand the biological significance of these gene regulators in G. arboreum, a bioinformatics analysis was performed on G. arboreum small RNAs produced from G. arboreum leaf, flower, and boll tissues. Consequently, 263 miRNAs derived from 353 precursors, including 155 conserved miRNAs (cs-miRNAs) and 108 novel lineage-specific miRNAs (ls-miRNAs). Along with miRNAs, 2,033 miRNA variants (isomiRNAs) were identified as well. Those isomiRNAs with variation at the 3'-miRNA end were expressed at the highest levels, compared to other types of variants. In addition, 755 pha-siRNAs derived 319 pha-siRNA gene transcripts (PGTs) were identified, and the potential pha-siRNA initiators were predicted. Also, 2,251 non-phased siRNAs were found as well, of which 1,088 appeared to be produced by so-called cis- or trans-cleavage of the PGTs observed at positions differing from pha-siRNAs. Of those sRNAs, 148 miRNAs/isomiRNAs and 274 phased/non-phased siRNAs were differentially expressed in one or more pairs of tissues examined. Target analysis revealed that target genes for both miRNAs and pha-siRNAs are involved a broad range of metabolic and enzymatic activities. We demonstrate that secondary siRNA production could result from initial cleavage of precursors by both miRNAs or isomiRNAs, and that subsequently produced phased and unphased siRNAs could result that also serve as triggers of a second round of both cis- and trans-cleavage of additional siRNAs, leading to the formation of complex sRNA regulatory networks mediating posttranscriptional gene silencing. Results from this study extended our knowledge on G. arboreum sRNAs and their biological importance, which would facilitate future studies on regulatory mechanism of tissue development in cotton and other plant species.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

miRVine: a microRNA expression atlas of grapevine based on small RNA sequencing

Background

miRNAs are the most abundant class of small non-coding RNAs, and they are involved in post-transcriptional regulations, playing a crucial role in the refinement of genetic programming during plant development. Here we present a comprehensive picture of miRNA regulation in Vitis vinifera L. plant during its complete life cycle. Furthering our knowledge about the post-transcriptional regulation of plant development is fundamental to understand the biology of such an important crop.

Results

We analyzed 70 small RNA libraries, prepared from berries, inflorescences, tendrils, buds, carpels, stamens and other samples at different developmental stages. One-hundred and ten known and 175 novel miRNAs have been identified and a wide grapevine expression atlas has been described. The distribution of miRNA abundance reveals that 22 novel miRNAs are specific to stamen, and two of them are, interestingly, involved in ethylene biosynthesis, while only few miRNAs are highly specific to other organs. Thirty-eight miRNAs are present in all our samples, suggesting a role in key regulatory circuit. On the basis of miRNAs abundance and distribution across samples and on the estimated correlation, we suggest that miRNA expression define organ identity. We performed target prediction analysis and focused on miRNA expression analysis in berries and inflorescence during their development, providing an initial functional description of the identified miRNAs.

Conclusions

Our findings represent a very extensive miRNA expression atlas in grapevine, allowing the definition of how the spatio-temporal distribution of miRNAs defines organ identity. We describe miRNAs abundance in specific tissues not previously described in grapevine and contribute to future targeted functional analyses. Finally, we present a deep characterization of miRNA involvement in berry and inflorescence development, suggesting a role for miRNA-driven hormonal regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1610-5) contains supplementary material, which is available to authorized users.

Manipulating microRNAs for improved biomass and biofuels from plant feedstocks

Petroleum-based fuels are nonrenewable and unsustainable. Renewable sources of energy, such as lignocellulosic biofuels and plant metabolite-based drop-in fuels, can offset fossil fuel use and reverse environmental degradation through carbon sequestration. Despite these benefits, the lignocellulosic biofuels industry still faces many challenges, including the availability of economically viable crop plants. Cell wall recalcitrance is a major economic barrier for lignocellulosic biofuels production from biomass crops. Sustainability and biomass yield are two additional, yet interrelated, foci for biomass crop improvement. Many scientists are searching for solutions to these problems within biomass crop genomes. MicroRNAs (miRNAs) are involved in almost all biological and metabolic process in plants including plant development, cell wall biosynthesis and plant stress responses. Because of the broad functions of their targets (e.g. auxin response factors), the alteration of plant miRNA expression often results in pleiotropic effects. A specific miRNA usually regulates a biologically relevant bioenergy trait. For example, relatively low miR156 overexpression leads to a transgenic feedstock with enhanced biomass and decreased recalcitrance. miRNAs have been overexpressed in dedicated bioenergy feedstocks such as poplar and switchgrass yielding promising results for lignin reduction, increased plant biomass, the timing of flowering and response to harsh environments. In this review, we present the status of miRNA-related research in several major biofuel crops and relevant model plants. We critically assess published research and suggest next steps for miRNA manipulation in feedstocks for increased biomass and sustainability for biofuels and bioproducts.

Rapid divergence and high diversity of miRNAs and miRNA targets in the Camelineae

MicroRNAs (miRNAs) are short RNAs involved in gene regulation through translational inhibition and transcript cleavage. After processing from imperfect fold-back structures, miRNAs are incorporated into RNA-induced silencing complexes (RISCs) before targeting transcripts with varying degrees of complementarity. Some miRNAs are evolutionarily deep-rooted, and sequence complementarity with their targets is maintained through purifying selection. Both Arabidopsis and Capsella belong to the tribe Camelineae in the Brassicaceae, with Capsella rubella serving as an outgroup to the genus Arabidopsis. The genome sequence of C. rubella has recently been released, which allows characterization of its miRNA complement in comparison with Arabidopsis thaliana and Arabidopsis lyrata. Through next-generation sequencing, we identify high-confidence miRNA candidates specific to the C. rubella lineage. Only a few lineage-specific miRNAs have been studied for evolutionary constraints, and there have been no systematic studies of miRNA target diversity within or divergence between closely related plant species. Therefore we contrast sequence variation in miRNAs and their targets within A. thaliana, and between A. thaliana, A. lyrata and C. rubella. We document a surprising amount of small-scale variation in miRNA-target pairs, where many miRNAs are predicted to have species-specific targets in addition to ones that are shared between species. Our results emphasize that the transitive nature of many miRNA-target pairs can be observed even on a relatively short evolutionary time-scale, with non-random occurrences of differences in miRNAs and their complements in the miRNA precursors, the miRNA* sequences.

Differential expression of miRNAs and their target genes in senescing leaves and siliques: insights from deep sequencing of small RNAs and cleaved target RNAs

MicroRNAs (miRNAs) are a class of small RNAs, which typically function by guiding cleavage of target mRNAs. They are known to play roles in a variety of plant processes including development, responses to environmental stresses and senescence. To identify senescence regulation of miRNAs in Arabidopsis thaliana, eight small RNA libraries were constructed and sequenced at four different stages of development and senescence from both leaves and siliques, resulting in more than 200 million genome-matched sequences. Parallel analysis of RNA ends libraries, which enable the large-scale examination of miRNA-guided cleavage products, were constructed and sequenced, resulting in over 750 million genome-matched sequences. These large datasets led to the identification a new senescence-inducible small RNA locus, as well as new regulation of known miRNAs and their target genes during senescence, many of which have established roles in nutrient responsiveness and cell structural integrity. In keeping with remobilization of nutrients thought to occur during senescence, many miRNAs and targets had opposite expression pattern changes between leaf and silique tissues during the progression of senescence. Taken together, these findings highlight the integral role that miRNAs may play in the remobilization of resources and alteration of cellular structure that is known to occur in senescence.

Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses

Arabidopsis microRNA expression regulation was studied in a wide array of abiotic stresses such as drought, heat, salinity, copper excess/deficiency, cadmium excess, and sulfur deficiency. A home-built RT-qPCR mirEX platform for the amplification of 289 Arabidopsis microRNA transcripts was used to study their response to abiotic stresses. Small RNA sequencing, Northern hybridization, and TaqMan® microRNA assays were performed to study the abundance of mature microRNAs. A broad response on the level of primary miRNAs (pri-miRNAs) was observed. However, stress response at the level of mature microRNAs was rather confined. The data presented show that in most instances, the level of a particular mature miRNA could not be predicted based on the level of its pri-miRNA. This points to an essential role of posttranscriptional regulation of microRNA expression. New Arabidopsis microRNAs responsive to abiotic stresses were discovered. Four microRNAs: miR319a/b, miR319b.2, and miR400 have been found to be responsive to several abiotic stresses and thus can be regarded as general stress-responsive microRNA species.

Endogenous target mimics down-regulate miR160 mediation of ARF10, -16, and -17 cleavage during somatic embryogenesis in Dimocarpus longan Lour

MicroRNA160 plays a critical role in plant development by negatively regulating the auxin response factors ARF10, -16, and -17. However, the ways in which miR160 expression is regulated at the transcriptional level, and how miR160 interacts with its targets during plant embryo development, remain unknown. Here, we studied the regulatory relationships among endogenous target mimics (eTMs), and miR160 and its targets, and their involvement in hormone signaling and somatic embryogenesis (SE) in Dimocarpus longan. We identified miR160 family members and isolated the miR160 precursor, primary transcript, and promoter. The promoter contained cis-acting elements responsive to stimuli such as light, abscisic acid, salicylic acid (SA) and heat stress. The pri-miR160 was down-regulated in response to SA but up-regulated by gibberellic acid, ethylene, and methyl jasmonate treatment, suggesting that pri-miR160 was associated with hormone transduction. Dlo-miR160a, -a(∗) and -d(∗) reached expression peaks in torpedo-shaped embryos, globular embryos and cotyledonary embryos, respectively, but were barely detectable in friable-embryogenic callus. This suggests that they have expression-related and functional diversity, especially during the middle and later developmental stages of SE. Four potential eTMs for miR160 were identified. Two of them, glucan endo-1,3-beta- glucosidase-like protein 2-like and calpain-type cysteine protease DEK1, were confirmed to control the corresponding dlo-miR160a(∗) expression level. This suggests that they may function to abolish the binding between dlo-miR160a(∗) and its targets. These two eTMs also participated in 2,4-D and ABA signal transduction. DlARF10, -16, and -17 targeting by dlo-miR160a was confirmed; their expression levels were higher in friable-embryogenic callus and incomplete compact pro-embryogenic cultures and responded to 2,4-D, suggesting they may play a major role in the early stages of longan SE dependent on 2,4-D. The eTMs, miR160, and ARF10, -16, and -17 exhibited tissue specificity in 'Sijimi' longan vegetative and reproductive organs, but were not significant negatively correlated. These results provide insights into the possible role of the eTM-miR160-ARF10-16-17 pathway in longan somatic embryo development.

The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling

Barley (Hordeum vulgare L.) Mla alleles encode coiled-coil (CC), nucleotide binding, leucine-rich repeat (NB-LRR) receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). How Mla or NB-LRR genes in grass species are regulated at post-transcriptional level is not clear. The microRNA family, miR9863, comprises four members that differentially regulate distinct Mla alleles in barley. We show that miR9863 members guide the cleavage of Mla1 transcripts in barley, and block or reduce the accumulation of MLA1 protein in the heterologous Nicotiana benthamiana expression system. Regulation specificity is determined by variation in a unique single-nucleotide-polymorphism (SNP) in mature miR9863 family members and two SNPs in the Mla miR9863-binding site that separates these alleles into three groups. Further, we demonstrate that 22-nt miR9863s trigger the biogenesis of 21-nt phased siRNAs (phasiRNAs) and together these sRNAs form a feed-forward regulation network for repressing the expression of group I Mla alleles. Overexpression of miR9863 members specifically attenuates MLA1, but not MLA10-triggered disease resistance and cell-death signaling. We propose a key role of the miR9863 family in dampening immune response signaling triggered by a group of MLA immune receptors in barley.

Plants rely on cell-surface and intracellular immune receptors to sense pathogen invasion and to mediate defense responses. However, uncontrolled activation of immune responses is harmful to plant growth and development. Small RNAs have recently been shown to fine-tune the expression of intracellular immune receptors and contribute to the regulation of defense signaling in dicot plants, while similar processes have not been well documented in monocot grain crops, such as barley and wheat. Here, we show that, in barley, some members of the miR9863 family target a subset of Mla alleles that confer race-specific disease resistance to the powdery mildew fungus. These miRNAs act on Mla transcripts by cleavage and translational repression. Production of a type of trans-acting small RNAs, designated as phasiRNAs, enhances the effects of miRNA regulation on Mla targets. We propose that Mla-mediated immune signaling is fine-tuned by the miRNAs at later stage of MLA activation to avoid overloading of immune responses in barley cells.

Is Argonaute 1 the only effective slicer of small RNA-mediated regulation of gene expression in plants?

To maintain normal growth and development in a plant, gene expression must be under strict surveillance. One of the post-transcriptional regulatory pathways involves small RNA (sRNA)-guided, Argonaute (AGO) protein complex-mediated target cleavages. MicroRNAs (miRNAs) are the well-known sRNA species participating in the cleavage-based regulation of gene expression. In plants, most miRNAs are incorporated into AGO1-associated silencing complexes. Thus, the AGO1 protein is considered to be the most important slicer for sRNA-mediated target cleavages. Previous phylogenetic analysis revealed that AGO1, AGO5, and AGO10 belonged to the same clade in Arabidopsis (Arabidopsis thaliana). In addition, experimental evidence pointed to the possibility that AGO2, AGO7, and AGO10 were implicated in specific miRNA-mediated regulatory pathways. To investigate the potential slicer activities of AGO2, AGO5, AGO7, and AGO10, a comprehensive search was performed for the sRNAs enriched in the four AGO proteins in Arabidopsis. A total of 3 499, 1 618, 4 632, and 63 sRNAs are enriched in AGO2, AGO5, AGO7, and AGO10, respectively. Interestingly, several miRNAs were found to be enriched in AGO2 or AGO5. Transcriptome-wide target identification based on degradome sequencing data uncovered that a number of sRNAs enriched in the four AGOs could perform target cleavages like AGO1-associated miRNAs in plants. Based on the above results, the opinion was put forward that not only AGO1, but also AGO2, AGO5, AGO7, and AGO10 might be essential for the sRNA-mediated regulation of gene expression in plants.

The functional scope of plant microRNA-mediated silencing

Deep sequencing has identified a complex set of plant miRNAs that potentially regulates many target genes of high complementarity. Furthermore, the discovery that many plant miRNAs work through a translational repression mechanism, along with the identification of noncanonical targets, has encouraged bioinformatic searches with less stringent parameters, identifying an even wider range of potential targets. Together, these findings suggest that any given plant miRNA family may regulate a highly diverse set of mRNAs. Here we present evolutionary, genetic, and mechanistic evidence that opposes this idea but instead suggests that families of sequence-related miRNAs regulate very few functionally related targets. We propose that complexities beyond complementarity impact plant miRNA target recognition, possibly explaining the current disparity between bioinformatic prediction and functional evidence.

Identification of SMG6 cleavage sites and a preferred RNA cleavage motif by global analysis of endogenous NMD targets in human cells

In metazoans, cleavage by the endoribonuclease SMG6 is often the first degradative event in non-sense-mediated mRNA decay (NMD). However, the exact sites of SMG6 cleavage have yet to be determined for any endogenous targets, and most evidence as to the identity of SMG6 substrates is indirect. Here, we use Parallel Analysis of RNA Ends to specifically identify the 5′ termini of decay intermediates whose production is dependent on SMG6 and the universal NMD factor UPF1. In this manner, the SMG6 cleavage sites in hundreds of endogenous NMD targets in human cells have been mapped at high resolution. In addition, a preferred sequence motif spanning most SMG6 cleavage sites has been discovered and validated by mutational analysis. For many SMG6 substrates, depletion of SMG6 resulted in the accumulation of decapped transcripts, an effect indicative of competition between SMG6-dependent and SMG6-independent NMD pathways. These findings provide key insights into the mechanisms by which mRNAs targeted by NMD are degraded.

Repression of cell proliferation by miR319-regulated TCP4

Leaf development has been extensively studied on a genetic level. However, little is known about the interplay between the developmental regulators and the cell cycle machinery--a link that ultimately affects leaf form and size. miR319 is a conserved microRNA that regulates TCP transcription factors involved in multiple developmental pathways, including leaf development and senescence, organ curvature, and hormone biosynthesis and signaling. Here, we analyze the participation of TCP4 in the control of cell proliferation. A small increase in TCP4 activity has an immediate impact on leaf cell number, by significantly reducing cell proliferation. Plants with high TCP4 levels have a strong reduction in the expression of genes known to be active in G2-M phase of the cell cycle. Part of these effects is mediated by induction of miR396, which represses Growth-Regulating Factor (GRF) transcription factors. Detailed analysis revealed TCP4 to be a direct regulator of MIR396b. However, we found that TCP4 can control cell proliferation through additional pathways, and we identified a direct connection between TCP4 and ICK1/KRP1, a gene involved in the progression of the cell cycle. Our results show that TCP4 can activate different pathways that repress cell proliferation.

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets, including new miRNA target-identification software

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.

miR396a-Mediated basic helix-loop-helix transcription factor bHLH74 repression acts as a regulator for root growth in Arabidopsis seedlings

miR396 targets seven GROWTH-REGULATING FACTOR (GRF) genes and the BASIC HELIX-LOOP-HELIX (bHLH) TRANSCRIPTION FACTOR 74 gene (bHLH74) in Arabidopsis. Previous research revealed that the miR396 target module regulates cell proliferation and plays a critical role in leaf development. However, no additional biological functions of miR396 have been investigated in detail. In this study, T-DNA insertion mutants and transgenic plants with altered levels of miR396 or its target genes were used to characterize the regulatory role of miR396 in root development. We found that AtMIR396a was the predominant source for miR396 accumulation in the roots of seedlings, and that the mir396a-1 mutant had longer roots than wild-type seedlings. Overexpression of AtMIR396a decreased the transcript levels of target genes such as GRF genes and bHLH74, and resulted in a shorter root phenotype. Furthermore, the bhlh74-1 mutant had shorter roots, whereas overexpression of an miR396-resistant form of bHLH74 (mbHLH74) had an enhanced root growth phenotype. Moreover, MIR396a regulated root growth by affecting the elongation zone. Taken together, these data indicate that miR396a-mediated bHLH74 repression helps regulate root growth in Arabidopsis seedlings.

miR-PREFeR: an accurate, fast and easy-to-use plant miRNA prediction tool using small RNA-Seq data

SUMMARY

Plant microRNA prediction tools that use small RNA-sequencing data are emerging quickly. These existing tools have at least one of the following problems: (i) high false-positive rate; (ii) long running time; (iii) work only for genomes in their databases; (iv) hard to install or use. We developed miR-PREFeR (miRNA PREdiction From small RNA-Seq data), which uses expression patterns of miRNA and follows the criteria for plant microRNA annotation to accurately predict plant miRNAs from one or more small RNA-Seq data samples of the same species. We tested miR-PREFeR on several plant species. The results show that miR-PREFeR is sensitive, accurate, fast and has low-memory footprint.

AVAILABILITY AND IMPLEMENTATION

https://github.com/hangelwen/miR-PREFeR

Sample sequencing of vascular plants demonstrates widespread conservation and divergence of microRNAs

Small RNAs are pivotal regulators of gene expression that guide transcriptional and post-transcriptional silencing mechanisms in eukaryotes, including plants. Here we report a comprehensive atlas of sRNA and miRNA from 3 species of algae and 31 representative species across vascular plants, including non-model plants. We sequence and quantify sRNAs from 99 different tissues or treatments across species, resulting in a data set of over 132 million distinct sequences. Using miRBase mature sequences as a reference, we identify the miRNA sequences present in these libraries. We apply diverse profiling methods to examine critical sRNA and miRNA features, such as size distribution, tissue-specific regulation and sequence conservation between species, as well as to predict putative new miRNA sequences. We also develop database resources, computational analysis tools and a dedicated website, http://smallrna.udel.edu/. This study provides new insights on plant sRNAs and miRNAs, and a foundation for future studies.

Seeing the forest for the trees: annotating small RNA producing genes in plants

A key goal in genomics is the complete annotation of the expressed regions of the genome. In plants, substantial portions of the genome make regulatory small RNAs produced by Dicer-Like (DCL) proteins and utilized by Argonaute (AGO) proteins. These include miRNAs and various types of endogenous siRNAs. Small RNA-seq, enabled by cheap and fast DNA sequencing, has produced an enormous volume of data on plant miRNA and siRNA expression in recent years. In this review, we discuss recent progress in using small RNA-seq data to produce stable and reliable annotations of miRNA and siRNA genes in plants. In addition, we highlight key goals for the future of small RNA gene annotation in plants.