Deep sequencing of grapevine flower and berry short RNA library for discovery of novel microRNAs and validation of precise sequences of grapevine microRNAs deposited in miRBase

Participation of miR165a in the Phytochrome Signal Transduction in Maize (Zea mays L.) Leaves under Changing Light Conditions

The involvement of the microRNA miR165a in the light-dependent mechanisms of regulation of target genes in maize (Zea mays) has been studied. The light-induced change in the content of free miR165a was associated with its binding by the AGO10 protein and not with a change in the rate of its synthesis from the precursor. The use of knockout Arabidopsis plants for the phytochrome A and B genes demonstrated that the presence of an active form of phytochrome B causes an increase in the level of the RNA-induced silencing miR165a complex, which triggers the degradation of target mRNAs. The two fractions of vesicles from maize leaves, P40 and P100 that bind miR165a, were isolated by ultracentrifugation. The P40 fraction consisted of larger vesicles of the size >0.170 µm, while the P100 fraction vesicles were <0.147 µm. Based on the quantitative PCR data, the predominant location of miR165a on the surface of extracellular vesicles of both fractions was established. The formation of the active form of phytochrome upon the irradiation of maize plants with red light led to a redistribution of miR165a, resulting in an increase in its proportion inside P40 vesicles and a decrease in P100 vesicles.

Assessment of 'Cabernet Sauvignon' Grape Quality Half-Véraison to Maturity for Grapevines Grown in Different Regions

Grapes are widely cultivated around the world and their quality has distinct regional characteristics. In this study, the qualitative characteristics of the 'Cabernet Sauvignon' grape variety in seven regions, from half-véraison to maturity, were analyzed comprehensively at physiological and transcriptional levels. The results indicated that the quality traits of 'Cabernet Sauvignon' grapes in different regions were significantly different with obvious regionality. Total phenols, anthocyanins, and titratable acids were the main factors of the regionality of berry quality, which were very sensitive to changes in the environment. It should be noted that the changes in titrating acids and total anthocyanin of berries vary greatly from half-véraison to maturity between regions. Moreover, the transcriptional analysis showed that the co-expressed genes between regions characterized the core transcriptome of berry development, while the unique genes of each region reflected the regionality of berries. The differentially expressed genes (DEGs) between half-véraison and maturity can be used to demonstrate that the environment of the regions could promote or inhibit gene expression. The functional enrichment suggested that these DEGs help to understand the interpretation of the plasticity of the quality composition of grapes according to the environment. Taken together, the information generated by this study could contribute to the development of viticultural practices aimed at making better use of native varieties for the development of wines with regional characteristics.

Advances in understanding cold tolerance in grapevine

Grapevine (Vitis ssp.) is a deciduous perennial fruit crop, and the canes and buds of grapevine should withstand low temperatures annually during winter. However, the widely cultivated Vitis vinifera is cold-sensitive and cannot survive the severe winter in regions with extremely low temperatures, such as viticulture regions in northern China. By contrast, a few wild Vitis species like V. amurensis and V. riparia exhibit excellent freezing tolerance. However, the mechanisms underlying grapevine cold tolerance remain largely unknown. In recent years, much progress has been made in elucidating the mechanisms, owing to the advances in sequencing and molecular biotechnology. Assembly of grapevine genomes together with resequencing and transcriptome data enable researchers to conduct genomic and transcriptomic analyses in various grapevine genotypes and populations to explore genetic variations involved in cold tolerance. In addition, a number of pivotal genes have been identified and functionally characterized. In this review, we summarize recent major advances in physiological and molecular analyses of cold tolerance in grapevine and put forward questions in this field. We also discuss the strategies for improving the tolerance of grapevine to cold stress. Understanding grapevine cold tolerance will facilitate the development of grapevines for adaption to global climate change.

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.

Tissue and regional expression patterns of dicistronic tRNA-mRNA transcripts in grapevine (Vitis vinifera) and their evolutionary co-appearance with vasculature in land plants

Transfer RNAs (tRNA) are crucial adaptor molecules between messenger RNA (mRNA) and amino acids. Recent evidence in plants suggests that dicistronic tRNA-like structures also act as mobile signals for mRNA transcripts to move between distant tissues. Co-transcription is not a common feature in the plant nuclear genome and, in the few cases where polycistronic transcripts have been found, they include non-coding RNA species, such as small nucleolar RNAs and microRNAs. It is not known, however, the extent to which dicistronic transcripts of tRNA and mRNAs are expressed in field-grown plants, or the factors contributing to their expression. We analysed tRNA-mRNA dicistronic transcripts in the major horticultural crop grapevine (Vitis vinifera) using a novel pipeline developed to identify dicistronic transcripts from high-throughput RNA-sequencing data. We identified dicistronic tRNA-mRNA in leaf and berry samples from 22 commercial vineyards. Of the 124 tRNA genes that were expressed in both tissues, 18 tRNA were expressed forming part of 19 dicistronic tRNA-mRNAs. The presence and abundance of dicistronic molecules was tissue and geographic sub-region specific. In leaves, the expression patterns of dicistronic tRNA-mRNAs significantly correlated with tRNA expression, suggesting that their transcriptional regulation might be linked. We also found evidence of syntenic genomic arrangements of tRNAs and protein-coding genes between grapevine and Arabidopsis thaliana, and widespread prevalence of dicistronic tRNA-mRNA transcripts among vascular land plants but no evidence of these transcripts in non-vascular lineages. This suggests that the appearance of plant vasculature and tRNA-mRNA occurred concurrently during the evolution of land plants.

Integrative Analysis of Gene Expression and miRNAs Reveal Biological Pathways Associated with Bud Paradormancy and Endodormancy in Grapevine

Transition of grapevine buds from paradormancy to endodormancy is coordinated by changes in gene expression, phytohormones, transcription factors, and other molecular regulators, but the mechanisms involved in transcriptional and post-transcriptional regulation of dormancy stages are not well delineated. To identify potential regulatory targets, an integrative analysis of differential gene expression profiles and their inverse relationships with miRNA abundance was performed in paradormant (long day (LD) 15 h) or endodormant (short day (SD), 13 h) Vitis riparia buds. There were 400 up- and 936 downregulated differentially expressed genes in SD relative to LD budsGene set and gene ontology enrichment analysis indicated that hormone signaling and cell cycling genes were downregulated in SD relative to LD buds. miRNA abundance and inverse expression analyses of miRNA target genes indicated increased abundance of miRNAs that negatively regulate genes involved with cell cycle and meristem development in endodormant buds and miRNAs targeting starch metabolism related genes in paradormant buds. Analysis of interactions between abundant miRNAs and transcription factors identified a network with coinciding regulation of cell cycle and epigenetic regulation related genes in SD buds. This network provides evidence for cross regulation occurring between miRNA and transcription factors both upstream and downstream of MYB3R1.

Characterization of the regulation mechanism of grapevine microRNA172 family members during flower development

BACKGROUND

