Comparative transcriptome analysis identifies differentially expressed genes between normal and late-blooming Siberian apricot

Abstract

Blooming date is an important trait in fruit tree species. Although several quantitative trait loci confirming blooming date were identified in Prunus spp., the molecular mechanism underlying it remains unclear. Arising from this, the transcriptomes of normal blooming and late-blooming Siberian apricot (P. sibirica L.) flower buds were analyzed using RNA-seq technology. A total of 68,855 unigenes were de novo assembled, among which 1204 were differentially expressed between normal and late blooming. Gene ontology enrichment analysis revealed that biological processes were enriched with metabolic processes. The catalytic-related gene transcripts between the two types of blooming were significantly changed in the molecular function. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that 156 genes were successfully annotated and 75 pathways enriched. Genes for gibberellin biosynthesis were up-regulated in normal blooming, whereas abscisic acid degradation-related genes were also up-regulated in normal blooming. Moreover, circadian rhythms related genes including EARLY FLOWERING 4, LATE ELONGATED HYPOCOTYL and CIRCANDIAN CLOCK ASSOCIATED1 were all up-regulated in normal blooming, indicating that circadian rhythms have a very important role in controlling blooming date. Furthermore, zinc finger protein CONSTANS-LIKE 12 was blasted onto the quantitative trait loci region on linkage group 4 in peach. However, changes in the abundance of key flowering genes such as SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1, FLOWERING LOCU T, LEAFY and FLOWERING LOCUS C were not significantly different, indicating that further investigation should explore the function of these genes on blooming date. The outcomes of this study will provide a valuable platform for further research on the molecular mechanism of blooming date in Prunus.

Keywords

Blooming date / Differentially expressed genes (DEGs) / Flowering time / Siberian apricot / Transcriptome

References

[1] Alburquerque N, García-Montiel F, Carrillo A, Burgos L. Chilling and heat requirements of sweet cherry cultivars and the relationship between altitude and the probability of satisfying the chill requirements. Environ Exp Bot, 2008, 64(2): 162-170.

[2] Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol, 1990, 215(3): 403-410.

[3] Amasino R. Seasonal and developmental timing of flowering. Plant J Cell Mol Biol, 2010 61 6 1001

[4] Anderson JL, Seeley SD. Bloom delay in deciduous fruits. Hortic Rev, 1993, 15: 97-144.

[5] Andrés F, Coupland G. The genetic basis of flowering responses to seasonal cues. Nat Rev Genet, 2012, 13(9): 627-639.

[6] Ashburner M, Ball CA, Blake JA. Gene ontology: tool for the unification of biology. Nat Genet, 2000, 25: 25-29.

[7] Bai S, Saito T, Sakamoto D, Ito A, Fujii H, Moriguchi T. Transcriptome analysis of japanese pear (Pyrus pyrifolia Nakai) flower buds transitioning through endodormancy. Plant Cell Physiol, 2013, 54(7): 1132-1151.

[8] Ballester J, Company RSI, Arus P, Vicente MCD. Genetic mapping of a major gene delaying blooming time in almond. Plant Breed, 2001, 120(3): 268-270.

[9] Barros PM, Gonçalves N, Saibo NJM, Oliveira MM. Cold acclimation and floral development in almond bud break: insights into the regulatory pathways. J Exp Bot, 2012, 63(12): 4585-4596.

[10] Blázquez MA, Weigel D. Integration of floral inductive signals in Arabidopsis. Nature, 2000, 404(6780): 889-892.

[11] Böhlenius H, Huang T, Charbonnel-Campaa L, Brunner AM, Jansson S, Strauss SH, Nilsson O. CO/FT regulatory module controls timing of flowering and seasonal growth cessation in trees. Science, 2006, 312(5776): 1040-1043.

[12] Boss PK, Bastow RM, Mylne JS, Dean C. Multiple pathways in the decision to flower: enabling, promoting, and resetting. Plant Cell, 2004, 16(Suppl): S18-S31.

[13] Bowman JL, Alvarez J, Weigel D, Meyerowitz EM, Smyth DR. Control of flower development in Arabidopsis thaliana by APETALA1 and interacting genes. Development, 1991, 119(3): 721-743.

[14] Campoy JA, Ruiz D, Egea J, Rees DJG, Celton JM, Martínez-Gómez P. Inheritance of flowering time in apricot (Prunus armeniaca L.) and analysis of linked quantitative trait loci (QTLs) using simple sequence repeat (SSR) markers. Plant Mol Biol Report, 2011, 29(2): 404-410.

[15] Campoy JA, Ruiz D, Allderman L, Cook N, Egea J. The fulfilment of chilling requirements and the adaptation of apricot (Prunus armeniaca l.) in warm winter climates: an approach in murcia (Spain) and the western cape (South Africa). Eur J Agron, 2012, 37(1): 43-55.

