Deep sequencing discovery of novel and conserved microRNAs in strawberry (Fragaria×ananassa)

Comparative Transcriptome Analysis Uncovers the Regulatory Roles of MicroRNAs Involved in Petal Color Change of Pink-Flowered Strawberry

The pink-flowered strawberry is popular in China due to its high ornamental value. In the present study, sRNAome, transcriptome, and degradome sequencing were performed to understand the functions of microRNAs (miRNAs) and their target genes during flower development in pink-flowered strawberry. Nine small RNA libraries and a mixed degradome library from flower petals at different developmental stages were constructed and sequenced. A total of 739 known miRNAs and 964 novel miRNAs were identified via small RNA sequencing, and 639 miRNAs were identified to cleave 2,816 target genes based on the degradome data. Additionally, 317 differentially expressed miRNAs among the various stages of flower development were identified, which regulated 2,134 differentially expressed target genes. These target genes were significantly enriched in the transcriptional regulation, phenylpropanoid biosynthesis, and plant hormone signal transduction pathways. Furthermore, integrated microRNAomic and transcriptomic analyses suggested that 98 miRNAs targeted several transcription factors, including MYBs (26), bHLHs (12), NACs (14), and SPLs (19), related to anthocyanin accumulation. In addition, 27 differentially expressed miRNAs might affect anthocyanin biosynthesis by regulating 23 targets involved in the hormone signal transduction pathway. The quantitative real-time PCR (qRT-PCR) analysis confirmed the expression changes of 21 miRNA-target pairs. Furthermore, the transient expression of candidate miRNAs was performed in the pink-flowered strawberry cultivar "Fenyun" at the bud stage. Introduction of FamiR156a, FamiR396e, and FamiR858_R-2 in the "Fenyun" increased flower color intensity, while transient expression of FamiR828a decreased flower color intensity. Overall, the present study uncovers the regulatory functions of microRNAs, including anthocyanin biosynthesis, hormone signaling, and regulation factors during flower development and coloration in pink-flowered strawberry. This work expands the knowledge of miRNAs affecting coloration in strawberry and provides rich resources for future functional studies.

Identification and verification of seed development related miRNAs in kernel almond by small RNA sequencing and qPCR

Many studies have investigated the role of miRNAs on the yield of various plants, but so far, no report is available on the identification and role of miRNAs in fruit and seed development of almonds. In this study, preliminary analysis by high-throughput sequencing of short RNAs of kernels from the crosses between almond cultivars 'Sefid' × 'Mamaee' (with small and large kernels, respectively) and 'Sefid' × 'P. orientalis' (with small kernels) showed that the expressions of several miRNAs such as Pdu-miR395a-3p, Pdu-miR8123-5p, Pdu-miR482f, Pdu-miR6285, and Pdu-miR396a were significantly different. These miRNAs targeted genes encoding different proteins such as NYFB-3, SPX1, PGSIP3 (GUX2), GH3.9, and BEN1. The result of RT-qPCR revealed that the expression of these genes showed significant differences between the crosses and developmental stages of the seeds, suggesting that these genes might be involved in controlling kernel size because the presence of these miRNAs had a negative effect on their target genes. Pollen source can influence kernel size by affecting hormonal signaling and metabolic pathways through related miRNAs, a phenomenon known as xenia.

Identification of miRNAs-mediated seed and stone-hardening regulatory networks and their signal pathway of GA-induced seedless berries in grapevine (V. vinifera L.)

BACKGROUND

Stone-hardening stage is crucial to the development of grape seed and berry quality. A significant body of evidence supports the important roles of MicroRNAs in grape-berry development, but their specific molecular functions during grape stone-hardening stage remain unclear.

RESULTS

Here, a total of 161 conserved and 85 species-specific miRNAs/miRNAs* (precursor) were identified in grape berries at stone-hardening stage using Solexa sequencing. Amongst them, 30 VvmiRNAs were stone-hardening stage-specific, whereas 52 exhibited differential expression profiles during berry development, potentially participating in the modulation of berry development as verified by their expression patterns. GO and KEGG pathway analysis showed that 13 VvmiRNAs might be involved in the regulation of embryo development, another 11 in lignin and cellulose biosynthesis, and also 28 in the modulation of hormone signaling, sugar, and proline metabolism. Furthermore, the target genes for 4 novel VvmiRNAs related to berry development were validated using RNA Ligase-Mediated (RLM)-RACE and Poly(A) Polymerase-Mediated (PPM)-RACE methods, and their cleavage mainly occurred at the 9th-11th sites from the 5' ends of miRNAs at their binding regions. In view of the regulatory roles of GA in seed embryo development and stone-hardening in grape, we investigated the expression modes of VvmiRNAs and their target genes during GA-induced grape seedless-berry development, and we validated that GA induced the expression of VvmiR31-3p and VvmiR8-5p to negatively regulate the expression levels of CAFFEOYL COENZYME A-3-O-METHYLTRANSFERASE (VvCCoAOMT), and DDB1-CUL4 ASSOCIATED FACTOR1 (VvDCAF1). The series of changes might repress grape stone hardening and embryo development, which might be a potential key molecular mechanism in GA-induced grape seedless-berry development. Finally, a schematic model of miRNA-mediated grape seed and stone-hardening development was proposed.

