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

What do microRNA concentrations tell us about the mechanical damage and storage period of strawberry fruits?

Although much research has been performed to investigate the effects of storage conditions, such as mechanical damage and storage period, on the morphological and physiological properties of strawberry fruits, almost all of them have considered severe stress conditions. Finding fruit characteristics that exert significant changes even toward mild and moderate stress conditions can help provide valuable information about the fruit quality during storage. This study aims to investigate various characteristics of strawberry fruits during storage to determine which type of fruit characteristics exert such significant changes toward stress conditions. Identical strawberry samples were subjected to mechanical loading at three levels (1, 2, and 3 N) and then stored at 6 °C for 13 days. Morphological and physiological features, as well as the concentration of several microRNAs involved in strawberry storage, were measured at three-day intervals. The effects of mechanical loading on morphological and physiological characteristics were not significant, while their effects were significant on miR-164, miR-167, and miR-399a. Moreover, while low correlation coefficients were observed between the fruit morphophysiological traits (< 0.6) toward storage conditions, high correlations were obtained between the concentrations of microRNAs. Instead of measuring the morphological and physiological characteristics of fruits, whose behavior is not generally specific toward the stresses, the results show that microRNA concentrations, which can be measured by an electrochemical biosensor, provide us with noteworthy information about fruit quality during storage. These small non-coding molecules exhibited remarkable responses even in mild and moderate stress conditions, making them reliable markers of fruit quality assessment.

A 24-nt miR9560 modulates the transporter gene BrpHMA2 expression in Brassica parachinensis

MicroRNAs (miRNAs) control gene expression in plant through transcript cleavage and translation inhibition. Recently, 24-nt miRNAs have been shown to direct DNA methylation at target sites, regulating the neighboring gene expression. Our study focused on miR9560, a 24-nt miRNA induced by cadmium (Cd) stress in Brassica rapa ssp. parachinensis (B. parachinensis). Phylogenetic analysis revealed miR9560 predominantly emerged in the Rosanae superorder and was conserved in Brassicaceae, with potential target sites adjacent to transporter family genes HMAs. RNA gel blotting showed that mature miR9560 was only detected in various Brassica crops roots after Cd stress. In B. parachinensis, miR9560's putative target site is upstream of BrpHMA2, an afflux-type Cd transporter. In a transient expression system of B. parachinensis protoplasts, the expression of miR9560 increased the DNA methylation upstream of BrpHMA2, reducing the transcription of BrpHMA2. This regulation was also observed in Arabidopsis wild-type protoplasts but not in the mutants dcl234 and ago4 with impairments in the RNA-dependent DNA methylation (RdDM) pathway. We deduced that miR9560 modulates BrpHMA2 expression via the RdDM pathway, potentially regulating Cd uptake and movement in B. parachinensis. Furthermore, this regulatory mechanism may extend to other Brassica plants. This study enhances our comprehension of 24-nt miRNAs role in regulating Cd accumulation within Brassica plants.

A smart multiplexed microRNA biosensor based on FRET for the prediction of mechanical damage and storage period of strawberry fruits

Today, measuring the concentration of various microRNAs in fruits has been introduced to model the storage conditions of agricultural products. However, there is a limiting factor in the extensive utilization of such techniques: the existing methods for measuring microRNA sequences, including PCR and microarrays, are time-consuming and expensive and do not allow for simultaneous measurement of several microRNAs. In this study, a biosensor based on the Förster resonance energy transfer (FRET) of fluorescence dyes that can lead to the hybridization of oligonucleotide probes labeled with such dyes by using an excitation wavelength has been used to simultaneously measure microRNAs. Three microRNA compounds, i.e., miRNA-164, miRNA-167, and miRNA-399a, which play significant roles in the postharvest characteristics of strawberry fruits were measured. The simultaneous measurement was performed using three fluorescence dyes which exert various emission wavelengths at 570, 596, and 670 nm. In the following, machine learning methods including artificial neural networks (ANNs) and support vector machines (SVMs), with hyperparameter values optimized ​with the help of metaheuristic optimization algorithms, were used to predict the amount of mechanical loading on strawberry fruits and their storage period having the microRNA concentrations. The results showed that the SVM with Gaussian kernel, which was optimized by the Harris hawks optimization, is capable of predicting the mechanical stress and storage period of strawberry fruits with a coefficient of determination (R2) of 0.89 and 0.92, respectively. The findings of this study reveal the application of combining FRET-based biosensors and machine learning methods in fruit storage quality assessment.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

The miRNA-mRNA Regulatory Modules of Pinus massoniana Lamb. in Response to Drought Stress

Masson pine (Pinus massoniana Lamb.) is a major fast-growing woody tree species and pioneer species for afforestation in barren sites in southern China. However, the regulatory mechanism of gene expression in P. massoniana under drought remains unclear. To uncover candidate microRNAs, their expression profiles, and microRNA-mRNA interactions, small RNA-seq was used to investigate the transcriptome from seedling roots under drought and rewatering in P. massoniana. A total of 421 plant microRNAs were identified. Pairwise differential expression analysis between treatment and control groups unveiled 134, 156, and 96 differential expressed microRNAs at three stages. These constitute 248 unique microRNAs, which were subsequently categorized into six clusters based on their expression profiles. Degradome sequencing revealed that these 248 differentially expressed microRNAs targeted 2069 genes. Gene Ontology enrichment analysis suggested that these target genes were related to translational and posttranslational regulation, cell wall modification, and reactive oxygen species scavenging. miRNAs such as miR482, miR398, miR11571, miR396, miR166, miRN88, and miRN74, along with their target genes annotated as F-box/kelch-repeat protein, 60S ribosomal protein, copper-zinc superoxide dismutase, luminal-binding protein, S-adenosylmethionine synthase, and Early Responsive to Dehydration Stress may play critical roles in drought response. This study provides insights into microRNA responsive to drought and rewatering in Masson pine and advances the understanding of drought tolerance mechanisms in Pinus.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Indian Himalayan natural Arabidopsis thaliana accessions with abolished miR158 levels exhibit robust miR173-initiated trans-acting cascade silencing

