Unique expression, processing regulation, and regulatory network of peach (Prunus persica) miRNAs

miRNA and Degradome Sequencing Reveal miRNA and Their Target Genes That May Mediate Shoot Growth in Spur Type Mutant "Yanfu 6"

The spur-type growth habit in apple trees is characterized by short internodes, increased number of fruiting spurs, and compact growth that promotes flowering and facilitates management practices, such as pruning. The molecular mechanisms responsible for regulating spur-type growth have not been elucidated. In the present study, miRNAs and the expression of their potential target genes were evaluated in shoot tips of "Nagafu 2" (CF) and spur-type bud mutation "Yanfu 6" (YF). A total of 700 mature miRNAs were identified, including 202 known apple miRNAs and 498 potential novel miRNA candidates. A comparison of miRNA expression in CF and YF revealed 135 differentially expressed genes, most of which were downregulated in YF. YF also had lower levels of GA, ZR, IAA, and ABA hormones, relative to CF. Exogenous applications of GA promoted YF shoot growth. Based on the obtained results, a regulatory network involving plant hormones, miRNA, and their potential target genes is proposed for the molecular mechanism regulating the growth of YF. miRNA164, miRNA166, miRNA171, and their potential targets, and associated plant hormones, appear to regulate shoot apical meristem (SAM) growth. miRNA159, miRNA167, miRNA396, and their potential targets, and associated plant hormones appear to regulate cell division and internode length. This study provides a foundation for further studies designed to elucidate the mechanism underlying spur-type apple architecture.

Biogenesis of diverse plant phasiRNAs involves an miRNA-trigger and Dicer-processing

It has been almost 30 years since RNA interference (RNAi) was shown to silence genes via double-stranded RNAs (dsRNAs) in Caenorhabditis elegans (Fire et al. 1998). 20-30-nucleotide (nt) small non-coding RNAs are a key element of the RNAi machinery. Recently, phased small interfering RNAs (phasiRNAs), small RNAs that are generated from a long RNA precursor at intervals of 21 to 26-nt, have been identified in plants and animals. In Drosophila, phasiRNAs are generated by the endonuclease, Zucchini (Zuc), in germlines. These phasiRNAs, known as one of PIWI-interacting RNAs (piRNAs), mainly repress transposable elements. Similarly, reproduction-specific phasiRNAs have been identified in the family Poaceae, although DICER LIKE (DCL) protein-dependent phasiRNA biogenesis in rice is distinct from piRNA biogenesis in animals. In plants, phasiRNA biogenesis is initiated when 22-nt microRNAs (miRNAs) cleave single-stranded target RNAs. Subsequently, RNA-dependent RNA polymerase (RDR) forms dsRNAs from the cleaved RNAs, and dsRNAs are further processed by DCLs into 21 to 24-nt phasiRNAs. Finally, the phasiRNAs are loaded to ARGONAUTE (AGO) proteins to induce RNA-silencing. There are diverse types of phasiRNA precursors and the miRNAs that trigger the biogenesis. Their expression patterns also differ among plant species, suggesting that species-specific combinations of these triggers dictate the spatio-temporal pattern of phasiRNA biogenesis during development, or in response to environmental stimuli.

Identification of anthocyanin biosynthesis related microRNAs in a distinctive Chinese radish (Raphanus sativus L.) by high-throughput sequencing

Anthocyanins are widely distributed water-soluble phytochemical pigments belonging to the flavonoid group. To date, limited knowledge is available about the regulatory roles of miRNAs in anthocyanin biosynthesis in plants. To identify the miRNAs associated with anthocyanin biosynthesis in radish, five small RNA (sRNA) libraries constructed from 'Xinlimei' radish roots at 11, 21, 44, 56 and 73 days (d) were examined using high-throughput sequencing technology. A total of 102.02 million (M) clean reads were generated, from which 483 known and 1415 novel miRNAs were identified. Combined with target prediction and annotation, 72 differentially expressed miRNAs (52 known and 20 novel miRNAs) were more likely to participate in anthocyanin biosynthesis. Several target genes for these miRNAs encode a few transcription factors, including Myb domain (MYB), basic helix-loop-helix (bHLH), WD40 repeat, squamosa promoter binding protein like (SPL), auxin response factor (ARF), ethylene insensitive 3 (EIN3), WRKY and MADS-box proteins. Furthermore, the expression patterns of some anthocyanin biosynthesis related miRNAs and their corresponding targets were validated by RT-qPCR. Based on the characterization of anthocyanin biosynthesis related miRNAs and their target genes, a putative miRNA-target module regulating anthocyanin biosynthesis was proposed. This study represents the first genome-wide identification of miRNAs associated with anthocyanin biosynthesis in radish, and provides insights into the molecular mechanisms underlying regulation of anthocyanin biosynthesis in radish and other crops.

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

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

Repression of MYBL2 by Both microRNA858a and HY5 Leads to the Activation of Anthocyanin Biosynthetic Pathway in Arabidopsis

