From Embryo to Adult: Low Temperatures Affect Phase Transitions of Allium sativum L. from Germination to Flowering

: Juvenile and vegetative adult shoot apical meristems (SAM) are actively involved in acquisition of flowering competence, while the embryonic SAM is regarded as “responsible” only for germination. Comparative analyses of imbibed and stratified seeds of garlic Allium sativum show that only stratified seedlings produced bulbs and flower stems at the end of the season. Since the seed morphology of stratified and non ‐ stratified seeds prior to sowing was similar, the differences are attributed to the molecular alterations in the embryonic SAM. Functional annotation analysis of 3000 differentially expressed genes suggests that seed imbibition reactivates the embryonic cell cycle, initiates the metabolism, and primes garlic seed germination. Stratification enhances DNA modifications, biosynthesis, cellular transport, and tissue development. Similar to vernalization of the vegetative buds, stratification of the embryonic SAM resulted in altered expression of meristem ‐ identity and flowering homologs. Phase transitions from seed germination to flowering and bulbing in A. sativum are tightly connected, and possibly associated with downregulation of specific flowering repressor(s). The embryonic SAM plays an important role not only in seed germination, but in the entire plant life cycle, providing the foundation for the genetic regulation of major functional shifts in metabolism and development.


Introduction
The regulation of plant developmental transitions and organ growth are carefully orchestrated by internal and environmental signals. The combined effect of multiple cues affects the timing of the progress from seed germination to completion of the life cycle [1]. Post-embryonic development proceeds through the successive juvenile, adult vegetative, and reproductive phases. Although morphologically distinct, these phases are linked by the same regulatory networks. It is generally accepted that plant transition from juvenility to adulthood is accompanied by the acquisition of reproductive competence, followed by morphological makeover [2]. However, the embryonic shoot apical meristem (SAM) lacks the competence to develop reproductive organs, and such a transition is possible only in the post-juvenile phase [3].
In temperate environments, seed dormancy release and germination often require a prolonged exposure to cold, wet conditions, called stratification [4]. Similarly, in order to be competent to flower, many temperate plants must undergo exposure to a prolonged cold period termed vernalization [5]. The co-existence of these processes in the life cycle suggest that natural selection for the adaptation to cold winter should act jointly in seeds and vegetative buds. Traditionally, most of the research on stratification and vernalization mechanisms was carried out separately, but recent reports have established associations between seed germination and flowering in Arabidopsis and other plants [6][7][8]. Since adaptation to long winters evolved several times in the course of plant evolution, various regulatory mechanisms were developed in different plant groups [1]. In general, the accumulation of chilling hours regulates changes in water status, hormonal balance, respiration, and carbohydrate mobilization [9], and consequent increases in plant sensitivity to photoperiodic changes [5].
Such mechanisms are especially important in geophytes-perennial herbs that survive unfavorable periods as underground dormant storage organs-such as bulbs, corms, and rhizomes. Geophytes were metaphorically named "plant computers" because their long-term environmentdependent "memory" strongly affects life cycle and phase transitions to flowering and storage-organ formation [10,11]. During harsh weather in the course of the annual cycle, adult geophytes rest as storage organs, but when the environment becomes favorable for growth, storage organs provide the energy required for re-sprouting and initial growth [12]. The physiology of economically important geophytes (e.g., onion, tulip, lily), has been thoroughly studied [10,13], yet there is only limited understanding of the molecular regulation of their life processes and phase transitions [14][15][16][17].
Many economically important geophytes are clonally propagated (e.g., potatoes, garlic, shallots, ornamental bulbs, and corms), therefore, their seed biology got only a little attention. Most springflowering geophytes depend on 6-8 weeks at 4-10 °C germination-inducing stratification, but the seed of geophytes from warm regions do not require low temperatures for germination [18]. In many Allium and Tulipa species that have evolved in Central Asia, seeds become dormant on maturation, and their germination requires cold stratification [19][20][21]. When grown from seeds, most genera undergo an extended period of juvenility of up to seven years before they can produce flowers [22]. In post-juvenile geophytes, vernalization is required for meristem transition, and then for both bulbing and flowering [9,10,23].
Garlic Allium sativum, an important bulbous geophyte, responds to thermoperiodic cues, thus the timing and performance of its annual life cycle depends on cue from the surrounding environment. In many garlic genotypes, vernalization is required for the transition of SAM from vegetative to reproductive state, for the growth of axillary buds and for bulbing [24]. All commercial garlic cultivars do not produce seeds and are propagated only vegetatively. Recent fertility restoration of bolting genotypes opened the way for garlic breeding, with the potential for the production of garlic from true seeds [24,25]. To accomplish this task, an understanding of garlic seed biology, especially the effect of environment on germination and seedling development is of high importance.
We studied the effect of stratification on seed germination and seedling development, as well as on plant morphology and architecture during the entire annual cycle of A. sativum. Under ambient Mediterranean conditions, only seedlings raised from stratified seeds underwent transition to the reproductive stage and formed bulbs. We therefore propose that cold stratification and vernalization share common mechanisms that affect phase transitions and reproductive success during the entire life cycle of the herbaceous perennial plant.