CONCLUSION

This work identified 30 stone-hardening stage-specific VvmiRNAs and 52 significant differential expression ones, and preliminary interpreted the potential molecular mechanism of GA-induced grape parthenocarpy. GA negatively manipulate the expression of VvCCoAOMT and VvDCAF1 by up-regulation the expression of VvmiR31-3p and VvmiR8-5p, thereby repressing seed stone and embryo development to produce grape seedless berries.

A roadmap for research in octoploid strawberry

The cultivated strawberry (Fragaria × ananassa) is an allo-octoploid species, originating nearly 300 years ago from wild progenitors from the Americas. Since that time the strawberry has become the most widely cultivated fruit crop in the world, universally appealing due to its sensory qualities and health benefits. The recent publication of the first high-quality chromosome-scale octoploid strawberry genome (cv. Camarosa) is enabling rapid advances in genetics, stimulating scientific debate and provoking new research questions. In this forward-looking review we propose avenues of research toward new biological insights and applications to agriculture. Among these are the origins of the genome, characterization of genetic variants, and big data approaches to breeding. Key areas of research in molecular biology will include the control of flowering, fruit development, fruit quality, and plant-pathogen interactions. In order to realize this potential as a global community, investments in genome resources must be continually augmented.

Identification of drought-responsive miRNAs in Hippophae tibetana using high-throughput sequencing

MicroRNAs (miRNAs) play an important role in abiotic stress response in plants. However, the total miRNA profiles (miRNome) and drought-responsive miRNAs in H. tibetana have not been identified. In this study, we present the first report on the miRNome profiles of H. tibetana by high-throughput sequencing technology. 116 known and 4 predicted novel miRNAs were all identified in six H. tibetana samples. Moreover, to reveal the drought-responsive miRNAs in H. tibetana, we compared the miRNA profiles of H. tibetana grown under water sufficiency and drought stress. The results showed that 39 known miRNAs were up-regulated, while 34 miRNAs were downregulated under drought stress. Moreover, the expression of two novel miRNAs (novel_mir_24 and novel_mir_87) showed notable changes in response to drought stress. The target genes of these differentially expressed miRNAs were mainly enriched in cellular process, metabolic process, cell part, and response to stimulus. The identified drought-responsive miRNAs might be used for improving drought tolerance in H. tibetana and other plateau plants.

Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences

Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKC^C -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

Sweet cherry fruit miRNAs and effect of high CO2 on the profile associated with ripening

157 known and 55 novel miRNAs were found in sweet cherry fruit. MiRNA target genes involved in fruit ripening and the differentially expressed miRNAs under CO₂ treatment were identified. MicroRNAs (miRNAs) are short non-coding RNAs and play important functions in many biological processes, including fruit ripening and senescence. In the current study, the high-throughput sequencing and bioinformatics methods were implemented to decipher the miRNAs landscape in sweet cherry fruit. A total of 157 known miRNAs belonging to 50 families and 55 putative novel miRNAs were found. Target genes of the miRNAs were predicted and genes involved in fruit ripening were found, including F-box proteins and TFs such as SPL, TCP, NAC, MYB, ARF and AP2/ERF. And these target genes were further confirmed by degradome sequencing. A regulatory network model was constructed to uncover the miRNAs and their targets involved in fruit ripening and senescence. Importantly, elevated carbon dioxide can significantly postpone the ripening and senescence of sweet cherry fruit and the differentially expressed miRNAs exposed to CO₂ were identified. These miRNAs included miR482j, miR6275, miR164, miR166, miR171, miR393, miR858, miR3627a, miR6284, miR6289 and miR7122b, and some of their functions were linked to fruit ripening. This study was the first report to profile miRNAs in sweet cherry fruit and it would provide more information for further study of miRNA roles in the ripening processes and their regulation mechanism underlying the effects of high carbon dioxide treatment on fruit ripening.

Verification of miRNAs in ginseng decoction by high-throughput sequencing and quantitative real-time PCR