Extensive studies in various plants show that the anthocyanin biosynthetic process is affected by environmental factors and regulated by many transcription factors through sophisticated regulatory networks. However, it remains largely unclear about the roles of microRNA in this process. Here, we demonstrate that miR858a is a positive regulator of anthocyanin biosynthesis in Arabidopsis seedlings. Overexpression of miR858a enhances the accumulation of anthocyanins, whereas the reduced miR858a activity results in low levels of anthocyanins in STTM858 transgenic plants. We found that miR858a inhibits the expression of MYBL2, a key negative regulator of anthocyanin biosynthesis, by translational repression. In addition, ELONGATED HYPOCOTYL 5 (HY5) was shown to directly bind the MYBL2 promoter and represses its expression via specific histone modifications. Interestingly, we found that miR858a exhibits light-responsive expression in an HY5-dependent manner. Together, these results delineate the HY5-MIR858a-MYBL2 loop as a cellular mechanism for modulating anthocyanin biosynthesis, suggesting that integration of transcriptional and posttranscriptional regulation is critical for governing proper anthocyanin accumulation in response to light and other environmental factors.

Identification of Known and Novel microRNAs and Their Targets in Peach (Prunus persica) Fruit by High-Throughput Sequencing

MicroRNAs (miRNAs) are a group of non-coding RNAs that have functions in post-transcriptional gene regulation in plants. Although the most important economic component of peach trees (Prunus persica) is the fruit, not much is known about miRNAs in this organ. In this study, miRNAs and their targets were identified and characterized from libraries of small RNAs of peach fruit through Solexa based-sequencing and bioinformatics approaches. A total of 557 known peach miRNAs belonging to 34 miRNA families were identified, and some of these miRNAs were found to be highly conserved in at least four other plant species. Using the most current criteria for miRNA annotation, 275 putative novel miRNAs were predicted, and the sequencing frequencies of these novel miRNAs were less than those of the conserved miRNAs. In total, 3959 and 1614 target genes for 349 known and 193 novel miRNAs, respectively, were predicted with the criteria that a single target gene can be targeted by different miRNAs and that a single miRNA can also have a large number of target genes. Three targets were even found to be targeted by 13 novel miRNAs that contained the same complete miRNA sequence at different locations and had different scaffolds. The proteins predicted to be targeted by the miRNAs identified in this study encompass a wide range of transcription factors and are involved in many biological processes and pathways, including development, metabolism, stress responses and signal transduction. A total of 115 and 101 target genes were identified to be cleaved by 60 known miRNAs and 27 novel miRNAs through degradome sequencing, respectively. These miRNAs induce cleavage of their targets precisely at the position between nucleotides 10 and 11 of the miRNA sequences from the 5’ to the 3’ end. Thirty conserved miRNAs and 19 novel miRNAs exhibited differential expression profiles in the peach, and the expression patterns of some miRNAs appeared to be tissue- or developmental stage-specific. The findings of this study provide an important basis for the analysis of miRNAs, their targets and the functions of these targets in peach fruit.

A naïve Bayesian classifier for identifying plant microRNAs

MicroRNAs (miRNAs) are important regulatory molecules in eukaryotic organisms. Existing methods for the identification of mature miRNA sequences in plants rely extensively on the search for stem-loop structures, leading to high false negative rates. Here, we describe a probabilistic method for ranking putative plant miRNAs using a naïve Bayes classifier and its publicly available implementation. We use a number of properties to construct the classifier, including sequence length, number of observations, existence of detectable predicted miRNA* sequences, the distribution of nearby reads and mapping multiplicity. We apply the method to small RNA sequence data from soybean, peach, Arabidopsis and rice and provide experimental validation of several predictions in soybean. The approach performs well overall and strongly enriches for known miRNAs over other types of sequences. By utilizing a Bayesian approach to rank putative miRNAs, our method is able to score miRNAs that would be eliminated by other methods, such as those that have low counts or lack detectable miRNA* sequences. As a result, we are able to detect several soybean miRNA candidates, including some that are 24 nucleotides long, a class that is almost universally eliminated by other methods.

In silico search and biological validation of microRNAs related to drought response in peach and almond

Plant responses to drought stress are regulated at the transcriptional and post-transcriptional levels through noncoding endogenous microRNAs. These microRNAs play key roles in gene expression, mainly by down-regulating target mRNAs. In this work, an in silico search and validation for microRNAs related to drought response in peach ('G.H. Hill'), almond ('Sefied') and an interspecific peach-almond hybrid ('GN 15') has been performed. We used qPCR to analyse the gene expression of several miRNAs described as being related to drought response in peach, including miR156, miR159, miR160, miR167, miR171, miR172, miR398, miR403, miR408, miR842 and miR2275 under mild and severe water deficit. These miRNAs were in silico selected on the basis of previous works, their conservation in plants and their drought response. qPCR analysis confirmed the implication of these miRNAs in the dehydration stress response in the three assayed genotypes. Comparison of miRNA expression patterns in the three evaluated genotypes indicated that the hybrid GN 15 showed higher expression levels of specific miRNAs which should be related to the observed drought tolerance. mRNA target transcripts of the miRNAs studied were predicted using the Rose database, which includes transcription factors that regulate plant growth and development. In addition, results showed that the promoter region contains responsive elements to hormone-mediated regulatory elements. Network analysis not only unravelled the interaction between miRNAs and their predicted gene targets but also highlighted the roles of miRNAs in response to drought stress.

Citrus psorosis virus 24K protein interacts with citrus miRNA precursors, affects their processing and subsequent miRNA accumulation and target expression