Grapevine (Vitis vinifera L.), which has important nutritional values and health benefits, is one of the most economically important fruit crops cultivated worldwide. Several studies showed a large number of microRNAs (VvmiRNAs) involved in the modulation of grape growth and development, and many VvmiRNA families have multiple members. However, the way by which various members from the same miRNA family work is unclear, particularly in grapes.

RESULTS

In this study, an important conserved VvmiR172 family (VvmiR172s) and their targets were set as a good example for elucidating the interaction degree, mechanism, and spatio-temporal traits of diverse members from the same miRNA family. miR-RACE and Stem-loop RT-PCR were used to identify the spatio-temporal expressions of various members of VvmiR172s; together with RLM-RACE, PPM-RACE, Western blot, transgenic technologies, their interaction degree, and regulation mechanism were further validated. The expression of VvmiR172c was significantly higher than that of VvmiR172a, b, and d and showed a positive correlation with the abundance of VvAP2 cleavage products. These findings indicated that VvmiR172c might be one of the main action factors of the VvmiR172 family in flower development. The ability of VvmiR172c to cleave target genes differed due to divergence in complementary degree with VvAP2 and expression levels of various members. In VvmiR172 transgenic lines, we observed that 35S::VvmiR172c resulted in the earliest and abundant flowering, indicating the strong function of VvmiR172c. In contrast, the non-significant phenotypic changes were detected in the VvAP2 transgenic lines. The qRT-PCR and Western bolt results demonstrated that VvmiR172c plays a major role in targeting VvAP2.

CONCLUSIONS

VvmiR172 up-regulated the expression of NtFT and decreased the expression of NtFLC. The up/down regulation of VvmiR172c was the most pronounced. The functions of four VvmiR172 members in grape differed, and miR172c had the strongest regulation on AP2.

Genome-wide identification and characterization of gibberellin metabolic and signal transduction (GA MST) pathway mediating seed and berry development (SBD) in grape (Vitis vinifera L.)

BACKGROUND

Grape is highly sensitive to gibberellin (GA), which is crucial during seed and berry development (SBD) either by itself or by interacting with other hormones, such as auxin, Abscisic acid (ABA), and Cytokinin (CK). However, no systematic analysis of GA metabolic and signal transduction (MST) pathway has been undertaken in grapevine.

RESULTS

In this study, total endogenous GA3 content significantly decreased during SBD, and a total of 48 known genes in GA metabolic (GAM; 31) and signal transduction (ST; 17) pathways were identified in this process. In the GAM pathway, out of 31 genes, VvGA20ox1-1, VvGA3ox4-1, and VvGA2ox1-1 may be the major factors interacting at the green-berry stage (GBS) accompanied with higher accumulation rate. GA biosynthesis was greater than GA inactivation at GBS, confirming the importance of seeds in GA synthesis. The visible correlation between endogenous GA3 content and gene expression profiles suggested that the transcriptional regulation of GA biosynthesis pathway genes was a key mechanism of GA accumulation at the stone-hardening stage (SHS). Interestingly, we observed a negative feedback regulation between VvGA3oxs-VvGAI1-4, VvGA2oxs-VvGAI1-4, and VvGID1B-VvGAI1-4 in maintaining the balance of GA3 content in berries. Moreover, 11 miRNAs may be involved in the modulation of GA MST pathway by mediating their target genes, such as VvGA3ox, VvGID1B, and VvGAMYB. Many genes in auxin, ABA, and CK MST pathways were further identified and found to have a special pattern in the berry, and the crosstalk between GA and these hormones may modulate the complex process during SBD through the interaction gene network of the multihormone pathway. Lastly, based on the expression characterization of multihormone MST pathway genes, a proposed model of the GA-mediated multihormone regulatory network during SBD was proposed.

CONCLUSIONS

Our results provided novel insights into GA-mediated regulatory networks during SBD in grape. The complexity of GA-mediated multihormone ST in SBD was also elucidated, thereby providing valuable information for future functional characterizations of specific genes in grape.

Transcriptomic and biochemical investigations support the role of rootstock-scion interaction in grapevine berry quality

Background

In viticulture, rootstock genotype plays a critical role to improve scion physiology, berry quality and to adapt grapevine (Vitis vinifera L.) to different environmental conditions. This study aimed at investigating the effect of two different rootstocks (1103 Paulsen - P - and Mgt 101–14 - M) in comparison with not grafted plants - NGC - on transcriptome (RNA-seq and small RNA-seq) and chemical composition of berry skin in Pinot noir, and exploring the influence of rootstock-scion interaction on grape quality. Berry samples, collected at veraison and maturity, were investigated at transcriptional and biochemical levels to depict the impact of rootstock on berry maturation.

Results

RNA- and miRNA-seq analyses highlighted that, at veraison, the transcriptomes of the berry skin are extremely similar, while variations associated with the different rootstocks become evident at maturity, suggesting a greater diversification at transcriptional level towards the end of the ripening process. In the experimental design, resembling standard agronomic growth conditions, the vines grafted on the two different rootstocks do not show a high degree of diversity. In general, the few genes differentially expressed at veraison were linked to photosynthesis, putatively because of a ripening delay in not grafted vines, while at maturity the differentially expressed genes were mainly involved in the synthesis and transport of phenylpropanoids (e.g. flavonoids), cell wall loosening, and stress response. These results were supported by some differences in berry phenolic composition detected between grafted and not grafted plants, in particular in resveratrol derivatives accumulation.

Conclusions

Transcriptomic and biochemical data demonstrate a stronger impact of 1103 Paulsen rootstock than Mgt 101–14 or not grafted plants on ripening processes related to the secondary metabolite accumulations in berry skin tissue. Interestingly, the MYB14 gene, involved in the feedback regulation of resveratrol biosynthesis was up-regulated in 1103 Paulsen thus supporting a putative greater accumulation of stilbenes in mature berries.

Understanding the molecular mechanisms underlying graft success in grapevine

BACKGROUND

Grafting is an intensive commercial practice required to protect the European grapevine against the Phylloxera pest. Rootstocks resistant to this pest are hybrids of American vine species with different levels of compatibility with European Vitis vinifera varieties. Aiming to understand what drives grafting compatibility in grapevine, a transcriptomic approach was used to search for master regulators of graft success. Two scion/rootstock combinations, with different levels of compatibility, were compared in a nursery-grafting context at two stages, at 21 and 80 days after grafting.

RESULTS

In the most compatible combination, an earlier and higher expression of genes signaling the metabolic and hormonal pathways as well as a reduced expression of genes of the phenolic metabolism and of the oxidative stress response was observed. At 80 days after grafting a higher expression of transcription factors regulating vascular maintenance, differentiation and proliferation was obtained in the most compatible combination. Moreover, lower expression levels of microRNAs potentially targeting important transcription factors related to plant development was observed in the more compatible combination when compared to the less compatible one.

CONCLUSION

In this context, a set of regulators was selected as potential expression markers for early prediction of a compatible grafting.

Understanding the molecular mechanisms underlying graft success in grapevine

BACKGROUND

