Skip to main content

Grapevine microRNAs responsive to exogenous gibberellin

Abstract

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 (GA3).

Results

Solexa technology was used to sequence small RNA libraries constructed from grapevine berries treated with GA3 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 GA3, 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 GA3, with 12 conserved and 16 novel Vvi-miRNAs, based on the analysis of qRT-PCR essays. There existed some consistency in expression levels of GA3-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 GA3-responsive miRNAs, among which 28 exhibited different expression profiles of response to GA3 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.

Background

MicroRNAs (miRNAs) are endogenous gene regulators distributed widely in plant genomes, and they play important roles in plant growth, development, signal transduction and response to environmental stimuli [15]. Identification of miRNAs is a key step in gaining insight into sRNA-based regulatory functions with many conserved miRNAs having been identified by traditional sequencing approaches such as the Sanger sequencing method [6]. However, most species-specific or tissue-specific miRNAs are hard to be detected probably because of their low accumulation and/or insufficient stringency of the sequencing approach [79]. The advent of new sequencing technologies could make it possible to mine even species-/tissue- specific miRNAs with low abundance, and they havebeen successfully used on Arabidopsis thaliana, Oryza sativa, Poplus tricocarpa, Medicago tuncatula, Gossypiumhirsutum, Zea mays, Arachis hypogaea L., Solanumlycopersicum, Citrus trifoliate, Vitis vinifera and Vitis amurensis Rupr. [1021].

Grapevine (Vitis vinifera L.) is one of the most economically important fruit crops worldwide and has nutritional and processing properties [22]. In recent years, sequencing of small RNA libraries from different grapevine cultivars or tissues has severally been reported where a large number of Vvi-miRNAs were identified [7, 1821]. Despite this, there are no reports on the study of response of Vvi-miRNAs to phytohormones. Phytohormones are important endogenous signals and regulators involved in plant growth and development [23], and they are classified into auxins, gibberellins(GA3), cytokinins, abscisic acid(ABA) and ethylene. All of these phytohormones act at low concentrations to regulate different aspects of plant growth and development to varying degrees [2325]. Among the hormones, GA3 play significant regulatoryroles in early berry expansion, berry set and berry ripening [26, 27]. Till now, it is still unclear on how GA3 participates in the regulation of the complicated developmental processes of grapevine berry.

It has been reported that phytohormone signaling pathways can be effectively controlled by modulation of positive and negative regulators during plant growth and development [4]. Among the modulators of phytohormones, miRNA was recently found to be a new growth regulator involved in phytohormone signaling [2830], with several studies showing the interactions between miRNAs and phytohormones in various plant responses. For instance, GA3 modulates the expression of miR159, while miR159 regulates the development of Arabidopsi s anthers and seeds by cleaving the GAMYB gene, during Arabidopsis anther development [28] and seed germination [29]. In strawberry, miR159 interacts with GAMYB during the course of receptacle development, and both of them act in a joint fashion to respond, in part, to changes in endogenous GA3 levels [30]. These studies strongly demonstrated that miRNAs are involved in the GA3 signaling process.

To account for the roles of miRNAs in response to GA3 during grapevine berry development, we constructed two small RNA libraries from mixed tissue samples of grapevine berries sprayed with GA3 (treatment) and with water (the control). The grapevine cultivar used for this study is ‘Summer Black’ (hybrid of V. vinifera?×?V. labrusca), an elite table grapevine cultivar native to Japan. After high-throughput sequencing, we identified a number of conserved and non-conserved Vvi-miRNAs responsive to GA3, and we further analyzed their potential role in mediation of GA3-induced regulation of grapevine berry growth and development. Further, qRT-PCR was utilized to analyze the expression of Vvi-miRNAs in different development stages of the grapevine berries subjected to exogenous GA3 and in exogenous GA3-free berries. Lastly, an attempt to elucidate the regulatory functions of Vvi-miRNAs being responsive to GA3 during grapevine berry development was done.

Results

Characterization of the Vvi-miRNAs from deep sequencing of grapevine sRNA libraries

To identify GA3-responsive miRNAs in grapevine berries, two small RNA libraries from grapevine berries treated with GA3 (GA3 treatment) and sprayed with water (the control) were constructed. Solexa, a high throughput sequencing technology, was employed to sequence these libraries, leading to a generation of 16,231,320 and 16,486,660 clean reads from GA3 treatment and the control libraries, respectively. All these clean reads were those from removal of adaptor, insert, polyA, and RNAs shorter than 18nt in length (Table 1). About 4,265,160 (GA3 treatment) and 4,326,915 (the control) clear reads could be mapped to the grapevine genome published in 2007 [31], and miRNA, tRNA, siRNA, snRNA, snoRNA, rRNA, repeat regions, exon and intron RNA reads were annotated. In addition, 9,075,238 and 9,207,588 un-annotated reads were used for prediction of new Vvi-miRNAs in GA3 treatment and the control grapevine berries, respectively (Table 1).