Small RNAs (sRNAs) such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) are short 20-24-nucleotide non-coding RNAs. They are key regulators of gene expression in plants and other organisms. Several 22-nucleotide miRNAs trigger biogenesis cascades of trans-acting secondary siRNAs, which are involved in various developmental and stress responses. Here we show that Himalayan Arabidopsis thaliana accessions having natural mutations in the miR158 locus exhibit robust cascade silencing of the pentatricopeptide repeat (PPR)-like locus. Furthermore, we show that these cascade sRNAs trigger tertiary silencing of a gene involved in transpiration and stomatal opening. The natural deletions or insertions in MIR158 led to improper processing of miR158 precursors, thereby blocking synthesis of mature miR158. Reduced miR158 levels led to increased levels of its target, a pseudo-PPR gene that is targeted by tasiRNAs generated by the miR173 cascade in other accessions. Using sRNA datasets derived from Indian Himalayan accessions, as well as overexpression and knockout lines of miR158, we show that absence of miR158 led to buildup of pseudo-PPR-derived tertiary sRNAs. These tertiary sRNAs mediated robust silencing of a gene involved in stomatal closure in Himalayan accessions lacking miR158 expression. We functionally validated the tertiary phasiRNA that targets NHX2, which encodes a Na⁺ -K⁺ /H⁺ antiporter protein, thereby regulating transpiration and stomatal conductance. Overall, we report the role of the miRNA-TAS-siRNA-pseudogene-tertiary phasiRNA-NHX2 pathway in plant adaptation.

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

BACKGROUND

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

RESULTS

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

CONCLUSIONS

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

Cultivar-specific miRNA-mediated RNA silencing in grapes

In-depth comparative degradome analysis of two domesticated grape cultivars with diverse secondary metabolite accumulation reveals differential miRNA-mediated targeting. Small (s)RNAs such as micro(mi)RNAs and secondary small interfering (si) often work as negative switches of gene expression. In plants, it is well known that miRNAs target and cleave mRNAs that have high sequence complementarity. However, it is not known if there are variations in miRNA-mediated targeting between subspecies and cultivars that have been subjected to vast genetic modifications through breeding and other selections. Here, we have used PAREsnip2 tool for analysis of degradome datasets derived from two contrasting domesticated grape cultivars having varied fruit color, habit and leaf shape. We identified several interesting variations in sRNA targeting using degradome and 5'RACE analysis between two contrasting grape cultivars that was further correlated using RNA-seq analysis. Several of the differences we identified are associated with secondary metabolic pathways. We propose possible means by which sRNAs might contribute to diversity in secondary metabolites and other development pathways between two domesticated cultivars of grapes.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development

MicroRNAs (miRNAs) are a class of 21-24 nucleotides (nt) noncoding small RNAs ubiquitously distributed across the plant kingdom. miR482/2118, one of the conserved miRNA superfamilies originating from gymnosperms, has divergent main functions in core-angiosperms. It mainly regulates NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT (NBS-LRR) genes in eudicots, functioning as an essential component in plant disease resistance; in contrast, it predominantly targets numerous long noncoding RNAs (lncRNAs) in monocot grasses, which are vital for plant reproduction. Usually, miR482/2118 is 22-nt in length, which can trigger the production of phased small interfering RNAs (phasiRNAs) after directed cleavage. PhasiRNAs instigated from target genes of miR482/2118 enhance their roles in corresponding biological processes by cis-regulation on cognate genes and expands their function to other pathways via trans activity on different genes. This review summarizes the origin, biogenesis, conservation, and evolutionary characteristics of the miR482/2118 superfamily and delineates its diverse functions in disease resistance, plant development, stress responses, etc.

PhasiRNAnalyzer: an integrated analyser for plant phased siRNAs

Phased siRNAs (phasiRNAs) are a class of small interfering RNAs (siRNAs) which play essential roles in plant development and defence. However, only a few phasiRNAs have been extensively studied due to the difficulties in identifying and characterizing plant phasiRNAs by plant biologists. Herein, we describe a comprehensive and multi-functional web server termed PhasiRNAnalyzer, which is able to identify all crucial components in plant phasiRNA's regulatory pathway (phase-initiator→PHAS gene→phasiRNA cluster→target gene). Currently, PhasiRNAnalyzer exhibits the following advantages: I) It is the most comprehensive platform which hosts 170 plant species with 256 genome data, 438 cDNA data and 271 degradome data. II) It can identify all crucial components in phasiRNA's regulatory pathway, and verify the interactions between phasiRNAs and their target genes based on degradome data. III) It can perform differential expression analysis of phasiRNAs on each PHAS gene locus between different samples conveniently. IV) It provides the user-friendly interfaces and introduces several improvements, primarily by making more accurate and efficient analysis when dealing with deep sequencing data. In summary, PhasiRNAnalyzer is a comprehensive and systemic phasiRNA analysis server with high sensitivity and efficiency. It can be freely accessed at https://cbi.njau.edu.cn/PPSA/.

Conserved and non-conserved triggers of 24-nucleotide reproductive phasiRNAs in eudicots