Grafting is an intensive commercial practice required to protect the European grapevine against the Phylloxera pest. Rootstocks resistant to this pest are hybrids of American vine species with different levels of compatibility with European Vitis vinifera varieties. Aiming to understand what drives grafting compatibility in grapevine, a transcriptomic approach was used to search for master regulators of graft success. Two scion/rootstock combinations, with different levels of compatibility, were compared in a nursery-grafting context at two stages, at 21 and 80 days after grafting.

RESULTS

In the most compatible combination, an earlier and higher expression of genes signaling the metabolic and hormonal pathways as well as a reduced expression of genes of the phenolic metabolism and of the oxidative stress response was observed. At 80 days after grafting a higher expression of transcription factors regulating vascular maintenance, differentiation and proliferation was obtained in the most compatible combination. Moreover, lower expression levels of microRNAs potentially targeting important transcription factors related to plant development was observed in the more compatible combination when compared to the less compatible one.

CONCLUSION

In this context, a set of regulators was selected as potential expression markers for early prediction of a compatible grafting.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.

miRVIT: A Novel miRNA Database and Its Application to Uncover Vitis Responses to Flavescence dorée Infection

Micro(mi)RNAs play crucial roles in plant developmental processes and in defense responses to biotic and abiotic stresses. In the last years, many works on small RNAs in grapevine (Vitis spp.) were published, and several conserved and putative novel grapevine-specific miRNAs were identified. In order to reorganize the high quantity of available data, we produced "miRVIT," the first database of all novel grapevine miRNA candidates characterized so far, and still not deposited in miRBase. To this aim, each miRNA accession was renamed, repositioned in the last version of the grapevine genome, and compared with all the novel and conserved miRNAs detected in grapevine. Conserved and novel miRNAs cataloged in miRVIT were then used for analyzing Vitis vinifera plants infected by Flavescence dorée (FD), one of the most severe phytoplasma diseases affecting grapevine. The analysis of small RNAs from healthy, recovered (plants showing spontaneous and stable remission of symptoms), and FD-infected "Barbera" grapevines showed that FD altered the expression profiles of several miRNAs, including those involved in cell development and photosynthesis, jasmonate signaling, and disease resistance response. The application of miRVIT in a biological context confirmed the effectiveness of the followed approach, especially for the identification of novel miRNA candidates in grapevine. miRVIT database is available at http://mirvit.ipsp.cnr.it. Highlights: The application of the newly produced database of grapevine novel miRNAs to the analysis of plants infected by Flavescence dorée reveals key roles of miRNAs in photosynthesis and jasmonate signaling.

Spatio-temporal expression of miRNA159 family members and their GAMYB target gene during the modulation of gibberellin-induced grapevine parthenocarpy

Grapevine, Vitis vinifera, is an important economic fruit crop that is highly sensitive to gibberellin (GA), and the exogenous application of GA can efficiently induce grapevine parthenocarpy. However, the molecular mechanisms underlying this process remain elusive. In this study, morphological changes during flower development in response to GA treatments were examined in the 'Zuijinxiang' cultivar. To obtain insights into the roles of miRNA159s in GA-induced grapevine parthenocarpy, VvmiR159a, VvmiR159b, VvmiR159c, and their target gene VvGAMYB were isolated, sequenced and characterized. Spatial-temporal expression analyses showed that VvmiR159c exhibited the highest expression levels at 4 d before flowering, followed by a gradual decrease, while VvGAMYB displayed an opposite pattern of expression with the lowest expression at the corresponding stage in response to GA treatment. A cleavage interaction between VvmiR159s and VvGAMYB and variations of their cleavage roles were confirmed in grapevine floral development. In addition, the potential roles of VvmiR159s in GA signaling were investigated through DELLA-protein repressors, indicating that GA-DELLA (SLR1)-VvmiR159c-VvGAMYB is the key signaling regulatory module in grapevine. Our findings provide novel insights into the GA-responsive roles of VvmiR159s in modulating grapevine floral development, which have important implications for the molecular breeding of high-quality seedless grapevine berry.

Comparative analysis of miRNA and mRNA abundance in determinate cucumber by high-throughput sequencing

Determinate cucumber is a type of short vines, fewer nodes, and terminal flowers, it is suitable for high-density planting and available harvesting in field cultivation, whereas the indeterminate cucumber is preferred to cultivate in greenhouses. However, many biotic or abiotic stresses could lead indeterminate cucumber to be determinate in greenhouse cultivation, which may decrease yield and fruit quality. Therefore, it is urgent and essential to investigate the key factors forming determinate and terminal flowering in cucumber. In this study, two close background inbred lines were selected and conducted the miRNA and mRNA high throughput sequencing. Interestingly, ethylene-associated miRNAs and mRNAs were intensively obtained, indicating that the plant hormone ethylene is a key factor impacting determinate and terminal flowering in cucumber. The ethylene metabolites analysis showed that significant higher ethylene was observed in determinate line than that in the indeterminate line. The RT-qPCR validation of ethylene related miRNAs Cas-miR172, Cas-miR396, and Cas-miR414 and their target mRNAs showed a significant negative correlation. These data suggested that ethylene-associated miRNAs might affect determinate and terminal flower phenotypes by regulating their target genes expression. This study not only provides a potential molecular mechanism for determinate formation in cucumber but also establishes a method to demonstrate important physiological processes through the comprehensive association of miRNA and mRNA high-throughput sequencing.

Cold-responsive miRNAs and their target genes in the wild eggplant species Solanum aculeatissimum

Background

Low temperature is an important abiotic stress in plant growth and development, especially for thermophilic plants. Eggplants are thermophilic vegetables, although the molecular mechanism of their response to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although the role of many plant miRNAs in facilitating chilling tolerance has been verified, little is known about the mechanisms of eggplant chilling tolerance.

Results

Here, we used high-throughput sequencing to extract the miRNA and target genes expression profiles of Solanum aculeatissimum (S. aculeatissimum) under low temperature stress at different time periods(0 h, 2 h, 6 h, 12 h, 24 h). Differentially regulated miRNAs and their target genes were analyzed by comparing the small RNA (sRNA) and miRBase 20.0 databases using BLAST or BOWTIE, respectively. Fifty-six down-regulated miRNAs and 28 up-regulated miRNAs corresponding to 220 up-regulated mRNAs and 94 down-regulated mRNAs, respectively, were identified in S. aculeatissimum. Nine significant differentially expressed miRNAs and twelve mRNAs were identified by quantitative Real-time PCR and association analysis, and analyzed for their GO function enrichment and KEGG pathway association.

Conclusions

In summary, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, which provides a theoretical basis for the further study of low temperature stress-related miRNAs and the regulation of cold-tolerance mechanisms of eggplant at the miRNA level.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4341-y) contains supplementary material, which is available to authorized users.

microRNAs and Their Targets in Apple (Malus domestica cv. "Fuji") Involved in Response to Infection of Pathogen Valsa mali