Table 1 Distribution of small RNAs among different categories in control and GA 3 Treated grapevine berries

The size distribution of all sRNAs was found to be uneven ranging from 18nt to 30nt long, with the majority being 19-25nt long (Figure 1). The sRNAs of 21 nt and 24nt formed two major classes, occupying 38.35% and 32.06% (Figure 1) of the total, respectively, an observation which is in agreement with some previous reports in grapevine and tomato [11, 16, 21, 31], but contrasts to those reported in Arabidopsis, rice and peanut [3235]. The two main peaks of sRNAs in this study were 24 and 21nt long, while the most common sRNAs in Taxus chinensis[35] and Citrus trifoliate[17] were those with 21 nt in length. These cases suggested that some differences might exist in the sRNA biogenesis pathways in various plants. In addition, analysis of the first nucleotide of 18-25nt long sRNAs indicated that many sRNAs started with a uridine (U) at their 5’-ends and most of them are 21 nt and 22 nt long, with the former being most outstanding in number (Figure 2). Similar to other plants, most miRNAs here were those with 21 and 22 nt in length and they also begin with a 5’uridine, which is one of the important characteristic features of miRNAs [10, 11, 36].

Figure 1
figure1

Size distribution of unique small RNA sequences from grapevine.

Figure 2
figure2

First nucleotide bias of 18-30nt sRNA tags.