Small RNAs play important roles in plant growth and development by modulating expression of genes and transposons. In many flowering plant species, male reproductive organs, the anthers, produce abundant phased small interfering RNAs (phasiRNAs). Two classes of reproductive phasiRNAs are generally known, mostly from monocots: (i) pre-meiotic 21-nucleotide (nt) phasiRNAs triggered by miR2118 and (ii) meiotic 24-nt phasiRNAs triggered by miR2275. Here, we describe conserved and non-conserved triggers of 24-nt phasiRNAs in several eudicots. We found that the abundant 24-nt phasiRNAs in the basal eudicot columbine (Aquilegia coerulea) are produced by the canonical trigger miR2275, as well as by other non-canonical triggers, miR482/2118 and miR14051. These triggering microRNAs (miRNAs) are localized in microspore mother cells and tapetal cells of meiotic and post-meiotic stage anthers. Furthermore, we identified a lineage-specific trigger (miR11308) of 24-nt phasiRNAs and an expanded number of 24-PHAS loci in wild strawberry (Fragaria vesca). We validated the presence of the miR2275-derived 24-nt phasiRNA pathway in rose (Rosa chinensis). Finally, we evaluated all eudicots that have been validated for the presence of 24-nt phasiRNAs as possible model systems in which to study the biogenesis and function of 24-nt phasiRNAs. We conclude that columbine (Aquilegia coerulea) would be a strong model because of its extensive number of 24-PHAS loci and its diversity of trigger miRNAs.

Pre-meiotic 21-nucleotide reproductive phasiRNAs emerged in seed plants and diversified in flowering plants

Plant small RNAs are important regulatory elements that fine-tune gene expression and maintain genome integrity by silencing transposons. Reproductive organs of monocots produce abundant phased, small interfering RNAs (phasiRNAs). The 21-nt reproductive phasiRNAs triggered by miR2118 are highly enriched in pre-meiotic anthers, and have been found in multiple eudicot species, in contrast with prior reports of monocot specificity. The 24-nt reproductive phasiRNAs are triggered by miR2275, and are highly enriched during meiosis in many angiosperms. Here, we report the widespread presence of the 21-nt reproductive phasiRNA pathway in eudicots including canonical and non-canonical microRNA (miRNA) triggers of this pathway. In eudicots, these 21-nt phasiRNAs are enriched in pre-meiotic stages, a spatiotemporal distribution consistent with that of monocots and suggesting a role in anther development. Although this pathway is apparently absent in well-studied eudicot families including the Brassicaceae, Solanaceae and Fabaceae, our work in eudicots supports an earlier singular finding in spruce, a gymnosperm, indicating that the pathway of 21-nt reproductive phasiRNAs emerged in seed plants and was lost in some lineages.

An Atlas of Genomic Resources for Studying Rosaceae Fruits and Ornamentals

Rosaceae, a large plant family of more than 3,000 species, consists of many economically important fruit and ornamental crops, including peach, apple, strawberry, raspberry, cherry, and rose. These horticultural crops are not only important economic drivers in many regions of the world, but also major sources of human nutrition. Additionally, due to the diversity of fruit types in Rosaceae, this plant family offers excellent opportunities for investigations into fleshy fruit diversity, evolution, and development. With the development of high-throughput sequencing technologies and computational tools, an increasing number of high-quality genomes and transcriptomes of Rosaceae species have become available and will greatly facilitate Rosaceae research and breeding. This review summarizes major genomic resources and genome research progress in Rosaceae, highlights important databases, and suggests areas for further improvement. The availability of these big data resources will greatly accelerate research progress and enhance the agricultural productivity of Rosaceae.

The long noncoding RNA FRILAIR regulates strawberry fruit ripening by functioning as a noncanonical target mimic

Long noncoding RNAs (lncRNAs) are emerging as important regulators in plant development, but few of them have been functionally characterized in fruit ripening. Here, we have identified 25,613 lncRNAs from strawberry ripening fruits based on RNA-seq data from poly(A)-depleted libraries and rRNA-depleted libraries, most of which exhibited distinct temporal expression patterns. A novel lncRNA, FRILAIR harbours the miR397 binding site that is highly conserved in diverse strawberry species. FRILAIR overexpression promoted fruit maturation in the Falandi strawberry, which was consistent with the finding from knocking down miR397, which can guide the mRNA cleavage of both FRILAIR and LAC11a (encoding a putative laccase-11-like protein). Moreover, LAC11a mRNA levels were increased in both FRILAIR overexpressing and miR397 knockdown fruits, and accelerated fruit maturation was also found in LAC11a overexpressing fruits. Overall, our study demonstrates that FRILAIR can act as a noncanonical target mimic of miR397 to modulate the expression of LAC11a in the strawberry fruit ripening process.

Evolution and diversification of reproductive phased small interfering RNAs in Oryza species

In grasses, two types of phased, small interfering RNAs (phasiRNAs) are expressed largely in young, developing anthers. They are 21 or 24 nucleotides (nt) in length and are triggered by miR2118 or miR2275, respectively. However, most of their functions and activities are not fully understood. We performed comparative genomic analysis of their source loci (PHAS) in five Oryza genomes and combined this with analysis of high-throughput sRNA and degradome datasets. In total, we identified 8216 21-PHAS and 626 24-PHAS loci. Local tandem and segmental duplications mainly contributed to the expansion and supercluster distribution of the 21-PHAS loci. Despite their relatively conserved genomic positions, PHAS sequences diverged rapidly, except for the miR2118/2275 target sites, which were under strong selection for conservation. We found that 21-nt phasiRNAs with a 5'-terminal uridine (U) demonstrated cis-cleavage at PHAS precursors, and these cis-acting sites were also variable among close species. miR2118 could trigger phasiRNA production from its own antisense transcript and the derived phasiRNAs might reversibly regulate miR2118 precursors. We hypothesised that successful initiation of phasiRNA biogenesis is conservatively maintained, while phasiRNA products diverged quickly and are not individually conserved. In particular, phasiRNA production is under the control of multiple reciprocal regulation mechanisms.