Sweet orange (Citrus sinensis), one of the most important fruit crops worldwide, may suffer from disease symptoms induced by virus infections, thus resulting in dramatic economic losses. Here, we show that the infection of sweet orange plants with two isolates of Citrus psorosis virus (CPsV) expressing different symptomatology alters the accumulation of a set of endogenous microRNAs (miRNAs). Within these miRNAs, miR156, miR167 and miR171 were the most down-regulated, with almost a three-fold reduction in infected samples. This down-regulation led to a concomitant up-regulation of some of their targets, such as Squamosa promoter-binding protein-like 9 and 13, as well as Scarecrow-like 6. The processing of miRNA precursors, pre-miR156 and pre-miR171, in sweet orange seems to be affected by the virus. For instance, virus infection increases the level of unprocessed precursors, which is accompanied by a concomitant decrease in mature species accumulation. miR156a primary transcript accumulation remained unaltered, thus strongly suggesting a processing deregulation for this transcript. The co-immunoprecipitation of viral 24K protein with pre-miR156a or pre-miR171a suggests that the alteration in the processing of these precursors might be caused by a direct or indirect interaction with this particular viral protein. This result is also consistent with the nuclear localization of both miRNA precursors and the CPsV 24K protein. This study contributes to the understanding of the manner in which a virus can alter host regulatory mechanisms, particularly miRNA biogenesis and target expression.

Small RNAs Add Zing to the Zig-Zag-Zig Model of Plant Defenses

Plant small RNAs play important roles in transcriptional and posttranscriptional regulation, with ongoing work demonstrating their functions in diverse pathways. Their roles in defense responses are a topic of active investigation, particularly the rich set of micro (mi)RNAs that target disease resistance genes such as nucleotide binding/leucine-rich repeat (NB-LRR) genes. The miRNA-NB-LRR interactions result in the production of phased, secondary small interfering (phasi)RNAs, and phasiRNAs function in both cis and trans to propagate negative regulatory effects across additional members of the target gene family. Yet, while phasiRNAs have the capacity to trigger targeted decay of specific targets, both in cis and trans, their functional relevance in NB-LRR regulation remains largely a matter of speculation.

Identification and expression analysis of genes related to calyx persistence in Korla fragrant pear

BACKGROUND

The objective of this study was to increase understanding about genetic mechanisms affecting calyx persistence in Korla fragrant pear (Pyrus brestschneideri Rehd). Flowers were collected at early bloom, full bloom, and late bloom. The RNA was extracted from the flowers and then combined according to calyx type. Transcriptome and digital gene expression (DGE) profiles of flowers, ovaries, and sepals with persistent calyx (SC_hua, SC_ep, and SC_zf, respectively) were compared with those of flowers, ovaries, and sepals with deciduous calyx (TL_hua, TL_ep, and TL_zf, respectively). Temporal changes in the expression of selected genes in floral organs with either persistent or deciduous calyx were compared using real-time quantitative PCR (qRT-PCR).

RESULTS

Comparison of the transcriptome sequences for SC_hua and TL_hua indicated 26 differentially expressed genes (DEGs) with known relationship to abscission and 10 DEGs with unknown function. We identified 98 MYB and 21 SPL genes from the assembled unigenes. From SC_zf vs TL_zf, we identified 21 DEGs with known relationship to abscission and 18 DEGs with unknown function. From SC_ep vs TL_ep, 12 DEGs with known relationship to abscission were identified along with 11 DEGs with unknown function. Ten DEGs were identified by both transcriptome sequencing and DGE sequencing.

CONCLUSIONS

More than 50 DEGs were observed that were related to calyx persistence in Korla fragrant pear. Some of the genes were related to cell wall degradation, plant hormone signal transduction, and stress response. Other DEGs were identified as zinc finger protein genes and lipid transfer protein genes. Further analysis showed that calyx persistence in Korla fragment pear was a metabolic process regulated by many genes related to cell wall degradation and plant hormones.

A highly specific microRNA-mediated mechanism silences LTR retrotransposons of strawberry

Small RNAs are involved in a plethora of functions in plant genomes. In general, transcriptional gene silencing is mediated by 24-nucleotide siRNAs and is required for maintaining transposable elements in a silenced state. However, microRNAs are not commonly associated with transposon silencing. In this study, we performed small RNA transcriptome and degradome analyses of the Rosaceae model plant Fragaria vesca (the woodland strawberry) at the genome-wide level, and identified miRNA families and their targets. We report a highly specific mechanism of LTR retrotransposon silencing mediated by an abundant, ubiquitously expressed miRNA (fve-miR1511) generated from a single locus. This miRNA specifically targets LTR retroelements, silencing them post-transcriptionally by perfectly pairing to the highly conserved primer binding site for methionyl initiator tRNA that is essential for reverse transcription. We investigated the possible origins of this miRNA, and present evidence that the pre-miR1511 hairpin structure probably derived from a locus coding for tRNA(iM) (et) through a single microinversion event. Our study shows that this miRNA targets retrotransposons specifically and constitutively, and contributes to features such as genome stability, size and architecture in a far more direct way than previously thought.

Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae

Background

Micro RNAs are a class of small non coding RNAs of 20–24 nucleotides transcribed as single stranded precursors from MIR gene loci. Initially described as post-transcriptional regulators involved in development, two decades ago, miRNAs have been proven to regulate a wide range of processes in plants such as germination, morphology and responses to biotic and abiotic stress. Despite wide conservation in plants, a number of miRNAs are lineage specific. We describe the first genome wide survey of Eucalyptus miRNAs based on high throughput sequencing.

Results

In addition to discovering small RNA sequences, MIR loci were mapped onto the reference genome and interspecific variability investigated. Sequencing was carried out for the two most world widely planted species, E. grandis and E. globulus. To maximize discovery, E. grandis samples were from BRASUZ1, the same tree whose genome provided the reference sequence. Interspecific analysis reinforces the variability in small RNA repertoire even between closely related species. Characterization of Eucalyptus small RNA sequences showed 95 orthologous to conserved miRNAs and 193 novel miRNAs. In silico target prediction confirmed 163 novel miRNAs and degradome sequencing experimentally confirmed several hundred targets. Experimental evidence based on the exclusive expression of a set of small RNAs across 16 species within Myrtaceae further highlighted variable patterns of conservation and diversity of these regulatory elements.