High throughput sequencing can verify a large number of known miRNAs together with the identification of novel specific miRNAs even with the low abundance in organisms. From two sRNA libraries in this study, we first searched for known Vvi-miRNAs by comparing our libraries with known miRNAs from other plant species in miRBase 19.0 (http://www.mirbase.org/). A total of 122 known Vvi-miRNAs were sequenced both in the control and GA3-treatment libraries, and they belonged to 27 conserved miRNA families according to a comparative genomics-based analysis in different plant species [37]. Eighteen of the 27 Vvi-miRNA families contained many members (Table 2), with four families (Vvi-miR169, Vvi-miR156, Vvi-miR166, Vvi-miR171 and Vvi-miR399) possessing 19, 8, 7, 8, and 7 members, respectively. Another nine Vvi-miRNA families (Vvi-miR162, Vvi-miR168, Vvi-miR390, Vvi-miR397, Vvi-miR408, Vvi-miR477, Vvi-miR479, Vvi-miR482, and Vvi-miR828) had only one member each. Among the known Vvi-miRNAs, the Vvi-miR166 family had the most abundant reads accounting for 63.2% of all the conserved miRNA reads. In this family, the number of Vvi-miR166h reads was over 500,000 in the two libraries, followed by Vvi-miR156, Vvi-miR168, Vvi-miR167, Vvi-miR479, Vvi-miR482 families, whose redundancies were more than several ten thousands. However, other miRNA families like Vvi-miR393, Vvi-miR394, Vvi-miR398 and Vvi-miR399 had only a few or tens of reads. Interestingly, the number and abundance of members of different Vvi-miRNA families between control and GA3 libraries were significant divergence (Table 2; Figure 3), which in turn could reflect the discrepancy in their potential functions during the development of grapevine berries responsive to gibberellin.

Table 2 Conserved Vv-miRNAs identified and normalized counts (NC) in control and GA 3 treated grapevine berries
Figure 3
figure3

Number of members from Vv-miRNA families in control and GA 3 treated grapevine berries.

Bioinformatics analysis of the sequencing data could identify novel Vvi-miRNAs based on the criteria of novel miRNA annotations developed by Meyers et al. (2008) [38]. The genomic sequences with flanking un-announced sequences were extracted. Among these sequences, 90 candidate Vvi-miRNAs were firstly uncovered, and the hairpin structures of their precursors’ could be predicted (Additional file 1). The miRNA* sequences available is of vital evidence supporting the release of miRNA duplex from the predicted foldback structure [16]. Among these 90 potential novel Vvi-miRNAs, 28 had their Vvi-miRNA* sequences were detected (Additional file 2), and the sequences of 45 novel potential Vvi-miRNAs started with a 5’ uridine, the important feature of miRNAs (Additional file 2). All these confirmed the existence of these novel candidate Vvi-miRNAs in grapevine. In this study, it was also noted that the novel Vvi-miRNAs had fewer reads than the conserved ones. 72% of the former had only several or dozens of copies except for Vvi-miRC20, Vvi-miRC36, Vvi-miRC43, Vvi-miRC44, Vvi-miRC45, Vvi-miRC46, Vvi-miRC50, Vvi-miRC51, Vvi-miRC52, Vvi-miRC71, and Vvi-miRC74 that were sequenced more than one thousand times (Additional file 2). This agreed with other previously reported results [10, 11, 36], where most novel species-specific miRNAs were usually expressed at lower levels than their conserved counterparts and were much more spatiotemporally expressed.

Discovery of miRNAs responsive to exogenous GA3 in grapevine

To identify Vvi-miRNAs responsive to exogenous GA3, we normalized the counts of reads sequenced based on the systematical analysis, and compared the members and the normalized counts (NC) of different members of Vvi-miRNA families between the control and GA3 libraries. The results revealed that seven known Vvi-miRNAs (Vvi-miR159b, Vvi-miR160f, Vvi-miR171e, Vvi-miR171f, Vvi-miR171h, Vvi-miR393a and Vvi-miR393b) were found only in the GA3 treated grapevine tissues not in the control, suggesting they were probably those induced by GA3; another 21 known Vvi-miRNAs (Vvi-miR160d, Vvi-miR166c, Vvi-miR169a, Vvi-miR169e, Vvi-miR169j, Vvi-miR169m, Vvi-miR169s, Vvi-miR319b, Vvi-miR394a, Vvi-miR395a, Vvi-miR395b, Vvi-miR395c, Vvi-miR395d, Vvi-miR395e, Vvi-miR395h, Vvi-miR395i, Vvi-miR395j, Vvi-miR399b, Vvi-miR403a, Vvi-miR403b and Vvi-miR535a) were strongly repressed in GA3 treatment for they could only be detected in the control but not in the GA3 treatment (Table 2).

Thirty novel Vvi-miRNAs (Vvi-miRC02, Vvi-miRC03, Vvi-miRC05, Vvi-miRC08, Vvi-miRC10, Vvi-miRC11, Vvi-miRC13, Vvi-miRC14, Vvi-miRC16, Vvi-miRC17, Vvi-miRC21-Vvi-miRC26, Vvi-miRC33, Vvi-miRC39, Vvi-miRC41, Vvi-miRC42, Vvi-miRC54, Vvi-miRC55, Vvi-miRC57, Vvi-miRC59, Vvi-miRC61, Vvi-miRC64, Vvi-miRC65, Vvi-miRC68 and Vvi-miRC72) were found to be responsive to exogenous GA3, for they could only be detected in GA3 treated grapevines (Table 2); while another 16 novel Vvi-miRNAs (from Vvi-miRC75 to Vvi-miRC90) were repressed by exogenous GA3, and found only in control but not in GA3 treatment (Table 2). These results demonstrate that these groups of Vvi-miRNAs could be responsive to GA3 treatment.

Further comparison of the normalized counts of Vvi-miRNAs identified in GA3 treated and the control plants indicated that many Vvi-miRNAs exhibited drastic variations (increases/decreases over several times) in expression frequencies. As shown in Table 2 and Table 3, the expression of 137 Vvi-miRNAs was strongly responsive to exogenous GA3, with 58 of them induced by GA3, showing conspicuously up-regulated expression under GA3 treatment (pointed as ↑ in Table 2 and 3).Conversely, 51 Vvi-miRNAs were repressed by exogenous GA3 (pointed as ↓ in Table 2 and 3), leading to drastic reduction in their expression levels in GA3-treated grapevines. An interesting revelation was that diverse members of the same Vvi-miRNA family could exhibit conspicuous discrepancy in their responses to exogenous GA3. This aspect is best exemplified by Vvi-miR166 family, where Vvi-miR166d/e/f/g had more than 246,000 reads in control samples, but possessed only several dozen or even a little more reads in the GA3 treatment samples, while Vvi-miR166b was detected 37,952 times in control and 176,156 times in GA3 treatment (Table 2). Vvi-miR166h, however, had no distinct variations between the control and GA3 treatment (Table 2). Similar situations were found in other Vvi-miRNA families like Vvi-miR156, Vvi-miR164, Vvi-miR167, Vvi-miR403 and Vvi-miR535 (Table 2), suggesting the Vvi-miRNAs responsive to exogenous GA3 application possess multiple aspects and functions during the development of grapevine berries.

Table 3 Comparison of normalized counts (NC) of novel Vv-miRNAs between GA 3 treated and control grapevine berries

Expression patterns of Vvi-miRNAs responsive to exogenous GA3 during grapevine berry development

Spatiotemporal expression of grapevine miRNAs could not only provide clues to their physiological functions, but also give fundamental evidence supporting the existence of the miRNAs in grapevine. In this study, 53 Vvi-miRNAs (27 conserved Vvi-miRNAs and 26 novel candidate Vvi-miRNAs) detected in grapevine berries treated with GA3 were subjected to qRT-PCR expression analysis as described in Section 2. This could also be applied to the analysis of the degree of response of these Vvi-miRNAs to GA3 treatments. We screened for the expression profiles of miRNAs responsive to GA3 treatments in the diverse stages of grapevine berries by qRT-PCR. The results showed that the expression levels of 15 Vvi-miRNAs (Vvi-miR169d, Vvi-miR319c, Vvi-miR393a, Vvi-miR396a, Vvi-miR398a, Vvi-miR399a, Vvi-miRC03, Vvi-miRC04, Vvi-miRC05, Vvi-miRC08, Vvi-miRC10, Vvi-miRC13, Vvi-miRC19, Vvi-miRC26 and Vvi-miRC37) were up-regulated by GA3, while 13 (Vvi-miR167a, Vvi-miR171d, Vvi-miR395f, Vvi-miR397a, Vvi-miR408, Vvi-miR482, Vvi-miRC06, Vvi-miRC14, Vvi-miRC15, Vvi-miRC24, Vvi-miRC27, Vvi-miRC30, Vvi-miRC38) were down-regulated by GA3, and seven (Vvi-miR156d, Vvi-miR166h, Vvi-miR390, Vvi-miR477, Vvi-miRC23, Vvi-miRC29, and Vvi-miRC36) were un-affected by GA3 (see Figure 4).

Figure 4
figure4

Expression patterns of miRNAs in grapevine berries underdifferent GA 3 treatments and control.

Further analysis of the expression results of 16 known Vvi-miRNAs responsive to GA3 both from qRT-PCR and high throughput sequencing revealed that the qRT-PCR expressions of 11 were consistent with the results from high throughput sequencing, which could not only demonstrated the reliability of these two technologies, but confirm the characteristics of these known GA3 responsive Vvi-miRNAs. On the other hand, it was discernible that the expression analysis results of 18 novel GA3 responsive Vvi-miRNAs from both qRT-PCR and high throughput sequencing were not of the same phenomenon, for only three (Vvi-miRC19, Vvi-miRC14 and Vvi-miRC29) was in the same way. This discrepancy in behavior of conserved and novel Vvi-miRNAs requires further research; while the consistency of the expression results, for the most conserved and a few novel Vvi-miRNAs, from the assays of the qRT-PCR and high throughput may provide some evidence supporting the fact about these Vvi-miRNAs responsive to GA3.

Characterization of potential target genes for novel Vvi-miRNAs responsive to GA3

To further comprehend the functions of these newly identified species-specific Vvi-miRNAs, it was essential to search for their target genes. Following the rules suggested by Schwab et al. (2005) [39], we searched the grapevine transcript database (http://www.genoscope.cns.fr/spip), and predicted a total of 117 putative target genes for 29 novel Vvi-miRNAs (Additional file 3). Among the 29 novel Vvi-miRNAs, 11 had multiple target genes, as exemplified by Vvi-miRC03 with 19 target genes, indicating these Vvi-miRNAs might possess comprehensive functions in grapevine. Based on orthologous functional annotation in other plants, these target genes seemed to be functionally involved in glucose metabolism (DNA glycosylase domain-containing protein), aromatic substance biosynthesis in berry (lipoxygenase, nigralipoxygenase), signal transduction (receptor-like protein kinase, serine-threonine protein kinase, protein phosphatase), stress resistance (ankyrin repeat-containing protein, disease resistance protein,cc-nbs-lrrresistance protein, domain-containing disease resistance protein, leucine-rich repeat receptor-like protein kinase), etc. (Figure 5; Additional file 3). We also noted that there were a large number of potential target genes with unknown functions for novel Vvi-miRNAs, especially those responsive to GA3 (Additional file 3). More in depth investigation on their functions would be essential for thorough understanding of the mechanisms of grapevine flower and fruit development and of the formation of berry quality.

Figure 5
figure5

Functional classification of target genes for novel Vv-miRNAs responsive to exogenous GA 3 .

Verification of potential Vvi-miRNA target genes using 5’-RLM-RACE

To verify the nature of potential miRNA targets and to study the Vvi-miRNAs’ regulation on their target genes, a modified RLM-RACE experiment was set up and used to map the cleavage sites in four predicted Vvi-miRNA target genes. Results showed that the cleavage sites of these four miRNA target genes GSVIVT01000639001, GSVIVT01026728001, GSVIVT01037667001 and GSVIVT01037667001 for Vvi-miRC15, Vvi-miRC23, Vvi-miRC60 and Vvi-miRC72 in this study is at the nucleotide that pairs with the 9th and/or 10th and/or 11th nucleotide of the corresponding miRNAs (Figure 6), consistent with previous related reports [7, 16, 18, 19, 40]. The GSVIVT01000639001, GSVIVT01026728001, GSVIVT01037667001 and GSVIVT01037667001 were confirmed as the true targets of Vvi-miR015, Vvi-miR023, Vvi-miR060 andVvi-miR072, respectively. Functional analysis indicated that GSVIVT01000639001, GSVIVT01026728001, GSVIVT01037667001 and GSVIVT01037667001 could be of cc-nbs-lrr resistance protein, nbs-lrrresistance protein, leucine-rich repeat family protein, and lipoxygenase (Additional file 3).

Figure 6
figure6

Verification of target genes for Vvi-miRNAs by RLM-5’-RACE.

Discussion

To verify the hypothesis that miRNAs play some role inregulating plant response to GA3[28], two small RNA libraries from grapevine berries treated with GA3 and without GA3 were constructed and used for high-throughput sequencing of sRNAs, from which a total of 212 Vvi-miRNAs (both known and novel), with 137 of them being firstly found to be responsive to GA3 in grapevine. Further comparison of our dataset from this work to Vvi-miRNAs reported earlier by us [20] revealed that since the grapevine materials used in these studies was the same grapevine cultivar, most of the Vvi-miRNAs could be observed in these two studies, while only a few members (such as Vvi-miR171f, Vvi-miR171h) could be discovered in one of both studies, which may be derived from the differences in the development stages of grapevine materials used in these studies or GA3 induced/depressed Vvi-miRNAs in this work.

The breakthroughs in sequencing technology have provided the most powerful tool for miRNA discovery, and one of the advantages of the thorough approach is its ability to reveal novel miRNAs. In this work, a large number of novel candidate Vvi-miRNAs were uncovered, with the miRNAs* for some novel Vvi-miRNAs being detected too. The identification of miRNAs*of the candidate miRNAs provides convincing evidence for consideration of these Vvi-miRNAs as authentic [38]. Importantly, the expression levels of some novel Vvi-miRNAs and their miRNAs* under GA3 treatments had much more diverse variation compared to the control. The best examples of this phenomenon are Vvi-miRC20/Vvi-miRC20*, Vvi-miRC34/Vvi-miRC34*, Vvi-miRC52/Vvi-miR C52* and Vvi-miRC53/Vvi-miRC53* (Figure 7). The expressions of Vvi-miRC20 and Vvi-miRC34 were up-regulated in GA3 treatments compared to the control, while the expression of their miRNAs* was down-regulated. On the contrary, the expressions of Vvi-miRC52 and Vvi-miRC53 were down-regulated in GA3 treatments, while those of their miRNAs* were up-regulated (Figure 7). The reasons for these discrepancies have not been clearly elucidated. The drastic variation in expression levels of the Vvi-miRNA* under GA3 treatments could indicate these miRNAs* were important regulatory genes like Vvi-miRNAs [41].

Figure 7
figure7

Comparison of expression modes for Vv-miRNAs and Vv-miRNAs*.

From high throughput sequencing and qRT-PCR analysis, it was further revealed that the conserve Vvi-miRNAs responsive to GA3 had the higher consistency between these two assays, while the novel ones showed some distinct discrepancy. This might be derived from the fact that the conserved Vvi-miRNAs possessed the higher conservation of development and function, while the novel ones had high specificity. Other, the samples used in the high throughput sequencing were the mixed materials from grapevine berries of several development stages, while those used for the qRT-PCR were of three stages berries (1 week after flowering (WAF), 5 WAF, 9 WAF). These could be the reasons explaining the expression levels of novel Vvi-miRNAs from qRT-PCR and high throughput sequencing had more apparent differences than those of conserved Vvi-miRNAs. In addition, the predication of potential target genes for novel Vvi-miRNAs responsive to exogenous GA3 and functional annotations of their orthologous target genes in other plants revealed that 20 genes were related to stress resistance. Whether or how these GA3 responsive novel Vvi-miRNAs are involved in the regulation of stress resistance can call for further research. The functions of one-third of the novel Vvi-miRNAs were unknown, indicating that more studies need to be performed on these novel miRNAs to elucidate their functions in the growth of grapevine.

Conclusions

Deep sequencing of short RNAs from grapevine berries in GA3 treatment and the control identified 137 GA3-responsive miRNAs, of which 28 exhibited different expression profiles of response to GA3 in the diverse developmental stages of grapevine berries. These Vvi-miRNAs identified might be involved in the grapevine berry development and response to various environments.

Methods

Plant material

Mixed samples of young berries (one week after flowering, WAF1) large berries (five week after flowering, WAF5), and old berries(nine week after flowering, WAF9) treated with 50 mg/lGA3, were collected in 2011 from four-year old ‘Summer Black’ grapevine (hybrids of V. vinifera and V. labrusca) trees grown under common field conditions at the Nanjing Agricultural University fruit farm, Nanjing, China. Each type of samples had three replicates during deep sequencing and qRT-PCR. All the samples were immediately frozen in liquid nitrogen and stored at −80°C until use.

Small RNA library construction and sequencing

Mixed Summer Black’ grapevine young berries (one week after flowering), large berries (five week after flowering after flowering), and old berries (nine week after flowering) treated with GA3, wereusedforRNAextraction. The total RNA samples were first extracted usingour modified CTAB method [20]. Isolation of small RNAs and preparation of small RNA libraries were performed based on the procedure of Wang et al. (2011) [20]. sRNAs were first separated from the total RNA by size fractionation with 15% TBE urea polyacrylamide gel (TBU) and small RNA regions correspondingto the 18–30 nucleotide bands in the marker lane were excised and recovered. The 18–30 nt small RNAs were 5’ and 3’ RNA adapter-ligated by T4 RNA ligase and at each step, length validated and purified by TBU electrophoretic separation. The adapter-ligated smallRNA was subsequently transcribed into cDNA by SuperScript II Reverse Transcriptase (Invitrogen) and PCR amplified using primers that annealed to the ends of the adapters. The amplified cDNA constructs were purified and recovered. 18 ng cDNA was loaded into the Illumina 1 G Genome Analyzer for sequencing.

Bioinformatics analysis and identification of Vvi-miRNA

To identify conserved and potential Vvi-miRNAs in grapevine, the raw sequences were processed as described by Sunkaret al. (2005) [8]. All sRNAs sequences from 18nt to 30nt obtained were removed from the vector sequences, then the modified sequences were further subjected to removal of rRNA, tRNA, snRNA, snoRNA and those containing the polyA tails, and finally the remaining sequences were compared against known plant Vvi-miRNAs in the miRBase [33]. Only the high matching (0-3mismatches) sequences were considered as conserved Vvi-miRNAs. To study potential Vvi-miRNAs precursor sequences, all sRNAs from grapevine were aligned against the grapevine genome and then the miRNA candidates were processed by miRCat (http://srna-tools.cmp.uea.ac.uk/) [33], using default parameters, to generate the secondary structures (Additional file 1).

qRT-PCR validation of miRNA expression

The template for RT-PCR was the miRNA-enriched library mentioned above. To amplify the Vvi-miRNAs from the reverse transcribed cDNAs, we used the Vvi-miRNA precise sequences as the forward primers and the mirRacer 3’Primer as the reverse primer [9] (Additional file 3). RT-PCR was conducted with the Rotor-Gene 3000 (Corbett Robotics, Australia) and the Rotor-Gene software version 6.1. For each reaction, 1 μL of diluted cDNA (equivalent to about 100 pg of total RNA) was mixed with 10 μL of 2X SYBR green reaction mix (SYBR® Green qRT-PCR Master Mix; Toyobo, Osaka, Japan), and 5 pmol each of the forward and the reverse primers were added in a final volume of 20 μL. The conditions for the PCR amplification were as follows: polymerase activation at 95°C for 1 min, then 95°C for 1 min, followed by 50 cycles of 95°C for 15 s, 95°C for 15 s, 60°C for 20 s, and 72°C for 20 s. The comparative quantification procedure was used to determine relative expression levels, and the homologous genes of the Arabidopsis 5.8S rRNA in grapevine berries was previously used as a reference gene in the qPCR detection of miRNAs in Arabidopsis[41]. The data were analyzed with an R2 above 0.998 using the LinRegPCR program [42].

Prediction of potential target mRNAs for Vvi-miRNAs

Target predictions were performed based on methods described by Schwab et al.(2005) [39]. Putative Vvi-miRNAs were first blasted against the grapevine unigene database on the Genoscope (http://www.genoscope.cns.fr/). BLASTn hits possessing less than four mismatches were chosen as the candidate targets, and then BLASTx was used to obtain their putative functions. The sequences of predicted targets and their functions are shown in Additional file 4.

Data access

The sRNA sequence data from this study have been submitted to Gene Expression Omnibus (GEO) under accession No. at website: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dlihnquimuscezm&acc=GSE3973.

References

  1. 1.

    Chen X: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303: 2022-2025. 10.1126/science.1088060.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Lauter N, Kampani A, Carlson S, Goebel M: Moose SP: microRNA172 downgulates glossy15 to promote vegetative phase change in maize. Proc Natl Acad Sci USA. 2005, 102: 9412-9417. 10.1073/pnas.0503927102.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  3. 3.

    Fu JH, Chiou TJ, Lin SI, Aung K, Zhu JK: A miRNA involved in phosphate-starvation response in Arabidopsis. Curr Biol. 2005, 15: 2038-2043. 10.1016/j.cub.2005.10.016.

    Article  Google Scholar 

  4. 4.

    Huq E: Degradation of negative regulators: a common theme in hormone and light signaling networks?. Trend Plant Sci. 2006, 11: 4-6.

    CAS  Article  Google Scholar 

  5. 5.

    Phillips J, Dalmay T, Bartels D: The role of small RNAs in abiotic stress. FEBS Letters. 2007, 581: 3592-3597. 10.1016/j.febslet.2007.04.007.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Carthew RW, Sontheimer EJ: Origins and mechanisms of miRNAs and siRNAs. Cell. 2009, 136: 642-655. 10.1016/j.cell.2009.01.035.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  7. 7.

    Carra A, Mica E, Gambino G, Pindo M, Moser C, Pe ME, Schubert A: Cloning and characterization of small non-coding RNAs from grape. Plant J. 2009, 59: 750-763. 10.1111/j.1365-313X.2009.03906.x.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Sunkar R, Girke T, Jain PK, Zhu JK: Cloning and characterization of microRNAs from rice. Plant Cell. 2005, 17: 1397-1411. 10.1105/tpc.105.031682.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  9. 9.

    Wang C, Shangguan LF, Nicholas KK, Wang XC, Han J, Song CN, Fang JG: Characterization of microRNAs identified in a table grapevine cultivar with validation of computationally predicted grapevine miRNAs by miRRACE. PLoSone. 2011, 6: e21259-10.1371/journal.pone.0021259.

    CAS  Article  Google Scholar 

  10. 10.

    Allen E, Xie Z, Gustafson AM, Sung GH, Spatafora JW, Carrington JC: Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsisthaliana. Nat Genet. 2004, 36: 1282-1290. 10.1038/ng1478.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Rajagopalan R, Vaucheret H, Trejo J: Bartel DP: A diverse and evolutionarily fluid set of microRNAs inArabidopsis thaliana. Genes Dev. 2006, 20: 3407-3425. 10.1101/gad.1476406.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  12. 12.

    Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate Vv-miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8: 25-10.1186/1471-2229-8-25.

    PubMed Central  PubMed  Article  Google Scholar 

  13. 13.

    Szittya G, Moxon S, Santos DM, Jing R, Fevereiro MP, Moulton V, Dalmay T: High-throughput sequencing of Medicagotruncatula short RNAs identifies eight new Vv-miRNA families. BMC Genomics. 2008, 9: 593-10.1186/1471-2164-9-593.

    PubMed Central  PubMed  Article  Google Scholar 

  14. 14.

    Ruan MB, Zhao YT, Meng ZH, Wang XJ, Yang WC: Conserved miRNA analysis in Gossypiumhirsutum through small RNA sequencing. Genomics. 2009, 94: 263-268. 10.1016/j.ygeno.2009.07.002.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Zhao CZ, Xia H, Frazier TP, Yao YY, Bi YP, Li YP, Li AQ, Li MJ, Zhang BH, Wang XJ: Deep sequencing identifies novel and conserved microRNAs in peanuts (Arachishypogaea L.). BMC Plant Biol. 2010, 10: 3-10.1186/1471-2229-10-3.

    PubMed Central  PubMed  Article  Google Scholar 

  16. 16.

    Moxon S, Jing R, Szittya G, Schwach F, RusholmePilcher RL, Moulton V, Dalmay T: Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res. 2008, 18: 1602-1609. 10.1101/gr.080127.108.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  17. 17.

    Song CN, Wang C, Zhang CQ, Korir NK, Yu HP, Ma ZQ, Fang JG: Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics. 2010, 11: 431-10.1186/1471-2164-11-431.

    PubMed Central  PubMed  Article  Google Scholar 

  18. 18.

    Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, Burgyan J: Identification of grapevine microRNAs and their targets using highthroughput sequencing and degradome analysis. The Plant Journal. 2010, 62: 960-976.

    CAS  PubMed  Google Scholar 

  19. 19.

    Mica E, Piccolo V, Delledonne M, Ferrarini A, Pezzotti M, Casati C, Fabbro CD, Valle G, Policriti A, Morgante M, Pesole G, Enrico MP, Horner DS: Correction: High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitisvinifera. BMC Genomics. 2010, 11: 109-10.1186/1471-2164-11-109.

    PubMed Central  PubMed  Article  Google Scholar 

  20. 20.

    Wang C, Wang X, Nicholas KK, Song C, Zhang C, Li X, Han J, Fang J: Deep sequencing of grapevine flower and berry short RNA library for discovery of new microRNAs and validation of precise sequences of grapevine microRNAs deposited in miRBase. PhysiologiaPlantarum. 2011, 143: 64-81.

    CAS  Google Scholar 

  21. 21.

    Wang C, Han J, Liu C, Kibet KN, Kayesh E, Shangguan L, Li X, Fang J: Identification of microRNAs from Amur grape (Vitisamurensis Rupr.) by deep sequencing and analysis of microRNA variations with bioinformatics. BMC Genomics. 2012, 13: 122-10.1186/1471-2164-13-122.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  22. 22.

    He PC: Grapevine [M]. 1998, BeiJing: Chinese Agricultural Press

    Google Scholar 

  23. 23.

    Leyser O: Plant hormones. Curr. Biol. 1998, 8: R5-R7. 10.1016/S0960-9822(98)70006-5.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Klee H: Hormones are in the air. PNAS. 2003, 100: 10144-10145. 10.1073/pnas.1934350100.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  25. 25.

    Abu-Zahra TR: Berry size of Thompson seedless as influenced by the application of Gibberellic acid and can girdling. Pak J Bot. 2010, 42: 1755-1760.

    Google Scholar 

  26. 26.

    Abu-Zahra TR, Salameh NM: Influence of Gibberellic acid and can girdling on berry size of black magic grape cultivar. Middle-East Journal of Scientific Research. 2012, 11: 718-722.

    CAS  Google Scholar 

  27. 27.

    Achard P, Herr A, Baulcombe DC, Harberd NP: Modulation of floral development by a gibberellin-regulated microRNA. Development. 2004, 131: 3357-3365. 10.1242/dev.01206.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Reyes JL, Chua NH: ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant Journal. 2007, 49 (4): 592-606. 10.1111/j.1365-313X.2006.02980.x.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Gsukasi F, Donaire L, Casanal A, Martinez-Priego L, Botella MA, Medina-Escobar N, Llave C, Valpuesta V: Two strawberry miR159 family members display developmental-specific expression patterns in fruit receptacle and cooperatively regulate Fa-GAMYB. New phytologist. 2012, 195: 47-57. 10.1111/j.1469-8137.2012.04134.x.

    Article  Google Scholar 

  30. 30.

    Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Fabbro CD, Alaux M, Gaspero GD, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Clainche IL, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pè ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P: French-Italian Public Consortium for Grapevine Genome Characterization: the grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Henderson IR, Zhang XY, Lu C, Johnson L, Blake C, Meyers BC, Green PJ, Jacobsen SE: Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nature Genet. 2006, 38: 721-725. 10.1038/ng1804.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Morin RD, Aksay G, Dolgosheina E, Ebhardt HA, Magrini V, Mardis ER, Sahinalp SC, Unrau PJ: Comparative analysis of the small RNA transcriptomes of Pinuscontorta and Oryza sativa. Genome Res. 2008, 18: 571-584. 10.1101/gr.6897308.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  33. 33.

    Wang C, Leng XP, Zhang YY, Kayesh E, Zhang YP, Sun X, Fang JG: Transcriptome-wide analysis of dynamic variations in regulation modes of grapevine microRNAs on their target genes during grapevine development. Plant Molecular Biology. 2014, 84 (3): 269-285. 10.1007/s11103-013-0132-2.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Qiu DY, Pan XP, Wilson IW, Li F, Liu M, Teng W, Zhang B: High throughput sequencing technology reveals that the taxoid elicitor methyl jasmonate regulates microRNA expression in Chinese yew (Taxuschinensis). Gene. 2009, 436: 37-44. 10.1016/j.gene.2009.01.006.

    CAS  PubMed  Article  Google Scholar 

  35. 35.

    Lindow M, Krogh A: Computational evidence for hundreds of non-conserved plant microRNAs. BMC Genomics. 2005, 6: 119-10.1186/1471-2164-6-119.

    PubMed Central  PubMed  Article  Google Scholar 

  36. 36.

    Zhang XR, Luo GG, Wang RH, Wang J, Himelrick DG: Growth and developmental responses of seeded and seedless grape berries to shoot girdling. J Am SocHortSci. 2003, 128: 316-323.

    Google Scholar 

  37. 37.

    Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao XF, Carrington JC, Chen XM, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi YJ, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20: 3186-3190. 10.1105/tpc.108.064311.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  38. 38.

    Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M: Weigel D:Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005, 8 (4): 517-527. 10.1016/j.devcel.2005.01.018.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Jagadeeswaran GR, Zheng Y, Sumathipala N, Jiang HB, Arrese EL, Soulages JL, Zhang WX, Sunkar R: Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development. BMC Genomics. 2010, 11: 52-10.1186/1471-2164-11-52.

    PubMed Central  PubMed  Article  Google Scholar 

  40. 40.

    Mallory AC, Bartel DP, Bartel B: MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell. 2005, 17: 1360-1375. 10.1105/tpc.105.031716.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  41. 41.

    Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005, 39: 519-525. 10.2144/000112010.

    CAS  PubMed  Article  Google Scholar 

  42. 42.

    Ramakers C, Ruijter JM, Deprez RH, Moorman AF: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neuroscience Lett. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

This research was supported by Project Funded by the Natural Science Foundation of China (NSFC) (No. 31301759), China Postdoctoral Science Foundation(2013 M531373), Postdoctoral Science Foundation of Jiangsu Province(1301116C), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the National Science Foundation of China (No. 60901053), and the Nanjing Agricultural University Youth Science and Technology Innovation Fund (KJ2013013).

Author information

Affiliations

Authors

Corresponding authors

Correspondence to Jinggui Fang or Chen Wang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JH and CW carried out the molecular genetic studies, participated in the sequence alignment. YY and XS carried out the RT-PCR arrays. XL and CS participated in the sequence alignment. CW performed the statistical analysis. JH and CW drafted the manuscript. JF and CW conceived of the study, and participated in its design and coordination. JF revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Han, J., Fang, J., Wang, C. et al. Grapevine microRNAs responsive to exogenous gibberellin. BMC Genomics 15, 111 (2014). https://doi.org/10.1186/1471-2164-15-111

Download citation

Keywords

  • Grapevine
  • Berry
  • microRNAs
  • Exogenous gibberellin
  • High throughput sequencing