- Research article
- Open Access
Transcriptome analysis of seed dormancy after rinsing and chilling in ornamental peaches (Prunus persica (L.) Batsch)
© The Author(s). 2016
- Received: 20 January 2016
- Accepted: 28 July 2016
- Published: 8 August 2016
Ornamental peaches cv. ‘Yaguchi’ (Prunus persica (L.) Batsch) can be propagated via seeds. The establishment of efficient seed treatments for early germination and seedling growth is required to shorten nursery and breeding periods. It is important, therefore, to identify potential candidate genes responsible for the effects of rinsing and chilling on seed germination. We hypothesized that longer rinsing combined with chilling of seeds can alter the genes expression in related to dormancy and then raise the germination rate in the peach. To date, most molecular studies in peaches have involved structural genomics, and few transcriptome studies of seed germination have been conducted. In this study, we investigated the function of key seed dormancy-related genes using next-generation sequencing to profile the transcriptomes involved in seed dormancy in peaches. De novo assembly and analysis of the transcriptome identified differentially expressed and unique genes present in this fruit.
De novo RNA-sequencing of peach was performed using the Illumina Miseq 2000 system. Paired-end sequence from mRNAs generated high quality sequence reads (9,049,964, 10,026,362 and 10,101,918 reads) from ‘Yaguchi’ peach seeds before rinsed (BR) and after rinsed for 2 or 7 days with a chilling period of 4 weeks (termed 2D4W and 7D4W), respectively. The germination rate of 7D4W was significantly higher than that of 2D4W. In total, we obtained 51,366 unique sequences. Differential expression analysis identified 7752, 8469 and 506 differentially expressed genes from BR vs 2D4W, BR vs 7D4W and 2D4W vs 7D4W libraries respectively, filtered based on p-value and an adjusted false discovery rate of less than 0.05. This study identified genes associated with the rinsing and chilling process that included those associated with phytohormones, the stress response and transcription factors. 7D4W treatment downregulated genes involved in ABA synthesis, catabolism and signaling pathways, which eventually suppressed abscisic acid activity and consequently promoted germination and seedling growth. Stress response genes were also downregulated by the 7D4W treatment, suggesting that this treatment released seeds from endodormancy. Transcription factors were upregulated by the BR and 2D4W treatment, suggesting that they play important roles in maintaining seed dormancy.
This work indicated that longer rinsing combined with chilling affects gene expression and germination rate, and identified potential candidate genes responsible for dormancy progression in seeds of ‘Yaguchi’ peach. The results could be used to develop breeding programs and will aid future functional genomic research in peaches and other fruit trees.
- Seed germination
- Ornamental peach
Peach (Prunus persica (L.) Batsch) is a deciduous tree of the rose family (Rosaceae). Peaches are not only widely planted as fruit trees, but also as ornamental plants such as garden trees and cut branches. Budding and grafting are common asexual propagation techniques for fruit trees except breeding purpose because fruit phenotype cannot be maintained via seed propagation. However, seed propagation is possible for ornamental peaches such as ‘Yaguchi’ and ‘Hokimomo’ because they have been repeatedly propagated via seed for a long time and thus genetically homogeneous. Early germination is required to shorten nursery and breeding periods.
Seed germination depends on a number of factors, including internal dormancy and the environment . Seed dormancy is observed in higher plants with distinct physiological and morphological characteristics in different species ; however, it is an undesirable characteristic in an agricultural crop where rapid germination and growth are required. Moist chilling plays an important role in providing the stimulus required to overcome dormancy, increase germination and produce normal seedlings in ‘GF305’ peaches , strawberry trees  and sweet cherries (Prunus avium) . Seeds of most Prunus species require a period of chilling to break seed dormancy [6–9]. In peaches, removal of the seed coat shortens the chilling periods needed to break dormancy and even improves the germination of non-chilled seeds [10, 11]. When the cold treatment is insufficient, seedlings show physiological dwarfing, which is considered a special case of embryo dormancy [12, 13]. These results implied that dormancy in peach seeds is caused by exogenous and endogenous dormancy associated with the seed covering layers and the embryo .
The ratio of the hormones abscisic acid (ABA) and gibberellic acid (GA) is considered a relevant factor regulating seed dormancy. Moist chilling induced an increase in GA levels in embryos of the European hazel (Corylus avellana), suggesting that GAs synthesized during cold treatment were responsible for breaking dormancy . ABA plays a key role in various aspects of plant growth. In the ABA biosynthetic pathway, 9-cis-epoxycarotenoid dioxygenase (NCED) is the key enzyme in ABA biosynthesis in higher plants . The pattern of PacNCED1 expression was coincident with that of ABA accumulation in sweet cherry fruit . Recent studies indicated that the key step of ABA inactivation is the hydroxylation of the 8′-methyl group of ABA in most plant tissues. ABA 8′-Hydroxylase is a key enzyme in the oxidative catabolism of ABA and is expressed throughout sweet cherry fruit development.
Germination commences with the uptake of water by imbibition by the dry seed, followed by embryo expansion. In our previous study , seeds after rinsing with running tap water for 2 days and chilling at 5 °C for more than 8 weeks showed decreased ABA contents in the embryonic axis and seed coat, which ultimately increased the uniform germination and normal growth in ‘Yaguchi’ peaches. Moreover, longer rinsing (about 8 days) increased germination rate and significantly increased the plant height in ‘Hokimomo’ peach . In Vitis vinifera, rinsing with running water for 8 days increased germination rate from 14 to 34 % compared to 4 days .
Therefore, we hypothesized that longer rinsing of seeds can alter the genes expression in related to dormancy and then raise the germination rate in ‘Yaguchi’ peach. In the present study, germination rate of seeds after rinsed for 2 days (2D) and 7 days (7D) and chilled for 4 weeks each (termed 2D4W and 7D4W, respectively) were determined. The aim of the present study was to perform transcriptome analysis among seeds of BR, 2D4W and 7D4W seeds and to gain an understanding of molecular mechanism during peach seed dormancy and to prove the hypothesis.
In recent years, RNA-sequencing (RNA-Seq) technology has become the most popular and powerful tool for transcriptome analysis. RNA-Seq is cheaper, more efficient, more sensitive and more accurate in generating transcriptome profiles than microarray analysis and other techniques [21–28]. We use RNA-Seq technology to identify and characterize the expression of the large number of genes, especially those differentially expressed during dormancy progression.
The germination rate of embryo in peach
De novo transcriptome assembly and assessments of expressed sequenced tags
Sequencing and de novo assembly
The numbers of reads
Raw NGS data
Quality filtered data
Mean quality value
Mean quality value
Transcriptome of BR (PE reads)
Transcriptome of 2D4W (PE reads)
Transcriptome of 7D4W (PE reads)
The summary of assembly
Min sequence length (bp)
Max sequence length (bp)
Average sequence length (bp)
N50 length (bp)
(G + C)s
Differentially expressed genes analysis (DEGs)
We compared the gene expressions among BR, 2D4W and 7D4W data and identified differentially expressed genes. Normalization was applied to the treatments to provide accurate differential expression rather than individual quantification. Total DEGs patterns for each sample are presented in Additional file 3. The significantly DEGs with FDRs < 0.05 obtained by comparing BR vs 2D4W, BR vs 7D4W and 2D4W vs 7D4W libraries were 7752, 8469 and 506 respectively (see Additional file 4). The distributions of DEGs were further analyzed using edgeR package. DEGs are visualized as an MA plot (log ratio versus abundance plot) for set of experiments (BR vs 2D4W, BR vs 7D4W and 2D4W vs 7D4W, see Additional file 5). The red dots highlight transcripts of positive and negative values of log2 Fold Change (logFC), indicating that the sequences were upregulated and downregulated genes in each treatment. The genes names, fold changes, and p values for up- and down-regulated DEGs in each treatment were listed in Additional file 4.
Functional annotation of DEGs
For annotation, the consensus sequences were first searched against the Swiss Institute of Bioinformatics databases (SwissProt) database  by local BLASTX (E-value cut-off was set at 1e-5) to search for the maximum number of similar genes. Among the 7752 DEGs in BR vs 2D4W, 3121 and 4631 were up- and down-regulated genes. From 8469 DEGs of BR vs 7D4W, 3290 and 5179 were up- and down-regulated. 377 and 129 DEGs were up- and down-regulated genes in 2D4W vs 7D4W (see Additional file 4). These genes were identified from BLAST nr database, SwissProt and Uniprot database and assigned with Kyoto Encyclopedia of Genes and Genomes (KEGG)  and Gene ontology (GO) terms in the biological process, molecular functions and cellular component categories .
Functional classification by KEGG
KEGG pathway database records networks of molecular interactions in cells, and variants of them, specific to particular organisms. Pathway-based analysis helps to understand the biological functions of gene products. Pathway information for all annotated sequences was obtained from KEGG pathway annotations. In total, 7752 DEGs of BR vs 2D4W, 8469 DEGs of BR vs 7D4W and 506 DEGs in 2D4W vs 7D4W, 2789, 2982 and 142 sequences, respectively were mapped to the reference pathway in KEGG (see Additional file 4). We focused our attention on carotenoid biosynthesis (see Additional file 6). In the pathways of carotenoid biosynthesis and plant hormone transduction, NCED (EC: 188.8.131.52), ABA 8′-hydroxylase (EC: 184.108.40.206), abscisic acid receptor PYR/PYL family (PYR/PYL), serine/threonine-protein kinase SRK2 (SNRK2), protein phosphatase 2C (PP2C) and ABA responsive element binding factor (ABF) were known to have important roles in the seed dormancy and germination .
GO enrichment analysis of DEGs
Quantitative real-time reverse transcription PCR (qRT-PCR) validation of DEGs from transcriptome analysis
Genes associated with phytohormones
Phytohormones, such as ABA and GA, play an important role in controlling seed dormancy and germination [35, 37, 38], and mediate responses to both biotic and abiotic stresses . ABA is a positive regulator of dormancy and a negative regulator of germination, whereas GA promotes seed germination . In particular, ABA has been termed as ‘dorman’ or ‘dormancy inductor’ .
In this study, pathway analysis of KEGG analysis showed upregulation in 2D4W of the unigenes encoding proteins involved in ABA biosynthesis and catabolism, including NCED1, ABA 8′-hydroxylase in carotenoid biosynthesis and PP2C in plant hormone transduction pathways (see Additional file 6). While only NCED1 and PP2C showed upreguration in BR. NCED is the key enzyme in ABA biosynthesis in higher plants , while ABA 8′-hydroxylase is the key enzyme in the oxidative catabolism of ABA and is expressed throughout fruit development . The expression of NCED1 is closely associated with the ABA level in imbibing dormant seeds of Brachypodium distachyon . Our data suggested that these genes could be involved in modulation of embryo dormancy through ABA synthesis and catabolism.
ABA signaling is important for stress responses, seed development and dormancy [42, 43]. The activity of PP2C is controlled by ABA receptors, which upon binding to ABA form a complex with, and inactivate PP2C. Inactivation of PP2C leads to upregulation of SnRK2, which activates downstream transcription factors, including the ABA responsive element (ABRE) binding factor (ABF), ABA insensitive5 (ABI5), ABI3 and ABI4, and thereby mediates seed responsiveness to ABA . Twenty-one DEGs were annotated as probable PP2C and were upregulated in the endodormancy state in flower bud of Japanese pear . In sorghum, the higher expression of an ABI5 homolog has been associated with higher levels of dormancy . In this study, KEGG pathway analysis showed that transcripts of a gene encoding PP2C were significantly higher in BR and 2D4W. Similarly, the qRT-PCR results showed upregulation of PP2C in BR and 2D4W and dramatic downregulation in 7D4W (Fig. 3h). In addition, ABI5-binding protein genes were also significantly upregulated in 2D4W and significantly downregulated in 7D4W (Fig. 3e and f).
These data suggested that the ABA signaling pathway was still activate in dry seed and 2D4W; ie short period rinsing, and it might be significant for maintaining ABA sensitivity and endodormancy of the embryo. Transcripts of genes encoding GA2-oxidase (GA2-ox) 8, which inactivates GA1, were significantly higher in 2D4W. GA catabolic genes such as GA2ox are important to control GA levels [47–50]. These results suggested that the ABA signaling pathway and GA inactivating genes play an important role in dormancy and germination and were downregulated by the longer rinsed treatment in ‘Yaguchi’ peach seeds.
Genes associated with stress response
Seed development and maturation are accompanied by increased desiccation tolerance . Seeds activate a series of mechanisms to respond to many biotic and abiotic stresses during germination when placed in water or their external environment changes. The EID1-like protein 3 (EDL3) shows high similarity to EID1 (Empfindlicher im dunkelroten Licht 1), an F-box protein . EID1 functions as a positive regulator in ABA-dependent signaling cascades that control seed germination, root growth, greening of etiolated seedlings and transition to flowering in Arabidopsis thaliana . EDL3 expression was induced under osmotic stress, high salinity, and ABA application . Our results showed that rinsing for 7D was effective to decrease the expression of EID1 gene both at 7D and 7D4W (Fig. 3a), indicating that longer rinsing decreased EID1 expression.
Dehydration-responsive element-binding proteins (DREBs) are transcriptional regulators of the APETALA2/Ethylene Responsive element-binding Factor (AP2/ERF) family that control the expressions of abiotic stress-related genes. Under conditions of mild heat stress, transgenic seeds overexpressing DREB2C showed delayed germination and increased ABA content compared with the untransformed wild-type in Arabidopsis . In this study, the DREB2CA loci were upregulated in BR, 2D and 2D4W (Fig. 3b), but downregulated in 7D4W, suggesting that 7D4W released seeds from endodormancy.
Furthermore, these results showed that transcripts of genes encoding the LEA D-34 protein or dehydrin (related to drought stress) there was higher expression of the genes encoding the LEA D-34 in BR and then downregulated after chilling (see Additional file 4 and Fig. 3c). LEA proteins are classified into at least five groups by the similarity in their amino acid sequences [55, 56] and are associated with dormancy transition . In the Norway spruce (Betula pubescens Ehrh.), the expression of certain dehydrin genes gradually decreased when approaching bud burst . In dormant Japanese pears, LEA genes are significantly higher in the deepest endodormancy period . Additionally, dehydrins are synthesized by cells in response to ABA  and most PmLEA genes are upregulated by ABA treatment . These results indicated that the LEA gene was induced to develop tolerance against drought in dry seeds without rinsing.
Genes associated with transcription factors (TFs)
Transcription factors (TFs) play a crucial role in plant development and stress response . One group of well-studied transcription factors involved in chilling responses is the AP2-EREBP family members, which have been subdivided into four major subfamilies: the DREB/CBF, AP2, ERF subfamilies . Of these, the DREB/CBF subfamily has been reported to play a major role in the early stages of the chilling response , as evidenced by studies in numerous species such as Anthurium  and rice , as well as tea plant [66, 67]. In the GO molecular function and biological process categories, sequence-specific DNA binding TF activity and regulation of transcription were significantly enriched in 2D4W compared with 7D4W (Fig. 2). In other words, 7D4W, ie, a longer rinsing treatment, downregulated the expression of several TFs, such as heat stress TF, AP2-like ethylene-responsive TF, DREB2C, ethylene-responsive TF (ERF TF) and the MYB TF family (see Additional file 4). These results suggested that genes of sequence-specific DNA binding TF activity are associated with the effects of the longer rinsing combined with chilling on breaking dormancy in ‘Yaguchi’ seeds.
MADS-box TFs play fundamental roles in plant development, such as floral organ and meristems identity determination and transition from vegetative to reproductive growth regulation. MADS-box TFs are involved in bud dormancy regulation in peaches and other species [68–71]. PpMADS13-1 and −2, isolated from the leaf bud of Japanese pear, appear to be upregulated towards endodormancy establishment and downregulated concomitant with endodormancy release [72, 73], similar to the peach DAM5 and DAM6 genes . In this study, chilling treatment effectively reduced the expression of LeMADS (Fig. 3j). Similarly, the expressions of DAM1 and DAM6 were decreased significantly by chilling treatment in peach seeds . These results indicated that MADS-box genes maintain seed dormancy and are repressed when seed dormancy is released by chilling.
This study indicated that longer rinsing and chilling affects gene expression and germination rate, and investigated key seed dormancy-related genes using next-generation sequencing to profile the transcriptomes of dormant ‘Yaguchi’ peach seeds. We identified genes associated with the effects of the rinsing and chilling process, including genes associated with phytohormones, the stress response and TFs. The 7D4W treatment downregulated genes involved in ABA synthesis; NCED1, catabolism; ABA 8′-hydroxylase, ABA signaling pathway; PP2C and ABI5 and GA inactivation; GA2-ox8, which eventually suppressed ABA activity and consequently promoting germination and seedling growth. Stress response genes (EID1, DREB2CA and LEA D-34) were also downregulated by 7D4W treatment, suggesting that it released seeds from endodormancy. TF genes such as AP2-like ethylene-responsive, DREB2Cs, ERF and LeMADS were upregulated in BR and 2D4W, suggesting these TFs play important roles in maintaining seed dormancy.
Fully ripened ‘Yaguchi’ peaches were collected from the peach flower garden of Koga city, Ibaraki prefecture, Japan, in middle October 2013 (36°10'N, 139°42'E). Voucher specimens were gathered from the leaf and seed specimens (Voucher specimen accession number: PPKYL0716-1; PPKYS0716-1) and identified by M. Yoshida and K. Yamane. After the skin and flesh were removed, endocarps were carefully cracked in a vice to remove the seeds. We hypothesized that longer rinsing combined with chilling of seeds can alter the genes expression in related to dormancy and then raise the germination rate in this peach. To prove this assumption, Seeds were rinsed for 2 or 7 days under running tap water and transferred into petri dishes containing two-layers of filter paper moistened with sterile water. The plastic petri dishes were sealed by parafilm and kept in the refrigerator at 5 °C for 2, 4, 6, and 8 weeks. The germination and height of seedlings were noted after 35 days. Three biological replicates of 10 embryonic axes of dry seed before and after rinsing and chilling were immediately frozen with liquid nitrogen and stored at −80 °C. For RNA extraction and qRT-PCR, BR; 2D; 7D; 2D4W and 7D4W were used. For library preparation, the sample from only BR, 2D4W and 7D4W were constructed due to these treatments were significantly different of germination rate (Fig. 1). Therefore, we would like to know how to difference of the differentially expression genes among these treatments.
Samples of 10 mg of frozen plant tissue were ground in liquid nitrogen to a fine powder. Total RNA was extracted from these tissues using the method described by . Ten embryonic axes were extracted with three replications. RNA integrity was evaluated with a 1.0 % agarose gel stained with ethidium bromide. Then, total RNA were quantified and examined for protein contamination (A260/A280) and reagent contamination (A260/A230) by using a NanoDrop ND-1000 spectrophotometer. The RNA integrity number determined by the Agilent Technologies 2100 Bioanalyzer (Agilent Technologies) together with high sensitivity DNA Lab chip kit . High quality total RNA samples give two distinct peaks and yield an RNA Integrity Number (RIN) value greater than 8. Total RNA (200 ng) from each biological triplicate was used for library preparation and sequenced.
Library preparation and Illumina sequencing
The libraries for RNA-Seq were constructed using a SureSelect Strand-Specific RNA Library prep kit (Agilent Technologies), according to manufacturer's instructions. Oligo-(dT) magnetic beads were used to isolate poly-(A) mRNA from total RNA, and mRNA was fragmented in fragmentation buffer. Using these short fragments as templates, random hexamer primers were used to synthesize first-strand cDNA. Second-strand cDNA was synthesized using buffer, dNTPs, RNaseH, and DNA polymerase I. The product will be double-standed cDNA (ds cDNA). ds cDNA fragments were purified from free nucleotides, enzymes, buffers, and RNA with a AMPure XP beads (Beckman Coulter Genomics). Perform end-repair on purified eluted ds cDNA and adding poly (A), then ligated to sequencing adapters. Adaptors were ligating to both ends of the ds cDNA. After purification via AMPure XP beads again, suitable fragments were enriched by polymerase chain reaction (PCR) amplification and index the adaptor-ligated cDNA library. Each adaptor had eight-nucleotide difference in adaptor sequenced. A different index for each library reaction was use allowed for pooling libraries later for sequenced. Finally, sent normalized and pooled libraries to sequencing facility for cluster generation on a MiSeq Sequencer (Illumina, USA), the paired-end library was prepared following the protocol of the Illumina MiSeq Reagent Kits v3 (150 cycles, paired-end). All three replicates of the BR, 2D4W and 7D4W samples were sequenced.
Pre-processing of raw short read sequences and de novo assembly
Quality of sequence reads was checked with FastQC v0.11.5. Raw reads were first processed using cutadapt version 1.8.1  to remove adapter sequences, low quality ends (quality scores < 30) and the last 76th base, while reads shorter than 50-bp were discarded. The high-quality short reads were assembled into unigenes using the Trinity program version: trinityrnaseq-2.0.6 . rRNA sequences were excluded from the unigenes by removing sequences matched to the SILVA rRNA database by the megablast program . The high-quality short reads were mapped to the rRNA-removed unigenes using Bowtie . Then, to estimate the expression levels of the transcripts, the number of uniquely mapped reads for each unigene was competed and normalized to RPKM values (reads per kilobase per million mapped reads) using RSEM method . The RPKM values were then compared pairwise as: BR/2D4W, BR/7D4W, and 2D4W/7D4W. To identify DEGs, the edgeR package  was used to compute the p-value and fold change. The p-value was used to identify genes expressed differentially between the paired treatments. Significantly DGEs were identified using a FDR (false discovery rate) threshold of ≤ 0.05 and a minimum two-fold change.
Validation of de novo assembly
To verify the accuracy of de novo assembly and annotation, we mapped reads on our assemblies and the Prunus_persica_v2.0.a1.all transcripts . Mapping statistics were obtained using Bowtie2 (V.2.2.5) package.
To predict the biological functions of unigenes, sequences were used for BLAST searches and annotation against an NCBI nr protein database (NCBI non-redundant sequence database). Consensus sequences were further aligned by BLASTX to protein databases such as Swiss-Prot (E-value cut-off was set at 1e-5) , KEGG and GO. KEGG pathway mapping for unigenes were done using the KEGG automatic annotation server (KAAS) . Putative coding regions were extracted from Trinity transcripts using TransDecoder . GO annotation of the transcriptome was performed using InterproScan .
Quantitative gene expression analysis
To validate the mRNA abundance of 10 genes that were significantly regulated by rinsing and chilling during the RNA-seq analysis, qRT-PCR was performed. Total RNA was isolated from embryonic axes and cDNAs were generated from the RNA samples as those used for the RNA-seq experiment. For each sample, 1 μg of total RNA was used with the QuantiTect® SYBR® Green RT-PCR Kit, according to the manufacturer’s protocol (Agilent). The cDNA was stored at −80 °C until use in the qRT-PCR analysis. The qRT-PCR assay was performed with three technical replicates using the Fast SYBR Green Master Mix (Life Technologies, Roche) on a Lightcycler 96 Real Time PCR system in a total volume of 15 μL. The PCR cycle comprised one 600 s cycle at 95 °C, followed by 45 cycles at 94 °C for 15 s, 60 °C for 30 s and 72 °C for 30 s. All amplified products were subjected to melt curve analysis. A negative control without a cDNA template was run with all analyses to evaluate the overall specificity. The reference gene ubiquitin was used to normalize the total amount of cDNA in each reaction. This gene was stably expressed in peach . Amplification efficiency and relative gene expression levels were calculated using the ΔΔC T and 2-ΔΔC T methods (C T; cycle threshold). The ΔC T value of each gene was calculated by subtracting the C T value of the endogenous control from the C T value of the target gene. Gene-specific primers were designed using primer-BLAST  (see Additional file 11).
ABA, Abscisic acid; DEGs, Differentially expressed genes; FC, Fold Change; FDR, False Discovery Rate; GA, Gibberelin; GO, Gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PE, Paired-end; qRT-PCR, quantitative real-time reverse transcription PCR; RNA-seq, RNA-sequencing; TF, Transcription Factor
The authors would like to thank Assoc. Prof. Xiaonan Xie for assistance with the ABA content analysis and Kitamura Rika for assistance with the RNA-seq analysis. This research was supported by United Graduate School of Agricultural Science, Tokyo University of Agriculture and Technology; Bioscience Education and Research Center, Utsunomiya University; Faculty of Agriculture, Utsunomiya University; Genome human resource development program and Japanese Government (Monbukagakusho) Scholarship.
Budget of the United Graduate School of Agricultural Science, Tokyo University of Agriculture and Technology.
Availability of data and materials
The data sets supporting the results of this article are available in the DDBJ Sequence Read Archive (DRA) (accession number DRA004674).
WK carried out all experiments and drafted the manuscript. WK and KY conceived the study, drafted and revised the manuscript. MI provides machine, TS, KI and TK gave advice on the experiments, data analysis. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Bewley JD, Black M. Seeds: Physiology of development and germination. NY: Plenum Press; 1994.View ArticleGoogle Scholar
- Finch-Savage WE, Leubner-Metzger G. Seed dormancy and the control of germination. New Phytol. 2006;171:501–23.View ArticlePubMedGoogle Scholar
- Martınez-Gomez P, Dicenta F. Mechanisms of dormancy in seeds of peach (Prunus persica (L.) Batsch) cv. ‘GF305’. Sci Hort. 2001;91:51–8.View ArticleGoogle Scholar
- Karam NS, AL-Salem MM. Breaking dormancy in Arbutus andrachna L. seeds by stratification and gibberellic acid. Seed Sci Technol. 2001;29:51–6.Google Scholar
- Jensen M, Eriksen EN. Development of primary dormancy in seeds of Prunus avium during maturation. Seed Sci Technol. 2001;29:307–20.Google Scholar
- Crocker W, Barton LV. After-ripening, germination and storage of certain rosaceous seeds. Contrib Boyce Thompson Inst. 1931;3(3):385–404.Google Scholar
- Grisez TJ. Prunus L. Cherry, peach, and plum: Seeds of woody plants in the United States.Handbook No.450. In: Schopmeyer CS, editor. Chapter 1 Seed biology. Washington DC: United States Department of Agriculture Forest Service; 1974. p. 658–673.Google Scholar
- Chang S, Werner DJ. Relationship of seed germination and respiration during stratification with cultivar chilling requirement in peach. Rev Neurol. 1983;139(10):42–5.Google Scholar
- Suszka B, Muller C, Bonnet-Masimbert M. Seeds of broadleaves from harvest to sowing. Paris: Institut National de la Recherche Agronomique (INRA); 1996.Google Scholar
- Tukey HB, Barrett MS. Approximate germination test for non-after-ripened peach seed. Plant Physiol. 1936;11:629–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Mehanna HT, Martin GC, Nishijima C. Effects of temperature, chemical treatments and endogenous hormone content on peach seed germination and subsequent seedling growth. Sci Hort. 1985;27:63–73.View ArticleGoogle Scholar
- Hartmann HT, Kester DE, Davies FT, Geneve RL. Principles and practices, 6th edition: Plant propagation. USA: Prentice Hall Inc.; 1997.Google Scholar
- Bewley JD, Black M. Physiology and biochemistry of seeds in relation to germination. 2nd ed. Berlin: Springer; 1982. p. 60–125.View ArticleGoogle Scholar
- Diaz DH, Martin GC. Peach seed dormancy in relation to endogenous inhibitors and applied growth substances. J Am Soc Hortic Sci. 1972;97:652–4.Google Scholar
- Williams PM, Bradbeer JW, Gaskin P, MacMillan J. Studies in seed dormancy. VIII. The identification and determination of gibberellins A1 and A9 in seeds of Corylus avellana L. Planta. 1974;117:101–8.View ArticlePubMedGoogle Scholar
- Tan BC, Schwartz SH, Zeevaart JAD, McCarty DR. Genetic control of abscisic acid biosynthesis in maize. Proc Natl Acad Sci U S A. 1997;94(22):12235–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Ren J, Sun L, Wu J, Zhao S, Wang C, Wang Y, et al. Cloning and expression analysis of cDNAs for ABA 8'-hydroxylase during sweet cherry fruit maturation and under stress conditions. J Plant Physiol. 2010;167(17):1486–93.View ArticlePubMedGoogle Scholar
- Pawasut A, Yamane K, Fujishige N, Yoneyama K, Yamaki Y, Honjo H. Influence of seed coat removal and chilling on abscisic acid content and germination in ornamental peach (Prunus persica Batsch) seeds. Journal of Horticultural Science & Biotechnology. 2010;85(3):248–52.View ArticleGoogle Scholar
- Yamane K, Burana C, Kisugi T, Xia X, Yoneyama K, Yamaki Y. Effects of rinsing and low temperature periods on the ABA contents, germination and seedling growth in ornamental peaches. Hort Res (Japan). 2011;10 Suppl 2:262. In Japanese.Google Scholar
- Kachru RB, Singh RN, Yadev IS. Physiological studies on dormancy in grape seeds (Vitis vinifera var. Black Muscat; II. On the effect of exogenous application of growth substances, low chilling temperature and subjection of the seeds to running water. Vitis. 1972;11:289–95.Google Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomic. Nat Rev Genet. 2009;10(1):57–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, et al. Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012;13(1):19.View ArticlePubMedPubMed CentralGoogle Scholar
- Xue J, Bao YY, Li B, Cheng YB, Peng ZY. Transcriptome analysis of the brown planthopper Nilaparvata lugens. PLoS One. 2010;5(12):e14233.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez-Gomez P, Crisosto CH, Bonghi C, Rubio M. New approaches to Prunus transcriptome analysis. Genetica. 2011;139(6):755–69.View ArticlePubMedGoogle Scholar
- Tang Q, Ma XJ, Mo CM, Wilson IW, Song C, Zhao H, et al. An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011;12:343.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu K, Xu Q, Da X, Guo F, Ding Y, Deng X. Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis). BMC Genomics. 2012;13(1):10.View ArticlePubMedPubMed CentralGoogle Scholar
- Schuster SC. Next-generation sequencing transforms today’s biology. Nat Methods. 2008;5:16–8.View ArticlePubMedGoogle Scholar
- Shendure J, Ji H. Next-generation DNA sequencing. Nat Biotechnol. 2008;26(10):1135–45.View ArticlePubMedGoogle Scholar
- Verde I, Abbott AG, Scalabrin S, Jung S, Shu S, et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013;45:487–94.View ArticlePubMedGoogle Scholar
- Alkio M, Jonas U, Declercq M, Van Nocker S, Knoche M. Transcriptional dynamics of the developing sweet cherry (Prunus avium L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes. Horticulture Res. 2014;1:11.View ArticleGoogle Scholar
- Schliesky S, Gowik U, Weber AP, Bräutigam A. RNA-seq assembly–are we there yet? Frontiers in plant science 3. 2012; doi: 10.3389/fpls.2012.00220.
- Bairoch A, Boeckmann B. The SWISS-PROT protein sequence data bank. Nucleic Acids Res. 1991;19 Suppl: 2247.Google Scholar
- Kanehisa MGS. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Consortium GO. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32 suppl 1:D258–61.View ArticleGoogle Scholar
- Kucera B, Cohn MA, Leubner-Metzger G. Plant hormone interactions during seed dormancy release and germination. Seed Sci Res. 2005;15:281–307.View ArticleGoogle Scholar
- Yano R, Takebayashi Y, Nambara E, Kamiya Y, Seo M. Combining association mapping and transcriptomics identify HD2B histone deacetylase as a genetic factor associated with seed dormancy in Arabidopsis thaliana. Plant J. 2013;74:815–28.View ArticlePubMedGoogle Scholar
- Finkelstein R, Reeves W, Ariizumi T, Steber C. Molecular aspects of seed dormancy. Annu Rev Plant Biol. 2008;59:387–415.View ArticlePubMedGoogle Scholar
- Nonogaki M, Sall K, Nambara E, Nonogaki H. Amplification of ABA biosynthesis and signaling through a positive feedback mechanism in seeds. Plant J. 2014;78(3):527–39.View ArticlePubMedGoogle Scholar
- Santner A, Estelle M. Recent advances and emerging trends in plant hormone signaling. Nature. 2009;459:1071–8.View ArticlePubMedGoogle Scholar
- Takemura Y, Sudo S, Ikeda T, Matsumoto K, Tamura F. Chilling induces bud endodormancy in Japanese pear ‘Gold Nijisseiki’. Hortic Res. 2011;10:87–92.View ArticleGoogle Scholar
- Barrero JM, Jacobsen JV, Talbot MJ, White RG, Swain SM, Garvin DF, Gubler F. Grain dormancy and light quality effects on germination in the model grass Brachypodium distachyon. New Phytol. 2012;193:376–86.View ArticlePubMedGoogle Scholar
- Hubbard KE, Nishimura N, Hitomi K, Getzoff ED, Schroeder JI. Early abscisic acid signal transduction mechanisms: newly discovered components and newly emerging questions. Genes Dev. 2010;24:1695–708.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakashima K, Fujita Y, Kanamori N, Katagiri T, Umezawa T, Kidokoro S, et al. Three Arabidopsis SnRK2 protein kinases, SRK2D/SnRK2. 2, SRK2E/SnRK2. 6/OST1 and SRK2I/SnRK2. 3, involved in ABA signaling are essential for the control of seed development and dormancy. Plant Cell Physiol. 2009;50:1345–63.View ArticlePubMedGoogle Scholar
- Nambara E, Okamoto M, Tatematsu K, Yano R, Seo M, et al. Abscisic acid and the control of seed dormancy and germination. Seed Sci Res. 2010;20:55–67.View ArticleGoogle Scholar
- Takemura Y, Kuroki K, Shida Y, Araki S, Takeuchi Y, Tanaka K, et al. Comparative transcriptome analysis of the less-dormant Taiwanese pear and the dormant Japanese pear during winter season. PLoS One. 2015;10(10):e0139595.View ArticlePubMedPubMed CentralGoogle Scholar
- Rodrı’guez MV, Mendiondo GM, Maskin L, Gudesblat GE, Iusem ND, et al. Expression of ABA signalling genes and ABI5 protein levels in imbibed Sorghum bicolor caryopses with contrasting dormancy and at different developmental stages. Ann Bot. 2009;104:975–85.View ArticleGoogle Scholar
- Reid JB, Ross JJ, Swain SM. Internode length in Pisum: a new, slender mutant with elevated levels of C19 gibberellins. Planta. 1992;188:462–7.View ArticlePubMedGoogle Scholar
- Ross JJ, Reid JB, Swain SM, Hasan O, Poole AT, Hedden P, et al. Genetic regulation of gibberellin deactivation in Pisum. Plant J. 1995;17:241–50.Google Scholar
- Thomas SG, Phillips AL, Hedden P. Molecular cloning and functional expression of gibberellin 2-oxidases, multifunctional enzymes involved in gibberellin deactivation. Proc Natl Acad Sci U S A. 1999;96:4698–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Sakamoto T, Kobayashi M, Itoh H, Tagiri A, Kayano T, Tanaka H, et al. Expression of a gibberellin 2-oxidase gene around the shoot apex is related to phase transition in rice. Plant Physiol. 2001;125:1508–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Finnie C, Melchior S, Roepstorff P, Svensson B. Proteome analysis of grain filling and seed maturation in barley. Plant Physiol. 2002;129(3):1308–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Marrocco K, Zhou Y, Bury E, Dieterle M, Funk M, Genschik P, et al. Functional analysis of EID1, an F-box protein involved in phytochrome A-dependent light signal transduction. Plant J. 2006;45(3):423–38.View ArticlePubMedGoogle Scholar
- Koops P, Pelser S, Ignatz M, Klose C, Marrocco-Selden K, Kretsch T. EDL3 is an F-box protein involved in the regulation of abscisic acid signalling in Arabidopsis thaliana. J Exp Bot. 2011;62(15):5547–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Je J, Chen H, Song C, Lim CO. Arabidopsis DREB2C modulates ABA biosynthesis during germination. Biochem Biophys Res Commun. 2014;452(1):91–8.View ArticlePubMedGoogle Scholar
- Cuming AC. Seed Proteins: LEA proteins. In: Shewry PR, Casey R, editors. Chapter 1 Seed protein. Dordrecht, Netherlands: Kluwer Academic Press. 1999. p. 753–780.Google Scholar
- Wise J. LEAping to conclusions: a computational reanalysis of the late embryogenesis abundant proteins and their possible roles. BMC Bioinformatics. 2003;4:52–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Ophir R, Pang X, Halaly T, Venkateswari J, Lavee S, Galbraith D, et al. Gene-expression profiling of grape bud response to two alternative dormancy-release stimuli expose possible links between impaired mitochondrial activity, hypoxia, ethylene-ABA interplay and cell enlargement. Plant Mol Biol. 2009;71:403–23.View ArticlePubMedGoogle Scholar
- Yakovlev IA, Asante DKA, Fossdal CG, Partanen J, Junttila O, Johnsen O. Dehydrins expression related to timing of bud burst in Norway spruce. Planta. 2008;228:459–72.View ArticlePubMedGoogle Scholar
- Chandler PM, Robertson M. Gene expression regulated by abscisic acid and its relation to stress tolerance. Annu Rev Plant Physiol Plant Mol Biol. 1994;47:113–41.View ArticleGoogle Scholar
- Du D, Zhang Q, Cheng T, Pan H, Yang W, Sun L. Genome-wide identification and analysis of late embryogenesis abundant (LEA) genes in Prunus mume. Mol Biol Rep. 2013;40(2):1937–46.View ArticlePubMedGoogle Scholar
- Lee BH, Henderson DA, Zhu JK. The Arabidopsis cold-responsive transcriptome and its regulation by ICE1. Plant Cell. 2005;17:3155–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Gilmour SJ, Zarka DG, Stockinger EJ, Salazar MP, Houghton JM, Thomashow MF. Low temperature regulation of the Arabidopsis CBF family of AP2 transcriptional activators as an early step in cold-induced COR gene expression. Plant J. 1998;16:433–42.View ArticlePubMedGoogle Scholar
- Cook D, Fowler S, Fiehn O, Thomashow MF. A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis. Proc Natl Acad Sci U S A. 2004;101:15243–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Tian DQ, Pan XY, Yu YM, Wang WY, Zhang F, Ge YY, et al. De novo characterization of the Anthurium transcriptome and analysis of its digital gene expression under cold stress. BMC Genomics. 2013;14:827.View ArticlePubMedPubMed CentralGoogle Scholar
- Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, et al. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought, high‐salt‐and cold‐responsive gene expression. Plant J. 2003;33:751–63.View ArticlePubMedGoogle Scholar
- Wang Y, Jiang CJ, Li YY, Wei CL, Deng WW. CsICE1 and CsCBF1: two transcription factors involved in cold responses in Camellia sinensis. Plant Cell Rep. 2012;31:27–34.View ArticlePubMedGoogle Scholar
- Zheng C, Zhao L, Wang Y, Shen J, Zhang Y, Jia S, et al. Integrated RNA-Seq and sRNASeq analysis identifies chilling and freezing responsive key molecular players and pathways in tea plant (Camellia sinensis). PLoS One. 2015;10(4):e0125031.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:536.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:169–79.View ArticlePubMedGoogle Scholar
- Li Z, Reighard GL, Abbott AG, Bielenberg DG. Dormancy associated MADS genes from the EVG locus of peach [Prunus persica (L.) Batsch] have distinct seasonal and photoperiodic expression patterns. J Exp Bot. 2009;60:3521–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, Akagi T, et al. Functional and expressional analyses of PmDAM genes associated with endodormancy in Japanese apricot. Plant Physiol. 2011;157:485–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishitani C, Saito T, Ubi BE, Shimizu T, Itai A, Saito T, et al. Transcriptome analysis of Pyrus pyrifolia leaf buds during transition from endodormancy to ecodormancy. Sci Hortic. 2012;147:49–55.View ArticleGoogle Scholar
- Saito T, Bai S, Ito A, Sakamoto D, Saito T, Ubi BE, et al. Expression and genomic structure of the dormancy-associated MADS box genes MADS13 in Japanese pears (Pyrus pyrifolia Nakai) that differ in their chilling requirement for endodormancy release. Tree Physiol. 2013;33:654–67.View ArticlePubMedGoogle Scholar
- Yamane H, Ooka T, Jotatsu H, Hosaka Y, Sasaki R, Tao R. Expressional regulation of PpDAM5 and PpDAM6, peach (Prunus persica) dormancy-associated MADS-box genes, by low temperature and dormancy breaking reagent treatment. J Exp Bot. 2011;62:3481–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Leida C, Conejero A, Arbona V, Gomez-Cadenas A, Llacer G, Badenes ML, Rios G. Chilling-dependent release of seed and bud dormancy in peach associates to common changes in gene expression. PLoS One. 2012;7(5):e35777.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang S, Puryear J, Cairney J. A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993;11:113–6.View ArticleGoogle Scholar
- Copois V, Bibeau F, Bascoul-Mollevi C, Salvetat N, Chalbos P, et al. Impact of RNA degradation on gene expression profile: Assessment of different methods to reliably determine RNA quality. J Biotechnol. 2007;127:549–59.View ArticlePubMedGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from highthroughput sequencing reads. EMBnet J. 2011;1(17):10–2.View ArticleGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.View ArticlePubMedGoogle Scholar
- Zhang Z, Schwartz S, Wagner L, Miller W. A greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7:203–14.View ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.View ArticlePubMedPubMed CentralGoogle Scholar
- Li B, Dewey C. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticlePubMedGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35 suppl 1:D61–5.View ArticlePubMedGoogle Scholar
- Zdobnov EM, Apweiler R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17(9):847–8.View ArticlePubMedGoogle Scholar
- Ye J, George C, Irena Z, Ioana C, Steve R, Thomas LM. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134.View ArticlePubMedPubMed CentralGoogle Scholar