Locally adaptive temperature response of vegetative growth in Arabidopsis thaliana

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

We investigated early vegetative growth of natural Arabidopsis thaliana accessions in cold, nonfreezing temperatures, similar to temperatures these plants naturally encounter in fall at northern latitudes. We found that accessions from northern latitudes produced larger seedlings than accessions from southern latitudes, partly as a result of larger seed size. However, their subsequent vegetative growth when exposed to colder temperatures was slower. The difference was too large to be explained by random population differentiation, and is thus suggestive of local adaptation, a notion that is further supported by substantial transcriptome and metabolome changes in northern accessions. We hypothesize that the reduced growth of northern accessions is an adaptive response and a consequence of reallocating resources toward cold acclimation and winter survival.

Editor's evaluation

The article combines genetic and phenotypic approaches to show convincing evidence of local adaptation in early vegetative growth of Arabidopsis lineages sampled from a wide range of locations. The authors show larger initial size and slower growth of northern accessions compared to southern accessions when exposed to cold temperatures, suggesting that northern accessions potentially reallocate resources for winter survival. This study is commendable for its scope and comprehensive analysis of local adaptation of a highly polygenic trait in a model weed.

https://doi.org/10.7554/eLife.77913.sa0

Introduction

Plants use a wide variety of life history strategies in adaptation to their local environment. These strategies have evolved to maximize fitness, but are constrained by trade-offs between components such as growth, survival, and reproduction (Lande, 1982; Stearns, 1992). While most life history studies investigate differences between species, there is also variation found within species, including in Arabidopsis thaliana, where life history variation has been linked to climate parameters (Estarague et al., 2022; Vasseur et al., 2018; Sartori et al., 2019). Less clear is how trade-offs are shaping this variation. In this article, we consider vegetative growth, a key component of life history, and use transcriptome and metabolome data to help explore potential trade-offs. We did this specifically in cold temperatures meant to simulate natural conditions in the northern regions of the species distribution.

Local adaptation studies in A. thaliana have found important roles for life history traits such as seed dormancy and flowering time (Takou et al., 2019). Temperature is a major regulator of these traits and local populations are adapted to their local climate (Martínez-Berdeja et al., 2020; Simpson and Dean, 2002; Hepworth et al., 2018). Plant growth is also affected by temperature, and previous studies have detected genetic variation underlying growth-related traits (Bac-Molenaar et al., 2015; Marchadier et al., 2019), as well as signals of polygenic adaptation (Wieters et al., 2021). How vegetative growth is adapted to local temperatures remains unclear, however.

Growth can be seen as the end sum of a vast number of physiological processes. All of these are genetically determined but can also be heavily influenced by environmental conditions (Bac-Molenaar et al., 2015; Fritz et al., 2018). Growth is therefore not only genetically a complex trait but also enormously plastic. The most straightforward environmental effect is when conditions are so adverse that growth reaches a physiological limit, making it impossible for the plant to grow any further. This is called ‘passive plasticity’ (Forsman, 2015; van Kleunen and Fischer, 2005). Yet, when survival is at stake, it may also be in the interest of the plant to actively inhibit growth upon deteriorating environmental conditions (Claeys and Inzé, 2013), called ‘active plasticity’ (Forsman, 2015; van Kleunen and Fischer, 2005). Since vegetative growth ultimately determines photosynthetic surface and thus energy input that can be invested in the next generation, it is directly related to fitness, and is typically in trade-off with survival. Allocation of resources towards either growth or survival is thus an important balance to keep, and plants are expected to be adapted to constantly perceiving and responding to specific environmental changes as cues for looming adverse conditions.

Cold acclimation is a well-studied mechanism in plants, in which decreasing temperatures induce freezing tolerance in preparation for winter (Thomashow, 1999; Hughes and Dunn, 1996). This temperature response is typically studied at 4°C, but has been observed in temperatures up to 12°C (Bond et al., 2011). The increased freezing tolerance is accomplished by changing membrane composition, producing cryoprotective polypeptides such as COR15A (Artus et al., 1996; Steponkus et al., 1998) and accumulating compatible solutes with cryoprotective properties such as raffinose, sucrose, and proline (Nanjo et al., 1999; Gilmour et al., 2000; Taji et al., 2002). Main regulators of cold acclimation are CBF1/DREB1b, CBF2/DREB1c, and CBF3/DREB1a, three AP2/ERF transcription factors, for which allelic variation in CBF2 has been linked to natural variation in freezing tolerance (Oakley et al., 2014; Park et al., 2018; Alonso-Blanco et al., 2005). Nothing is known, however, about whether natural variation in freezing tolerance regulation also influences a trade-off with growth responses to cold temperatures.

Here, we investigated the role of growth in adaptation to cold temperatures by comparing vegetative growth of 249 accessions (Figure 1) grown in daily maximum temperatures of 16 and 6°C for a period of 3 weeks following seedling establishment (Figure 1—figure supplement 1). Rosette growth of each plant was measured twice a day during temperature treatments using automated phenotyping. The experiment generated rosette growth estimates at a high temporal resolution in two ecologically realistic temperature conditions in a wide set of accessions, allowing us to look for patterns of local adaptation.

Figure 1 with 1 supplement see all

Download asset Open asset

Geographic origin of the 249 accessions.

Map color shows winter temperature (mean temperature of coldest quarter). Accessions are colored according to subpopulation (1001 Genomes Consortium, 2016). Accessions from the warmest and coldest regions are from Greece and the Himalayas, respectively.

Figure 1—source data 1 List of all 249 accessions with indication of the 8 accessions used for the transcriptome analysis.: https://cdn.elifesciences.org/articles/77913/elife-77913-fig1-data1-v1.csv
Download elife-77913-fig1-data1-v1.csv

Results

Estimating plant growth

Our highly replicated experiment yielded dense (two measurements per day) time-series growth data for over 7000 individual plants (5 replicate plants × 249 accessions × 2 treatments × 3 replicate experiments). These data were used to model plant growth and estimate growth parameters for further analysis. Unlimited growth should be exponential, but plant growth is known to slow down with increasing size, and therefore a power-law function, $\frac{d M}{d t} = r M^{β}$ , with $β < 1$ is typically a better fit than a pure exponential function (for which $β = 1$ — in the equation, M is the size, r is the growth rate, and $β$ is a scaling factor that allows rate of size increase to change with size). Growth according to a power-law function typically describes early stages of plant growth especially well (Paine et al., 2012), and our rosette size measurements were no exception (see ‘Materials and methods’). To calculate the rosette size from a power-law function at a given time point, only three parameters are required: the initial size ( $M_{0})$ , growth rate ( $r)$ , and $β$ . Note that ( $M_{0})$ is the rosette size at the start of the temperature treatment 14 days after stratification (Figure 1—figure supplement 1) and is thus not affected by the temperature treatment. We used a nonlinear mixed model to obtain estimates for the initial size, growth rates, and $β$ . Accession was added as fixed effect for initial size and growth rate, temperature and accession × temperature interactions were added as fixed effects for growth rate only. $β$ was considered to be constant over accessions and temperatures. The ‘temperature response’ of the growth rate was calculated for each accession as the slope between the growth rate at 16 and 6°C. As expected, all accessions grew faster when it was warmer. The observed phenotypic variation (Figure 2—figure supplement 1) is to a large extent explained by genetic variation; broad-sense heritabilities are 0.41 for initial size, and 0.57 and 0.32 for growth rate at 16 and 6°C, respectively.

Growth parameters correlate with the environment of origin