Secondary siRNAs in Plants: Biosynthesis, Various Functions, and Applications in Virology

The major components of RNA silencing include both transitive and systemic small RNAs, which are technically called secondary sRNAs. Double-stranded RNAs trigger systemic silencing pathways to negatively regulate gene expression. The secondary siRNAs generated as a result of transitive silencing also play a substantial role in gene silencing especially in antiviral defense. In this review, we first describe the discovery and pathways of transitivity with emphasis on RNA-dependent RNA polymerases followed by description on the short range and systemic spread of silencing. We also provide an in-depth view on the various size classes of secondary siRNAs and their different roles in RNA silencing including their categorization based on their biogenesis. The other regulatory roles of secondary siRNAs in transgene silencing, virus-induced gene silencing, transitivity, and trans-species transfer have also been detailed. The possible implications and applications of systemic silencing and the different gene silencing tools developed are also described. The details on mobility and roles of secondary siRNAs derived from viral genome in plant defense against the respective viruses are presented. This entails the description of other compatible plant-virus interactions and the corresponding small RNAs that determine recovery from disease symptoms, exclusion of viruses from shoot meristems, and natural resistance. The last section presents an overview on the usefulness of RNA silencing for management of viral infections in crop plants.

PhasiRNAs in Plants: Their Biogenesis, Genic Sources, and Roles in Stress Responses, Development, and Reproduction

Phased secondary small interfering RNAs (phasiRNAs) constitute a major category of small RNAs in plants, but most of their functions are still poorly defined. Some phasiRNAs, known as trans-acting siRNAs, are known to target complementary mRNAs for degradation and to function in development. However, the targets or biological roles of other phasiRNAs remain speculative. New insights into phasiRNA biogenesis, their conservation, and their variation across the flowering plants continue to emerge due to the increased availability of plant genomic sequences, deeper and more sophisticated sequencing approaches, and improvements in computational biology and biochemical/molecular/genetic analyses. In this review, we survey recent progress in phasiRNA biology, with a particular focus on two classes associated with male reproduction: 21-nucleotide (accumulate early in anther ontogeny) and 24-nucloetide (produced in somatic cells during meiosis) phasiRNAs. We describe phasiRNA biogenesis, function, and evolution and define the unanswered questions that represent topics for future research.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Identification and Characterization of Stress-Responsive TAS3-Derived TasiRNAs in Melon

Small interfering RNAs (siRNA) are key regulators of gene expression that play essential roles in diverse biological processes. Trans-acting siRNAs (tasiRNAs) are a class of plant-endogenous siRNAs that lead the cleavage of nonidentical transcripts. TasiRNAs are usually involved in fine-tuning development. However, increasing evidence supports that tasiRNAs may be involved in stress response. Melon is a crop of great economic importance extensively cultivated in semiarid regions frequently exposed to changing environmental conditions that limit its productivity. However, knowledge of the precise role of siRNAs in general, and of tasiRNAs in particular, in regulating the response to adverse environmental conditions is limited. Here, we provide the first comprehensive analysis of computationally inferred melon-tasiRNAs responsive to two biotic (viroid-infection) and abiotic (cold treatment) stress conditions. We identify two TAS3-loci encoding to length (TAS3-L) and short (TAS3-S) transcripts. The TAS candidates predicted from small RNA-sequencing data were characterized according to their chromosome localization and expression pattern in response to stress. The functional activity of cmTAS genes was validated by transcript quantification and degradome assays of the tasiRNA precursors and their predicted targets. Finally, the functionality of a representative cmTAS3-derived tasiRNA (TAS3-S) was confirmed by transient assays showing the cleavage of ARF target transcripts.

Conserved and novel roles of miR164-CUC2 regulatory module in specifying leaf and floral organ morphology in strawberry

MicroRNAs (miRNAs) are a kind of short noncoding RNA (20-24 nt), playing versatile roles in plant growth and development. Strawberry generates leaves and flowers with unique features. However, few miRNAs have been functionally characterised in strawberry, especially for their developmental regulation. Here, we identified one ethyl methanesulfonate (EMS) mutant, deeply serrated (des), in the woodland strawberry Fragaria vesca that has wrinkled leaves with deeper serrations, serrated petals and deformed carpels. The causative mutation occurs in the 19^th nucleotide of the FvemiR164a mature sequence. Overexpressing FveMIR164A rescued the phenotypes of des/fvemir164a except the petal serrations. Furthermore, we identified two allelic mutants of FveCUC2a, one target of FvemiR164a, which developed leaves with smooth margins and fused leaflets. Phenotypes of the double mutant fvemir164a fvecuc2a indicated that the two genes act linearly in leaf and carpel development, but synergistically in the development of other floral organs and inflorescence architecture. This work demonstrates the conserved and novel roles of the miR164-CUC2 module in leaf and flower development in different plant species, and reveals that the 19^th nucleotide of FvemiR164a is important for its processing.

miR828 and miR858 regulate VvMYB114 to promote anthocyanin and flavonol accumulation in grapes