Panax ginseng C. A. Meyer is a precious traditional Chinese medicine that has been clinically used for over thousands of years. In general, ginseng needs to be prepared to ginseng decoction before taking it. MicroRNAs are a class of small (18–24 nt), single-stranded molecules that regulate gene expression at the post-transcriptional level. Considering that ginseng miRNAs may be bioactive compounds, we used Illumina high-throughput sequencing and quantitative real-time PCR (qRT-PCR) to validate the existence of miRNAs in fresh ginseng decoction which have been boiled at high temperature. Our previous studies have demonstrated that there are several miRNAs in fresh ginseng. The roots of fresh Panax ginseng were prepared according to routine methods, from which miRNAs were extracted and sequenced. A total of 43 miRNAs were identified from water decoction by Illumina high-throughput sequencing, belonging to 71 miRNA families. The target genes of these miRNAs were predicted by sequencing, and were annotated by GO, KEGG and Nr databases. The functions of these target genes mainly included plant hormone signal transduction, transcription regulation, macromolecular metabolism and auxin signaling. Nine highly expressed miRNAs (miR159, miR167, miR396, miR166, miR168, miR156, miR165, miR162 and miR394) were verified by qRT-PCR, and the results of Illumina high-throughput sequencing and qRT-PCR were consistent. Results from this study indicate that miRNAs remained stable in P. ginseng after high-temperature boiling. Additionally, Illumina high-throughput sequencing was superior in the acquisition of higher amount of small RNAs.

High-throughput sequencing and differential expression analysis of miRNAs in response to Brassinosteroid treatment in Arabidopsis thaliana

Brassinosteroids are a class of phytohormones that play crucial roles in improving stress tolerance in plants. Many biochemical and physiological changes in response to abiotic stress are related to regulation of gene expression and accumulation of associated proteins. MicroRNAs (miRNAs) are class of small non-coding RNAs that regulate gene expression post-transcriptionally. Roles of these regulatory RNAs in brassinosteroid (BR) signalling have however remained elusive. In this study using high-throughput small RNA sequencing method, we present a comprehensive compilation of BR-induced differentially expressed microRNAs in root and shoots of Arabidopsis thaliana seedlings. We identified 229 known miRNAs belonging to 102 families and 27 novel miRNAs that express in response to exogenous BR treatment. Out of 102 families, miRNAs belonging to known 48 families and out of 27 novel miRNAs, 23 were observed to be differentially expressed in response to BR treatment. Among the conserved miRNAs, all members of miR169 were observed to be downregulated in both shoot and root samples. While, auxin-responsive factors were predicted to be direct targets of some novel miRNAs that are upregulated in shoots and suppressed in roots. The BR-responsive tissue-specific miRNome characterized in this study can be used as a starting point by investigators for functional validation studies that will shed light on the underlying molecular mechanism of BR-mediated stress tolerance at the level of post-transcriptional gene regulation.

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.

Comparative transcriptome analysis uncovers the regulatory functions of long noncoding RNAs in fruit development and color changes of Fragaria pentaphylla

To investigate the molecular mechanism underlying fruit development and color change, comparative transcriptome analysis was employed to generate transcriptome profiles of two typical wild varieties of Fragaria pentaphylla at three fruit developmental stages (green fruit stage, turning stage, and ripe fruit stage). We identified 25,699 long noncoding RNAs (lncRNAs) derived from 25,107 loci in the F. pentaphylla fruit transcriptome, which showed distinct stage- and genotype-specific expression patterns. Time course analysis detected a large number of differentially expressed protein-coding genes and lncRNAs associated with fruit development and ripening in both of the F. pentaphylla varieties. The target genes downregulated in the late stages were enriched in terms of photosynthesis and cell wall organization or biogenesis, suggesting that lncRNAs may act as negative regulators to suppress photosynthesis and cell wall organization or biogenesis during fruit development and ripening of F. pentaphylla. Pairwise comparisons of two varieties at three developmental stages identified 365 differentially expressed lncRNAs in total. Functional annotation of target genes suggested that lncRNAs in F. pentaphylla may play roles in fruit color formation by regulating the expression of structural genes or regulatory factors. Construction of the regulatory network further revealed that the low expression of Fra a and CHS may be the main cause of colorless fruit in F. pentaphylla.

Preharvest Ultraviolet C Treatment Affected Senescence of Stored Strawberry Fruit with a Potential Role of MicroRNAs in the Activation of the Antioxidant System

Recent studies presented preharvest ultraviolet C (UV-C) as an environmentally friendly approach for the management of horticultural crop diseases. The effect of this approach on quality preservation during postharvest storage has not yet been investigated. Strawberry fruit harvested from plants grown with supplemental UV-C were stored at room temperature for 72 h, and their postharvest shelf-life biochemical indicators were evaluated. The involvement of microRNAs (miRNAs) in the activation of UV-C-induced antioxidant systems was investigated. Preharvest UV-C contributed to the preservation of sugar and organic acid and reduced overall lipid peroxidation in strawberry fruit during storage. We found that miR159 and miR398 were downregulated by preharvest UV-C and that their respective targets were upregulated at the early stage of storage with enhancement of the activity of antioxidant enzymes. The initial burst of H₂O₂ and O₂^• - suggested that preharvest UV-C primed the fruit in an antioxidative activated state via reactive-oxygen-species-mediated feedback control with post-transcriptional involvement of miRNAs.