[16] Castède S, Campoy JA, García JQ, Le DL, Lafargue M, Barreneche T. Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements. New Phytol, 2014, 202(2): 703-715.

[17] Coen ES, Meyerowitz EM. The war of the whorls: genetic interactions controlling flower development. Nature, 1991, 353(6339): 31-37.

[18] Conesa A, Götz S, Garcíagómez JM, Terol J, Talón M, Robles M. Blast2go: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics, 2005, 21(18): 3674-3676.

[19] Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol, 2005, 139(1): 5-17.

[20] Dirlewanger E, Moing A, Rothan C, Svanella L, Pronier V, Guye A. Mapping QTLs controlling fruit quality in peach (Prunus persica (L.) Batsch). Theor Appl Genet, 1999, 98(1): 18-31.

[21] Dodd AN, Salathia N, Hall A, Kévei E, Tóth R, Nagy F. Plant circadian clocks increase photosynthesis, growth, survival, and competitive advantage. Science, 2005, 309(5734): 630-633.

[22] Dogramaci M, Foley ME, Chao WS, Christoffers MJ, Anderson JV. Induction of endodormancy in crown buds of leafy spurge (Euphorbia esula L.) implicates a role for ethylene and cross-talk between photoperiod and temperature. Plant Mol Biol, 2013, 81(6): 577-593.

[23] Doyle MR, Davis SJ, Bastow RM, Mcwatters HG, Kozmabognár L, Nagy F. The ELF4 gene controls circadian rhythms and flowering time in Arabidopsis thaliana. Nature, 2002, 419(6902): 74-77.

[24] Erez A. Bud dormancy; phenomenon, problems and solutions in the tropics and subtropics. Temperate fruit crops in warm climates, 2000, Dordrecht: Springer 17 48

[25] Eriksson ME, Millar AJ. The circadian clock. A plant’s best friend in a spinning world. Plant Physiol, 2003, 132(2): 732-738.

[26] Fan S, Bielenberg DG, Zhebentyayeva TN, Reighard GL, Okie WR, Holland D. Mapping quantitative trait loci associated with chilling requirement, heat requirement and bloom date in peach (Prunus persica). New Phytol, 2010, 185(4): 917-930.

[27] Ferrandiz C, Gu Q, Martienssen R, Yanofsky MF. Redundant regulation of meristem identity and plant architecture by fruitfull, APETALA1 and cauliflower. Development, 2000, 127(4): 725-734.

[28] Gong DY, Ho CH. The Siberian High and climate change over middle to high latitude Asia. Theoret Appl Climatol, 2002, 72(1): 1-9.

[29] Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol, 2011, 29(7): 644-652.

[30] Guedon Y, Legave JM. Analyzing the time-course variation of apple and pear tree dates of flowering stages in the global warming context. Ecol Model, 2008, 219(1–2): 189-199.

[31] Gumus M, Kasifoglu S. Performance and emission evaluation of a compression ignition engine using a biodiesel (apricot seed kernel oil methyl ester) and its blends with diesel fuel. Biomass Bioenerg, 2010, 34(1): 134-139.

[32] Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc, 2013, 8(8): 1494-1512.

[33] Hardtke CS, Okamoto H, Stoop-Myer C, Deng XW. Biochemical evidence for ubiquitin ligase activity of the Arabidopsis COP1 interacting protein 8 (CIP8). Plant J Cell Mol Biol, 2002, 30(4): 385-394.

[34] Hur J, Ahn J. The change of first-flowering date over South Korea projected from downscaled IPCC AR5 simulation: peach and pear. Int J Climatol, 2015, 35(8): 1926-1937.

[35] Jofuku KD, Den Boer BG, Van Montagu M, Okamuro JK. Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell, 1994, 6(9): 1211-1225.

[36] Jyothi T, Suttle JC, Chao WS, Horvath DP, Anderson JV. Transcriptome analysis identifies novel responses and potential regulatory genes involved in seasonal dormancy transitions of leafy spurge (Euphorbia Esula L.). BMC Genom, 2008 9 1 536

[37] Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M. KEGG for linking genomes to life and the environment. Nucl Acids Res, 2008 36 Database issue D480

[38] Karlberg A, Englund M, Petterle A, Molnar G, Sjödin A, Bako L. Analysis of global changes in gene expression during activity-dormancy cycle in hybrid aspen apex. Plant Biotechnol, 2010, 27(1): 1-16.

[39] Klie M, Debener T. Identification of superior reference genes for data normalisation of expression studies via quantitative PCR in hybrid roses (Rosa hybrida). BMC Res Notes, 2011 4 1 518

[40] Kobayashi Y, Kaya H, Goto K, Iwabuchi M, Araki T. A pair of related genes with antagonistic roles in mediating flowering signals. Science, 1999, 286(5446): 1960-1962.