MicroRNAs are a class of non-coding small RNAs involved in the negative regulation of gene expression, which play critical roles in developmental and metabolic pathways. Studies in several plants have identified a few microRNAs and other small RNAs that target regulators of the phenylpropanoid metabolic pathway called the MYB transcription factors. However, it is not well understood how sRNA-mediated regulation of MYBs influences the accumulation of specific secondary metabolites. Using sRNA sequencing, degradome analysis, mRNA sequencing, and proteomic analysis, we establish that grape lines with high anthocyanin content express two MYB-targeting microRNAs abundantly, resulting in the differential expression of specific MYB proteins. miR828 and miR858 target coding sequences of specific helix motifs in the mRNA sequences of MYB proteins. Targeting by miR828 caused MYB RNA decay and the production of a cascade of secondary siRNAs that depend on RNA-dependent RNA polymerase 6. MYB suppression and cascade silencing was more robust in grape lines with high anthocyanin content than in a flavonol-rich grape line. We establish that microRNA-mediated silencing targeted the repressor class of MYBs to promote anthocyanin biosynthesis in grape lines with high anthocyanins. We propose that this process regulates the expression of appropriate MYBs in grape lines to produce specific secondary metabolites.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

24-nt reproductive phasiRNAs are broadly present in angiosperms

Small RNAs are key regulators in plant growth and development. One subclass, phased siRNAs (phasiRNAs) require a trigger microRNA for their biogenesis. In grasses, two pathways yield abundant phasiRNAs during anther development; miR2275 triggers one class, 24-nt phasiRNAs, coincident with meiosis, while a second class of 21-nt phasiRNAs are present in premeiotic anthers. Here we report that the 24-nt phasiRNA pathway is widely present in flowering plants, indicating that 24-nt reproductive phasiRNAs likely originated with the evolutionary emergence of anthers. Deep comparative genomic analyses demonstrated that this miR2275/24-nt phasiRNA pathway is widely present in eudicots plants, however, it is absent in legumes and in the model plant Arabidopsis, demonstrating a dynamic evolutionary history of this pathway. In Solanaceae species, 24-nt phasiRNAs were observed, but the miR2275 trigger is missing and some loci displaying 12-nt phasing. Both the miR2275-triggered and Solanaceae 24-nt phasiRNAs are enriched in meiotic stages, implicating these phasiRNAs in anther and/or pollen development, a spatiotemporal pattern consistent in all angiosperm lineages that deploy them.

Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening

The perception and signal transduction of the plant hormone abscisic acid (ABA) are crucial for strawberry fruit ripening, but the underlying mechanism of how ABA regulates ripening-related genes has not been well understood. By employing high-throughput sequencing technology, we comprehensively analyzed transcriptomic and miRNA expression profiles simultaneously in ABA- and nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker)-treated strawberry fruits with temporal resolution. The results revealed that ABA regulated many genes in different pathways, including hormone signal transduction and the biosynthesis of secondary metabolites. Transcription factor genes belonging to WRKY and heat shock factor (HSF) families might play key roles in regulating the expression of ABA inducible genes, whereas the KNOTTED1-like homeobox protein and Squamosa Promoter-Binding-like protein 18 might be responsible for ABA-downregulated genes. Additionally, 20 known and six novel differentially expressed miRNAs might be important regulators that assist ABA in regulating target genes that are involved in versatile physiological processes, such as hormone balance regulation, pigments formation and cell wall degradation. Furthermore, degradome analysis showed that one novel miRNA, Fa_novel6, could degrade its target gene HERCULES1, which likely contributed to fruit size determination during strawberry ripening. These results expanded our understanding of how ABA drives the strawberry fruit ripening process as well as the role of miRNAs in this process.

Integrated analysis of high-throughput sequencing data shows abscisic acid-responsive genes and miRNAs in strawberry receptacle fruit ripening

The perception and signal transduction of the plant hormone abscisic acid (ABA) are crucial for strawberry fruit ripening, but the underlying mechanism of how ABA regulates ripening-related genes has not been well understood. By employing high-throughput sequencing technology, we comprehensively analyzed transcriptomic and miRNA expression profiles simultaneously in ABA- and nordihydroguaiaretic acid (NDGA, an ABA biosynthesis blocker)-treated strawberry fruits with temporal resolution. The results revealed that ABA regulated many genes in different pathways, including hormone signal transduction and the biosynthesis of secondary metabolites. Transcription factor genes belonging to WRKY and heat shock factor (HSF) families might play key roles in regulating the expression of ABA inducible genes, whereas the KNOTTED1-like homeobox protein and Squamosa Promoter-Binding-like protein 18 might be responsible for ABA-downregulated genes. Additionally, 20 known and six novel differentially expressed miRNAs might be important regulators that assist ABA in regulating target genes that are involved in versatile physiological processes, such as hormone balance regulation, pigments formation and cell wall degradation. Furthermore, degradome analysis showed that one novel miRNA, Fa_novel6, could degrade its target gene HERCULES1, which likely contributed to fruit size determination during strawberry ripening. These results expanded our understanding of how ABA drives the strawberry fruit ripening process as well as the role of miRNAs in this process.

Updated annotation of the wild strawberry Fragaria vesca V4 genome