Plant Material and Growth Conditions
True hybrid garlic seeds (Allium sativum L. 96 × 87) were harvested in July 2016 and 2018, from plants grown in a living collection at the Agricultural Research Organization (ARO), Volcani Center, Israel. Upon ripening, seeds were threshed, cleaned, and sorted by weight using a Zigzag type 1 Selecta Gravity Seed Separator (The Netherlands), and stored in paper bags in an open-sided shed.
Stratification experiments, including both morpho-physiological and transcriptome analyses, were performed in the 2016/17 and 2018/19 growing seasons. Seeds of the same batches were used in the comparative experiments and sown at the same day.

Pre-Sowing Treatments and Germination Test
The seeds of the same batch were first disinfected by immersion in 6% sodium hypochlorite for 10 min. Thereafter, one group of seeds were dipped in 0.2% Merpan fungicide (Captan) solution (Makhteshim Chemicals Ltd., Beer Sheva, Israel) solution for 10 min followed by 8 h imbibition in sterilized tap water. For stratification, disinfected seeds were mixed with moist vermiculite, and stratified in darkness at 4 °C for eight weeks. Moisture was maintained by controlled watering with 0.2% Merpan.
Treated seeds were sown in Petri dishes on MS medium mixed with 0.2% sucrose and 0.08% agar, in three replicates, 20 per replicate. Each Petri dish was wrapped in polystyrene linear lowdensity polyethylene (LLDPE) film and then wrapped in aluminum foil, for dark incubation at 20 °C. Germination percentage was recorded every four days for 50 days, and was determined by when the radicle reached a length of 4 mm.

Growth and Development In Situ
Two hundred treated seeds were sown in 40 L plastic boxes containing 80% ground coconut husk: 10% Styrofoam: 10% compost (Even Ari, Israel) in triplicates. The boxes were placed in an insect-proof 50-mesh screen-house providing 30% shade, at ARO. Regular fertigation consisted of "Shefer" liquid fertilizer (N:P:K = 59:35:94 g L −1 , Dshanim, Haifa, Israel). Temperature and photoperiod records during plant growth are presented in Supplementary Figure S1. Emerging seedlings were counted when the cotyledon was visible above ground. Records were taken once every four days, for 60 days. At the end of each growing period, all plants were harvested, and measurements taken on bulb and neck diameter in cm, in order to calculate the bulbing index (BI = bulb diameter/neck diameter) [26].

