Comparative transcriptome analysis of nonchilled, chilled, and late-pink bud reveals flowering pathway genes involved in chilling-mediated flowering in blueberry

Comparative transcriptome analysis of nonchilled, chilled, and late-pink bud reveals flowering... Background: Blueberry cultivars require a fixed quantity of chilling hours during winter endo-dormancy for vernalization. In this study, transcriptome analysis using RNA sequencing data from nonchilled, chilled, and late pink buds of southern highbush blueberry ‘Legacy’ was performed to reveal genes associated with chilling accumulation and bud break. Results: Fully chilled ‘Legacy’ plants flowered normally whereas nonchilled plants could not flower. Compared to nonchilled flower buds, chilled flower buds showed differential expression of 89% of flowering pathway genes, 86% of MADS-box genes, and 84% of cold-regulated genes. Blueberry orthologues of FLOWERING LOCUS T (FT) did not show a differential expression in chilled flower buds (compared to nonchilled flower bud) but were up-regulated in late-pink buds (compared to chilled flower bud). Orthologoues of major MADS-box genes were significantly up- regulated in chilled flower buds and down-regulated in late-pink buds. Functional orthologues of FLOWERING LOCUS C (FLC) were not found in blueberry. Orthologues of Protein FD (FD), TERMINAL FLOWER 1 (TFL1),and LEAFY (LFY) were down-regulated in chilled flower buds and in late-pink buds compared to nonchilled flower bud. Conclusions: The changes from nonchilled to chilled and chilled to late-pink buds are associated with transcriptional changes in a large number of differentially expressed (DE) phytohormone-related genes and DE flowering pathway genes. The profile of DE genes suggests that orthologues of FT, FD, TFL1, LFY, and MADS-box genes are the major genes involved in chilling-mediated blueberry bud-break. The results contribute to the comprehensive investigation of the vernalization-mediated flowering mechanism in woody plants. Keywords: Chilling requirement, Cold hardness, Flowering time control, Freezing tolerance, Vaccinium corymbosum, Vernalization, Woody plant Background Climate change in the last 40 years has caused earlier Winter dormancy (endo-dormancy) is essential for de- shifts in the onset of the growing season for trees (e.g., ciduous fruit crop survival [1, 2]. Under inductive low 2.3 days/decade in temperate Europe) and increased temperatures in the fall, deciduous woody fruit and nut temperature fluctuation [4]. Early onset of the growing crops are acclimated to develop freezing tolerance; season causes insufficient chilling hours and prevents meanwhile, accumulation of effective chilling hours is bud-break in fruit trees. Increased temperature fluctu- stimulated [3]. Sufficient chilling accumulation gives ation during plant bloom turns early season frosts into a plants full vernalization, which is a prerequisite for danger, with freezing injuries to flowers and young fruits bud-break in the spring. [5]. Plant breeding to manipulate chilling requirements and develop improved freeze tolerant cultivars are con- sidered to be long-term solutions to mitigate reduced * Correspondence: songg@msu.edu Plant Biotechnology Resource and Outreach Center, Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA © The Author(s). 2018 Open Access This 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. Song and Chen BMC Plant Biology (2018) 18:98 Page 2 of 13 winter chill, decrease freezing damages, and secure de- natural winter conditions flowered normally (Fig. 1). Tran- ciduous fruit production [6]. scriptome comparison of the chilled to nonchilled flower Seasonal flowering plays a significant role in a plant’slife buds using Trinity and a non-annotated transcriptome cycle and is controlled by a network of flowering pathway reference (Reftrinity) of tetraploid blueberry in GenBank genes [7–9]. FLOWERING LOCUS C (FLC) is a key (Accession number: SRX2728597) [19, 20] revealed 37,000 regulator in the vernalization pathway of winter-annual differentially expressed genes (DEGs) and 47,000 DE Arabidopsis thaliana ecotypes [7]. In winter wheat and transcripts. barley, VERNALIZATION2 (VRN2) is a major regulator of To conduct transcriptome analysis, Reftrinity (180,000 vernalization-mediated flowering [8]. FLC and VRN2 ana- genes and 250,000 isoforms) was annotated using Trino- logs are recorded in peach (Prunus persica). These analogs tate [20]. The annotation resulted in 14,000 known are the DORMANCY ASSOCIATED MADS-box (DAM) genes from 30,000 genes and 55,000 isoforms of Reftri- genes which are a cluster of six MADS-box transcription nity. With this annotated reference, 64% of known blue- factors. The loss of all or part of the DAMs resulted in the berry genes and isoforms showed differential expression non-vernalized peach evergrowing mutant [10, 11]. The in the comparison of chilled to nonchilled flower buds DAM genes are considered alternatives to FLC in regulat- (herein referring to DEGs/DE transcripts in chilled ing vernalization-mediated chilling requirement and flow- flower buds). Chilling affected expression of numerous ering [10, 12]. However, the DAM genes show high genes simultaneously in blueberry flower buds. similarities to A. thaliana AGAMOUS-LIKE 24 (AGL24) and SHORT VEGETATIVE PHASE (SVP)genes [12, 13]. Effect of chilling accumulation on flowering pathway Additionally, functional analysis of DAMs to reveal their genes roles in chilling-mediated flowering through reverse gen- Differential expression was detected in chilled blueberry etics has not been reported in peach. To date, neither a flower buds for 89% of flowering pathway orthologues of functional FLC-LIKE nor a VRN2-LIKE gene has been known A. thaliana flowering pathway genes (Table 1; verified in Vaccinium plants [14]. Fig. 1). One of the top two orthologues of FT (VcFT: Blueberries and cranberries are the most important hereafter Vc before any A. thaliana gene refers to the Vaccinium fruits due to their high antioxidant and blueberry orthologue gene) showed a slight decrease anti-inflammatory capacities [15]. Deciphering the while the other showed a slight increase. Similar results mechanism of vernalization/chilling-mediated flowering were observed for CONSTANS-LIKE2 orthologues will facilitate molecular breeding of blueberry cultivars (Table 1). for low-chilling requirement. To investigate flowering Orthologues of seven major MADS-box genes in flower- responses under nonconducive conditions, functional ing pathway genes of blueberry buds, APETALA1 (AP1), analysis of a blueberry FLOWERING LOCUS T (VcFT) FRUITFULL (FUL), SUPPRESSOR OF OVEREXPRESSION gene has been conducted in highbush blueberry (Vacci- OF CONSTANS 1 (SOC1), CAULIFLOWER (CAL), FLC, nium corymbosum L.) [16–18]. Overexpression of VcFT AGL24,and SVP, were significantly up-regulated (Table 1). (about 2900-fold increase in leaf tissues) caused continu- The orthologues of Protein FD (FD), TERMINAL FLOWER ous and precocious flowering in in vitro shoots and in 1(TFL1),and LEAFY (LFY), and ACTIN-RELATED PRO- one-year old ‘Aurora’ plants [16]. However, Overexpres- TEIN6 (ARP6) were down-regulated genes. Both VcFLC sion of VcFT was not capable of fulfilling all chilling re- and VcSVP showed a contrasting response of FLC to quirements in blueberry [17, 18]. Over 80% of the flower vernalization in A. thaliana. The down-regulated VcFD and buds in two- and three-year old VcFT-overexpressing up-regulated VcSVP support the DE VcFT results [7]. How- plants could not bloom under greenhouse conditions ever, decreased expression of VcLFY contradicts the de- without a chilling period. To discover vernalization/chil- creased VcTFL1 and increased VcAP1, VcAGL24,and ling-responsive genes, transcriptome analyses were con- VcSOC1 expressions. The unchanged expression of VcFT ducted with blueberry flower buds. The profile of and decreased expression of VcLFY orthologues could pre- differentially expressed (DE) genes will facilitate our un- vent the flowering of chilled flower buds during plant derstanding of the role of vernalization/chilling in blue- vernalization. berry bud-break. The MADS-box transcription factor and flowering re- pressor, FLC, and FRIGIDA (FRI) are the major genes Results regulating the vernalization of A. thaliana. FRI activates Identification of DE transcripts in chilled flower buds FLC, but vernalization represses FLC [7]. Both VcFRI and Nonchilled flower buds from the southern highbush blue- VcFLC in blueberry are present in the annotated Reftri- berry ‘Legacy’ were grown in a heated greenhouse through nity. VcFRI showed a high similarity to FRI (e = − 98) while the winter season. These buds did not flower in the fol- two VcFLC transcripts showed a lower similarity to FLC lowing spring. In contrast, flower buds fully chilled under (e = − 31). These two VcFLC transcripts were annotated to Song and Chen BMC Plant Biology (2018) 18:98 Page 3 of 13 Fig. 1 Effect of vernalization of blueberry flowering. a Nonchilled flower buds could not bloom. b Early bloom of fully chilled flower buds. c Percentages of differentially expressed (DE) genes and transcripts in chilled flower buds (in comparison to nonchilled flower buds). *DE genes including A. thaliana genes and genes annotated to other species SEP1 of A. thaliana and MADS6 of rice, respectively late-pink buds to those in chilled flower buds (hereafter (Table 1). In chilled flower buds, VcFRI did not show dif- referring to DEGs/DE transcripts in late-pink buds) re- ferential expression whereas expression of VcFLC was in- sulted in 28,000 DE isoforms, which were annotated to creased about 25-fold for the top-two VcFLC candidate 11,000 known genes. genes (Table 1). When all 26 potential VcFLC candidate genes (e < − 20) were included in the analysis, four candi- Expression of flowering pathway genes in late-pink buds dates showed a down-regulation with a fold-change of 0.6– The major DE flower pathway genes in chilled flower buds 0.8. The other 22 candidate genes showed an up-regulation were compared to those in late-pink buds (Additional file with an average of 4.9-fold (Additional file 1:Table S2).The 2: Table S1). The major flowering pathway genes VcFT, inconsistency between VcFLC and FLC response to chill- VcSOC1, VcAP1,VcFUL, VcFLC, VcSVP,and VcLFY showed ing/vernalization suggests that a different or more compli- decreased expression in late-pink buds. In A. thaliana, FT, cated mechanism is involved in vernalization-mediated SOC1, AP1, FUL,and LFY promote flowering while SVP flowering in tetraploid, blueberry plants. Additionally, all and FLC are flowering repressors [7, 21]. The expression other DEGs involved in known vernalization-mediated changes of VcFD and VcTFL1 in late-pink buds were flowering pathway were up-regulated in chilled flower buds similar to those in chilled flower buds (Table 1). The de- with the exception of VcARP6 (c49456_g2, Log FC = − 7.7) creased VcFD expression was related, at least in part, to (Additional file 2:Table S1). VcFT down-regulation. Decreased VcTFL1 expression in late-pink buds was associated with a down-regulation of Identification of DE transcripts in late-pink buds VcAP1 and VcLFY. Additionally, decreased VcSVP and Fully chilled blueberry flower buds remained dormant VcFLC expression were associated with a decrease in VcFT under chilling conditions in January until continuous expression. warm conditions in April drove dormancy release and bloom. To investigate the effect of dormancy release on Expression of MADS-box genes gene expression, RNA sequencing data was obtained Blueberry DE MADS-box and DAMs genes were identi- from late-pink buds. Comparative transcripts analysis in fied using A. thaliana MADS-box genes and Japanese Song and Chen BMC Plant Biology (2018) 18:98 Page 4 of 13 Table 1 Differentially expressed major flowering pathway genes in both chilled flower buds (CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: Log (LPB/CB) 2 2 Arabidopsis_ID|Vc_ LogFC for LogFC for FDR for FDR for e-value_to_blast Gene_name Annotation_by_ Query_ transcript_ID chilled buds late-pink buds chilled buds late-pink buds Trinotate PmDAMs AT5G60910.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 9E-44 FUL, AGL8 AGL11_ARATH AT1G69120.1|c92021_g1_i1 −0.6 −1.8 4E-02 5E-15 1E-80 AP1, AGL7 AGL8_SOLTU AT5G60910.1|c88116_g1_i1 −0.4 0.8 4E-02 1E-10 5E-91 FUL, AGL8 AGL8_SOLTU AT1G69120.1|c88116_g8_i1 3.2 −5.4 1E-16 2E-59 4E-94 AP1, AGL7 AGL8_SOLTU AT1G69120.1|c81830_g1_i1 3.5 −2.1 2E-42 4E-24 2E-52 AP1, AGL7 AGL9_PETHY AT1G69120.1|c88116_g7_i1 3.6 −2.7 3E-33 7E-41 3E-90 AP1, AGL7 AP1_SINAL AT5G15840.1|c77980_g1_i1 −0.5 0.7 2E-03 3E-10 1E-66 CO, FG COL2_ARATH AT5G15840.1|c76265_g1_i1 2.3 −2.5 2E-09 6E-14 4E-30 CO, FG COL2_ARATH AT5G60910.1|c79125_g1_i1 −0.3 −1.1 4E-02 2E-13 2E-35 FUL, AGL8 DEFA_ANTMA AT5G60910.1|c67980_g1_i1 3.4 −4.6 7E-46 2E-70 5E-36 FUL, AGL8 DEFA_ANTMA AT5G61850.1|c96427_g2_i1 −2.4 −5.7 4E-13 4E-24 7E-102 LFY, LFY3 FLO_ANTMA AT5G61850.1|c96427_g2_i2 −2.1 −6.7 2E-13 9E-31 3E-98 LFY, LFY3 FLO_ANTMA AT1G65480.1|c84088_g2_i3 −0.9 −4.5 2E-02 7E-12 5E-95 FT HD3A_ORYSJ AT1G65480.1|c84088_g2_i5 1.1 −4.9 1E-03 9E-28 1E-94 FT HD3A_ORYSJ AT1G26310.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 2E-56 CAL, AGL10 MADS6_ORYSJ AT2G45660.1|c89673_g3_i1 2.6 −4.0 1E-09 3E-15 3E-53 SOC1, AGL20 SOC1_ARATH AT2G45660.1|c86010_g2_i1 3.9 −4.4 1E-08 2E-09 2E-54 SOC1, AGL20 SOC1_ARATH AT2G22540.1|c91377_g1_i7 1.0 −1.6 3E-02 1E-03 5E-69 SVP, AGL22 SVP_ARATH AT2G22540.1|c90829_g2_i1 1.2 −1.2 1E-03 5E-04 4E-77 SVP, AGL22 SVP_ARATH AT4G24540.1|c91377_g1_i7 1.0 −1.6 3E-02 1E-03 1E-39 AGL24 SVP_ARATH AT5G10140.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 2E-31 FLC, FLF, MADS6_ORYSJ AGL25 AT5G10140.1|c88116_g6_i1 4.1 −1.9 5E-91 1E-26 7E-31 FLC, FLF, SEP1_ARATH AGL25 AT4G18960.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 1E-72 MAF2, AGL31 AGL11_ARATH PmDAM1 AT2G14210.1|c80388_g1_i2 3.9 3.2 4E-02 5E-09 3E-79 MAF2, AGL31 MAD23_ORYSJ PmDAM1 AT5G60910.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 7E-54 MAF2, AGL31 MADS6_ORYSJ PmDAM1 AT4G22950.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 1E-42 MAF3, AGL70 AGL11_ARATH PmDAM1 AT4G18960.1|c88293_g3_i1 4.7 −4.7 5E-21 5E-20 5E-59 MAF4 AG_TOBAC AT5G15800.2|c88116_g6_i1 4.1 −1.9 5E-91 1E-26 2E-87 MAF4 SEP1_ARATH PmDAM1 AT1G24260.2|c68983_g1_i1 2.5 0.0 3E-06 0E + 00 3E-42 MAF5, AGL68 . AT5G15800.2|c81830_g1_i1 3.5 −2.1 2E-42 4E-24 3E-89 MAF5, AGL68 AGL9_PETHY PmDAM1 #N/A: No differential expression apricot (Prunus mume) DAMs (PmDAMs) (Additional VcSOC1 and VcFLC could have multiple functions in file 1: Table S2; Table 1). Orthologues of 62 A. thaliana blueberry. MADS-box gene were identified in blueberry (unpub- Three DE PmDAM orthologues (VcPmDAM1, lished data). DE orthologues of 53 and 44 MADS-box VcPmDAM2, and VcPmDAM5) were identified and showed genes were detected in chilled flower buds and late-pink high similarities to four A. thaliana MADS-box genes, buds, respectively. These orthologues include the major MADS AFFECTING FLOWERING 2 (MAF2), MAF4, flowering pathway genes VcFLC, VcSOC1, VcSVP, VcAP1, MAF5 and FOREVER YOUNG FLOWER.In A. thaliana, VcFUL, VcCAL,and VcAGL24 (Table 1). The annotated MAF2, MAF4,and MAF5 are FLC paralogs. MAF2, MAF5 VcSOC1 (c86010_g2_i1) showed high similarities to 25 and FLC are down-regulated and MAF5 is up-regulated A. thaliana MADS-box genes. Similarly, the annotated during vernalization [22]. In contrast, the blueberry ortho- VcFLC homologues were similar to 23 A. thaliana logues VcFLC, VcMAF1, VcMAF2, VcMAF4, and VcMAF5 MADS-box genes (Table 1). The results suggest that were up-regulated while VcMAF3 was repressed in chilled Song and Chen BMC Plant Biology (2018) 18:98 Page 5 of 13 flower buds (Table 1; Additional file 1: Table S2). Add- phyotohormone orthologues than the chilled buds itionally, three DE VcPmDAMs were annotated to the (Fig. 2). TheDE phytohormonegenes suggest thepo- homologues VcSOC1, VcSVP, VcAP1,and VcSEP1.The tential involvement of these phytohormones during up-regulated VcPmDAM1 homologues were the only DE chilling and flowering. orthologue in chilled flower buds. In late-pink buds, 75% Due to the tetraploid nature of ‘Legacy’, orthologues of of DE VcPmDAM1 homologues and all DE VcPmDAM5 each A. thaliana gene used for query often have more homologues were down-regulated while DE VcPmDAM2 than one homologue (Fig. 2). Thus, an average of the Fold Change was up-regulated (Additional file 1:Table S2).These blue- changes (Log ) for all the DE transcripts that berry MADS-box genes showed significant changes in re- are orthologues derived from a single A. thaliana query sponse to both chilling and flower bud breaking (Table 1). gene was used to represent the overall change of each However, during vernalization, the responses of VcFLC, phytohormone-related gene (Fig. 3). Increased expres- VcMAFs, and VcPmDAMs diverges from FLC’s response sion of ABA1, ABA2 and NINE-CIS-EPOXYCAROTE- to vernalization in A. thaliana. NOID DIOXYGENASE 3 (NCED3) in the ABA biosynthesis pathway were seen in chilled flower buds. Response of phytohormone-related genes in chilled and ABA1 and ABA2 continued to increase and NCED3 de- late-pink buds creased in late-pink buds (Fig. 3). The increased expres- For both chilled buds and late-pink buds, DE transcripts sion of these orthologues indicates a potential increase showed high similarities to the pathway genes for five in ABA biosynthesis during vernalization. Regardless of major phytohormones (Additional file 3: Table S3). Over decreased ABA1 expression, increased NCED3 expres- 50% of the DE blueberry orthologues were related to sion suggests that there is an increase in ABA biosyn- abscisic acid (ABA), ethylene, auxin, and gibberellin thesis during floral bud break (Fig. 3). (GA) genes while 25% were related to cytokinin genes DE orthologues of ethylene signaling pathway genes (Fig. 2). The late-pink bud showed more DE were all up-regulated in chilled buds and down-regulated in flowering buds (Fig. 3). These orthologues are consid- ered regulators for freezing tolerance in A. thaliana. Indole-3-acetic acid (IAA) and GA biosynthesis pathway orthologues were up-regulated in chilled flower buds and decreased in late-pink buds (Fig. 3). GIBBERELLIN 3-BETA-DIOXYGENASE 2 (GA3OX2) in the GA pathway and AUX1 in the IAA pathway were major DEGs with high expression changes. The DE ARABIDOPSIS RE- SPONSE REGULATORS (ARRs) orthologues included two A-type, two B-type and five ARR-like genes in chilled buds and flowering buds. One B-type orthologue (ARR10) was suppressed only in chilled buds (Fig. 3). Gene networks of DEGs in chilled flower buds and late-pink buds Over-represented Gene Ontology (GO) terms (P <0.05) were grouped to visualize gene networks of the annotated DE transcripts using the GOslim_Plant as the selected GO file and A. thaliana annotation as the reference. The DE transcripts were classified in 70 and 73 over-represented GO terms for chilled flower buds and late-pink buds, respectively (Fig. 4). The over-represented GO terms for chilled flower buds and late-pink buds were identical except for two GO terms (Fig. 4), suggesting that the same tran- Fig. 2 Response of phytohormone-related genes and transcripts in scripts responded to temperature changes in these buds. chilled flower buds (vs. nonchilled flower buds) and late-pink buds (vs. Thedifferencein “biological_process” was two additional chilled flower buds). a Percentage of differentially expressed (DE) orthologues of A. thaliana genes (The number of DE A. thaliana genes over-represented GO terms (GO:0007610-behavior and ÷ total number of A. thaliana genes × 100). b Percentage of DE GO:040029-regulation of gene expression, epigenetic) transcripts of transcripts of blueberry phytohormone-related genes found in late-pink buds but not in chilled flower (The number of DE transcripts ÷ total number of transcripts × 100) buds. (Fig. 4). Song and Chen BMC Plant Biology (2018) 18:98 Page 6 of 13 Fig. 3 Average fold changes (Log FC) of differentially expressed homologues for each of the phytohormone-related orthologue of A. thaliana in chilled flower buds (CB) [vs. nonchilled flower buds (NB)] and late-pink buds (LPB) [vs. chilled flower buds (CB)]. LogFC for chilled buds: Log (CB/ NB). LogFC for late-pink buds: Log (LPB/CB). a Abscisic acid biosynthesis pathway genes [46]. b Ethylene biosynthesis and signaling pathway genes [45, 47]. c Gibberellin biosynthesis pathway genes [48]. d Two-component ARABIDOPSIS RESPONSE REGULATORS (ARR) [49, 50]. e Auxin biosynthesis pathway genes [51]. The bars represent standard deviation Over-represented GO terms in chilled flower buds terms related to growth, response to stress, and revealed the impact of chilling on GO terms in three reproduction (Fig. 4). Thegenenetwork basedon categories (Fig. 4). The over-represented GO terms over-represented GO terms facilitate our understand- in “biological_process” revealed that vernalization/ ing of the role DEGs in both chilled flower buds and chilling affected expression of genes in multiple GO late-pink buds (Fig. 4). Song and Chen BMC Plant Biology (2018) 18:98 Page 7 of 13 cd Fig. 4 (See legend on next page.) Song and Chen BMC Plant Biology (2018) 18:98 Page 8 of 13 (See figure on previous page.) Fig. 4 Gene networks of differentially expressed genes (DEGs) in chilled flower buds and late-pink buds. DEGs in chilled flower buds were identified in comparison to nonchilled flower buds while DEGs in late-pink buds were identified in comparison to chilled flower buds. The ontology file of GOSlim_Plants in BiNGO was used to identify over-represented GO terms (P < 0.05). a Comparison of the gene network in chilled flower buds to late-pink buds; white nodes and black edges are present in both gene networks; red nodes and edges are present only in the chilled buds; and green nodes and edges are present only in the late-pink buds. The number in each circle is a GO identity number. A gene network in chilled flower buds (b “Biological_process” c “Cellular component” d “Molecular function”). I, II, and III in b show GO terms related to stress, plant growth, and reproduction, respectively Validation of the expression of selected genes For blueberries, expressed sequence tags have been gener- In chilled flower buds and late-pink buds, qRT-PCR were ated from blueberry flower buds [25]. However, compara- used to validate DE transcripts of VcFD, VcTFL1,and tive transcriptome analyses of different stage floral buds VcARP6 (Fig. 5). The results suggested high-reliability of have not been documented. the RNA-seq data. The roles of VcFT, VcLFY and VcARP6 in vernalization- Discussion mediated blueberry flowering Transcriptome analysis is an effective approach to study Overexpression of VcFT (expression level > 2000-fold flowering pathway genes [23, 24]. Using this approach, in leaf tissues) resulted in precocious flowering [16, DEGs in response to vernalization have been identified in 17]. However, the high expression level did not com- Japanese pear (Pyrus pyrifolia Nakai) and oriental lily [24]. pletely reverse the need for chilling for normal plant Fig. 5 Comparison of RNA sequencing and qRT-PCR analysis of three differentially expressed genes in (a) chilled flower buds (compared to nonchilled flower buds) (b) and late-pink buds (compared to chilled flower buds). Eukaryotic translation initiation factor 3 subunit H is the internal control. Log fold-change was calculated by -ΔΔCt = − [(Ct – Ct ) – (Ct – Ct ) ]. Average Log fold-change ± standard deviation of three 2 GOI nom tissue 1 GOI nom tissue 2 2 biological replicates. Significant average fold-change determined using a Student’s t-test is denoted. An asterisk (*) indicates p <0.001 Song and Chen BMC Plant Biology (2018) 18:98 Page 9 of 13 Fig. 6 Response of major flowering pathway genes in chilled flower buds (compared to nonchilled flower buds). The relationships among the listed genes are drawn according to the diagram for A. thaliana by Fornara et al. 2010 [7], although not all DE genes of blueberry align perfectly with the correlations proposed for A. thaliana. All the listed genes in this diagram showed down-regulation in late-pink buds (compared to chilled flower buds) (Additional file 2: Table S1) flowering which suggests that chilling requirement is The interaction of FD and FT promotes flowering in not replaceable by VcFT manipulation. When the re- A. thaliana while TFL1 is a negative regulator of FT [27, sults were aligned to the flowering pathway of A. 28]. In this study, reduced VcFD and VcTFL1 expres- thaliana [7], the acquired chill did not change VcFT sions did not changed VcFT expression in chilled flower expression in blueberry flower buds (Fig. 6). This re- buds. However, increased VcFD and decreased VcTFL1 sult was also observed in woody pear but not in in late-pink buds were associated with increased VcFT. herbaceous lily [24]. The inactive VcFT expression in TFL1 was considered a repressor of both LFY and AP1 response to chilling may be one major reason that in A. thaliana until recent evidence suggest that TFL1 chilled flower buds remain dormant prior to exposure transcription was suppressed by AP1 but promoted by to bud-breaking temperatures since VcFT increased in LFY [29, 30]. For chilled blueberry flower buds, in- late-pink buds (Table 1). creased expression of VcAP1 and decreased expression Overexpression of VcFT in leaves promoted expression of VcLFY were associated with decreased expression of of downstream genes VcSOC1, VcFUL, VcAP1, and VcTFL1. This result is consistent with the recent report VcLFY [17]. In this study, expression of VcSOC1, VcFUL, about the interactions among these three genes [30]. and VcAP1 were up-regulated but VcLFY was repressed During flower bud break, repressed VcTFL1 associated regardless of VcFT expression in flower buds (Table 1: with decreased VcLFY and VcAP1 supports the theory Fig. 5). The expression of VcFT-downstream genes was that TFL1 represses LFY [30]. Although some DEGs of regulated independently of VcFT in chilled flower buds. flowering pathway genes in blueberry match the pro- Additionally, repressed VcLFY response in chilled flower posed interactions in A. thaliana [7], VcFD and VcTFL1 buds is similar to results observed in grapevine (Vitis vi- seem to be playing a more significant role in blueberry nifera)[26]. These results suggest that VcLFY supression (Fig. 4). may play a role in chilling-mediated flowering by main- FLC interacts with SVP and both are repressed during taining bud dormancy before bud break. vernalization in A. thaliana [19]. In chilled blueberry Song and Chen BMC Plant Biology (2018) 18:98 Page 10 of 13 flower buds, both VcFLC and VcSVP homologues tolerance, and chilling-mediated flowering in blueberries. showed decreased expression but increased in late-pink For example, chilled blueberry buds showed higher buds (Table 1). In A. thaliana, ARP6 activates FLC, freezing tolerance than nonchilled buds and flower tis- MAF4, and MAF5, which are all repressors of plant sue [37]. Increased DE orthologues of ethylene genes in flowering in vernalization pathway [31]. In blueberry, the chilled blueberry buds were responsible for the en- decreased VcARP6 in chilled flower buds was not associ- hanced freezing tolerance in chilled buds while decreased ated with decreased expression of VcFLC, VcMAF4,or expression of DE ethylene orthologues in late-pink buds VcMAF5.However,increased VcARP6 (c49456_g2, Log FC reduced freezing tolerance (Fig. 3). = 8.1) was associated with an increase in these genes in late-pink buds (Table 1;Additional file 2: Table S1). Due Conclusions to ARP6’srole in A. thaliana vernalization, DE VcARP6 The changes from nonchilled to chilled and chilled to may significantly contribute to chilling-mediated flowering late-pink buds are associated with transcriptional of blueberryflowerbuds. changes in a large number of DE phytohormone-related genes and DE flowering pathway genes. The DE flower- Expression of blueberry MADS-box genes in chilled ing pathway genes suggest that orthologues of FT, FD, flower buds and late-pink buds TFL1, LFY, and MADS-box genes are the major genes The major flowering pathway genes SOC1, FLC, AP1, involved in chilling-mediated blueberry bud-break. The FUL, SVP, and AGL24 are MADS-box genes encoding DE phytohormone genes reveal the potential roles of MIKC (classical MIKC) proteins [7]. Similar to A. thali- phytohormone genes in cold acclimation, dormancy, ana, multiple blueberry MADS-box genes are present freezing tolerance, and chilling-mediated flowering in and activated at different flowering stages. VcSOC1, blueberries. The results contribute to the comprehensive VCAP1 and VcFUL are responsive to VcFT overexpres- investigation of the chilling-mediated flowering mechan- sion [17]. Additionally, constitutively expressed VcSOC1 ism in woody plants. or Keratin-like (K) domain of VcSOC1 promoted blue- berry flowering [32]. In this study, the functional ortho- Method logues of FLC and AGL24 were not detected in Plant materials blueberry, suggesting that the vernalization/chilling-me- The tetraploid southern highbush blueberry ‘Legacy’ needs diated flowering pathway of blueberry is different from over 800 chilling units (CU) for normal flowering. Twelve A. thaliana. VcSVP showed differential expression in 4-year old ‘Legacy’ plants were obtained through micro- chilled and late-pink buds (Table 1). propagation of in vitro cultured shoots. All plants were In woody fruit crops, functional FLC orthologues have grown in 4-gal pots in a secured courtyard under natural not been identified. The peach DAM genes mimic FLC light conditions at Michigan State University, East Lansing, response in A. thaliana under dormancy [10]. The DAM Michigan (latitude 42.701847, longitude − 84.482170). The genes are the orthologues of A. thaliana AGL24 and average low and high temperatures in January 2016 SVP genes [12, 33]. In this study, the DE DAMs showed were − 11 °C and − 1.8 °C, respectively (http://www.u- similarity to several MADS-box genes (VcAP1, VcSVP, sclimatedata.com/climate/east-lansing/michigan/united VcSOC1, and VcSPL3) (Table 1). Therefore, it is possible -states/usmi0248). In September 2015, six plants were that the interaction of multiple MADS-box genes moved to a heated greenhouse with a 12-h photo- co-regulates chilling-mediated flowering in blueberry as period and a minimum temperature of 23 °C in order well as other woody plants. to keep the plants from any chilling hour accumula- tion. The remaining six plants were kept in the se- Response of phytohormone genes during vernalization cured courtyard. In November, three plants were and devernalization selected from the greenhouse and 30–50 flower buds Phytohormones are involved in plant flowering and dor- were harvested per plant. These flower buds did not mancy. In A. thaliana, cold acclimation, dormancy, and receive any chilling temperatures and were labeled as plant flowering are affected by phytohormone gene ex- nonchilled flower buds. At the end of January 2016, pression [7, 34, 35]. The gibberellin pathway interacts three plants were selected from the courtyard and with the flowering pathway through SOC1 [7, 36]. Ethyl- 30–50 flower buds were harvested per plants. These ene signaling pathway genes are considered regulators flower buds experienced natural chilling conditions for freezing tolerance in A. thaliana. In this study, DE through mid-winter and were labeled as chilled flower phytohormone genes were identified in both chilled buds buds. In April, 20–30 flower buds per plant were ob- and late-pink buds of blueberry (Fig. 2; Additional file 3: tained from a second harvest of the same three plants Table S3). These DE phytohormone genes reveal their in the courtyard. These flower buds experienced nat- potential roles in cold acclimation, dormancy, freezing ural chilling conditions and began to flower in early Song and Chen BMC Plant Biology (2018) 18:98 Page 11 of 13 spring. The buds selected were at early-pink-stage Identification of the selected pathway genes and were labeled as late-pink buds. All tissues col- Representative protein sequences of selected genes of A. lected were frozen immediately in liquid nitrogen and thaliana were downloaded from the TAIR server (https:// stored at − 80 °C. Three plants for each bud type www.arabidopsis.org/tools/bulk/sequences/index.jsp). The were used as the three biological replicates for tran- retrieved sequences were used to search the blueberry scriptome analysis. transcriptome reference (refTrinity) using the tblastn command of BLAST+. The resultant transcripts that show e-value lower than − 20 were used to screen the DE tran- RNA preparation, sequencing, and de novo transcriptome script list of nonchilled floral buds. assembly The blueberry floral genes identified in the previous Total RNA of each blueberry sample (from individual study [17] were used to analyze flowering pathway plants) was isolated from 200 mg of bud tissues using genes. The pathway genes of major phytohormones a separate CTAB method [38] and was purified using (gibberellin [41], abscisic acid [42], cytokinin/Arabi- RNeasy Mini Kit (Qiagen, Valencia, CA, USA). All RNA dopsis Responsive Regulator [43], indole-3-acetic acid samples were purified using On-Column DNase digestion [44], and ethylene [45]) in A. thaliana were retrieved with the RNase-free DNase Set (Qiagen). The integrity of from TAIR_10 server based on published gene iden- the RNA samples was assessed using the Agilent RNA 6000 tities (Additional file 3: Table S3). Additionally, se- Pico Kit (Agilent Technologies, Inc., Germany). All samples quences of A. thaliana MADS-box proteins were had an RNA quality score greater than 8.0 prior to submis- used to analyze blueberry MADS-box genes. Percent- sion for sequencing (100-bp pair end reads) using the Illu- ages of DE phyotohormone genes were calculated mina HiSeq2500 platform at the Research Technology based either on the number of orthologues to A. Support Facility at Michigan State University (East Lansing, thaliana genes or on the number of DE blueberry Michigan, USA). The FastQC program (www.bioinforma- transcripts. tics.babraham.ac.uk/projects/fastqc/) was used to assess the quality of sequencing reads for the per base quality scores RT-PCR of DE transcripts ranging from 30 to 40. Reliability of DE genes or transcripts identified through RNA-seq was evaluated through qRT-PCR analysis of six Differential expression analysis and transcriptome selected transcripts (Additional file 4: Table S4). These annotation transcripts are from the representative DE genes in RNA-seq reads of three biological replicates for non- auxin, ethylene, cytokinin, and gibberellin pathways. chilled, chilled, and late-pink buds were analyzed. Two They have high fold changes (> 2) and sequence specifi- technical replicates were sequenced for each biological city (based on alignment result of different isoforms) for replicate and were combined together for analysis. The PCR amplification. Eukaryotic translation initiation paired reads, two sets for each biological replicate, were factor 3 subunit H was the internal control (Additional aligned to the transcriptome reference Reftrinity file 4: Table S4). developed for ‘Legacy’ [17] and the abundance of each The same RNA samples used for RNA-sequencing, read was estimated using the Trinity command including samples of three biological replicates, were “align_and_estimate_abundance.pl”. The Trinity com- used for cDNA preparation. Reverse transcription of mand “run_DE_analysis.pl –method edgeR” was used RNA to cDNA was performed using SuperScript II re- for differential expression analysis. The DE transcripts verse transcriptase (Invitrogen, Carlsbad, CA, USA). The with false discovery rate (FDR) values below 0.05 were resulting cDNA of one micro gram of RNA was diluted used for further analyses. Comparison of transcriptome (volume 1: 4) in water and a 1 μl/sample (25 ng) was in chilled to nonchilled flower buds resulted in DE used for PCR reactions. transcripts/genes in chilled flower buds. Comparison of Integrated DNA Technologies, Inc. (https://www.idtdna. transcriptome in late-pink buds to chilled flower buds com/Primerquest/Home/Index) provided the online tool resulted in DE transcripts/genes in late-pink buds. DE for primer design and synthesized the primers (Additional transcripts in chilled buds were annotated using Trinota- file 4: Table S4). Three qRT-PCR analyses were performed te_v2.0 (https://trinotate.github.io). on an Agilent Technologies Stratagene Mx3005P (Agilent Technologies, Santa Clara, CA) using the SYBR Green sys- Gene network construction tem (Life Technologies, Carlsbad, CA). In each 25 μl reac- Annotated transcripts were imported to Cytoscape 3.5.0 tion mixture, 25 ng of cDNA, 200 nM of primers, and under BiNGO’s default parameters with selected ontol- 12.5 μl of 2× SYBR Green master mix were included. The ogy file ‘GOSlim_Plants’ and selected organism A. thali- reaction conditions for all primer pairs were 95 °C for ana [39, 40]. 10 min, 40 cycles of 30 s at 95 °C, 60 s at 60 °C and 60 s at Song and Chen BMC Plant Biology (2018) 18:98 Page 12 of 13 72 °C, and followed by one cycle of 60 s at 95 °C, 30 s at Received: 5 January 2018 Accepted: 15 May 2018 55 °C and 30 s at 95 °C. The specificity of the amplification reaction for each primer pair was determined by the melt- ing curve. Transcript levels within samples were normalized References 1. Anderson JV. Advances in plant dormancy. Switzerland: Springer to the eukaryotic translation initiation factor 3 subunit H. International Publishing; 2015. Fold changes were calculated by -ΔΔCt = − [(Ct – GOI 2. Zinn KE, Tunc-Ozdemir M, Harper JF. Temperature stress and plant sexual Ct ) – (Ct – Ct ) ](n =3). nom tissue 1 GOI nom tissue 2 reproduction: uncovering the weakest links. J Exp Bot. 2010;61(7):1959–68. 3. Ouellet F, Charron J-B. Cold acclimation and freezing tolerance in plants. In: eLS. Chichester: John Wiley & Sons, Ltd; 2013. Additional files 4. Chuine IC, Bonhomme M, Legave J-M, De Cortázar-atauri I, Charrier G, Lacointe A, Améglio T. Can phenological models predict tree phenology accurately in the future? The unrevealed hurdle of endodormancy break. Additional file 1: Table S2. DE MADS-box genes in chilled flower buds Glob Chang Biol. 2016;22:17. (CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in 5. Luedeling E, Girvetz EH, Semenov MA, Brown PH. Climate change affects ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: winter chill for temperate fruit and nut trees. PLoS One. 2011;6(5):e20155. Log (LPB/CB). Except #N/A (no differential expression), all the rest are DE 6. Atkinson CJ, Brennan RM, Jones HG. Declining chilling and its impact on genes. (XLSX 23 kb) temperate perennial crops. Environ Exp Bot. 2013;91:48–62. Additional file 2: Table S1. DE floral genes in chilled flower buds (CB) 7. Fornara F, de Montaigu A, Coupland G. SnapShot: control of flowering in [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in Arabidopsis. Cell. 2010;141(3):550, 550.e1–2. ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: 8. Greenup A, Peacock WJ, Dennis ES, Trevaskis B. The molecular biology of Log (LPB/CB). #N/A: no differential expression. (XLSX 140 kb) seasonal flowering-responses in Arabidopsis and the cereals. Ann Bot. 2009; Additional file 3: Table S3. DE phytohormones in chilled flower 103(8):1165–72. buds(CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs 9. Higgins JA, Bailey PC, Laurie DA. Comparative genomics of flowering time CB) in ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink 2 pathways using Brachypodium distachyon as a model for the temperate buds: Log (LPB/CB). Except #N/A (no differential expression), all the rest 2 grasses. PLoS One. 2010;5(4):e10065. are DE genes. (XLSX 401 kb) 10. Bielenberg DG, Wang Y, Li ZG, Zhebentyayeva T, Fan SH, Reighard GL, Scorza R, Abbott AG. Sequencing and annotation of the evergrowing locus Additional file 4: Table S4. Primers used in this study. (DOCX 54 kb) 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 Genet Genomes. 2008;4(3):495–507. Abbreviations 11. Wang Y, Georgi LL, Reighard GL, Scorza R, Abbott AG. Genetic mapping of ABA: abscisic acid; DE: Differentially expressed; FDR: False discovery rate; the evergrowing gene in peach [Prunus persica (L.) Batsch]. J Hered. 2002; GA: gibberellin; GO: Gene ontology; IAA: indole-3-acetic acid; qRT- 93(5):352–8. PCR: Quantitative reverse transcriptase polymerase chain reaction 12. Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, Akagi T, Tao R. Functional and expressional analyses of PmDAM genes associated with Acknowledgements endodormancy in Japanese apricot. Plant Physiol. 2011;157(1):485–97. The authors would thank Dr. Wayne H. Loescher for reviewing this 13. Jimenez S, Reighard GL, Bielenberg DG. Gene expression of DAM5 and manuscript, Dr. Jeff Landgraf and Mr. Kevin Carr at Michigan State University DAM6 is suppressed by chilling temperatures and inversely correlated with Research Technology Support Facility for RNA sequencing. This research is bud break rate. Plant Mol Biol. 2010;73(1–2):157–67. partially supported by AgBioResearch Project GREEEN of Michigan State 14. Wilkie JD, Sedgley M, Olesen T. Regulation of floral initiation in horticultural University (http://www.canr.msu.edu/research/plant-agriculture/ trees. J Exp Bot. 2008;59(12):3215–28. project_greeen/). 15. Ehlenfeldt MK, Prior RL. Oxygen radical absorbance capacity (ORAC) and phenolic and anthocyanin concentrations in fruit and leaf tissues of highbus blueberry. J Agr Food Chem. 2001;49(5):2222–7. Funding 16. Song GQ, Walworth A, Zhao DY, Jiang N, Hancock JF. The Vaccinium This research was partially supported by AgBioResearch of Michigan State corymbosum FLOWERING LOCUS T-like gene (VcFT): a flowering activator University (http://agbioresearch.msu.edu/programs/info/project_greeen). reverses photoperiodic and chilling requirements in blueberry. Plant Cell Rep. 2013;32(11):1759–69. Availability of data and materials 17. Walworth AE, Chai B, Song GQ. Transcript profile of flowering regulatory Our blueberry transcriptome reference Reftrinity has been deposited in genes in VcFT-overexpressing blueberry plants. PLoS One. 2016;11(6): GenBank (Accession number: SRX2728597). Datasets from the current study e0156993. are available from the corresponding author on request. 18. Gao X, Walworth AE, Mackie C, Song GQ. Overexpression of blueberry FLOWERING LOCUS T is associated with changes in the expression of Authors’ contributions phytohormone-related genes in blueberry plants. Hortic Res. 2016;3:16053. GS conceived and supervised the study; QC and GS conducted the 19. Song GQ, Gao X. Transcriptomic changes reveal gene networks responding experiments; GS analyzed the data and wrote the manuscript. Both authors to the overexpression of a blueberry DWARF AND DELAYED FLOWERING 1 read and approved the manuscript. gene in transgenic blueberry plants. BMC Plant Biol. 2017;17(1):106. 20. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, et al. De novo transcript sequence Ethics approval and consent to participate reconstruction from RNA-seq using the trinity platform for reference Not applicable. generation and analysis. Nat Protoc. 2013;8(8):1494–512. 21. Mateos JL, Madrigal P, Tsuda K, Rawat V, Richter R, Romera-Branchat M, Competing interests Fornara F, Schneeberger K, Krajewski P, Coupland G. Combinatorial The authors declare that they have no competing interests. activities of SHORT VEGETATIVE PHASE and FLOWERING LOCUS C define distinct modes of flowering regulation in Arabidopsis. Genome Biol. 2015;16:31. Publisher’sNote 22. Ratcliffe OJ, Kumimoto RW, Wong BJ, Riechmann JL. Analysis of the Springer Nature remains neutral with regard to jurisdictional claims in published Arabidopsis MADS AFFECTING FLOWERING gene family: MAF2 prevents maps and institutional affiliations. Vernalization by Short periods of cold. Plant Cell. 2003;15(5):1159–69. Song and Chen BMC Plant Biology (2018) 18:98 Page 13 of 13 23. Wen Z, Guo W, Li J, Lin H, He C, Liu Y, Zhang Q, Liu W. Comparative 47. Corbineau F, Xia Q, Bailly C, El-Maarouf-Bouteau H. Ethylene, a key factor in Transcriptomic analysis of Vernalization- and Cytokinin-induced floral the regulation of seed dormancy. Front Plant Sci. 2014;5:539. transition in Dendrobium nobile. Sci Rep. 2017;7:45748. 48. Yamauchi Y, Ogawa M, Kuwahara A, Hanada A, Kamiya Y, Yamaguchi S. 24. Li W, Liu X, Lu Y. Transcriptome comparison reveals key candidate genes in Activation of gibberellin biosynthesis and response pathways by low response to vernalization of oriental lily. BMC Genomics. 2016;17:664. temperature during imbibition of Arabidopsis thaliana seeds. Plant Cell. 2004;16(2):367–78. 25. Rowland LJ, Alkharouf N, Darwish O, Ogden EL, Polashock JJ, Bassil NV, Main D. Generation and analysis of blueberry transcriptome sequences from 49. Muller B, Sheen J. Advances in cytokinin signaling. Science. 2007; leaves, developing fruit, and flower buds from cold acclimation through 318(5847):68–9. 50. Greenham K, McClung CR. Integrating circadian dynamics with physiological deacclimation. BMC Plant Biol. 2012;12:46. processes in plants. Nat Rev Genet. 2015;16(10):598–610. 26. Carmona MJ, Cubas P, Martinez-Zapater JM. VFL, the grapevine 51. Velasquez SM, Barbez E, Kleine-Vehn J, Estevez JM. Auxin and cellular FLORICAULA/LEAFY ortholog, is expressed in meristematic regions elongation. Plant Physiol. 2016;170(3):1206–15. independently of their fate. Plant Physiol. 2002;130(1):68–77. 27. Abe M, Kobayashi Y, Yamamoto S, Daimon Y, Yamaguchi A, Ikeda Y, Ichinoki H, Notaguchi M, Goto K, Araki T. FD, a bZIP protein mediating signals from the floral pathway integrator FT at the shoot apex. Science. 2005;309(5737): 1052–6. 28. Wigge PA, Kim MC, Jaeger KE, Busch W, Schmid M, Lohmann JU, Weigel D. Integration of spatial and temporal information during floral induction in Arabidopsis. Science. 2005;309(5737):1056–9. 29. Liljegren SJ, Gustafson-Brown C, Pinyopich A, Ditta GS, Yanofsky MF. Interactions among APETALA1, LEAFY, and TERMINAL FLOWER1 specify meristem fate. Plant Cell. 1999;11(6):1007–18. 30. Goslin K, Zheng BB, Serrano-Mislata A, Rae L, Ryan PT, Kwasniewska K, Thomson B, O'Maoileidigh DS, Madueno F, Wellmer F, et al. Transcription factor interplay between LEAFY and APETALA1/CAULIFLOWER during floral initiation. Plant Physiol. 2017;174(2):1097–109. 31. Deal RB, Kandasamy MK, McKinney EC, Meagher RB. The nuclear actin- related protein ARP6 is a pleiotropic developmental regulator required for the maintenance of FLOWERING LOCUS C expression and repression of flowering in Arabidopsis. Plant Cell. 2005;17(10):2633–46. 32. Song GQ, Walworth A, Zhao DY, Hildebrandt B, Leasia M. Constitutive expression of the K-domain of a Vaccinium corymbosum SOC1-like (VcSOC1-K) MADS-box gene is sufficient to promote flowering in tobacco. Plant Cell Rep. 2013;32(11):1819–26. 33. Jimenez S, Lawton-Rauh AL, Reighard GL, Abbott AG, Bielenberg DG. Phylogenetic analysis and molecular evolution of the dormancy associated MADS-box genes from peach. BMC Plant Biol. 2009;9:81. 34. Shi Y, Ding Y, Yang S. Cold signal transduction and its interplay with phytohormones during cold acclimation. Plant Cell Physiol. 2015;56(1):7–15. 35. Kendall SL, Hellwege A, Marriot P, Whalley C, Graham IA, Penfield S. Induction of dormancy in Arabidopsis summer annuals requires parallel regulation of DOG1 and hormone metabolism by low temperature and CBF transcription factors. Plant Cell. 2011;23(7):2568–80. 36. El-Showk S, Ruonala R, Helariutta Y. Crossing paths: cytokinin signalling and crosstalk. Development. 2013;140(7):1373–83. 37. Walworth AE, Rowland LJ, Polashock JJ, Hancock JF, Song GQ. Overexpression of a blueberry-derived CBF gene enhances cold tolerance in a southern highbush blueberry cultivar. Mol Breed. 2012;30(3):1313–23. 38. Zamboni A, Pierantoni L, De Franceschi P. Total RNA extraction from strawberry tree (Arbutus unedo) and several other woodyplants. Iforest. 2008;1:122–5. 39. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11): 2498–504. 40. Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9. 41. Regnault T, Daviere JM, Heintz D, Lange T, Achard P. The gibberellin biosynthetic genes AtKAO1 and AtKAO2 have overlapping roles throughout Arabidopsis development. Plant J. 2014;80(3):462–74. 42. Xiong L, Zhu JK. Regulation of abscisic acid biosynthesis. Plant Physiol. 2003; 133(1):29–36. 43. Hwang D, Chen HC, Sheen J. Two-component signal transduction pathways in Arabidopsis. Plant Physiol. 2002;129(2):500–15. 44. Normanly J, Bartel B. Redundancy as a way of life - IAA metabolism. Curr Opin Plant Biol. 1999;2(3):207–13. 45. Wang KL, Li H, Ecker JR. Ethylene biosynthesis and signaling networks. Plant Cell. 2002;14(Suppl):S131–51. 46. Dong T, Park Y, Hwang I. Abscisic acid: biosynthesis, inactivation, homoeostasis and signalling. Essays Biochem. 2015;58:29–48. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png BMC Plant Biology Springer Journals

Comparative transcriptome analysis of nonchilled, chilled, and late-pink bud reveals flowering pathway genes involved in chilling-mediated flowering in blueberry

Free
13 pages

Loading next page...
 
/lp/springer_journal/comparative-transcriptome-analysis-of-nonchilled-chilled-and-late-pink-mdXQhSeUIx
Publisher
Springer Journals
Copyright
Copyright © 2018 by The Author(s).
Subject
Life Sciences; Plant Sciences; Agriculture; Tree Biology
eISSN
1471-2229
D.O.I.
10.1186/s12870-018-1311-8
Publisher site
See Article on Publisher Site

Abstract

Background: Blueberry cultivars require a fixed quantity of chilling hours during winter endo-dormancy for vernalization. In this study, transcriptome analysis using RNA sequencing data from nonchilled, chilled, and late pink buds of southern highbush blueberry ‘Legacy’ was performed to reveal genes associated with chilling accumulation and bud break. Results: Fully chilled ‘Legacy’ plants flowered normally whereas nonchilled plants could not flower. Compared to nonchilled flower buds, chilled flower buds showed differential expression of 89% of flowering pathway genes, 86% of MADS-box genes, and 84% of cold-regulated genes. Blueberry orthologues of FLOWERING LOCUS T (FT) did not show a differential expression in chilled flower buds (compared to nonchilled flower bud) but were up-regulated in late-pink buds (compared to chilled flower bud). Orthologoues of major MADS-box genes were significantly up- regulated in chilled flower buds and down-regulated in late-pink buds. Functional orthologues of FLOWERING LOCUS C (FLC) were not found in blueberry. Orthologues of Protein FD (FD), TERMINAL FLOWER 1 (TFL1),and LEAFY (LFY) were down-regulated in chilled flower buds and in late-pink buds compared to nonchilled flower bud. Conclusions: The changes from nonchilled to chilled and chilled to late-pink buds are associated with transcriptional changes in a large number of differentially expressed (DE) phytohormone-related genes and DE flowering pathway genes. The profile of DE genes suggests that orthologues of FT, FD, TFL1, LFY, and MADS-box genes are the major genes involved in chilling-mediated blueberry bud-break. The results contribute to the comprehensive investigation of the vernalization-mediated flowering mechanism in woody plants. Keywords: Chilling requirement, Cold hardness, Flowering time control, Freezing tolerance, Vaccinium corymbosum, Vernalization, Woody plant Background Climate change in the last 40 years has caused earlier Winter dormancy (endo-dormancy) is essential for de- shifts in the onset of the growing season for trees (e.g., ciduous fruit crop survival [1, 2]. Under inductive low 2.3 days/decade in temperate Europe) and increased temperatures in the fall, deciduous woody fruit and nut temperature fluctuation [4]. Early onset of the growing crops are acclimated to develop freezing tolerance; season causes insufficient chilling hours and prevents meanwhile, accumulation of effective chilling hours is bud-break in fruit trees. Increased temperature fluctu- stimulated [3]. Sufficient chilling accumulation gives ation during plant bloom turns early season frosts into a plants full vernalization, which is a prerequisite for danger, with freezing injuries to flowers and young fruits bud-break in the spring. [5]. Plant breeding to manipulate chilling requirements and develop improved freeze tolerant cultivars are con- sidered to be long-term solutions to mitigate reduced * Correspondence: songg@msu.edu Plant Biotechnology Resource and Outreach Center, Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA © The Author(s). 2018 Open Access This 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. Song and Chen BMC Plant Biology (2018) 18:98 Page 2 of 13 winter chill, decrease freezing damages, and secure de- natural winter conditions flowered normally (Fig. 1). Tran- ciduous fruit production [6]. scriptome comparison of the chilled to nonchilled flower Seasonal flowering plays a significant role in a plant’slife buds using Trinity and a non-annotated transcriptome cycle and is controlled by a network of flowering pathway reference (Reftrinity) of tetraploid blueberry in GenBank genes [7–9]. FLOWERING LOCUS C (FLC) is a key (Accession number: SRX2728597) [19, 20] revealed 37,000 regulator in the vernalization pathway of winter-annual differentially expressed genes (DEGs) and 47,000 DE Arabidopsis thaliana ecotypes [7]. In winter wheat and transcripts. barley, VERNALIZATION2 (VRN2) is a major regulator of To conduct transcriptome analysis, Reftrinity (180,000 vernalization-mediated flowering [8]. FLC and VRN2 ana- genes and 250,000 isoforms) was annotated using Trino- logs are recorded in peach (Prunus persica). These analogs tate [20]. The annotation resulted in 14,000 known are the DORMANCY ASSOCIATED MADS-box (DAM) genes from 30,000 genes and 55,000 isoforms of Reftri- genes which are a cluster of six MADS-box transcription nity. With this annotated reference, 64% of known blue- factors. The loss of all or part of the DAMs resulted in the berry genes and isoforms showed differential expression non-vernalized peach evergrowing mutant [10, 11]. The in the comparison of chilled to nonchilled flower buds DAM genes are considered alternatives to FLC in regulat- (herein referring to DEGs/DE transcripts in chilled ing vernalization-mediated chilling requirement and flow- flower buds). Chilling affected expression of numerous ering [10, 12]. However, the DAM genes show high genes simultaneously in blueberry flower buds. similarities to A. thaliana AGAMOUS-LIKE 24 (AGL24) and SHORT VEGETATIVE PHASE (SVP)genes [12, 13]. Effect of chilling accumulation on flowering pathway Additionally, functional analysis of DAMs to reveal their genes roles in chilling-mediated flowering through reverse gen- Differential expression was detected in chilled blueberry etics has not been reported in peach. To date, neither a flower buds for 89% of flowering pathway orthologues of functional FLC-LIKE nor a VRN2-LIKE gene has been known A. thaliana flowering pathway genes (Table 1; verified in Vaccinium plants [14]. Fig. 1). One of the top two orthologues of FT (VcFT: Blueberries and cranberries are the most important hereafter Vc before any A. thaliana gene refers to the Vaccinium fruits due to their high antioxidant and blueberry orthologue gene) showed a slight decrease anti-inflammatory capacities [15]. Deciphering the while the other showed a slight increase. Similar results mechanism of vernalization/chilling-mediated flowering were observed for CONSTANS-LIKE2 orthologues will facilitate molecular breeding of blueberry cultivars (Table 1). for low-chilling requirement. To investigate flowering Orthologues of seven major MADS-box genes in flower- responses under nonconducive conditions, functional ing pathway genes of blueberry buds, APETALA1 (AP1), analysis of a blueberry FLOWERING LOCUS T (VcFT) FRUITFULL (FUL), SUPPRESSOR OF OVEREXPRESSION gene has been conducted in highbush blueberry (Vacci- OF CONSTANS 1 (SOC1), CAULIFLOWER (CAL), FLC, nium corymbosum L.) [16–18]. Overexpression of VcFT AGL24,and SVP, were significantly up-regulated (Table 1). (about 2900-fold increase in leaf tissues) caused continu- The orthologues of Protein FD (FD), TERMINAL FLOWER ous and precocious flowering in in vitro shoots and in 1(TFL1),and LEAFY (LFY), and ACTIN-RELATED PRO- one-year old ‘Aurora’ plants [16]. However, Overexpres- TEIN6 (ARP6) were down-regulated genes. Both VcFLC sion of VcFT was not capable of fulfilling all chilling re- and VcSVP showed a contrasting response of FLC to quirements in blueberry [17, 18]. Over 80% of the flower vernalization in A. thaliana. The down-regulated VcFD and buds in two- and three-year old VcFT-overexpressing up-regulated VcSVP support the DE VcFT results [7]. How- plants could not bloom under greenhouse conditions ever, decreased expression of VcLFY contradicts the de- without a chilling period. To discover vernalization/chil- creased VcTFL1 and increased VcAP1, VcAGL24,and ling-responsive genes, transcriptome analyses were con- VcSOC1 expressions. The unchanged expression of VcFT ducted with blueberry flower buds. The profile of and decreased expression of VcLFY orthologues could pre- differentially expressed (DE) genes will facilitate our un- vent the flowering of chilled flower buds during plant derstanding of the role of vernalization/chilling in blue- vernalization. berry bud-break. The MADS-box transcription factor and flowering re- pressor, FLC, and FRIGIDA (FRI) are the major genes Results regulating the vernalization of A. thaliana. FRI activates Identification of DE transcripts in chilled flower buds FLC, but vernalization represses FLC [7]. Both VcFRI and Nonchilled flower buds from the southern highbush blue- VcFLC in blueberry are present in the annotated Reftri- berry ‘Legacy’ were grown in a heated greenhouse through nity. VcFRI showed a high similarity to FRI (e = − 98) while the winter season. These buds did not flower in the fol- two VcFLC transcripts showed a lower similarity to FLC lowing spring. In contrast, flower buds fully chilled under (e = − 31). These two VcFLC transcripts were annotated to Song and Chen BMC Plant Biology (2018) 18:98 Page 3 of 13 Fig. 1 Effect of vernalization of blueberry flowering. a Nonchilled flower buds could not bloom. b Early bloom of fully chilled flower buds. c Percentages of differentially expressed (DE) genes and transcripts in chilled flower buds (in comparison to nonchilled flower buds). *DE genes including A. thaliana genes and genes annotated to other species SEP1 of A. thaliana and MADS6 of rice, respectively late-pink buds to those in chilled flower buds (hereafter (Table 1). In chilled flower buds, VcFRI did not show dif- referring to DEGs/DE transcripts in late-pink buds) re- ferential expression whereas expression of VcFLC was in- sulted in 28,000 DE isoforms, which were annotated to creased about 25-fold for the top-two VcFLC candidate 11,000 known genes. genes (Table 1). When all 26 potential VcFLC candidate genes (e < − 20) were included in the analysis, four candi- Expression of flowering pathway genes in late-pink buds dates showed a down-regulation with a fold-change of 0.6– The major DE flower pathway genes in chilled flower buds 0.8. The other 22 candidate genes showed an up-regulation were compared to those in late-pink buds (Additional file with an average of 4.9-fold (Additional file 1:Table S2).The 2: Table S1). The major flowering pathway genes VcFT, inconsistency between VcFLC and FLC response to chill- VcSOC1, VcAP1,VcFUL, VcFLC, VcSVP,and VcLFY showed ing/vernalization suggests that a different or more compli- decreased expression in late-pink buds. In A. thaliana, FT, cated mechanism is involved in vernalization-mediated SOC1, AP1, FUL,and LFY promote flowering while SVP flowering in tetraploid, blueberry plants. Additionally, all and FLC are flowering repressors [7, 21]. The expression other DEGs involved in known vernalization-mediated changes of VcFD and VcTFL1 in late-pink buds were flowering pathway were up-regulated in chilled flower buds similar to those in chilled flower buds (Table 1). The de- with the exception of VcARP6 (c49456_g2, Log FC = − 7.7) creased VcFD expression was related, at least in part, to (Additional file 2:Table S1). VcFT down-regulation. Decreased VcTFL1 expression in late-pink buds was associated with a down-regulation of Identification of DE transcripts in late-pink buds VcAP1 and VcLFY. Additionally, decreased VcSVP and Fully chilled blueberry flower buds remained dormant VcFLC expression were associated with a decrease in VcFT under chilling conditions in January until continuous expression. warm conditions in April drove dormancy release and bloom. To investigate the effect of dormancy release on Expression of MADS-box genes gene expression, RNA sequencing data was obtained Blueberry DE MADS-box and DAMs genes were identi- from late-pink buds. Comparative transcripts analysis in fied using A. thaliana MADS-box genes and Japanese Song and Chen BMC Plant Biology (2018) 18:98 Page 4 of 13 Table 1 Differentially expressed major flowering pathway genes in both chilled flower buds (CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: Log (LPB/CB) 2 2 Arabidopsis_ID|Vc_ LogFC for LogFC for FDR for FDR for e-value_to_blast Gene_name Annotation_by_ Query_ transcript_ID chilled buds late-pink buds chilled buds late-pink buds Trinotate PmDAMs AT5G60910.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 9E-44 FUL, AGL8 AGL11_ARATH AT1G69120.1|c92021_g1_i1 −0.6 −1.8 4E-02 5E-15 1E-80 AP1, AGL7 AGL8_SOLTU AT5G60910.1|c88116_g1_i1 −0.4 0.8 4E-02 1E-10 5E-91 FUL, AGL8 AGL8_SOLTU AT1G69120.1|c88116_g8_i1 3.2 −5.4 1E-16 2E-59 4E-94 AP1, AGL7 AGL8_SOLTU AT1G69120.1|c81830_g1_i1 3.5 −2.1 2E-42 4E-24 2E-52 AP1, AGL7 AGL9_PETHY AT1G69120.1|c88116_g7_i1 3.6 −2.7 3E-33 7E-41 3E-90 AP1, AGL7 AP1_SINAL AT5G15840.1|c77980_g1_i1 −0.5 0.7 2E-03 3E-10 1E-66 CO, FG COL2_ARATH AT5G15840.1|c76265_g1_i1 2.3 −2.5 2E-09 6E-14 4E-30 CO, FG COL2_ARATH AT5G60910.1|c79125_g1_i1 −0.3 −1.1 4E-02 2E-13 2E-35 FUL, AGL8 DEFA_ANTMA AT5G60910.1|c67980_g1_i1 3.4 −4.6 7E-46 2E-70 5E-36 FUL, AGL8 DEFA_ANTMA AT5G61850.1|c96427_g2_i1 −2.4 −5.7 4E-13 4E-24 7E-102 LFY, LFY3 FLO_ANTMA AT5G61850.1|c96427_g2_i2 −2.1 −6.7 2E-13 9E-31 3E-98 LFY, LFY3 FLO_ANTMA AT1G65480.1|c84088_g2_i3 −0.9 −4.5 2E-02 7E-12 5E-95 FT HD3A_ORYSJ AT1G65480.1|c84088_g2_i5 1.1 −4.9 1E-03 9E-28 1E-94 FT HD3A_ORYSJ AT1G26310.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 2E-56 CAL, AGL10 MADS6_ORYSJ AT2G45660.1|c89673_g3_i1 2.6 −4.0 1E-09 3E-15 3E-53 SOC1, AGL20 SOC1_ARATH AT2G45660.1|c86010_g2_i1 3.9 −4.4 1E-08 2E-09 2E-54 SOC1, AGL20 SOC1_ARATH AT2G22540.1|c91377_g1_i7 1.0 −1.6 3E-02 1E-03 5E-69 SVP, AGL22 SVP_ARATH AT2G22540.1|c90829_g2_i1 1.2 −1.2 1E-03 5E-04 4E-77 SVP, AGL22 SVP_ARATH AT4G24540.1|c91377_g1_i7 1.0 −1.6 3E-02 1E-03 1E-39 AGL24 SVP_ARATH AT5G10140.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 2E-31 FLC, FLF, MADS6_ORYSJ AGL25 AT5G10140.1|c88116_g6_i1 4.1 −1.9 5E-91 1E-26 7E-31 FLC, FLF, SEP1_ARATH AGL25 AT4G18960.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 1E-72 MAF2, AGL31 AGL11_ARATH PmDAM1 AT2G14210.1|c80388_g1_i2 3.9 3.2 4E-02 5E-09 3E-79 MAF2, AGL31 MAD23_ORYSJ PmDAM1 AT5G60910.1|c77146_g1_i1 4.9 −3.1 2E-57 9E-37 7E-54 MAF2, AGL31 MADS6_ORYSJ PmDAM1 AT4G22950.1|c77424_g2_i2 1.0 2.2 6E-03 3E-19 1E-42 MAF3, AGL70 AGL11_ARATH PmDAM1 AT4G18960.1|c88293_g3_i1 4.7 −4.7 5E-21 5E-20 5E-59 MAF4 AG_TOBAC AT5G15800.2|c88116_g6_i1 4.1 −1.9 5E-91 1E-26 2E-87 MAF4 SEP1_ARATH PmDAM1 AT1G24260.2|c68983_g1_i1 2.5 0.0 3E-06 0E + 00 3E-42 MAF5, AGL68 . AT5G15800.2|c81830_g1_i1 3.5 −2.1 2E-42 4E-24 3E-89 MAF5, AGL68 AGL9_PETHY PmDAM1 #N/A: No differential expression apricot (Prunus mume) DAMs (PmDAMs) (Additional VcSOC1 and VcFLC could have multiple functions in file 1: Table S2; Table 1). Orthologues of 62 A. thaliana blueberry. MADS-box gene were identified in blueberry (unpub- Three DE PmDAM orthologues (VcPmDAM1, lished data). DE orthologues of 53 and 44 MADS-box VcPmDAM2, and VcPmDAM5) were identified and showed genes were detected in chilled flower buds and late-pink high similarities to four A. thaliana MADS-box genes, buds, respectively. These orthologues include the major MADS AFFECTING FLOWERING 2 (MAF2), MAF4, flowering pathway genes VcFLC, VcSOC1, VcSVP, VcAP1, MAF5 and FOREVER YOUNG FLOWER.In A. thaliana, VcFUL, VcCAL,and VcAGL24 (Table 1). The annotated MAF2, MAF4,and MAF5 are FLC paralogs. MAF2, MAF5 VcSOC1 (c86010_g2_i1) showed high similarities to 25 and FLC are down-regulated and MAF5 is up-regulated A. thaliana MADS-box genes. Similarly, the annotated during vernalization [22]. In contrast, the blueberry ortho- VcFLC homologues were similar to 23 A. thaliana logues VcFLC, VcMAF1, VcMAF2, VcMAF4, and VcMAF5 MADS-box genes (Table 1). The results suggest that were up-regulated while VcMAF3 was repressed in chilled Song and Chen BMC Plant Biology (2018) 18:98 Page 5 of 13 flower buds (Table 1; Additional file 1: Table S2). Add- phyotohormone orthologues than the chilled buds itionally, three DE VcPmDAMs were annotated to the (Fig. 2). TheDE phytohormonegenes suggest thepo- homologues VcSOC1, VcSVP, VcAP1,and VcSEP1.The tential involvement of these phytohormones during up-regulated VcPmDAM1 homologues were the only DE chilling and flowering. orthologue in chilled flower buds. In late-pink buds, 75% Due to the tetraploid nature of ‘Legacy’, orthologues of of DE VcPmDAM1 homologues and all DE VcPmDAM5 each A. thaliana gene used for query often have more homologues were down-regulated while DE VcPmDAM2 than one homologue (Fig. 2). Thus, an average of the Fold Change was up-regulated (Additional file 1:Table S2).These blue- changes (Log ) for all the DE transcripts that berry MADS-box genes showed significant changes in re- are orthologues derived from a single A. thaliana query sponse to both chilling and flower bud breaking (Table 1). gene was used to represent the overall change of each However, during vernalization, the responses of VcFLC, phytohormone-related gene (Fig. 3). Increased expres- VcMAFs, and VcPmDAMs diverges from FLC’s response sion of ABA1, ABA2 and NINE-CIS-EPOXYCAROTE- to vernalization in A. thaliana. NOID DIOXYGENASE 3 (NCED3) in the ABA biosynthesis pathway were seen in chilled flower buds. Response of phytohormone-related genes in chilled and ABA1 and ABA2 continued to increase and NCED3 de- late-pink buds creased in late-pink buds (Fig. 3). The increased expres- For both chilled buds and late-pink buds, DE transcripts sion of these orthologues indicates a potential increase showed high similarities to the pathway genes for five in ABA biosynthesis during vernalization. Regardless of major phytohormones (Additional file 3: Table S3). Over decreased ABA1 expression, increased NCED3 expres- 50% of the DE blueberry orthologues were related to sion suggests that there is an increase in ABA biosyn- abscisic acid (ABA), ethylene, auxin, and gibberellin thesis during floral bud break (Fig. 3). (GA) genes while 25% were related to cytokinin genes DE orthologues of ethylene signaling pathway genes (Fig. 2). The late-pink bud showed more DE were all up-regulated in chilled buds and down-regulated in flowering buds (Fig. 3). These orthologues are consid- ered regulators for freezing tolerance in A. thaliana. Indole-3-acetic acid (IAA) and GA biosynthesis pathway orthologues were up-regulated in chilled flower buds and decreased in late-pink buds (Fig. 3). GIBBERELLIN 3-BETA-DIOXYGENASE 2 (GA3OX2) in the GA pathway and AUX1 in the IAA pathway were major DEGs with high expression changes. The DE ARABIDOPSIS RE- SPONSE REGULATORS (ARRs) orthologues included two A-type, two B-type and five ARR-like genes in chilled buds and flowering buds. One B-type orthologue (ARR10) was suppressed only in chilled buds (Fig. 3). Gene networks of DEGs in chilled flower buds and late-pink buds Over-represented Gene Ontology (GO) terms (P <0.05) were grouped to visualize gene networks of the annotated DE transcripts using the GOslim_Plant as the selected GO file and A. thaliana annotation as the reference. The DE transcripts were classified in 70 and 73 over-represented GO terms for chilled flower buds and late-pink buds, respectively (Fig. 4). The over-represented GO terms for chilled flower buds and late-pink buds were identical except for two GO terms (Fig. 4), suggesting that the same tran- Fig. 2 Response of phytohormone-related genes and transcripts in scripts responded to temperature changes in these buds. chilled flower buds (vs. nonchilled flower buds) and late-pink buds (vs. Thedifferencein “biological_process” was two additional chilled flower buds). a Percentage of differentially expressed (DE) orthologues of A. thaliana genes (The number of DE A. thaliana genes over-represented GO terms (GO:0007610-behavior and ÷ total number of A. thaliana genes × 100). b Percentage of DE GO:040029-regulation of gene expression, epigenetic) transcripts of transcripts of blueberry phytohormone-related genes found in late-pink buds but not in chilled flower (The number of DE transcripts ÷ total number of transcripts × 100) buds. (Fig. 4). Song and Chen BMC Plant Biology (2018) 18:98 Page 6 of 13 Fig. 3 Average fold changes (Log FC) of differentially expressed homologues for each of the phytohormone-related orthologue of A. thaliana in chilled flower buds (CB) [vs. nonchilled flower buds (NB)] and late-pink buds (LPB) [vs. chilled flower buds (CB)]. LogFC for chilled buds: Log (CB/ NB). LogFC for late-pink buds: Log (LPB/CB). a Abscisic acid biosynthesis pathway genes [46]. b Ethylene biosynthesis and signaling pathway genes [45, 47]. c Gibberellin biosynthesis pathway genes [48]. d Two-component ARABIDOPSIS RESPONSE REGULATORS (ARR) [49, 50]. e Auxin biosynthesis pathway genes [51]. The bars represent standard deviation Over-represented GO terms in chilled flower buds terms related to growth, response to stress, and revealed the impact of chilling on GO terms in three reproduction (Fig. 4). Thegenenetwork basedon categories (Fig. 4). The over-represented GO terms over-represented GO terms facilitate our understand- in “biological_process” revealed that vernalization/ ing of the role DEGs in both chilled flower buds and chilling affected expression of genes in multiple GO late-pink buds (Fig. 4). Song and Chen BMC Plant Biology (2018) 18:98 Page 7 of 13 cd Fig. 4 (See legend on next page.) Song and Chen BMC Plant Biology (2018) 18:98 Page 8 of 13 (See figure on previous page.) Fig. 4 Gene networks of differentially expressed genes (DEGs) in chilled flower buds and late-pink buds. DEGs in chilled flower buds were identified in comparison to nonchilled flower buds while DEGs in late-pink buds were identified in comparison to chilled flower buds. The ontology file of GOSlim_Plants in BiNGO was used to identify over-represented GO terms (P < 0.05). a Comparison of the gene network in chilled flower buds to late-pink buds; white nodes and black edges are present in both gene networks; red nodes and edges are present only in the chilled buds; and green nodes and edges are present only in the late-pink buds. The number in each circle is a GO identity number. A gene network in chilled flower buds (b “Biological_process” c “Cellular component” d “Molecular function”). I, II, and III in b show GO terms related to stress, plant growth, and reproduction, respectively Validation of the expression of selected genes For blueberries, expressed sequence tags have been gener- In chilled flower buds and late-pink buds, qRT-PCR were ated from blueberry flower buds [25]. However, compara- used to validate DE transcripts of VcFD, VcTFL1,and tive transcriptome analyses of different stage floral buds VcARP6 (Fig. 5). The results suggested high-reliability of have not been documented. the RNA-seq data. The roles of VcFT, VcLFY and VcARP6 in vernalization- Discussion mediated blueberry flowering Transcriptome analysis is an effective approach to study Overexpression of VcFT (expression level > 2000-fold flowering pathway genes [23, 24]. Using this approach, in leaf tissues) resulted in precocious flowering [16, DEGs in response to vernalization have been identified in 17]. However, the high expression level did not com- Japanese pear (Pyrus pyrifolia Nakai) and oriental lily [24]. pletely reverse the need for chilling for normal plant Fig. 5 Comparison of RNA sequencing and qRT-PCR analysis of three differentially expressed genes in (a) chilled flower buds (compared to nonchilled flower buds) (b) and late-pink buds (compared to chilled flower buds). Eukaryotic translation initiation factor 3 subunit H is the internal control. Log fold-change was calculated by -ΔΔCt = − [(Ct – Ct ) – (Ct – Ct ) ]. Average Log fold-change ± standard deviation of three 2 GOI nom tissue 1 GOI nom tissue 2 2 biological replicates. Significant average fold-change determined using a Student’s t-test is denoted. An asterisk (*) indicates p <0.001 Song and Chen BMC Plant Biology (2018) 18:98 Page 9 of 13 Fig. 6 Response of major flowering pathway genes in chilled flower buds (compared to nonchilled flower buds). The relationships among the listed genes are drawn according to the diagram for A. thaliana by Fornara et al. 2010 [7], although not all DE genes of blueberry align perfectly with the correlations proposed for A. thaliana. All the listed genes in this diagram showed down-regulation in late-pink buds (compared to chilled flower buds) (Additional file 2: Table S1) flowering which suggests that chilling requirement is The interaction of FD and FT promotes flowering in not replaceable by VcFT manipulation. When the re- A. thaliana while TFL1 is a negative regulator of FT [27, sults were aligned to the flowering pathway of A. 28]. In this study, reduced VcFD and VcTFL1 expres- thaliana [7], the acquired chill did not change VcFT sions did not changed VcFT expression in chilled flower expression in blueberry flower buds (Fig. 6). This re- buds. However, increased VcFD and decreased VcTFL1 sult was also observed in woody pear but not in in late-pink buds were associated with increased VcFT. herbaceous lily [24]. The inactive VcFT expression in TFL1 was considered a repressor of both LFY and AP1 response to chilling may be one major reason that in A. thaliana until recent evidence suggest that TFL1 chilled flower buds remain dormant prior to exposure transcription was suppressed by AP1 but promoted by to bud-breaking temperatures since VcFT increased in LFY [29, 30]. For chilled blueberry flower buds, in- late-pink buds (Table 1). creased expression of VcAP1 and decreased expression Overexpression of VcFT in leaves promoted expression of VcLFY were associated with decreased expression of of downstream genes VcSOC1, VcFUL, VcAP1, and VcTFL1. This result is consistent with the recent report VcLFY [17]. In this study, expression of VcSOC1, VcFUL, about the interactions among these three genes [30]. and VcAP1 were up-regulated but VcLFY was repressed During flower bud break, repressed VcTFL1 associated regardless of VcFT expression in flower buds (Table 1: with decreased VcLFY and VcAP1 supports the theory Fig. 5). The expression of VcFT-downstream genes was that TFL1 represses LFY [30]. Although some DEGs of regulated independently of VcFT in chilled flower buds. flowering pathway genes in blueberry match the pro- Additionally, repressed VcLFY response in chilled flower posed interactions in A. thaliana [7], VcFD and VcTFL1 buds is similar to results observed in grapevine (Vitis vi- seem to be playing a more significant role in blueberry nifera)[26]. These results suggest that VcLFY supression (Fig. 4). may play a role in chilling-mediated flowering by main- FLC interacts with SVP and both are repressed during taining bud dormancy before bud break. vernalization in A. thaliana [19]. In chilled blueberry Song and Chen BMC Plant Biology (2018) 18:98 Page 10 of 13 flower buds, both VcFLC and VcSVP homologues tolerance, and chilling-mediated flowering in blueberries. showed decreased expression but increased in late-pink For example, chilled blueberry buds showed higher buds (Table 1). In A. thaliana, ARP6 activates FLC, freezing tolerance than nonchilled buds and flower tis- MAF4, and MAF5, which are all repressors of plant sue [37]. Increased DE orthologues of ethylene genes in flowering in vernalization pathway [31]. In blueberry, the chilled blueberry buds were responsible for the en- decreased VcARP6 in chilled flower buds was not associ- hanced freezing tolerance in chilled buds while decreased ated with decreased expression of VcFLC, VcMAF4,or expression of DE ethylene orthologues in late-pink buds VcMAF5.However,increased VcARP6 (c49456_g2, Log FC reduced freezing tolerance (Fig. 3). = 8.1) was associated with an increase in these genes in late-pink buds (Table 1;Additional file 2: Table S1). Due Conclusions to ARP6’srole in A. thaliana vernalization, DE VcARP6 The changes from nonchilled to chilled and chilled to may significantly contribute to chilling-mediated flowering late-pink buds are associated with transcriptional of blueberryflowerbuds. changes in a large number of DE phytohormone-related genes and DE flowering pathway genes. The DE flower- Expression of blueberry MADS-box genes in chilled ing pathway genes suggest that orthologues of FT, FD, flower buds and late-pink buds TFL1, LFY, and MADS-box genes are the major genes The major flowering pathway genes SOC1, FLC, AP1, involved in chilling-mediated blueberry bud-break. The FUL, SVP, and AGL24 are MADS-box genes encoding DE phytohormone genes reveal the potential roles of MIKC (classical MIKC) proteins [7]. Similar to A. thali- phytohormone genes in cold acclimation, dormancy, ana, multiple blueberry MADS-box genes are present freezing tolerance, and chilling-mediated flowering in and activated at different flowering stages. VcSOC1, blueberries. The results contribute to the comprehensive VCAP1 and VcFUL are responsive to VcFT overexpres- investigation of the chilling-mediated flowering mechan- sion [17]. Additionally, constitutively expressed VcSOC1 ism in woody plants. or Keratin-like (K) domain of VcSOC1 promoted blue- berry flowering [32]. In this study, the functional ortho- Method logues of FLC and AGL24 were not detected in Plant materials blueberry, suggesting that the vernalization/chilling-me- The tetraploid southern highbush blueberry ‘Legacy’ needs diated flowering pathway of blueberry is different from over 800 chilling units (CU) for normal flowering. Twelve A. thaliana. VcSVP showed differential expression in 4-year old ‘Legacy’ plants were obtained through micro- chilled and late-pink buds (Table 1). propagation of in vitro cultured shoots. All plants were In woody fruit crops, functional FLC orthologues have grown in 4-gal pots in a secured courtyard under natural not been identified. The peach DAM genes mimic FLC light conditions at Michigan State University, East Lansing, response in A. thaliana under dormancy [10]. The DAM Michigan (latitude 42.701847, longitude − 84.482170). The genes are the orthologues of A. thaliana AGL24 and average low and high temperatures in January 2016 SVP genes [12, 33]. In this study, the DE DAMs showed were − 11 °C and − 1.8 °C, respectively (http://www.u- similarity to several MADS-box genes (VcAP1, VcSVP, sclimatedata.com/climate/east-lansing/michigan/united VcSOC1, and VcSPL3) (Table 1). Therefore, it is possible -states/usmi0248). In September 2015, six plants were that the interaction of multiple MADS-box genes moved to a heated greenhouse with a 12-h photo- co-regulates chilling-mediated flowering in blueberry as period and a minimum temperature of 23 °C in order well as other woody plants. to keep the plants from any chilling hour accumula- tion. The remaining six plants were kept in the se- Response of phytohormone genes during vernalization cured courtyard. In November, three plants were and devernalization selected from the greenhouse and 30–50 flower buds Phytohormones are involved in plant flowering and dor- were harvested per plant. These flower buds did not mancy. In A. thaliana, cold acclimation, dormancy, and receive any chilling temperatures and were labeled as plant flowering are affected by phytohormone gene ex- nonchilled flower buds. At the end of January 2016, pression [7, 34, 35]. The gibberellin pathway interacts three plants were selected from the courtyard and with the flowering pathway through SOC1 [7, 36]. Ethyl- 30–50 flower buds were harvested per plants. These ene signaling pathway genes are considered regulators flower buds experienced natural chilling conditions for freezing tolerance in A. thaliana. In this study, DE through mid-winter and were labeled as chilled flower phytohormone genes were identified in both chilled buds buds. In April, 20–30 flower buds per plant were ob- and late-pink buds of blueberry (Fig. 2; Additional file 3: tained from a second harvest of the same three plants Table S3). These DE phytohormone genes reveal their in the courtyard. These flower buds experienced nat- potential roles in cold acclimation, dormancy, freezing ural chilling conditions and began to flower in early Song and Chen BMC Plant Biology (2018) 18:98 Page 11 of 13 spring. The buds selected were at early-pink-stage Identification of the selected pathway genes and were labeled as late-pink buds. All tissues col- Representative protein sequences of selected genes of A. lected were frozen immediately in liquid nitrogen and thaliana were downloaded from the TAIR server (https:// stored at − 80 °C. Three plants for each bud type www.arabidopsis.org/tools/bulk/sequences/index.jsp). The were used as the three biological replicates for tran- retrieved sequences were used to search the blueberry scriptome analysis. transcriptome reference (refTrinity) using the tblastn command of BLAST+. The resultant transcripts that show e-value lower than − 20 were used to screen the DE tran- RNA preparation, sequencing, and de novo transcriptome script list of nonchilled floral buds. assembly The blueberry floral genes identified in the previous Total RNA of each blueberry sample (from individual study [17] were used to analyze flowering pathway plants) was isolated from 200 mg of bud tissues using genes. The pathway genes of major phytohormones a separate CTAB method [38] and was purified using (gibberellin [41], abscisic acid [42], cytokinin/Arabi- RNeasy Mini Kit (Qiagen, Valencia, CA, USA). All RNA dopsis Responsive Regulator [43], indole-3-acetic acid samples were purified using On-Column DNase digestion [44], and ethylene [45]) in A. thaliana were retrieved with the RNase-free DNase Set (Qiagen). The integrity of from TAIR_10 server based on published gene iden- the RNA samples was assessed using the Agilent RNA 6000 tities (Additional file 3: Table S3). Additionally, se- Pico Kit (Agilent Technologies, Inc., Germany). All samples quences of A. thaliana MADS-box proteins were had an RNA quality score greater than 8.0 prior to submis- used to analyze blueberry MADS-box genes. Percent- sion for sequencing (100-bp pair end reads) using the Illu- ages of DE phyotohormone genes were calculated mina HiSeq2500 platform at the Research Technology based either on the number of orthologues to A. Support Facility at Michigan State University (East Lansing, thaliana genes or on the number of DE blueberry Michigan, USA). The FastQC program (www.bioinforma- transcripts. tics.babraham.ac.uk/projects/fastqc/) was used to assess the quality of sequencing reads for the per base quality scores RT-PCR of DE transcripts ranging from 30 to 40. Reliability of DE genes or transcripts identified through RNA-seq was evaluated through qRT-PCR analysis of six Differential expression analysis and transcriptome selected transcripts (Additional file 4: Table S4). These annotation transcripts are from the representative DE genes in RNA-seq reads of three biological replicates for non- auxin, ethylene, cytokinin, and gibberellin pathways. chilled, chilled, and late-pink buds were analyzed. Two They have high fold changes (> 2) and sequence specifi- technical replicates were sequenced for each biological city (based on alignment result of different isoforms) for replicate and were combined together for analysis. The PCR amplification. Eukaryotic translation initiation paired reads, two sets for each biological replicate, were factor 3 subunit H was the internal control (Additional aligned to the transcriptome reference Reftrinity file 4: Table S4). developed for ‘Legacy’ [17] and the abundance of each The same RNA samples used for RNA-sequencing, read was estimated using the Trinity command including samples of three biological replicates, were “align_and_estimate_abundance.pl”. The Trinity com- used for cDNA preparation. Reverse transcription of mand “run_DE_analysis.pl –method edgeR” was used RNA to cDNA was performed using SuperScript II re- for differential expression analysis. The DE transcripts verse transcriptase (Invitrogen, Carlsbad, CA, USA). The with false discovery rate (FDR) values below 0.05 were resulting cDNA of one micro gram of RNA was diluted used for further analyses. Comparison of transcriptome (volume 1: 4) in water and a 1 μl/sample (25 ng) was in chilled to nonchilled flower buds resulted in DE used for PCR reactions. transcripts/genes in chilled flower buds. Comparison of Integrated DNA Technologies, Inc. (https://www.idtdna. transcriptome in late-pink buds to chilled flower buds com/Primerquest/Home/Index) provided the online tool resulted in DE transcripts/genes in late-pink buds. DE for primer design and synthesized the primers (Additional transcripts in chilled buds were annotated using Trinota- file 4: Table S4). Three qRT-PCR analyses were performed te_v2.0 (https://trinotate.github.io). on an Agilent Technologies Stratagene Mx3005P (Agilent Technologies, Santa Clara, CA) using the SYBR Green sys- Gene network construction tem (Life Technologies, Carlsbad, CA). In each 25 μl reac- Annotated transcripts were imported to Cytoscape 3.5.0 tion mixture, 25 ng of cDNA, 200 nM of primers, and under BiNGO’s default parameters with selected ontol- 12.5 μl of 2× SYBR Green master mix were included. The ogy file ‘GOSlim_Plants’ and selected organism A. thali- reaction conditions for all primer pairs were 95 °C for ana [39, 40]. 10 min, 40 cycles of 30 s at 95 °C, 60 s at 60 °C and 60 s at Song and Chen BMC Plant Biology (2018) 18:98 Page 12 of 13 72 °C, and followed by one cycle of 60 s at 95 °C, 30 s at Received: 5 January 2018 Accepted: 15 May 2018 55 °C and 30 s at 95 °C. The specificity of the amplification reaction for each primer pair was determined by the melt- ing curve. Transcript levels within samples were normalized References 1. Anderson JV. Advances in plant dormancy. Switzerland: Springer to the eukaryotic translation initiation factor 3 subunit H. International Publishing; 2015. Fold changes were calculated by -ΔΔCt = − [(Ct – GOI 2. Zinn KE, Tunc-Ozdemir M, Harper JF. Temperature stress and plant sexual Ct ) – (Ct – Ct ) ](n =3). nom tissue 1 GOI nom tissue 2 reproduction: uncovering the weakest links. J Exp Bot. 2010;61(7):1959–68. 3. Ouellet F, Charron J-B. Cold acclimation and freezing tolerance in plants. In: eLS. Chichester: John Wiley & Sons, Ltd; 2013. Additional files 4. Chuine IC, Bonhomme M, Legave J-M, De Cortázar-atauri I, Charrier G, Lacointe A, Améglio T. Can phenological models predict tree phenology accurately in the future? The unrevealed hurdle of endodormancy break. Additional file 1: Table S2. DE MADS-box genes in chilled flower buds Glob Chang Biol. 2016;22:17. (CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in 5. Luedeling E, Girvetz EH, Semenov MA, Brown PH. Climate change affects ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: winter chill for temperate fruit and nut trees. PLoS One. 2011;6(5):e20155. Log (LPB/CB). Except #N/A (no differential expression), all the rest are DE 6. Atkinson CJ, Brennan RM, Jones HG. Declining chilling and its impact on genes. (XLSX 23 kb) temperate perennial crops. Environ Exp Bot. 2013;91:48–62. Additional file 2: Table S1. DE floral genes in chilled flower buds (CB) 7. Fornara F, de Montaigu A, Coupland G. SnapShot: control of flowering in [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs CB) in Arabidopsis. Cell. 2010;141(3):550, 550.e1–2. ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink buds: 8. Greenup A, Peacock WJ, Dennis ES, Trevaskis B. The molecular biology of Log (LPB/CB). #N/A: no differential expression. (XLSX 140 kb) seasonal flowering-responses in Arabidopsis and the cereals. Ann Bot. 2009; Additional file 3: Table S3. DE phytohormones in chilled flower 103(8):1165–72. buds(CB) [vs nonchilled flower buds (NB)] and late-pink buds (LPB) (vs 9. Higgins JA, Bailey PC, Laurie DA. Comparative genomics of flowering time CB) in ‘Legacy’. LogFC for chilled buds: Log (CB/NB). LogFC for late-pink 2 pathways using Brachypodium distachyon as a model for the temperate buds: Log (LPB/CB). Except #N/A (no differential expression), all the rest 2 grasses. PLoS One. 2010;5(4):e10065. are DE genes. (XLSX 401 kb) 10. Bielenberg DG, Wang Y, Li ZG, Zhebentyayeva T, Fan SH, Reighard GL, Scorza R, Abbott AG. Sequencing and annotation of the evergrowing locus Additional file 4: Table S4. Primers used in this study. (DOCX 54 kb) 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 Genet Genomes. 2008;4(3):495–507. Abbreviations 11. Wang Y, Georgi LL, Reighard GL, Scorza R, Abbott AG. Genetic mapping of ABA: abscisic acid; DE: Differentially expressed; FDR: False discovery rate; the evergrowing gene in peach [Prunus persica (L.) Batsch]. J Hered. 2002; GA: gibberellin; GO: Gene ontology; IAA: indole-3-acetic acid; qRT- 93(5):352–8. PCR: Quantitative reverse transcriptase polymerase chain reaction 12. Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, Akagi T, Tao R. Functional and expressional analyses of PmDAM genes associated with Acknowledgements endodormancy in Japanese apricot. Plant Physiol. 2011;157(1):485–97. The authors would thank Dr. Wayne H. Loescher for reviewing this 13. Jimenez S, Reighard GL, Bielenberg DG. Gene expression of DAM5 and manuscript, Dr. Jeff Landgraf and Mr. Kevin Carr at Michigan State University DAM6 is suppressed by chilling temperatures and inversely correlated with Research Technology Support Facility for RNA sequencing. This research is bud break rate. Plant Mol Biol. 2010;73(1–2):157–67. partially supported by AgBioResearch Project GREEEN of Michigan State 14. Wilkie JD, Sedgley M, Olesen T. Regulation of floral initiation in horticultural University (http://www.canr.msu.edu/research/plant-agriculture/ trees. J Exp Bot. 2008;59(12):3215–28. project_greeen/). 15. Ehlenfeldt MK, Prior RL. Oxygen radical absorbance capacity (ORAC) and phenolic and anthocyanin concentrations in fruit and leaf tissues of highbus blueberry. J Agr Food Chem. 2001;49(5):2222–7. Funding 16. Song GQ, Walworth A, Zhao DY, Jiang N, Hancock JF. The Vaccinium This research was partially supported by AgBioResearch of Michigan State corymbosum FLOWERING LOCUS T-like gene (VcFT): a flowering activator University (http://agbioresearch.msu.edu/programs/info/project_greeen). reverses photoperiodic and chilling requirements in blueberry. Plant Cell Rep. 2013;32(11):1759–69. Availability of data and materials 17. Walworth AE, Chai B, Song GQ. Transcript profile of flowering regulatory Our blueberry transcriptome reference Reftrinity has been deposited in genes in VcFT-overexpressing blueberry plants. PLoS One. 2016;11(6): GenBank (Accession number: SRX2728597). Datasets from the current study e0156993. are available from the corresponding author on request. 18. Gao X, Walworth AE, Mackie C, Song GQ. Overexpression of blueberry FLOWERING LOCUS T is associated with changes in the expression of Authors’ contributions phytohormone-related genes in blueberry plants. Hortic Res. 2016;3:16053. GS conceived and supervised the study; QC and GS conducted the 19. Song GQ, Gao X. Transcriptomic changes reveal gene networks responding experiments; GS analyzed the data and wrote the manuscript. Both authors to the overexpression of a blueberry DWARF AND DELAYED FLOWERING 1 read and approved the manuscript. gene in transgenic blueberry plants. BMC Plant Biol. 2017;17(1):106. 20. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, et al. De novo transcript sequence Ethics approval and consent to participate reconstruction from RNA-seq using the trinity platform for reference Not applicable. generation and analysis. Nat Protoc. 2013;8(8):1494–512. 21. Mateos JL, Madrigal P, Tsuda K, Rawat V, Richter R, Romera-Branchat M, Competing interests Fornara F, Schneeberger K, Krajewski P, Coupland G. Combinatorial The authors declare that they have no competing interests. activities of SHORT VEGETATIVE PHASE and FLOWERING LOCUS C define distinct modes of flowering regulation in Arabidopsis. Genome Biol. 2015;16:31. Publisher’sNote 22. Ratcliffe OJ, Kumimoto RW, Wong BJ, Riechmann JL. Analysis of the Springer Nature remains neutral with regard to jurisdictional claims in published Arabidopsis MADS AFFECTING FLOWERING gene family: MAF2 prevents maps and institutional affiliations. Vernalization by Short periods of cold. Plant Cell. 2003;15(5):1159–69. Song and Chen BMC Plant Biology (2018) 18:98 Page 13 of 13 23. Wen Z, Guo W, Li J, Lin H, He C, Liu Y, Zhang Q, Liu W. Comparative 47. Corbineau F, Xia Q, Bailly C, El-Maarouf-Bouteau H. Ethylene, a key factor in Transcriptomic analysis of Vernalization- and Cytokinin-induced floral the regulation of seed dormancy. Front Plant Sci. 2014;5:539. transition in Dendrobium nobile. Sci Rep. 2017;7:45748. 48. Yamauchi Y, Ogawa M, Kuwahara A, Hanada A, Kamiya Y, Yamaguchi S. 24. Li W, Liu X, Lu Y. Transcriptome comparison reveals key candidate genes in Activation of gibberellin biosynthesis and response pathways by low response to vernalization of oriental lily. BMC Genomics. 2016;17:664. temperature during imbibition of Arabidopsis thaliana seeds. Plant Cell. 2004;16(2):367–78. 25. Rowland LJ, Alkharouf N, Darwish O, Ogden EL, Polashock JJ, Bassil NV, Main D. Generation and analysis of blueberry transcriptome sequences from 49. Muller B, Sheen J. Advances in cytokinin signaling. Science. 2007; leaves, developing fruit, and flower buds from cold acclimation through 318(5847):68–9. 50. Greenham K, McClung CR. Integrating circadian dynamics with physiological deacclimation. BMC Plant Biol. 2012;12:46. processes in plants. Nat Rev Genet. 2015;16(10):598–610. 26. Carmona MJ, Cubas P, Martinez-Zapater JM. VFL, the grapevine 51. Velasquez SM, Barbez E, Kleine-Vehn J, Estevez JM. Auxin and cellular FLORICAULA/LEAFY ortholog, is expressed in meristematic regions elongation. Plant Physiol. 2016;170(3):1206–15. independently of their fate. Plant Physiol. 2002;130(1):68–77. 27. Abe M, Kobayashi Y, Yamamoto S, Daimon Y, Yamaguchi A, Ikeda Y, Ichinoki H, Notaguchi M, Goto K, Araki T. FD, a bZIP protein mediating signals from the floral pathway integrator FT at the shoot apex. Science. 2005;309(5737): 1052–6. 28. Wigge PA, Kim MC, Jaeger KE, Busch W, Schmid M, Lohmann JU, Weigel D. Integration of spatial and temporal information during floral induction in Arabidopsis. Science. 2005;309(5737):1056–9. 29. Liljegren SJ, Gustafson-Brown C, Pinyopich A, Ditta GS, Yanofsky MF. Interactions among APETALA1, LEAFY, and TERMINAL FLOWER1 specify meristem fate. Plant Cell. 1999;11(6):1007–18. 30. Goslin K, Zheng BB, Serrano-Mislata A, Rae L, Ryan PT, Kwasniewska K, Thomson B, O'Maoileidigh DS, Madueno F, Wellmer F, et al. Transcription factor interplay between LEAFY and APETALA1/CAULIFLOWER during floral initiation. Plant Physiol. 2017;174(2):1097–109. 31. Deal RB, Kandasamy MK, McKinney EC, Meagher RB. The nuclear actin- related protein ARP6 is a pleiotropic developmental regulator required for the maintenance of FLOWERING LOCUS C expression and repression of flowering in Arabidopsis. Plant Cell. 2005;17(10):2633–46. 32. Song GQ, Walworth A, Zhao DY, Hildebrandt B, Leasia M. Constitutive expression of the K-domain of a Vaccinium corymbosum SOC1-like (VcSOC1-K) MADS-box gene is sufficient to promote flowering in tobacco. Plant Cell Rep. 2013;32(11):1819–26. 33. Jimenez S, Lawton-Rauh AL, Reighard GL, Abbott AG, Bielenberg DG. Phylogenetic analysis and molecular evolution of the dormancy associated MADS-box genes from peach. BMC Plant Biol. 2009;9:81. 34. Shi Y, Ding Y, Yang S. Cold signal transduction and its interplay with phytohormones during cold acclimation. Plant Cell Physiol. 2015;56(1):7–15. 35. Kendall SL, Hellwege A, Marriot P, Whalley C, Graham IA, Penfield S. Induction of dormancy in Arabidopsis summer annuals requires parallel regulation of DOG1 and hormone metabolism by low temperature and CBF transcription factors. Plant Cell. 2011;23(7):2568–80. 36. El-Showk S, Ruonala R, Helariutta Y. Crossing paths: cytokinin signalling and crosstalk. Development. 2013;140(7):1373–83. 37. Walworth AE, Rowland LJ, Polashock JJ, Hancock JF, Song GQ. Overexpression of a blueberry-derived CBF gene enhances cold tolerance in a southern highbush blueberry cultivar. Mol Breed. 2012;30(3):1313–23. 38. Zamboni A, Pierantoni L, De Franceschi P. Total RNA extraction from strawberry tree (Arbutus unedo) and several other woodyplants. Iforest. 2008;1:122–5. 39. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11): 2498–504. 40. Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005;21(16):3448–9. 41. Regnault T, Daviere JM, Heintz D, Lange T, Achard P. The gibberellin biosynthetic genes AtKAO1 and AtKAO2 have overlapping roles throughout Arabidopsis development. Plant J. 2014;80(3):462–74. 42. Xiong L, Zhu JK. Regulation of abscisic acid biosynthesis. Plant Physiol. 2003; 133(1):29–36. 43. Hwang D, Chen HC, Sheen J. Two-component signal transduction pathways in Arabidopsis. Plant Physiol. 2002;129(2):500–15. 44. Normanly J, Bartel B. Redundancy as a way of life - IAA metabolism. Curr Opin Plant Biol. 1999;2(3):207–13. 45. Wang KL, Li H, Ecker JR. Ethylene biosynthesis and signaling networks. Plant Cell. 2002;14(Suppl):S131–51. 46. Dong T, Park Y, Hwang I. Abscisic acid: biosynthesis, inactivation, homoeostasis and signalling. Essays Biochem. 2015;58:29–48.

Journal

BMC Plant BiologySpringer Journals

Published: May 31, 2018

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off