The diploid strawberry Fragaria vesca serves as an ideal model plant for cultivated strawberry (Fragaria × ananassa, 8x) and the Rosaceae family. The F. vesca genome was initially published in 2011 using older technologies. Recently, a new and greatly improved F. vesca genome, designated V4, was published. However, the number of annotated genes is remarkably reduced in V4 (28,588 genes) compared to the prior annotations (32,831 to 33,673 genes). Additionally, the annotation of V4 (v4.0.a1) implements a new nomenclature for gene IDs (FvH4_XgXXXXX), rather than the previous nomenclature (geneXXXXX). Hence, further improvement of the V4 genome annotation and assigning gene expression levels under the new gene IDs with existing transcriptome data are necessary to facilitate the utility of this high-quality F. vesca genome V4. Here, we built a new and improved annotation, v4.0.a2, for F. vesca genome V4. The new annotation has a total of 34,007 gene models with 98.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCOs). In this v4.0.a2 annotation, gene models of 8,342 existing genes are modified, 9,029 new genes are added, and 10,176 genes possess alternatively spliced isoforms with an average of 1.90 transcripts per locus. Transcription factors/regulators and protein kinases are globally identified. Interestingly, the transcription factor family FAr-red-impaired Response 1 (FAR1) contains 82 genes in v4.0.a2 but only two members in v4.0.a1. Additionally, the expression levels of all genes in the new annotation across a total of 46 different tissues and stages are provided. Finally, miRNAs and their targets are reanalyzed and presented. Altogether, this work provides an updated genome annotation of the F. vesca V4 genome as well as a comprehensive gene expression atlas with the new gene ID nomenclature, which will greatly facilitate gene functional studies in strawberry and other evolutionarily related plant species.

Genome-wide identification and expression analysis of the SPL gene family in woodland strawberry Fragaria vesca

SQUAMOSA promoter-binding protein-like (SPL) is a class of plant-specific transcription factors that play critical roles in regulating plant growth and development. However, little systematic research on SPL genes has been conducted in strawberry. In this study, 14 SPL genes were identified in the genome of woodland strawberry (Fragaria vesca), one of the model plants of the family Rosaceae. Chromosome localization analysis indicated that the 14 FvSPL genes were unevenly distributed on six chromosomes. Phylogenetic analysis indicated that the FvSPL proteins could be clustered into six groups (G1 to G6). Genes with similar structure were classified into the same group, implying their functional redundancy. In addition, nine out of the 14 FvSPL genes, belonging to G1, G2, and G5, were found to be the putative targets of FvmiR156 genes. Expression analysis indicated FvSPL genes exhibited highly diverse expression patterns in the tissues and organs examined. The transcript levels of most FvmiR156-targeted FvSPL genes in fruit were lower than those non-miR156-targeted genes. In addition, the expression of the FvmiR156-targeted FvSPL genes decreased during fruit ripening, whereas the expression of FvmiR156 genes increased in fruit during this process. The results provide a foundation for future functional analysis of FvSPL genes in strawberry growth and development.

Biogenesis and regulatory hierarchy of phased small interfering RNAs in plants

Several varieties of small RNAs including microRNAs (miRNAs) and small interfering RNAs (siRNAs) are generated in plants to regulate development, genome stability and response to adverse environments. Phased siRNA (phasiRNA) is a type of secondary siRNA that is processed from a miRNA-mediated cleavage of RNA transcripts, increasing silencing efficiency or simultaneously suppressing multiple target genes. Trans-acting siRNAs (ta-siRNAs) are a particular class of phasiRNA produced from noncoding transcripts that silence targets in trans. It was originally thought that 'one-hit' and 'two-hit' models were essential for processing distinct TAS precursors; however, a single hit event was recently shown to be sufficient at triggering all types of ta-siRNAs. This review discusses the findings about biogenesis, targeting modes and regulatory networks of plant ta-siRNAs. We also summarize recent advances in the generation of other phasiRNAs and their possible biological benefits to plants.

Despacito: the slow evolutionary changes in plant microRNAs

MicroRNAs (miRNAs) are key regulators of gene expression. A handful of miRNAs are broadly conserved in land plants, while the majority are lineage specific; this review describes the processes by which new miRNAs are hypothesized to have emerged. Two major models describe miRNA origins, firstly, de novo emergence via inverted duplication of target gene fragments, and secondly, the expansion and neofunctionalization of existing miRNA families. The occasional acquisition of target sites by previously un-targeted genes adds further dynamism to the process by which miRNAs may shift roles during evolution. Additional factors guiding miRNA evolution include functional constraints on their length and the importance of precursor conservation that is observed in regions above or below the mature miRNA duplex; these regions represent recognition sites for components of biogenesis machinery and direct precursor processing. Insights into the mechanisms of miRNA emergence and divergence are important for understanding plant genome evolution and the impact of miRNA regulatory networks.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

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

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

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

FRG3, a Target of slmiR482e-3p, Provides Resistance against the Fungal Pathogen Fusarium oxysporum in Tomato

The vast majority of plant disease resistance (R) genes encode nucleotide binding site-leucine-rich repeat (NBS-LRR) proteins, which specifically determine the plant immune response and have been demonstrated to be targets of several microRNA (miRNA) families. The fungus Fusarium oxysporum f. sp. lycopersici (FOL) causes vascular wilt disease in tomato worldwide. Here, we explored a possible role for FGR3 in tomato defense against FOL. FRG3 is a predicted NBS-LRR like gene that is targeted by slmiR482e-3p, a member of slmiR482 miRNA family. Northern blot data demonstrated that all seven members of the slmiR482 family were regulated in diverse ways after infection by FOL. The ability of FRG3 to be regulated by slmiR482e-3p was confirmed at the transcript level by co-expression studies in Nicotiana benthamiana. A virus-induced gene silencing (VIGS) approach revealed that FRG3 confers resistance to the Motelle tomato cultivar. Taken together, our study has identified a novel R gene, FRG3, which is targeted by slmiR482e-3p at the transcript level, and is necessary for resistance to tomato wilt disease in planta.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Genome re-annotation of the wild strawberry Fragaria vesca using extensive Illumina- and SMRT-based RNA-seq datasets