[41] Lang GA, Early JD, Martin GC, Darnell RL. Endo-, para-, and ecodormancy: physiological terminology and classification for dormancy research. Hortscience, 1987, 22: 271.

[42] Lavee S. Dormancy and bud break in warm climates; considerations of growth regulator involvement. Acta Hortic, 1973, 22: 371-377.

[43] Li XM, Qin ZQ, Tu JF, Yang FC, Zhu HY, Liu XQ. Study on chilling requirement of different pear cultivars. J Henan Agric Sci, 2011, 40(7): 126-129.

[44] Liu G, Zhang L, Wang T, Zhao X. Observation of flower bud differentiation and phenology in Siberian apricot (Prunus sibirica L.). J Shenyang Agric Univ, 1996, 27(1): 95-96. (in Chinese)

[45] Luedeling E, Brown PH. A global analysis of the comparability of winter chill models for fruit and nut trees. Int J Biometeorol, 2011, 55(3): 411-421.

[46] Mandel MA, Gustafsonbrown C, Savidge B, Yanofsky MF. Molecular characterization of the Arabidopsis floral homeotic gene APETALA1. Nature, 1992, 360(6401): 273-277.

[47] McClung CR. Circadian rhythms in plants. Annu Rev Plant Biol, 2001, 52(1): 139-162.

[48] McClung CR. Plant circadian rhythms. Plant Cell, 2006, 18(4): 792-803.

[49] Meehl GA, Stocker TF, Collins WD, Friedlingstein P, Gaye T, Gregory JM, Kitoh A, Knutti R, Murphy JM, Noda A, Raper SCB, Watterson IG, Weaver AJ, Zhao ZC Solomon S, Qin D, Manning M Global climate projections. Climate change: the physical science basis contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change, 2007, Cambridge: Cambridge University Press 1 99

[50] Moon J, Lee H, Kim M, Lee I. Analysis of flowering pathway integrators in Arabidopsis. Plant Cell Physiol, 2005, 46(2): 292-299.

[51] Mornya P, Cheng F. The levels of hormone and carbohydrate in autumn and non-autumn flowering tree peonies. Can J Plant Sci, 2011, 91(6): 991-998.

[52] Naor A, Flaishman M, Stern R, Moshe A, Erez A. Temperature effects on dormancy completion of vegetative buds in apple. J Am Soc Hortic Sci, 2003, 128(5): 636-641.

[53] Ohme-Takagi M, Shinshi H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell, 1995, 7(2): 173-182.

[54] Okie WR, Blackburn B. Increasing chilling reduces heat requirement for floral budbreak in peach. Hortscience, 2011, 46(2): 245-252. A Publication of the American Society for Horticultural Science)

[55] Olsen JE. Light and temperature sensing and signaling in induction of bud dormancy in woody plants. Plant Mol Biol, 2010, 73(1–2): 37-47.

[56] Onouchi H, Igeño MI, Périlleux C, Graves K, Coupland G. Mutagenesis of plants overexpressing constans demonstrates novel interactions among Arabidopsis flowering-time genes. Plant Cell, 2000, 12(6): 885-900.

[57] Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet, 2011, 12(2): 87-98.

[58] Porri A, Torti S, Romera-Branchat M, Coupland G. Spatially distinct regulatory roles for gibberellins in the promotion of flowering of Arabidopsis under long photoperiods. Development, 2012, 139(12): 2198-2209.

[59] Quero-Garcia J, Ivancic A, Lebot V. Taro and cocoyam. Root and tuber crops, 2010, New York: Springer 149 172

[60] Quilot B, Wu BH, Kervella J, Génard M, Foulongne M, Moreau K. QTL analysis of quality traits in an advanced backcross between Prunus persica cultivars and the wild relative species P. davidiana. Theor Appl Genet, 2004, 109(4): 884-897.

[61] Reinoso H, Luna V, Pharis RP, Bottini R. Dormancy in peach (Prunus persica) flower buds. V. Anatomy of bud development in relation to phenological stage. Can J Bot, 2002, 80(6): 656-663.

[62] Richardson EA. A model for estimating the completion of rest for ‘redhaven’ and ‘elberta’ peach trees. HortScience, 1974, 9: 331-332.

[63] Rieu I, Ruizrivero O, Fernandezgarcia N, Griffiths J, Powers SJ, Gong F. The gibberellin biosynthetic genes ATGA20ox1 and ATGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J, 2008, 53(3): 488-504.

[64] Rinne PL, Welling A, Vahala J, Ripel L, Ruonala R, Kangasjärvi J. Chilling of dormant buds hyperinduces FLOWERING LOCUS T and recruits GA-inducible 1,3-β-glucanases to reopen signal conduits and release dormancy in Populus. Plant Cell, 2011, 23(1): 130-146.