Conclusions

The description of miRNAs in Eucalyptus contributes to scientific knowledge of this vast genre, which is the most widely planted hardwood crop in the tropical and subtropical world, adding another important element to the annotation of Eucalyptus grandis reference genome.

Electronic supplementary material

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

Transition from two to one integument in Prunus species: expression pattern of INNER NO OUTER (INO), ABERRANT TESTA SHAPE (ATS) and ETTIN (ETT)

While gymnosperm ovules have one integument, in most angiosperms two integuments surround the ovules. Unitegmic ovules have arisen independently several times during the evolution of angiosperms, but the ultimate genetic cause of the presence of a single integument remains elusive. We compared species of the genus Prunus that have different numbers of integuments: bitegmic species, such as Prunus armeniaca (apricot) and Prunus persica (peach), and unitegmic species, such as Prunus incisa, analyzing the expression pattern of genes that are involved in integument development in Arabidopsis thaliana: INNER NO OUTER (INO), ABERRANT TESTA SHAPE (ATS) and ETTIN (ETT). Bitegmic and unitegmic species showed similar INO expression patterns, indicative of the conservation of an outer integument. However, expression of ETT, which occurs in the boundary of the outer and inner integuments, was altered in unitegmic ovules, which showed lack of ETT expression. These results strongly suggest that the presence of a single integument could be attributable to the amalgamation of two integuments and support the role of ETT in the fusion of the outer and inner integuments in unitegmic ovules, a situation that could be widespread in other unitegmic species of angiosperms.

Novel and Recently Evolved MicroRNA Clusters Regulate Expansive F-BOX Gene Networks through Phased Small Interfering RNAs in Wild Diploid Strawberry

The wild strawberry (Fragaria vesca) has recently emerged as an excellent model for cultivated strawberry (Fragaria × ananassa) as well as other Rosaceae fruit crops due to its short seed-to-fruit cycle, diploidy, and sequenced genome. Deep sequencing and parallel analysis of RNA ends were used to identify F. vesca microRNAs (miRNAs) and their target genes, respectively. Thirty-eight novel and 31 known miRNAs were identified. Many known miRNAs targeted not only conserved mRNA targets but also developed new target genes in F. vesca. Significantly, two new clusters of miRNAs were found to collectively target 94 F-BOX (FBX) genes. One of the miRNAs in the new cluster is 22 nucleotides and triggers phased small interfering RNA production from six FBX genes, which amplifies the silencing to additional FBX genes. Comparative genomics revealed that the main novel miRNA cluster evolved from duplications of FBX genes. Finally, conserved trans-acting siRNA pathways were characterized and confirmed with distinct features. Our work identified novel miRNA-FBX networks in F. vesca and shed light on the evolution of miRNAs/phased small interfering RNA networks that regulate large gene families in higher plants.

Differential mRNA Accumulation upon Early Arabidopsis thaliana Infection with ORMV and TMV-Cg Is Associated with Distinct Endogenous Small RNAs Level

Small RNAs (sRNAs) play important roles in plant development and host-pathogen interactions. Several studies have highlighted the relationship between viral infections, endogenous sRNA accumulation and transcriptional changes associated with symptoms. However, few studies have described a global analysis of endogenous sRNAs by comparing related viruses at early stages of infection, especially before viral accumulation reaches systemic tissues. An sRNA high-throughput sequencing of Arabidopsis thaliana leaf samples infected either with Oilseed rape mosaic virus (ORMV) or crucifer-infecting Tobacco mosaic virus (TMV-Cg) with slightly different symptomatology at two early stages of infection (2 and 4dpi) was performed. At early stages, both viral infections strongly alter the patterns of several types of endogenous sRNA species in distal tissues with no virus accumulation suggesting a systemic signaling process foregoing to virus spread. A correlation between sRNAs derived from protein coding genes and the associated mRNA transcripts was also detected, indicating that an unknown recursive mechanism is involved in a regulatory circuit encompassing this sRNA/mRNA equilibrium. This work represents the initial step in uncovering how differential accumulation of endogenous sRNAs contributes to explain the massive alteration of the transcriptome associated with plant-virus interactions.

Small tandem target mimic-mediated blockage of microRNA858 induces anthocyanin accumulation in tomato

Our work strongly suggests that microRNA858 regulates anthocyanin biosynthesis in tomato by modulating the expression of two R2R3 MYB transcription factors, underscoring the importance of microRNAs in the gene regulatory network controlling specialized metabolism in plants. The biological functions of microRNA858 (miR858), a recently identified small RNA, are not well understood. Here, we identified miR858 as a negative regulator of anthocyanin biosynthesis in tomato (Solanum lycopersicum). RNA ligase-mediated 5'RACE cleavage assay showed that miR858 mediates the cleavage of SlMYB7-like and SlMYB48-like transcripts in tomato. Expression analysis revealed an inverse correlation between the accumulation of miR858 and its target SlMYB7-like mRNA, in different tissues of tomato. Subsequently, a small tandem target mimic construct for blocking miR858 (STTM858) was generated and transformed into tomato. The majority of endogenous miR858 was blocked in STTM858 over-expressing tomato plants, whereas SlMYB7-like transcripts increased significantly. Concomitantly, upregulated expression was detected for several anthocyanin biosynthetic genes, including PAL, CHS, DFR, ANS and 3GT. As a result, anthocyanins were highly accumulated in young seedlings, leaves, stems and leaf buds of the transgenic plants under normal growth conditions. In addition, over-expression of STTM858 in tomato activated another MYB transcription factor, SlMYB48, implicating the possible involvement of SlMYB48 in anthocyanin biosynthesis.

