
==== Front
Theor Appl Genet
Theor Appl Genet
TAG. Theoretical and Applied Genetics. Theoretische Und Angewandte Genetik
0040-5752
1432-2242
Springer Berlin Heidelberg Berlin/Heidelberg

38761248
4643
10.1007/s00122-024-04643-8
Original Article
Low frequency of the wild-type freezing-tolerance LsCBF7 allele among lettuce population suggests a negative selection during domestication and breeding
http://orcid.org/0000-0002-7398-9476
Park Sunchung Sunchung.park@usda.gov

1
http://orcid.org/0000-0002-1066-7920
Shi Ainong 2
http://orcid.org/0000-0002-6135-9944
Mou Beiquan 3
1 https://ror.org/02d2m2044 grid.463419.d 0000 0001 0946 3608 U.S. Department of Agriculture, Agricultural Research Service, Beltsville, MD 20705 USA
2 grid.411017.2 0000 0001 2151 0999 Horticulture Dept, University of Arkansas, Fayetteville, AR 72701 USA
3 https://ror.org/02d2m2044 grid.463419.d 0000 0001 0946 3608 U.S. Department of Agriculture, Agricultural Research Service, Salinas, CA 93905 USA
Communicated by Diane E. Mather.

18 5 2024
18 5 2024
2024
137 6 13523 1 2024
3 5 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Key message

Sustainable winter production in lettuce requires freezing tolerant varieties. This study identified a wild-type allele of LsCBF7 that could contribute to freezing tolerance improvement in lettuce.

Abstract

Lettuce is one of the most consumed vegetables globally. While ideally grown in 13–21 °C, its cultivation extends into winter in milder climates. However, occasional freezing temperatures can significantly reduce yields. Therefore, the development of freezing-tolerant lettuce varieties has become a long-term goal of lettuce breeding programs. Despite its significance, our understanding of freezing tolerance in lettuce remains limited. Plants have evolved a coping mechanism against freezing, known as cold acclimation, whereby they can increase freezing tolerance when pre-exposed to low nonfreezing temperatures. The CBF pathway is well-known for its central role in cold acclimation. Previously, we identified 14 CBF genes in lettuce and discovered that one of them, LsCBF7, had a loss-of-function mutation. In this study, we uncovered that accessions from colder regions carried the wild-type allele of LsCBF7 and this allele likely contributed to increased freezing tolerance, with 14% of the lettuce population carrying this allele. Interestingly, in wild lettuce (L. serriola) that is considered a progenitor of cultivated lettuce, this wild-type allele was much more common, with a frequency of 90%. This finding suggests that this wild-type allele may have undergone negative selection during the domestication or breeding of lettuce. Our data strongly indicate that this allele could be linked to early bolting, an undesirable trait in lettuce, which may have driven the negative selection. While this wild-type allele shows promise for improving freezing tolerance in lettuce, it is crucial to decouple it from the early bolting trait to fully harness its potential in lettuce breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00122-024-04643-8.

http://dx.doi.org/10.13039/100005826 Agricultural Marketing Service 16SCCMAR0001 Mou Beiquan http://dx.doi.org/10.13039/100005825 National Institute of Food and Agriculture 2016-68004-24931 2015-51181-24283 2017-51181-26830 2021-51181-35903 Mou Beiquan issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Freezing temperature is a major environmental factor that restricts the geographical distribution of plants, leading to growth inhibition and yield loss in crops (Sakai and Larcher 1987; Yadav et al. 2020). To counter this freezing stress, plants have evolved diverse coping mechanisms. One such mechanism involves the enhancement of freezing tolerance when plants are exposed to nonfreezing low temperatures, a process known as cold acclimation (Thomashow 1999). While the precise mechanisms that underlie the enhancement of freezing tolerance through cold acclimation are not fully understood, many studies have suggested that C-repeat binding factors (CBFs) play a critical role in this process (Thomashow 2010; Park et al. 2015, 2018). In Arabidopsis thaliana (hereafter referred to as Arabidopsis), where the CBF pathway has been extensively studied, this pathway involves the action of three CBF genes encoding closely related AP2/ERF family transcription factors (Park et al. 2023). These three CBF genes are rapidly induced within minutes in response to low temperatures (Stockinger et al. 1997). They subsequently activate over 100 downstream target genes associated with metabolic, biochemical, and physiological processes that collectively enhance freezing tolerance (Gilmour et al. 2004; Cook et al. 2004; Kaplan et al. 2007; Li et al. 2020).

One important process during cold acclimation is the accumulation of raffinose family oligosaccharides (Taji et al. 2002; Pennycooke et al. 2003). This raffinose family plays a critical role in freezing tolerance by providing osmotic protection to plants, safeguarding chloroplasts, and stabilizing photosystem II against freezing damage (Knaupp et al. 2011). Galactinol synthase (GOLS), the enzyme catalyzing the initial and rate-limiting step in raffinose biosynthesis, plays a regulatory role in partitioning carbon between sucrose and raffinose family oligosaccharides (Saravitz et al. 1987). This gene is well-documented for its crucial role in freezing tolerance in various plant species including alfalfa (Cunningham et al. 2003), cucumber (Dai et al. 2022), Ammopiptanthus (Liu et al. 2020), Medicago falcata (Zhuo et al. 2013), Petunia (Zhang et al. 2022), and Arabidopsis (Taji et al. 2002).

In the Arabidopsis genome, there are at least seven GOLS-like genes, but only three of them, AtGOLS1, AtGOLS2, and AtGOLS3, respond to abiotic stresses (Taji et al. 2002). Among these, AtGOLS1 and AtGOLS2 were induced in response to water deficit, salinity, and heat stress by heat shock transcription factors (Panikulangara et al. 2004; Schramm et al. 2008; Nishizawa-Yokoi et al. 2009), while AtGOL3 transcripts highly accumulated exclusively in response to cold stress. This cold-induction appears to be specifically regulated by the action of CBF transcription factors (Taji et al. 2002; Park et al. 2015). CBF-overexpressing plants induced only AtGOLS3 among the seven family genes at warm temperatures, and the induction of AtGOLS3 by low temperatures was significantly impaired in CBF loss-of-function mutants (Zhao et al. 2016; Park et al. 2018). Moreover, of the 14 first-wave transcription factors (CBF1-3, MYB73, CRF2, RAV1, CRF3, ERF5, DEAR1, MYB44, CZF1, ZAT10, ZF and HSFC1) that are rapidly induced during cold acclimation, none was capable of inducing AtGOLS3 when overexpressed (Park et al. 2015). These findings indicate that the cold-induction of AtGOLS3 relies solely on the CBF regulatory module in Arabidopsis, with AtGOLS3 as a primary target of the CBF transcription factors.

The conservation of the CBF pathway is widespread among higher plants including lettuce, with many studies identifying CBF orthologous genes and highlighting their crucial roles in freezing tolerance (Owens et al. 2002; Skinner et al. 2005; Badawi et al. 2007; Carvallo et al. 2011; Hadi et al. 2011). The CBF genes have been identified as a quantitative trait locus (QTL) associated with freezing tolerance in crop species, making it a promising candidate for improving freezing tolerance in crops (Francia et al. 2007; Dumont et al. 2009; Li et al. 2013; Tayeh et al. 2013; Sieber et al. 2016; Adhikari et al. 2021). In fact, in natural populations, varying degrees of freezing tolerance were attributed to genetic variations in the CBF pathway. For example, Arabidopsis accessions originating from cooler northern latitudes generally displayed greater freezing tolerance compared to those from warmer southern latitudes (McKhann et al. 2008; Zhen and Ungerer 2008; Ågren and Schemske 2012). This clinal variation was mapped to the CBF locus in Arabidopsis, underscoring the pivotal role of CBF in adaptation to local climates (Ågren et al. 2013). Interestingly, the effect of CBF in local adaptation is reciprocal. Gain-of-function CBF variants are favored in colder regions, while they confer a disadvantage in warmer regions (Lee et al. 2024). Likewise, loss-of-function CBF alleles provide better fitness in warmer regions (Ågren and Schemske 2012; Ågren et al. 2013). Consistently, multiple studies in Arabidopsis found that accessions adapted to warmer regions carried nonfunctional CBF alleles, potentially contributing to adaptive evolution in the population (Kang et al. 2013; Gehan et al. 2015; Monroe et al. 2016). Although the precise mechanism underlying the diminished fitness associated with increased CBF activity in warmer climates remains elusive, it is postulated that erroneous activation of CBF, potentially triggered by fluctuating temperatures with a rare risk of freezing, could divert valuable resources away from growth-related processes in these warmer environments (Zhen et al. 2011; Oakley et al. 2014).

Lettuce (Lactuca sativa) is a leafy vegetable that offers numerous health benefits, including dietary fiber, minerals, and vitamins (Mou 2008). For year-round production, lettuce continues to be cultivated through the colder months in regions with mild winter temperatures, such as Imperial Valley, California, and Yuma, Arizona, which are major winter production areas in the USA (Geisseler and Horwath 2014). However, brief exposure to frost can cause lettuce leaves to blister and peel, making them vulnerable to decay and plant pathogens. This can significantly reduce lettuce yield and quality (Luna 2013; Boling 2022). Therefore, enhancing freezing tolerance has become a crucial, long-term objective in lettuce breeding for winter production.