If growth rates are locally adaptive, they may reflect the environment of origin of each accession. To investigate this, we correlated our estimated growth rates with climate data. The climate variables showing the strongest correlations with the different growth parameters were linked to winter temperatures (Figure 2—figure supplement 2), also when correlations were corrected for population structure (Figure 2—figure supplement 2). In particular, the mean temperature during the coldest quarter (henceforth referred to as ‘winter temperature’) was most strongly correlated with our parameter estimates, and we focus on it in what follows.

Initial size

Accessions from colder climates generally had higher initial rosette size (M₀), 2 weeks after germination, than accessions from warmer climates ( $r = - 0.39$ ), but then grew more slowly during the temperature experiment – regardless of temperature regime (Figure 2, Figure 2—figure supplement 3). Because during the first 2 weeks the plants were growing at 21°C, it is impossible to disentangle early growth from growth at a warmer temperature.

Figure 2 with 6 supplements see all

Download asset Open asset

Correlations of growth parameters with winter temperature.

(A) Initial size. (B) Growth rate at 16°C. (C) Growth rate at 6°C. (D) Temperature response of growth rate. Colors indicate genetically defined subpopulations of the accessions (1001 Genomes Consortium, 2016).

One reason for this pattern is likely to be the differences in seed size between accessions. Using unpublished seed size measurements for a subset of 123 Swedish accessions from previous experiments, we found that seed size is positively correlated with initial size ( $r = 0.28$ ; Figure 2—figure supplement 5A), and also with winter temperature ( $r = - 0.75$ ; Figure 2—figure supplement 5B), at least for the subset of 123 Swedish accessions. In a random-effect model, winter temperature explained 32.7% of the variation in initial size, whereas seed size explained 11.9%. Winter temperature is still significantly associated with initial size when seed size is taken into account ( $p - v a l u e < 1 e - 04$ ). The precise role of seed size in initial growth is surely a topic that would benefit from further studies, but for the purpose of this study, it is clear that seed size alone cannot explain the geographic pattern we observe for M₀ and that there must be a role for variation in growth rate during the very initial phases of seedling growth.

Growth rates

While the initial sizes correlate negatively with winter temperature, we observed the opposite relation for the growth rates. Despite being larger initially, accessions from colder climates grew more slowly during both the 16°C ( $r = 0.33$ ) and 6°C ( $r = 0.28$ ) treatments (Figure 2). The higher growth rates of accessions from warmer climates prove that resources are not limiting, suggesting that the northern lines are actively inhibiting their growth, and that growing slower may be beneficial in colder climates, perhaps in preparation for winter. Accessions from colder climates were also less sensitive to the temperature experiment in the sense that the temperature response of the growth rate increased with winter temperature of origin ( $r = 0.28$ ; Figure 2D). Even though accessions from the Asian and north Swedish subpopulations were more variable in their growth rate temperature response (Figure 2—figure supplement 4A), the correlations still hold when removing either Asian or northern and southern Swedish subpopulations (Figure 2—figure supplement 6B, Figure 2—figure supplement 4D), and when looking specifically within the northern and southern Swedish subpopulations (Figure 2—figure supplement 6).

Cold acclimation response

Just like the observed geographic pattern of the growth rates, metabolite measurements taken at the final day of our experiment and presented in an earlier publication (Weiszmann et al., 2020) showed clear differences between accessions from cold and warm regions, and many of these differences involved metabolites with a known role in cold acclimation. Since the transcriptomic component of cold acclimation is well studied, we analyzed the expression profiles of 251 previously described cold-acclimation genes (Figure 3—source data 1) in eight accessions that were representative in terms of their growth and metabolome profiles (Figure 3—figure supplement 1). The selected genes are described in the literature as being differentially expressed upon exposure to cold, and their expression is under control of at least one of the known transcription factors regulating cold acclimation: CBF1, CBF2, CBF3, HSFC1 (Park et al., 2015), or ZAT12 (Vogel et al., 2005). In our experiment, expression of these genes is likewise more affected by temperature than expected by chance (Figure 3; $χ^{2}$ -test: $p - v a l u e < 0.001$ ) and separates the two accessions from the coldest region (northern Sweden) from the rest in the 16°C treatment, and the three accessions from the warmest regions (Spain and Azerbaijan) from the rest in the 6°C treatment. Expression of different subsets of the selected cold-acclimation genes shows clear correlations with winter temperature of origin (Figure 3—figure supplement 2, Figure 3—figure supplement 3). In particular, the genes that were previously found to be upregulated upon cold exposure showed higher expression in accessions from cold climates (Figure 3—figure supplement 4). Since the expression of these cold-acclimation genes has been linked to the strength of cold acclimation in previous experiments (Park et al., 2015; Vogel et al., 2005), these accessions likely differ in their ability to cope with freezing temperatures upon cold treatment.

Figure 3 with 5 supplements see all

Download asset Open asset

Expression of 251 previously described cold-acclimation genes.

Expression is shown as the gene-wide z-scores of the normalized counts. The z-scores allow for grouping genes with a similar expression behavior over the different accessions in both temperatures. The top bar indicates winter temperature (°C) for each accession’s origin. Both dendrograms along y-axis and x-axis, respectively, show hierarchical clustering of genes, and of accessions in both temperatures.

Figure 3—source data 1 Cold-acclimation genes and their expression cluster membership as shown in Figure 3.: https://cdn.elifesciences.org/articles/77913/elife-77913-fig3-data1-v1.csv
Download elife-77913-fig3-data1-v1.csv

Growth is polygenic and shows signs of local adaptation

We used genome-wide association to investigate the genetic architecture underlying variation for the different growth parameters. As expected, these traits appear to be highly polygenic, and there were no genome-wide significant associations (Figure 4—figure supplement 1). The strongest association was found for overall growth rate at 16°C (Figure 4A). Inflated significance levels after correcting for population structure are consistent with what we would expect from a polygenic trait (Figure 4B) and were also observed for the other traits, except for growth rate at 6°C (Figure 4—figure supplement 1). Plausible candidates within 10 kb of the most significant SNP (chr5: 23334281; $- \log_{1} 0 (p - v a l u e) = 6.85$ ) include CIPK21 and MYB36. CIPK21 encodes a CBL-interacting protein kinase that is upregulated in cold conditions and is involved in the salt and osmotic stress response (Pandey et al., 2015). MYB36 is a key regulator of root endodermal differentiation (Liberman et al., 2015). Slightly more distant, 22 kb away, is $C O L 5$ , encoding a transcription factor that is part of the gene network that is regulated by AN3, a regulator of cell proliferation in leaf growth (Vercruyssen et al., 2014).

Figure 4 with 1 supplement see all

Download asset Open asset

Genome-wide association study (GWAS) results for the growth rate at 16°C.

(A) Manhattan plot showing the significance of the association between the phenotype and each of the tested SNPs ( $M A F > 10 %$ ). The Bonferroni-corrected threshold is shown with a dashed red line. (B) QQ-plot showing the relation between observed and expected $- \log 10 (p - v a l u e)$ distributions. Red line shows the observed relationship. The gray line and band show the expected relationship under the null hypothesis of no differentiation between both distributions.

To test for potential polygenic adaptation, we compared the phenotypic divergence to the expected neutral genome-wide genetic divergence. This can be done using a $Q_{S T} - F_{S T}$ test (Prout and Barker, 1993; Whitlock, 2008; Spitze, 1993); however, this test is not well suited for species with complex population structure, and so we used a variation that was designed to detect adaptive differentiation for traits measured in structured GWAS panels (Josephs et al., 2019). Instead of looking at divergences between predefined populations, this method uses principal components (PCs) of the genetic relatedness matrix as axes of potential adaptive differentiation. Adaptive differentiation is then detected as a correlation between the focal phenotype and any of these relatedness PCs that is significantly different than expected under neutrality.