miRNAs are important regulators involving in plant-pathogen interactions. However, their roles in apple tree response to Valsa canker pathogen (Valsa mali, Vm) infection were poorly understood. In this study, we constructed two miRNA libraries using the twig bark tissues of apple tree (Malus domestica Borkh. cv. "Fuji") inoculated with Vm (IVm) and PDA medium (control, BMd). Among all detected miRNAs, 23 miRNAs were specifically isolated from BMd and 39 miRNAs were specifically isolated from IVm. Meanwhile, the expression of 294 miRNAs decreased; and another 172 miRNAs showed an increased expression trend in IVm compared with that in BMd. Furthermore, two degradome sequencing libraries were also constructed to identify the target genes of these miRNAs. In total, 353 differentially expressed miRNAs between IVm and BMd were detected to be able to target 1,077 unigenes with 2,251 cleavage sites. Based on GO and KEGG analysis, these genes were found to be mainly related to transcription regulation and signal transduction. In addition, we selected 17 miRNAs and 22 corresponding target genes to screen the expression profiles when apple twigs were infected by Vm. The expression trends of most miRNAs/target genes were consist with the results of deep sequencing. Many of them may involve in the apple twig-Vm interaction by inducing/reducing their expression. What's more, miRNAs and their target genes regulate the apple twig-Vm interaction by forming many complicated regulation networks rather than one to one model. It is worth that a conserved miRNAs mdm-miR482b, which was down regulated in IVm compared with BMd, has 14 potential target genes, most of which are disease resistance related genes. This indicates that mdm-miR482b may play important roles in apple twig response to Vm. More important, the feedback regulation of sRNA pathway in apple twig is also very complex, and play critical role in the interaction between apple twig and Vm based on the results of expression analysis. In all, the results will provide insights into the crucial functions of miRNAs in the woody plant, apple tree-Vm interaction.

Ectopic expression of CSD1 and CSD2 targeting genes of miR398 in grapevine is associated with oxidative stress tolerance

MicroRNAs (miRNAs) are endogenous small RNAs of -21 nucleotides that play an important role in diverse plant physiological processes at the post-transcriptional level by directing mRNA cleavage or translational inhibition. Previous studies have indicated that down-regulation of miR398 in response to oxidative stress allows up-regulation of the two target genes, cytosolic CSD1 and chloroplastic CSD2 (copper/zinc superoxide dismutase), resulting in protecting the plants to tolerate oxidative stress. In this study, we provide evidence that grapevine miR398 (Vv-miR398), by regulating the expression of its target genes, VvCSD1 and VvCSD2, mediates responses of grapevine to copper (Cu) stress which have been magnified due to increase in Cu-containing pesticide application. The expression of Vv-miR398 was inhibited by different concentrations of Cu stress; on the other hand, there was a steady increase in the activity of VvCSD1 and VvCSD2 genes. The function of VvCSD1 and VvCSD2 under Cu stress was thoroughly examined by overexpressing the use of the VvCSD1 and VvCSD2 in transgenic tobacco (Nicotiana tabacum). We found that both the overexpressed transgenic lines had lower Cu sensitivity and higher Cu tolerance compared with the wild type. In addition, lower levels of ROS and higher levels of SOD activities were accumulated in the transgenic lines in comparison with the wild type under the higher Cu conditions. Furthermore, these transgenic tobacco lines also recorded a higher UV and salt tolerance than the WT plants. These results suggested that overexpressing the VvCSDs will enhance the ROS-scavenging systems and protect the plant against more oxidative damage. Also, more investigations in this line are needed that would provide significant improvements in our understanding the resistance of fruit crops to environmental stress.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Characterization of miR061 and its target genes in grapevine responding to exogenous gibberellic acid

MicroRNAs (miRNAs), as an important growth regulator, are also involved in gibberellic acid (GA) signaling, revealing much relationship between miRNAs and GA in various plant responses. Grape is highly sensitive to GA₃, which plays a significant regulatory role in regulation of flower development, berry expansion, berry set, berry ripening, and seedlessness induction; further, it was found that grapevine miR061 (VvmiR061) is a GA₃ responsive miRNA. In this study, grapevine REV (VvREV) and HOX32 (VvHOX32), two target genes of VvmiR061, were predicted, verified, and cloned; homologous conservation was analyzed in various plants. The expression profiles of both VvmiR061 and its target genes (VvREV and VvHOX32) under GA₃ treatment were detected by qRT-PCR during grapevine flower and berry development. Results revealed that GA₃ treatment has upregulated the transcription of VvREV and VvHOX32, while it downregulated the expression of VvmiR061. The function of VvmiR061 in cleaving target genes VvREV and VvHOX32 was diminished by GA₃ treatment during flower developmental process. The results of this study exhibited the importance of VvmiR061 in regulating flower development and GA₃ signaling pathway and also contributed some to the knowledge of small RNA-mediated regulation in grape.

Combined analysis of mRNA and miRNA identifies dehydration and salinity responsive key molecular players in citrus roots

Citrus is one of the most economically important fruit crops around world. Drought and salinity stresses adversely affected its productivity and fruit quality. However, the genetic regulatory networks and signaling pathways involved in drought and salinity remain to be elucidated. With RNA-seq and sRNA-seq, an integrative analysis of miRNA and mRNA expression profiling and their regulatory networks were conducted using citrus roots subjected to dehydration and salt treatment. Differentially expressed (DE) mRNA and miRNA profiles were obtained according to fold change analysis and the relationships between miRNAs and target mRNAs were found to be coherent and incoherent in the regulatory networks. GO enrichment analysis revealed that some crucial biological processes related to signal transduction (e.g. ‘MAPK cascade’), hormone-mediated signaling pathways (e.g. abscisic acid- activated signaling pathway’), reactive oxygen species (ROS) metabolic process (e.g. ‘hydrogen peroxide catabolic process’) and transcription factors (e.g., ‘MYB, ZFP and bZIP’) were involved in dehydration and/or salt treatment. The molecular players in response to dehydration and salt treatment were partially overlapping. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-seq and sRNA-seq analysis. This study provides new insights into the molecular mechanisms how citrus roots respond to dehydration and salt treatment.

Small RNA-Sequencing Links Physiological Changes and RdDM Process to Vegetative-to-Floral Transition in Apple

Transition from vegetative to floral buds is a critical physiological change during flower induction that determines fruit productivity. Small non-coding RNAs (sRNAs) including microRNAs (miRNAs) and small interfering RNAs (siRNAs) are pivotal regulators of plant growth and development. Although the key role of sRNAs in flowering regulation has been well-described in Arabidopsis and some other annual plants, their relevance to vegetative-to-floral transition (hereafter, referred to floral transition) in perennial woody trees remains under defined. Here, we performed Illumina sequencing of sRNA libraries prepared from vegetative and floral bud during flower induction of the apple trees. A large number of sRNAs exemplified by 33 previously annotated miRNAs and six novel members display significant differential expression (DE) patterns. Notably, most of these DE-miRNAs in floral transition displayed opposite expression changes in reported phase transition in apple trees. Bioinformatics analysis suggests most of the DE-miRNAs targeted transcripts involved in SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) gene regulation, stress responses, and auxin and gibberellin (GA) pathways, with further suggestion that there is an inherent link between physiological stress response and metabolism reprogramming during floral transition. We also observed significant changes in 24 nucleotide (nt) sRNAs that are hallmarks for RNA-dependent DNA methylation (RdDM) pathway, suggestive of the correlation between epigenetic modifications and the floral transition. The study not only provides new insight into our understanding of fundamental mechanism of poorly studied floral transition in apple and other woody plants, but also presents important sRNA resource for future in-depth research in the apple flowering physiology.