The genome of the wild diploid strawberry species Fragaria vesca, an ideal model system of cultivated strawberry (Fragaria × ananassa, octoploid) and other Rosaceae family crops, was first published in 2011 and followed by a new assembly (Fvb). However, the annotation for Fvb mainly relied on ab initio predictions and included only predicted coding sequences, therefore an improved annotation is highly desirable. Here, a new annotation version named v2.0.a2 was created for the Fvb genome by a pipeline utilizing one PacBio library, 90 Illumina RNA-seq libraries, and 9 small RNA-seq libraries. Altogether, 18,641 genes (55.6% out of 33,538 genes) were augmented with information on the 5′ and/or 3′ UTRs, 13,168 (39.3%) protein-coding genes were modified or newly identified, and 7,370 genes were found to possess alternative isoforms. In addition, 1,938 long non-coding RNAs, 171 miRNAs, and 51,714 small RNA clusters were integrated into the annotation. This new annotation of F. vesca is substantially improved in both accuracy and integrity of gene predictions, beneficial to the gene functional studies in strawberry and to the comparative genomic analysis of other horticultural crops in Rosaceae family.

The Emergence, Evolution, and Diversification of the miR390-TAS3-ARF Pathway in Land Plants

In plants, miR390 directs the production of tasiRNAs from TRANS-ACTING SIRNA3 (TAS3) transcripts to regulate AUXIN RESPONSIVE FACTOR (ARF) genes, critical for auxin signaling; these tasiRNAs are known as tasiARFs. To understand the evolution of this miR390-TAS3-ARF pathway, we characterized homologs of these three genes from thousands of plant species, from bryophytes to angiosperms. We found the lower-stem region of MIR390 genes, critical for accurate DICER-LIKE1 processing, is conserved in sequence in seed plants. We propose a model for the transition of functional tasiRNA sequences in TAS3 genes occurred at the emergence of vascular plants, in which the two miR390 target sites of TAS3 genes showed distinct pairing patterns. Based on the cleavability of miR390 target sites and the distance between target site and tasiARF, we inferred a potential bidirectional processing mechanism exists for some TAS3 genes. We also demonstrated a tight mutual selection between tasiARF and its target genes and that ARGONAUTE7, the partner of miR390, was specified later than other factors in the pathway. All these data illuminate the evolutionary path of the miR390-TAS3-ARF pathway in land plants and demonstrate the significant variation that occurs in this functionally important and archetypal regulatory circuit.

The evolution of microRNAs in plants

MicroRNAs (miRNAs) are a central player in post-transcriptional regulation of gene expression and are involved in numerous biological processes in eukaryotes. Knowledge of the origins and divergence of miRNAs paves the way for a better understanding of the complexity of the regulatory networks that they participate in. The biogenesis, degradation, and regulatory activities of miRNAs are relatively better understood, but the evolutionary history of miRNAs still needs more exploration. Inverted duplication of target genes, random hairpin sequences and small transposable elements constitute three main models that explain the origination of miRNA genes (MIR). Both inter-species and intra-species divergence of miRNAs exhibits functional adaptation and adaptation to changing environments in evolution. Here we summarize recent progress in studies on the evolution of MIR and related genes.

PMS1T, producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice

Phased small-interfering RNAs (phasiRNAs) are a special class of small RNAs, which are generated in 21- or 24-nt intervals from transcripts of precursor RNAs. Although phasiRNAs have been found in a range of organisms, their biological functions in plants have yet to be uncovered. Here we show that phasiRNAs generated by the photopheriod-sensetive genic male sterility 1 (Pms1) locus were associated with photoperiod-sensitive male sterility (PSMS) in rice, a germplasm that started the two-line hybrid rice breeding. The Pms1 locus encodes a long-noncoding RNA PMS1T that was preferentially expressed in young panicles. PMS1T was targeted by miR2118 to produce 21-nt phasiRNAs that preferentially accumulated in the PSMS line under long-day conditions. A single nucleotide polymorphism in PMS1T nearby the miR2118 recognition site was critical for fertility change, likely leading to differential accumulation of the phasiRNAs. This result suggested possible roles of phasiRNAs in reproductive development of rice, demonstrating the potential importance of this RNA class as regulators in biological processes.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

Comparative Analysis and Identification of miRNAs and Their Target Genes Responsive to Salt Stress in Diploid and Tetraploid Paulownia fortunei Seedlings

Salt stress is a global environmental problem that affects plant growth and development. Paulownia fortunei is an adaptable and fast-growing deciduous tree native to China that is environmentally and economically important. MicroRNAs (miRNAs) play important regulatory roles in growth, development, and stress responses in plants. MiRNAs that respond to biotic stresses have been identified; however, how miRNAs in P. fortunei respond to salt stress has not yet been reported. To identify salt-stress-responsive miRNAs and predict their target genes, four small RNA and four degradome libraries were constructed from NaCl-treated and NaCl-free leaves of P. fortunei seedlings. The results indicated that salt stress had different physiological effects on diploid and tetraploid P. fortunei. We detected 53 conserved miRNAs belonging to 17 miRNA families and 134 novel miRNAs in P. fortunei. Comparing their expression levels in diploid and tetraploid P. fortunei, we found 10 conserved and 10 novel miRNAs that were significantly differentially expressed under salt treatment, among them eight were identified as miRNAs probably associated with higher salt tolerance in tetraploid P. fortunei than in diploid P. fortunei. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the target genes of the conserved and novel miRNAs. The expressions of 10 differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on P. fortunei miRNAs and their target genes under salt stress. The results provided information at the physiological and molecular levels for further research into the response mechanisms of P. fortunei to salt stress.

FragariaCyc: A Metabolic Pathway Database for Woodland Strawberry Fragaria vesca