Comparative Transcriptome Profiling Analysis of Red- and White-Fleshed Strawberry (Fragaria�ananassa) Provides New Insight into the Regulation of the Anthocyanin Pathway

Anthocyanins are water-soluble pigments in plants. They confer both economic and healthy profits for humans. To gain a deeper insight into the regulation of anthocyanin biosynthesis in octoploid strawberry (Fragaria�ananassa; Fa), a widely consumed economically important fruit, we performed comparative transcriptomic analysis of red- and white-fleshed strawberry cultivars in two ripening stages. In total, 365,455 non-redundant transcripts were assembled from the RNA sequencing (RNAseq) data. Of this collection, 377 were annotated as putative anthocyanin-related transcripts. Differential expression analysis revealed that 57 anthocyanin biosynthesis transcripts were down-regulated, and 89 transcription factors (TFs) were either down- or up-regulated under anthocyanin deficiency. Additionally, amongst the 50,601 putative long non-coding RNAs (lncRNAs) identified here, 68 lncRNAs were differentially expressed and co-expressed with differentially expressed anthocyanin-related mRNAs; 2,070 co-expressing lncRNA-mRNA pairs were generated. Expression profile analysis revealed that it was the limited expression of FaF3'H (flavonoid 3'-hydroxylase) that blocked the cyanidin 3-glucoside accumulation in the two investigated strawberry cultivars. This was further supported by a transient overexpression experiment with FaMYB10. The down-regulated lncRNAs might participate in anthocyanin regulation by acting as targets for microRNAs (miRNAs). The level of competitive intensity in miRNA and lncRNA for the same mRNA targets was probably lower in the white-fleshed strawberries, which can release the repression effect of the mRNAs in red-fleshed strawberry as a result. This study for the first time presents lncRNAs related to anthocyanins in strawberries, provides new insights into the anthocyanin regulatory network and also lays the foundation for identifying new anthocyanin regulators in strawberry.

Expansion of Capsicum annuum fruit is linked to dynamic tissue-specific differential expression of miRNA and siRNA profiles

Small regulatory RNAs, such as microRNAs (miRNAs) and small interfering RNAs (siRNAs) have emerged as important transcriptional and post-transcriptional regulators controlling a wide variety of physiological processes including fruit development. Data are, however, limited for their potential roles in developmental processes determining economically important traits of crops. The current study aimed to discover and characterize differentially expressed miRNAs and siRNAs in sweet pepper (Capsicum annuum) during fruit expansion. High-throughput sequencing was employed to determine the small regulatory RNA expression profiles in various fruit tissues, such as placenta, seed, and flesh at 28 and 40 days after anthesis. Comparative differential expression analyses of conserved, already described and our newly predicted pepper-specific miRNAs revealed that fruit expansion is accompanied by an increasing level of miRNA-mediated regulation of gene expression. Accordingly, ARGONAUTE1 protein, the primary executor of miRNA-mediated regulation, continuously accumulated to an extremely high level in the flesh. We also identified numerous pepper-specific, heterochromatin-associated 24-nt siRNAs (hetsiRNAs) which were extremely abundant in the seeds, as well as 21-nt and 24-nt phased siRNAs (phasiRNAs) that were expressed mainly in the placenta and the seeds. This work provides comprehensive tissue-specific miRNA and siRNA expression landscape for a developing pepper fruit. We identified several novel, abundantly expressing tissue- and pepper-specific small regulatory RNA species. Our data show that fruit expansion is associated with extensive changes in sRNA abundance, raising the possibility that manipulation of sRNA pathways may be employed to improve the quality and quantity of the pepper fruit.

Identification and characterization of known and novel microRNAs in strawberry fruits induced by Botrytis cinerea

MicroRNAs are endogenous small non-coding RNAs that negatively regulate mRNAs, mainly at the post-transcriptional level, and play an important role in resistance response of plants. To date, there are few reports on resistance response of strawberry miRNAs to pathogens. In this study, using high-throughput sequencing, 134 conserved and 35 novel miRNAs were identified in six libraries within the treatment of Botrytis cinerea. A total 497 potential target genes were predicted using Fragaria vesca genome. Most of the differential expressed miRNAs in strawberry fruits were up-regulated in early libraries and down-regulated in late libraries. PIRL, the target gene of miR5290a, showed the opposite expressed trend compared with miR5290 from T1 to T3 libraries, and functional analysis of the PIRL gene shows that it has obvious resistance to B. cinerea in the strawberry fruits with overexpressed PIRL gene. We speculate that miR5290a negatively regulates its target gene PIRL to increase resistance to pathogen infection, and further analysis of PIRL function is meaningful for studying the plant-pathogen relationship and improving strawberry fruit quality and yield.

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.