Differential expression of miRNAs and associated gene targets in grapevine leafroll-associated virus 3-infected plants

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs (sRNA) that play an essential role in the regulation of target mRNAs expressed during plant development and in response to stress. MicroRNA expression profiling has helped to identify miRNAs that regulate a range of processes, including the plant's defence response to pathogens. In this study, differential miRNA expression in own-rooted Vitis vinifera cv. Cabernet Sauvignon plants infected with grapevine leafroll-associated virus 3 was investigated with microarrays and next-generation sequencing (NGS) of sRNA and mRNA. These high-throughput approaches identified several differentially expressed miRNAs. Four miRNAs, identified by both approaches, were validated by stemloop RT-PCRs. Three of the predicted targets of the differentially expressed miRNAs were also differentially expressed in the transcriptome data of infected plants, and were validated by RT-qPCR. Identification of these miRNAs and their targets can lead to a better understanding of host-pathogen interactions involved in grapevine leafroll disease and the identification of possible targets for virus resistance.

Small RNA and degradome profiling reveals important roles for microRNAs and their targets in tea plant response to drought stress

Tea (Camellia sinensis) is a popular beverage worldwide. Drought stress (DS) is a major constraint on the growth, yield and quality of tea plants. MicroRNAs (miRNAs) play important roles in plant responses to DS. We constructed eight small RNA libraries from the drought-tolerant 'Ningzhou 2' (NZ2) and drought-susceptible 'Zhuyeqi' (ZYQ) cultivars during four stages [control (CK), the fourth day of DS, the eighth day of DS and after recovery (RC)]. A total of 268 conserved and 62 novel miRNAs were identified using small RNA sequencing. In total, 139 (52.9%) and 96 (36.0%) conserved miRNAs were differentially expressed during the four stages (P ≤ 0.05) in NZ2 and ZYQ, respectively. A total of 814 predicted target genes were identified as differentially regulated by 199 miRNAs through degradome sequencing. Among them, 201 and 218 genes were specific to the NZ2 and ZYQ cultivars, respectively, and 395 were common to both cultivars. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed the biological roles of these targets and showed that some of the targets responded to DS in a stress- and cultivar-dependent manner. Correlated expression patterns between miRNA and their targets showed that specific miRNAs target the miRNA effector Argonaute 1 (AGO1), drought signaling-related receptors and enzymes, transcription factors, and other structural and functional proteins. The predicted regulatory networks provide insights into a potential miRNA-mediated regulatory mechanism. These results will contribute to the breeding of drought-tolerant tea plants and to elucidating miRNA regulation in response to drought.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Study on Expression Modes and Cleavage Role of miR156b/c/d and its Target Gene Vv-SPL9 During the Whole Growth Stage of Grapevine

miR156 regulates the expression of its target SPL (PROMOTER BINDING-LIKE) genes during flower and fruit development, diverse developmental stage transitions, especially from vegetative to reproductive growth phases, by cleaving the target mRNA SPL of one plant-specific transcription factor. However, systematic reports on grapevine have yet to be presented. Here, the precise sequence of miR156 (vvi-miR156b/c/d) in grapevine "Takatsuma" was cloned with a previously cloned grapevine SPL (Vv-SPL9). Expression profiles in 18 grapevine tissues were identified through stem-loop RT-PCR. The interaction mode between vvi-miR156b/c/d and Vv-SPL9 was further validated by detecting the cleavage site and cleavage products of 3'- and 5'-ends via an integrated approach of 5'-RLM-RACE (RNA ligase-mediated 5'-rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3'-rapid amplification of cDNA ends), and qRT-PCR (real time reverse transcriptase-polymerase chain reaction). The variation in their cleavage roles in the whole growth stage of grapevine was also systematically investigated. Results showed that vvi-miR156b/c/d exhibited typical temporal-spatial-specific expression levels. The expression levels were higher in vegetative organs, such as leaf, than in reproductive organs, such as tendrils, flowers, and berries. A significant variation was observed during vegetative-to-reproductive transition. The expression patterns of Vv-SPL9 showed the opposite trends with those of vvi-miR156b. We confirmed that the cleavage site was at the 10th site of vvi-miR156b/c/d complementary to Vv-SPL9 in "Takatsuma" grapevine. We also identified the temporal-spatial variation of the cleavage products. This variation can indicate the regulatory function of miR156 on SPL in grapevines. Our findings provide further insights into the functions of vvi-miR156b/c/d and its target Vv-SPL9, and also help enrich our knowledge of small RNA-mediated regulation in grapevine.

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

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

Genomic profiling of exogenous abscisic acid-responsive microRNAs in tomato (Solanum lycopersicum)

BACKGROUND

Plant microRNAs (miRNAs) are involved in various biological pathways and stress responses as negative regulators at the posttranscriptional level. Abscisic acid (ABA) is a key signaling molecule that mediates plant stress response by activating many stress-related genes. Although some miRNAs in plants are previously identified to respond to ABA, a comprehensive profile of ABA-responsive miRNAs has not yet been elucidated.

RESULTS

Here, we identified miRNAs responding to exogenous application of ABA, and their predicted target genes in the model plant organism tomato (Solanum lycopersicum). Deep sequencing of small RNAs from ABA-treated and untreated tomatoes revealed that miRNAs can be up- or down-regulated upon treatment with ABA. A total of 1067 miRNAs were detected (including 365 known and 702 candidate novel miRNAs), of those, 416 miRNAs which had an abundance over two TPM (transcripts per million) were selected for differential expression analysis. We identified 269 (180 known and 89 novel) miRNAs that respond to exogenous ABA treatment with a change in expression level of |log2FC|≥0.25. 136 of these miRNAs (90 known and 46 novel) were expressed at significantly different levels |log2FC|≥1 between treatments. Furthermore, stem-loop RT-PCR was applied to validate the RNA-seq data. Target prediction and analysis of the corresponding ABA-responsive transcriptome data uncovered that differentially expressed miRNAs are involved in condition stress and pathogen resistance, growth and development. Among them, approximately 90 miRNAs were predicted to target transcription factors and pathogen resistance genes. Some miRNAs had functional overlap in biotic and abiotic stress. Most of these miRNAs were down-regulated following exposure to exogenous ABA, while their related target genes were inversely up-regulated, which is consistent with their negative regulatory role in gene expression.

CONCLUSIONS

