- Research article
- Open Access
Comparative RNA-seq based transcriptomic analysis of bud dormancy in grape
BMC Plant Biology volume 17, Article number: 18 (2017)
Bud dormancy is an important biological phenomenon of perennial plants that enables them to survive under harsh environmental circumstances. Grape (Vitis vinifera) is one of the most grown fruit crop worldwide; however, underlying mechanisms involved in grape bud dormancy are not yet clear. This work was aimed to explore the underlying molecular mechanism regulating bud dormancy in grape.
We have performed transcriptome and differential transcript expression analyses of “Shine Muscat” grape buds using the Illumina RNA-seq system. Comparisons of transcript expression levels among three stages of dormancy, paradormancy (PD) vs endodormancy (ED), summer buds (SB) vs ED and SB vs PD, resulted in the detection of 8949, 9780 and 3938 differentially expressed transcripts, respectively. Out of approximately 78 million high-quality generated reads, 6096 transcripts were differentially expressed (log2 ratio ≥ 1, FDR ≤ 0.001). Grape reference genome was used for alignment of sequence reads and to measure the expression level of transcripts. Furthermore, findings obtained were then compared using two different databases; Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), to annotate the transcript descriptions and to assign a pathway to each transcript. KEGG analysis revealed that secondary metabolites biosynthesis and plant hormone signaling was found most enriched out of the 127 total pathways. In the comparisons of the PD vs ED and SB vs ED stages of grape buds, the gibberellin (GA) and abscisic acid (ABA) pathways were found to be the most enriched. The ABA and GA pathways were further analyzed to observe the expression pattern of differentially expressed transcripts. Transcripts related to the PP2C family (ABA pathway) were found to be up-regulated in the PD vs ED comparison and down-regulated in the SB vs ED and SB vs PD comparisons. GID1 family transcripts (GA pathway) were up-regulated while DELLA family transcripts were down-regulated during the three dormancy stages. Differentially expressed transcripts (DEGs) related to redox activity were abundant in the GO biological process category. RT-qPCR assay results for 12 selected transcripts validated the data obtained by RNA-seq.
At this stage, taking into account the results obtained so far, it is possible to put forward a hypothesis for the molecular mechanism underlying grape bud dormancy, which may pave the way for ultimate improvements in the grape industry.
Grape (Vitis vinifera) is the most widely grown fruit crop globally. The area under grape cultivation is approximately 7.8 million hectares with a production of about 67.5 million tons. The berries are categorized mainly into table grapes (fresh) and wine grapes (wine), as well as for several value-added products . China is the leading grape-producing country, accounting for 14% of the global grape production .
There are several developmental and metabolic processes that occur in the buds and twigs of grape plants during the winter period. These processes include enzyme synthesis, respiration, cell division, photosynthesis, growth stimulator production and growth inhibitor down-regulation. Dormancy is a controlling mechanism that enables woody perennials to adapt seasonal environmental changes and thus affects the following season’s vegetative growth and fruit production. Currently, global warming has a substantial influence on winter chilling accumulation and dormancy release of fruit trees . To ensure sustainable fruit production, it is necessary to investigate the underlying genetic factors responsible for controlling dormancy . Extended dormancy is a key hindrance for the large scale fruit production, including grape, in warm or mild winter regions under temperate and subtropical climates [5, 6]. Several studies have been conducted to determine the association between natural and chemical-induced ED, analyze gene expression during long and short photoperiods, and identify the transcript profile of bud development and signaling of bud dormancy break in grape [7–10]. Dormancy is generally classified into three main types: paradormancy (PD), endodormancy (ED), and ecodormancy (ECD) . PD is the plant growth suspension initiated by factors outside the meristem. It is essentially the effect of one organ on another and involves the dominance of apical buds. ED is regulated by internal growth inhibitors, even under favorable conditions; without exposure to cold temperature for a specific duration (chilling requirement), endodormant buds (EDBs) cannot initiate growth. Exposure to low temperature (2–9 °C) shifts the ED state of the plant to ECD. ECDBs can break and grow when exposed to suitable growth conditions . When EDB’s chilling requirement are fulfilled, the ED is released. EDBs steadily transition to the ECD state, especially under adverse environmental conditions. Summer buds (SB), which are green in color and small in size and grow on one side of winter buds that have no scales, can be observed after dormancy release during the new growth period and remain active for a short time during the transition from dormancy release to early summer dormancy. Like other perennial deciduous fruit plants, grape undergoes a characteristic dormant period during its growth cycle. In southeast China, grape buds fulfill their chilling requirement in the end of February and blossom in following spring. Inadequate cold accumulation hours during this period lead to irregular flowering, which consequently decreases fruit production.
The investigations have been made on dormancy at physiological as well as molecular levels in different deciduous fruits. MADS-box (DAM) genes associated with dormancy-have been isolated to investigate their expression pattern in some fruit plants during dormancy [12, 13]. For example, DAM1 through DAM6 have been identified in peach and Japanese apricot [14, 15], while MADS13-1, MADS13-2, MADS13-3, PpMADS1 and PpMADS2 were found in Japanese pear and Chinese white pear (Suli) [16, 17]. The expression profile of these genes during the induction and release of endodormancy indicated that DAMs serve as dose-dependent inhibitors of bud break . Additionally, several other genes are involved in the complex molecular network regulating dormancy in deciduous plants. Therefore, segregating single gene is not sufficient for illuminating underlying molecular processed associated with bud dormancy .
Recently, the next-generation sequencing (NGS) technology has uplifted the transcriptomic by allowing the RNA-sequencing using cDNA libraries on a large scale. RNA-seq is a highly efficient and modern tool that involves deep sequencing technologies to generate millions of short cDNA reads which is considerably more efficient than microarray analysis . In previous studies, RNA-seq was successfully applied to investigate dormancy based on direct sequencing of cDNAs in several woody plants using 454-pyrosequencing technology . Moreover, in another study the transcriptomic analysis revealed the dormancy-related regulatory pathways involving photoperiod, hormones and circadian clocks [20–22]. Although previous studies have investigated the physiological as well as the molecular mechanism of bud dormancy using the transcriptomic approach in deciduous fruits as well as other crops [13, 16, 23], no attempt has yet been made to study grape bud dormancy at the transcriptomic level.
This study was undertaken to investigate underlying molecular processes regulating bud dormancy in grape and to develop robust foundation for molecular research. RNA-seq technology was used to categorize and characterize the expression profile of differentially expressed genes (DEGs) during three different grape bud dormancy stages. This novel transcriptome and transcript expression profiling data generated through RNA-seq will offer an improved understanding of underlying molecular process of bud dormancy and will pave the way to identifying key genes involved in dormancy for the ultimate improvement of table grape industry.
Analysis of RNA-seq libraries
In this study, three cDNA libraries constructed from grape buds during three different stages were sequenced and generated 79.6 million sequence reads. After elimination of low-quality reads and adaptor sequences, 78.5 million clean reads (98.5% of the generated data) were recorded, which were then mapped to the reference genome of grape using HISAT . Furthermore, out of high-quality reads generated from the three samples, uniquely mapped reads were 73.28 to78%, while total mapped reads were 75.16 to 79.33% (Table 1).
Differential expression analysis of transcripts
To understand and interpret the results of the RNA-seq experiment, the differential expression patterns of transcripts were analyzed among the three different bud dormancy stages. From three different libraries, differential expression analysis identified 943 to 7596 transcripts with significant expression changes (p ≤ 0.05 and fold change ≥ 2). The different expression patterns among the three stages revealed that the maximum differences (7596 down-regulated transcripts and 2184 up-regulated transcripts) were examined between the SB and ED stages. In contrast, in the PD vs ED comparison, 2969 transcripts were up-regulated and 5980 were down-regulated, while in the SB vs PD comparison, 943 transcripts were up-regulated and 2995 were down-regulated. Whereas, in comparison between SB and ED stages, the maximum number of 1280 distinctive transcripts was observed, while fraction of unique transcripts were identified in the PD vs ED (1048) and SB vs PD (453) comparisons. Among these, 70 transcripts were commonly up-regulated and 565 transcripts were down-regulated in all three stages of dormancy (Fig. 1, Additional files 1, 2 and 3).
Cluster analysis of DEGs
A cluster analysis of transcript expression patterns with functional enrichment was performed using familiar log ratio values for the transcript expression analysis. The transcripts were arranged into three groups, SB vs PD, SB vs ED and PD vs ED. In the SB vs PD group, 969 transcripts (24.70%) were up-regulated and 2953 transcripts (75.29%) were down-regulated, while in the SB vs ED and PD vs ED groups, 2152 transcripts (54.86%) and 2907 transcripts (74.12%) were up-regulated and 1770 transcripts (45.13%) and 1015 transcripts (25.87%) were down-regulated, respectively. Split plots are shown for each cluster with the data presented as the means of the standard deviation of the RMKM expression values. The cluster analysis grouped up-regulated and down-regulated transcripts separately. A majority of transcripts were up-regulated; while a smaller number of transcripts were down-regulated (Fig. 2 and Additional file 4).
GO and KEGG analysis of DEGs
Gene Ontology based enrichment analysis was carried out using a threshold value (p-value ≤ 0.05) to evaluate the major biological functions of DEGs that are further classified into three main categories such as, cellular component (CC), molecular function (MF) and biological process (BP). BP category contained the majority of GO annotations (26,989; 42.15%) followed by MF (21,686; 33.87%) and CC (15,352; 23.97%). The major subcategories along with the analysis of all the transcripts among the three different stages of bud dormancy are shown in Fig. 3. The PD vs ED, SB vs ED and SB vs PD comparisons represent 26,434 (41.28%), 27,559 (43.04%) and 10,034 (15.67%) transcripts, respectively, of the total 64,027 transcripts annotated in GO major categories. A total of 15,352 transcripts were categorized as CC, with 6669 (43.44%) recognized in the PD vs ED comparison, 6642 (43.26%) in the SB vs ED comparison and 2041 (13.29%) in the SB vs PD comparison. Transcripts associated with the CC subcategories integral component of membrane (595; 8.92%, 632; 9.51%, 215; 10.53%) and nucleus (510; 7.64%, 500; 7.52%, 197; 9.65%) were identified in the PD vs ED, SB vs ED and SB vs PD comparisons, respectively. A total of 26,989 transcripts were categorized as BP, with 10,999 (40.75%) identified in the PD vs ED comparison, 11,582 (42.91%) in the SB vs ED comparison and 4408 (16.33%) in the SB vs PD comparison. Transcripts associated with the BP subcategories oxidation-reduction process (667; 6.06%, 712; 6.14%, 288; 2.48%) and metabolic process (534; 4.85%, 551; 4.75%, 199; 4.51%) were recognized in the PD vs ED, SB vs ED and SB vs PD comparisons, respectively. A total of 21,686 transcripts were categorized as MF, with 8766 (40.42%) identified in the PV vs ED comparison, 9335 (43.04%) in the SB vs ED comparison and 3585 (16.53%) in the SB vs PD comparison. Transcripts associated with the MF subcategories ATP binding (653; 7.44%, 730; 7.82%, 277; 7.72%) and DNA binding (282; 3.21%, 290; 3.10%, 128; 3.57%) were recognized in the PD vs ED, SB vs ED and SB vs PD comparisons, respectively (Table 2). A sum of 13,740 DEGs were allocated to 127 pathways (Additional files 5, 6 and 7). Based on KEGG analysis, biosynthesis of secondary metabolites with 1504 transcripts was the most enriched pathway, followed by plant hormone signal transduction (659 transcripts) and ribosome (299 transcripts) in three different dormancy stages (Fig. 4).
Transcripts related to plant hormone signal transduction and secondary metabolism pathways
In the present study, 1504 transcripts linked secondary metabolism pathways were identified in three dormancy stages. Out of which, 482 and 1022 were up and down-regulated during all three stages of dormancy. 10,312 DEGs were annotated in plant hormone signaling pathways, of which the ABA, gibberellin (GA), and ethylene signaling pathways were further analyzed. Sixteen transcripts were annotated as protein phosphatase 2C (PP2C) transcripts, out of which, 11 were up-regulated in the PD vs ED comparison. A large quantity of transcripts abundance of a gene annotated as serine/threonine-protein kinase (SnRK2) was lower in the PD vs ED comparison. In GA-responsive pathway, six out of the total 16 transcripts encoding DELLA proteins were found to be down-regulated in the PD vs ED comparison, while five transcripts were up-regulated in the SB vs ED comparison. In the ethylene response pathway, two transcripts annotated as ethylene response receptor (ETR) were down-regulated in the PD vs ED comparison, while three ETR transcripts were down-regulated in the SB vs PD comparison (Tables 3 and 4). Moreover, differential expression of genes involved in plant hormone signaling pathways was also identified. In the auxin biosynthesis pathway, four out of 15 transcripts encoding Aux-1 proteins showed up-regulation in the PD vs ED comparison. In the zeatin biosynthesis (cytokinin) pathway, 14 transcripts encoding CRE1 proteins were identified, with one transcript up-regulated in the PD vs ED comparison and 13 transcripts were down-regulated in the SB vs ED comparison.
Validation of DEGs by RT-qPCR
Twelve DEGs were chosen for RT-qPCR analysis to verify the precision and reproducibility of the transcriptome analysis results. In each case, the qRT-PCR assay results closely related to the transcript levels assessment by the RNA-seq analysis (Fig. 5).
Grape, being one of the most important fruit crops, is globally consumed fresh as well as in the form of several value-added products . Dormancy is a very complex and highly programmed mechanism used by perennial plants to cope with unfavorable environmental conditions. The beginning of dormancy requires sensing and development of regular environmental signals . In grape, a shorter photoperiod and low temperatures cause the alteration of buds into ED [26, 27]. Dormancy can be generally categorized into three dormant states, ED (growth suspension by factors outside the meristem), ED (growth inhibition by internal bud signals) and ECD (growth inhibition by momentary adverse ecological situations) . The molecular and physiological aspects of bud dormancy in grape have been previously examined in several studies [7–10]. This is first ever report on application of RNA-seq technique to classify a large number of transcripts from grape buds of different dormancy stages. Using a transcriptomic approach, we observed that the number and expression profiles of DEGs differed during dormancy stages. A sum of 8949, 9780 and 3983 transcripts were differentially expressed in the PD vs ED, SB vs ED and SB vs PD comparisons, respectively. Transcripts with a like expression patterns might be functionally correlated during bud dormancy. A cluster analysis of DEGs during three comparative dormancy stages was carried out to know the expression pattern of the 11,766 transcripts that were differentially expressed during dormancy stages. The cluster analysis revealed that the most of transcripts were up-regulated while a relatively smaller number of transcripts were down-regulated. Our findings revealed that a number of DEGs were highly expressed in SB vs ED than in the other two stages of dormancy. Previous studies showed that gene activity in black current was minimum at early stages of dormancy and maximum at the moment of bud break . In our study, very high transcript activity in SB vs ED as well as very low activity in SB vs PD was likely due to growth-conducive conditions or signaling from other plants. Additionally, using KEGG analysis, we found that these DEGs belonged to several pathways. Substantial variations were noticed in five pathways, secondary metabolites biosynthesis, ribosome, starch and sucrose metabolism in PD vs ED and SB vs ED stages, while secondary metabolites biosynthesis, signaling of plant hormone and flavonoid biosynthesis pathways were represented in SB vs PD stage of dormancy. Our findings were in consensus with previous work on Chinese pear, in which comparison of transcriptomic analysis between ED and ECD during the whole dormancy cycle showed substantial variations in five KEGG pathways, plant-pathogen interaction, metabolism of ether lipid, ribosome, endocytosis in glycerophospholipid and metabolic pathways [16, 17]. Enriched GO terms recognized in our study, oxidation-reduction process, hormone metabolism and jasmonic acid stimulus, were also in agreement with previous reports .
Oxidative stress is proposed to be an important process involved in ED release . Consistent with this perspective, H2O2 has been reported to be a signaling factor increasing the expression of genes related with release of ED . An increase in H2O2 levels take place earlier to release ED in grape buds, proposed that H2O2 could be a signal molecule that triggers gene expression for release of ED. Recent researches have figured out the key role of hydrogen cynamide and calcium signaling in bud break of Perlette grapevines . The higher expression of calcium signaling-related genes corresponds with the optimum bud break potentiation in V. riparia, additionally proposing a key role for calcium in the transition from ED to ECD . A significantly down-regulated group of 130 genes was identified during the alteration from ED to ECD at chilling accumulation time in grape and in leafy spurge, and included proline-rich protein, glutathione S-transferase, peroxidase,, serine decarboxylase, thaumatin, serine carboxy peptidase and xyloglucan endo-transglycosylase [12, 33]. Our data demonstrated the up-regulated expression of catalase along with down-regulated expression of some peroxidase genes among all three dormancy stages. Down-regulation of peroxidase genes and up-regulation of catalase genes could enhance or decrease the H2O2, thus increase release of ED. Therefore, further investigation into the relationship between activity of catalase and levels of H2O2 after ED is required. Generally, metabolic networks are controlled by hormone function and signaling. The involvement of ABA to maintain and promote bud dormancy in woody plants has been projected [34–36]. A gradual decreas of ABA contents during ED to ECD have been reported in leafy spurge and pear buds previously [29, 37] and peaked in poplar after few weeks of short days . Moreover, an ABA related transcript has showed down-regulation during the chilling phase essential for ED release in grape . Similar to these findings, our study showed higher ABA expression in the PD vs ED comparison, while lower expression was observed in the SB vs ED comparison. Based on previous reports, we speculate that ABA might play acrucial role in initiation and maintenance of ED in grape.
Gibberellin (GA) are plant hormones that control several growth processes including seed germination; stem elongation, growth regulation and dormancy. Previous reports have depicted the involvement of GA in bud break, and an increase in GA levels has been considered to be essential for ED release . GA signaling via GID1 receptors is essential for seed germination in Arabidopsis . Five transcripts in the GID1 and DELLA families were identified and validated by qRT-PCR in the present study. These transcripts also showed different expression patterns during the three dormancy stages. GID1 family transcripts were up-regulated while DELLA family transcripts were down-regulated during the three dormancy stages. Overall, these results suggested that GA was not associated with release of ED activities, with the exception of bud burst initiation.
Basipetally transported auxin is considered as a key signal regulating PD. Cytokinin synthesis is inhibited by auxin. Several genes have been identified in Arabidopsis and pea which involved in auxin-regulated growth inhibition . Cytokinin and auxin signaling have been identified in regulation of PD; however, their involvements in ED are not yet clear . The auxin and cytokinin-responsive transcripts are differentially expressed as plants alteration from PD to ED . In our study, transcripts related to signaling pathways of cytokinin and auxin showed lower expression in all three stages of dormancy. Based on previous studies, we speculate that auxin and cytokinin might be associated in PD and ED regulation of grape buds.
The functional category of identified transcription factors was significantly enriched in the transcript expression profile of this comparative study. Among these identified transcription factors, within the AP2-like transcription factor family, ERF subfamily with two transcripts was significantly enriched , while many of them can regulate the ethylene responses during dormancy and similar responses of ERF-like transcription factor have also been reported in poplar . In fact, potato, leafy spurge, and poplar all exhibited the momentary peak in ethylene or ethylene perception that is linked with ED induction as verified by several studies on similar aspect [37, 38, 41]. Another finding on leafy spurge showed contradictory results during PD as revealed by microarray analysis; at least ten ethylene responsive genes were highly induced but were repressed during Ed and ECD . In our study, transcripts related to ethylene signaling pathway showed synchronized expression patterns, with higher ETR levels in SB vs PD and lower levels of CTR1-like transcripts in PD vs ED. Based on our results, we suggest that ethylene signaling might be involved in endodormancy release.
As stated above, the results obtained so far allow for the development of a hypothesis regarding the molecular mechanism underlying bud dormancy. By comparing the transcriptomes among three stages, the potential contribution of various pathways in this method became evident. This work implicated several pathways, including plant hormone signaling as well as secondary metabolites biosynthesis. Further confirmation of most enriched pathways and DEGs will be the major emphasis of future studies.
Shine Muscat, the most popular table grape cultivar in Japan  and China due to its aroma and good taste, was used as the plant material in this study. Four-year-old grape plants were spaced at 6 m × 3 m apart under a rain shelter covered with polyvinyl film and supplemented with drip irrigation at Nanjing Agricultural University Vineyard located in Tangshan Valley, Nanjing, Jiangsu province, China. During the sampling period, plants were not pruned or chemically treated. Buds were harvested on February 02, (ED stage), May 19 (SB) and August 08 (PD stage) in 2015.
The dormancy stages of grape buds prior to constructing gene expression profile were defined as ED, SB, and PD. The growth in the ED stage is stopped due to low chilling exposure and factors within the meristem, while in the PD stage, plant growth is suspended due to factors outside the meristem. SB grows on one side of winter buds having no scales. No bud break was noticed on shoots sampled on 2nd February. These buds were considered to be in ED phase and the collected buds were designated endodormant buds (EDB). The bud samples collected on 19th May and 8th August were designated summer buds (SB) and paradormant buds (PDB), respectively. The samples were instantly frozen in liquid nitrogen and then kept at −80 °C until RNA extraction.
Preparation of RNA-seq libraries
Total RNA was extracted using Foregene RNA isolation kit (Foregene Co.Ltd, China) according to manufacturer’s instructions. RNA quality was checked with a 2200 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). Total RNA extracted from the three samples collected per dormancy stage was pooled into three sample stages. From each sample, to isolate poly (A) mRNA, 10 μg of total RNA was used to prepare Illumina RNA-seq libraries. From three biological replicates for each stage, each library was pooled by mixing equal quantities of RNA. An insert size of 200 bp was used for sequencing of each library using the Illumina HiSeqTM 2000 system following the manufacturer’s protocol.
Mapping of reads to the reference genome and gene annotation
The raw sequence data were filtered by removing adaptor sequences, low quality reads with more than 10% anonymous nucleotides (N) and 50% bases of quality value ≤5 by using hierarchical indexing for spliced alignment of transcripts (HISAT)  and standard parameters used for mapping (−−phred64 --n-ceil -q “L, 0, 0.05” -I 100 -X 1000 -t -p 6 --no-una) prior to mapping against a reference grape genome database. Clean reads were mapped to the Vitis vinifera reference genome (Assembly accession = GCF_000003745.3; Assembly version = 12X; http://ftp.ncbi.nlm.nih.gov/genomes/Vitis_vinifera/Assembled_chromosomes/seq/ vinifera) using the mapping software HISAT (version 0.1.6). For our data, a read length ≥ 100 bp was used and included reads mapped to the reference genome with ≤ two mismatches . Reads that failed to be mapped were cleaned and mapped to the genome again until a match was found (Fig. 6).
GO analysis and gene expression evaluation from RNA seq
To compare gene expression levels among three samples, the relative transcript level of each expressed transcript was normalized and calculated to the reads per kilobase of exon model per million mapped reads (RPKM) values . For all RPKM values of each transcript, the cutoff value was determined for shaping gene transcriptional activity based on a 95% confidence threshold. To obtain GO annotations, Blast2GO program was used (version 2.3.5) (https://www.blast2go.com/) for all the transcripts . Further, we performed GO enrichment analysis using GO seq  to classify genes or their products into terms (molecular function, biological process and cellular component) that are helpful in understanding the biological functions of the genes.
Differentially expressed genes (DEGs) and cluster analysis during the three stages of dormancy
DEG seq  and DEG seq2  were used to detect the differentially expressed genes. The p-value threshold was determined by FDR to account for multiple tests of significance. In this study, FDR threshold ≤ 0.001 and fold change ≥ 2 were adopted to observe the significance of the transcript expression differences . For pathway analysis, all DEGs were mapped to terms in KEGG database and then looked for significantly enriched pathway terms compared to the background genome. KEGG pathways fulfilling the criterion of a Bonferroni  corrected p-value ≤0.05 were defined as significantly enriched in DEGs. Cluster analyses of gene expression patterns in PD vs ED, SB vs ED and SB vs PD comparisons were performed using R package pheatmap . The sequences obtained from the Illumina sequencing were deposited in the NCBI Sequence Read Archive (accession number, GSE77119).
Real-time quantitative PCR (RT-qPCR) analysis of DEGs
Twelve genes were selected for validation using quantitative real-time PCR. Primer pairs were designed using Beacon Designer software (Premier Biosoft, version 7.0), which are listed in (Additional file 8). The qPCR reaction was performed in a total volume of 20 μl, containing 1 μl of diluted cDNA, 0.6 μl of reverse and forward primers, 7.4 μl of ddH2O, 0.4 μl of ROX and 10 μl of the PCR master mix (Thermo Fisher Scientific, Waltham, MA, USA). According to the standard protocol of the ABI 7300 system, the amplification program was performed as follows: 30 s at 95 °C, followed by 40 cycles of 5 s at 95 °C for and 30 s at 60 °C. To verify the formation of single peaks and to exclude the possibility of primer dimer and non-specific product formation, a melt curve (15 s at 95 °C, 60 s at 60 °C, and 15 s at 95 °C) was generated by the end of each PCR reaction. All reactions were performed in triplicate, including the non-template control reactions. In addition, the threshold cycles (Ct) of the triplicate reactions for each tested gene were averaged, and then the values were normalized to that of the control V. vinifera Actin gene (accession number XM_010659103) .
Arginine and proline metabolism
Amino sugar and nucleotide sugar metabolism
Biosynthesis of secondary metabolites
Carbon fixation in photosynthetic organisms
Cysteine and methionine metabolism
Carbohydrate metabolic process
Circadian rhythm plant
Differentially expressed genes
Flavone and flavonol biosynthesis
Fructose and mannose metabolism
Integral component of membrane
Metal ion binding
Porphyrin and chlorophyll biosynthesis
Porphyrin and chlorophyll metabolism
Plant hormone signal transduction
Protein processing in endoplasmic reticulum
Protein serine/threonine kinase activity
Regulation of transcription, DNA-templated
Sequence-specific DNA binding transcription factor activity
Starch and sucrose metabolism
Zinc ion binding
Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Ferrarini A. The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24(9):3489–505.
Food and agriculture organization of the united nations statistic devision Available at http://faostat3.fao.org/download/Q/QC/E. Accessed 28 Dec 2016.
Luedeling E, Brown PH. A global analysis of the comparability of winter chill models for fruit and nut trees. Int J Biometeorol. 2011;55(3):411–21.
Campoy JA, Ruiz D, Egea J. Dormancy in temperate fruit trees in a global warming context: a review. Sci Hortic. 2011;130(2):357–72.
George AP, Nissen RJ, Baker JA. Low chill peach and nectarine cultivars. Qld Agric J. 1986;112:27–33.
Erez A. Chemical control of budbreak. HortSci. 1987;22(6):1240–3.
Or E, Vilozny I, Eyal Y, Ogrodovitch A. The transduction of the signal for grape bud dormancy breaking induced by hydrogen cyanamide may involve the SNF-like protein kinase GDBRPK. Plant Mol Biol. 2000;43(4):483–94.
Vergara R, Pérez FJ. Similarities between natural and chemically induced bud-endodormancy release in grapevine Vitis vinifera L. Sci Hortic. 2010;125(4):648–53.
Sreekantan L, Mathiason K, Grimplet J, Schlauch K, Dickerson JA, Fennell AY. Differential floral development and gene expression in grapevines during long and short photoperiods suggests a role for floral genes in dormancy transitioning. Plant Mol Biol. 2010;73(1–2):191–205.
Díaz-Riquelme J, Grimplet J, Martínez-Zapater JM, Carmona MJ. Transcriptome variation along bud development in grapevine (Vitis vinifera L.). BMC Plant Biol. 2012;12(1):181.
Lang GA, Early JD, Martin GC, Darnell RL. Endo-, para-, and ecodormancy: physiological terminology and classification for dormancy research. HortSci. 1987;22(3):371–37.
Mathiason K, He D, Grimplet J, Venkateswari J, Galbraith DW, Or E, Fennell A. Transcript profiling in Vitis riparia during chilling requirement fulfillment reveals coordination of gene expression patterns with optimized bud break. Funct Integr Genomics. 2009;9(1):81–96.
Zhu Y, Li Y, Xin D, Chen W, Shao X, Wang Y, Guo W. RNA-Seq-based transcriptome analysis of dormant flower buds of Chinese cherry (Prunus pseudocerasus). Gene. 2015;555(2):362–76.
Jiménez S, Lawton-Rauh AL, Reighard GL, Abbott AG, Bielenberg DG. Phylogenetic analysis and molecular evolution of the dormancy associated MADS-box genes from peach. BMC Plant Biol. 2009;9(1):1.
Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, Akagi T, Tao R. Functional and expressional analyses of PmDAM genes associated with endodormancy in Japanese apricot. Plant Physiol. 2011;157(1):485–97.
Liu G, Li W, Zheng P, Xu T, Chen L, Liu D, Teng Y. Transcriptomic analysis of ‘Suli’pear (Pyrus pyrifolia white pear group) buds during the dormancy by RNA-Seq. BMC Genomics. 2012;13(1):1.
Bai S, Saito T, Sakamoto D, Ito A, Fujii H, and Moriguchi T. Transcriptome analysis of Japanese pear (Pyrus pyrifolia Nakai) flower buds transitioning through endodormancy. Plant Cell Physiol. 2013;pct067. doi:10.1093/pcp/pct067.
Annadurai RS, Jayakumar V, Mugasimangalam RC, Katta MA, Anand S, Gopinathan S, Rao SN. Next generation sequencing and de novo transcriptome analysis of Costus pictus D. Don, a non-model plant with potent anti-diabetic properties. BMC Genomics. 2012;13(1):663.
Habu T, Yamane H, Igarashi K. 454-pyrosequencing of the transcriptome in leaf and flower buds of Japanese apricot (Prunus mume Sieb. et Zucc.) at different dormant stages. Jpn Soc Hortic Sci. 2012;81(3):239–50.
Böhlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, Strauss SH, Nilsson O. CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science. 2006;312(5776):1040–3.
Olsen JE. Light and temperature sensing and signaling in induction of bud dormancy in woody plants. Plant Mol Biol. 2010;73(1–2):37–47.
Doğramacı M, Foley ME, Chao WS, Christoffers MJ, Anderson JV. Induction of endodormancy in crown buds of leafy spurge (Euphorbia esula L.) implicates a role for ethylene and cross-talk between photoperiod and temperature. Plant Mol Biol. 2013;81(6):577–93.
Ueno S, Klopp C, Leplé JC, Derory J, Noirot C, Léger V, Le Provost G. Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing. BMC Genomics. 2013;14(1):236.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Olsen JE. Mechanisms of dormancy regulation Proc, Xth IS on plant bioregulation in Frui. In: Webster AD, Ramírez H, editors. Acta Hort. 2006. p. 727–30.
Fennell A, Hoover E. Photoperiod influences growth, bud dormancy, and cold acclimation in Vitis labruscana and V. riparia. J Am Soc Hort Sci. 1991;116(2):270–3.
Kühn N, Ormeño-Núñez J, Jaque-Zamora G, Pérez FJ. Photoperiod modifies the diurnal expression profile of VvPHYA and VvPHYB transcripts in field-grown grapevine leaves. J Plant Physiol. 2009;166(11):1172–80.
Hedley PE, Russell JR, Jorgensen L, Gordon S, Morris JA, Hackett CA, Brennan R. Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrum L.). BMC Plant Biol. 2010;10(1):202.
Horvath DP, Chao WS, Suttle JC, Thimmapuram J, Anderson JV. Transcriptome analysis identifies novel responses and potential regulatory genes involved in seasonal dormancy transitions of leafy spurge (Euphorbia esula L.). BMC Genomics. 2008;9(1):536.
Horvath DP, Sung S, Kim D, Chao W, Anderson J. Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge. Plant Mol Biol. 2010;73(1–2):169–79.
Pérez FJ, Vergara R, Rubio S. H2O2 is involved in the dormancy-breaking effect of hydrogen cyanamide in grapevine buds. Plant Growth Regul. 2008;55(2):149–55.
Pang X, Halaly T, Crane O, Keilin T, Keren-Keiserman A, Ogrodovitch A, Or E. Involvement of calcium signalling in dormancy release of grape buds. J Exp Bot. 2007;58(12):3249–62.
Horvath DP. Microarray results corresponding to Geo Seriesaccession number GSE8849 link at davidhorvath.com.whois. Accessed 19 Oct 2016.
Arora R, Rowland LJ, Tanino K. Induction and release of bud dormancy in woody perennials: a science comes of age. HortSci. 2003;38(5):911–21.
Horvath DP, Anderson JV, Chao WS, Foley ME. Knowing when to grow: signals regulating bud dormancy. Trends Plant Sci. 2003;8(11):534–40.
Rohde A, Ruttink T, Hostyn V, Sterck L, Van Driessche K, Boerjan W. Gene expression during the induction, maintenance, and release of dormancy in apical buds of poplar. J Exp Bot. 2007;58(15–16):4047–60.
Tamura F, Tanabe K, Banno K. Effect of chilling treatment on intensity of bud dormancy, respiration and endogenousgrowth regulators in Japanese pear ‘Nijisseiki’. Jpn Soc Hortic Sci. 1992;60:763–9 (in Japanese with English abstract).
Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, et al. A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007;19:2370–90.
Voegele A, Linkies A, Müller K, Leubner-Metzger G. Members of the gibberellin receptor gene family GID1 (GIBBERELLIN INSENSITIVE DWARF1) play distinct roles during Lepidium sativum and Arabidopsis thaliana seed germination. J Exp Bot. 2011;62(14):5131–47.
Beveridge CA, Symons GM, Turnbull CG. Auxin inhibition of decapitation-induced branching is dependent on graft-transmissible signals regulated by genes Rms1 and Rms2. Plant Physiol. 2000;123(2):689–98.
Horvath D. Common mechanisms regulate flowering and dormancy. Plant Sci. 2009;177(6):523–31.
Dietz K, Vogel MO, Viehhauser A. AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma. 2010;245:3–14.
Yamada M, Yamane H, Sato A, Hirakawa N, Iwanami H, Yoshinaga K, Nakajima I. New grape cultivar ‘Shine Muscat’. Bull Natl Inst Fruit Tree Sci. 2008;7:21–38.
Balwierz PJ, Carninci P, Daub CO, Kawai J, Hayashizaki Y, Van Belle W, Beisel C. Van Nimwegen E Methods foranalyzing deep sequencing expression data: constructing the human and mouse promoterome with deepCAGE data. Genome Biol. 2009;10:R79. doi:10.1186/gb-2009-10-7-r79.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5((7):621–8.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Method Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.
Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):1–21.
Benjamini Y, and Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001:1165–1188.
Bland JM, Altman DG. Multiple significance tests: the Bonferroni method. BMJ. 1995;310(6973):170.
Schefe JH, Lehmann KE, Buschmann IR, Unger T, Funke-Kaiser H. Quantitative real-time RT-PCR data analysis: current concepts and the novel “gene expression’s C T difference” formula. J Mol Med. 2006;84(11):901–10.
We are thankful Dr. Syed Tahir Ata.ul.Karim and Dr. Muhammad Faheem for their technical assistance regarding critical review of manuscript.
The present work was supported by the China National ‘948’ key project [2011; G28] and China Agriculture Research System (CARS-30).
Availability of data and materials
All supporting data can be found within the manuscript and its additional files.
Conceived and design the experiment: JMT and MKR; Performed the experiments: MKR and LS; Writing of the manuscript: JMT and MKR. Analyzed the data: MKR. CL, MF and WW. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Differentially expressed genes between paradormancy vs endodormancy. (XLSX 417 kb)
Differentially expressed genes between summer buds vs endodormancy. (XLSX 448 kb)
Differentially expressed genes between summer buds vs paradormancy. (XLSX 187 kb)
Up and down regulated differentially expressed genes in cluster analysis. (XLSX 326 kb)
Differentially expressed genes involved in KEGG pathway between summer buds vs para dormancy. (XLSX 13 kb)
Differentially expressed genes involved in KEGG pathway between summer buds vs endo dormancy. (XLSX 13 kb)
Differentially expressed genes involved in KEGG pathways between paradormancy vs endodormancy. (XLSX 24 kb)
Genes and primer pairs used for quantitative real-time PCR. (DOCX 16 kb)