[65] Robinson MD, Mccarthy DJ, Smyth GK. Edger: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 2010, 26(1): 139-140.

[66] Rohde A, Ruttink T, Hostyn V, Sterck L, Driessche KV, Boerjan W. Gene expression during the induction, maintenance, and release of dormancy in apical buds of poplar. J Exp Bot, 2007, 58(15–16): 4047-4060.

[67] Ruiz D, Campoy JA, Egea J. Chilling and heat requirements of apricot cultivars for flowering. Environ Exp Bot, 2007, 61(3): 254-263.

[68] 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.

[69] Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-Sommer Z, Yanofsky MF, Coupland G. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science, 2000, 288(5471): 1613-1616.

[70] Sánchez-Pérez R, Dicenta F, Martínez-Gómez P. Inheritance of chilling and heat requirements for flowering in almond and QTL analysis. Tree Genet Genomes, 2012, 8(2): 379-389.

[71] Suárez-López P, Wheatley K, Robson F, Onouchi H, Valverde F, Coupland G. Constans mediates between the circadian clock and the control of flowering in Arabidopsis. Nature, 2001, 410(6832): 1116-1120.

[72] Suttle JC. Involvement of ethylene in potato microtuber dormancy. Plant Physiol, 1998, 118(3): 843-848.

[73] Tabuenca MC. Chilling requirements in almond. Spanish. Anal Estacion Exp Aula Dei, 1972, 11: 325-329.

[74] Valverde F, Mouradov A, Soppe W, Ravenscroft D, Samach A, Coupland G. Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science, 2004, 303(5660): 1003-1006.

[75] Van Verk MC, Hickman R, Pieterse CMJ, Van Wees SCM. RNA-seq: revelation of the messengers. Trends Plant Sci, 2013, 18(4): 175-179.

[76] Verde I, Quarta R, Cedrola C, Dettori MT. QTL analysis of agronomic traits in a BC1 peach population. Acta Hort, 2002, 592: 291-297.

[77] Wang L. Evaluation of Siberian Apricot (Prunus sibirica, L.) germplasm variability for biodiesel properties. J Am Oil Chem Soc, 2012, 89(9): 1743-1747.

[78] Wang LR, Zhu GR, Fang WC, Zuo QY. Estimating models of the chilling requirement for peach. Acta Hortic Sin, 2003, 30(4): 379-383.

[79] Wang Z, Gerstein M, Snyder M. RNA-seq: a revolutionary tool for transcriptomics. Nat Rev Genet, 2009, 10(1): 57-63.

[80] Wilkins O, Nahal H, Foong J, Provart NJ, Campbell MM. Expansion and diversification of the Populus R2R3-MYB family of transcription factors. Plant Physiol, 2009, 149(2): 981-993.

[81] Winter D, Vinegar B, Nahal H, Ammar R, Wilson GV, Provart NJ. An “electronic fluorescent pictograph” browser for exploring and analyzing large-scale biological data sets. PLoS ONE, 2007 2 8 e718

[82] Wu J, Seng S, Sui J, Vonapartis E, Luo X, Gong B. Gladiolus hybridus ABSCISIC ACID INSENSITIVE 5 (GhABI5) is an important transcription factor in ABA signaling that can enhance Gladiolus corm dormancy and Arabidopsis seed dormancy. Front Plant Sci, 2015, 6(e0166): 749-753.

[83] Yamaguchi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol, 2008, 59: 225-251.

[84] Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z. Wego: a web tool for plotting go annotations. Nucl Acids Res, 2006, 34(Web Server issue): W293-W297.

[85] Yoo SK, Chung KS, Kim J, Lee JH, Hong SM, Yoo SJ. CONSTANS activates SUPPERSSOR OF OVEREXPRESSION OF CONSTAN 1 through FLOWERING LOCUS T to promote flowering in Arabidopsis. Plant Physiol, 2005, 139(2): 770-778.

[86] Zhang G, Ming C, Chen X, Xu Z, Shan G, Li LC. Phylogeny, gene structures, and expression patterns of the ERF gene family in soybean (Glycine max L.). J Exp Bot, 2008, 59(15): 4095-4107.

[87] Zhebentyayeva T, Fan S, Olukolu B, Barakat A, Leida C, Badenes ML, Bielenberg D, Reighard G, Okie W, Abbott AG (2010) From genetics to epigenetics in control of chilling requirements and flowering time in peach. In: 5th international Rosaceae genomics conference, November 2010, Stellenbosch (South Africa)

[88] Zhu G, Ye N, Zhang J. Glucose-induced delay of seed germination in rice is mediated by the suppression of ABA catabolism rather than an enhancement of ABA biosynthesis. Plant Cell Physiol, 2009, 50(3): 644-651.