Prunus transcription factors: breeding perspectives

Many plant processes depend on differential gene expression, which is generally controlled by complex proteins called transcription factors (TFs). In peach, 1533 TFs have been identified, accounting for about 5.5% of the 27,852 protein-coding genes. These TFs are the reference for the rest of the Prunus species. TF studies in Prunus have been performed on the gene expression analysis of different agronomic traits, including control of the flowering process, fruit quality, and biotic and abiotic stress resistance. These studies, using quantitative RT-PCR, have mainly been performed in peach, and to a lesser extent in other species, including almond, apricot, black cherry, Fuji cherry, Japanese apricot, plum, and sour and sweet cherry. Other tools have also been used in TF studies, including cDNA-AFLP, LC-ESI-MS, RNA, and DNA blotting or mapping. More recently, new tools assayed include microarray and high-throughput DNA sequencing (DNA-Seq) and RNA sequencing (RNA-Seq). New functional genomics opportunities include genome resequencing and the well-known synteny among Prunus genomes and transcriptomes. These new functional studies should be applied in breeding programs in the development of molecular markers. With the genome sequences available, some strategies that have been used in model systems (such as SNP genotyping assays and genotyping-by-sequencing) may be applicable in the functional analysis of Prunus TFs as well. In addition, the knowledge of the gene functions and position in the peach reference genome of the TFs represents an additional advantage. These facts could greatly facilitate the isolation of genes via QTL (quantitative trait loci) map-based cloning in the different Prunus species, following the association of these TFs with the identified QTLs using the peach reference genome.

A dynamic evolutionary and functional landscape of plant phased small interfering RNAs

Background

Secondary, phased small interfering RNAs (phasiRNAs) derived from protein-coding or noncoding loci (PHAS) are emerging as a new type of regulators of gene expression in plants. However, the evolution and function of these novel siRNAs in plant species remain largely unexplored.

Results

We systematically analyzed PHAS loci in 23 plant species covering major phylogenetic groups spanning alga, moss, gymnosperm, basal angiosperm, monocot, and dicot. We identified over 3,300 PHAS loci, among which ~1,600 were protein-coding genes. Most of these PHAS loci were novel and clade- or species-specific and showed distinct expression patterns in association with particular development stages, viral infection, or abiotic stresses. Unexpectedly, numerous PHAS loci produced phasiRNAs from introns or exon–intron junction regions. Our comprehensive analysis suggests that phasiRNAs predominantly regulate protein-coding genes from which they are derived and genes from the same families of the phasiRNA-deriving genes, in contrast to the dominant trans-regulatory mode of miRNAs. The stochastic occurrence of many PHAS loci in the plant kingdom suggests their young evolutionary origins.

Conclusions

Our study discovered an unprecedented diversity of protein-coding genes that produce phasiRNAs in a wide variety of plants, and set a kingdom-wide foundation for investigating the novel roles of phasiRNAs in shaping phenotype diversities of plants.

Electronic supplementary material

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

The liverwort Pellia endiviifolia shares microtranscriptomic traits that are common to green algae and land plants

Liverworts are the most basal group of extant land plants. Nonetheless, the molecular biology of liverworts is poorly understood. Gene expression has been studied in only one species, Marchantia polymorpha. In particular, no microRNA (miRNA) sequences from liverworts have been reported. Here, Illumina-based next-generation sequencing was employed to identify small RNAs, and analyze the transcriptome and the degradome of Pellia endiviifolia. Three hundred and eleven conserved miRNA plant families were identified, and 42 new liverwort-specific miRNAs were discovered. The RNA degradome analysis revealed that target mRNAs of only three miRNAs (miR160, miR166, and miR408) have been conserved between liverworts and other land plants. New targets were identified for the remaining conserved miRNAs. Moreover, the analysis of the degradome permitted the identification of targets for 13 novel liverwort-specific miRNAs. Interestingly, three of the liverwort microRNAs show high similarity to previously reported miRNAs from Chlamydomonas reinhardtii. This is the first observation of miRNAs that exist both in a representative alga and in the liverwort P. endiviifolia but are not present in land plants. The results of the analysis of the P. endivifolia microtranscriptome support the conclusions of previous studies that placed liverworts at the root of the land plant evolutionary tree of life.

Identification of novel miRNAs and miRNA expression profiling in wheat hybrid necrosis

MicroRNAs (miRNAs) play essential roles in a vast array of biological processes, including growth and development, defense against viral infection, and responses to environmental changes in plant. Wheat hybrid necrosis is an interesting genetic phenomenon observed frequency and it is lethal or semi lethal, resulting in gradual death or loss of productivity. However, the molecular basis and mechanisms associated with hybrid necrosis in wheat are still not well understood. Here, we report the population and expression profiles of miRNAs in wheat hybrid necrosis. We identified a total of 57 conserved miRNA families as well as 182 putative novel miRNAs. Expression profiling revealed that expression of 49 known miRNAs and 165 novel miRNAs was changed in hybrid necrosis. And the expression levels of some miRNAs and their predicated targets have been confirmed by qRT-PCR. These results indicate that these miRNAs, especially miR159, miR166, miR167 and miR5072 could be involved in the extensive regulation of gene expression in response to hybrid necrosis.