Adaptive differentiation was detected for initial size and for growth rate at 16°C and its temperature response. These traits show adaptive differentiation along different genetic axes (Figure 5). Initial size shows significant adaptive differentiation along PC6 ( $p - v a l u e < 0.05$ ), whereas growth rate at 16°C and its temperature response showed significant adaptive differentiation along PC5 ( $p - v a l u e s < 0.05$ ). Adaptive differentiation was not significant along the other axes of genetic differentiation (PC1–4, PC7–10). The adaptive differentiation for initial size along PC6 seems to stem from higher initial sizes in Swedish accessions compared to central European accessions. The adaptive differentiation along PC5 seems to be driven by the lower growth rate temperature responses in Asian and northern Swedish accessions in contrast to higher growth rates in a subset of southern Swedish accessions. The accessions in our set that come from northern Sweden and Asia hail from the coldest climates. Thus, these results suggest adaptive differentiation driven by adaptation to cold winters. Given the seemingly strong influence from the Asian accessions, we repeated the analysis without them. Also in this analysis we detected significant adaptive differentiation ( $p - v a l u e s < 0.05$ ) for initial size, growth rate at 16°C, and its temperature response (Figure 5—figure supplement 1).

Figure 5 with 1 supplement see all

Download asset Open asset

Adaptive differentiation of initial size, growth rate at 16°C, and the temperature response of growth rate along different axes of genetic differentiation.

Plots represent the phenotypes and axes of genetic differentiation for which we detected significant adaptive differentiation; initial size and PC6 (A), growth rate in 16ºC and PC5 (B), and the growthrate's temperature response and PC5 (C). Accessions are colored according to their respective admixture groups, as specified in 1001 Genomes Consortium, 2016. The gray ribbon represents the expected correlation between phenotype and axis of genetic differentiation under neutrality with a 90% confidence interval. The neutral expectation is based on axes of genetic differentiation within populations (see ‘Materials and methods’ and Josephs et al., 2019 for further details). The blue line represents the observed correlation between phenotype and axis of genetic differentiation. Percentages refer to the genetic variation explained by the respective principal component.

Discussion

This study explores natural variation of rosette growth in nonfreezing temperatures. We detect genetic variation for the different growth parameters, and environmental correlations that suggest local adaptation. GWAS analyses reveal, not surprisingly, a polygenic trait architecture. We speculated that the slower growth measured in accessions from colder climates reflects relocation of resources from growth towards cold acclimation. Both metabolome and gene expression data are consistent with accessions from colder climates preparing for a harsh winter. In our temperature experiment, we see that the growth of northern lines is affected less than southern lines by switching from 16°C to 6°C.

Our conclusion that slower growth is likely adaptive in populations facing fiercer winters is in line with recent results of Wieters et al., 2021, who concluded that the reduced growth in northern lines was adaptive and not a consequence of an accumulation of deleterious mutations at the species border. If slower growth were indeed a consequence of accumulated deleterious mutations, we would expect to see slower growth also during the initial seedling establishment, which we measured here as the initial size. On the contrary, we saw a fast seedling establishment for accessions from colder regions. We speculate that the fast seedling establishment is a potential adaptation for short growth seasons, which often coincide with colder climates (high latitude or high altitude). This fast seedling establishment seems to be partly supported by larger seeds. These larger seeds may provide more nutrients to initiate faster seedling establishment, while this is of less importance in warmer climates with longer growth seasons. Further work is needed to disentangle initial growth from seed size effects and confirm that there is a causal relationship between seed size and fast seedling establishment, whether this is due to seed nutrient storage, and whether it is adaptive.

The adaptation of growth to local climates is likely to be influenced by a trade-off with cold acclimation. General growth-survival trade-offs have long been observed and are described in general ecological strategy schemes such as Grime’s C-S-R triangle (Grime, 1979) and the leaf–height–seed scheme (Westoby, 1998). Specific trade-offs between growth and cold/frost survival were observed for wheat (Hayes and Aamodt, 1927; Quisenberry, 1931), alfalfa (Castonguay et al., 2006), Dactylis glomerata (Bristiel et al., 2018), and multiple tree species (Koehler et al., 2012; Loehle, 1998; Molina-Montenegro et al., 2012; Savage and Cavender-Bares, 2013). Here, we observed higher expression of genes involved in cold acclimation in accessions from colder regions. This is clearest at 6°C, but is also happening at 16°C, suggesting that activation of cold acclimation is stronger in accessions from cold climates and may be triggered more easily. More generally, it establishes that geographic differentiation exists in the expression of known cold-response genes. Accessions from warm climates may instead activate cold acclimation as a stress response rather than as a preventive measure. Even though this is based on a limited set of 8 accessions, metabolome measurements in all 249 accessions lead to the same conclusion. Metabolites involved in cold acclimation such as raffinose, sucrose, and proline were found in higher concentrations in accessions from colder climates (Weiszmann et al., 2020). We believe that accessions from colder environments are relocating more energy and resources from growth towards preparations for upcoming freezing temperatures, which is a clear example of active plasticity (Forsman, 2015; van Kleunen and Fischer, 2005). Even though we have no direct survival measurements, we speculate that this results in stronger cold acclimation and consequently increased freezing tolerance in the accessions from colder regions. Indeed, accessions originating from colder environments show increased freezing tolerance upon cold acclimation (Zhen et al., 2011; Zuther et al., 2012; Hannah et al., 2006; Horton et al., 2016). This fits with observations of northern and colder regions favoring slower growing, more stress-tolerant plants (Vasseur et al., 2018; Estarague et al., 2022). Also, biogeographic studies in A. thaliana found that winter temperatures are a major determinant of suitable habitats for this species (Hoffmann, 2002; Yim et al., 2022), and reciprocal transplant experiments detected an important role for freezing tolerance in fitness variation in northern sites (Ågren and Schemske, 2012). The high variability we observed in our data does, however, show that there is more at play than selection for cold resistance alone. What these factors are we can only speculate about. Phenotypes are shaped by a mixture of neutral and adaptive processes, with a plethora of trade-offs between traits. Investigating phenotypes at different organismal scales in specific and realistic environments will further elucidate how phenotypes and, ultimately, life history strategies are shaped. We speculate that the reduced growth plasticity observed for accessions from colder climates is due to the stronger growth reduction at 16°C in these accessions compared to accessions from warmer climates. They may well be anticipating winter, whereas accessions from warmer climates do not, and hence show a stronger difference between 16 and 6°C. Although both metabolite and gene expression data suggest an involvement of cold acclimation, we note that it is impossible to rule out that accessions from colder climates grow slower due to reduced resource efficiency as they become larger or, for example, increase leaf thickness (Adams et al., 2016). Further work is needed to understand the mechanism underlying the growth response observed here.

There is strong evidence from QTL mapping that genetic variation in the CBF2 gene is one of the drivers for adaptation to freezing stress (Oakley et al., 2014; Gehan et al., 2015). Here, we looked at growth phenotypes and did not detect associations with the CBF loci. In the transcriptome analysis, we did pick up a role for CBF and other known cold-acclimation genes. The most significant locus detected in our GWAS analysis (for growth rate at 16°C) lies in the vicinity of COL5, a gene that is part of a leaf growth regulatory network (Vercruyssen et al., 2014) and whose expression is induced by both cold treatment and CBF1, CBF2, and CBF3 overexpression (Park et al., 2015). It is however unclear what its exact regulatory role in growth in cold conditions might be. In summary, we detected adaptive differentiation for growth between accessions from warm and cold climates. Our transcriptome data and previous metabolome data suggest that resources are relocated from growth to cold acclimation in accessions from colder regions. This allows these accessions to be fully prepared for the coming of winter.