Exogenous ABA application influences the composition and expression level of tomato miRNAs. ABA mainly down-regulates miRNAs that their target genes involve in abiotic stress adaption and disease resistance. ABA might increase expression of stress-related genes via miRNA-mediated posttranscriptional regulation, and our results indicate that ABA treatment has the potential to improve both abiotic stress tolerance and pathogen resistance. This study presents a comprehensive profile of ABA-regulated miRNAs in the tomato, and provides a robust database for further investigation of ABA regulatory mechanisms.

New insights of medicinal plant therapeutic activity-The miRNA transfer

MicroRNA (miRNA) has become the spotlight of the biomedical research around the world and is considered to be a major post-transcriptional gene regulator. This small, endogenous RNA (21-25 nucleotides long) plays an important role by targeting specific mRNAs in plants, animals and humans. Herbal medicine has been used for thousands of years, however little is known about its molecular mechanism of action. Since the discovery of plant miRNA in human tissue and sera after ingestion, the connection between the two kingdoms is presented under a new perspective. Forward pharmacology, such as miRNA therapeutics could be the next best step toward identifying novel therapeutic options involving medicinal plants. Besides reporting the latest findings regarding the cross-kingdom transfer of miRNA and its therapeutic application, this review can inform further investigations that could lead to a modern definition of herbal medicine.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Transcriptome-Wide Identification of miRNAs and Their Targets from Typha angustifolia by RNA-Seq and Their Response to Cadmium Stress

MicroRNAs (miRNAs) play important roles in plant responses to environmental stress. In this work, we used high-throughput sequencing to analyze transcriptome and small RNAs (sRNAs) in Typha angustifolia under cadmium (Cd) stress. 57,608,230 raw reads were obtained from deep sequencing of a pooled cDNA library. Sequence assembly and analysis yielded 102,473 unigenes. We subsequently sequenced two sRNA libraries from T. angustifolia with or without Cd exposure respectively. Based on transcriptome data of T. angustifolia, we catalogued and analyzed the sRNAs, resulting in the identification of 114 conserved miRNAs and 41 novel candidate miRNAs in both small RNA libraries. In silico analysis revealed 764 targets for 89 conserved miRNAs and 21 novel miRNAs. Statistical analysis on sequencing reads abundance and experimental validation revealed that 4 conserved and 6 novel miRNAs showed specific expression. Combined with function of target genes, these results suggested that miRNAs might play a role in plant Cd stress response. This study provided the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in T. angustifolia, which provide a framework for further analysis of miRNAs and their role in regulating plant responses to Cd stress.

Integrated RNA-Seq and sRNA-Seq Analysis Identifies Chilling and Freezing Responsive Key Molecular Players and Pathways in Tea Plant (Camellia sinensis)

Tea [Camellia sinensis (L) O. Kuntze, Theaceae] is one of the most popular non-alcoholic beverages worldwide. Cold stress is one of the most severe abiotic stresses that limit tea plants’ growth, survival and geographical distribution. However, the genetic regulatory network and signaling pathways involved in cold stress responses in tea plants remain unearthed. Using RNA-Seq, DGE and sRNA-Seq technologies, we performed an integrative analysis of miRNA and mRNA expression profiling and their regulatory network of tea plants under chilling (4℃) and freezing (-5℃) stress. Differentially expressed (DE) miRNA and mRNA profiles were obtained based on fold change analysis, miRNAs and target mRNAs were found to show both coherent and incoherent relationships in the regulatory network. Furthermore, we compared several key pathways (e.g., ‘Photosynthesis’), GO terms (e.g., ‘response to karrikin’) and transcriptional factors (TFs, e.g., DREB1b/CBF1) which were identified as involved in the early chilling and/or freezing response of tea plants. Intriguingly, we found that karrikins, a new group of plant growth regulators, and β-primeverosidase (BPR), a key enzyme functionally relevant with the formation of tea aroma might play an important role in both early chilling and freezing response of tea plants. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-Seq and sRNA-Seq analysis. This is the first study to simultaneously profile the expression patterns of both miRNAs and mRNAs on a genome-wide scale to elucidate the molecular mechanisms of early responses of tea plants to cold stress. In addition to gaining a deeper insight into the cold resistant characteristics of tea plants, we provide a good case study to analyse mRNA/miRNA expression and profiling of non-model plant species using next-generation sequencing technology.

Identification and characterization of microRNAs from tree peony (Paeonia ostii) and their response to copper stress

MicroRNAs (miRNAs) are a class of non-coding, small RNAs recognized as important regulators of gene expression. Although plant miRNAs have been extensively studied in model systems, less is known in other plants with limited genome sequence data, including Paeonia ostii. In this work, we used high-throughput sequencing to identify conserved and nonconserved miRNAs and other short RNAs in Paeonia ostii under control and copper stressed condition. 102 previously known plant miRNAs were identified and classified into 89 families according to their gene sequence identity. Some miRNAs were highly conserved in the plant kingdom suggesting that these miRNA play important and conserved roles. Combined our transcriptome sequencing data of Paeonia ostii under same conditions, 34 novel potential miRNAs were identified. The potential targets of the identified known and novel miRNAs were also predicted based on sequence homology search. Comparing the two libraries, it was observed that 12 conserved miRNAs and 18 novel miRNAs showed significantly changes in response to copper stress. Some of the new identified potential miRNAs might be involved in Paeonia ostii-specific regulating mechanisms under copper stress. These results provide a framework for further analysis of miRNAs and their role in regulating Paeonia ostii response to copper stress.

Identification and characterization of microRNAs from Chinese pollination constant non-astringent persimmon using high-throughput sequencing

BACKGROUND

microRNAs (miRNAs) have been shown to play key roles in regulating gene expression at post-transcriptional level, but miRNAs associated with natural deastringency of Chinese pollination-constant nonastringent persimmon (CPCNA) have never been identified.

RESULTS

In this study, two small RNA libraries established using 'Eshi No. 1' persimmon (Diospyros kaki Thunb.; CPCNA) fruits collected at 15 and 20 weeks after flowering (WAF) were sequenced through Solexa platform in order to identify miRNAs involved in deastringency of persimmon. A total of 6,258,487 and 7,634,169 reads were generated for the libraries at 15 and 20 WAF, respectively. Based on sequence similarity and hairpin structure prediction, 236 known miRNAs belonging to 65 miRNA families and 33 novel miRNAs were identified using persimmon transcriptome data. Sixty one of the characterized miRNAs exhibited pronounced difference in the expression levels between 15 and 20 WAF, 17 up-regulated and 44 down-regulated. Expression profiles of 12 conserved and 10 novel miRNAs were validated by stem loop qRT-PCR. A total of 198 target genes were predicted for the differentially expressed miRNAs, including several genes that have been reported to be implicated in proanthocyanidins (PAs, or called tannin) accumulation. In addition, two transcription factors, a GRF and a bHLH, were experimentally confirmed as the targets of dka-miR396 and dka-miR395, respectively.

CONCLUSIONS