miRNA alteration is an important mechanism in sugarcane response to low-temperature environment

BACKGROUND

Cold is a major abiotic stress limiting the production of tropical and subtropical crops in new production areas. Sugarcane (Saccharum spp.) originates from the tropics but is cultivated primarily in the sub-tropics where it frequently encounters cold stress. Besides regulating plant growth, miRNAs play an important role in environmental adaption.

RESULTS

In this study, a total of 412 sugarcane miRNAs, including 261 known and 151 novel miRNAs, were obtained from 4 small RNA libraries through the Illumina sequencing method. Among them, 62 exhibited significant differential expression under cold stress, with 34 being upregulated and 28 being downregulated. The expression of 13 miRNAs and 12 corresponding targets was validated by RT-qPCR, with the majority being consistent with the sequencing data. GO and KEGG analysis indicated that these miRNAs were involved in stress-related biological pathways. To further investigate the involvement of these miRNAs in tolerance to abiotic stresses, sugarcane miR156 was selected for functional analysis. RT-qPCR revealed that miR156 levels increased in sugarcane during cold, salt and drought stress treatments. Nicotiana benthamiana plants transiently overexpressing miR156 exhibited better growth status, lower ROS levels, higher anthocyanin contents as well as the induction of some cold-responsive genes, suggesting its positive role in the plant cold stress response.

CONCLUSIONS

This study provides a global view of the association of miRNA expression with the sugarcane response to cold stress. The findings have enriched the present miRNA resource and have made an attempt to verify the involvement of miR156 in plant response to cold stress.

Genome-wide identification of Hami melon miRNAs with putative roles during fruit development

MicroRNAs represent a family of small endogenous, non-coding RNAs that play critical regulatory roles in plant growth, development, and environmental stress responses. Hami melon is famous for its attractive flavor and excellent nutritional value, however, the mechanisms underlying the fruit development and ripening remains largely unknown. Here, we performed small RNA sequencing to investigate the roles of miRNAs during Hami melon fruit development. Two batches of flesh samples were collected at four fruit development stages. Small RNA sequencing yielded a total of 54,553,424 raw reads from eight libraries. 113 conserved miRNAs belonging to 30 miRNA families and nine novel miRNAs comprising nine miRNA families were identified. The expression of 42 conserved miRNAs and three Hami melon-specific miRNAs significantly changed during fruit development. Furthermore, 484 and 124 melon genes were predicted as putative targets of 29 conserved and nine Hami melon-specific miRNA families, respectively. GO enrichment analysis were performed on target genes, "transcription, DNA-dependent", "rRNA processing", "oxidation reduction", "signal transduction", "regulation of transcription, DNA-dependent", and "metabolic process" were the over-represented biological process terms. Cleavage sites of six target genes were validated using 5' RACE. Our results present a comprehensive set of identification and characterization of Hami melon fruit miRNAs and their potential targets, which provide valuable basis towards understanding the regulatory mechanisms in programmed process of normal Hami fruit development and ripening. Specific miRNAs could be selected for further research and applications in breeding practices.

New technologies accelerate the exploration of non-coding RNAs in horticultural plants

Non-coding RNAs (ncRNAs), that is, RNAs not translated into proteins, are crucial regulators of a variety of biological processes in plants. While protein-encoding genes have been relatively well-annotated in sequenced genomes, accounting for a small portion of the genome space in plants, the universe of plant ncRNAs is rapidly expanding. Recent advances in experimental and computational technologies have generated a great momentum for discovery and functional characterization of ncRNAs. Here we summarize the classification and known biological functions of plant ncRNAs, review the application of next-generation sequencing (NGS) technology and ribosome profiling technology to ncRNA discovery in horticultural plants and discuss the application of new technologies, especially the new genome-editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) systems, to functional characterization of plant ncRNAs.

Identification of miRNAs and their targets through high-throughput sequencing and degradome analysis in male and female Asparagus officinalis

BACKGROUND

MicroRNAs (miRNAs), a class of non-coding small RNAs (sRNAs), regulate various biological processes. Although miRNAs have been identified and characterized in several plant species, miRNAs in Asparagus officinalis have not been reported. As a dioecious plant with homomorphic sex chromosomes, asparagus is regarded as an important model system for studying mechanisms of plant sex determination.

RESULTS