Despite its importance, our understanding of freezing tolerance in lettuce is limited. A recent study by Park et al (2020) unveiled 14 CBF family genes in lettuce through a genome-wide search based on the reference genome constructed from the ‘Salinas’ cultivar (Reyes-Chin-Wo et al. 2017). Intriguingly, among these genes, one lettuce CBF gene, LsCBF7, was identified to have a loss-of-function mutation. This finding suggested that the mutant allele could potentially underlie freezing tolerance variation among the lettuce population. In this study, we investigated the allelic variation of the LsCBF7 gene and its contribution to freezing tolerance in lettuce and explored the frequency distribution of different LsCBF7 gene alleles among cultivated lettuce and wild lettuce populations. Furthermore, we addressed the potential of the wild-type allele as a promising candidate for enhancing freezing tolerance in lettuce. Our study sheds light on the intricate interplay between genetic variation and genetic trade-off, offering valuable insights for future breeding and improvement efforts in lettuce.

Materials and methods

Plant material and growth conditions

The lettuce genotypes used in this study consisted of 578 accessions maintained at the USDA-ARS, in Salinas, CA, USA, including crisphead, butterhead, leaf, romaine, primitive, and wild lettuce (Table S3).The lettuce plants were grown in pots with potting mix soil in a controlled growth chamber, where conditions were maintained at 20 °C with a 16 h photoperiod, except for the bolting time measurements. For bolting time measurements, plants were grown in a growth chamber with a day temperature of 33 °C and a night temperature of 19 °C. This temperature regime closely resembled the conditions used in the bolting experiment conducted by Park et al (2021). The plants were periodically rotated using a complete randomized block design. The light intensity was maintained in the range of 350 to 400 μmol m−2 s−1. For cold acclimation, plants were exposed to 4 °C for 24 h or 7 days with a light intensity of 100 μmol m−2 s−1.

Identification of GOLS genes in lettuce through phylogenetic analysis

To identify GOLS genes and determine their orthology, a comprehensive approach involving protein databases from 15 plant species including lettuce was undertaken. Among these, nine species were from the Asterid and six were from the Rosid clade (Table S1). The protein databases except for lettuce were obtained from the NCBI database. The lettuce protein database (genome v8: id37106) was sourced from https://genomevolution.org/coge. Initially, the protein sequences of Arabidopsis GOLS family genes (i.e., AtGOLS1-7) were used as queries and the protein databases were searched using the BLASTP method with an e-value threshold of < 1e−20. This search resulted in the selection of 100 genes, which were subsequently subjected to phylogenetic analysis (Table S2).

