Transcriptomic analysis of ‘Suli’ pear (Pyrus pyrifolia white pear group) buds during the dormancy by RNA-Seq
© Liu et al.; licensee BioMed Central Ltd. 2012
Received: 20 June 2012
Accepted: 7 December 2012
Published: 12 December 2012
Bud dormancy is a critical developmental process that allows perennial plants to survive unfavorable environmental conditions. Pear is one of the most important deciduous fruit trees in the world, but the mechanisms regulating bud dormancy in this species are unknown. Because genomic information for pear is currently unavailable, transcriptome and digital gene expression data for this species would be valuable resources to better understand the molecular and biological mechanisms regulating its bud dormancy.
We performed de novo transcriptome assembly and digital gene expression (DGE) profiling analyses of ‘Suli’ pear (Pyrus pyrifolia white pear group) using the Illumina RNA-seq system. RNA-Seq generated approximately 100 M high-quality reads that were assembled into 69,393 unigenes (mean length = 853 bp), including 14,531 clusters and 34,194 singletons. A total of 51,448 (74.1%) unigenes were annotated using public protein databases with a cut-off E-value above 10-5. We mainly compared gene expression levels at four time-points during bud dormancy. Between Nov. 15 and Dec. 15, Dec. 15 and Jan. 15, and Jan. 15 and Feb. 15, 1,978, 1,024, and 3,468 genes were differentially expressed, respectively. Hierarchical clustering analysis arranged 190 significantly differentially-expressed genes into seven groups. Seven genes were randomly selected to confirm their expression levels using quantitative real-time PCR.
The new transcriptomes offer comprehensive sequence and DGE profiling data for a dynamic view of transcriptomic variation during bud dormancy in pear. These data provided a basis for future studies of metabolism during bud dormancy in non-model but economically-important perennial species.
Keywords‘Suli’ pear (Pyrus pyrifolia white pear group) Transcriptome Bud dormancy RNA-Seq
Dormancy is a complex phase of plant development that is necessary for survival under unfavorable environmental conditions. According to Lang, dormancy is a temporary suspension of visible growth of any plant structures containing meristems and can be divided into five well-defined phases: paradormancy, endodormancy, ecodormancy, and the two transitional phases between para- and endodormancy and endo- and ecodormancy. Dormancy transitions are regulated by short photoperiods and/or chilling temperatures. In pear low temperatures have been proven to control dormancy. Several excellent studies have investigated the physiological and molecular mechanisms of bud-dormancy transitions in perennial woody and herbaceous plants, including leafy spurge (Euphorbia esula Linn.)[3–5], poplar (Populus spp.)[6, 7], peach (Prunus persica [L.] Batsch)[8–12], apple (Malus × domestica Borkh)[13, 14], Japanese apricot (Prunus mume Siebold Zucc.)[15, 16], chestnut (Castanea sativa Mill.)[17–19], grape (Vitis vinifera linn.)[20, 21], raspberry (Rubus idaeus Linn.), kiwifruit (Actinidia spp.)[23, 24], and blackcurrant (Ribes nigrum Linn.). These results suggested that bud dormancy involves many biochemical pathways related to photoperiod, temperature, circadian clocks, water, energy, reactive oxygen species, and hormones. Several genes involved in dormancy transition have been identified, providing useful references for studying perennial plant dormancy. However, many unresolved questions remain about how many pathways are involved, how they interact, and significant differences in dormancy transition and regulation among species, genotypes, tissues, and environments.
Pears (Pyrus spp.) are among the world's most important perennial deciduous fruit trees and have a key feature of transition from growth to dormancy during their annual growth cycles. Studies of pear dormancy have focused mainly on the physiological level, including respiration, carbohydrate[27, 28] and protein metabolism, and chilling requirements. To date, few studies of pear dormancy at the molecular level have been conducted. Ubi et al. (2010) isolated two dormancy-associated MADS-box (DAM) genes and studied their expression patterns during the seasonal endodormancy transition phases in Japanese pear (Pyrus pyrifolia). Although these data highlighted the potential of molecular research to understand dormancy in this crop, they were insufficient to elucidate the molecular regulation mechanism. Furthermore, with global warming, many fruit trees, including pears, growing in warm areas have suffered from inadequate winter chill and showed advanced springtime and delayed autumnal phenologies, uneven foliation and blossoming, and long blooming periods, which are unfavorable for sustainable pear production[32–37]. Therefore, understanding the molecular mechanisms of pear dormancy transitions will greatly assist programs to breed cultivars with lower chilling requirements and to develop agronomic measures to cope with insufficient winter chill.
Traditionally, researchers have studied target nucleotide sequences by cloning, sequencing and comparing with known sequences, annotating their functions, then verifying their functions using tools such as RNAi, microarrays, and genetic transformation. These methods are very useful, but characterizing a large number of genes in a single experiment is difficult, especially with respect to specific genes. In recent years, RNA-sequencing (RNA-Seq) technology based on pyro- sequencing has become the most popular and powerful tool for species that lack reference genome information. RNA-Seq is less costly, more efficient, and allows faster gene discovery and more sensitive and accurate profiles of the transcriptome than microarray analysis or other techniques[38–45]. To better understand the molecular mechanisms of bud dormancy transition in pear trees, we used RNA-Seq technology to identify and characterize the expression of a large number of genes, especially those differentially expressed during dormancy progression.
The aim of the present study was to gain an understanding of molecular mechanisms during pear bud dormancy and establish a sound foundation for future molecular studies. We sequenced cDNA libraries from lateral flower buds of ‘Suli’ pear (Pyrus pyrifolia white pear group) from endo to ecodormancy stages using Illumina deep-sequencing technology. Approximately 100 M high-quality reads were obtained and assembled into 69,393 unigenes. Furthermore, four digital gene expression (DGE) libraries were constructed to compare gene expression patterns in different dormancy states using an upgraded digital gene expression system. The assembled and annotated transcriptome sequences and gene expression profiles will provide valuable resources for the identification of pear genes involved in bud dormancy.
Dormancy status of lateral flower buds in pears
Illumina sequencing and de novo assembly
Summary of the 'Suli' pear transcriptome during bud dormancy based on the RNA-Seq data
Total number of reads
Total base pairs (bp)
Average read length (bp)
Total number of contigs
Mean length of contigs (bp)
Total number of unigenes
Mean length of unigenes (bp)
Sequences with E-value < 10-5
Annotation of predicted proteins
There were 29,500 unigenes (approximately 42.5%) that mapped into 128 KEGG pathways (Additional file4). The maps with highest unigene representation were metabolic pathways (KO 01100; 7,043 unigenes, 25.1%), followed by biosynthesis of secondary metabolites (KO 01110; 3,349 unigenes, 11.4%), plant-pathogen interactions (KO 04626; 1,531 unigenes, 5.2%), and plant hormone signal transduction (KO 04075; 1,357 unigenes, 4.6%).
DGE library sequencing and mapping sequences to the reference transcriptome database
Summary of read numbers based on the RNA-Seq data from 'Suli' pear during bud dormancy
Total Clean Reads
Total Mapped Reads
Total Mapped Reads / Total Clean Reads
Perfect Match Reads
Perfect Match Reads / Total Clean Reads
≤2 bp Mismatch Reads
≤2 bp Mismatch Reads / Total Clean Reads
Unique Match Reads
Unique Match Reads / Total Clean Reads
Multi-Position Match Reads
Multi-Position Match Reads / Total Clean Reads
Total-Unmapped Reads / Total Clean Reads
Changes in gene expression profiles among dormancy stages
We analyzed the functions of the most differentially-expressed genes with the expression fold (log2Ratio ≥ 2) and false discovery rate (FDR ≤ 10-5) as the threshold in the Nov. 15-VS-Dec. 15, Dec. 15-VS-Jan. 15, and Jan. 15-VS-Feb. 15 comparisons (Additional file13) based on NCBI annotations. In the Nov. 15-VS-Dec. 15 comparison, only two of the 20 most up-regulated genes (stachyose synthase precursor and phytosulfokines precursor) had defined functions in the NCBI database, and only one (cytochrome P450) of the twenty most down-regulated genes had defined functions. Comparing the Dec. 15 and Jan. 15 libraries, three of the twenty most up-regulated genes (i.e., WRKY transcription factor 12, transcriptional factor TINY, blue copper protein precursor) and one (non-specific lipid-transfer protein 8-like) of the twenty most down-regulated genes had defined functions. Between the Jan. 15 and Feb. 15 libraries, three of the twenty most up-regulated genes had defined functions, including genes encoding LRR receptor-like serine, phd finger protein, and arabinogalactan-protein; none of the most down-regulated genes had defined functions.
Functional classification of DEGs during dormancy stages
In the Nov. 15-VS-Dec. 15, Dec. 15-VS-Jan. 15, and Jan. 15-VS-Feb. 15 comparisons, 666, 393, and 1,223 DEGs mapped to 104, 92, and 117 KEGG pathways, respectively. Of the 666 DEGs in the Nov. 15-VS-Dec. 15 comparison, 599 (89.9%) mapped to 15 pathways (Additional file14). Remarkably, specific enrichment of unigenes was observed for 15 pathways involved in metabolic processes, such as phenylpropanoid biosynthesis, ether lipid metabolism, ribosome, cutin, suberine and wax biosynthesis, stilbenoid, diarylheptanoid and gingerol biosynthesis, endocytosis, glycerophospholipid metabolism, oxidative phosphorylation, zeatin biosynthesis, flavonoid biosynthesis, starch and sucrose metabolism, fructose and mannose metabolism, other glycan degradation, and pentose and glucuronate interconversions. Comparing the Dec. 15 and Jan. 15 libraries, 343 (87.3%) DEGs were identified in 6 pathways: ribosome, ether lipid metabolism, glycerophospholipid metabolism, endocytosis, plant-pathogen interaction, and metabolic pathways (Additional file14). In the Jan. 15-VS-Feb. 15 comparison, 793 (64.8%) DEGs were significantly enriched in 15 pathways (Additional file14).
Clustering analysis of DEGs during dormancy stages
The second largest group (Group 4) contained 56 (29.5%) genes, from CL 12859 to CL 9922, that were down-regulated from Nov. 15 to Jan. 15 then up-regulated between Jan. 15 and Feb. 15. This group mainly included genes encoding proteins associated with photosynthesis metabolism, such as chlorophyll A/B binding protein, photosystem I reaction center subunit III, photosystem I reaction center subunit IV A, plastocyanin A, chloroplast oxygen-evolving enhancer protein 1, ribulose bisphosphate carboxylase, cytochrome P450, and magnesium-protoporphyrin IX monomethyl ester cyclase. Additionally, zinc finger protein, SPL domain class transcription factor, and basic helix-loop-helix domain-containing protein were encoded.
The third largest group (Group 5) contained 30 (15.8%) genes, from CL 10 to unigene 27694, that were up-regulated from Nov. 15 to Jan. 15, then down-regulated between Jan. 15 and Feb. 15. This group mainly included genes encoding proteins associated with oxidation-reduction reaction and resistance, such as blue copper protein, 2-oxoglutarate-dependent dioxygenase, cytochrome P450, heat shock protein, glycine rich protein, and dehydrin 1.
In Group 6, 17 genes (8.9%), from CL 13025 to unigene 26541, were up-regulated between Nov. 15 and Dec. 15, but were down-regulated between Dec. 15 and Jan. 15 and between Jan. 15 and Feb. 15. This group mainly involved in genes encoding dehydration-responsive element-binding protein, galactinol synthase 1, and GA2-oxidase.
In Group 7, 18 genes (9.5%), from CL 1316 to unigene 12417, were down-regulated between Nov. 15 and Dec. 15 and between Jan. 15 and Feb. 15 but were up-regulated between Dec. 15 and Jan. 15. These genes encoded proteins including fructose-bisphosphate aldolase cytoplasmic isozyme, tyrosine-protein phosphatase SIW14, glycerol-3-phosphate transporter, NAC domain class transcription factor, WRKY domain class transcription factor, and AP2 domain class transcription factor.
In Group 3, seven genes were down-regulated from Nov. 15 to Dec. 15, but were up-regulated between Dec. 15 and Feb. 15. Finally, two genes (Group 1) were up-regulated from Nov. 15 to Dec. 15, but were down-regulated between Jan. 15 and Feb. 15. Of these nine genes, only 2-aminoethanethiol dioxygenase isoform 2 and transferase were definitely annotated.
Gene expression analysis and Q-PCR validation
Phylogenetic analysis of dormancy-associated MADS-box (DAM) genes and their expression variations
In this study, using RNA-Seq technology, a ‘Suli’ pear cDNA library and four DGE libraries (from samples collected on Nov. 15, Dec. 15, Jan. 15, and Feb. 15, 2010–2011) were constructed and used to screen DEGs during dormancy. Surprisingly, we obtained 69,393 unique sequences, and of which 51,448 could be annotated, in total 48,725 genes, including 14,531 clusters and 34,194 singletons. Until 23 May, 2012, there were only 4,339 expressed sequence tags (ESTs) and 2,837 nucleotide sequences of Pyrus plants deposited in GenBank. Several recent studies have used traditional EST analyses to study dormancy in other species. Horvath et al. (2008) and Mazzitelli et al. (2007) identified nearly 1,000 genes and 327 ESTs associated with bud dormancy in leafy spurge and raspberry, respectively, using cDNA microarrays[5, 22]. Leida et al. (2010) identified nearly 400 ESTs associated with bud dormancy in peach by constructing four subtracted-cDNA libraries. These ESTs participated in different metabolic pathways related to photoperiod, temperature, circadian clocks, water, energy, reactive oxygen species, and hormones[5, 22, 46]. To our knowledge, this is the first report to use RNA-Seq technique to identify large numbers of genes involved in different metabolic pathways in pear bud dormancy. Compared with traditional EST analyses, RNA-Seq was less expensive, more efficient, and allowed faster gene discovery in bud dormancy studies.
Through RNA-seq analysis, we found that the numbers and expression profiles of DEGs differed at different times during dormancy. A total of 1,978, 1,024, and 3,468 genes were differentially expressed between Nov. 15 and Dec. 15, Dec. 15 and Jan. 15, and Jan. 15 and Feb. 15, respectively. These results showed that the number of DEGs was fewer in the endodormant stage (Nov. 15 and Dec. 15) than in the ecodormant stage (Jan. 15 and Feb. 15), increased with endodormancy-release to reach a maximum by Feb. 15. Hedley et al. (2010) reported that gene activity was lowest in the early stages of dormancy and peaked around the time of bud break in blackcurrant (Ribes nigrum L.). By analyzing KEGG pathways, we found DEGs that participated in several different pathways. Some pathways (such as starch and sucrose metabolism, circadian rhythm, and flavonoid biosynthesis) had been previously correlated to bud break in other species[5, 17, 22, 46], and some like phenylpropanoid biosynthesis, stilbenoid, diarylheptanoid and gingerol biosynthesis, zeatin biosynthesis, ether lipid metabolism, endocytosis, and glycerophospholipid metabolism were associated with bud break for the first time in this study. These data may suggest new research directions for understanding bud dormancy.
Some of the genes found in this work had been previously identified in other perennial plants. The DAM genes, widely described in perennial species such as leafy spurge, peach[9–12], raspberry, Japanese apricot, kiwifruit, and Japanese pear, are candidates for internal factors controlling endodormancy. In this study, we also found two DAM genes, and phylogenetic analysis revealed that CL 1161.contig2 was more closely related to PpMADS13-1 of Pyrus pyrifolia, whereas CL 1161.contig5 was more similar to PpMADS13-2 of Pyrus pyrifolia.Changes in the expression of CL 1161.contig2 and CL 1161.contig5 decreased with endodormancy release in lateral flower buds were consistent with the findings of earlier work comparing PpMADS13-1 and PpMADS13-2 gene expression in lateral leaf buds of Japanese pear, PpDAM5 and PpDAM6 in lateral vegetative buds and lateral flower buds of peach[10, 11], and all PmDAMs (PmDAM1 PmDAM6) in lateral vegetative buds of Japanese apricot. Our study, along with previous studies, suggested that DAM genes might play significant roles in the regulation of bud dormancy in ‘Suli’ pear.
The accumulation of dehydrin (DHN) is known to be associated with freezing tolerance in plants. Recent studies have reported that the accumulation of DHNs in woody plants correlates with seasonal transitions in dormancy and cold acclimation stages during winter[16, 48], but characterizations of DHN genes expressed during dormancy are limited. Yakovlev et al. (2008) found a reduction in the transcript levels of most of the 15 DHNs that they cloned as Norway spruce neared bud-burst. Garcia-Bañuelos et al. (2009) reported that transcripts of apple DHN were highly expressed in bark and bud tissues when the plant was dormant and cold hardy in midwinter, but were not expressed in early spring when cold hardiness was lost and buds were growing. Recently, several studies have identified DHN genes that were activated by ABA and C- repeat binding factor (CBF) in response to abiotic stresses[14, 50–52]. Intriguingly, in leafy spurge, ABA levels were elevated during endodormancy but dropped following the transition to ecodormancy. Horvath et al. (2008) found that CBF genes involved in cold regulated were up-regulated during the transition from para- to endo-dormancy. In the present study, one gene (CL 9148) encoding dehydrin showed significantly higher expression during the transition from endo- to eco-dormancy; thereafter, the expression level of this gene rapidly decreased, as indicated by DGE analysis and Q-PCR data. Based on previous studies, we speculated that DHN genes may act as signals and offer some protection for ‘Suli’ pear after the end of endodormancy, when pear often encounters unfavorable environmental conditions, such as cold. Therefore, more attention should be paid to ABA and CBF, which activate DHN genes, in future studies of transcriptional regulation related to the pear dormancy process.
Generally, sugar transport is thought to occur via H+/sugar symports that depend on a pH gradient generated by a plasma membrane H+-ATPase. Gevaudant et al. (2001) examined expression of the four H+-ATPase genes and reported that the levels of three H+-ATPase gene mRNAs increased, whereas the level of one H+-ATPase gene decreased in vegetative buds of peach trees after dormancy release. Mazzitelli et al. (2007) demonstrated that the plasma membrane H+-ATPase gene was highly expressed during the dormancy transition. Our results showed that the expression of a plasma membrane H+-ATPase gene (CL 1729) was up-regulated in pear buds during the endodormant maintenance period, down-regulated during endodormancy-release, and then up-regulated again. The expression patterns of plasma membrane H+-ATPase in ‘Suli’ pear were different from those of peach and raspberry buds, perhaps due to species-level or tissue-level differences.
In the present study, some genes encoding galactinol synthase (CL 9475), plastocyanin A (CL 1227), chlorophyll A/B binding protein (CL 7129, CL 514, CL 4948, CL 9178, CL 9961, CL 3279), and S-adenosylemethionine decarboxylase (CL 2122) were differentially expressed. Of these, chlorophyll A/B binding protein and S-adenosylemethionine decarboxylase were previously reported in other perennial plants[46, 55]. The expression levels of these genes changed significantly during the dormancy process. Thus, these genes may play roles in the regulation of bud dormancy in ‘Suli’ pear.
In addition, differentially-regulated transcription factors were identified in this study, including AP2 (unigene 1554), Zn-finger (unigene 16353), NAC (CL 7187), WRKY (unigene 19749), SPL (CL 3589), and bHLH (CL 12548). Of these, AP2, Zn-finger, and NAC were previously reported in leafy spurge and peach. Based on DGE analysis, the expression levels of the genes encoding these transcription factors significantly changed during dormancy in ‘Suli’ pear.
Although the molecular mechanisms associated with dormancy transitions in pear trees remain largely unknown, the present transcriptome analysis provided valuable information that could facilitate future studies on the detailed molecular functions of genes related to pear bud dormancy.
We obtained transcriptome data that provided comprehensive sequencing and DGE profiling data for a dynamic view of transcriptomic variation during the dormancy stage in pear. Physiological processes such as phenylpropanoid biosynthesis, stilbenoid, diarylheptanoid and gingerol biosynthesis, zeatin biosynthesis, ether lipid metabolism, endocytosis, glycerophospholipid metabolism, photosynthesis, phenylalanine metabolism, starch and sucrose metabolism were all differentially regulated during bud dormancy. Approximately 190 genes involved in many metabolic processes were significant differentially regulated during bud dormancy. Genes related to bud dormancy and their expression profiles at four time-points during dormancy were analyzed further. This offered new insights into the molecular mechanisms underlying pear bud dormancy. This provided a relatively complete molecular platform for future research on the progression of pear bud dormancy. To our knowledge, this work is the first to study pear bud dormancy using RNA-Seq.
‘Suli’ pear cultivars grafted onto Pyrus betulaefolia Bunge rootstocks were obtained from the Dangshan Germplasm Resources Center (Dangshan County, Anhui Province, China). Pear trees were 10 years old and were considered to be in the adult phase. Trees used in the experiment were not pruned or chemically treated during the experimental period. All samples were collected from the same trees at each stage. The current season’s growth shoots were collected on Nov. 15, Dec. 15, Jan. 15, Jan. 25and Feb. 15 from Nov. 2010 through Feb. 2011. Lateral flower bud samples were stored immediately in liquid N2 and then at −80°C until RNA extraction after picking.
Dormancy status of lateral flower buds
The dormancy status of lateral flower buds at four collection dates (Nov. 15, Dec. 15, Jan. 15, and Feb. 15) was estimated by evaluating excised shoots from field-grown trees. To measure the percentage of bud break, 12 current season’s growth shoots, with lengths of about 60 cm and bearing apical buds and 10–12 lateral flower buds, were collected and placed in water in 500 mL vials in a phytotron kept at day/night 25 ± 1/18 ± 1°C, with a 12-h photoperiod of white light (320 μ photon mol m-2 s-1) and 75% humidity. The water in the vials was changed and the basal ends of the shoots were cut every 2–3 d. After 21 d, the dormancy status was valued by percentage bud break; the beginning of bud break was defined as green leaf tips enclosing visible flowers. Lateral flower buds of shoots with bud break percentages < 50% were considered to remain in the endodormant stage.
RNA extraction, library preparation and RNA-seq
Thirty shoots (three biological replicates, with 10 shoots per replicate) with lengths of about 60 cm and bearing 10–12 lateral flower buds were sampled at each stage (Nov. 15, Dec. 15, Jan. 15, Jan. 25, and Feb. 15). The buds were sampled from three biological replications of each stage and produced an independent pool. Total RNA was extracted from lateral flower buds using the pine tree extraction protocol. The transcriptome assembly library was pooled by mixing equal quantities of RNA of five dormancy stages. The four DGE libraries consisted of separate RNA extracts from buds of four different dormancy stages, i.e., Nov. 15, Dec. 15, Jan. 15, and Feb. 15. Each library was pooled by mixing equal quantities of RNA from three biological replications for each stage. Each pool was sequenced once technically since the RNA-seq data are highly replicable with relatively little technical variation. The following protocols were performed by staff at the Beijing Genome Institute (BGI; Shenzhen, China). RNA integrity was confirmed with a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). 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 (≈200 bp) 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. Short double-stranded cDNA fragments were purified with QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands), resolved with EB buffer for end reparation and adding poly (A), then ligated to sequencing adapters. After purification via agarose gel electrophoresis, suitable fragments were enriched by PCR amplification before library sequencing using Illumina HiSeq™ 2000 (San Diego, CA, USA).
De novo assembly and function annotation
Raw sequence data in fastq format were filtered to remove reads containing adaptors, reads with more than 5% unknown nucleotides, and low-quality reads with more than 20% bases of quality value ≤ 10. Only clean reads were used in the following analysis. The sequences from the Illumina sequencing were deposited in the NCBI Sequence Read Archive (Accession SRX147917). First, transcriptome de novo assembly was carried out by BGI using the short-read assembly program Trinity with the following parameters: min_contig_length = 100, min_glue = 2, group_pairs_distance = 250, path_reinforcement_distance = 75, bfly_opts = '-V 5 --edge-thr = 0.05 --stderr', min_kmer_cov = 2. Meanwhile, all reads of approximately 100 M (i.e. transcriptome sequencing reads and RNA-Seq reads) were mapped back to the apple genome reference (ftp://ftp.jgi-psf.org/pub/JGI_data/phytozome/v8.0/Mdomestica/assembly/Mdomestica_196.fa.gz) to identify continuous gene regions using SOAPsplice software (Release 1.9; http://soap.genomics.org.cn/soapsplice.html). Secondly, we realigned all the transcripts onto the reference genome. When more than one transcript were placed in one gene region and they each other had an overlap less than 24 bp, we connected them into a longer transcript. A total of 31,727 transcripts assembled by Trinity were connected into 19,309 transcripts. Then the redundancy of unigenes was removed by TGICL (v.2.1) with options ‘-l 40 -v 25’. Finally based on gene family clustering, the unigenes were divided into two classes: clusters and singletons. The former was prefixed with 'CL' and the latter with 'unigene'. The id number of each unigene followed this prefix. In a cluster, the similarity between unigenes was more than 70%.
Blastx alignment (E-value < 10-5) between unigenes and protein databases such as nr, Swiss-Prot, KEGG, and GO was performed, and the best-aligning results were used to determine the sequence direction of unigenes. When different databases conflicted, the results were prioritized in the order: nr, Swiss-Prot, KEGG, and GO. When a unigene did not align to any of the databases, ESTScan was used to decide its sequence direction. GO annotation was analyzed by Blast2GO software (v.2.5.0). KEGG pathway annotation was performed using Blastall software against the KEGG database. The assembled sequences could be searched using the Gene-ID listed in Additional materials (Additional file3).
Analysis and mapping of DGE reads
Raw sequence data in fastq format were filtered to remove reads containing adaptors, reads with more than 10% unknown nucleotides, and low-quality reads with more than 50% bases of quality value ≤ 5. The sequences from the DGE analysis were deposited in the NCBI Sequence Read Archive under accession numbers SRX148326 (Nov. 15), SRX148327 (Dec. 15), SRX148328 (Jan. 15), and SRX148329 (Feb. 15). Clean reads were mapped to our transcriptome reference database using the short oligonucleotide analysis package SOAPaligner/soap2, allowing mismatches of no more than two bases. The unique mapped reads were used in subsequent analyses. For gene expression analysis, the number of unique-match reads was calculated and then normalized to RPKM (reads per kb per million reads). The RPKM method eliminates the influence of different gene lengths and sequencing discrepancies on the calculation of gene expression, so that the calculated gene expressions can be directly compared among samples.
Evaluation of DGE libraries
We used FDR ≤ 0.001, the absolute value of log2Ratio ≥ 1, and the RPKM value of each gene for either sample over 10 as the thresholds to judge the significance of gene expression differences, where log2Ratio indicates the degree of differential expression between two samples and was the ratio of RPKM values for the treatment and control samples. This analysis found genes with significantly differential expression among samples prior to GO function and KEGG pathway analyses.
Clustering of gene expression profiles
Hierarchical cluster analysis of 190 gene expression patterns was performed with cluster and Java Treeview software. The log2Ratio for each gene was used for the hierarchical clustering analysis.
A phylogenetic tree was constructed based on the nucleotide sequences of two unigenes and 18 additional DAM genes. The tree was generated with MEGA (v. 4.0.1) software, using the neighbor-joining method.
Real-time quantitative RT-PCR analysis
Total RNA used for Q-PCR analysis was extracted from lateral flower buds of four different dormancy stages, i.e., Nov. 15, Dec. 15, Jan. 15, and Feb. 15, using three biological replicates of about 300 buds. Total RNA was extracted as described above, genomic DNA was removed with DNase I, and total RNA concentration was measured. First-strand cDNA was synthesized from 4 μg of DNA-free RNA using the Revert Aid™ First Strand cDNA Synthesis Kit (Fermentas, Glen Burnie, MD, USA) according to the manufacturer’s instructions. The cDNA was used as the template for Q-PCR. Primer sequences (designed using primer 3,http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) are listed in Additional file16. The Q-PCR mixture (20 μl total volume) contained 10 μl of SYBR Premix ExTaq™ (Takara, Kyoto, Japan), 0.4 μl of each primer (10 μM), 2 μl of cDNA, and 7.2 μl of RNase-free water. The reactions were performed on a LightCycler 1.5 instrument (Roche, Basel, Switzerland) according to the manufacturer’s instructions. The two-step Q-PCR program began with 30 seconds at 95°C, followed by 40 cycles of 95°C for 5 seconds and 60°C for 20 seconds. Template-less controls for each primer pair were included in each run. The specificity of Q-PCR primers was confirmed by melting curve and sequencing of Q-PCR products. Expression was calculated as 2-∆ ∆ Ct and normalized to that of the actin gene (PpActin, JN684184), and data were managed with the Data Processing System (DPS, v. 7.05, Zhejiang University, Hangzhou, China).
Beijing Genome Institute
Basic local alignment search tool
Digital gene expression
Data processing system
Expressed sequence tag
False discovery rate
The Kyoto encyclopedia of genes and genomes database
KEGG Orthology ids
Molecular evolutionary genetics analysis
NCBI non-redundant database
Real-time quantitative PCR
Reads per kb per million reads
Short oligonucleotide analysis package.
Acknowledgements and funding
This research was financed by the earmarked fund for Modern Agro-industry Technology Research System (nycytx-29).
We acknowledge the Dangshan Germplasm Resources Center for providing materials.
- Lang G, Early J, Martin G, Darnell R: Endo-, para-, and ecodormancy: physiological terminology and classification for dormancy research. Hortscience. 1987, 22: 371-377.Google Scholar
- Heide O, Prestrud A: Low temperature, but not photoperiod, controls growth cessation and dormancy induction and release in apple and pear. Tree Physiol. 2005, 25 (1): 109-114.View ArticlePubMedGoogle Scholar
- Chao WS, Serpe MD, Anderson JV, Gesch RW, Horvath DP: Sugars, hormones, and environment affect the dormancy status in underground adventitious buds of leafy spurge (Euphorbia esula). Weed Science. 2006, 54 (1): 59-68.View ArticleGoogle Scholar
- Horvath DP, Chao WS, Anderson JV: Molecular analysis of signals controlling dormancy and growth in underground adventitious buds of leafy spurge. Plant Physiol. 2002, 128 (4): 1439-1446.PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
- Bohlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S: CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science. 2006, 312 (5776): 1040-View ArticlePubMedGoogle Scholar
- Ruttink T, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, Boerjan W, Rohde A: A molecular timetable for apical bud formation and dormancy induction in poplar. Plant Cell. 2007, 19 (8): 2370-2390.PubMed CentralView ArticlePubMedGoogle Scholar
- Yooyongwech S, Sugaya S, Sekozawa Y, Gemma H: Differential adaptation of high- and low-chill dormant peaches in winter through aquaporin gene expression and soluble sugar content. Plant Cell Rep. 2009, 28 (11): 1709-1715.View ArticlePubMedGoogle Scholar
- Leida C, Conesa A, Llácer G, Badenes ML, Ríos G: Histone modifications and expression of DAM6 gene in peach are modulated during bud dormancy release in a cultivar-dependent manner. New Phytol. 2012, 193 (1): 67-80.View ArticlePubMedGoogle Scholar
- Yamane H, Ooka T, Jotatsu H, Sasaki R, Tao R: Expression analysis of PpDAM5 and PpDAM6 during flower bud development in peach (Prunus persica). Sci Hortic. 2011, 129 (4): 844-848.View ArticleGoogle Scholar
- Yamane H, Ooka T, Jotatsu H, Hosaka Y, Sasaki 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 (10): 3481-3488.PubMed CentralView ArticlePubMedGoogle Scholar
- Bielenberg DG, Wang Y, Li ZG, Zhebentyayeva T, Fan SH, Reighard GL, Scorza R, Abbott AG: Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. Tree Genetics & Genomes. 2008, 4 (3): 495-507.View ArticleGoogle Scholar
- Brunel N, Leduc N, Poupard P, Simoneau P, Mauget JC, Viemont JD: KNAP2, a class I KN1-like gene is a negative marker of bud growth potential in apple trees (Malus domestica [L.] Borkh.). J Exp Bot. 2002, 53 (378): 2143-2149.View ArticlePubMedGoogle Scholar
- Wisniewski M, Norelli J, Bassett C, Artlip T, Macarisin D: Ectopic expression of a novel peach (Prunus persica) CBF transcription factor in apple (Malus × domestica) results in short-day induced dormancy and increased cold hardiness. Planta. 2011, 233 (5): 971-983.View ArticlePubMedGoogle Scholar
- 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-497.PubMed CentralView ArticlePubMedGoogle Scholar
- Yamane H, Kashiwa Y, Kakehi E, Yonemori K, Mori H, Hayashi K, Iwamoto K, Tao R, Kataoka I: Differential expression of dehydrin in flower buds of two Japanese apricot cultivars requiring different chilling requirements for bud break. Tree Physiol. 2006, 26 (12): 1559-1563.View ArticlePubMedGoogle Scholar
- Ramos A, Pérez-Solís E, Ibáñez C, Casado R, Collada C, Gómez L, Aragoncillo C, Allona I: Winter disruption of the circadian clock in chestnut. Proc Natl Acad Sci U S A. 2005, 102 (12): 7037-PubMed CentralView ArticlePubMedGoogle Scholar
- Santamaria ME, Rodriguez R, Canal MJ, Toorop PE: Transcriptome analysis of chestnut (Castanea sativa) tree buds suggests a putative role for epigenetic control of bud dormancy. Ann Bot. 2011, 108: 485-498.PubMed CentralView ArticlePubMedGoogle Scholar
- Ibañez C, Ramos A, Acebo P, Contreras A, Casado R, Allona I, Aragoncillo C: Overall alteration of circadian clock gene expression in the chestnut cold response. PLoS One. 2008, 3 (10): e3567-PubMed CentralView ArticlePubMedGoogle Scholar
- Pacey-Miller T, Scott K, Ablett E, Tingey S, Ching A, Henry R: Genes associated with the end of dormancy in grapes. Funct Integr Genomics. 2003, 3 (4): 144-152.View ArticlePubMedGoogle Scholar
- Ophir R, Pang X, Halaly T, Venkateswari J, Lavee S, Galbraith D, Or E: 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 (4): 403-423.View ArticlePubMedGoogle Scholar
- Mazzitelli L, Hancock RD, Haupt S, Walker PG, Pont SDA, McNicol J, Cardle L, Morris J, Viola R, Brennan R: Co-ordinated gene expression during phases of dormancy release in raspberry (Rubus idaeus L.) buds. J Exp Bot. 2007, 58 (5): 1035-1045.View ArticlePubMedGoogle Scholar
- Richardson AC, Walton EF, Meekings JS, Boldingh HL: Carbohydrate changes in kiwifruit buds during the onset and release from dormancy. Sci Hortic. 2010, 124 (4): 463-468.View ArticleGoogle Scholar
- Wu R-M, Walton EF, Richardson AC, Wood M, Hellens RP, Varkonyi-Gasic E: Conservation and divergence of four kiwifruit SVP-like MADS-box genes suggest distinct roles in kiwifruit bud dormancy and flowering. J Exp Bot. 2012, 63 (2): 797-807.PubMed CentralView ArticlePubMedGoogle Scholar
- Hedley PE, Russell JR, Jorgensen L, Gordon S, Morris JA, Hackett CA, Cardle L, Brennan R: Candidate genes associated with bud dormancy release in blackcurrant (Ribes nigrum L.). BMC Plant Biol. 2010, 10: 202-PubMed CentralView ArticlePubMedGoogle Scholar
- Bi L, Zhang Y, Poudyal B: Effects of growth regulators on the respiration metabolism of pear buds during dormant period. Frontiers of Agriculture in China. 2011, 5 (1): 45-50.View ArticleGoogle Scholar
- Zimmerma RH, Faust M: Pear Bud Metabolism - Seasonal Changes in Glucose Utilization. Plant Physiol. 1969, 44 (9): 1273-1276.View ArticleGoogle Scholar
- Marafon AC, Citadin I, Amarante L, Herter FG, Hawerroth FJ: Chilling privation during dormancy period and carbohydrate mobilization in Japanese pear trees. Scientia Agricola. 2011, 68 (4): 462-468.View ArticleGoogle Scholar
- Tamura F, Tanabe K, Itai A, Tanaka H: Protein changes in the flower buds of Japanese pear during breaking of dormancy by chilling or high-temperature treatment. J Am Soc Hortic Sci. 1998, 123 (4): 532-536.Google Scholar
- Kretzschmar A, Brighenti L, Rufato L, Pelizza T, Silveira F, Miquelutti D, Faoro I: Chilling Requirement for Dormancy Bud Break in European Pear. Acta Hortic. 2011, 909: 85-87.View ArticleGoogle Scholar
- Ubi BE, Sakamoto D, Ban Y, Shimada T, Ito A, Nakajima I, Takemura Y, Tamura F, Saito T, Moriguchi T: Molecular cloning of dormancy-associated MADS-box gene homologs and their characterization during seasonal endodormancy transitional phases of Japanese pear. J Am Soc Hortic Sci. 2010, 135 (2): 174-182.Google Scholar
- Cannell M, Smith R: Climatic warming, spring budburst and forest damage on trees. J Appl Ecol. 1986, 177-191.Google Scholar
- Young E: Timing of high temperature influences chilling negation in dormant apple trees. J Am Soc Hortic Sci. 1992, 117 (2): 271-273.Google Scholar
- Luedeling E, Zhang M, Girvetz EH: Climatic changes lead to declining winter chill for fruit and nut trees in California during 1950–2099. PLoS One. 2009, 4 (7): e6166-PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson J, Horvath D, Chao W, Foley M: Bud Dormancy in Perennial Plants: A Mechanism for Survival. Dormancy and Resistance in Harsh Environments. Edited by: Lubzens E, Cerda J, Clark M. 2010, Heidelberg: Springer Berlin, 69-90. vol. 21View ArticleGoogle Scholar
- Campoy J, Ruiz D, Egea J: Dormancy in temperate fruit trees in a global warming context: A review. Sci Hortic. 2011, 130 (2): 357-372.View ArticleGoogle Scholar
- Luedeling E, Girvetz EH, Semenov MA, Brown PH: Climate change affects winter chill for temperate fruit and nut trees. PLoS One. 2011, 6 (5): e20155-PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10 (1): 57-63.PubMed CentralView ArticlePubMedGoogle Scholar
- Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, Allan AC, Ferguson IB, Chen KS: Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012, 13: 19-PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez-Gomez P, Crisosto CH, Bonghi C, Rubio M: New approaches to Prunus transcriptome analysis. Genetica. 2011, 139 (6): 755-769.View ArticlePubMedGoogle Scholar
- Tang Q, Ma XJ, Mo CM, Wilson IW, Song C, Zhao H, Yang YF, Fu W, Qiu DY: An efficient approach to finding Siraitia grosvenorii triterpene biosynthetic genes by RNA-seq and digital gene expression analysis. BMC Genomics. 2011, 12: 343-PubMed CentralView ArticlePubMedGoogle 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-PubMed CentralView ArticlePubMedGoogle Scholar
- Schuster SC: Next-generation sequencing transforms today’s biology. Nat Methods. 2008, 5: 16-18.View ArticlePubMedGoogle Scholar
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26 (10): 1135-1145.View ArticlePubMedGoogle Scholar
- Leida C, Terol J, Marti G, Agusti M, Llacer G, Badenes ML, Rios G: Identification of genes associated with bud dormancy release in Prunus persica by suppression subtractive hybridization. Tree Physiol. 2010, 30 (5): 655-666.View ArticlePubMedGoogle Scholar
- Kosová K, Vítámvás P, Prášil IT: The role of dehydrins in plant response to cold. Biol Plantarum. 2007, 51 (4): 601-617.View ArticleGoogle Scholar
- Yakovlev I, Asante D, Fossdal C, Partanen J, Junttila O, Johnsen Ø: Dehydrins expression related to timing of bud burst in Norway spruce. Planta. 2008, 228 (3): 459-472.View ArticlePubMedGoogle Scholar
- Garcia-Bañuelos M, Gardea A, Winzerling J, Vazquez-Moreno L: Characterization of a Midwinter-Expressed Dehydrin (DHN) Gene from Apple Trees (Malus domestica). Plant Mol Biol Rep. 2009, 27 (4): 476-487.View ArticleGoogle Scholar
- Pulla RK, Kim YJ, Kim MK, Senthil KS, In JG, Yang DC: Isolation of a novel dehydrin gene from Codonopsis lanceolata and analysis of its response to abiotic stresses. BMB Rep. 2008, 41 (4): 338-343.View ArticlePubMedGoogle Scholar
- Mingeot D, Dauchot N, Van Cutsem P, Watillon B: Characterisation of two cold induced dehydrin genes from Cichorium intybus L. Mol Biol Rep. 2009, 36 (7): 1995-2001.View ArticlePubMedGoogle Scholar
- Weiss J, Egea-Cortines M: Transcriptomic analysis of cold response in tomato fruits identifies dehydrin as a marker of cold stress. J Appl Genet. 2009, 50 (4): 311-319.View ArticlePubMedGoogle Scholar
- Améglio T, Guilliot A, Lacointe A, Julien JL, Alves G, Valentin V, Pétel G: Water relations in winter: effect on bud break of walnut tree. Dormancy in plants: from whole plant behaviour to cellular control. Edited by: Viémont JD, Crabbé J. 2000, Wallingford, UK: CAB International, 109-120.View ArticleGoogle Scholar
- Gevaudant F, Petel G, Guilliot A: Differential expression of four members of the H+-ATPase gene family during dormancy of vegetative buds of peach trees. Planta. 2001, 212 (4): 619-626.View ArticlePubMedGoogle Scholar
- Horvath DP, Anderson JV, Soto-Suarez M, Chao WS: Transcriptome analysis of leafy spurge (Euphorbia esula) crown buds during shifts in well-defined phases of dormancy. Weed Science. 2006, 54 (5): 821-827.View ArticleGoogle Scholar
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993, 11 (2): 113-116.View ArticleGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517.PubMed CentralView ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29 (7): 644-652.PubMed CentralView ArticlePubMedGoogle Scholar
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999, 138-148.Google Scholar
- Li R, Yu C, Li Y, Lam T-W, Yiu S-M, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967.View ArticlePubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628.View ArticlePubMedGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.PubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57 (1): 289-300.Google Scholar
- de Hoon MJL, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20 (9): 1453-1454.View ArticlePubMedGoogle Scholar
- Saldanha AJ: Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599.View ArticlePubMedGoogle Scholar
- Zhang D, Yu B, Bai J, Qian M, Shu Q, Su J, Teng Y: Effects of high temperatures on UV-B/visible irradiation induced postharvest anthocyanin accumulation in ‘Yunhongli No. 1’ (Pyrus pyrifolia Nakai) pears. Sci Hortic. 2012, 134: 53-59.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.