Materials and methods

Plant growth and phenotyping

Request a detailed protocol

Seeds of 249 natural accessions (Figure 1—source data 1) of A. thaliana described in the 1001 genomes project (1001 Genomes Consortium, 2016) were sown on sieved (6 mm) substrate (Einheitserde ED63). Pots were filled with 71.5 g ± 1.5 g of soil to ensure homogenous packing. The prepared pots were all covered with blue mats (Junker et al., 2014) to enable a robust performance of the high-throughput image analysis algorithm. Seeds were stratified (4 days at 4°C in darkness) after which they germinated and left to grow for 2 weeks at 21°C (relative humidity: 55%; light intensity: $160 μ m o l m^{- 2} s^{- 1}$ ; 14 hr light). The temperature treatments were started by transferring the seedlings to either 6 or 16°C. To simulate natural conditions, temperatures fluctuated diurnally between 16–21°C, 0.5–6°C, and 8–16°C for the 21°C initial growth conditions and the 6 and 16°C treatments, respectively (Figure 1—figure supplement 1). Light intensity was kept constant at $160 μ m o l m^{- 2} s^{- 1}$ throughout the experiment. Relative humidity was set at 55% but in colder temperatures it rose uncontrollably to maximum 95%. Daylength was 9 hr during the 16 and 6°C treatments. Each temperature treatment was repeated in three independent experiments. Five replicate plants were grown for every genotype per experiment. Plants were randomly distributed across the growth chamber with an independent randomization pattern for each experiment. During the temperature treatments (14–35 days after stratification), plants were photographed twice a day (1 hr after/before lights switched on/off), using an RGB camera (IDS uEye UI-548xRE-C; 5MP) mounted to a robotic arm. Rosette areas were extracted from the plant images using Lemnatec OS (LemnaTec GmbH, Aachen, Germany) software. Plant growth profiles were visually inspected, and datapoints with smaller rosette areas than earlier time points (negative growth) were discarded from further analyses. At 35 days after stratification, whole rosettes were harvested, immediately frozen in liquid nitrogen, and stored at –80°C until further analysis.

Nonlinear modeling

Request a detailed protocol

Nonlinear modeling was used to describe plant growth in a minimum number of parameters. In a first step, we constructed a simple nonlinear model with plant size being explained by either the exponential (Equation 1a, Equation 1b) or the power-law function (Equation 2a, Equation 2b), with individual plant as a random effect for each of the model parameters; M₀, $r$ , and $β$ . With $β$ being only present in the power-law model. Models were constructed using the nlsList and nlme functions from the nlme package (3.1.152; Pinheiro et al., 2021) for R (4.0.3; R Development Core Team, 2017). Exponential and power-law SelfStart functions were used from Paine et al., 2012. Based on Akaike information criterion and likelihood ratio test generated by the anova function (Table 1), we decided to use the power-law model for further analyses.

\frac{d M}{d t} = r M

M_{t} = M_{0} e^{r t}

\frac{d M}{d t} = r M^{β}

M_{t} = {(M_{0}^{1 - β} + r t (1 - β))}^{1 / (1 - β)}

Table 1

ANOVA table for the comparison between the exponential and power-law model with degrees of freedom (df), Akaike information criterion (AIC), Bayesian information criterion (BIC), and log-likelihood (logLik) for each model.

The likelihood ratio statistic (L.ratio) and p-value are given for the likelihood ratio test that was used to compare these models.

	Model	df	AIC	BIC	logLik	Test	L.ratio	p-Value
Exponential	1	6	–668551.6	–668488.1	334281.8		-	-
Power-law	2	10	–870605.0	–870499.1	435312.5	1 vs. 2	202.061	<0.001

In a second step, we constructed a model with fixed effects for the different power-law parameters. For initial size (M₀), we added accession as fixed effect. Temperature treatment only started from the initial time point onwards, and thus could not have an effect on the initial plant size. The growth rate, on the other hand, should be affected by temperature; therefore, we included accession, temperature, and their interaction as fixed effects for growth rate ( $r$ ). No fixed effects were added for $β$ . The idea here is that it is an adjustment factor for decreasing growth rates (when $β < 1$ ) with increasing plant sizes, which is general for plant growth, or at least for our data in this case. Individuals nested within experiment were added as random effects for each of the model parameters. The correlation structure intrinsic to measuring the same individuals over time was accounted for by adding the first-order continuous autoregressive correlation structure (corCAR1). The estimated fixed effects of this model were then used to obtain initial size estimates for each accession and growth rate estimates for each accession in both temperatures. These estimates were used for all further analyses apart from broad-sense heritability calculations (see below). For each accession, we calculated the growth rate response as the slope between the growth rate at 6°C and the growth rate at 16°C. The slope was obtained from linear regression with the lm function in R (4.0.3; R Development Core Team, 2017) using temperature as an ordered categorical variable (6°C < 16°C).

Climate correlations

Request a detailed protocol

The different phenotypes were correlated with each of the different (bio)climate variables downloaded from https://www.worldclim.org (Fick and Hijmans, 2017). Correlations were calculated as Pearson’s correlations using the cor function in R (4.0.3; R Development Core Team, 2017). Population structure may confound the correlation between phenotype and climate. Therefore, we included additional phenotype–climate correlations with correction for population structure (Figure 2—figure supplement 3). For the population–structure-corrected correlations, we used a mixed-effects model as implemented in the lmekin function from the coxme (2.2.16; Therneau, 2020) package with phenotype as dependent variable, climate variable as fixed effects, and the kinship matrix as random effect. The kinship matrix was based on the SNPs from the 1001 genomes consortium (1001 Genomes Consortium, 2016) and was calculated using ’mixmogam’ (https://github.com/bvilhjal/mixmogam; Vilhjalmsson, 2019) based on Kang et al., 2010. In this analysis, phenotype and climate variables were standardized, so that regression coefficients were comparable to correlation coefficients. Even though the strength and significance of the correlations weaken upon population structure correction, the growth parameters still demonstrate the same pattern, being most strongly correlated with winter temperatures.

Seed size correlations

Request a detailed protocol

