Transcriptomic analysis of the underground renewal buds during dormancy transition and release in 'Hangbaishao' peony (Paeonia lactiflora)

Abstract

Paeonia lactiflora is one of the most famous species of herbaceous peonies with gorgeous flowers. Bud dormancy is a crucial developmental process that allows P. lactiflora to survive unfavorable environmental conditions. However, little information is available on the molecular mechanism of the bud dormancy in P. lactiflora. We performed de novo transcriptome sequencing using the Illumina RNA sequencing platform for the underground renewal buds of P. lactiflora 'Hangbaishao' to study the molecular mechanism underlying its bud dormancy transition (the period from endodormancy to ecodormancy) and release (the period from ecodormancy to bud elongation and sprouting). Approximately 300 million high-quality clean reads were generated and assembled into 207,827 (mean length = 828 bp) and 51,481 (mean length = 1250 bp) unigenes using two assembly methods named "Trinity" and "Trinity+PRICE", respectively. Based on the data obtained by the latter method, 32,316 unigenes were annotated by BLAST against various databases. Approximately 1,251 putative transcription factors were obtained, of which the largest number of unique transcripts belonged to the basic helix-loop-helix protein (bHLH) transcription factor family, and five of the top ten highly expressed transcripts were annotated as dehydrin (DHN). A total of 17,705 simple sequence repeat (SSR) motifs distributed in 13,797 sequences were obtained. The budbreak morphology, levels of indole-3-acetic acid (IAA) and abscisic acid (ABA), and activities of guaiacol peroxidase (POD) and catalase (CAT) were observed. The expression of 20 interested unigenes, which annotated as DHN, heat shock protein (HSP), histone, late elongated hypocotyl (LHY), and phytochrome (PHY), and so on, were also analyzed. These studies were based on morphological, physiological, biochemical, and molecular levels and provide comprehensive insight into the mechanism of dormancy transition and release in P. lactiflora. Transcriptome dataset can be highly valuable for future investigation on gene expression networks in P. lactiflora as well as research on dormancy in other non-model perennial horticultural crops of commercial significance.

Fig 1. Ornamental organs and experimental samples of P. lactiflora ‘Hangbaishao’.

(A) The pink purple flower in the early florescence stage. (B) The pink white flower in the end of florescence stage. (C) The purple red young leaves. (D) The purple red ripe fruits. (E) The underground renewal buds sampled for deep sequencing and transcriptome research.

Fig 2. Change in the daily average temperatures (℃) in sampling base from November 23, 2012 to March 11, 2013.

Fig 3. The length distribution (bp) of ‘Hangbaishao’ unigenes.

(A) The length distribution of ‘Hangbaishao’ unigenes assembled by “Trinity+PRICE”. (B) Comparison of the length distribution of unigenes assembled by “Trinity” and “Trinity+PRICE”.

Fig 4. Characteristics of homology search of ‘Hangbaishao’ unigenes against the NR database.

(A) Venn diagram of number of unigenes annotated by BLASTx with an E-value threshold of 10^-5 against protein databases. The numbers in the circles indicate the number of unigenes annotated by single or multiple databases. (B) E-value distribution of the top BLASTx hits for each unigene. (C) Similarity distribution of the top BLASTx hits for each unigene. (D) Species distribution of the top BLASTx hits for all homologous sequences.

Fig 5. COG function classification of ‘Hangbaishao’ unigenes.

Fig 6. GO classification of ‘Hangbaishao’ unigenes.

Fig 7. The proportions of sequences related to carbohydrate metabolism, plant hormones and antioxidant enzymes.

(A) The proportion of sequences related to carbohydrate metabolism. (B) The proportion of sequences related to five internal plant hormones. (C) The proportion of sequences related to four antioxidant enzymes.

Fig 8. The number of DEGs derived from 15 pairwise comparisons of six treatments.

Fig 9. Expression of 20 interested unigenes related to the bud dormancy transition or release.