Predicted trans-acting siRNAs in the human brain

Endogenous small non-coding RNAs play pivotal roles in regulating gene expression in eukaryotes. Many studies have investigated the function and molecular mechanism of microRNAs in the development and disease of various organisms via mRNA repression of protein-coding genes. Recent findings indicate microRNAs might trigger the generation of trans-acting small interfering RNAs (ta-siRNAs). The interaction among different types of small RNA molecules reveals an even more complicated and elaborate pattern of RNA regulation during gene expression than previously thought. We developed a method for mining ta-siRNA sequences and evaluated the performance of our novel method using data from Arabidopsis thaliana. Additionally, using small RNA and degradome data for the human brain, we identified 155 small RNAs that satisfied ta-siRNA characteristics. The DRAXIN and ATCAY genes, which are preferentially expressed in the human brain, were predicted to be the targets of 12 potential ta-siRNAs.

Genome-wide identification of vegetative phase transition-associated microRNAs and target predictions using degradome sequencing in Malus hupehensis

Background

A long juvenile period between germination and flowering is a common characteristic among fruit trees, including Malus hupehensis (Pamp.) Rehd., which is an apple rootstock widely used in China. microRNAs (miRNAs) play an important role in the regulation of phase transition and reproductive growth processes.

Results

M. hupehensis RNA libraries, one adult and one juvenile phase, were constructed using tree leaves and underwent high-throughput sequencing. We identified 42 known miRNA families and 172 novel miRNAs. We also identified 127 targets for 25 known miRNA families and 168 targets for 35 unique novel miRNAs using degradome sequencing. The identified miRNA targets were categorized into 58 biological processes, and the 123 targets of known miRNAs were associated with phase transition processes. The KEGG analysis revealed that these targets were involved in starch and sucrose metabolism, and plant hormone signal transduction. Expression profiling of miRNAs and their targets indicated multiple regulatory functions in the phase transition. The higher expression level of mdm-miR156 and lower expression level of mdm-miR172 in the juvenile phase leaves implied that these two small miRNAs regulated the phase transition. mdm-miR160 and miRNA393, which regulate genes involved in auxin signal transduction, could also be involved in controlling this process. The identification of known and novel miRNAs and their targets provides new information on this regulatory process in M. hupehensis, which will contribute to the understanding of miRNA functions during growth, phase transition and reproduction in woody fruit trees.

Conclusions

The combination of sRNA and degradome sequencing can be used to better illustrate the profiling of hormone-regulated miRNAs and miRNA targets involving complex regulatory networks, which will contribute to the understanding of miRNA functions during growth, phase transition and reproductive growth in perennial woody fruit trees.

Electronic supplementary material

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

Genome-wide discovery and analysis of phased small interfering RNAs in Chinese sacred lotus

Phased small interfering RNA (phasiRNA) generating loci (briefly as PHAS) in plants are a novel class of genes that are normally regulated by microRNAs (miRNAs). Similar to miRNAs, phasiRNAs encoded by PHAS play important regulatory roles by targeting protein coding transcripts in plant species. We performed a genome-wide discovery of PHAS loci in Chinese sacred lotus and identified a total of 106 PHAS loci. Of these, 47 loci generate 21 nucleotide (nt) phasiRNAs and 59 loci generate 24 nt phasiRNAs, respectively. We have also identified a new putative TAS3 and a putative TAS4 loci in the lotus genome. Our results show that some of the nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance proteins and MYB transcription factors potentially generate phasiRNAs. Furthermore, our results suggest that some large subunit (LSU) rRNAs can derive putative phasiRNAs, which is potentially resulted from crosstalk between small RNA biogenesis pathways that are employed to process rRNAs and PHAS loci, respectively. Some of the identified phasiRNAs have putative trans-targets with less than 4 mismatches, suggesting that the identified PHAS are involved in many different pathways. Finally, the discovery of 24 nt PHAS in lotus suggests that there are 24 nt PHAS in dicots.

MicroRNAs in fruit trees: discovery, diversity and future research directions

Since the first description of microRNAs (miRNAs) 20 years ago, the number of miRNAs identified in different eukaryotic organisms has exploded, largely due to the recent advances in DNA sequencing technologies. Functional studies, mostly from model species, have revealed that miRNAs are major post-transcriptional regulators of gene expression in eukaryotes. In plants, they are implicated in fundamental biological processes, from plant development and morphogenesis, to regulation of plant pathogen and abiotic stress responses. Although a substantial number of miRNAs have been identified in fruit trees to date, their functions remain largely uncharacterised. The present review aims to summarise the progress made in miRNA research in fruit trees, focusing specifically on the economically important species Prunus persica, Malus domestica, Citrus spp, and Vitis vinifera. We also discuss future miRNA research prospects in these plants and highlight potential applications of miRNAs in the on-going improvement of fruit trees.

Extending the sRNAome of apple by next-generation sequencing