Two independent sRNA libraries from male and female asparagus plants were sequenced with Illumina sequencing, thereby generating 4.13 and 5.88 million final clean reads, respectively. Both libraries predominantly contained 24-nt sRNAs, followed by 21-nt sRNAs. Further analysis identified 154 conserved miRNAs, which belong to 26 families, and 39 novel miRNA candidates seemed to be specific to asparagus. Comparative profiling revealed that 63 miRNAs exhibited significant differential expression between male and female plants, which was confirmed by real-time quantitative PCR analysis. Among them, 37 miRNAs were significantly up-regulated in the female library, whereas the others were preferentially expressed in the male library. Furthermore, 40 target mRNAs representing 44 conserved and seven novel miRNAs were identified in asparagus through high-throughput degradome sequencing. Functional annotation showed that these target mRNAs were involved in a wide range of developmental and metabolic processes.

CONCLUSIONS

We identified a large set of conserved and specific miRNAs and compared their expression levels between male and female asparagus plants. Several asparagus miRNAs, which belong to the miR159, miR167, and miR172 families involved in reproductive organ development, were differentially expressed between male and female plants, as well as during flower development. Consistently, several predicted targets of asparagus miRNAs were associated with floral organ development. These findings suggest the potential roles of miRNAs in sex determination and reproductive developmental processes in asparagus.

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.

The Influence of Genotype and Environment on Small RNA Profiles in Grapevine Berry

Understanding the molecular mechanisms involved in the interaction between the genetic composition and the environment is crucial for modern viticulture. We approached this issue by focusing on the small RNA transcriptome in grapevine berries of the two varieties Cabernet Sauvignon and Sangiovese, growing in adjacent vineyards in three different environments. Four different developmental stages were studied and a total of 48 libraries of small RNAs were produced and sequenced. Using a proximity-based pipeline, we determined the general landscape of small RNAs accumulation in grapevine berries. We also investigated the presence of known and novel miRNAs and analyzed their accumulation profile. The results showed that the distribution of small RNA-producing loci is variable between the two cultivars, and that the level of variation depends on the vineyard. Differently, the profile of miRNA accumulation mainly depends on the developmental stage. The vineyard in Riccione maximizes the differences between the varieties, promoting the production of more than 1000 specific small RNA loci and modulating their expression depending on the cultivar and the maturation stage. In total, 89 known vvi-miRNAs and 33 novel vvi-miRNA candidates were identified in our samples, many of them showing the accumulation profile modulated by at least one of the factors studied. The in silico prediction of miRNA targets suggests their involvement in berry development and in secondary metabolites accumulation such as anthocyanins and polyphenols.

Identification of microRNAs and their targets in Finger millet by high throughput sequencing

MicroRNAs are short non-coding RNAs which play an important role in regulating gene expression by mRNA cleavage or by translational repression. The majority of identified miRNAs were evolutionarily conserved; however, others expressed in a species-specific manner. Finger millet is an important cereal crop; nonetheless, no practical information is available on microRNAs to date. In this study, we have identified 95 conserved microRNAs belonging to 39 families and 3 novel microRNAs by high throughput sequencing. For the identified conserved and novel miRNAs a total of 507 targets were predicted. 11 miRNAs were validated and tissue specificity was determined by stem loop RT-qPCR, Northern blot. GO analyses revealed targets of miRNA were involved in wide range of regulatory functions. This study implies large number of known and novel miRNAs found in Finger millet which may play important role in growth and development.

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.

miRNA expression during prickly pear cactus fruit development

miRNAs are a class of small non-coding RNAs that regulate gene expression. They are involved in the control of many developmental processes, including fruit development. The increasing amount of information on miRNAs, on their expression, abundance, and conservation between various species, provides a new opportunity to study the role of miRNAs in non-model plant species. In this work, we used a combination of Northern blot and tissue print hybridization analysis to identify conserved miRNAs expressed during prickly pear cactus (Opuntia ficus indica) fruit development. Comparative profiling detected the expression of 34 miRNAs, which were clustered in three different groups that were associated with the different phases of fruit development. Variation in the level of miRNA expression was observed. Gradual expression increase of several miRNAs was observed during fruit development, including miR164. miR164 was selected for stem-loop RT-PCR and for a detailed spatial-temporal expression analysis. At early floral stages, miR164 was mainly localized in meristematic tissues, boundaries and fusion zones, while it was more homogenously expressed in fruit tissues. Our results provide the first evidence of miRNA expression in the prickly pear cactus and provide the basis for future research on miRNAs in Opuntia. Moreover, our analyses suggest that miR164 plays different roles during prickly pear cactus fruit development.

De novo transcriptome assembly, gene annotation, marker development, and miRNA potential target genes validation under abiotic stresses in Oenanthe javanica

