- Research article
- Open Access
Crosstalk in the darkness: bulb vernalization activates meristem transition via circadian rhythm and photoperiodic pathway
BMC Plant Biology volume 20, Article number: 77 (2020)
Geophytes possess specialized storage organs - bulbs, tubers, corms or rhizomes, which allow their survival during unfovarable periods and provide energy support for sprouting and sexual and vegetative reproduction. Bulbing and flowering of the geophyte depend on the combined effects of the internal and external factors, especially temperature and photoperiod. Many geophytes are extensively used in agriculture, but mechanisms of regulation of their flowering and bulbing are still unclear.
Comparative morpho-physiological and transcriptome analyses and quantitative validation of gene expression shed light on the molecular regulation of the responses to vernalization in garlic, a typical bulbous plant. Long dark cold exposure of bulbs is a major cue for flowering and bulbing, and its interactions with the genetic makeup of the individual plant dictate the phenotypic expression during growth stage. Photoperiod signal is not involved in the initial nuclear and metabolic processes, but might play role in the later stages of development, flower stem elongation and bulbing. Vernalization for 12 weeks at 4 °C and planting in November resulted in flower initiation under short photoperiod in December–January, and early blooming and bulbing. In contrast, non-vernalized plants did not undergo meristem transition. Comparisons between vernalized and non-vernalized bulbs revealed ~ 14,000 differentially expressed genes.
Low temperatures stimulate a large cascades of molecular mechanisms in garlic, and a variety of flowering pathways operate together for the benefit of meristem transition, annual life cycle and viable reproduction results.The circadian clock appears to play a central role in the transition of the meristem from vegetative to reproductive stage in bulbous plant, serving as integrator of the low-temperature signals and the expression of the genes associated with vernalization, photoperiod and meristem transition. The reserved photoperiodic pathway is integrated at an upstream point, possibly by the same receptors. Therefore, in bulb, low temperatures stimulate cascades of developmental mechanisms, and several genetic flowering pathways intermix to achieve successful sexual and vegetative reproduction.
Perennial plant species carefully time their flowering with seasonal changes to ensure and maximize reproductive success. Flowering is induced, initiated and promoted both by internal and environmental cues, such as physiological age, gibberellin synthesis, photoperiod, and/or temperatures. The combined effect of two or more of these agents results in an integrated regulatory network that controls flowering time and quality [4, 89].
Despite the considerable genetic differences between and within species [9, 10, 18, 29], there is a great physiological similarity among plants from the temperate zone as many of them depend on vernalization for floral induction. Accumulation of chilling hours regulates essential changes in water status, in hormonal balance, in respiration and in carbohydrate mobilization , with a consequent increase in plant receptivity to changes in day length in the following spring .
Dependence on low temperatures evolved several times in the history of plants, with the consequent involvement of a number of regulatory mechanisms [4, 10, 40, 74]. For instance, Arabidopsis remains vegetative when FRIGIDA (FRI) upregulates the floral repressor FLOWERING LOCUS C (FLC), which in turn downregulates the floral integrators FT, FD, and SOC1. Cold induction downregulates FLC, thus enabling the expression of floral integrators and consequently the meristem transition from the vegetative to the reproductive state [44, 60, 79, 80]. An alternative regulating mechanism is common in the monocot cereals, where FT-homolog VERNALIZATION 3 (VRN3) is repressed by VRN2. Following cold induction, VRN1 represses VRN2 expression in the leaves, thus enabling the expression of VRN3 and the consequent transition of the meristem . In bulb onion, FT-like genes control the initiation of both bulbing and flowering . Under short photoperiod, high AcFT4 expression inhibits bulb formation by repressing AcFT1. This genotype-specific inhibitory effect gradually weakens when days elongate, and thereafter the inductive photoperiod downregulates AcFT4, with the consequent expression of bulb-promoting AcFT1. On the other hand, low temperatures in storage and/or in the field promote the up-regulation of AcFT2, which encodes for flowering in the spring/summer. Vernalization effects on the upregulation of FT-like sequences (LiFTL) were also reported for Lilium longiflorum [47, 49]. In tulip, TgFT2 is considered to act as florigen, whereas TgFT1 and TgFT3 may have bulb-specific functions .
Unlike vernalization, the molecular mechanism of the photoperiodic pathway evolved in early times and is well conserved in plants [4, 91]. It involves stabilizing/destabilizing balance between photoreceptors, with the consequent timely induction and initiation of flowering by the expression of CONSTANS (CO), which is regulated by both, GIGANTEA (GI) and the circadian rhythm [84, 91].
In many plant species, flowering induction and initiation require sequential and combined effects of vernalization and photoperiod. Hence a ‘memory’ of the vernalization effect is maintained by epigenetic mechanisms [8, 86]. In Arabidopsis, these chromatin modifications suppress the floral repressor FLC [16, 86], while in cereals they upregulate the floral activator VRN1 [67, 90]. It is generally accepted that for both mono- and dicots, FTs homologs act as the main genetic floral integrator of various pathways .
Geophytes are defined as plants with underground storage organs - bulbs, tubers, corms or rhizomes. In climates with environmental extremes, storage organs assist in survival during unfovarable periods and for providing energy support for sprouting and sexual and vegetative reproduction. Their life cycle consists of complex environment-dependent annual sequences, where induction and initiation of flowering and bulbing are tightly controlled by the surrounding conditions . Many geophyte species are extensively used in agriculture, and deep knowledge of their physiology is readily available, while main mechanisms of molecular regulation is still unclear .
The annual development cycle of typical bulbous geophyte garlic (Allium sativum L.), is well adapted to the environmental changes common in its center of origin in Central Asia . Hence, both bulbing and flowering are regulated by the combined effects of temperature and photoperiod [38, 58, 71, 72, 82, 94]. In bolting genotypes, low temperatures induce the transition of the meristem from vegetative to the reproductive state, with the consequent differentiation of inflorescence, scape elongation, and growth of axillary buds. In the early summer, when temperatures increase and days get longer, blooming occurs and mature compound bulbs enter dormancy.
Pre-planting cold treatment and/or winter temperatures significantly affect the developmental regulation and phenotypic expression of growing garlic plants, their bulbing and bulb quality, as well as its reproductive development . Extended vernalization strongly expedite vegetative growth, early meristem transition and bulbing, thus resulting in the development of only a few axillary buds (cloves) and consequently in small bulbs and inflorescences [5, 7, 26, 76, 95].
It is commonly accepted that photoperiod is the dominant environmental factor that controls stem elongation and bulbing in both onion and garlic [38, 70]. However, our recent findings show that garlic pre-planting storage at 4 °C for 12 weeks resulted in photoperiod-independent meristem transition already in mid-winter, in December–January, when day length is shortest . The present study aims to identify the molecular mechanisms involved in the regulation of responses to low temperature in a flowering garlic genotype. Furthermore, by comparing gene expression in plants exposed to long and short vernalization with non-vernalized plants, we obtained molecular evidence for flowering induction during the bulb cold storage.
Plants produced from the non-vernalized (NV) bulbs maintained a vegetative habitus and continuously produced new leaves throughout the growing period, up to 20 leaves in June . In these plants, the apical meristem did not shift to floral initiation, apical dominance did not weaken and lateral meristems did not differentiate into leaves or cloves. Dry weight of aboveground parts reached 17 g, as compared with 11–12 and 3-4 g for the plants raised after short (SV) and long vernalization (LV), respectively (Fig. 1a).
In comparison, plants raised from vernalized bulbs lost their apical dominance a number of weeks after planting, developed reproductive organs and cloves and yielded dormant bulbs (Fig. 1b). In plants raised after SV, the short vernalization did not suffice for induction of meristem transition prior to planting, but at the end of March, 140 days after planting (DAP), apical dominance weakened and reproductive meristem was observed. These plants produced 17–19 leaves prior to flowering in May, 185 DAP, and developed approximately 160 flowers per inflorescence. On harvest, 230 DAP, bulbs had 6–7 cloves, and mean bulb weight was 30 g (Fig. 1).
In LV plants, meristem transition was visible already early in January ~ 50 DAP, after the development of 7–8 leaves (Fig. 1). The inflorescences produced 100–120 flowers, and the bulbs had only 2–3 cloves. Bulb maturation occurred on 200 DAP, and mean bulb weight was 8 g (Fig. 1). Therefore, the cold units accumulated during the long vernalization sufficed for initiation of both flowering and bulbing in clone #87.
Global transcriptome analysis
RNA samples taken from clove apical meristems were used for the construction of pre-planting transcriptomes. The high-throughput parallel RNA-Seq libraries served for building cDNA libraries with 36–42 million 100-bp pair-end reads. Quality trimming and filtration resulted in 33–39 million clean reads. Together with the fertile garlic data (NCBI bioproject PRJNA243415; ) these reads were assembled using Trinity software, thus generated 112,388 trinity ‘genes’ with an average length of 1.178 bp and N50 of 1812 bp. About 90% of the clean reads were mapped to the A. sativum transcriptome catalog.
Comparisons between pairs of treatments revealed that only 4% of the detected ~ 14,000 differentially expressed genes (DEGs) were present in all comparisons. Comparisons LV/NV and LV/SV shared the highest number of common DEGs (38%) (Fig. 2a). Long cold exposure resulted in the highest number of differences in transcription activity, in both up- and down-regulated genes. Moreover, in all comparisons, the number of down-regulated genes was considerably higher than that of up-regulated ones (Fig. 2b).
Functional GO analysis of the 547 DEGs (4% of the total number) shared by all comparisons (Fig. 2a) revealed three main domains of activity - RNA methylation, chromatin organization, and biosynthesis of aromatic compounds (Fig. 3). These metabolic processes play general developmental roles towards the end of garlic dormancy, sprouting and meristem identity affected by vernalization.
Using GO algorithm, we analyzed the biological process functions of up- and down-regulated genes, in all comparisons (Fig. 4). For SV/NV comparison, the GO enrichment showed over-representation of the domains associated with hormone transport, response to temperature stimulus, and flavonoid metabolism (Fig. 4a, c), whereas chromatin organization, DNA packaging, regulation of gene expression, metabolic processes and biosynthesis, were down-regulated (Fig. 4b, d). The LV/NV comparison revealed strong over-represented domains associated with photosynthesis, carbon fixation, response to biotic stimulus and cytokinin, plastid organization, and light reaction (Fig. 4a, c). In contrast, function domains involved in cellular catabolism, defense response, and RNA metabolism were clearly down-regulated in LV cloves (Fig. 4b, d). Finally, the LV/SV comparison analysis showed upregulation of 2947 genes (Fig. 2a), associated with biosynthesis, chromatin organization, DNA modification and methylation, carbohydrate metabolism, photosynthesis and response to light, and biological processes related to photoperiodism, photosynthesis, response to external stimuli and auxin influx were abundant (Fig. 4a, c). The cellular component organization was extremely over-represented, whereas genes coding for a response to temperature stimulus and to RNA methylation were down-regulated (Fig. 4b).
Data mining for functional analysis
Our results indicate that differential growth performance of plants is associated with expression of genes involved in meristem transition. To verify this assumption, we performed data-mining in the DEGs list, searching for gene annotations in the databases using the keywords: flowering, meristem transition, photoperiod, circadian clock and vernalization. Further, we analysed the literature data on specific gene families or genes with higher complexity, e.g., FT family, in order to compare garlic with model and non-model plant species and assess the meristem identitiy pathways.
Consequently, we compiled a list of 74 DEGs associated with determination of meristem identity (Additional file 5: Table S2).
Gene co-expression networks (GCNs)
Gene co-expression networks were constructed as the integration of multiple expression datasets by connecting genes with similar patterns of expression across treatments to reveal the specific biological context of positive and negative gene linkages in the apical meristem. In positive GCN representation, at least three main modules were identified, and co-expression was found between genes involved in vernalization and photoperiod pathways, and in those annotated as genes involved in meristem transition (Fig. 5). Vernalization-associated genes (marked green) are present in all three modules and are closely linked with genes annotated to photoperiod response (marked blue) and to meristem identity genes (marked purple). In general, genes associated with a light reception, photoperiod response and circadian clock (Additional file 5: Table S2) represent a major part of the positive GCN and allied with the meristem transition network. A garlic homolog of SUPPRESSOR OF FRIGIDA 4 (SUF4), associated in Arabidopsis with activation of FLC expression, was found to be centrally positioned in the positive GCN. Hence, it seems to play a conserved role in integrating flowering-related pathways (Fig. 5).
Our analyses also revealed three FLOR1 (FLR1) homologs, associated with long photoperiod pathway of meristem transition. FLR1 homologs with slightly different sequences were present in each of the major modules, suggesting their possible involvement in flowering through alternative splicing. SALT TOLERANCE (STO), strongly associated with photoperiod response and flowering transition, is another gene of interest. It is located in the center of the largest GCN module, connecting two sub-modules.
The FT family of genes (FT1, FT2, FT4, and FT3/5) is centrally located in this network, but with no direct links connect between them. FT2 co-expresses with meristem-transition genes CCR2, EMB, COR27, while FT1 and FT4 co-express with genes involved in photoperiod response (CKA1, CKA2, LHY, and ADG1).
In the negative GCN, a vernalization-associated gene AGAMOUS-like 19 (AGL19) negatively correlates with photoperiodic and meristem identity genes (Additional file 3: Figure S3).
Differential expression of the “flowering” genes
For differential expression of selected candidate genes that are involved in flowering, we combined the results of the in silico analyses of 74 genes of interest (Additional file 5: Table S2) with qRT-PCR analysis of 15 putative genes, annotated for their involvement in various flowering pathways.
The genes whose expression is associated with vernalization are clustered in three large groups (Fig. 6a). One that contains homologs of SUF4, VIN3, and EMF2, shows low expression in NV and SV tissues but a strong expression in LV plants. A second cluster consists of genes slightly upregulated in LV, and the third contains VRN1, SOC1, TFL2 whose over-representation is clear in NV and SV plants, but they are downregulated in LV plants. qRT-PCR expression analysis of three genes assigned to this group, PRMT5, VIN3, and AGL19, show consistent over-representation in vernalized bulbs (Fig. 6b).
In silico analyses of genes associated with the photoperiodic pathway generated four clusters of the differential expression (Fig. 7a). One cluster includes genes over-represented in LV plants, with exceedingly high over-expression of VOZ1, known for its involvement in positive regulation of long-day photoperiodism, cold acclimation and flowering. Another cluster consists of up-regulated genes, especially FLR1.2, in SV plants, while LV resulted in down-regulation of these genes. The same cluster includes also a group of circadian rhythm genes (CKA1, CKA2, LHY, CCA, CUL4) that are over-represented following both NV and SV treatments and down-regulated after long vernalization. The other two clusters show moderate expression differences between treatments. Some are slightly down-regulated after prolonged vernalization (e.g., CKB1, ATAN11, CDF3, ADG1) (Fig. 7a). qRT-PCR expression analyses of four specific genes STO, LHY, JMJ18, and FKF, show their overexpression mainly in SV, and STO and JMJ18 responded similarly also to LV (Fig. 7b).
Microscopic examination of apical meristems on the day storage treatment ended clearly show a vegetative appearence, irrespective of treatment [data not shown]. Transcriptome analyses, however, indicates high expression of numerous genes, annotated as active in meristem identity. The heatmap consists of at least four distinct clusters of genes (Fig. 8a). Most DEGs associated with meristem transition from vegetative to the reproductive state were down-regulated after LV, but AGAL2, AHL29, LD and FZL in one cluster were over-represented after LV treatment. Another cluster consisted of down-regulated genes in NV and LV plants, while SV resulted in a considerable up-regulation. A third cluster represents minor up-regulation at NV, while the fourth cluster contains over-represented genes in NV and SV bulbs (ARF2, MBR1, J3, FT1, FD2) that are remarkably down-regulated after long cold storage (Fig. 8a). It should be noted that RT-qPCR method has higher variance in comparison with RNAseq results and, therefore, the expression results obtained by these methods might vary. The qRT-PCR measurements of FT1 and FT4 expression were not always consistent with the in silico results, but the relative expression of CAL and AP1 confirmed the calculated trends (Fig. 8b).
KEGG analysis of the circadian rhythm genes
The endogenous circadian oscillator provides plants with information on the daily changes and thereby controls a number of developmental processes. The analysis of the specific DEGs involved in the circadian rhythm (ko04712) using KEGG algorithm, reveals that the expression of at least half of the number of genes presented in the circadian rhythm is affected by temperature (Fig. 9). Thus, despite the complete darkness in storage, morning components LATE ELONGATED HYPOCOTYL (LHY) and CIRCADIAN CLOCK ASSOCIATED1 (CCA1) showed relatively high expression at the end of NV and SV and low expression on completion of the LV treatment.
Flowering promoting genes FLAVIN-BINDING KELCH REPEAT F-box 1(FKF1) and FT were over-expressed in SV compared to NV and LV treatments. The expression of PHYTOCHROM B (PHYB) was not significantly affected by treatments, while strong expression of PHYTOCHROM INTERACTING FACTOR 3 (PIF3) was evident in LV plants.
Long vernalization activates meristem transition in garlic
Many geophyte species, including alliaceous crops, evolved in winter-cold climates and require vernalization for successful completion of their annual development and flowering [34, 39]. In onion and garlic, photoperiod is also considered the significant environmental cue for the transition from active growth to dormancy period [58, 70]. Our data support the findings from Ben Michael et al.  that long vernalization alone satisfies the needs for apical meristem transition - the first and obligatory step for flowering and bulbing in garlic. Here we studied the effect of cold treatment on the genetic interplay involved in meristem transition to active growth and flowering.
Our results show the dominant environmental cue of low temperature for triggering meristem transition (Fig. 1), while photoperiod might play only a regulatory role at the later stages of plant development . It should be noted, that in onion low temperatures promote flowering and long photoperiod induces bulbing . In garlic, however, cold satisfied the requirements for both processes . Therefore, the molecular control of meristem transition and further plant development might vary between A. cepa and A. sativum.
Epigenetics is involved in meristem shifting in garlic
Transcriptome analysis revealed significant modifications in epigenetic and chromatin-based processes in the apical meristems in association with pre-planting storage conditions (Fig. 3). It is widely accepted that methylation plays a central role in adaptation to environmental stress and in plants’ “winter memory” [28, 42]. Our study shows significant methylation activity in all comparisons in both, vernalized and non-vernalized bulbs, and these pathways should be further investigated.
Vernalization in dark triggers light reaction and photoperiod-associated genes
The meristem transcriptional changes coded for major modifications in two important metabolic routes: (1) dormancy release and initiation of active growth, and (2) meristem transition from vegetative to reproductive state, with the consequent end of apical dominance and outgrowth of axillary buds. Previous physiological and molecular studies suggested that the temperature-regulated dormancy release is associated with hormone transport, and carbon, peptide and flavonoids’ metabolism [59, 63, 66, 93]. We now provide evidence that for the same physiological change, at the transcription level, numerous genes are involved in cell division, differentiation, and growth (Figs. 2b and 3). Garlic meristems reacted to vernalization by intense re-organization of a large range of cellular components and metabolic processes (Fig. 4). In dark-stored vernalized meristems, we found a large number of GO terms associated with light reaction and photosynthesis. We thus conclude that the tagged genes annotated as involved in “photosynthetic activity” are basically regulated here by low temperature. Similar results were obtained for the evergreen orchardgrass, where vernalization and photoperiod genes were involved and co-expressed during meristem transition . Leeggangers et al.  reported that GO terms of photomorphogenesis, photoperiodism, response to light stimulus, and circadian rhythm were upregulated in dormant tulip bulbs during floral induction by high summer temperatures. They offered the explanation that the aboveground drying leaves may receive light or photoperiodic signals. These cues are then passed on to the meristem with the consequent flower induction effect. In our experiments, however, leaves and roots were all removed prior to the bulbs storage in the dark, hence light induction of photoperiodic genes was unlikely. We thus conclude that upregulation of the above genes is exclusively cold dependent.
Gene co-expressing during meristem transition
In plants, meristem transition from the vegetative to reproductive state is coded by a multifaceted genetic network involving photoperiod, gibberellins, vernalization and autonomous signaling modules [2, 21, 61, 84]. The reception of photoperiod cue provides reliable information on seasonal changes and is thus well-conserved , whereas vernalization requirements and their genetic make-up vary considerably among species. The latter variation stems from the fact that acquisition of winter memory evolved independently in several plant taxa for better adaptation to local climate cycles . Thus, considerable differences exist in vernalization perception between Arabidopsis and related species (Brassicaceae; ), cereals (Poaceae; ), Rosaceae  and others. We employed garlic transcriptome catalog  in a comparative study with a large number of genes and pathways. These analyses of flowering genes revealed a few modules in gene co-expression network (Fig. 5). While garlic transcripts annotated as vernalization genes took only minor parts in this network, they were centrally located. Vernalization genes co-expressed with numerous genes associated with light reaction and with those involved in the photoperiodic control of flowering (Fig. 5). Therefore, we propose that the molecular cascade of vernalization includes, to some extent, genes coding for the photoperiodic pathway. The latter are homologs to the Arabidopsis STO (major role in response to light and photoperiodism, ), LHY (circadian rhythm and photoperiodism, ), FKF1 (regulation of circadian clock-dependent transition to flowering time; ), VOZ1 (promotion flowering downstream of phytochrome B (phyB), ) and CKA1 and CKA2 (circadian rhythm, photomorphogenesis, ).
Treatment effects on differential expression of vernalization and photoperiod genes were evident (Fig. 8). A large number of those genes were downregulated in LV plants (Figs. 6 and 7), thus indicating that, on a molecular level, a transition to the reproductive stage had already occurred by the end of the LV treatment. Indeed, the differential expression of the meristem-identity genes supports this assumption (Fig. 8).
Flowering signal integrators, such as the Arabidopsis TFL1, play an important role as a hub connecting between the photoperiod and temperature coding genes during meristem transition . In strawberry, FvTFL1 integrates the photoperiod and temperature signals to repress flowering . In garlic, gene co-expression analysis pointed out a few genes with a putative key role in the flowering process. Hence, SUF4 links between the three modules that positively co-express with numerous flowering genes (Fig. 5), while a vernalization-associated gene AGL19 negatively correlates with photoperiodic and meristem identity genes and seems to have a considerable repressing effect on these genes (Additional file 3: Figure S3). Therefore, SUF4 and AGL19 are likely to serve as floral integrators in garlic (and other geophytes), but further research is needed to determine their actual role in floral initiation.
The onset of flower development in Arabidopsis is controlled by the interplay between LEAFY (LFY), APETALA1 (AP1) and its paralog CAULIFLOWER (CAL) . We have already shown that garlic homolog of LFY (gaLFY) acts as a key-gene in the floral transition, as well as in inflorescence and flower differentiation . Our present findings show only a weak expression of CAL and AP1 in non-vernalized plants and a clear upregulation following vernalization (Fig. 9). We thus propose that the latter genes interact with gaLFY during floral transition in garlic.
FT genes are involved in meristem transition
In addition to the important and conserved role of FT as the ‘florigen’, other genes of this family serve in many plant species as “master” regulators of an array of functions. In geophytes, they are also involved in bulbing/tuberization [1, 45] and are associated with vernalization, photoperiod, and heat response [47, 65]. In garlic, six FT homologs were identified , including AsFT2, a homolog of onion ‘florigen’, whose expression well correlates with flowering initiation in garlic [45, 76].
Flowering and bulbing in onion, bunching onion (A. fistulosum) and many other geophytes are two competitive processes on both physiological and molecular levels [11, 17, 45, 70]. Our present results show that similar to onion, in garlic low temperatures affect down-regulation of AsFT4 and up-regulation of AsFT2. Differences in FT expression between short- and long-day onion varieties were already shown [56, 57]. Since in vernalized garlic bulbing is also possible under short day, FT genes in this species can express and function in different manner, and further functional characterization of this family in various garlic genotypes will help to understand their role in florogenesis and bulbing.
Circadian clock acts as integrator between temperature and meristem transition
Plant circadian clock, an endogenous timing system based on a cellular self-sustained oscillation , can be adjusted by light cycles and temperature, and may thus serve as an integration point for the two major environment cues [13, 23]. Indeed, our data show a strong expression of genes associated with the internal clock cycle, e.g., CKA1/2, LHY, ELF3, PIF3 in vernalized garlic (Fig. 9), as previously demonstrated in chestnut and other plant species [23, 30]. The involvement of circadian rhythm in low-temperature induced florogenesis was also confirmed in lily , orchard grass (Dactylis glomerata; ), radish (Raphanus sativus; [53, 64]), Chinese cabbage (Brassica rapa ssp. pekinensis; ), as well as in bud dormancy break in trees [51, 68].
We argue that during dark cold storage of bulbs, the circadian clock integrates the low-temperature signal and consequently induces the expression of genes associated with vernalization, photoperiod and meristem transition (Fig. 10). Indeed, recent reports demonstrated that Phytochrome B (PHYB), a major light sensor of the circadian rhythm, acts as thermosensor [32, 48, 69], while EARLY FLOWERING 3 (ELF3) transmits thermal information . Similarly, in garlic, homolog of PhyB may act as thermosensor, conveying temperatures information to the circadian clock via an interface with Phytochrome Interacting Factors (PIFs) . Moreover, cold induced the differential expression of at least 14 circadian clock genes, tightly connected to the core of the photoperiodic pathway (e.g., LHY, CCA1, FKF1, CDF1, and PRRs; ) (Fig. 9). We conclude that molecular control of vernalization in garlic integrates into the reserved photoperiodic pathway in an upstream point, possibly by the same receptors.
Long dark cold exposure of bulbs is a major cue for flowering and bulbing, and its interactions with the genetic makeup of the individual plant dictate the phenotypic expression during growth stage. Low temperatures stimulate large cascades of molecular mechanisms in garlic, and a variety of flowering mechanisms, including photoperiodic pathway, operate together via circadian clock for the benefit of meristem transition, annual life cycle and viable reproduction results.
Plant materials, storage treatments, and growth conditions
Flowering hardneck garlic #87, selected in 2004 from a segregating seedling population and then clonally propagated at the Agricultural Research Organization (ARO), The Volcani Center, Israel, was used in our experiments (Additional file 1: Figure S1). Plants were grown in 2015–2016 at the Experimental Farm, the Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot, Israel (latitude 31.9°N).
In July–August 2015, freshly harvested bulbs were stored under ambient conditions (22–35 °C), in a roofed open-shed. The experiment was designed to ensure that all bulbs are planted at the same time to avoid seasonal variations.In August healthy intact uniform sized bulbs were placed in a temperature-controlled dark chamber at 4 °C for 12 weeks (LV, long vernalization), and in October for 4 weeks (SV, short vernalization). Control bulbs (NV, no vernalization) were kept in the original shed where ambient temperatures ranged between 20 and 30 °C. On November 11, 2015, when all cold treatments ended, uniform-sized cloves from healthy intact bulbs were sampled for RNA extraction and concomitantly intact cloves were planted in the soil, in a 30% shaded screenhouse, in triplicate randomized block design, 30 cloves per replication, at a 50 plants/m2 density. Indoors air temperature records and photoperiod for the entire growth period are presented in Additional file 2: Figure S2. Standard agriculture practices were applied throughout, yet for each treatment, fertigation [“Shefer” liquid fertilizer (N:P:K = 59:35:94 g L-1), Dshanim, Israel] discontinued when the aboveground parts died back, in May–June 2016, and bulbs were selectively harvested upon maturation. Twenty intact growing plants per replication (60 plants/treatment) served for phenology and morphology studies. On full bloom, floral buds and flowers were detached from the receptacle and counted in ten inflorescences per treatment. After harvest, foliage leaves and floral stalk were detached from the bulb, for fresh and dry weight measurements. The weight and structure of cured bulbs were recorded 3 weeks after harvest. Data on plant development were subjected to one-way analysis of variance (ANOVA) with the Tukey-Kramer test using JMP software .
Tissues sampling for RNA/DNA analyses and extraction procedures
When storage treatments ended, three biological replications of apical buds were sampled and immediately dipped in liquid nitrogen before storage at − 80 °C. Total RNA was extracted using the CTAB protocol , and extract quality was assessed using the Agilent 2100 Bioanalyzer. Only extract with a minimum RNA integrated number value of seven were used in the following studies.
Library preparation and sequencing were performed by the Genome Center, Life Sciences and Engineering, Technion, Israel. Six libraries of paired-end RNA-seq of 100 nucleotides, were prepared for analysis using Illumina Hiseq 2500 and Trueseq protocols. Raw reads were subjected to a cleaning procedure with FASTX Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html, version 0.0.13.2) as follows:
trimming read-end nucleotides with quality scores < 30 using fastq_quality_trimmer;
read pairs were discarded if either one had less than 70% base pairs with quality score ≤ 30 using fastq_quality_filter.
Following processing and cleaning, a total of approximately 230 million cleaned paired-end reads were assembled de novo, using Trinity software (Trinityrnaseq version v2.3.2, ) with default parameters by Trimmomatic option.
The resulting data were compiled with previously published transcriptome catalog of fertile garlic . The new transcriptome catalog now contains a total of 112,388 trinity genes with an average length of 1.178 bp and N50 of 1812 bp. The data have been deposited in the NCBI Sequence Read Archive (SRA) database as a bioproject PRJNA566287.
Abundance estimation and differential expression analysis
The cleaned reads from each library were aligned with the newly assembled transcriptome, using the Bowtie2 aligner . The abundance estimation was calculated using the Trinity protocol  and the Expectation-Maximization method (RSEM). The normalized expression was calculated by TMM (trimmed mean of M values) normalization and by FPKM (fragments per kilobase of transcript per million mapped reads).
Comparisons between each pair of samples and differential expression analyses of the sequence count data were performed by Bioconductor edgeR package in the R environment . Differential expression was defined as a value larger than the two-fold difference in transcription expression with a false discovery-corrected statistical significance below 0.05. The differentially expressed transcripts that had a minimum of one pairwise comparison between treatments were examined, using cluster analysis. Based on the means of two replications, hierarchical clustering of transcripts and samples was performed and clusters were extracted using R scripts.
The assembled transcriptome was used for a search of the NCBI non-redundant (nr) protein database, employing the DIAMOND software . Homologous sequences were identified within the Swiss-Prot database using the BLASTX tool  at E-value threshold of 10− 5. The results were exported to Blast2GO version 4.0  for gene ontology (GO) assignments. The KAAS tool (KEGG Automatic Annotation Server; http://www.genome.jp/tools/kaas/) was used for KEGG orthology and KEGG pathway assignments. Gene ontology enrichment analyses were carried out using Blast2GO  software based on Fisher’s Exact Test  with multiple testing correction of false discovery rates (FDR; ). The threshold was set as FDR with corrected P-value smaller than 0.05. The ReviGO web server was used for visualization of the GO terms in a semantic similarity-based scatterplot (http://revigo.irb.hr; ).
Real-time quantitative PCR validation
For gene expression validation, 15 genes were selected upon their putative function in flower induction and meristem transition. The RNA samples employed for the construction of libraries and transcriptome analyses were further analyzed. Two micrograms RNA samples from each of the total RNA obtained from three samples/storage treatment, served for the synthesis of cDNA using the High-Capacity Reverse Transcription Kit (Applied Biosystems, Foster City, CA). Primers were designed using primer3.0 web-tool (Additional file 4: Table S1). Amplicon identities were validated by sequencing with Sanger sequencing (Life Sciences Core Facilities, Weizmann Institute of Science, Israel).
Gene expression was validated using the Fluidigm® system. The cDNA samples along with non-template control (NTC) were analyzed by Biological Services, Weizmann Institute for Science, Israel. Each of the three technical repetitions per sample was analyzed twice for all sets of primers by using the 48.48 Dynamic Array™ IFCs, the BioMark™ HD System (Fluidigm, South San Francisco, CA) and GE Fast PCR + Melt v2 Protocol (https://www.fluidigm.com/documents). Data analysis was performed using Fluidigm Real-time PCR software. Garlic homologs of actin (AY821677) and tubulin (AY148156) served as reference genes for data normalization and calculation of relative amounts of mRNA in the studied samples. All sequences were deposited in the garlic transcriptome catalog (NCBI bioproject PRJNA566287). The expression was quantified by the ∆∆CT method  and presented in a log scale (2-ΔΔCT).
Availability of data and materials
The datasets generated during the current study are submitted to the NCBI repository, bioproject PRJNA566287.
Days after planting
Abelenda JA, Cruz-Oró E, Franco-Zorrilla JM, Prat S. Potato StCONSTANS-like1 suppresses storage organ formation by directly activating the FT-like StSP5G repressor. Curr Biol. 2016;26:872–81.
Albani MC, Coupland G. Comparative analysis of flowering in annual and perennial plants. Curr Top Dev Biol. 2010;91:323–48.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Andrés F, Coupland G. The genetic basis of flowering responses to seasonal cues. Nat Rev Gen. 2012;13:627–39.
Bandara MS, Krieger K, Slinkard AE, Tanino KK. Pre-plant chilling requirements for cloving of spring-planted garlic. Can J Plant Sci. 2000;80:379–84.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Methodol. 1995;57:289–300.
Ben Michael T, Shemesh-Mayer E, Kimhi S, Gershberg C, Forer I, de Ávila VT, Rabinowitch HD, Kamenetsky-Goldstein R. Temporal and spatial effect of low pre-planting temperatures on plant architecture and flowering in bolting garlic. Sci Hort. 2018;242:69–75.
Berry S, Dean C. Environmental perception and epigenetic memory: mechanistic insight through FLC. Plant J. 2015;83:133–48.
Bloomer RH, Dean C. Fine-tuning timing: natural variation informs the mechanistic basis of the switch to flowering in Arabidopsis thaliana. J Exp Bot. 2017;8:5439–52.
Bouché F, Woods DP, Amasino RM. Winter memory throughout the plant kingdom: different paths to flowering. Plant Physiol. 2017;73:27–35.
Brewster JL. Onions and other vegetable Alliums. 2nd ed. Wallingford: CABI; 2008.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.
Casal JJ, Qüesta JI. Light and temperature cues: multitasking receptors and transcriptional integrators. New Phytol. 2018;217:1029–34.
Chang S, Puryear J, Cairney J. A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993;1:113–6.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Dennis ES, Peacock WJ. Epigenetic regulation of flowering. Curr Opin Plant Biol. 2007;10:520–7.
Dong Y, Cheng Z, Meng H, Liu H, Wu C, Khan AR. The effect of cultivar, sowing date and transplant location in field on bolting of welsh onion (Allium fistulosum L.). BMC Plant Biol. 2013;13:154. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2229-13-154.
Duncan S, Holm S, Questa J, Irwin J, Grant A, Dean C. Seasonal shift in timing of vernalization as an adaptation to extreme winter. eLife. 2015;4:e06620. https://0-doi-org.brum.beds.ac.uk/10.7554/eLife.06620.
Ezer D, Jung JH, Lan H, et al. The evening complex coordinates environmental and endogenous signals in Arabidopsis. Nat Plants. 2017;3:17087. https://0-doi-org.brum.beds.ac.uk/10.1038/nplants.2017.87.
Feng G, Huang L, Li J, Wang J, Xu L, Pan L, Zhao X, Wang X, Huang T, Zhang X. Comprehensive transcriptome analysis reveals distinct regulatory programs during vernalization and floral bud development of orchardgrass (Dactylis glomerata L.). BMC Plant Biol. 2017. https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-017-1170-8.
Fornara F, de Montaigu A, Coupland G. SnapShot: control of flowering in Arabidopsis. Cell. 2010;141:550.
Fujiwara S, Oda A, Yoshida R, et al. Circadian clock proteins LHY and CCA1 regulate SVP protein accumulation to control flowering in Arabidopsis. Plant Cell. 2008;20:2960–71.
Gil KE, Park CM. Thermal adaptation and plasticity of the plant circadian clock. New Phytol. 2019;221:1215–29.
Goslin K, Zheng B, Serrano-Mislata A, et al. Transcription factor interplay between LEAFY and APETALA1/CAULIFLOWER during floral initiation. Plant Physiol. 2017;174:1097–109.
Grabherr MG, Haas BJ, Yassour M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Guevara-Figueroa T, López-Hernández L, Lopez MG, Hurtado MD, Vázquez-Barrios ME, Guevara-Olvera L, González RG, Rivera-Pastrana DM, Torres-Robles H, Mercado-Silva EM. Conditioning garlic “seed” cloves at low temperature modifies plant growth, sugar, fructan content, and sucrose sucrose fructosyl transferase (1-SST) expression. Sci Hort. 2015;189:150–8.
Haas BJ, Papanicolaou A, Yassour M, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.
He Y, Li Z. Epigenetic environmental memories in plants: establishment, maintenance, and reprogramming. Trends Genet. 2018;34:856–66.
Hepworth J, Dean C. Flowering Locus C’s lessons: conserved chromatin switches underpinning developmental timing and adaptation. Plant Physiol. 2015;168:1237–45.
Ibañez C, Ramos A, Acebo P, Contreras A, Casado R, Allona I, Aragoncillo C. Overall alteration of circadian clock gene expression in the chestnut cold response. PLoS One. 2008. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0003567.
JMP 13 Pro. Cary, NC: SAS Institute Inc. 2016.
Jung JH, Domijan M, Klose C, et al. Phytochromes function as thermosensors in Arabidopsis. Science. 2016;354:886–9.
Jung WY, Lee A, Moon JS, Kim YS, Cho HS. Genome-wide identification of flowering time genes associated with vernalization and the regulatory flowering networks in Chinese cabbage. Plant Biotechnol Rep. 2018;12:347–63.
Kamenetsky R. Biodiversity of geophytes: phytogeography, morphology and survival strategies. In: Kamenetsky R, Okubo H, editors. Ornamental geophytes: from basic science to sustainable production. Boca Raton: CRC Press, Taylor and Francis Group; 2012. p. 57–76.
Kamenetsky R, Zaccai M, Flaishman MA. Florogenesis. In: Kamenetsky R, Okubo H, editors. Ornamental geophytes: from basic science to sustainable production. Boca Raton: CRC Press, Taylor and Francis Group; 2012. p. 197–232.
Kamenetsky R. Garlic: botany and horticulture. In: Horticultural reviews-Westport then New York, vol. 33; 2007. p. 123–72.
Kamenetsky R, Faigenboim A, Shemesh-Mayer E, et al. Integrated transcriptome catalogue and organ-specific profiling of gene expression in fertile garlic (Allium sativum L.). BMC Genomics. 2015. https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-015-1212-2.
Kamenetsky R, London Shafir I, Zemah H, Barzilay A, Rabinowitch HD. Environmental control of garlic growth and florogenesis. J Am Soc Hortic Sci. 2004;129:144–51.
Khodorova NV, Boitel-Conti M. The role of temperature in the growth and flowering of geophytes. Plants. 2013;2:699–711.
Kim DH, Doyle MR, Sung S, Amasino RM. Vernalization: winter and the timing of flowering in plants. Annu Rev Cell Dev Bi. 2009;25:277–99.
Kurokura T, Mimida N, Battey NH, Hytönen T. The regulation of seasonal flowering in the Rosaceae. J Exp Bot. 2013;64:4131–41.
Lämke J, Bäurle I. Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants. Genome Biol. 2017. https://0-doi-org.brum.beds.ac.uk/10.1186/s13059-017-1263-6.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009. https://0-doi-org.brum.beds.ac.uk/10.1186/gb-2009-10-3-r25.
Lee H, Suh SS, Park E, Cho E, Ahn JH, Kim SG, Lee JS, Kwon YM, Lee I. The AGAMOUS-LIKE 20 MADS domain protein integrates floral inductive pathways in Arabidopsis. Gen Dev. 2000;14:2366–76.
Lee R, Baldwin S, Kenel F, McCallum J, Macknight R. FLOWERING LOCUS T genes control onion bulb formation and flowering. Nat Commun. 2013. https://0-doi-org.brum.beds.ac.uk/10.1038/ncomms3884.
Leeggangers HA, Nijveen H, Bigas JN, Hilhorst HW, Immink RG. Molecular regulation of temperature-dependent floral induction in Tulipa gesneriana. Plant Physiol. 2017;173:1904–19.
Leeggangers HA, Rosilio-Brami T, Bigas-Nadal J, Rubin N, et al. Tulipa gesneriana and Lilium longiflorum PEBP genes and their putative roles in flowering time control. Plant Cell Physiol. 2018;59:90–106.
Legris M, Klose C, Burgie ES, Rojas CC, Neme M, Hiltbrunner A, Wigge PA, Schäfer E, Vierstra RD, Casal JJ. Phytochrome B integrates light and temperature signals in Arabidopsis. Science. 2016;354:897–900.
Li W, Liu X, Lu Y. Transcriptome comparison reveals key candidate genes in response to vernalization of Oriental lily. BMC Genomics. 2016. https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-2955-0.
Li F, Sun J, Wang D, Bai S, Clarke AK, Holm M. The B-box family gene STO (BBX24) in Arabidopsis thaliana regulates flowering time in different pathways. PLoS One. 2014. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0087544.
Liu G, Li W, Zheng P, Xu T, Chen L, Liu D, Hussain S, Teng Y. Transcriptomic analysis of ‘Suli’pear (Pyrus pyrifolia white pear group) buds during the dormancy by RNA-Seq. BMC Genomics. 2012. https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2164-13-700.
Liu X, Wang Q, Gu J, Lü Y. Vernalization of Oriental hybrid lily ‘Sorbonne’: changes in physiology metabolic activity and molecular mechanism. Mol Biol Rep. 2014;41:6619–34.
Liu C, Wang S, Xu W, Liu X. Genome-wide transcriptome profiling of radish (Raphanus sativus L.) in response to vernalization. PLoS One. 2017. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0177594.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25:402–8.
Lu SX, Liu H, Knowles SM, Li J, Ma L, Tobin EM, Lin C. A role for protein kinase casein kinase2 α-subunits in the Arabidopsis circadian clock. Plant Physiol. 2011;157:1537–45.
Lyngkhoi F, Khar A, Mangal M, Gaikwad AB, Thirunavukkarasu N. Expression analysis and association of bulbing to FLOWERING LOCUS T (FT) gene in short day onion (Allium cepa L.). Indian J Genet Pl Br. 2019;79:77–81.
Manoharan RK, Suk J, Han H, Vijayakumar H, Subramani B, Thamilarasan SK, Park JI, Nou IS. Molecular and functional characterization of FLOWERING LOCUS T homologs in Allium cepa. Molecules. 2016;21:217. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules21020217.
Mathew D, Forer Y, Rabinowitch HD, Kamenetsky R. Effect of long photoperiod on the reproductive and bulbing processes in garlic (Allium sativum L.) genotypes. Environ Exp Bot. 2011;71:166–73.
Maurya JP, Bhalerao RP. Photoperiod-and temperature-mediated control of growth cessation and dormancy in trees: a molecular perspective. Ann Bot. 2017;120:351–60.
Michaels SD, Amasino RM. Loss of FLOWERING LOCUS C activity eliminates the late-flowering phenotype of FRIGIDA and autonomous pathway mutations but not responsiveness to vernalization. Plant Cell. 2001;13:935–41.
Mouradov A, Cremer F, Coupland G. Control of flowering time: interacting pathways as a basis for diversity. Plant Cell. 2002;14:111–30.
Muranaka T, Oyama T. Monitoring circadian rhythms of individual cells in plants. J Plant Res. 2018;131:15–21.
Née G, Xiang Y, Soppe WJ. The release of dormancy, a wake-up call for seeds to germinate. Curr Opin Plant Biol. 2017;35:8–14.
Nie S, Li C, Xu L, Wang Y, Huang D, Muleke EM, Sun X, Xie Y, Liu L. De novo transcriptome analysis in radish (Raphanus sativus L.) and identification of critical genes involved in bolting and flowering. BMC Genomics. 2016. https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-2633-2.
Noy Porat T, Cohen D, Mathew D, Eshel A, Kamenetsky R, Flaishman MA. Turned on by heat: differential expression of FT and LFY-like genes in Narcissus tazetta during floral transition. J Exp Bot. 2013;64:3273–84.
Okubo H, Sochacki D. Botanical and horticultural aspects of major ornamental geophytes. In: Kamenetsky R, Okubo H, editors. Ornamental geophytes: from basic science to sustainable production. Boca Raton: CRC Press, Taylor and Francis Group; 2012. p. 77–117.
Oliver SN, Finnegan EJ, Dennis ES, Peacock WJ, Trevaskis B. Vernalization-induced flowering in cereals is associated with changes in histone methylation at the VERNALIZATION1 gene. Proc Nat A Sci. 2009;106:8386–91.
Penfield S. Temperature perception and signal transduction in plants. New Phytol. 2008;179:615–28.
Qiu Y, Li M, Jean R, Moore CM, Chen M. Daytime temperature is sensed by phytochrome B in Arabidopsis through a transcriptional activator HEMERA. Nat Commun. 2019. https://0-doi-org.brum.beds.ac.uk/10.1038/s41467-018-08059-z.
Rabinowitch HD. Physiology of flowering. In: Rabinowitch HD, Brewster JL, editors. Onions and allied crops, vol. 1. Boca Raton: CRC Press; 1990. p. 113–34.
Rahim MA, Fordham R. Effect of storage temperature on the initiation and development of garlic cloves (Allium sativum L.). Sci Hort. 1988;7:25–38.
Rahim MA, Fordham R. Environmental manipulation for controling bulbing in garlic. Acta Hortic. 2001;555:181–8. https://0-doi-org.brum.beds.ac.uk/10.17660/ActaHortic.2001.555.27.
Rantanen M, Kurokura T, Jiang P, Mouhu K, Hytönen T. Strawberry homologue of TERMINAL FLOWER 1 integrates photoperiod and temperature signals to inhibit flowering. Plant J. 2015;82:163–73.
Ream TS, Woods DP, Amasino RM. The molecular basis of vernalization in different plant groups. Cold Spring Harb Symp Quant Biol. 2012;7:105–15 Cold Spring Harbor Laboratory Press.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Rohkin-Shalom S, Gillett D, Zemach H, Kimhi S, Forer I, Zutahy Y, Tam Y, Teper-Bamnolker P, Kamenetsky R, Eshel D. Storage temperature controls the timing of garlic bulb formation via shoot apical meristem termination. Planta. 2015;242:951–62.
Rotem N, David-Schwartz R, Peretz Y, Sela I, Rabinowitch HD, Flaishman M, Kamenetsky R. Flower development in garlic: the ups and downs of gaLFY expression. Planta. 2011;233:1063–72.
Sakamoto T, Kimura S. Plant temperature sensors. Sensors. 2018. https://0-doi-org.brum.beds.ac.uk/10.3390/s18124365.
Samach A, Onouchi H, Gold SE, Ditta GS, Schwarz-Sommer Z, Yanofsky MF, Coupland G. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science. 2000;288:1613–6.
Searle I, He Y, Turck F, Vincent C, Fornara F, Kröber S, Amasino RA, Coupland G. The transcription factor FLC confers a flowering response to vernalization by repressing meristem competence and systemic signaling in Arabidopsis. Genes Dev. 2006;20:898–912.
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:2498–504.
Shemesh-Mayer E, Ben-Michael T, Rotem N, Rabinowitch HD, Doron-Faigenboim A, Kosmala A, Perlikowski D, Sherman A, Kamenetsky R. Garlic (Allium sativum L.) fertility: transcriptome and proteome analyses provide insight into flower and pollen development. Front Plant Sci. 2015. https://0-doi-org.brum.beds.ac.uk/10.3389/fpls.2015.00271.
Shim JS, Kubota A, Imaizumi T. Circadian clock and photoperiodic flowering in Arabidopsis: CONSTANS is a hub for signal integration. Plant Physiol. 2017;173:5–15.
Song YH, Shim JS, Kinmonth-Schultz HA, Imaizumi T. Photoperiodic flowering: time measurement mechanisms in leaves. Annu Rev Plant Biol. 2015;66:441–64.
Strasser B, Alvarez MJ, Califano A, Cerdán PD. A complementary role for ELF3 and TFL1 in the regulation of flowering time by ambient temperature. Plant J. 2009;58:629–40.
Sung S, Amasino RM. Vernalization and epigenetics: how plants remember winter. Curr Opin Plant Biol. 2004. https://0-doi-org.brum.beds.ac.uk/10.1016/j.pbi.2003.11.010.
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0021800.
Takato I, Kay SA. Photoperiodic control of flowering: not only by coincidence. Trends Plant Sci. 2006;11:550–8.
Teotia S, Tang G. To bloom or not to bloom: role of microRNAs in plant flowering. Mol Plant. 2015;8:359–77.
Trevaskis B, Hemming MN, Dennis ES, Peacock WJ. The molecular basis of vernalization-induced flowering in cereals. Trends Plant Sci. 2007;12:352–7.
Turck F, Fornara F, Coupland G. Regulation and identity of florigen: FLOWERING LOCUS T moves center stage. Annu Rev Plant Biol. 2008;59:573–94.
Upton GJ. Fisher’s exact test. J R Stat Soc A Stat Soc. 1992;155:395–402.
Wang Y, Zhao H, Wang Y, Yu S, Zheng Y, Chan Z. Comparative physiological and metabolomic analyses reveal natural variations of tulip in response to storage temperatures. Planta. 2019;249:1379–90.
Wu C, Wang M, Cheng Z, Meng H. Response of garlic (Allium sativum L.) bolting and bulbing to temperature and photoperiod treatments. Biol Open. 2016;5:507–18.
Wu C, Wang M, Dong Y, Cheng Z, Meng H. Growth, bolting and yield of garlic (Allium sativum L.) in response to clove chilling treatment. Sci Hort. 2015;194:43–52.
Yasui Y, Kohchi T. VASCULAR PLANT ONE-ZINC FINGER1 and VOZ2 repress the FLOWERING LOCUS C clade members to control flowering time in Arabidopsis. Biosci Biotechnol Biochem. 2014;78:1850–5.
This study was partly funded by Classeed Ltd. The funding body did not participate in the design of the study, the collection, analysis, and interpretation of data, and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Flowchart depicting the experimental design, including pre-planting treatments, morpho-physiological studies, transcriptome analyses and validation of candidate genes in #87 garlic.
Day length and mean daily temperatures from October 2015 to September 2016, in Rehovot, Israel.
Negative co-expression of garlic genes associated with vernalization (green), photoperiod (blue) pathway and meristem transition (purple). Data from three vernalization treatments were analyzed using the network-drawing software Cytoscape . Pearson correlation value higher than 0.9.
List of primers designed for qRT-PCR validation. Garlic homologs of actin (AY821677), tubulin (AY148156), and NADH (TRINITY_DN98625_c0_g3) served as reference genes for data normalization and calculation of relative amounts of mRNA in the studied samples.
List of 74 genes associated with meristem transition, photoperiodic pathway and vernalization, their relevant GO terms, and short description according to UniProt database and literature survey.
About this article
Cite this article
Ben Michael, T.E., Faigenboim, A., Shemesh-Mayer, E. et al. Crosstalk in the darkness: bulb vernalization activates meristem transition via circadian rhythm and photoperiodic pathway. BMC Plant Biol 20, 77 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-020-2269-x
- Allium sativum
- Reproductive meristem
- Low temperature