The global importance of apple as a fruit crop necessitates investigations into molecular aspects of the processes that influence fruit quality and yield, including plant development, fruit ripening and disease resistance. In order to study and understand biological processes it is essential to recognise the range of molecules, which influence these processes. Small non-coding RNAs are regulatory agents involved in diverse plant activities, ranging from development to stress response. The occurrence of these molecules in apple leaves was studied by means of next-generation sequencing. 85 novel microRNA (miRNA) gene loci were predicted and characterized along with known miRNA loci. Both cis- and trans-natural antisense transcript pairs were identified. Although the trans-overlapping regions were enriched in small RNA (sRNA) production, cis-overlaps did not seem to agree. More than 150 phased regions were also identified, and for a small subset of these, potential miRNAs that could initiate phasing, were revealed. Repeat-associated siRNAs, which are generated from repetitive genomic regions such as transposons, were also analysed. For this group almost all available repeat sequences, associated with the apple genome and present in Repbase, were found to produce siRNAs. Results from this study extend our current knowledge on apple sRNAs and their precursors significantly. A rich molecular resource has been created and is available to the research community to serve as a baseline for future studies.

Recent advancements to study flowering time in almond and other Prunus species

Flowering time is an important agronomic trait in almond since it is decisive to avoid the late frosts that affect production in early flowering cultivars. Evaluation of this complex trait is a long process because of the prolonged juvenile period of trees and the influence of environmental conditions affecting gene expression year by year. Consequently, flowering time has to be studied for several years to have statistical significant results. This trait is the result of the interaction between chilling and heat requirements. Flowering time is a polygenic trait with high heritability, although a major gene Late blooming (Lb) was described in "Tardy Nonpareil." Molecular studies at DNA level confirmed this polygenic nature identifying several genome regions (Quantitative Trait Loci, QTL) involved. Studies about regulation of gene expression are scarcer although several transcription factors have been described as responsible for flowering time. From the metabolomic point of view, the integrated analysis of the mechanisms of accumulation of cyanogenic glucosides and flowering regulation through transcription factors open new possibilities in the analysis of this complex trait in almond and in other Prunus species (apricot, cherry, peach, plum). New opportunities are arising from the integration of recent advancements including phenotypic, genetic, genomic, transcriptomic, and metabolomics studies from the beginning of dormancy until flowering.

Genome-wide identification of non-coding RNAs interacted with microRNAs in soybean

A wide range of RNA species interacting with microRNAs (miRNAs) form a complex gene regulation network and play vital roles in diverse biological processes. In this study, we performed a genome-wide identification of endogenous target mimics (eTMs) for miRNAs and phased-siRNA-producing loci (PHAS) in soybean with a focus on those involved in lipid metabolism. The results showed that a large number of eTMs and PHAS genes could be found in soybean. Additionally, we found that lipid metabolism related genes were potentially regulated by 28 miRNAs, and nine of them were potentially further regulated by a number of eTMs with expression evidence. Thirty-three miRNAs were found to trigger production of phasiRNAs from 49 PHAS genes, which were able to target lipid metabolism related genes. Degradome data supported miRNA- and/or phasiRNA-mediated cleavage of genes involved in lipid metabolism. Most eTMs for miRNAs involved in lipid metabolism and phasiRNAs targeting lipid metabolism related genes showed a tissue-specific expression pattern. Our bioinformatical evidences suggested that lipid metabolism in soybean is potentially regulated by a complex non-coding network, including miRNAs, eTMs, and phasiRNAs, and the results extended our knowledge on functions of non-coding RNAs.

Parallel analysis of RNA ends enhances global investigation of microRNAs and target RNAs of Brachypodium distachyon

BACKGROUND

The wild grass Brachypodium distachyon has emerged as a model system for temperate grasses and biofuel plants. However, the global analysis of miRNAs, molecules known to be key for eukaryotic gene regulation, has been limited in B. distachyon to studies examining a few samples or that rely on computational predictions. Similarly an in-depth global analysis of miRNA-mediated target cleavage using parallel analysis of RNA ends (PARE) data is lacking in B. distachyon.

RESULTS

B. distachyon small RNAs were cloned and deeply sequenced from 17 libraries that represent different tissues and stresses. Using a computational pipeline, we identified 116 miRNAs including not only conserved miRNAs that have not been reported in B. distachyon, but also non-conserved miRNAs that were not found in other plants. To investigate miRNA-mediated cleavage function, four PARE libraries were constructed from key tissues and sequenced to a total depth of approximately 70 million sequences. The roughly 5 million distinct genome-matched sequences that resulted represent an extensive dataset for analyzing small RNA-guided cleavage events. Analysis of the PARE and miRNA data provided experimental evidence for miRNA-mediated cleavage of 264 sites in predicted miRNA targets. In addition, PARE analysis revealed that differentially expressed miRNAs in the same family guide specific target RNA cleavage in a correspondingly tissue-preferential manner.

CONCLUSIONS

B. distachyon miRNAs and target RNAs were experimentally identified and analyzed. Knowledge gained from this study should provide insights into the roles of miRNAs and the regulation of their targets in B. distachyon and related plants.

Identification of miRNAs and their target genes in peach (Prunus persica L.) using high-throughput sequencing and degradome analysis

MicroRNAs play critical roles in various biological and metabolic processes. The function of miRNAs has been widely studied in model plants such as Arabidopsis and rice. However, the number of identified miRNAs and related miRNA targets in peach (Prunus persica) is limited. To understand further the relationship between miRNAs and their target genes during tissue development in peach, a small RNA library and three degradome libraries were constructed from three tissues for deep sequencing. We identified 117 conserved miRNAs and 186 novel miRNA candidates in peach by deep sequencing and 19 conserved miRNAs and 13 novel miRNAs were further evaluated for their expression by RT-qPCR. The number of gene targets that were identified for 26 conserved miRNA families and 38 novel miRNA candidates, were 172 and 87, respectively. Some of the identified miRNA targets were abundantly represented as conserved miRNA targets in plant. However, some of them were first identified and showed important roles in peach development. Our study provides information concerning the regulatory network of miRNAs in peach and advances our understanding of miRNA functions during tissue development.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