Statistical Analysis
A completely randomized one-way analysis of variance (ANOVA) with mean separation by the Tukey-Kramer test was applied using the general linear model (GLM) procedure SAS software package JMP version 12 Pro (SAS Institute Inc., Cary, NC, USA, 1989USA, -2007. Logistic regression was used to distinguish differences in germination and emergence percentages between treatments.

Tissue Sampling for RNA/DNA Analyses and Extraction Procedures
Two replications of 50 seeds each, either dry, stratified, or imbibed, were sampled for nucleic acid analyses prior to sowing. The seeds were dipped in liquid nitrogen and stored at −80 °C.
Total RNA was extracted using the CTAB protocol [27], and extract quality was assessed using an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Only extracts with a minimum RNA integrated value of seven were used for analysis. This resulted in six libraries. The pairs of replication libraries were compared and the data was averaged for further evaluation.

Transcriptome Assembly
Library preparation and sequencing were performed at the Genome Center, Life Sciences and Engineering, Technion, Israel Institute of Technology, Haifa, Israel. Six libraries of paired-end RNA sequences of 100 nucleotides were prepared for analysis using Illumina Hiseq2500 and Trueseq protocols. Raw reads were subjected to a cleaning procedure with the FASTX Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html, version 0.0.13.2) as follows: (1) Trimming read-end nucleotides with quality scores <30 using fastq_quality_trimmer; (2) Read pairs were discarded if either one had less than 70% base pairs with quality score ≤ 30 as determined using fastq_quality_filter.
Following processing and cleaning, a total of ca. 263 million cleaned paired-end reads were assembled de novo, using Trinity software (Trinityrnaseq version v2.3.2) [28] with default parameters and a Trimmomatic option [29]. The resulting data were compiled using the transcriptome catalog of fertile garlic [26]. The newly assembled catalog contained a total of 120,319 Trinity 'genes' with an average length of 909.6 bp and N50 of 1481 bp. The entire catalog has been deposited in the NCBI (National Center for Biotechnology Information) Sequence Read Archive (SRA) database as bioproject PRJNA647152.

Abundance Estimation and Differential Expression Analysis
The cleaned reads from each library were aligned with the newly assembled transcriptome, using the Bowtie2 aligner [30]. The abundance estimation was calculated using the Trinity protocol [31] 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 the Bioconductor edgeR package in the R environment [32]. Differential expression was defined as a value larger than a two-fold difference in transcription expression with a false discovery-corrected statistical significance below 0.05. The differentially expressed transcripts with 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.