To determine orthology among the selected genes, their protein sequences were aligned using ClustalW2 with default parameters, and further refinement of the alignment was performed, if necessary, using BioEdit (Hall 1999). A phylogenetic tree was generated based on the alignment using the neighbor-joining (NJ) method implemented in MEGA 11 (Tamura et al. 2021) with the parameters of Jones-Taylor-Thornton model, uniform rates among sites, and pairwise deletion of gaps. The resulting tree was visualized using FigTree (version 1.4.4) (http://tree.bio.ed.ac.uk/software/figtree).

Freezing tolerance assay

To assess the freezing tolerance of the plants, electrolyte leakage assays were performed as described by Park et al (2015). The experiments were conducted using 28 day old plants. For cold acclimation, plants were subjected to 4 °C for 7 days under a 16 h photoperiod. Leaf disks were obtained by using a 1.5 mm diameter punch from the 4th or 5th leaf of two different plants and were placed in glass tubes. Each experimental condition was replicated three times. The tubes were then placed in a freezing chamber (Percival model LT-41VL, https://www.percival-scientific.com). To initiate the freezing process, the tubes were lowered to a temperature of − 2 °C in the dark. After an hour of equilibrium at − 2 °C, the leaf discs were frozen by adding ice-chips into the tubes. Following an additional 30 min, the tubes for the − 2 °C test were removed from the freezing chamber and placed on ice, and the temperature in the chamber was gradually lowered by 1 °C every 30 min until reaching − 12 °C. Once the specified temperatures were achieved, the tubes for each test temperature were taken out and thawed overnight on ice at 4 °C in the dark. To measure electrolytes leaked from leaf disks, 10 ml of deionized water was added into the tubes, which were then gently agitated for 4 h. After measuring the conductivity of the solutions, the solutions were transferred to new tubes, and the leaf discs were frozen at − 70 °C overnight. The next day, the preserved solutions were poured back to the frozen leaf disks and shaken for 4 h. The conductivity of the solutions was measured once again, and the ratio between the initial and subsequent conductivities was calculated.

LsCBF7 allele genotyping

SNP genotyping for the LsCBF7 alleles was performed for 578 accessions (Table S3), with 441 of these accessions previously genotyped across the genome using the genotyping-by-sequencing (GBS) method (Park et al. 2021). The SNP genotype was conducted using the rhAmp SNP Genotyping System (Integrated DNA Technologies, https://www.idtdna.com). Primers were designed using the rhAmp-genotyping design tool (Table S4). Each genotyping assay consisted of 5 µL of rhAmp Genotyping Master Mix, 0.25 µL of rhAmp Reporter Mix, 0.5 µL of rhAmp SNP Assay, 2.25 µL of water, and 2 µL of template genomic DNA (10 ng). PCR reactions were conducted on a LightCycler 480 (Roche) following a cycling condition: enzyme activation for 10 min at 95 °C, followed by 40 cycles of amplification (95 °C for 10 s, 60 °C for 30 s, 68 °C for 20 s). The fluorescent signals were analyzed using Light Cycler 480 software according to the manufacturer’s instructions. The LsCBF7 allele regions from at least 5 different accessions for each allele were sequenced to confirm the SNP.

Quantitative real-time PCR

Total RNA was extracted from leaf tissues of 18-day-old lettuce seedlings using the RNeasy Plant Mini kits (Qiagen, http://www.qiagen.com/). cDNA synthesis was carried out with 200 ng of total RNA and random primers using the Reverse Transcription System (Promega, https://www.promega.com). Quantitative real-time PCR (qRT-PCR) was conducted using the fast SYBR Green master mix (Life Technologies, http://www.lifetechnologies.com). Two housekeeping genes, eukaryotic translation initiation factor 2A (EIF2a; Ls6g95581) and isopentenyl diphosphate isomerase 2 (IPP2; Ls2g17540), were used as reference genes following the method described by Park et al (2020). The primer sequences used for qRT‐PCR are listed in Table S4. Relative expression values were calculated using the ddCt method, where the average of the two reference genes served as a baseline, and normalized values were expressed as fold-changes relative to the control treatment.

RNA-seq analysis

Leaf tissues were collected from cultivar ‘Salinas’ and ‘PI 284702’ plants exposed to 4 °C for 0 h, 24 h and 7 days. Total RNA was isolated using Qiagen Plant RNeasy kit (https://www.qiagen.com) and submitted to Novogene Corporation (https://en.novogene.com/) for RNA-seq. Sequencing was performed in a 150 bp paired-end format using the Illumina HiSeq platform (http://www.illumina.com). The RNA-seq reads were aligned to the Lactuca sativa reference genome (version 8) using STAR v2.5.2 (Dobin et al. 2013). The resulting alignments (BAM files) were processed to count gene-level reads using featureCounts (Liao et al. 2014). Differential expression analysis was conducted using the edgeR package in the R environment (Robinson et al. 2010). To mitigate the influence of lowly expressed genes, only those with a minimum of 0.5 reads per million (CPM) in at least two biological replicates were included in the analysis.

Our primary interest was in the differential gene expression resulting from the interaction between LsCBF7 alleles and the cold-acclimation treatment. To investigate the interaction between genotype and cold-acclimation treatment on gene expression, we employed a generalized linear model in edgeR. Differences were expressed as log2-fold changes (FC), where the subscript ‘cold’ represents either 24 h or 7 days of cold treatment, as follows:FCcold.PI284702-FCcold.Salinas=log2PI284702coldPI284702warm-log2SalinascoldSalinaswarm

Genes with a two-fold change (log2 = 1) or greater, and a false discovery rate of 0.01, were considered differentially cold-responsive between Salinas and PI 284702. To identify genes analysespotentially regulated by LsCBF7 in response to low temperatures, we further refined this gene set to include only those induced in PI 284702 during cold acclimation. Hierarchical clustering analyzes were performed using the hcluster method from the R package amap (R core team 2023). The resulting clusters were visualized with Treeview (http://rana.lbl.gov/EisenSoftware.htm). The RNA-seq data for PI 284702 have been deposited in the Gene Expression Omnibus under accession number GSE241604, and the RNA-seq data for Salinas was obtained from GSE134012.

Bolting time measurement

To assess bolting time, a randomized complete-block design was employed, with at least five replicates for each individual accession. Bolting time was quantified as the number of days from the initial sowing to the point when the main stem reached a length of 5 cm.

Genetic diversity between LsCBF7 allelic groups and L. serriola accessions

The GBS-based SNP data for the 441 lettuce accessions was obtained from the European Variation Archive (Project Number: PRJEB40369) (Park et al. 2021). To investigate the genetic divergence between L. sativa accessions with the LsCBF7 wild-type allele or the mutant allele, and L. serriola accessions, pairwise fixation index (FST) values between the groups were calculated using the method of Weir and Hill (2002) implemented in the R package BEDASSLE (Bradburd 2013). FST serves as a metric for assessing the variation in allele frequencies between populations, ranging between 0 and 1. Genetic diversity within each population was measured using two methods: (1) the average expected heterozygosity for all loci in a population (Nei 1973) and (2) the frequency of shared alleles among all individuals in a population (Chakraborty and Jin 1993).

Results

Genetic variation in LsCBF7 gene in lettuce

In our earlier research, we identified a family of fourteen CBF genes in lettuce through a comparative phylogenetic analysis (Park et al. 2020). Among these genes, LsCBF7 was found to have a premature stop within the AP2 DNA binding domain, resulting in a truncated and non-functional protein (Fig. 1). The expression analysis indicated that this mutant allele of LsCBF7 was significantly induced in response to low temperature (4 °C), implying a role in cold stress responses in lettuce (Park et al. 2020). Natural variation in CBF genes has been reported in multiple plants, often correlated with variation in freezing tolerance (Francia et al. 2007; Dumont et al. 2009; Knox et al. 2010; Tayeh et al. 2013; Gehan et al. 2015; Sieber et al. 2016). These observations prompted us to investigate whether this mutant allele is associated with variation in the freezing tolerance of lettuce. Our hypothesis was that lettuce accessions adapted to colder regions might have wild-type alleles of LsCBF7, and the wild-type alleles could confer enhanced freezing tolerance to the lettuce accessions. To address this hypothesis, we sequenced the genomic DNA of the LsCBF7 coding region for multiple accessions originating from colder climate, such as Sweden and Russia. Interestingly, our analysis identified a wild-type allele of LsCBF7 within one of the accessions from Sweden, namely PI 284702 (Fig. 1). A specific single base difference (A-to-G) at position 210 bp within the coding region of LsCBF7 resulted in a codon encoding tryptophan (W). This alteration led to the synthesis of a functional LsCBF7 protein.Fig. 1 Alignment of nucleotide and amino acid sequences for the LsCBF7 coding regions in the ‘Salinas’ and PI 284702 accessions. The predicted amino acid sequence is represented in single letter code, with the stop codon represented by an asterisk. Differences between the two accessions are highlighted in yellow; the A-to-G change at the 210th nucleotide results in a stop codon (*) changing to tryptophan (W). Red lines indicate the AP2 DNA-binding domain (color figure online)

To assess the potential impact of the putatively functional wild-type LsCBF7 allele on freezing tolerance, we conducted a freezing tolerance assay for PI 284702 (carrying the wild-type allele of LsCBF7) and ‘Salinas’ plants (carrying the mutant LsCBF7 allele). The freezing tolerance assay was performed with and without cold acclimation at 4 °C for 7 days to assess the influence of cold acclimation on freezing tolerance in lettuce. Our electrolyte leakage assays showed minimal disparities between non-acclimated PI 284702 and ‘Salinas’ plants, both exhibiting an EL50 of around –4 °C (the temperature at which freezing damage results in leakage of 50% of the total electrolytes) (Fig. 2a). However, after cold acclimation, the PI 284702 plants exhibited significantly greater freezing tolerance compared to the ‘Salinas’ plants. Specifically, the EL50 of cold-acclimated PI 284702 plants improved to − 8 °C, while ‘Salinas’ plants showed an EL50 of approximately − 6 °C (Fig. 2b). Therefore, the cold acclimation contributed to 4 °C enhancement in freezing tolerance (from − 4 to − 8 °C) for PI 284702 plants, whereas ‘Salinas’ plants exhibited 2 °C improvement (from − 4 to − 6 °C). The increase in freezing tolerance due to cold acclimation was reduced by about 50% in the ‘Salinas’ plants. These results suggested that the loss of CBF function might account for the 2 °C reduction in freezing tolerance, potentially underlying the variance in freezing tolerance between these two accessions.Fig. 2 The PI 284702 accession exhibits greater freezing tolerance compared to the Salinas accession. Electrolyte leakage assays were conducted on non-acclimated (NA) plants grown at 20 °C and on cold-acclimated (ACC) plants (4 °C for 7 days). Error bars indicate standard error of three biological replicates. An asterisk indicates statistically significant differences between two accessions (P < 0.05, student t-test)

Cold transcriptome analysis reveals LsGOLS1 gene as a primary target of LsCBF7

In response to low temperatures, CBF transcription factors play a pivotal role by upregulating an array of downstream target genes, which can orchestrate vital metabolic and biochemical processes that lead to the enhancement of freezing tolerance in plants (Gilmour et al. 2004; Cook et al. 2004; Kaplan et al. 2007). Therefore, the potential consequences of the loss-of-function mutation in LsCBF7 on cold acclimation can be reflected in the expression levels of its downstream target genes. To investigate the impact of the LsCBF7 mutation on the cold-induced expression of its targets, we performed RNA-seq experiments for PI 284702 and ‘Salinas’ plants subjected to low temperatures (4 °C) for 24 h and 7 days. The cold transcriptome analysis identified 15 genes that displayed significant differential expression between the two cultivars (Fig. 3; Table S5). Among these, one gene appeared to have substantial sequence similarity to a well-known CBF target—galactinol synthase, a gene widely recognized for its activation by CBF transcription factors in Arabidopsis (Park et al. 2015, 2018). In Arabidopsis, the GOLS genes belong to a multigene family with seven members of genes, with AtGOLS3 being the only gene significantly activated by CBFs in response to cold stress (Taji et al. 2002; Gehan et al. 2015).Fig. 3 Heatmap for the differentially expressed genes during cold acclimation in Salinas and PI 284702. It shows the log2-fold change in transcript levels for the upregulated genes in PI 284702 plants relative to ‘Salinas’ plants after exposure to low temperature (4 °C) for 24 h or 7 days. Transcript levels in cold conditions were normalized by those in warm conditions (20 °C). The color scale represents log2-fold change

To examine whether a similar regulatory relationship exists in lettuce, we first conducted a comprehensive search for GOLS genes within the lettuce genome. Our analysis indicated the presence of two GOLS genes in lettuce. Subsequently, we determined the orthologous relationships between the GOLS genes in Arabidopsis and lettuce through a comparative phylogenetic analysis. Our phylogenetic analysis revealed three ancestral groups within the GOLS gene family in higher plants, designated as I, II, and III (Fig. 4). All three groups contained genes from both Asterid and Rosid clades of flowering plants, indicating that ancestral genes of these groups predate the separation of Asterid and Rosid clades. The two lettuce GOLS genes, named LsGOLS1 and LsGOLS2, were assigned to groups I and III, respectively, while at least one of the GOLS genes from Arabidopsis was assigned to each of the three groups: AtGOLS2 and AtGOLS3 belonged to group I, AtGOLS1 to group II, and ATGOLS4, ATGOLS5, ATGOLS6, and ATGOLS7 belonged to group III. Therefore, based on the phylogenetic tree, LsGOLS1 (Ls8g23120) was closely related to both AtGOLS2 and AtGOLS3 that are downstream target of AtCBFs. On the other hand, LsGOLS2 (Ls9g201) was orthologous to four other Arabidopsis genes: AtGOLS4–7 (Table S2).Fig. 4 Phylogenetic analysis of 100 GOLS genes from 15 species, including lettuce and Arabidopsis. The tree was constructed based on the protein sequences using the neighbor-joining method. The nine species from the Asterid clade, including lettuce are marked with blue circles, and the six species from the Rosid clade, including Arabidopsis are marked with red circles. GOLS genes only from Lettuce and Arabidopsis are annotated in the tree (refer to Table S1 for genes from other species). Three ancestral groups are denoted I, II, and III

This orthologous relationship aligned with our RNA-seq results, where only LsGOLS1 (Ls8g23120) was identified as cold-induced in PI 284702 (Fig. 3). Subsequent qRT-PCR corroborated this, confirming significant cold induction of LsGOLS1 and minimal, if any, induction of LsGOLS2 in response to low temperatures. The degree of cold induction for LsGOLS1 was significantly higher in PI 284702 compared to ‘Salinas’, whereas LsGOLS2 exhibited no significant difference in cold induction between the two cultivars (Fig. 5). This result strongly supported the idea that LsGOLS1 is likely a direct target of LsCBF7, as in Arabidopsis, and the loss-of-function LsCBF7 mutation dampened the cold-induced expression of LsGOLS1 in ‘Salinas’ plants. Given the important role of GOLS genes in freezing tolerance, this finding supported the conclusion that the loss-of-function LsCBF7 mutation led to reduced freezing tolerance in ‘Salinas’ plants.Fig. 5 Expression of lettuce GOLS genes in response to low temperatures. Expression levels of LsGOLS1 and LsGOLS2 were determined by qRT-PCR in plants grown at 20 °C (0 h) and exposed to 4 °C for 4 h, 24 h and 7 days. Error bars represent standard error of three biological replicates

LsCBF7 allele frequency varies among lettuce horticultural types and wild lettuce

As the accession (PI 284702) with the wild-type LsCBF7 allele exhibited greater freezing tolerance compared to the accession (Salinas) with the mutant allele, we sought to investigate further, with a larger set of accessions, whether genetic variation in the LsCBF7 gene could contribute to difference in freezing tolerance. To this end, we genotyped LsCBF7 allele polymorphisms for 578 accessions, including four main horticultural types—115 crisphead, 118 romaine, 116 leaf, 136 butterhead—as well as 67 wild lettuce L. serriola accessions (Table S3). Our data revealed that the mutant LsCBF7 allele was predominant within the lettuce population, representing 86% (438) of the sampled L. sativa accessions, while only 14% (72) of the accessions carried the wild-type allele. Interestingly, the frequency of the wild-type allele varied among the four main horticultural types: 25% of leaf type accessions, 15% of romaine type, 10% of butterhead type, and 1% of crisphead type (Table 1). Furthermore, we observed a striking difference in the frequency of the LsCBF7 alleles among L. serriola accessions, the species that were commonly believed to be the progenitor of domesticated lettuce, L. sativa. In L. serriola, 90% (60) of accessions had the wild-type LsCBF7 allele, while only 9% (6) had the mutant allele, and one accession was heterozygous (Fig. 6). This significant discrepancy (binomial test, P value < 2.2e−16) suggested that the wild-type LsCBF7 allele may have been negatively selected during lettuce domestication or breeding programs.Table 1 Distribution of LsCBF7 allele genotype among four horticulturals types of lettuce

Horticultural type	Allele genotype*	% (Wt/total) (%)	
Wt	Mt	Ht	
Leaf	29	87	0	25	
Romaine	18	99	1	15	
Butterhead	14	122	0	10	
Crisphead	1	114	0	1	
*Wt represents wild-type; Mt, mutant; Ht, heterozygous

Fig. 6 Distribution of different LsCBF7 allele among cultivated lettuce (L. sativa) and wild lettuce (L. serriola). The numbers of accessions are shown on the pie chart, along with their respective percentages in parentheses

LsCBF7 allelic variation likely contributes to freezing tolerance variation among lettuce population

To investigate the association between genetic variation in the LsCBF7 gene and variation in freezing tolerance among lettuce accessions, we initially examined the cold-induced expression of LsGOLS genes in the additional accessions with different LsCBF7 genotypes. We randomly chose five accessions, each for the wild-type allele and the mutant allele of LsCBF7, and assessed the expression levels of LsGOLS1 and LsGOLS2 in plants exposed to a temperature of 4 °C for 24 h (Fig. 7). Consistent with our earlier findings from PI 284702 and ‘Salinas’ plants, the accessions with the wild-type LsCBF7 allele exhibited significant induction of LsGOLS1 expression following 24 h of cold treatment. In contrast, those carrying the mutant allele showed little or no induction of LsGOLS1. The expression levels of LsGOLS2, however, did not significantly differ between the two genotype groups (Fig. S1). These observations provided compelling evidence that LsGOLS1 is a primary downstream target of the LsCBF7 and its activation in response to low temperatures is dependent on LsCBF7 function.Fig. 7 LsGOLS1 expression in response to low temperatures in accessions carrying mutant or wild-type LsCBF7 alleles. Five accessions were randomly selected for each mutant (mt lscbf7) and wild-type (WT LsCBF7) allele group to assess LsGOLS1 expression. ‘Salinas’ and PI 284702 were included as controls. Accession identification numbers are shown on the x-axis. Gene expression was determined by qRT-PCR in plants grown at 20 °C (0 h) and exposed to 4 °C for 24 h. Error bars represent the standard error of three biological replicates

We further assessed the degree of freezing tolerance in these two genotype groups of accessions using electrolyte leakage assay. Different horticultural types have distinct growth characteristics, such as leaf texture and growth rate, which could potentially affect the assays. To minimize this possible bias, we selected accessions from the same horticultural types and tested them at a temperature range of − 5 to − 8 °C, which appeared to yield the most reliable results in electrolyte leakage assays at the earlier test with PI 284702 and ‘Salinas’ accessions (Fig. 2). In all six comparisons, the accessions carrying the wild-type LsCBF7 allele showed greater freezing tolerance following cold acclimation compared to those with the mutant allele (Fig. 8). These results further substantiated the idea that the wild-type LsCBF7 allele contributed to an increase in freezing tolerance in lettuce.Fig. 8 Loss-of-function mutation of LsCBF7 impairs freezing tolerance in lettuce plants. Electrolyte leakage freeze tests were conducted on cold-acclimated plants exposed to 4 °C for 7 days. Accessions from the butterhead type (a), leaf type (b), and romaine type (c) carrying either mutant (mt) or wild-type (WT) LsCBF7 allele were compared. Error bar indicates the standard error of three biological replicates. Significance was tested by Student’s t-test (**< 0.01, *< 0.05)

The LsCBF7 wild-type allele is associated with early bolting

Our genotyping analysis revealed a stark difference in the frequency of the wild-type LsCBF7 allele between L. sativa and L. serriola accessions (Fig. 6). Given that lettuce (L. sativa) is widely believed to be domesticated from L. serriola, this discrepancy led us to hypothesize that the wild-type allele might have been subjected to negative selection during domestication process or breeding efforts. Notably, we observed a strong negative correlation between bolting time and the wild-type allele frequency for the four lettuce horticultural types, with R2 of 0.96 (Fig. S2). We compared the average bolting time of the four horticultural types obtained from Park et al. (2021), with their LsCBF7 wild-type allele frequencies. The leaf type with the shortest average bolting time exhibited the highest frequency of the wild-type allele. The romaine type with the second-shortest bolting time showed the second-highest frequency, followed by butterhead and crisphead (Fig. S2). This finding suggested a potential association between the wild-type allele and early bolting trait.

Early bolting is generally considered an undesirable trait in lettuce cultivation, as it reduces the biomass of the leafy vegetable and leads to a bitter taste. Therefore, this association could be a driving force for the negative selection of the wild-type allele. To explore this idea, we referred to a phylogenetic tree previously constructed based on 186,000 SNP data with 441 accessions (Park et al. 2021), which were also genotyped for the LsCBF7 alleles in this study (Table S3). Using the neighbor-joining tree as a guide, we selected eight pairs of genetically closely related accessions that carried different LsCBF7 allele genotypes (Fig. S3). This approach was designed to minimize the influence of other genes on the bolting trait. We then grew these paired accessions under identical conditions and monitored their phenotypes. Our observations revealed that accessions carrying the wild-type LsCBF7 allele exhibited earlier bolting compared to those with mutant alleles in seven out of eight pairs, while one pair showed similar bolting time between the different genotypes (Fig. S4). This finding strongly supported the notion that the wild-type LsCBF7 allele is associated with the undesirable trait of early bolting, which could lead to the negative selection of the wild-type allele during domestication or breeding programs.

Genetic relationship between the two LsCBF7 allelic groups of L. sativa and wild lettuce (L. serriola)

To gain insights into the genetic relationship between the two LsCBF7 allelic groups of L. sativa and wild lettuce (L. serriola), we examined genetic similarity between these three groups and genetic diversity within the two allelic groups of L. sativa. The extent of genetic similarity was estimated using pairwise FST values, where lower FST indicates greater genetic similarity (Weir and Hill 2002). When the two allelic groups of L. sativa were compared to L. serriola accessions, the wild-type LsCBF7 accessions showed a lower FST value (0.242) than the mutant accessions (0.263). When FST was estimated only for the chromosome 9 region where LsCBF7 gene is located, difference between the two groups became more pronounced, with wild-type accessions displaying FST of 0.205 and mutant accessions having 0.286. These results indicated that accessions with the wild-type LsCBF7 allele were genetically closer to L. serriola than those with the mutant allele, particularly in the context of the chromosome 9 region.

In addition, we estimated genetic diversity within the two allelic groups of L. sativa. The extent of genetic diversity was estimated as expected heterozygosity (Hexp) (Nei 1973) and shared allele frequency (SAF) (Gao and Martin 2009), where higher Hexp and SAF indicates greater diversity within a population. The wild-type accessions exhibited higher Hexp (0.28) and SAF (0.54) values than the mutant accessions with Hexp of 0.26 and SAF of 0.50. When this estimation was limited to the chromosome 9 region, the genetic diversity for the wild-type accessions was even more pronounced. The wild-type accessions showed Hexp of 0.28 and SAF of 0.56, while the mutant accessions displayed Hexp of 0.23 and SAF of 0.44. These results indicated that accessions with the wild-type allele were genetically more diverse compared to accessions with the mutant allele.

Discussion

The CBF pathway plays a prominent role in cold acclimation and has been identified as a freezing tolerance QTL in many plant species (Francia et al. 2007; Dumont et al. 2009; Li et al. 2013; Tayeh et al. 2013; Sieber et al. 2016; Adhikari et al. 2021). Previously, we discovered a loss-of-function mutation in one of the lettuce CBF family genes (Park et al. 2020). In this study, we further investigated the potential impact of this mutant allele on freezing tolerance. Our results demonstrated that cold acclimation in lettuce can increase freezing tolerance by 4 °C, and the mutant allele might cause the reduction of this increase by 50% compared to its wild-type allele. Additionally, we unveiled variations in the frequency of the wild-type allele among horticultural types. Notably, these variations displayed a strong negative correlation with the average bolting time of these lettuce types, with the higher frequency of the wild-type allele leading to earlier bolting. This observation suggested an association between the wild-type allele and early bolting trait. This notion was further substantiated by bolting time tests conducted on genetically similar pairs of accessions but with different LsCBF7 alleles. Furthermore, we found a striking difference in the prevalence of the wild-type allele between cultivated lettuce (L.sativa) and wild lettuce (L. serriola). Only 14% of the L.sativa population carried the wild-type allele, while 90% of the L. serriola population carried the wild-type allele. Given that L. sativa is derived from L. serriola, our findings suggested that the wild-type allele underwent negative selection during the domestication of lettuce, despite its positive effect on freezing tolerance. We proposed that the negative selection of the wild-type allele might be driven by its association with early bolting. These results provide insights into how the genetic trade-off between freezing tolerance and early bolting, whether arising from linkage drag or pleiotropic effects, has played a role in the genetic make-up of lettuce during its domestication and breeding process.

Among the 14 CBF genes in lettuce, twelve were significantly induced in response to cold stress, suggesting their roles in freezing tolerance (Park et al. 2020). Nonetheless, the intriguing observation in this study is that a single CBF mutation in lettuce may significantly reduce the cold-induction of LsGOLS1. In Arabidopsis, the three CBFs are known to be functionally redundant, capable of inducing a similar set of target genes when overexpressed (Gilmour et al. 2004; Park et al. 2015). However, there is also evidence that they are not entirely redundant. Single or double mutations in the three Arabidopsis CBF genes still reduced freezing tolerance, although the effect was weaker compared to the triple CBF mutation (Zhao et al. 2016; Park et al. 2018), indicating that each of the three CBF transcription factors has its own non-redundant function. Consistently, CBF2 appears to have a more prominent influence on freezing tolerance than the other CBF genes in Arabidopsis. CBF2 overexpression activated the downstream target genes more strongly (Park et al. 2015), and double mutations including cbf2 (cbf1/cbf2 and cbf2/cbf3) showed a greater reduction in freezing tolerance than other double mutations without cbf2 (cbf1/cbf3) (Zhao et al. 2016). Therefore, these observations suggest that there may be some degree of non-redundancy among Arabidopsis CBF genes. In addition, when natural CBF mutations were surveyed across natural Arabidopsis populations, cbf2 mutants were restricted to much warmer regions than other CBF mutants (Monroe et al. 2016), indicating that CBF2 mutation has a more significant impact on freezing tolerance. Taken together, these observations support the notion that LsCBF7, like CBF2 in Arabidopsis, may play a more prominent role among the lettuce CBF family, and its mutation alone can lead to the significant reduction in the cold-induced expression of LsGOLS1.

Our results revealed intriguing variations in the frequency of the CBF wild-type allele across different horticultural lettuce types. Particularly in crisphead type, only one accession carried the wild-type LsCBF7 allele, while at least 10–25% of other horticultural types—butterhead, romaine, and leaf—possessed the wild-type allele. This distinction in the crisphead type may be attributed to a population bottleneck effect. Crisphead lettuce, the most recent development among the four types, originated from the French Batavia types in the twentieth century (de Vries 1997). Since then, most modern crisphead varieties have been developed in the United States from three cultivars (Calmar, Salinas, and Vanguard) as parental lines (Mikel 2007). Consistently, a study by Park et al. (2021) highlighted higher homogeneity within the crisphead type due to the extensive use of a small number of parents in a short period. Therefore, it is conceivable that the crisphead type may have originated from parents already possessing the mutant LsCBF7 allele, contributing to this frequency pattern. In addition, positive selection for delayed bolting may also play a role in maintaining the near-absence of the wild-type allele in the crisphead type. Compact head formation is a distinguishing feature of crisphead lettuce, and early bolting disrupts this head formation, resulting in more severe yield loss compared to other types. Therefore, the crisphead type is under stronger selection pressure for delayed bolting. In fact, the crisphead type exhibited the most pronounced delay in bolting among the four types (Park et al. 2021). Given the association of the wild-type allele with early bolting, this selection pressure would prevent the introduction of the wild-type allele into the crisphead population.

While the wild-type allele is dominant among L. serriola population, the presence of the mutant allele in some wild accessions raised a question about its origin. There are two possible scenarios: (1) the mutation occurred in L. serriola and was inherited by cultivated lettuce (L. sativa), or (2) the mutation occurred in cultivated lettuce and transferred back to L. serriola, as the two species are cross-compatible. We propose that the latter scenario is more likely. In nature, L. serriola plants typically germinate in late fall and overwinter in a rosette stage (Prince et al. 1978). Therefore, the wild-type LsCBF7 allele would provide a survival advantage during cold winter months, which subsequently allows for a high frequency of the wild-type allele within the population.

Considering the predominance of the wild-type allele in the L. serriola population, it is plausible that earlier domesticated lettuce would also carry the wild-type allele. Lettuce plants are typically cultivated during the warm growing season, avoiding the cold winter. Therefore, the beneficial effect of the wild-type allele on cultivated lettuce is not as pronounced as it is on wild lettuce, and the selection pressure on the wild-type allele would be relaxed. This relaxation may have allowed for the occurrence of the mutant allele, as observed in natural Arabidopsis populations adapted to warmer climates (Kang et al. 2013; Gehan et al. 2015; Monroe et al. 2016). Furthermore, the wild-type CBF allele might be less advantageous than the mutant allele in warmer regions that frequently experience temperatures low enough to induce CBF genes, but rarely experiencing freezing temperatures. False activation of CBF transcription factors could consume resources that could otherwise be utilized for growth, with the rare risk of freezing (Zhen et al. 2011; Oakley et al. 2014). Nevertheless, studies involving natural Arabidopsis populations suggest that the fitness benefit gained by the mutant allele may not be sufficient to account for its prevalence among cultivated lettuce in this study. For example, a survey for adaptive CBF mutations in warmer regions by Monroe et al. (2016) found an excess of nonsynonymous polymorphisms in CBF genes, but only 6% of the population carried nonfunctional mutations, with the functional CBF alleles still prevailing in the population. While previous studies (Ågren et al. 2013; Lee et al. 2024) found limited phenotypic effects of CBF loss-of-function mutations in warmer regions, primarily on fecundity, we observed a significant difference in bolting time between accessions carrying the wild-type and mutant alleles. This finding suggests that there could be a more influential factor contributing to the prevalence of the mutant allele within the lettuce population, such as delayed bolting.

Our bolting experiments (Figs. S2, S4) suggested that the mutant allele was associated with delayed bolting. Early bolting could be crucial for the survival and proliferation of wild lettuce, as it allows plants to produce seeds before the onset of unfavorable environmental conditions. On the other hand, delayed bolting in cultivated lettuce is a desirable trait with multiple benefits including increased yield and better quality, and this trait has been continuously selected in lettuce breeding programs (de Vries 1997; Hartman et al. 2013; Han et al. 2021). Therefore, during domestication or breeding, the association of the mutant allele with the delayed bolting trait accelerated the genetic fixation of the mutant allele among the modern lettuce population.

Our genetic diversity and similarity analyzes among the three genetic groups—L. serriola, accessions with the wild-type allele, and accessions with the mutant allele—also supported the idea that the wild-type allele was inherited from L. serriola and the mutant allele might be under positive selection in lettuce population. Our FST analysis indicated that accessions with the wild-type allele were genetically closer to L. serriola than those with the mutant allele (Table S6), making it more likely that the wild-type allele is passed down from L. serriola rather than the mutant allele. Moreover, the mutant accessions are genetically more homogeneous than the wild-type accessions, according to our Hexp and SAF analyzes (Table S6), suggesting that the mutant allele might be under strong selection, thereby resulting in reduced genetic diversity.

Although future empirical work will be needed to determine the nature of the association between the wild-type allele and early bolting—whether it is due to physical linkage or arises from a pleiotropic function of the wild-type allele—the utilization of this beneficial allele for improving freezing tolerance in lettuce would require decoupling the undesirable trait from the wild-type allele. Modern plant breeding techniques, such as CRISPR (clustered regularly interspaced short palindromic repeats), could also be utilized to modify the mutant allele and restore the cold tolerance trait in lettuce.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file 1 (DOCX 327 kb)

Supplementary file 2 (XLSX 58 kb)

Acknowledgements

We appreciate Dr. Lyndel W. Meinhardt for critical review of the manuscript. This research used resources provided by the SCINet project of the USDA Agricultural Research Service, ARS project number 0500-00093-001-00-D. Mention of a trade name, proprietary product, or vendor does not constitute an endorsement, guarantee, or warranty by the USDA and does not imply its approval to the exclusion of other products or vendors that may be suitable.

Author contribution statement

B.M. and A.S. obtained funding for this investigation. S.P. and B.M. conceptualized the study and designed experiments. S.P conducted the experiments, analyzed the results with the help of B.M., and wrote the first draft of the manuscript. The results were validated by B.M, and A.S, both of whom also contributed to writing the manuscript. All authors read and approved the final manuscript.

Funding

This research was supported by grants from the USDA-AMS Specialty Crop Multistate Program (Award No. 16SCCMAR0001), the USDA-NIFA Agriculture and Food Research Initiative (Award No. 2016-68004-24931), and the USDA-NIFA Specialty Crop Research Initiative (Award No. 2015-51181-24283, 2017-51181-26830, and 2021-51181-35903).

Data availability

The RNA-seq data are available in the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo/) under accession number GSE241604. All relevant data are included in the manuscript and the Supporting Information files.

Declarations

Conflict of interest

The authors have no relevant financial or non-financial interests to disclose.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Change history

9/23/2024

A Correction to this paper has been published: 10.1007/s00122-024-04727-5
==== Refs
References

Adhikari L Makaju SO Lindstrom OM Missaoui AM Mapping freezing tolerance QTL in alfalfa: based on indoor phenotyping BMC Plant Biol 2021 21 403 10.1186/s12870-021-03182-4 34488630
Adhikari L, Makaju SO, Lindstrom OM, Missaoui AM (2021) Mapping freezing tolerance QTL in alfalfa: based on indoor phenotyping. BMC Plant Biol 21:403. 10.1186/s12870-021-03182-434488630
Ågren J Schemske DW Reciprocal transplants demonstrate strong adaptive differentiation of the model organism Arabidopsis thaliana in its native range New Phytol 2012 194 1112 1122 10.1111/j.1469-8137.2012.04112.x 22432639
Ågren J, Schemske DW (2012) Reciprocal transplants demonstrate strong adaptive differentiation of the model organism Arabidopsis thaliana in its native range. New Phytol 194:1112–1122. 10.1111/j.1469-8137.2012.04112.x22432639
Ågren J Oakley CG McKay JK Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana Proc Natl Acad Sci 2013 110 21077 21082 10.1073/pnas.1316773110 24324156
Ågren J, Oakley CG, McKay JK et al (2013) Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana. Proc Natl Acad Sci 110:21077–21082. 10.1073/pnas.131677311024324156
Badawi M Danyluk J Boucho B The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs Mol Genet Genom 2007 277 533 554 10.1007/s00438-006-0206-9
Badawi M, Danyluk J, Boucho B et al (2007) The CBF gene family in hexaploid wheat and its relationship to the phylogenetic complexity of cereal CBFs. Mol Genet Genom 277:533–554. 10.1007/s00438-006-0206-9
Boling E (2022) Lettuce crops bouncing back after early January ice. In: The packer. https://www.thepacker.com/news/industry/lettuce-crops-bouncing-back-after-early-january-ice. Accessed 6 Nov 2023
Bradburd G (2013) Package ‘BEDASSLE.’ In: Compr. R Arch. Netw. https://cran.r-project.org/web//packages/BEDASSLE/BEDASSLE.pdf. Accessed 6 Nov 2023
Carvallo MA Pino M-T Jeknić Z A comparison of the low temperature transcriptomes and CBF regulons of three plant species that differ in freezing tolerance: Solanum commersonii, Solanum tuberosum, and Arabidopsis thaliana J Exp Bot 2011 62 3807 3819 10.1093/jxb/err066 21511909
Carvallo MA, Pino M-T, Jeknić Z et al (2011) A comparison of the low temperature transcriptomes and CBF regulons of three plant species that differ in freezing tolerance: Solanum commersonii, Solanum tuberosum, and Arabidopsis thaliana. J Exp Bot 62:3807–3819. 10.1093/jxb/err06621511909
Chakraborty R Jin L Pena SDJ Chakraborty R Epplen JT Jeffreys AJ A unified approach to study hypervariable polymorphisms: statistical considerations of determining relatedness and population distances DNA fingerprinting: state of the science 1993 Basel Birkhäuser 153 175
Chakraborty R, Jin L (1993) A unified approach to study hypervariable polymorphisms: statistical considerations of determining relatedness and population distances. In: Pena SDJ, Chakraborty R, Epplen JT, Jeffreys AJ (eds) DNA fingerprinting: state of the science. Birkhäuser, Basel, pp 153–175
Cook D Fowler S Fiehn O Thomashow MF A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis Proc Natl Acad Sci 2004 101 15243 15248 10.1073/pnas.0406069101 15383661
Cook D, Fowler S, Fiehn O, Thomashow MF (2004) A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis. Proc Natl Acad Sci 101:15243–15248. 10.1073/pnas.040606910115383661
Cunningham SM Nadeau P Castonguay Y Raffinose and stachyose accumulation, galactinol synthase expression, and winter injury of contrasting alfalfa germplasms Crop Sci 2003 43 562 570 10.2135/cropsci2003.5620
Cunningham SM, Nadeau P, Castonguay Y et al (2003) Raffinose and stachyose accumulation, galactinol synthase expression, and winter injury of contrasting alfalfa germplasms. Crop Sci 43:562–570. 10.2135/cropsci2003.5620
Dai H Zhu Z Wang Z Galactinol synthase 1 improves cucumber performance under cold stress by enhancing assimilate translocation Hortic Res 2022 9 uhab063 10.1093/hr/uhab063 35048123
Dai H, Zhu Z, Wang Z et al (2022) Galactinol synthase 1 improves cucumber performance under cold stress by enhancing assimilate translocation. Hortic Res 9:uhab063. 10.1093/hr/uhab06335048123
de Vries IM Origin and domestication of Lactuca sativa L. Genet Resour Crop Evol 1997 44 165 174 10.1023/A:1008611200727
de Vries IM (1997) Origin and domestication of Lactuca sativa L. Genet Resour Crop Evol 44:165–174. 10.1023/A:1008611200727
Dobin A Davis CA Schlesinger F STAR: ultrafast universal RNA-seq aligner Bioinformatics 2013 29 15 21 10.1093/bioinformatics/bts635 23104886
Dobin A, Davis CA, Schlesinger F et al (2013) STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29:15–21. 10.1093/bioinformatics/bts63523104886
Dumont E Fontaine V Vuylsteker C Association of sugar content QTL and PQL with physiological traits relevant to frost damage resistance in pea under field and controlled conditions Theor Appl Genet 2009 118 1561 1571 10.1007/s00122-009-1004-7 19322559
Dumont E, Fontaine V, Vuylsteker C et al (2009) Association of sugar content QTL and PQL with physiological traits relevant to frost damage resistance in pea under field and controlled conditions. Theor Appl Genet 118:1561–1571. 10.1007/s00122-009-1004-719322559
Francia E Barabaschi D Tondelli A Fine mapping of a HvCBF gene cluster at the frost resistance locus Fr-H2 in barley Theor Appl Genet 2007 115 1083 1091 10.1007/s00122-007-0634-x 17763839
Francia E, Barabaschi D, Tondelli A et al (2007) Fine mapping of a HvCBF gene cluster at the frost resistance locus Fr-H2 in barley. Theor Appl Genet 115:1083–1091. 10.1007/s00122-007-0634-x17763839
Gao X Martin ER Using allele sharing distance for detecting human population stratification Hum Hered 2009 68 182 191 10.1159/000224638 19521100
Gao X, Martin ER (2009) Using allele sharing distance for detecting human population stratification. Hum Hered 68:182–191. 10.1159/00022463819521100
Gehan MA Park S Gilmour SJ Natural variation in the C-repeat binding factor cold response pathway correlates with local adaptation of Arabidopsis ecotypes Plant J 2015 84 682 693 10.1111/tpj.13027 26369909
Gehan MA, Park S, Gilmour SJ et al (2015) Natural variation in the C-repeat binding factor cold response pathway correlates with local adaptation of Arabidopsis ecotypes. Plant J 84:682–693. 10.1111/tpj.1302726369909
Geisseler D, Horwath WR (2014) Lettuce production in California. In: Fertil. Res. Educ. Program. https://growyourhealthgardening.com/wp-content/uploads/2020/02/lettuce_production_ca.pdf. Accessed 6 Nov 2023
Gilmour SJ Fowler SG Thomashow MF Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities Plant Mol Biol 2004 54 767 781 10.1023/B:PLAN.0000040902.06881.d4 15356394
Gilmour SJ, Fowler SG, Thomashow MF (2004) Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol Biol 54:767–781. 10.1023/B:PLAN.0000040902.06881.d415356394
Hadi F Gilpin M Fuller MP Identification and expression analysis of CBF/DREB1 and COR15 genes in mutants of Brassica oleracea var. botrytis with enhanced proline production and frost resistance Plant Physiol Biochem 2011 49 1323 1332 10.1016/j.plaphy.2011.08.013 22000056
Hadi F, Gilpin M, Fuller MP (2011) Identification and expression analysis of CBF/DREB1 and COR15 genes in mutants of Brassica oleracea var. botrytis with enhanced proline production and frost resistance. Plant Physiol Biochem 49:1323–1332. 10.1016/j.plaphy.2011.08.01322000056
Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In: Nucleic acids symposium series. Oxford, pp 95–98
Han R Truco MJ Lavelle DO Michelmore RW A composite analysis of flowering time regulation in lettuce Front Plant Sci 2021 12 632708 10.3389/fpls.2021.632708 33763095
Han R, Truco MJ, Lavelle DO, Michelmore RW (2021) A composite analysis of flowering time regulation in lettuce. Front Plant Sci 12:632708. 10.3389/fpls.2021.63270833763095
Hartman Y Hooftman DAP Eric Schranz M van Tienderen PH QTL analysis reveals the genetic architecture of domestication traits in Crisphead lettuce Genet Resour Crop Evol 2013 60 1487 1500 10.1007/s10722-012-9937-0
Hartman Y, Hooftman DAP, Eric Schranz M, van Tienderen PH (2013) QTL analysis reveals the genetic architecture of domestication traits in Crisphead lettuce. Genet Resour Crop Evol 60:1487–1500. 10.1007/s10722-012-9937-0
Kang J Zhang H Sun T Natural variation of C-repeat-binding factor (CBFs) genes is a major cause of divergence in freezing tolerance among a group of Arabidopsis thaliana populations along the Yangtze River in China New Phytol 2013 199 1069 1080 10.1111/nph.12335 23721132
Kang J, Zhang H, Sun T et al (2013) Natural variation of C-repeat-binding factor (CBFs) genes is a major cause of divergence in freezing tolerance among a group of Arabidopsis thaliana populations along the Yangtze River in China. New Phytol 199:1069–1080. 10.1111/nph.1233523721132
Kaplan F Kopka J Sung DY Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content Plant J 2007 50 967 981 10.1111/j.1365-313X.2007.03100.x 17461790
Kaplan F, Kopka J, Sung DY et al (2007) Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. Plant J 50:967–981. 10.1111/j.1365-313X.2007.03100.x17461790
Knaupp M Mishra KB Nedbal L Heyer AG Evidence for a role of raffinose in stabilizing photosystem II during freeze–thaw cycles Planta 2011 234 477 486 10.1007/s00425-011-1413-0 21533754
Knaupp M, Mishra KB, Nedbal L, Heyer AG (2011) Evidence for a role of raffinose in stabilizing photosystem II during freeze–thaw cycles. Planta 234:477–486. 10.1007/s00425-011-1413-021533754
Knox AK Dhillon T Cheng H CBF gene copy number variation at frost resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals Theor Appl Genet 2010 121 21 35 10.1007/s00122-010-1288-7 20213518
Knox AK, Dhillon T, Cheng H et al (2010) CBF gene copy number variation at frost resistance-2 is associated with levels of freezing tolerance in temperate-climate cereals. Theor Appl Genet 121:21–35. 10.1007/s00122-010-1288-720213518
Lee G Sanderson BJ Ellis TJ A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation Proc Natl Acad Sci 2024 121 e2317461121 10.1073/pnas.2317461121 38289961
Lee G, Sanderson BJ, Ellis TJ et al (2024) A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation. Proc Natl Acad Sci 121:e2317461121. 10.1073/pnas.231746112138289961
Li L Liu X Xie K qLTG-9, a stable quantitative trait locus for low-temperature germination in rice (Oryza sativa L.) Theor Appl Genet 2013 126 2313 2322 10.1007/s00122-013-2137-2 23748708
Li L, Liu X, Xie K et al (2013) qLTG-9, a stable quantitative trait locus for low-temperature germination in rice (Oryza sativa L.). Theor Appl Genet 126:2313–2322. 10.1007/s00122-013-2137-223748708
Li T Zhang Y Liu Y Raffinose synthase enhances drought tolerance through raffinose synthesis or galactinol hydrolysis in maize and Arabidopsis plants J Biol Chem 2020 295 8064 8077 10.1074/jbc.RA120.013948 32366461
Li T, Zhang Y, Liu Y et al (2020) Raffinose synthase enhances drought tolerance through raffinose synthesis or galactinol hydrolysis in maize and Arabidopsis plants. J Biol Chem 295:8064–8077. 10.1074/jbc.RA120.01394832366461
Liao Y Smyth GK Shi W featureCounts: an efficient general purpose program for assigning sequence reads to genomic features Bioinformatics 2014 30 923 930 10.1093/bioinformatics/btt656 24227677
Liao Y, Smyth GK, Shi W (2014) featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923–930. 10.1093/bioinformatics/btt65624227677
Liu Y Zhang L Meng S Expression of galactinol synthase from Ammopiptanthus nanus in tomato improves tolerance to cold stress J Exp Bot 2020 71 435 449 10.1093/jxb/erz450 31616940
Liu Y, Zhang L, Meng S et al (2020) Expression of galactinol synthase from Ammopiptanthus nanus in tomato improves tolerance to cold stress. J Exp Bot 71:435–449. 10.1093/jxb/erz45031616940
Luna N (2013) Lettuce prices jump due to cold snap. In: Orange Cty. Regist. https://www.ocregister.com/2013/01/23/lettuce-prices-jump-due-to-cold-snap/. Accessed 6 Nov 2023
McKhann HI Gery C Bérard A Natural variation in CBF gene sequence, gene expression and freezing tolerance in the Versailles core collection of Arabidopsis thaliana BMC Plant Biol 2008 8 105 10.1186/1471-2229-8-105 18922165
McKhann HI, Gery C, Bérard A et al (2008) Natural variation in CBF gene sequence, gene expression and freezing tolerance in the Versailles core collection of Arabidopsis thaliana. BMC Plant Biol 8:105. 10.1186/1471-2229-8-10518922165
Mikel MA Genealogy of contemporary north American lettuce HortScience 2007 42 489 493 10.21273/HORTSCI.42.3.489
Mikel MA (2007) Genealogy of contemporary north American lettuce. HortScience 42:489–493. 10.21273/HORTSCI.42.3.489
Monroe JG McGovern C Lasky JR Adaptation to warmer climates by parallel functional evolution of CBF genes in Arabidopsis thaliana Mol Ecol 2016 25 3632 3644 10.1111/mec.13711 27247130
Monroe JG, McGovern C, Lasky JR et al (2016) Adaptation to warmer climates by parallel functional evolution of CBF genes in Arabidopsis thaliana. Mol Ecol 25:3632–3644. 10.1111/mec.1371127247130
Mou B Prohens J Nuez F Lettuce Vegetables I: asteraceae, brassicaceae, chenopodicaceae, and cucurbitaceae 2008 New York, NY Springer 75 116
Mou B (2008) Lettuce. In: Prohens J, Nuez F (eds) Vegetables I: asteraceae, brassicaceae, chenopodicaceae, and cucurbitaceae. Springer, New York, NY, pp 75–116
Nei M Analysis of gene diversity in subdivided populations Proc Natl Acad Sci U S A 1973 70 3321 3323 10.1073/pnas.70.12.3321 4519626
Nei M (1973) Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci U S A 70:3321–3323. 10.1073/pnas.70.12.33214519626
Nishizawa-Yokoi A Yoshida E Yabuta Y Shigeoka S Analysis of the regulation of target genes by an Arabidopsis heat shock transcription factor, HsfA2 Biosci Biotechnol Biochem 2009 73 890 895 10.1271/bbb.80809 19352026
Nishizawa-Yokoi A, Yoshida E, Yabuta Y, Shigeoka S (2009) Analysis of the regulation of target genes by an Arabidopsis heat shock transcription factor, HsfA2. Biosci Biotechnol Biochem 73:890–895. 10.1271/bbb.8080919352026
Oakley CG Ågren J Atchison RA Schemske DW QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs Mol Ecol 2014 23 4304 4315 10.1111/mec.12862 25039860
Oakley CG, Ågren J, Atchison RA, Schemske DW (2014) QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs. Mol Ecol 23:4304–4315. 10.1111/mec.1286225039860
Owens CL Thomashow MF Hancock JF Iezzoni AF CBF1 orthologs in sour cherry and strawberry and the heterologous expression of CBF1 in strawberry J Am Soc Hortic Sci 2002 127 489 494 10.21273/JASHS.127.4.489
Owens CL, Thomashow MF, Hancock JF, Iezzoni AF (2002) CBF1 orthologs in sour cherry and strawberry and the heterologous expression of CBF1 in strawberry. J Am Soc Hortic Sci 127:489–494. 10.21273/JASHS.127.4.489
Panikulangara TJ Eggers-Schumacher G Wunderlich M Galactinol synthase1. A novel heat shock factor target gene responsible for heat-induced synthesis of raffinose family oligosaccharides in Arabidopsis Plant Physiol 2004 136 3148 3158 10.1104/pp.104.042606 15466240
Panikulangara TJ, Eggers-Schumacher G, Wunderlich M et al (2004) Galactinol synthase1. A novel heat shock factor target gene responsible for heat-induced synthesis of raffinose family oligosaccharides in Arabidopsis. Plant Physiol 136:3148–3158. 10.1104/pp.104.04260615466240
Park S Lee C-M Doherty CJ Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network Plant J Cell Mol Biol 2015 82 193 207 10.1111/tpj.12796
Park S, Lee C-M, Doherty CJ et al (2015) Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J Cell Mol Biol 82:193–207. 10.1111/tpj.12796
Park S Gilmour SJ Grumet R Thomashow MF CBF-dependent and CBF-independent regulatory pathways contribute to the differences in freezing tolerance and cold-regulated gene expression of two Arabidopsis ecotypes locally adapted to sites in Sweden and Italy PLoS ONE 2018 13 e0207723 10.1371/journal.pone.0207723 30517145
Park S, Gilmour SJ, Grumet R, Thomashow MF (2018) CBF-dependent and CBF-independent regulatory pathways contribute to the differences in freezing tolerance and cold-regulated gene expression of two Arabidopsis ecotypes locally adapted to sites in Sweden and Italy. PLoS ONE 13:e0207723. 10.1371/journal.pone.020772330517145
Park S Shi A Mou B Genome-wide identification and expression analysis of the CBF/DREB1 gene family in lettuce Sci Rep 2020 10 5733 10.1038/s41598-020-62458-1 32235838
Park S, Shi A, Mou B (2020) Genome-wide identification and expression analysis of the CBF/DREB1 gene family in lettuce. Sci Rep 10:5733. 10.1038/s41598-020-62458-132235838
Park S Kumar P Shi A Mou B Population genetics and genome-wide association studies provide insights into the influence of selective breeding on genetic variation in lettuce Plant Genome 2021 14 e20086 10.1002/tpg2.20086 33629537
Park S, Kumar P, Shi A, Mou B (2021) Population genetics and genome-wide association studies provide insights into the influence of selective breeding on genetic variation in lettuce. Plant Genome 14:e20086. 10.1002/tpg2.2008633629537
Park S Shi A Meinhardt LW Mou B Genome-wide characterization and evolutionary analysis of the AP2/ERF gene family in lettuce (Lactuca sativa) Sci Rep 2023 13 21990 10.1038/s41598-023-49245-4 38081919
Park S, Shi A, Meinhardt LW, Mou B (2023) Genome-wide characterization and evolutionary analysis of the AP2/ERF gene family in lettuce (Lactuca sativa). Sci Rep 13:21990. 10.1038/s41598-023-49245-438081919
Pennycooke JC Jones ML Stushnoff C Down-regulating α-galactosidase enhances freezing tolerance in transgenic petunia Plant Physiol 2003 133 901 909 10.1104/pp.103.024554 14500789
Pennycooke JC, Jones ML, Stushnoff C (2003) Down-regulating α-galactosidase enhances freezing tolerance in transgenic petunia. Plant Physiol 133:901–909. 10.1104/pp.103.02455414500789
Prince SD Marks MK Carter RN Induction of flowering in wild lettuce (Lactuca serriola L.) New Phytol 1978 81 265 277 10.1111/j.1469-8137.1978.tb02632.x
Prince SD, Marks MK, Carter RN (1978) Induction of flowering in wild lettuce (Lactuca serriola L.). New Phytol 81:265–277. 10.1111/j.1469-8137.1978.tb02632.x
R core Team (2023) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Retrieved from https://www.R-project.org. In: Lang. Environ. Stat. Comput. Vienna Austria R Found. Stat. Comput. https://www.r-project.org/. Accessed 23 Jun 2023
Reyes-Chin-Wo S Wang Z Yang X Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce Nat Commun 2017 8 14953 10.1038/ncomms14953 28401891
Reyes-Chin-Wo S, Wang Z, Yang X et al (2017) Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce. Nat Commun 8:14953. 10.1038/ncomms1495328401891
Robinson MD McCarthy DJ Smyth GK edgeR: a Bioconductor package for differential expression analysis of digital gene expression data Bioinformatics 2010 26 139 140 10.1093/bioinformatics/btp616 19910308
Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139–140. 10.1093/bioinformatics/btp61619910308
Sakai A Larcher W Sakai A Larcher W Regional distribution of plants and their adaptive responses to low temperatures Frost survival of plants: responses and adaptation to freezing stress 1987 Berlin, Heidelberg Springer 174 234
Sakai A, Larcher W (1987) Regional distribution of plants and their adaptive responses to low temperatures. In: Sakai A, Larcher W (eds) Frost survival of plants: responses and adaptation to freezing stress. Springer, Berlin, Heidelberg, pp 174–234
Saravitz DM Pharr DM Carter TE Jr Galactinol synthase activity and soluble sugars in developing seeds of four soybean genotypes 1 Plant Physiol 1987 83 185 189 10.1104/pp.83.1.185 16665199
Saravitz DM, Pharr DM, Carter TE Jr (1987) Galactinol synthase activity and soluble sugars in developing seeds of four soybean genotypes 1. Plant Physiol 83:185–189. 10.1104/pp.83.1.18516665199
Schramm F Larkindale J Kiehlmann E A cascade of transcription factor DREB2A and heat stress transcription factor HsfA3 regulates the heat stress response of Arabidopsis Plant J 2008 53 264 274 10.1111/j.1365-313X.2007.03334.x 17999647
Schramm F, Larkindale J, Kiehlmann E et al (2008) A cascade of transcription factor DREB2A and heat stress transcription factor HsfA3 regulates the heat stress response of Arabidopsis. Plant J 53:264–274. 10.1111/j.1365-313X.2007.03334.x17999647
Sieber A-N Longin CFH Leiser WL Würschum T Copy number variation of CBF-A14 at the Fr-A2 locus determines frost tolerance in winter durum wheat Theor Appl Genet 2016 129 1087 1097 10.1007/s00122-016-2685-3 26883046
Sieber A-N, Longin CFH, Leiser WL, Würschum T (2016) Copy number variation of CBF-A14 at the Fr-A2 locus determines frost tolerance in winter durum wheat. Theor Appl Genet 129:1087–1097. 10.1007/s00122-016-2685-326883046
Skinner JS von Zitzewitz J Szűcs P Structural, functional, and phylogenetic characterization of a large CBF gene family in barley Plant Mol Biol 2005 59 533 551 10.1007/s11103-005-2498-2 16244905
Skinner JS, von Zitzewitz J, Szűcs P et al (2005) Structural, functional, and phylogenetic characterization of a large CBF gene family in barley. Plant Mol Biol 59:533–551. 10.1007/s11103-005-2498-216244905
Stockinger EJ Gilmour SJ Thomashow MF Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit Proc Natl Acad Sci U S A 1997 94 1035 1040 10.1073/pnas.94.3.1035 9023378
Stockinger EJ, Gilmour SJ, Thomashow MF (1997) Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc Natl Acad Sci U S A 94:1035–1040. 10.1073/pnas.94.3.10359023378
Taji T Ohsumi C Iuchi S Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana Plant J 2002 29 417 426 10.1046/j.0960-7412.2001.01227.x 11846875
Taji T, Ohsumi C, Iuchi S et al (2002) Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J 29:417–426. 10.1046/j.0960-7412.2001.01227.x11846875
Tamura K Stecher G Kumar S MEGA11: molecular evolutionary genetics analysis version 11 Mol Biol Evol 2021 38 3022 3027 10.1093/molbev/msab120 33892491
Tamura K, Stecher G, Kumar S (2021) MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol 38:3022–3027. 10.1093/molbev/msab12033892491
Tayeh N Bahrman N Sellier H A tandem array of CBF/DREB1 genes is located in a major freezing tolerance QTL region on Medicago truncatula chromosome 6 BMC Genom 2013 14 814 10.1186/1471-2164-14-814
Tayeh N, Bahrman N, Sellier H et al (2013) A tandem array of CBF/DREB1 genes is located in a major freezing tolerance QTL region on Medicago truncatula chromosome 6. BMC Genom 14:814. 10.1186/1471-2164-14-814
Thomashow MF PLANT COLD ACCLIMATION: freezing tolerance genes and regulatory mechanisms Annu Rev Plant Physiol Plant Mol Biol 1999 50 571 599 10.1146/annurev.arplant.50.1.571 15012220
Thomashow MF (1999) PLANT COLD ACCLIMATION: freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol 50:571–599. 10.1146/annurev.arplant.50.1.57115012220
Thomashow MF Molecular basis of plant cold acclimation: insights gained from studying the CBF cold response pathway Plant Physiol 2010 154 571 577 10.1104/pp.110.161794 20921187
Thomashow MF (2010) Molecular basis of plant cold acclimation: insights gained from studying the CBF cold response pathway. Plant Physiol 154:571–577. 10.1104/pp.110.16179420921187
Weir BS Hill WG Estimating F-statistics Annu Rev Genet 2002 36 721 750 10.1146/annurev.genet.36.050802.093940 12359738
Weir BS, Hill WG (2002) Estimating F-statistics. Annu Rev Genet 36:721–750. 10.1146/annurev.genet.36.050802.09394012359738
Yadav S Modi P Dave A Effect of abiotic stress on crops Sustain Crop Prod 2020 3 5 16
Yadav S, Modi P, Dave A et al (2020) Effect of abiotic stress on crops. Sustain Crop Prod 3:5–16
Zhang H Sun Z Feng S The C2H2-type zinc finger protein PhZFP1 regulates cold stress tolerance by modulating galactinol synthesis in Petunia hybrida J Exp Bot 2022 73 6434 6448 10.1093/jxb/erac274 35726094
Zhang H, Sun Z, Feng S et al (2022) The C2H2-type zinc finger protein PhZFP1 regulates cold stress tolerance by modulating galactinol synthesis in Petunia hybrida. J Exp Bot 73:6434–6448. 10.1093/jxb/erac27435726094
Zhao C Zhang Z Xie S Mutational evidence for the critical role of CBF transcription factors in cold acclimation in Arabidopsis Plant Physiol 2016 171 2744 2759 10.1104/pp.16.00533 27252305
Zhao C, Zhang Z, Xie S et al (2016) Mutational evidence for the critical role of CBF transcription factors in cold acclimation in Arabidopsis. Plant Physiol 171:2744–2759. 10.1104/pp.16.0053327252305
Zhen Y Ungerer MC Relaxed selection on the CBF/DREB1 regulatory genes and reduced freezing tolerance in the southern range of Arabidopsis thaliana Mol Biol Evol 2008 25 2547 2555 10.1093/molbev/msn196 18775899
Zhen Y, Ungerer MC (2008) Relaxed selection on the CBF/DREB1 regulatory genes and reduced freezing tolerance in the southern range of Arabidopsis thaliana. Mol Biol Evol 25:2547–2555. 10.1093/molbev/msn19618775899
Zhen Y Dhakal P Ungerer MC Fitness benefits and costs of cold acclimation in Arabidopsis thaliana Am Nat 2011 178 44 52 10.1086/660282 21670576
Zhen Y, Dhakal P, Ungerer MC (2011) Fitness benefits and costs of cold acclimation in Arabidopsis thaliana. Am Nat 178:44–52. 10.1086/66028221670576
Zhuo C Wang T Lu S A cold responsive galactinol synthase gene from Medicago falcata (MfGolS1) is induced by myo-inositol and confers multiple tolerances to abiotic stresses Physiol Plant 2013 149 67 78 10.1111/ppl.12019 23253102
Zhuo C, Wang T, Lu S et al (2013) A cold responsive galactinol synthase gene from Medicago falcata (MfGolS1) is induced by myo-inositol and confers multiple tolerances to abiotic stresses. Physiol Plant 149:67–78. 10.1111/ppl.1201923253102