Oenanthe javanica is an aquatic perennial herb with known medicinal properties and an edible vegetable with high vitamin and mineral content. The understanding of the biology of O. javanica is limited by the absence of information on its genome, transcriptome, and small RNA. In this study, transcriptome sequencing and small RNA sequencing were performed to annotate function genes, develop SSR markers and analyze potential target genes of miRNAs in O. javanica. All reads with total nucleotides number of 1,440,321,408 bp were assembled into 58,072 transcripts and 40,208 unigenes. A total of 1,233 SSRs were identified from O. javanica. Generated unigenes were aligned against seven databases and annotated with functions. A total of 29 potential targets were predicted. Expression of 10 miRNAs and their corresponding target genes under abiotic stresses (heat, cold, salinity, and drought) was validated. All ten miRNAs were confirmed to response to abiotic stresses. A pair of miRNA and its target gene was found. This study can serve as a valuable resource for future studies on O. javanica, which may focus on novel gene discovery, SSR development, gene mapping, and miRNA-affected processes and pathways. This can promote the development of the useful medicinal properties of O. javanica in medical science.

Fragaria: a genus with deep historical roots and ripe for evolutionary and ecological insights

The cultivated strawberry, Fragaria ×ananassa, is one of the youngest domesticated plants. Its 18th century origin via hybridization in Europe between the North American F. virginiana and the South American F. chiloensis was documented by the botanist Antoine Nicolas Duchesne. His 1766 "Natural History of Strawberries" is an extraordinary work that integrates fundamental discoveries on the biology, ecology, and phylogeny of Fragaria with applied information on cultivation and ethnobotanical uses, serving as an inspiration for current research in the genus. Fragaria species exhibit the full range of sexual systems in the gynodioecy pathway from hermaphroditism to dioecy (and back again), as well as variation in self-compatibility, and evidence of sex chromosomes with female heterogamety. The genus is also characterized by interspecific hybridization and polyploidy, with a natural range of ploidy levels from diploids to decaploids. This biological diversity, combined with the availability of genomic resources and the ease of growing and experimenting with the plants, makes Fragaria a very attractive system for ecological and evolutionary genomics. The goal of this review is to introduce Fragaria as a model genus and to provide a roadmap for future integrative research. These research directions will deepen our understanding of the ecological and evolutionary context that shaped the ancestors of the cultivated strawberry, not only providing information that can be applied to efforts to shape the future of this important fruit crop but also our understanding of key transitions in plant evolution.

Advances in identification and validation of plant microRNAs and their target genes

Developments in the field of molecular biology and genetics, such as microarray, gene transfer and discovery of small regulatory RNAs, have led to significant advances in plant biotechnology. Among the small RNAs, microRNAs (miRNAs) have elicited much interest as key post-transcriptional regulators in eukaryotic gene expression. Advances in genome and transcriptome sequencing of plants have facilitated the generation of a huge wealth of sequence information that can find much use in the discovery of novel miRNAs and their target genes. In this review, we present an overview of the developments in the strategies and methods used to identify and study miRNAs, their target genes and the mechanisms by which these miRNAs interact with their target genes since the discovery of the first miRNA. The approaches discussed include both reverse and forward genetics. We observed that despite the availability of advanced methods, certain limitations ranging from the cost of materials, equipment and personnel to the availability of genome sequences for many plant species present a number of challenges for the development and utilization of modern scientific methods for the elucidation and development of miRNAs in many important plant species.

Identification of the conserved and novel miRNAs in Mulberry by high-throughput sequencing

miRNAs are a class of non-coding endogenous small RNAs. They play vital roles in plant growth, development, and response to biotic and abiotic stress by negatively regulating genes. Mulberry trees are economically important species with multiple uses. However, to date, little is known about mulberry miRNAs and their target genes. In the present study, three small mulberry RNA libraries were constructed and sequenced using high-throughput sequencing technology. Results showed 85 conserved miRNAs belonging to 31 miRNA families and 262 novel miRNAs at 371 loci. Quantitative real-time PCR (qRT-PCR) analysis confirmed the expression pattern of 9 conserved and 5 novel miRNAs in leaves, bark, and male flowers. A total of 332 potential target genes were predicted to be associated with these 113 novel miRNAs. These results provide a basis for further understanding of mulberry miRNAs and the biological processes in which they are involved.

Grafting-responsive miRNAs in cucumber and pumpkin seedlings identified by high-throughput sequencing at whole genome level