Functional Annotations
The assembled transcriptome was used for a search of the NCBI non-redundant (nr) protein database, employing DIAMOND software [33]. Homologous sequences were identified within the Swiss-Prot database using the BLASTX tool [34] at E-value threshold of 10 −5 . The results were exported to Blast2GO version 4.0 [35] for gene ontology (GO) assignments.
The KAAS tool (KEGG Automatic Annotation Server http://www.genome.jp/tools/kaas) was used for KEGG ontology and KEGG pathway assignments. Gene ontology enrichment analyses were carried out using Blast2GO [35] software based on Fisher's Exact Test [36] with multiple testing corrections of false discovery rates (FDR) [37]. The threshold was set as FDR with corrected P-value smaller than 0.05. The REViGO web server, which sorts out long lists of gene ontology terms and summarizes them by removing redundant GO terms, was used for visualization of the GO terms in a semantic similarity-based scatterplot (http://revigo.irb.hr) [38].

Effect of Stratification on Seed Germination and Seedling Phenology
On the day of sowing, there were no histological and morphological differences in seed structure between imbibed and stratified seeds. The seeds had normal embryos, with no deformations or degenerated tissues ( Figure 1A). The embryonic SAM was visible ( Figure 1B), and in some seeds from both groups initial differentiation of the first leaf primordium was evident ( Figure 1C). In comparison with imbibed seeds, stratification clearly enhanced germination and emergence, in both Petri dishes and the planting mixture ( Figure 2). During juvenile development, stratification had no effect on leaf number, but it significantly affected leaf elongation. Leaf length in two-monthold seedlings raised from imbibed seeds was 7.7 ± 0.6 cm, compared with 10.9 ± 1.0 cm in stratified seedlings (data not shown). Five-month-old plants grown from imbibed seeds produced approximately 7-9 prostrate leaves, compared with 7-8 erect leaves in plants from stratified seeds ( Figure 3A). The latter formed bulbs and flower stems, while the apical meristem of the plants raised from non-stratified seeds remained vegetative and new leaf production continued ( Figure 3B). No bulbing, branching, or flowering were evident in these plants throughout the experiment.
Early in June, the plants were harvested for morphological evaluation. All plants raised from stratified seeds produced mature dormant bulbs up to 4-6 cm in diameter with closed necks, thus ranked as BI = 3 or 4+ (Figure 4). No bulbs were formed by any plant raised from imbibed seeds.

Global Transcriptome Analysis of Garlic Seeds
The high-throughput parallel RNA-Seq libraries extracted from dry, imbibed, and stratified garlic seeds, yielded 19.8 million 100-bp pair-end reads. Following quality trimming and filtration 16.7 million clean reads were generated and then assembled with the published garlic transcriptome catalog [39] (NCBI bioproject PRJNA243415) using Trinity software. The combined data produced a new comprehensive catalog of about 120,000 Trinity 'genes' with a mean length of 909.6 bp and N50 of 1481 bp.
The global gene expression profiles show sharp increases in the number of genes expressed as a result of both imbibition and stratification, compared to dry seeds ( Figure 5A). The total number of expressed genes was relatively low in dry seeds, about 43,000 transcripts, out of which approximately 40% is shared with the imbibed and stratified seeds. Seed imbibition and stratification led to a doubling of the total number of transcripts to about 84,000 and 72,000, respectively, but only half of the expressed genes were common to both imbibed and stratified seeds ( Figure 5A). In the comparison between imbibed and stratified seeds, 2933 differentially expressed genes (DEGs) were found ( Figure 5B), of which more than half, 1511, were upregulated in imbibed seeds compared to stratified seeds ( Figure 5C).
Functional analysis of the 2933 DEGs revealed enrichment of 24 and 160 biological processes in imbibed and in stratified seeds, respectively ( Figure 6). The main enriched terms in imbibed seeds were related to lipid storage, cellular biosynthesis, and to regulation of DNA endoreduplication ( Figure 6A). In stratified seeds, the biological processes with significantly enriched terms included DNA replication, cell cycle, glucose metabolism, meristem development, epigentic chromatin silencing, photosynthesis carbon fixation, and more ( Figure 6B). These terms are all associated with enhanced plant metabolism and growth.

Differential Expression of the Genes Associated with Seed Germination
Seed stratification significantly enhanced germination, emergence, and seedling vigor ( Figure  2). This physiological activity is associated, on the transcriptome level, with intensive gene expression ( Figure 6B). Hence, differential expression of genes associated with breaking dormancy, embryo development, sugar transport, and cell proliferation processes was expected. Data mining for functional annotations narrowed the list of 2933 DEGs between imbibed and stratified seeds to 38 DEGs (Supplementary Table S1). The patterns of expression of the individual gene functions were compared between dry seeds, imbibed seeds, and stratified seeds (Figure 7). The homologs of genes associated with seed dormancy and with germination are grouped in three clusters ( Figure 7A). With a single exception, the lowest expression of all DEGs was recorded in dry seeds. Genes in cluster I ( Figure 7A, violet) have relatively low expression in both dry and imbibed seeds and are over-expressed in stratified seeds. This cluster includes homologs of TRYPTOPHAN AMINOTRANSFERASE-RELATED PROTEIN 2 (TAR2), which is involved in the auxin signaling pathway that regulates embryo and seedling growth, PHOTOTROPIN 1A (PHOT1A), a blue light photoreceptor responsible for phototropism during germination and PHYTOCROME A (PHYA), which enhances germination even under low light. HEAT SHOCK PROTEIN 90.2 (HSP90.2) is also presented in this cluster.
Cluster II ( Figure 7A, orange) has only the single gene RGA-like1 (RGL1), which is involved in repression of gibberellin-induced germination. This gene exhibits a unique pattern in that it is highly expressed in dry seeds, poorly expressed in imbibed seeds, and barely expressed in stratified seeds, corresponding inversely to the germinability of the respective seed lots.
The largest cluster III ( Figure 7A, green) consists of genes upregulated only in imbibed seeds. Many members of the family of LATE EMBYOGENESIS ABUNDANT PROTEINS (LEA), associated with seed development and stress tolerance, share this cluster. Notably, DELAY OF GERMINATION1 (DOG1), known to regulate dormancy and juvenility, and GLUTELIN PRECURSOR 4 (GLUP4) involved in storage protein trafficking in seeds, are also overrepresented only in imbibed seeds.
Genes associated with cell proliferation and division, are clustered into two groups. Cluster I ( Figure 7B Figure 7C). Enrichment in carbon metabolism and glycolysis/gluconeogenesis in stratified seeds was confirmed by KEGG analysis (Supplementary Figure S2).

Differential Expression of Meristem Identity Genes in Germinating Seeds
Pre-sowing seed treatments greatly affected plant phenotype throughout the growth period. From seeding to maturation, plants from both treatments experienced the same environment, hence the differences in meristem transition are attributed to the pre-sowing embryonic meristematic response to seed treatment. Data mining for meristem identity genes, including those associated with meristem transition, vernalization, photoperiod, circadian clock and flowering, resulted in 26 DEGs between imbibed and stratified seeds (Supplementary Table S1). The two pre-sowing treatments had a significant effect on their activity (Figure 8). Cluster II (Figure 8, orange) is over-represented only in imbibed seeds. Two of these genes act as repressors in the photoperiod pathway of Arabidopsis, e.g., CONSTANS-LIKE 9-like (COL9) and JUMONJI 14 (JMJ14). AGAMOUS-LIKE14 (AGL14), a key player in at the complex gene regulatory network underlying SAM transitions promoting flowering in response to multiple signals, is also included in this module.
Several major components of the circadian clock are differentially expressed in response to the pre-sowing treatments. PHYTOCHROME A (PHYA) and CRYPTOCHROME2 (CRY), as well as oscillator output CYCLIC DOF FACTOR1 (CDF1) and the morning component of the circadian clock LATE ELONGATED HYPOCOTYL (LHY) were upregulated in stratified seeds. On the other hand, CIRCADIAN CLOCK ASSOCIATED1 (CCA1), usually co-expressed with LHY, showed the opposite pattern-it was highly expressed in imbibed seeds.

Discussion
Following germination, seedlings undergo several developmental phases, and the apical meristem gradually acquires competence for flowering. It is generally accepted that this process begins in the juvenile stage and is regulated by the miR156-SPL-miR172 cascade [2,3]. The plant's life cycle, however, is not a simple sequence of separate developmental events, but a chain of interactive linked processes, regulated by short-and long-term mechanisms. In this context, the role of the embryonic meristem in phase transitions and further development of the adult plant is still not clear.
Our results show that germination, flowering, and bulbing of A. sativum are tightly connected and are affected by similar environmental cues. Seed stratification promoted not only germination and emergence (Figures 1 and 2), but also had a long-term effect on meristem transition in adult plants (Figure 9).

Seed Germination is Enhanced by Low Temperatures
Similar to other Allium species [40], when seed moisture content decreases upon maturation, garlic seeds become dormant [20]. Gene expression in dry seeds is rather low, but one specific gene, a homolog of RGL1, was highly expressed. RGL1 belongs to the DELLA group (GAI, RGA, RGL1, and RGL2), and participates in coding for large multiprotein complexes that repress transcription of GAinducible genes. In Arabidopsis it enhances RGL2 repression of germination and is downregulated by imbibition [41]. We have now provided evidence of a significant decrease in RGL1 expression by cold stratification, in garlic (Figure 7). We therefore propose RGL1 as a possible marker for seed dormancy and embryo meristem activation, but its role in the genetic control of seed germination has to be further investigated.
Imbibition of the mature seed reactivates a variety of metabolic systems, such as mitochondrial biogenesis and the synthesis of new components, required for renewed cell expansion and cell division [42]. In onion, imbibition prompted the swelling of storage tissues, elongation of embryo cells, and initiation of sugar metabolism [43]. Our histological observations showed structural similarities between imbibed and stratified garlic seeds, including their embryonic meristems producing a single leaf primordium prior to germination (Figure 1). However, this morphological similarity did not correlate with germination and emergence. As in many species from temperate regions [44], cold treatment induced uniform and rapid germination (Figure 2). Water uptake doubled the number of total expressed genes, compared with dry seeds ( Figure 5). Similar to other plant species [45], imbibition facilitates the initiation of carbohydrate and lipid metabolism. Dry seeds harbor most of their cells' activities at the G1-phase, and imbibition reactivates meristematic cells into the cell cycle in preparation for seedling establishment [46]. We found that DNA endoreduplication and DNA metabolism processes are enriched in imbibed garlic seeds, and that active metabolism restarts prior to germination sensu stricto ( Figure 6A).
The expression of DOG1, a homolog of the major dormancy gene, especially drew our attention due to its dual role in seed germination [47] and in flowering control [48]. It is also involved in juvenility regulation governed by the miR156-SPL-miR172 cascade and the age-dependent flowering pathway [49]. In garlic, however, seed imbibition led to DOG1 over-expression, while after stratification its expression decreased. This DOG1 expression may precede the initiation of epigenetic silencing [50], or act to delay the germination of imbibed seeds, enabling pre-germination processes to complete.
Following an eight-week cold period, metabolic and developmental activities increased considerably ( Figure 6B). These were associated with cell cycle, sugar metabolism (e.g., FBA), cell division, and proliferation (e.g., CDC2, MCM2), as expected when storage substances gradually break down to support metabolic processes during embryo development and growth [51,52].
We conclude that initially seed imbibition reactivates the embryonic cell cycle, initiates metabolism (KINB1 and PV42A) and stress tolerance (as evidence by the large group of LEA genes, Figure 7A), and primes garlic seed in support of germination. Further, wet cold stratification serves two functions: it provides an ample supply of water to the dry seeds and the chilling initiates a number of biological pathways. These include enhanced carbon metabolism, desired DNA modifications, biosynthesis, cellular transport, and tissue development (Figures 6-8). On a physiological level, this dynamic is translated into rapid and synchronized radicle protrusion, followed by seedling emergence.

Embryonic Meristem Stratification Affects Phase Transitions in Adult Plant
Information on the long-term effect of seed treatments on the entire annual cycle and especially on phase transition is scarce. Most research is focused on the acquisition of flowering competence by juvenile and vegetative adult meristems [2,3,53], while embryonic SAM is regarded as "responsible" only for germination. We found that a long period of moist chilling of garlic seeds contributes to meristem transition from vegetative to reproductive state, even without vernalization of the adult plant. Phenotypic differences between seedlings produced from imbibed and stratified seeds were obvious after three months of growth (Figure 3), and only plants raised from stratified seeds formed flower stems and mature bulbs six to seven months after sowing (Figure 4). Bulb formation was independent of plant size (Figure 4), and the dramatic differences between plants raised from imbibed vs. stratified seed should thus be attributed to meristematic changes on the molecular level, initiated by the stratification of the embryonic SAM.
In adult plants, meristem transition from the vegetative to the reproductive state is coded by a multifaceted gene network involving photoperiod, gibberellins, vernalization, and autonomous signaling modules [54,55]. In the vernalization pathway, core transcription factors such as FLC in Brassicaceae, and VRN1 and VRN2 in cereals act as main repressors of flower transition [56][57][58]. However, since climate dependence evolved several times in the history of plants, the occurrence of different regulatory mechanisms in various plant groups is not surprising [1,5]. Thus, in Beta vulgaris the FLC homolog does not play a major role in the control of vernalization response [59]. The discovery of alternative mechanisms of flowering in wheat, maize, and rice, suggest that monocots may also be FLC-independent for flower initiation [58]. Bulbous monocots probably evolved different FLC-dependent pathways, which modulate the activity of other genes that might act as coldregulated flowering repressors [60]. It is possible, therefore, that the effect of cold on garlic embryonic and adult meristems is associated with downregulation of a specific repressor(s) of meristem transition to flowering. We found that two homologs of the known flowering repressors COL9, a repressor of CO and FT in Arabidopsis [61], and JMJ14, whichrepresses FT and TSF [62], were expressed in imbibed seeds, but downregulated by cold stratification. The same patterns were found in a few other meristem-identity genes (e.g., LD, FIP, and AGL14). Perhaps one or more of these genes are putative FLC ortholog(s) that act as repressors of meristem phase transition in garlic (Figure 9).
Further research, employing more sensitive tools (e.g., transformation systems or genome sequencing), is required to find specific flower repressors in geophyte species. Figure 9. Effect of cold period on seeds (stratification) or bulbs (vernalization) of Allium sativum. Putative flower repressor(s) (FLC orthologs) are downregulated by cold in both embryonic and adult meristems. In seedling development, acquisition of flowering competence depends on the miR156-miR172 cascade [2] and meristem enlargement, while in bud development the meristem is ready for phase transition to flowering [63].
In adult garlic plants, the molecular cascade of vernalization includes, to some extent, genes coding for the photoperiodic pathway [63]. The latter are homologs of Arabidopsis STO (which has a 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). In lily, homologs of flowering genes FT, VRN1, VRN2, SOC1, and genes of circadian rhythm were differentially expressed in vernalized bulbs, in comparison with non-vernalized [64], while during floral induction in tulip bulbs 57 homologs of Arabidopsis flowering genes were generated, including those associated with photoperiodism [65]. Similar to the results found in adult geophyte plants, our results for stratified garlic seeds show abundant processes related to meristem development, photosynthesis, light reaction, and over-representation of genes involved in photoperiodic flowering response (e.g., CRY2, PRR9, FBH4), as well as downregulation of photoperiodic repressors (e.g., JMJ14, COL9) (Figures 6 and 9).
We conclude that in garlic the processes associated with chilling-dependent meristem transition, such as increased sugar metabolism, molecular SAM activation, and an upsurge in photoperiodassociated genes, are similar for both seeds and buds/bulbs [63]. Similar phenomena were reported for peach [6], almonds [66], Persian walnut [67], and other fruit trees. Although initial genetic work on model plants identified independent molecular mechanisms for stratification [68,69] and vernalization [70,71], recent studies suggest that overlapping molecular pathways may indeed coordinate seed dormancy and vernalization-induced flowering [3,49,[72][73][74].
Embryonic and adult SAMs differ in function, size, and genetic regulation. In garlic clove, prior to vernalization, the adult SAM has already produced several leaf primordia, and is morphologically ready for the flower transition [63]. In contrast, the embryonic SAM is very small (Figure 1) and has to reach maturity (formation of 7-8 leaves and more [75]) prior to achieving flowering competence ( Figure 9). It is obvious that a juvenile garlic plant has to reach a critical mass and physiological age prior to acquisition of flower competence, and the pivotal role of miR156-miR172 interplay in this process has to be further investigated.
We propose that the embryonic SAM functions not only in seed germination and emergence, but plays an important role in the entire plant life cycle, providing the foundation for the molecular regulation of major functional transitions. The similarities we have highlighted between seed and bud chilling requirements and responses shed light on the evolution of cold-mediated plant development. Future research should address shared genetic pathways involved in the response to cold at the two separate yet critical life stages. The study of molecular mechanisms that regulate the embryonic meristem, including the epigenetic memory of the cold winter experienced by seeds, will help to elucidate successful flowering and bulb formation in adult geophytes.
Supplementary Materials: The following are available online at www.mdpi.com/2073-4395/10/11/1651/s1, Figure S1: Daylength and average air temperatures (based on 8 measures/day) in Bet Dagan, Israel during growth season 2015-2016; Figure S2: Pathways of glycolysis and gluconeogenesis in imbibed and stratified garlic seeds as revealed by KEGG algorithm osa00010. Over-representation in stratified seeds is marked in red, overrepresentation in imbibed seeds is marked in blue; Table S1: List of 61 homologous genes associated with seed dormancy and germination, proliferation, glycolysis and glucose transport and meristem identity in seeds of Allium sativum. Short description according to UniProt database and literature survey. Funding: This research was partly supported by Classeed LTD.