Taken together, the present data unraveled several important miRNAs in persimmon. Among them, miR395p-3p and miR858b may regulate bHLH and MYB, respectively, which are influenced by SPL under the control of miR156j-5p and in turn regulate the structural genes involved in PA biosynthesis. In addition, dka-miR396g and miR2911a may regulate their target genes associated with glucosylation and insolubilization of tannin precursors. All of these miRNAs might play key roles in the regulation of (de)astringency in persimmon fruits under normal development conditions.

Identification of cold-inducible microRNAs in grapevine

Low temperature is one of the most important environmental factors that limits the geographical distribution and productivity of grapevine. However, the molecular mechanisms on how grapevine responds to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although miRNAs and their targets have been identified in several Vitis species, their participation during cold accumulation in grapevine remains unknown. In this study, two small RNA libraries were generated from micropropagated 'Muscat Hamburg' (V. vinifera) plantlets under normal and low temperatures (4°C). A total of 163 known miRNAs and 67 putative novel miRNAs were detected from two small RNA libraries by Solexa sequencing. Forty-four cold-inducible miRNAs were identified through differentially expressed miRNAs (DEMs) analysis; among which, 13 belonged to upregulated DEMs while 31 belonged downregulated DEMs. The expression patterns of the 13 DEMs were verified by real-time RT-PCR analysis. The prediction of the target genes for DEMs indicated that miRNA may regulate transcription factors, including AP2, SBP, MYB, bHLH, GRAS, and bZIP under cold stress. The 5'-RLM RACE were conducted to verify the cleavage site of predicted targets. Seven predicted target genes for four known and three novel vvi-miRNAs showed specific cleavage sites corresponding to their miRNA complementary sequences. The expression pattern of these seven target genes revealed negative correlation with the expression level of the corresponding vvi-miRNAs. Our results indicated that a diverse set of miRNAs in V. vinifera are cold-inducible and may play an important role in cold stress response.

Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play versatile roles in plants, functioning in processes such as growth, development and stress responses. Chilling is a common abiotic stress that seriously affects plants growth and development. Recently, chilling-responsive miRNAs have been detected in several plant species. However, little is known about the miRNAs in the model plant tomato. 'LA1777' (Solanum habrochaites) has been shown to survive chilling stress due to its various characteristics.

RESULTS

Here, two small RNA libraries and two degradome libraries were produced from chilling-treated (CT) and non-chilling-treated (NT) leaves of S. habrochaites seedlings. Following high-throughput sequencing and filtering, 161 conserved and 236 novel miRNAs were identified in the two libraries. Of these miRNAs, 192 increased in the response to chilling stress while 205 decreased. Furthermore, the target genes of the miRNAs were predicted using a degradome sequencing approach. It was found that 62 target genes were cleaved by 42 conserved miRNAs, while nine target genes were cleaved by nine novel miRNAs. Additionally, nine miRNAs and six target genes were validated by quantitative real-time PCR (qRT-PCR). Target gene functional analysis showed that most target genes played positive roles in the chilling response, primarily by regulating the expression of anti-stress proteins, antioxidant enzyme and genes involved in cell wall formation.

CONCLUSIONS

Tomato is an important model plant for basic biological research. In this study, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, and the target genes were analyzed by degradome sequencing. The work helps identify chilling-responsive miRNAs in tomato and increases the number of identified miRNAs involved in chilling stress. Furthermore, the work provides a foundation for further study of the regulation of miRNAs in the plant response to chilling stress.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Transient expression assays in grapevine: a step towards genetic improvement

In the past few years, the usefulness of transient expression assays has continuously increased for the characterization of unknown gene function and metabolic pathways. In grapevine (Vitis vinifera L.), one of the most economically important fruit crops in the world, recent systematic sequencing projects produced many gene data sets that require detailed analysis. Due to their rapid nature, transient expression assays are well suited for large-scale genetic studies. Although genes and metabolic pathways of any species can be analysed by transient expression in model plants, a need for homologous systems has emerged to avoid the misinterpretation of results due to a foreign genetic background. Over the last 10 years, various protocols have thus been developed to apply this powerful technology to grapevine. Using cell suspension cultures, somatic embryos, leaves or whole plantlets, transient expression assays enabled the study of the function, regulation and subcellular localization of genes involved in specific metabolic pathways such as the biosynthesis of phenylpropanoids. Disease resistance genes that could be used for marker-assisted selection in conventional breeding or for stable transformation of elite cultivars have also been characterized. Additionally, transient expression assays have proved useful for shaping new tools for grapevine genetic improvement: synthetic promoters, silencing constructs, minimal linear cassettes or viral vectors. This review provides an update on the different tools (DNA constructs, reporter genes, vectors) and methods (Agrobacterium-mediated and direct gene transfer methods) available for transient gene expression in grapevine. The most representative results published thus far are then described.

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.

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

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

Identification and comparative analysis of the Tegillarca granosa haemocytes microRNA transcriptome in response to Cd using a deep sequencing approach

Background

MicroRNAs (miRNAs) are endogenous non-coding small RNAs (sRNAs) that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. To identify miRNAs regulated by heavy metal stress, we constructed two sRNA libraries for the blood clam Tegillarca granosa: one for organisms exposed to toxic levels of cadmium (Cd) and one for a control group.

Results

Sequencing of the two libraries and subsequent analysis revealed 215 conserved and 39 new miRNAs. Most of the new miRNAs in T. granosa were up- or down-regulated in response to Cd exposure. There were significant differences in expression between the Cd and control groups for 16 miRNAs. Of these, five miRNAs were significantly up-regulated and 11 were significantly down-regulated in the Cd stress library. Potential targets were predicted for the 16 differential miRNAs in pre-miRNAs identified according to sequence homology. Some of the predicted miRNA targets are associated with regulation of the response to stress induced by heavy metals. Five differentially expressed miRNAs (Tgr-nmiR-8, Tgr-nmiR-21, Tgr-miR-2a, Tgr-miR-10a-5p, and Tgr-miR-184b) were validated by qRT-PCR.

Conclusion

Our study is the first large-scale identification of miRNAs in T. granosa haemocytes. Our findings suggest that some miRNAs and their target genes and pathways may play critical roles in the responses of this species to environmental heavy metal stresses.

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.

Grapevine microRNAs responsive to exogenous gibberellin

Background

MicroRNAs (miRNAs), involving in various biological and metabolic processes, have been discovered and analyzed in quite a number of plants species, such as Arabidopsis, rice and other plants. However, there have been few reports about grapevine miRNAs in response to gibberelline (GA₃).

Results

Solexa technology was used to sequence small RNA libraries constructed from grapevine berries treated with GA₃ and the control. A total of 122 known and 90 novel grapevine miRNAs (Vvi-miRNAs) were identified. Totally, 137 ones were found to be clearly responsive to GA₃, among which 58 were down-regulated, 51 were up-regulated, 21 could only be detected in the control, and seven were only detected in the treatment. Subsequently, we found that 28 of them were differentially regulated by GA₃, with 12 conserved and 16 novel Vvi-miRNAs, based on the analysis of qRT-PCR essays. There existed some consistency in expression levels of GA₃-responsive Vvi-miRNAs between high throughput sequencing and qRT-PCR essays. In addition, 117 target genes for 29 novel miRNAs were predicted.