We used the seeds produced by Kerdaffrec et al., 2016 and limited our measurements to the set of 123 Swedish accessions that overlapped with our growth dataset. After seed stratification for four days at 4°C in darkness, mother plants were grown for 8 weeks at 4°C under long-day conditions (16 hr light; 8 hr dark) to ensure proper vernalization. Temperature was raised to 21°C (light) and 16°C (dark) for flowering and seed ripening. Seeds were kept in darkness at 16°C and 30% relative humidity, from the harvest until seed size measurements. For each genotype, three replicates were pooled and about 200–300 seeds were sprinkled on 12 × 12 cm², transparent Petri dishes. Image acquisition was performed as described in Exposito-Alonso et al., 2018 by scanning dishes on a cluster of eight Epson V600 scanners. The resulting 1200 dpi .tiff images were analyzed with the ImageJ software (2.1.0/1.53c). Images were converted to eight-bit binary images and thresholded with the setAutoThreshold("Defaultdark”) command, and seed area was measured in mm² by running the Analyse Particles command (inclusion parameters: $s i z e = 0.04 - 0.25 c i r c u l a r i t y = 0.70 - 1.00$ ). All scripts used for image processing are available at https://github.com/vevel/seed_size; Kerdaffrec, 2022. The variance decomposition for initial size into variance explained by winter temperature and seed size was done with a random-effect model where initial size was explained by winter temperature and seed size, using the lmer and VarCorr functions from the lme4 package (1.1.27.1; Bates et al., 2015) in R (4.0.3; R Development Core Team, 2017). The seed size-corrected correlation between initial size and winter temperature was estimated with the lme function from the nlme package (3.1.152; Pinheiro et al., 2021) in R (4.0.3; R Development Core Team, 2017), and the correction for seed size was done by including seed size as a random effect.

Transcriptome profiling

Request a detailed protocol

35 days after stratification, rosette tissue of all plants were harvested and flash frozen in liquid nitrogen. Random samples from each replicate experiment for both temperatures were taken for eight accessions to profile the transcriptome with RNA-sequencing. The eight accessions were selected to represent the climatic variation in the full panel (Figure 3—figure supplement 1, Figure 1—source data 1). Total RNA was extracted using the KingFisher Duo Prime System (Thermo Fisher Scientific) together with a high-performance RNA bead isolation kit (Molecular Biology Service, VBC Core Facilities, Vienna). To determine the quantity of RNA, we used Fluorometer Qubit 4 (Invitrogen) and Qubit RNA BR Kit (Invitrogen). For each sample, 1 µg of total RNA was treated with the poly(A) RNA Selection Kit (Lexogen) and eluted in 12 µl of Nuclease-Free Water. Libraries were prepared according to the manufacturer’s protocol in NEBNext Ultra II Directional RNA Library Prep Kit (New England Biolabs) and individually indexed with NEBNext Multiplex Oligos for Illumina (New England Biolabs). The quantity and quality of each amplified library were analyzed by using Fragment Analyzer (Agilent) and HS NGS Fragment Kit (Agilent). Libraries were sequenced with an Illumina HiSeq2500 in paired-end mode with read length of 125 bp. Sequencing was performed by the Next Generation Sequencing Facility at Vienna BioCenter Core Facilities (VBCF), member of the Vienna BioCenter (VBC), Austria. Samples were distributed over four independent libraries. This was due to failed samples that needed replacement. Detailed info on which samples belong to which library is listed in SRA (https://www.ncbi.nlm.nih.gov/sra/PRJNA807069). Gene expression was quantified by using quasi-mapping in salmon (1.2.1; Patro et al., 2017). The salmon indices were built separately for each accession as we incorporated the SNP variation from the 1001 Genomes Consortium, 2016 into the reference transcriptome. The heatmap was built using pheatmap (1.0.12; Kolde, 2019) in R (4.0.3; R Development Core Team, 2017). The clustering was done by complete clustering on Euclidean distances for both rows and columns. The gene clusters were defined by cutting the dendrogram in seven groups. With a -test, we tested for overrepresentation of a temperature effect on the expression of the 251 selected cold-acclimation genes ( $d f = 1$ ) compared to the remaining 18,784 background genes. We used the chisq.test function in R (4.0.3; R Development Core Team, 2017). Differential expression analysis was conducted with the DESeq2 package (1.30.0; Love et al., 2014) in R (4.0.3; R Development Core Team, 2017). A full model was used, with expression depending on $r e p l i c a t e + a c c e s s i o n + t e m p e r a t u r e + r e p l i c a t e : t e m p e r a t u r e + a c c e s s i o n : t e m p e r a t u r e$ , after which significance of each model coefficient was defined with a negative binomial Wald test. Differential expression for each accession was then extracted by specifying the respective contrasts using the lfcShrink function in DESeq2 with the adaptive shrinkage estimator (Stephens, 2017). Genes were considered differentially expressed when the adjusted p-value was <0.05.

Metabolome profiling

Request a detailed protocol

Besides transcriptome profiling of eight accessions, we also conducted metabolome profiling on all 249 accessions. Samples for metabolome measurements were taken from the same experiments described in this study and, just like the transcriptome samples, were taken 35 days after stratification. Results and detailed methodology are described in Weiszmann et al., 2020.

Broad-sense heritabilities

Request a detailed protocol

Broad-sense heritabilities ( $H^{2}$ ) were calculated as the ratio between phenotypic variation explained by genotype ( $V_{g}$ ) and the total phenotypic variation ( $V_{p}$ ), which is the sum of $V_{g}$ and phenotypic variation explained by environment ( $V_{e}$ ). These variances were obtained from a mixed model by including accession as a random effect (estimate for $V_{g}$ ). Because our accession estimates were corrected for experiment effects, we removed the variance explained by experiment by including experiment as a fixed effect. In order to estimate the variance within each accession ( $V_{e}$ ), the dependent variables in this model were the growth parameter estimates for each individual plant, in contrast to the estimates for each accession that were used in all other analyses. These individual plant estimates were obtained from the same model, but took the random effects into account. The variance explained by accession was then taken as an estimate for $V_{g}$ , the residual variance was taken as $V_{e}$ . For initial size, we calculated heritability over all experiments, and growth rate heritabilities were calculated for each temperature independently. The mixed model was constructed with the lmer function in the lme4 package (1.1.27.1; Bates et al., 2015) in R (4.0.3; R Development Core Team, 2017).

Genome-wide association mapping

Request a detailed protocol

Genome-wide association mapping was done for each of the growth parameters in both temperatures and also the temperature response for the growth rate. We used a mixed model with phenotype as dependent variable, genotype as fixed effect, and genetic relatedness as random factor. Nonimputed SNPs obtained from the 1001 Genomes Consortium, 2016 were used as genotypes. This model was run in GEMMA (0.98.3; Zhou and Stephens, 2012), with kinship matrix calculated as the centered relatedness matrix, as implemented in GEMMA.

Testing for adaptive differentiation

Request a detailed protocol

Adaptive differentiation was tested with the method described by Josephs et al., 2019 and the accompanying quaint package (0.0.0.9; https://github.com/emjosephs/quaint; Josephs, 2020) in R (4.0.3; R Development Core Team, 2017). The kinship matrix was calculated using the make_k function in the quaint package. Genetic PCs were then calculated from the eigen decomposition of the kinship matrix. Adaptive differentiation of each phenotype along the first 10 PCs was tested with the calcQpc function in the quaint package. PCs 11–248 were used to build the expected phenotypic differentiation under neutrality.

Supplementary information

Request a detailed protocol

Scripts can be found at https://github.com/picla/growth_16C_6C/; Clauw, 2022. Scripts for seedsize analysis can be found at https://github.com/vevel/seed_size; Kerdaffrec, 2022.

All RNA-sequencing were uploaded to SRA under http://www.ncbi.nlm.nih.gov/bioproject/807069. All generated phenotyping data are filed under 10.5281/zenodo.6076948.

Data availability

Scripts can be found at https://github.com/picla/growth_16C_6C/, copy archived at swh:1:rev:f32793d07a068a6c49c876536f684ffebdcd9e6b. Scripts for seedsize analysis can be found in https://github.com/vevel/seed_size, copy archived at swh:1:rev:748523031aa64601720328ab49e02009a6fb70da. All RNA-sequencing were uploaded to SRA under http://www.ncbi.nlm.nih.gov/bioproject/807069. All generated phenotyping data are filed under https://doi.org/10.5281/zenodo.6076948.

The following data sets were generated

1. Clauw P
(2022) NCBI BioProject
ID PRJNA807069. Arabidopsis transcriptome in 16C and 6C.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA807069

References

1. 1001 Genomes Consortium
(2016) 1,135 genomes reveal the global pattern of polymorphism in Arabidopsis thaliana
Cell 166:481–491.

https://doi.org/10.1016/j.cell.2016.05.063
- PubMed
- Google Scholar
(2016) Habitat temperature and precipitation of Arabidopsis thaliana ecotypes determine the response of foliar vasculature, photosynthesis, and transpiration to growth temperature
Frontiers in Plant Science 7:1026.

https://doi.org/10.3389/fpls.2016.01026
- PubMed
- Google Scholar
1. Ågren J
2. Schemske DW
(2012) Reciprocal transplants demonstrate strong adaptive differentiation of the model organism Arabidopsis thaliana in its native range
The New Phytologist 194:1112–1122.

https://doi.org/10.1111/j.1469-8137.2012.04112.x
- PubMed
- Google Scholar
(2005) Genetic and molecular analyses of natural variation indicate CBF2 as a candidate gene for underlying a freezing tolerance quantitative trait locus in Arabidopsis
Plant Physiology 139:1304–1312.

https://doi.org/10.1104/pp.105.068510
- PubMed
- Google Scholar
1. Artus NN
2. Uemura M
3. Steponkus PL
4. Gilmour SJ
5. Lin C
6. Thomashow MF
(1996) Constitutive expression of the cold-regulated Arabidopsis thaliana COR15a gene affects both chloroplast and protoplast freezing tolerance
PNAS 93:13404–13409.

https://doi.org/10.1073/pnas.93.23.13404
- PubMed
- Google Scholar
(2015) Genome-wide association mapping of growth dynamics detects time-specific and general quantitative trait loci
Journal of Experimental Botany 66:5567–5580.

https://doi.org/10.1093/jxb/erv176
- PubMed
- Google Scholar
1. Bates D
2. Mächler M
3. Bolker B
4. Walker S
(2015) Fitting linear Mixed-Effects models using lme4
Journal of Statistical Software 67:1–48.

https://doi.org/10.18637/jss.v067.i01
- Google Scholar
(2011) The low temperature response pathways for cold acclimation and vernalization are independent
Plant, Cell & Environment 34:1737–1748.

https://doi.org/10.1111/j.1365-3040.2011.02370.x
- PubMed
- Google Scholar
1. Bristiel P
2. Gillespie L
3. Østrem L
4. Balachowski J
5. Violle C
6. Volaire F
7. Field K
(2018) Experimental evaluation of the robustness of the growth–stress tolerance trade‐off within the perennial grass Dactylis glomerata
Functional Ecology 32:1944–1958.

https://doi.org/10.1111/1365-2435.13112
- Google Scholar
(2006)
In Advances in Agronomy

203–265, Alfalfa winter hardiness: A research retrospective and integrated perspective, In Advances in Agronomy, Academic Press.
- Google Scholar
1. Claeys H
2. Inzé D
(2013) The agony of choice: how plants balance growth and survival under water-limiting conditions
Plant Physiology 162:1768–1779.

https://doi.org/10.1104/pp.113.220921
- PubMed
- Google Scholar
Software
1. Clauw P
(2022) growth_16C_6C, version swh:1:rev:f32793d07a068a6c49c876536f684ffebdcd9e6b
GitHub.

https://github.com/picla/growth_16C_6C
1. Estarague A
2. Vasseur F
3. Sartori K
4. Bastias CC
5. Cornet D
6. Rouan L
7. Beurier G
8. Exposito-Alonso M
9. Herbette S
10. Bresson J
11. Vile D
12. Violle C
(2022) Into the range: a latitudinal gradient or a center-margins differentiation of ecological strategies in Arabidopsis thaliana?
Annals of Botany 129:343–356.

https://doi.org/10.1093/aob/mcab149
- PubMed
- Google Scholar
1. Exposito-Alonso M
2. Becker C
3. Schuenemann VJ
4. Reiter E
5. Setzer C
6. Slovak R
7. Brachi B
8. Hagmann J
9. Grimm DG
10. Chen J
11. Busch W
12. Bergelson J
13. Ness RW
14. Krause J
15. Burbano HA
16. Weigel D
(2018) The rate and potential relevance of new mutations in a colonizing plant lineage
PLOS Genetics 14:e1007155.

https://doi.org/10.1371/journal.pgen.1007155
- PubMed
- Google Scholar
1. Fick SE
2. Hijmans RJ
(2017) WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas
International Journal of Climatology 37:4302–4315.

https://doi.org/10.1002/joc.5086
- Google Scholar
1. Forsman A
(2015) Rethinking phenotypic plasticity and its consequences for individuals, populations and species
Heredity 115:276–284.

https://doi.org/10.1038/hdy.2014.92
- PubMed
- Google Scholar
1. Fritz MA
2. Rosa S
3. Sicard A
(2018) Mechanisms underlying the environmentally induced plasticity of leaf morphology
Frontiers in Genetics 9:478.

https://doi.org/10.3389/fgene.2018.00478
- PubMed
- Google Scholar
1. Gehan MA
2. Park S
3. Gilmour SJ
4. An C
5. Lee CM
6. Thomashow MF
(2015) Natural variation in the C-repeat binding factor cold response pathway correlates with local adaptation of Arabidopsis ecotypes
The Plant Journal 84:682–693.

https://doi.org/10.1111/tpj.13027
- PubMed
- Google Scholar
(2000) Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation
Plant Physiology 124:1854–1865.

https://doi.org/10.1104/pp.124.4.1854
- PubMed
- Google Scholar
1. Grime JP
(1979) Primary strategies in plants
Transactions of the Botanical Society of Edinburgh 43:151–160.

https://doi.org/10.1080/03746607908685348
- Google Scholar
1. Hannah MA
2. Wiese D
3. Freund S
4. Fiehn O
5. Heyer AG
6. Hincha DK
(2006) Natural genetic variation of freezing tolerance in Arabidopsis
Plant Physiology 142:98–112.

https://doi.org/10.1104/pp.106.081141
- PubMed
- Google Scholar
1. Hayes HK
2. Aamodt OS
(1927)
Inheritance of winter hardiness and growth habit in crosses of marquis with minhardi and minturki wheats

Journal of Agricultural Research 35:223–236.
- Google Scholar
1. Hepworth J
2. Antoniou-Kourounioti RL
3. Bloomer RH
4. Selga C
5. Berggren K
6. Cox D
7. Collier Harris BR
8. Irwin JA
9. Holm S
10. Säll T
11. Howard M
12. Dean C
(2018) Absence of warmth permits epigenetic memory of winter in Arabidopsis
Nature Communications 9:639.

https://doi.org/10.1038/s41467-018-03065-7
- PubMed
- Google Scholar
1. Hoffmann MH
(2002) Biogeography of Arabidopsis thaliana (l.) heynh. (brassicaceae)
Journal of Biogeography 29:125–134.

https://doi.org/10.1046/j.1365-2699.2002.00647.x
- Google Scholar
(2016) The genetic architecture of freezing tolerance varies across the range of Arabidopsis thaliana
Plant, Cell & Environment 39:2570–2579.

https://doi.org/10.1111/pce.12812
- PubMed
- Google Scholar
1. Hughes MA
2. Dunn MA
(1996) The molecular biology of plant acclimation to low temperature
Journal of Experimental Botany 47:291–305.

https://doi.org/10.1093/jxb/47.3.291
- Google Scholar
(2019) Detecting adaptive differentiation in structured populations with genomic data and common gardens
Genetics 211:989–1004.

https://doi.org/10.1534/genetics.118.301786
- PubMed
- Google Scholar
Software
1. Josephs E
(2020) quaint, version swh:1:rev:748dd884e3b9cbb8a5a42932ec50b6c7f0cef979
GitHub.

https://github.com/emjosephs/quaint
(2014) Optimizing experimental procedures for quantitative evaluation of crop plant performance in high throughput phenotyping systems
Frontiers in Plant Science 5:770.

https://doi.org/10.3389/fpls.2014.00770
- PubMed
- Google Scholar
1. Kang HM
2. Sul JH
3. Service SK
4. Zaitlen NA
5. Kong SY
6. Freimer NB
7. Sabatti C
8. Eskin E
(2010) Variance component model to account for sample structure in genome-wide association studies
Nature Genetics 42:348–354.

https://doi.org/10.1038/ng.548
- PubMed
- Google Scholar
1. Kerdaffrec E
2. Filiault DL
3. Korte A
4. Sasaki E
5. Nizhynska V
6. Seren Ü
7. Nordborg M
(2016) Multiple alleles at a single locus control seed dormancy in Swedish Arabidopsis
eLife 5:e22502.

https://doi.org/10.7554/eLife.22502
- PubMed
- Google Scholar
Software
1. Kerdaffrec E
(2022) vevel/seed_size, version swh:1:rev:748523031aa64601720328ab49e02009a6fb70da
GitHub.

https://github.com/vevel/seed_size
(2012) Evidence for a freezing tolerance-growth rate trade-off in the live oaks (Quercus series Virentes) across the tropical-temperate divide
The New Phytologist 193:730–744.

https://doi.org/10.1111/j.1469-8137.2011.03992.x
- PubMed
- Google Scholar
Software
1. Kolde R
(2019) pheatmap: Pretty heatmaps, version 1.0.12
CRAN.

https://CRAN.R-project.org/package=pheatmap
1. Lande R
(1982) A quantitative genetic theory of life history evolution
Ecology 63:607–615.

https://doi.org/10.2307/1936778
- Google Scholar
(2015) MYB36 regulates the transition from proliferation to differentiation in the Arabidopsis root
PNAS 112:12099–12104.

https://doi.org/10.1073/pnas.1515576112
- PubMed
- Google Scholar
1. Loehle C
(1998) Height growth rate tradeoffs determine northern and southern range limits for trees
Journal of Biogeography 25:735–742.

https://doi.org/10.1046/j.1365-2699.1998.2540735.x
- Google Scholar
1. Love MI
2. Huber W
3. Anders S
(2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
Genome Biology 15:550.

https://doi.org/10.1186/s13059-014-0550-8
- PubMed
- Google Scholar
1. Marchadier E
2. Hanemian M
3. Tisné S
4. Bach L
5. Bazakos C
6. Gilbault E
7. Haddadi P
8. Virlouvet L
9. Loudet O
(2019) The complex genetic architecture of shoot growth natural variation in Arabidopsis thaliana
PLOS Genetics 15:e1007954.

https://doi.org/10.1371/journal.pgen.1007954
- PubMed
- Google Scholar
(2020) Functional variants of DOG1 control seed chilling responses and variation in seasonal life-history strategies in Arabidopsis thaliana
PNAS 117:2526–2534.

https://doi.org/10.1073/pnas.1912451117
- PubMed
- Google Scholar
(2012) The trade-off between cold resistance and growth determines the Nothofagus pumilio treeline
Plant Ecology 213:133–142.

https://doi.org/10.1007/s11258-011-9964-5
- Google Scholar
(1999) Antisense suppression of proline degradation improves tolerance to freezing and salinity in Arabidopsis thaliana
FEBS Letters 461:205–210.

https://doi.org/10.1016/s0014-5793(99)01451-9
- PubMed
- Google Scholar
(2014) QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs
Molecular Ecology 23:4304–4315.

https://doi.org/10.1111/mec.12862
- PubMed
- Google Scholar
1. Paine CET
2. Marthews TR
3. Vogt DR
4. Purves D
5. Rees M
6. Hector A
7. Turnbull LA
(2012) How to fit nonlinear plant growth models and calculate growth rates: an update for ecologists
Methods in Ecology and Evolution 3:245–256.

https://doi.org/10.1111/j.2041-210X.2011.00155.x
- Google Scholar
1. Pandey GK
2. Kanwar P
3. Singh A
4. Steinhorst L
5. Pandey A
6. Yadav AK
7. Tokas I
8. Sanyal SK
9. Kim BG
10. Lee SC
11. Cheong YH
12. Kudla J
13. Luan S
(2015) Calcineurin b-like protein-interacting protein kinase cipk21 regulates osmotic and salt stress responses in arabidopsis
Plant Physiology 169:780–792.

https://doi.org/10.1104/pp.15.00623
- PubMed
- Google Scholar
1. Park S
2. Lee CM
3. Doherty CJ
4. Gilmour SJ
5. Kim Y
6. Thomashow MF
(2015) Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network
The Plant Journal 82:193–207.

https://doi.org/10.1111/tpj.12796
- PubMed
- Google Scholar
(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.

https://doi.org/10.1371/journal.pone.0207723
- PubMed
- Google Scholar
1. Patro R
2. Duggal G
3. Love MI
4. Irizarry RA
5. Kingsford C
(2017) Salmon provides fast and bias-aware quantification of transcript expression
Nature Methods 14:417–419.

https://doi.org/10.1038/nmeth.4197
- PubMed
- Google Scholar
Software
1. Pinheiro J
2. Bates D
3. DebRoy S
4. Sarkar D
5. Team RC
(2021) Linear and nonlinear mixed effects models, version 3.1-158
Nlme.

https://svn.r-project.org/R-packages/trunk/nlme/
1. Prout T
2. Barker JS
(1993) F statistics in Drosophila buzzatii: selection, population size and inbreeding
Genetics 134:369–375.

https://doi.org/10.1093/genetics/134.1.369
- PubMed
- Google Scholar
Book
1. Quisenberry KS
(1931)
Inheritance of Winter Hardiness, Growth Habit, and Stem-Rust Reaction in Crosses Between Minhardi Winter and H-44 Spring Wheats

U.S. Department of Agriculture.
- Google Scholar
Software
1. R Development Core Team
(2017) R: A Language and Environment for Statistical Computing
R Foundation for Statistical Computing, Vienna, Austria.

https://www.r-project.org/index.html
1. Sartori K
2. Vasseur F
3. Violle C
4. Baron E
5. Gerard M
6. Rowe N
7. Ayala-Garay O
8. Christophe A
9. Jalón Lg
10. Masclef D
11. Harscouet E
12. Granado MDR
13. Chassagneux A
14. Kazakou E
15. Vile D
(2019) Leaf economics and slow-fast adaptation across the geographic range of Arabidopsis thaliana
Scientific Reports 9:10758.

https://doi.org/10.1038/s41598-019-46878-2
- PubMed
- Google Scholar
1. Savage JA
2. Cavender-Bares J
(2013) Phenological cues drive an apparent trade-off between freezing tolerance and growth in the family Salicaceae
Ecology 94:1708–1717.

https://doi.org/10.1890/12-1779.1
- PubMed
- Google Scholar
1. Simpson GG
2. Dean C
(2002) Arabidopsis, the Rosetta stone of flowering time?
Science 296:285–289.

https://doi.org/10.1126/science.296.5566.285
- PubMed
- Google Scholar
1. Spitze K
(1993) Population structure in Daphnia obtusa: quantitative genetic and allozymic variation
Genetics 135:367–374.

https://doi.org/10.1093/genetics/135.2.367
- PubMed
- Google Scholar
Book
1. Stearns SC
(1992)
The Evolution of Life Histories

Oxford University Press.
- Google Scholar
1. Stephens M
(2017) False discovery rates: a new deal
Biostatistics 18:275–294.

https://doi.org/10.1093/biostatistics/kxw041
- PubMed
- Google Scholar
(1998) Mode of action of the COR15a gene on the freezing tolerance of Arabidopsis thaliana
PNAS 95:14570–14575.

https://doi.org/10.1073/pnas.95.24.14570
- PubMed
- Google Scholar
1. Taji T
2. Ohsumi C
3. Iuchi S
4. Seki M
5. Kasuga M
6. Kobayashi M
7. Yamaguchi-Shinozaki K
8. Shinozaki K
(2002) Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana
The Plant Journal 29:417–426.

https://doi.org/10.1046/j.0960-7412.2001.01227.x
- PubMed
- Google Scholar
(2019) Linking genes with ecological strategies in Arabidopsis thaliana
Journal of Experimental Botany 70:1141–1151.

https://doi.org/10.1093/jxb/ery447
- PubMed
- Google Scholar
Software
1. Therneau TM
(2020) Mixed effects cox models, version 2.2-17
Coxme.

https://CRAN.R-project.org/package=coxme
1. Thomashow MF
(1999) Plant cold acclimation: Freezing tolerance genes and regulatory mechanisms
Annual Review of Plant Physiology and Plant Molecular Biology 50:571–599.

https://doi.org/10.1146/annurev.arplant.50.1.571
- PubMed
- Google Scholar
1. van Kleunen M
2. Fischer M
(2005) Constraints on the evolution of adaptive phenotypic plasticity in plants
The New Phytologist 166:49–60.

https://doi.org/10.1111/j.1469-8137.2004.01296.x
- PubMed
- Google Scholar
1. Vasseur F
2. Sartori K
3. Baron E
4. Fort F
5. Kazakou E
6. Segrestin J
7. Garnier E
8. Vile D
9. Violle C
(2018) Climate as a driver of adaptive variations in ecological strategies in Arabidopsis thaliana
Annals of Botany 122:935–945.

https://doi.org/10.1093/aob/mcy165
- PubMed
- Google Scholar
1. Vercruyssen L
2. Verkest A
3. Gonzalez N
4. Heyndrickx KS
5. Eeckhout D
6. Han S-K
7. Jégu T
8. Archacki R
9. Van Leene J
10. Andriankaja M
11. De Bodt S
12. Abeel T
13. Coppens F
14. Dhondt S
15. De Milde L
16. Vermeersch M
17. Maleux K
18. Gevaert K
19. Jerzmanowski A
20. Benhamed M
21. Wagner D
22. Vandepoele K
23. De Jaeger G
24. Inzé D
(2014) ANGUSTIFOLIA3 binds to SWI/SNF chromatin remodeling complexes to regulate transcription during Arabidopsis leaf development
The Plant Cell 26:210–229.

https://doi.org/10.1105/tpc.113.115907
- PubMed
- Google Scholar
Software
1. Vilhjalmsson BJ
(2019) mixmogam, version 766b889
GitHub.

https://github.com/bvilhjal/mixmogam
(2005) Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis
The Plant Journal 41:195–211.

https://doi.org/10.1111/j.1365-313X.2004.02288.x
- PubMed
- Google Scholar
Preprint
1. Weiszmann J
2. Clauw P
3. Jagoda J
4. Reichardt-Gomez I
5. Koemeda S
6. Jez J
7. Nordborg M
8. Walther D
9. Nägele T
10. Weckwerth W
(2020) Plasticity of the primary metabolome in 241 cold grown Arabidopsis thaliana accessions and its relation to natural habitat emperature
bioRxiv.

https://doi.org/10.1101/2020.09.24.311092
- Google Scholar
1. Westoby M
(1998) A leaf-height-seed (LHS) plant ecology strategy scheme
Plant and Soil 199:213–227.

https://doi.org/10.1023/A:1004327224729
- Google Scholar
1. Whitlock MC
(2008) Evolutionary inference from QST
Molecular Ecology 17:1885–1896.

https://doi.org/10.1111/j.1365-294X.2008.03712.x
- PubMed
- Google Scholar
1. Wieters B
2. Steige KA
3. He F
4. Koch EM
5. Ramos-Onsins SE
6. Gu H
7. Guo YL
8. Sunyaev S
9. de Meaux J
(2021) Polygenic adaptation of rosette growth in Arabidopsis thaliana
PLOS Genetics 17:e1008748.

https://doi.org/10.1371/journal.pgen.1008748
- PubMed
- Google Scholar
Preprint
1. Yim C
2. Bellis ES
3. DeLeo VL
4. Gamba D
5. Muscarella R
6. Lasky JR
(2022) Climate biogeography of Arabidopsis thaliana: Linking distribution models, individual performance, and life historyclimate biogeography of Arabidopsis thaliana
bioRxiv.

https://doi.org/10.1101/2022.03.06.483202
- Google Scholar
(2011) Fitness benefits and costs of cold acclimation in Arabidopsis thaliana
The American Naturalist 178:44–52.

https://doi.org/10.1086/660282
- PubMed
- Google Scholar
1. Zhou X
2. Stephens M
(2012) Genome-wide efficient mixed-model analysis for association studies
Nature Genetics 44:821–824.

https://doi.org/10.1038/ng.2310
- PubMed
- Google Scholar
1. Zuther E
2. Schulz E
3. Childs LH
4. Hincha DK
(2012) Clinal variation in the non-acclimated and cold-acclimated freezing tolerance of Arabidopsis thaliana accessions
Plant, Cell & Environment 35:1860–1878.

https://doi.org/10.1111/j.1365-3040.2012.02522.x
- PubMed
- Google Scholar

Article and author information

Author details

Pieter Clauw

Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna BioCenter, Vienna, Austria

Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Writing - original draft

For correspondence
pieter.clauw@gmi.oeaw.ac.at

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9677-8727
Envel Kerdaffrec

Department of Biology, University of Fribourg, Fribourg, Switzerland

Contribution
Resources

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-8667-6850
Joanna Gunis

Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna BioCenter, Vienna, Austria

Contribution
Investigation

Competing interests
No competing interests declared
Ilka Reichardt-Gomez

Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany

Contribution
Investigation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7964-9352
Viktoria Nizhynska

Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna BioCenter, Vienna, Austria

Contribution
Investigation

Competing interests
No competing interests declared
Stefanie Koemeda

Plant Sciences Facility, Vienna BioCenter Core Facilities GmbH, Vienna, Austria

Contribution
Investigation

Competing interests
No competing interests declared
Jakub Jez

Plant Sciences Facility, Vienna BioCenter Core Facilities GmbH, Vienna, Austria

Contribution
Supervision, Investigation

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6481-4383
Magnus Nordborg

Gregor Mendel Institute of Molecular Plant Biology, Austrian Academy of Sciences, Vienna BioCenter, Vienna, Austria

Contribution
Conceptualization, Supervision, Funding acquisition, Project administration, Writing - review and editing

For correspondence
magnus.nordborg@gmi.oeaw.ac.at

Competing interests
Reviewing editor, eLife

"This ORCID iD identifies the author of this article:" 0000-0001-7178-9748

Funding

No external funding was received for this work.

Acknowledgements

We thank past and current members of the Nordborg group for their help in setting up and harvesting these experiments. Thanks also to Daniele Filiault for her valuable comments on the manuscript.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.