FragariaCyc is a strawberry-specific cellular metabolic network based on the annotated genome sequence of Fragaria vesca L. ssp. vesca, accession Hawaii 4. It was built on the Pathway-Tools platform using MetaCyc as the reference. The experimental evidences from published literature were used for supporting/editing existing entities and for the addition of new pathways, enzymes, reactions, compounds, and small molecules in the database. To date, FragariaCyc comprises 66 super-pathways, 488 unique pathways, 2348 metabolic reactions, 3507 enzymes, and 2134 compounds. In addition to searching and browsing FragariaCyc, researchers can compare pathways across various plant metabolic networks and analyze their data using Omics Viewer tool. We view FragariaCyc as a resource for the community of researchers working with strawberry and related fruit crops. It can help understanding the regulation of overall metabolism of strawberry plant during development and in response to diseases and abiotic stresses. FragariaCyc is available online at http://pathways.cgrb.oregonstate.edu.

Comparative Analysis of miRNAs and Their Target Transcripts between a Spontaneous Late-Ripening Sweet Orange Mutant and Its Wild-Type Using Small RNA and Degradome Sequencing

Fruit ripening in citrus is not well-understood at the molecular level. Knowledge of the regulatory mechanism of citrus fruit ripening at the post-transcriptional level in particular is lacking. Here, we comparatively analyzed the miRNAs and their target genes in a spontaneous late-ripening mutant, "Fengwan" sweet orange (MT) (Citrus sinensis L. Osbeck), and its wild-type counterpart ("Fengjie 72-1," WT). Using high-throughput sequencing of small RNAs and RNA degradome tags, we identified 107 known and 21 novel miRNAs, as well as 225 target genes. A total of 24 miRNAs (16 known miRNAs and 8 novel miRNAs) were shown to be differentially expressed between MT and WT. The expression pattern of several key miRNAs and their target genes during citrus fruit development and ripening stages was examined. Csi-miR156k, csi-miR159, and csi-miR166d suppressed specific transcription factors (GAMYBs, SPLs, and ATHBs) that are supposed to be important regulators involved in citrus fruit development and ripening. In the present study, miRNA-mediated silencing of target genes was found under complicated and sensitive regulation in citrus fruit. The identification of miRNAs and their target genes provide new clues for future investigation of mechanisms that regulate citrus fruit ripening.

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.

Post-transcriptional gene silencing in plants: a double-edged sword

In plants, post-transcriptional gene silencing (PTGS) protects the genome from foreign genes and restricts the expression of certain endogenous genes for proper development. Here, we review the recent progress about how the unwanted PTGS is avoided in plants. As a decision-making step of PTGS, aberrant transcripts from most endogenous coding genes are strictly sorted to the bidirectional RNA decay pathways in cytoplasm but not to the short interference RNA (siRNA)-mediated PTGS, with the exception of a few development-relevant endogenous siRNA-producing genes. We also discuss a finely balanced PTGS threshold model that plants fully take advantage of the power of PTGS without self-harm.

Global identification and analysis of long non-coding RNAs in diploid strawberry Fragaria vesca during flower and fruit development

Background

Long non-coding RNAs (lncRNAs) are a new class of regulatory molecules with roles in diverse biological processes. While much effort has been invested in the analysis of lncRNAs from established plant models Arabidopsis, maize, and rice, almost nothing is known about lncRNAs from fruit crops, including those in the Rosaceae family.

Results

Here, we present a genome-scale identification and characterization of lncRNAs from a diploid strawberry, Fragaria vesca, based on rich RNA-seq datasets from 35 different flower and fruit tissues. 5,884 Fve-lncRNAs derived from 3,862 loci were identified. These lncRNAs were carefully cataloged based on expression level and whether or not they contain repetitive sequences or generate small RNAs. About one fourth of them are termed high-confidence lncRNAs (hc-lncRNAs) because they are expressed at a level of FPKM higher than 2 and produce neither small RNAs nor contain repetitive sequence. To identify regulatory interactions between lncRNAs and their potential protein-coding (PC) gene targets, pairs of lncRNAs and PC genes with positively or negatively correlated expression trends were identified based on their expression; these pairs may be candidates of cis- or trans-acting lncRNAs and their targets. Finally, blast searches within plant species indicate that lncRNAs are not well conserved.

Conclusions

Our study identifies a large number of tissue-specifically expressed lncRNAs in F. vesca, thereby highlighting their potential contributions to strawberry flower and fruit development and paving the way for future functional studies.

Electronic supplementary material

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

Multilocus Sex Determination Revealed in Two Populations of Gynodioecious Wild Strawberry, Fragaria vesca subsp. bracteata

Gynodioecy, the coexistence of females and hermaphrodites, occurs in 20% of angiosperm families and often enables transitions between hermaphroditism and dioecy. Clarifying mechanisms of sex determination in gynodioecious species can thus illuminate sexual system evolution. Genetic determination of gynodioecy, however, can be complex and is not fully characterized in any wild species. We used targeted sequence capture to genetically map a novel nuclear contributor to male sterility in a self-pollinated hermaphrodite of Fragaria vesca subsp. bracteata from the southern portion of its range. To understand its interaction with another identified locus and possibly additional loci, we performed crosses within and between two populations separated by 2000 km, phenotyped the progeny and sequenced candidate markers at both sex-determining loci. The newly mapped locus contains a high density of pentatricopeptide repeat genes, a class commonly involved in restoration of fertility caused by cytoplasmic male sterility. Examination of all crosses revealed three unlinked epistatically interacting loci that determine sexual phenotype and vary in frequency between populations. Fragaria vesca subsp. bracteata represents the first wild gynodioecious species with genomic evidence of both cytoplasmic and nuclear genes in sex determination. We propose a model for the interactions between these loci and new hypotheses for the evolution of sex determining chromosomes in the subdioecious and dioecious Fragaria.