Grafting is an important agricultural technique widely used for improving growth, yields and tolerance of crops to abiotic and biotic stresses. As one type of endogenous, non-coding small RNAs, microRNAs (miRNAs) regulate development and responsiveness to biotic and abiotic stresses by negatively mediating expression of target genes at the post-transcriptional level. However, there have been few detailed studies to evaluate the role of miRNAs in mediation of grafting-induced physiological processes in plants. Cucumis sativus and Cucurbita moschata are important vegetables worldwide. We constructed eight small RNA libraries from leaves and roots of seedlings that were grafted in the following four ways: (1) hetero-grafting, using cucumber as scion and pumpkin as rootstock; (2) hetero-grafting, with pumpkin as scion and cucumber as rootstock; (3) auto-grafting of cucumbers and (4) auto-grafting of pumpkins. High-throughput sequencing was employed, and more than 120 million raw reads were obtained. We annotated 112 known miRNAs belonging to 40 miRNA families and identified 48 new miRNAs in the eight libraries, and the targets of these known and novel miRNAs were predicted by bioinformatics. Grafting led to changes in expression of most miRNAs and their predicted target genes, suggesting that miRNAs may play significant roles in mediating physiological processes of grafted seedlings by regulating the expression of target genes. The potential role of the grafting-responsive miRNAs in seedling growth and long-distance transport of miRNA was discussed. These results are useful for functional characterization of miRNAs in mediation of grafting-dependent physiological processes.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

miRNAs expression profile in bast of ramie elongation phase and cell wall thickening and end wall dissolving phase

MicroRNAs (miRNAs) are 20-24 nucleotide long non-coding RNAs that play critical regulatory roles during plant development, organ morphogenesis, and cell fate determination and differentiation. In this study, miRNA microarray chips were used to explore the expression profile of ramie miRNAs between the bast of fiber elongation phase and those of cell wall thickening and end wall dissolving phase. There are 150 and 148 credible miRNAs in the bast of fiber elongation phase and cell wall thickening and end wall dissolving phase, respectively. These miRNAs distributed in 27 species and mainly concentrated in nine species. Analysis showed that 51 miRNAs were differentially expressed: 27 up-regulated (miR166, miR172, miR396, miR482, miR894 and miR2911 families) and 24 down-regulated (miR156, miR159, miR164, miR319 and miR1450 families) in the bast of fiber elongation phase compared with the bast of cell wall thickening and end wall dissolving phase. To further confirm our results, we examined the expression of three miRNAs (zma-miR172b*, pvu-miR482 and vvi-172a) by quantitative real-time reverse transcriptase-PCR. Our results will provide a molecular basis for future research miRNA function on ramie genetics and breeding.

Deep sequencing discovery of novel and conserved microRNAs in wild type and a white-flesh mutant strawberry

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression by base pairing to mRNA target sequences, and play crucial roles in plant development and stress responses. The knowledge on post-transcriptional regulation by miRNAs in strawberry is rather limited so far. In order to understand the role of miRNA in the molecular control during strawberry fruit development, small RNA libraries were constructed from fruits at the turning stage of strawberry cultivar 'Sachinoka' and its white-flesh mutant by using the Solexa platform. One hundred and twenty conserved miRNAs belonging to 27 miRNA families and 33 putative novel strawberry miRNAs were identified in both libraries. Their target genes were predicted using the Fragaria vesca genome. Nine of all miRNAs showed significant expression differences between two types of samples. Four miRNAs were up-regulated and five were down-regulated in white-flesh mutant. The sequencing results were partially validated by quantitative RT-PCR. Among them, the expression of miR399a shows the biggest change between the two samples. The prediction of its target gene showed that miR399 may play an important role in phosphate homeostasis of strawberry fruits. Furthermore, we deduce that the expression of miR399 has negative correlation with the content of sugars.

Identification and characterization of novel and conserved microRNAs in radish (Raphanus sativus L.) using high-throughput sequencing

MicroRNAs (miRNAs) are endogenous, non-coding, small RNAs that play significant regulatory roles in plant growth, development, and biotic and abiotic stress responses. To date, a great number of conserved and species-specific miRNAs have been identified in many important plant species such as Arabidopsis, rice and poplar. However, little is known about identification of miRNAs and their target genes in radish (Raphanus sativus L.). In the present study, a small RNA library from radish root was constructed and sequenced using the high-throughput Solexa sequencing. Through sequence alignment and secondary structure prediction, a total of 545 conserved miRNA families as well as 15 novel (with their miRNA* strand) and 64 potentially novel miRNAs were identified. Quantitative real-time PCR (qRT-PCR) analysis confirmed that both conserved and novel miRNAs were expressed in radish, and some of them were preferentially expressed in certain tissues. A total of 196 potential target genes were predicted for 42 novel radish miRNAs. Gene ontology (GO) analysis showed that most of the targets were involved in plant growth, development, metabolism and stress responses. This study represents a first large-scale identification and characterization of radish miRNAs and their potential target genes. These results could lead to the further identification of radish miRNAs and enhance our understanding of radish miRNA regulatory mechanisms in diverse biological and metabolic processes.