Conclusions

Deep sequencing of short RNAs from grapevine berries treated with GA3 and the control identified 137 GA₃-responsive miRNAs, among which 28 exhibited different expression profiles of response to GA₃ in the diverse developmental stages of grapevine berries. These identified Vvi-miRNAs might be involved in the grapevine berry development and response to environmental stresses.

Transcriptome-wide analysis of dynamic variations in regulation modes of grapevine microRNAs on their target genes during grapevine development

MicroRNAs (miRNAs) play critical regulatory roles mainly through cleaving their target mRNAs or repressing gene translation during plant development. Grapevines are among the most economically important fruit crops with available whole genome sequences. Studies on grapevine miRNAs (Vv-miRNAs) are also widely available. However, studies on the regulation mode of Vv-miRNAs on their target mRNAs during grapevine development have not been studied well, especially at the transcriptome-wide level. Here, six small RNA and mRNA libraries from various grapevine tissues were constructed for Illumina and Degradome sequencing. Subsequently, we systematically analyzed the spatiotemporal variations in the regulation of the target genes of regulation of Vv-miRNAs. In total, 242 known and 132 novel Vv-miRNAs and 193 target mRNAs were identified, including 103 target mRNAs for known and 90 target mRNAs for novel miRNAs, were validated in one or more of the tissues examined. More than 50 % of novel miRNAs were expressed exclusively in the flowers and berries, where they cleaved their target genes in a tissue-specific manner, especially, the breadth of their cleavage sites in flower tissues. Moreover, six novel miRNAs in berries responded to exogenous gibberellin and/or ethylene under a quantitative real time RT-PCR analysis, which confirmed their regulatory functions during berry development. Up to 93.6 % of the known miRNAs were highly conserved in various tissues, where their expression levels exhibited dynamic variations during grapevine development. Significantly, some Vv-miRNA families had one key member that acted as the main regulator of their target genes during grapevine development.

Cloning and characterization of miRNAs and their targets, including a novel miRNA-targeted NBS-LRR protein class gene in apple (Golden Delicious)

MicroRNA (miRNA) has emerged as an important regulator of gene expression in plants. 146 miRNAs were identified from apple (Malus domestica cv. Golden Delicious) by bioinformatic analysis and RNA library sequencing. From these, 135 were conserved and 11 were novel miRNAs. Target analysis predicted one of the novel miRNAs, Md-miRLn11 (Malus domestica microRNA Ln11), targeted an apple nucleotide-binding site (NBS)-leucine-rich repeat (LRR) class protein coding gene (Md-NBS). 5' RACE assay confirmed the ability of Md-miRLn11 to cleave Md-NBS at the 11-12-nt position. Analysis of the expression of Md-miRLn11 and Md-NBS during the optimum invasion period in 40 apple varieties showed that the expression of Md-NBS gene in resistant varieties is higher than in susceptible varieties, with an inverse pattern for Md-miRLn11. Seedlings from the resistant apple variety 'JiGuan' were used to carry out an Agrobacterium infiltration assay, and then inoculated with the apple leaf spot disease. The result showed a clear decline of disease resistance in JiGuan apples. In contrast, the susceptible variety 'FuJi' infiltrated with the Md-NBS gene showed a significant increase in disease resistance. Based on the above results, we propose that Md-miRLn11 regulates Md-NBS gene expression in particular under the condition of pathogen infection, and that the Md-miRLn11 targeting P-loop site may regulate many NBS-LRR protein class genes in woody plants.

Identification of potential miRNAs and their targets in Vriesea carinata (Poales, Bromeliaceae)

The miRNAs play important roles in regulation of gene expression at the post-transcriptional level. A small RNA and RNA-seq of libraries were constructed to identify miRNAs in Vriesea carinata, a native bromeliad species from Brazilian Atlantic Rainforest. Illumina technology was used to perform high throughput sequencing and data was analyzed using bioinformatics tools. We obtained 2,191,509 mature miRNAs sequences representing 54 conserved families in plant species. Further analysis allowed the prediction of secondary structures for 19 conserved and 16 novel miRNAs. Potential targets were predicted from pre-miRNAs by sequence homology and validated using RTqPCR approach. This study provides the first identification of miRNAs and their potential targets of a bromeliad species.

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

In plants, microRNAs (miRNAs) play a critical role in post-transcriptional gene regulation and have been shown to control many genes involved in various biological and metabolic processes. There have been extensive studies to discover miRNAs and analyze their functions in model plant species, such as Arabidopsis and rice. Deep sequencing technologies have facilitated identification of species-specific or lowly expressed as well as conserved or highly expressed miRNAs in plants. In this research, we used Solexa sequencing to discover new miRNAs in cultivated strawberry (Fragaria×ananassa). A total of 23,282 ,309 reads representing 22,500 ,402 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from strawberry fruit tissues. On the basis of sequence similarity and hairpin structure prediction, we found that 156,639 reads representing 153 sequences have good matches to known miRNAs. We also identified 37 novel miRNA candidates. These sequences had not been previously described in other plant species. Potential target genes were predicted for the majority of and novel miRNAs. These results show that regulatory miRNAs exist in the agriculturally important cultivated strawberry and may play an important role in its growth, development and response to disease.

Characterization of target mRNAs for grapevine microRNAs with an integrated strategy of modified RLM-RACE, newly developed PPM-RACE and qPCRs

MicroRNAs (miRNAs) regulate target gene expression by mediating target gene cleavage or inhibition of translation at transcriptional and post-transcriptional levels in higher plants. Until now, many grapevine microRNAs (Vv-miRNAs) have been identified and quite a number of miRNA target genes were also verified by various analysis. However, global interaction of miRNAs with their target genes still remained to perform more research. We reported experimental validation of a number of miRNA target genes in table grapevine that had been previously identified by bioinformatics in our earlier studies. To verify more predicted target genes of Vv-miRNAs and elucidate the modes by which these Vv-miRNAs work on their target genes, 31 unverified potential target genes for 18 Vv-miRNAs were experimentally verified by a new integrated strategy employing a modified 5'-RLM-RACE (RNA ligase-mediated 5' rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3' rapid amplification of cDNA ends) and qRT-PCRs of cleavage products. The results showed that these Vv-miRNAs negatively regulated expression of their target messenger RNAs (mRNAs) through guiding corresponding target mRNA cleavage, of which about 94.4% Vv-miRNAs cleaved their target mRNAs mainly at the tenth nucleotide of 5'-end of miRNAs. Expression levels of both miRNAs and their target mRNAs in eight tissues exhibited inverse relationships, and expressions both of cleaved targets and miRNAs indicated a cleavage mode of Vv-miRNAs on their target genes. Our results confirm the importance of Vv-miRNAs in grapevine growth and development, and suggest more study on Vv-miRNAs and targets can enrich the knowledge of miRNA mediated-regulation in grapevine.