Computational identification of conserved microRNAs and their putative targets in the Hypericum perforatum L. flower transcriptome

MicroRNAs (miRNAs) have recently emerged as important regulators of gene expression in plants. Many miRNA families and their targets have been extensively studied in model species and major crops. We have characterized mature miRNAs along with their precursors and potential targets in Hypericum to generate a comprehensive list of conserved miRNA families and to investigate the regulatory role of selected miRNAs in biological processes that occur in the flower. St. John's wort (Hypericum perforatum L., 2n = 4x = 32), a medicinal plant that produces pharmaceutically important metabolites with therapeutic activities, was chosen because it is regarded as an attractive model system for the study of apomixis. A computational in silico prediction of structure, in combination with an in vitro validation, allowed us to identify 7 pre-miRNAs, including miR156, miR166, miR390, miR394, miR396, and miR414. We demonstrated that H. perforatum flowers share highly conserved miRNAs and that these miRNAs potentially target dozens of genes with a wide range of molecular functions, including metabolism, response to stress, flower development, and plant reproduction. Our analysis paves the way toward identifying flower-specific miRNAs that may differentiate the sexual and apomictic reproductive pathways.

Phased, secondary, small interfering RNAs in posttranscriptional regulatory networks

Plant genomes are the source of large numbers of small RNAs, generated via a variety of genetically separable pathways. Several of these pathways converge in the production of phased, secondary, small interfering RNAs (phasiRNAs), originally designated as trans-acting small interfering RNAs or tasiRNAs. PhasiRNA biogenesis requires the involvement of microRNAs as well as the cellular machinery for the production of siRNAs. PhasiRNAs in Arabidopsis thaliana have been well described for their ability to function in trans to suppress target transcript levels. Plant genomic data from an expanding set of species have demonstrated that Arabidopsis is relatively sparing in its use of phasiRNAs, while other genomes contain hundreds or even thousands of phasiRNA-generating loci. In the dicots, targets of those phasiRNAs include several large or conserved families of genes, such as those encoding disease resistance proteins or transcription factors. Suppression of nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes by small RNAs is particularly unusual because of a high level of redundancy. In this review, we discuss plant phasiRNAs and the possible mechanistic significance of phasiRNA-based regulation of the NB-LRRs.

Integrative analysis of miRNA and mRNA profiles in response to ethylene in rose petals during flower opening

MicroRNAs play an important role in plant development and plant responses to various biotic and abiotic stimuli. As one of the most important ornamental crops, rose (Rosa hybrida) possesses several specific morphological and physiological features, including recurrent flowering, highly divergent flower shapes, colors and volatiles. Ethylene plays an important role in regulating petal cell expansion during rose flower opening. Here, we report the population and expression profiles of miRNAs in rose petals during flower opening and in response to ethylene based on high throughput sequencing. We identified a total of 33 conserved miRNAs, as well as 47 putative novel miRNAs were identified from rose petals. The conserved and novel targets to those miRNAs were predicted using the rose floral transcriptome database. Expression profiling revealed that expression of 28 known (84.8% of known miRNAs) and 39 novel (83.0% of novel miRNAs) miRNAs was substantially changed in rose petals during the earlier opening period. We also found that 28 known and 22 novel miRNAs showed expression changes in response to ethylene treatment. Furthermore, we performed integrative analysis of expression profiles of miRNAs and their targets. We found that ethylene-caused expression changes of five miRNAs (miR156, miR164, miR166, miR5139 and rhy-miRC1) were inversely correlated to those of their seven target genes. These results indicate that these miRNA/target modules might be regulated by ethylene and were involved in ethylene-regulated petal growth.

MicroRNA superfamilies descended from miR390 and their roles in secondary small interfering RNA Biogenesis in Eudicots

Trans-acting small interfering RNAs (tasiRNAs) are a major class of small RNAs performing essential biological functions in plants. The first reported tasiRNA pathway, that of miR173-TAS1/2, produces tasiRNAs regulating a set of pentatricopeptide repeat (PPR) genes and has been characterized only in Arabidopsis thaliana to date. Here, we demonstrate that the microRNA (miRNA)-trans-acting small interfering RNA gene (TAS)-pentatricopeptide repeat-containing gene (PPR)-small interfering RNA pathway is a highly dynamic and widespread feature of eudicots. Nine eudicot plants, representing six different plant families, have evolved similar tasiRNA pathways to initiate phased small interfering RNA (phasiRNA) production from PPR genes. The PPR phasiRNA production is triggered by different 22-nucleotide miRNAs, including miR7122, miR1509, and fve-PPRtri1/2, and through distinct mechanistic strategies exploiting miRNA direct targeting or indirect targeting through TAS-like genes (TASL), one-hit or two-hit, or even two layers of tasiRNA-TASL interactions. Intriguingly, although those miRNA triggers display high sequence divergence caused by the occurrence of frequent point mutations and splicing shifts, their corresponding MIRNA genes show pronounced identity to the Arabidopsis MIR173, implying a common origin of this group of miRNAs (super-miR7122). Further analyses reveal that super-miR7122 may have evolved from a newly defined miR4376 superfamily, which probably originated from the widely conserved miR390. The elucidation of this evolutionary path expands our understanding of the course of miRNA evolution, especially for relatively conserved miRNA families.