# Leaf economics and slow-fast adaptation across the geographic range of Arabidopsis thaliana

Kevin Sartori, François Vasseur, Cyrille Violle, Etienne Baron, Marianne Gerard, Nick P Rowe, Oscar Ayala-Garay, Ananda Christophe, Laura Garcia De Jalón, Diane Masclef, et al.

► To cite this version:

Kevin Sartori, François Vasseur, Cyrille Violle, Etienne Baron, Marianne Gerard, et al. Leaf economics and slow-fast adaptation across the geographic range of Arabidopsis thaliana. Scientific Reports, 2019, 9 (1), pp.10758. 10.1038/s41598-019-46878-2. hal-02411727

HAL Id: hal-02411727

https://hal.science/hal-02411727

Submitted on 16 Dec 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Leaf economics and slow-fast adaptation across the geographic range of *Arabidopsis thaliana*

Received: 20 August 2018

Accepted: 1 July 2019

Published online: 24 July 2019

Kevin Sartori<sup>1</sup>, François Vasseur<sup>1,2</sup>, Cyrille Violle<sup>1</sup>, Etienne Baron<sup>1</sup>, Marianne Gerard<sup>1</sup>, Nick Rowe<sup>3</sup>, Oscar Ayala-Garay<sup>2,4</sup>, Ananda Christophe<sup>1</sup>, Laura Garcia de Jalón<sup>1</sup>, Diane Masclef<sup>2</sup>, Erwan Harscouet<sup>1</sup>, Maria del Rey Granado<sup>1</sup>, Agathe Chassagneux<sup>1,5</sup>, Elena Kazakou<sup>1,6</sup> & Denis Vile<sup>2</sup>

Life history strategies of most organisms are constrained by resource allocation patterns that follow a 'slow-fast continuum'. It opposes slow growing and long-lived organisms with late investment in reproduction to those that grow faster, have earlier and larger reproductive effort and a short longevity. In plants, the Leaf Economics Spectrum (LES) depicts a leaf-level trade-off between the rate of carbon assimilation and leaf lifespan, as stressed in functional ecology from interspecific comparative studies. However, it is still unclear how the LES is connected to the slow-fast syndrome. Interspecific comparisons also impede a deep exploration of the linkage between LES variation and adaptation to climate. Here, we measured growth, morpho-physiological and life-history traits, at both the leaf and whole-plant levels, in 378 natural accessions of *Arabidopsis thaliana*. We found that the LES is tightly linked to variation in whole-plant functioning, and aligns with the slow-fast continuum. A genetic analysis further suggested that phenotypic differentiation results from the selection of different slow-fast strategies in contrasted climates. Slow growing and long-lived plants were preferentially found in cold and arid habitats while fast growing and short-lived ones in more favorable habitats. Our findings shed light on the role of the slow-fast continuum for plant adaptation to climate. More broadly, they encourage future studies to bridge functional ecology, genetics and evolutionary biology to improve our understanding of plant adaptation to environmental changes.

Investigations of the genetic and phenotypic differentiation of plant lineages along environmental gradients is crucial for the understanding of plant evolutionary responses to current and future climate variations. As plants cannot simultaneously optimize competing eco-physiological functions, an important question is how plant adaptation occurs under the influence of major trade-offs between traits.

The slow-fast continuum is a pervasive trade-off between resource allocation to growth, reproduction and survival, spread across the tree of life<fnref n="1" />. The slow end of this continuum is characterized by slow growing, long-lived species and low reproductive output, while species at the fast end reach reproductive maturity faster and produce more offsprings. In plants, the leaf economics spectrum (LES hereafter)²⁻<fnref n="4" /> is thought to reflect the physiological basis of the slow-fast continuum<fnref n="4" />. The LES arrays plant species along a continuum of leaf trait syndromes going from short-lived leaves with fast metabolism to the reverse syndrome³. Core LES traits include leaf dry mass per area (LMA), leaf lifespan (LLS) and net photosynthetic rate per mass unit ($A_{max}$)³,⁵⁻⁷. LES traits are widely used in comparative ecology to infer whole-plant ecological strategies<fnref n="4" />,⁸⁻¹². However, the extent to which leaf-level resource economics reflects whole-plant physiology, performance, and ultimately fitness, is still under debate¹³. Many processes can lead to a mismatch between LES and whole-plant functioning¹⁴, including the impact of self-shading among leaves and resource allocation patterns, such as carbon investment in non-photosynthetic tissues¹⁵,¹⁶. To gain insights into the robustness of the slow-fast continuum at different organizational levels, we need

<footnote n="1">¹Univ Montpellier, CNRS, EPHE, IRD, Univ Paul Valéry Montpellier 3, Montpellier, France. ²Univ Montpellier, INRA, Montpellier SupAgro, LEPSE, Montpellier, France. ³Univ Montpellier, CIRAD, CNRS, INRA, IRD, Montpellier, France.</footnote>

<footnote n="4">⁴Programa de Recursos Genéticos y Productividad (RGP)-Fisiología Vegetal, Colegio de Postgraduados, 56230, Texcoco, Mexico. ⁵Office National de la Chasse et de la Faune sauvage, DRE Unité, Ongulés sauvages, Birieux, France.</footnote>

<footnote n="6">⁶Univ Montpellier, INRA, Montpellier SupAgro, Montpellier, France. Kevin Sartori and François Vasseur contributed equally. Correspondence and requests for materials should be addressed to K.S. (email: kevinfrsartori@gmail.com)</footnote>
to examine how LES traits scale up to plant level resource-use strategies, life history and performance. However, it remains difficult to compare individual performance across species with different growth forms, phenology and dispersal strategies since cross-species comparisons impede a clear linkage between physiological and adaptive trade-offs<sup>17-19</sup>.

The LES has been associated with differences in the ability of plants to adapt to more or less harsh environmental conditions<sup>4,12,20,21</sup>; species displaying high photosynthetic, respiration and growth rates, are short-lived, with thin and nitrogen-rich leaves are preferentially found in nutrient-rich and/or growth-suitable climatic conditions. Those species are qualified as acquisitive species in contrast to conservative ones that exhibit the opposite set of traits. Despite these observations, functional ecology has no tools to test for adaptation, and empirical evidences of the adaptive value of being at one end or the other of the continuum in a given environment remain scarce (see ref. 22 for a review). Furthermore, sampling procedure in field observation studies often impedes to disentangle the effects of plasticity vs. genetic differentiation on the emergence of the LES<sup>23</sup>. Thus, comparative studies looking for plant adaptation are at best incomplete<sup>22-24</sup>, and the role of selection in shaping the LES and driving adaptation to diverse environments is hardly understood. To fill this gap, intraspecific studies are encouraged since they can take benefit from tools developed in population ecology and genetics<sup>22,25,26</sup>. The LES has started to be analyzed at the intraspecific level, with contrasting findings depending on the studied organism and type of study<sup>23,27-32</sup>. LES relationships appeared consistent with cross-species ones when using species with broad environmental niche spectra<sup>31-33</sup> and/or broad phenotypic variability<sup>23</sup>, but inconsistent when using species with narrow phenotypic (and genetic) diversity<sup>34</sup>. Genetic differentiation of LES strategies has been demonstrated among populations of *Helianthus anomalus* along a 400 km rainfall gradient<sup>35</sup>. However, the question whether LES diversifies because of adaptation to climate among lineages spanning large geographic distribution remains open. Overall, we still miss a comprehensive understanding of within-species LES variation and the subsequent insights they can provide to well-described interspecific patterns from an evolutionary perspective.

Within-species genetic and phenotypic diversity is driven by natural selection, genetic drift, mutation and migration. The measure of $F_{ST}$ statistics among populations, i.e. a comparison of genetic or phenotypic differences among populations, is commonly used to distinguish neutral *versus* adaptive processes. At the genetic level, neutral loci are characterized by low $F_{ST}$ values, i.e. low between-population differentiation due to the combined effect of migration, mutation and drift, while high $F_{ST}$ values reflect genetic divergence as a consequence of directional selection in contrasted environments. From a phenotypic standpoint, adaptive divergence relies on the $Q_{ST}$ statistics, which was developed, by analogy to $F_{ST}$, to evaluate the population structure for quantitative traits. $Q_{ST}$ values above neutral $F_{ST}$ are interpreted as a signature of diversifying selection on the underlying trait. For instance, $Q_{ST}-F_{ST}$ comparisons have been successfully used in *Campanula rotundifolia*, *Arrhenatherum elatius*, *Quercus oleoides* and *Arabidopsis thaliana* to investigate the role of selection in the diversification of life-history traits, growth strategies and drought resistance among lineages at both local and global scales<sup>36-39</sup>. This method is expected to be particularly powerful in model species where modern genomics have given access to high-throughput genotyping<sup>40,41</sup>. By comparing genetic and phenotypic differentiation between populations or lineages, $Q_{ST}-F_{ST}$ provides a powerful tool to infer adaptation in polygenic quantitative traits such as LES traits<sup>41</sup>.

The species *Arabidopsis thaliana* has been widely used in molecular biology, cell biology and quantitative genetics. Thanks to the efforts to characterize the genetic diversity in this species<sup>42-45</sup>, it is also a model in population dynamics<sup>46</sup> and evolutionary ecology<sup>47</sup>. For instance, the genetic determinism of *A. thaliana* life history has been extensively studied, notably with the discovery of genes that control major developmental transitions such as flowering time (see ref. 48 for a review). Allelic variation in these genes appears to be adaptive to climatic and altitudinal gradients<sup>49</sup>. A recent study supports the hypothetic link between life history variation and the LES in *A. thaliana*, highlighted by strong genetic correlations between these traits<sup>28</sup>. However, this analysis was performed on recombinant inbred lines used for genetic mapping. Made of artificial crosses, they preclude examining the relationships between LES and the natural environment. Interestingly, *A. thaliana* has recently gained a renewed interest in functional ecology and biogeography<sup>50-52</sup>, notably due to the large panel of natural accessions that have been collected from contrasting climates, and genotyped at high density (e.g.<sup>42-45</sup>). As genetic data in *A. thaliana* allow an unprecedented large-scale analysis of genetic variation among populations and lineages, this species is promising to investigate the extent of intraspecific diversity and its role for adaptation to contrasted climates.

In this study, we explored the evolutionary bases of intraspecific leaf and plant trait variation using a pan-European collection of 378 natural *A. thaliana* accessions from the RegMap panel<sup>44</sup>. Specifically, we investigated whether plant adaptation to various climates is associated with genetic differentiation along the LES and the slow-fast continuum. To test this hypothesis, we first examined how the LES shapes phenotypic diversity across contrasted genotypes of *A. thaliana*, and tested whether LES traits scale up to plant level resource-use strategies and life history. Next, we took benefit from the large genomic information available in *A. thaliana* to evaluate to what extend phenotypic differences among lineages are attributable to adaptive processes such as adaptation to contrasted climates using $Q_{ST}-F_{ST}$ comparisons.

## Results

**Geographic clustering of *A. thaliana* lineages.** Following the Whittaker's biome classification method, two climatic variables, mean annual rainfall (MAR) and mean annual temperature (MAT), were used for study site classification into the major biome types. The range of biomes experienced by the sampled genotypes covers woodlands-shrublands and the less humid part of temperate forests, with a few extremes in boreal forests and deserts like environments (Fig. 1b). Using the 250 K SNPs data available from Horton *et al.*<sup>44</sup>, we performed a genetic clustering of the genotype set, which revealed the presence of five genetic lineages. These lineages were moderately differentiated (mean $F_{ST}$ = 0.11), geographically (Fig. 1a) as well as in the Whittaker's biome classification (Fig. 1b). The analysis revealed the existence of two genetic groups exclusively located in France in our sample (French 1 and French 2 hereafter) of 76 and 49 genotypes, respectively. Among the 20 genotypes of the
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location and climatic conditions of the genotype collecting sites. (a) Distribution of the 378 natural genotypes used in this study. The small points represent the collecting sites of genotypes and bigger points give the number of collecting sites overlapped at these positions. The colors represent the five genetic groups: Admixed (grey), French 1 (brown), French 2 (orange), Swedish (purple), Central Europe (dark blue), Western Europe (light blue). (b) Mean annual rainfall (MAR) and mean annual temperature (MAT) for the sites where genotypes were collected, in relation to major biome types of the world following Whittaker's classification. 1–9: Tundra, Boreal forest, Temperate Grassland Desert, Woodland Shrubland, Temperate Forest, Temperate Rain Forest, Tropical Forest Savana, Tropical Rain Forest, and Desert.</figcaption></figure>

third group, seven were defined as North Swedish in the 1001 genomes dataset⁴⁵. Consistently, the 20 "Swedish" genotypes, although not all in Sweden (Fig. 1a), were mainly located in cold environments and woodland-shrubland in Whittaker's classification (Fig. 1b). We considered the 76 genotypes from group 4 as "Central European" (Fig. 1a), typically living at intermediate temperatures and rainfall (Fig. 1b). 83 genotypes composed the group 5, all located in Western Europe (Fig. 1a), in a range of relatively warm environments with intermediate rainfall.

**Leaf economics of *A. thaliana*.** Assimilation rate ($A_{mass}$) was the most variable trait among the leaf economics traits in our dataset (15-fold; from 40.9 to 608.9 µmol g⁻¹ s⁻¹) while leaf mass per area (LMA) and leaf lifespan (LLS) varied 5 and 3.5 fold (from 18.7 to 101 g m⁻², and from 15 to 53.5 days), respectively. In regard to interspecific meta-analyses, variation in $A_{mass}$ was wide (from 5 to 660 µmol g⁻¹ s⁻¹³) and was relatively weak for LMA (from 14 to 1,500 g m⁻²³) and for LLS (from 0.9 to 288 months³). Pairwise regression revealed strong correlations between traits, independently of the genetic structure of the sample ($r²_{Amass-LLS}$ = 0.32, $r²_{Amass-LMA}$ = 0.73, $r²_{LMA-LLS}$ = 0.38, all p-value < 0.01). The 3-dimension space shaped by traits covariations in *A. thaliana* was comparable to the interspecific spectrum previously observed³: genotypes are ranked from low $A_{mass}$ and high LMA and LLS, toward high $A_{mass}$ and low LMA and LLS (Fig. 2). A principal component analysis (PCA) showed that 78% of the covariation between these three traits was explained by a single Principal Component (PC1; Fig. S1a). Hereafter, we assigned a position along the LES for each genotype according to its score on PC1. $A_{mass}$ was highly negatively correlated with PC1 ($r$ = −0.90) while LMA and LLS were positively correlated with PC1 ($r$ = 0.91 and 0.79, respectively). Thus, high and low PC1 values are representative of genotypes located at the conservative and acquisitive side of the LES, respectively.

**From the Leaf Economics Spectrum to the plant slow-fast continuum.** Trait measurements at the plant level revealed that assimilation rate was again the most variable trait with a 68-fold variation (from 8.4 to 578.1 µmol g⁻¹ s⁻¹), while plant mass per area and age of maturity (AM) both varied 5 fold (from 17.7 to 85.4 g m⁻² and 22 to 111 d, respectively). Standardized major axis (SMA) regressions between traits measured at the leaf and plant levels were all significant. Leaf and whole plant-level LMA were highly correlated ($r$ = 0.89; $P$ < 0.001; Fig. S1e) and the slope was close to, but significantly different from 1 (95% Confidence Interval slope = [1.07, 1.18]), as well as for leaf-level and whole plant-level net assimilation rate ($r$ = 0.78; $P$ < 0.001; Fig. S1f, 95% CI slope = [0.78, 0.9]). Similarly, LLS and AM were significantly correlated with a slope below 1 (slope = 0.68 [0.55; 0.65], $r$ = 0.67, $P$ < 0.001; Fig. S1g). A single principal component explained 86% of the trait covariation at the plant level (Fig. S1b) and was highly correlated with PC1 at the leaf level ($r$ = 0.87, $P$ < 0.001; Fig. S1h). Furthermore, life history and performance at the plant level co-varied in parallel with this leaf level phenotypic variation. Relative growth rate (RGR) and age at maturity (AM) were negatively correlated (Fig. 3) in our dataset, consistently with the slow-fast continuum pattern: fast growing genotypes reproduced earlier than slow growing ones. As for the LES, we performed a PCA with slow-fast traits and we assigned a position along the continuum for each genotype according to its score on PC1. The positions of the genotypes along the leaf economics spectrum and the slow-fast continuum were tightly correlated (Fig. 4a) and the slope did not differ
<figure><img src="image_3.png" /><figcaption>**Figure 2.** The leaf economics spectrum in *A. thaliana*. Three-way relationships among the main leaf economics traits: A<sub>mass</sub>, mass based assimilation rate (µmol CO<sub>2</sub> g<sup>-1</sup> s<sup>-1</sup>); LMA, leaf mass per area (g m<sup>-2</sup>); LLS, leaf lifespan (days). Each point represents a different genotype.</figcaption></figure>

significantly when taking the kinship matrix of the genotypes as a covariate. In other words, slow growing and late reproducing genotypes have low leaf carbon assimilation rate and long leaf lifespan, whereas fast growing and early reproducing genotypes have high leaf assimilation rate and short leaf lifespan.

**Adaptation cues of the slow-fast continuum.** Consistently with the original description of the RegMap genotypes⁴⁴, the genetic differentiation among *A. thaliana* populations was substantial. Pairwise comparisons based on the whole genetic data indicated a strong genetic divergence between French 2 and the other groups (F<sub>ST</sub> > 0.19). French 1 group was genetically closer to Western European group (F<sub>ST</sub> = 0.13) than French 2 group (F<sub>ST</sub> = 0.20). Interestingly, North Swedish lines showed strong phenotypic and genetic differentiation with other lineages (Fig. S4). We evaluated the significance of the differentiation of the slow-fast strategies performing Q<sub>ST</sub>-F<sub>ST</sub> comparisons for both RGR and AM. Heritability was high for both traits (h<sup>2</sup><sub>AM</sub> > 0.95, h<sup>2</sup><sub>RGR</sub> ~ 0.64). P<sub>ST</sub> is a good estimator of Q<sub>ST</sub> when trait heritability (h<sup>2</sup>) is high and the phenotypic between-populations component of variance genetically determined (c) is close to h<sup>2</sup> (see the Methods section). For the null hypothesis (H0: c = h<sup>2</sup> = 1), P<sub>ST</sub> values were high and significantly higher than the neutral F<sub>ST</sub> for both traits (P<sub>ST AM</sub> ~ 0.35, P<sub>ST RGR</sub> ~ 0.3; see Table 1). However, further analyses revealed that the estimation of P<sub>ST</sub> were highly sensitive to the ratio between the among-population (c) and overall heritability (h<sup>2</sup>) of the trait. On the other hand, we showed that all LES traits, exhibited modest and non-significant P<sub>ST</sub> values (Table 1) under the null assumption. Finally, we performed a P<sub>ST</sub> analysis on PCA scores for the two set of traits. Despite a lower heritability (h<sup>2</sup><sub>SFC</sub> ~ 0.59, h<sup>2</sup><sub>LES</sub> > 0.95), the slow-fast continuum P<sub>ST</sub> was high and significant compared to the leaf economics spectrum P<sub>ST</sub> (Fig. 4b) (P<sub>ST SFC</sub> ~ 0.35, critical c/h<sup>2</sup> ~ 0.35; P<sub>ST LES</sub> ~ 0.12), although both were highly sensitive to c/h<sup>2</sup>. This suggests that slow-fast traits behaved like outlier variants that diverged among lineages due to the effect of diversifying selection. Our data do not support the same hypothesis for LES traits, despite a tight coordination with slow-fast traits.

**Climatic drivers of *A. thaliana* phenotypes.** We investigated whether annual rainfall and temperature, two key climatic variables in plant functional biogeography, explain the position of the genotypes along the slow-fast continuum and leaf economics spectrum axes. We extracted 19 climatic variables at the collecting sites of the genotypes on CHELSA website (www.chelsa-climate.org/) related to temperature and precipitation, their temporal variability and extremes. Consistently with the population structure analysis, only the slow-fast continuum axis was significantly correlated with MAR (r = -0.18, P < 0.01) and MAT (r = -0.16, P < 0.01) (Fig. 5). The correlation with the 17 other CHELSA variables revealed the same pattern: the slow-fast continuum axis and the related traits were more often significantly correlated with climatic variables (Table S1). Therefore, we investigated whether climatic variables at the collecting sites of the genotypes can predict the *A. thaliana* slow-fast strategies. Firstly, stepwise regressions revealed that position on the slow-fast continuum was best predicted by a subset of climatic variables, including the mean annual temperature and rainfall, and variables related to temperature variance and extremes throughout the year. We evaluated the accuracy of the model using a repeated cross-validation method, which revealed that phenotypes were more accurately predicted by the reduced model (r<sup>2</sup> ~ 0.26) than the formal model (r<sup>2</sup> ~ 0.20). Extrapolating prediction of phenotypes across Europe from climate variables, we showed that slow strategies, characterized by slow growth and late reproduction, were favored in North Europe and Central East of Spain and in the highest European reliefs (Fig. 6). In addition, fast strategies characterized by fast growth and early reproduction were found in Central Europe and near the coasts.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** The slow-fast continuum in *A. thaliana*. Relationship between relative growth rate (mm² mm⁻² d⁻¹) and age at maturity (days). Each point represents a different genotype. The solid line represents the output of a linear model, the slopes do not significantly differ while taking into account the genetic relatedness of genotypes or not.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** The Leaf Economics Spectrum and the slow-fast continuum covariate in *A. thaliana* and are differentiated among populations. (a) Correlation between the position of each genotype along the leaf economics spectrum and the slow-fast continuum. Colors represent the five genetic groups (see Fig. 1). (b) Phenotypic differentiation (P<sub>ST</sub>, black vertical lines) and their confidence intervals (horizontal brackets) relatively to the distribution of genetic differentiation (F<sub>ST</sub>), its median (red dashed line), 90<sup>th</sup> quantile (light blue) and 95<sup>th</sup> quantile (dark blue).</figcaption></figure>

<table><caption>**Table 1.** The phenotypic differentiation of *A. thaliana* regarding slow-fast and LES traits is genetically determined and structured among regional populations. Abbreviations: h², narrow-sense heritability; P<sub>ST</sub>, phenotypic differentiation statistic; CI, confidence interval; c, among-population heritability.</caption><thead><tr><th></th><th>LES</th><th>SFC</th><th>A<sub>mass</sub></th><th>LLS</th><th>LMA</th><th>AM</th><th>RGR</th></tr></thead><tbody><tr><td>h<sup>2</sup></td><td>0.99</td><td>0.59</td><td>0.70</td><td>0.66</td><td>0.88</td><td>0.99</td><td>0.64</td></tr><tr><td>P<sub>ST</sub></td><td>0.12</td><td>0.35</td><td>0.12</td><td>0.02</td><td>0.19</td><td>0.35</td><td>0.3</td></tr><tr><td>CI 95%</td><td>[-0.02;0.28]</td><td>[0.07;0.7]</td><td>[-0.01;0.29]</td><td>[-0.05;0.09]</td><td>[0.01;0.42]</td><td>[0.07;0.71]</td><td>[0.03;0.62]</td></tr><tr><td>Critical c/h<sup>2</sup></td><td>NA</td><td>1.25</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>1.49</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Slow-fast continuum, Leaf Economics Spectrum, and their relationships with climate. Position of genotypes along the slow-fast continuum as a function of (a) mean annual temperature (MAT) and (b) mean annual rainfall (MAR) and position of genotypes along the Leaf Economics Spectrum as a function of (c) MAT and (d) MAR. Linear regressions are represented by black lines when significant.</figcaption></figure>

Discussion

The comparison of multiple species based on a few traits is the historical approach of functional ecology⁵³. While fruitful¹², such an approach impedes a deeper investigation of how evolutionary forces and trade-offs operate together to shape the observable phenotypic diversity²⁴,²⁶. Notably, several trait-trait covariations have been discussed in functional ecology in the light of trade-off theories. One of the most prominent phenotypic pattern discussed in the last decades, the so-called Leaf Economics Spectrum (LES), is thought to reflect a trade-off between metabolic rate and lifespan at the leaf level³,⁵,⁵⁴,⁵⁵. Plant species that exhibit long-lived leaves have been referred as resource conservative species. They optimize long-term carbon gain and extended nutrient residence time, as well as nutrient use efficiency⁵⁶. By contrast species with short-lived leaves sacrifice nutrient retention to maximize the rate of carbon fixation. The LES is expected to reflect an adaptive trade-off between fast and slow growth strategies across plant species⁴. Two assumptions underline this assertion: (i) the negative correlation between leaf photosynthetic rate and leaf lifespan is translated into a negative correlation between plant growth rate and the duration of the life cycle, (ii) particular combinations of slow-fast traits are selected in different environments. Both assumptions are difficult to test at the interspecific level. This has generated a living debate about the evolutionary causes of the LES²²,²⁸,⁵⁷⁻⁶⁰. Taking benefit from a large collection of sequenced genotypes in a model species, our results show that LES traits are correlated with slow-fast strategies at the plant level, and that trait divergence among genetic lineages is non-neutral. This supports the idea that plant populations evolve different slow-fast strategies along with different LES traits in order to adapt to contrasting climates.

We showed that LES trait correlations in *A. thaliana* follow the interspecific pattern³,⁵: individuals that invest a large amount of biomass per unit leaf area have a lower leaf assimilation rate and a longer leaf lifespan than plants that invest less biomass per unit leaf area. Moreover, the economics spectrum is still detected when scaling from leaf to whole-plant traits. This gives strong support to the idea that, at least in an herbaceous species such as *A. thaliana*, a trait value obtained on a single leaf using a standardized method, reflects the average phenotypic value expressed by all the leaves of an individual plant⁶¹,⁶². Furthermore, our results showed that the ranking of genotypes was similar along the leaf economics spectrum and the slow-fast continuum, suggesting that carbon economy at the leaf level is connected to the slow-fast strategies at the plant level. Up to now, functional ecology has favored interspecific comparisons, focusing on species trait means⁶³, with the perspective of generalization
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Slow strategy is favored in Mediterranean and Scandinavian regions while fast strategy is selected in the center of *A. thaliana* distribution. Map representing the prediction of slow-fast strategies based on climatic data: darker regions endure combination of climatic variables associated with slower strategies as opposed to faster strategies in lighter regions.</figcaption></figure>

and prediction of the whole biota functioning⁶⁴. Conversely, model species and annual ones in particular, have often been considered as extreme and atypical strategies when examining the phenotypic space of the worldwide flora¹¹. Here we showed how the eco-physiological examination of these species enrich our interpretation of interspecific trait covariations. Our results thus encourage future studies analyzing intraspecific trait covariations using model species such as *A. thaliana* in plant functional ecology to further advance our understanding of their underlying origins and mechanisms. Nonetheless, it is also true that our findings can be specific to the relatively simple organization and functioning of an annual rosette species. In particular, the correlations between leaf-level and whole-plant traits are presumably strongly variable among species. This relation is notably expected to be weaker in woody species because of the varying proportion of non-photosynthesizing tissues⁶⁵. This result has a consequence for mass-based traits: when a given leaf trait increases by one unit, the same trait at the plant level increases by less than one unit. Our results illustrate this statement: leaf A<sub>mass</sub> correlated with plant A<sub>mass</sub> with a slope below of 1. Similarly, leaf lifespan varied less across genotypes than plant age at maturity. As a result, there is a room for a decoupling between leaf and plant life history, even though the ranking among genotypes is globally conserved at the two organizational scales. Further explorations of how much leaf-level trade-offs and plant functioning are coordinated in herbaceous and woody species are needed.

Despite the autogamous nature of *A. thaliana* which makes it peculiar in eco-evolutionary studies, its large geographic distribution allows to explore multiple trait-environment relationships at large scales. These relationships are still hardly known and quantified for most functional traits and in most biomes¹²,⁶⁴ because of the myriad of species-specific responses that can blur the general interspecific trends. Again, our results illustrate how intraspecific comparisons can help elucidating the mechanisms underlying these relationships, notable putative adaptive ones. Using the Q<sub>ST</sub>-F<sub>ST</sub> comparisons, we brought evidence that *A. thaliana* regional populations are structured according to slow-fast traits, although better characterization of the genetic determinism of these phenotypes is needed. We reinforced this assertion identifying climatic factors related to this structure. The predicted distribution of slow-fast strategies across Europe revealed differential selection between roughly Norway, Sweden and Spain on one side, and central and Western Europe on the other. Selection for slow genotypes toward higher latitude in *A. thaliana*, specifically in North Swedish genotypes, is supported by previous findings on flowering time⁶⁶,⁶⁷. More surprisingly however, our results suggest that similar trait combinations representative of slow strategies are selected in two contrasted climates: Spain and Scandinavia, which are at the opposite edges of the *A. thaliana* latitudinal range. This clustering of *A. thaliana* genotypes echoes a recent study showing fixation of drought-related alleles in both Scandinavian and Spanish *A. thaliana* populations⁶⁸. If we consider together the absence of significant effect of the kinship matrix on trait-trait relationships tested, the globally low average differentiation among genetic groups (F<sub>ST</sub> = 0.11), and the phenotypic similarity observed at two distant locations, our results suggest that the genetic determinism of slow strategies as well as phenotypic differentiation could have occurred by convergence through adaptive processes. Thus, slow strategies could be selected in response to environmental stress in regions from nonetheless very different climates: low average temperature at Scandinavian sites and Mediterranean climate at Spanish sites. Interspecific studies at global scale revealed a negative relationship between conservative strategies and rainfall³,⁶⁹, possibly linked to a higher investment in cell wall complex macromolecules to face drought stress⁷⁰. Large-scale interspecific studies also reported a bias toward acquisitive strategies with increasing temperature in herbaceous species⁷¹,⁷². Similarly, theoretical frameworks suggest that drought and cold favor slow growing individuals in environments limited by water or temperature⁷³. Together, this suggests a general selection pressure for slow strategies in stressful environments, as opposed to selection for fast strategies in non-stressing environments²³. A promising avenue for future studies is to perform reciprocal transplant experiments to test the effect of drought and high temperatures on *A. thaliana* genotypes distributed along the slow-fast continuum.
Overall, our results suggest that slow-fast strategies are differentially selected in contrasted climates. Despite strong coordination of strategies from leaf level to plant level, slow-fast traits were significantly differentiated across populations and were associated with climatic conditions at the collecting sites, while significance was lacking for LES traits. By contrast, previous experiments under controlled conditions reported significant effects of temperature and soil water content on both whole plant and leaf traits in *A. thaliana*<sup>74,75</sup>. Complex interactions between climatic variables could change their relationships with phenotype in natural environments. This echoes the long standing search for global pattern of covariation between climate and LES traits in functional ecology and biogeography<sup>3,71,76</sup>. In addition to the expected anisometric relationship between plant and leaf traits (slope ≠1), room is left for leaf level traits to desynchronize from individual strategy. This suggests a potential for leaf traits to adapt to microclimatic conditions, including those at the canopy level. More broadly, our results corroborate the weak predictive power of leaf-level traits only in functional ecology when searching for adaption to particular climatic variable combination<sup>12,13</sup>. Indeed, even for an annual herb such as *A. thaliana*, whole-plant traits are more likely to capture adaptation to the environment compared to organ-level traits.

Using a model species, with large collections of well-characterized genetic material, appears particularly successful to go deeper into the evolutionary underpinning of major eco-physiological trade-offs, such as the LES and the slow-fast continuum. Combined with global climatic data, our findings notably revealed that adaptation to cold or dry habitats tends to favor slow strategies in *A. thaliana*. Next steps will be to merge approaches, and fully benefit from what a model species can provide both to genetics and ecophysiology. For instance, the climatic cues detected here despite the lack of climate data precision, is encouraging for the future of functional biogeography<sup>64</sup>. There is also evidence that the connection between functional trait and environmental adaptation requires a better characterization of plant fitness through demographic measures<sup>26</sup>. Comparative studies integrating demographic approach at population level are promising to understand how selection and macro-ecological gradients shape the evolutionary responses of plants to climate variation<sup>24,26</sup>.

## Materials and Methods

**Plant material.** We used a total of 378 natural genotypes of *A. thaliana* L. Heynh sampled from the worldwide lines of the RegMap population (http://bergelson.uchicago.edu/wp-content/uploads/2015/04/Justins-360-lines.xls), which were genotyped for 250 K bi-allelic SNPs<sup>44</sup>. Seeds used in the present study were originally collected and multiplied by Fabrice Roux (UMR LIPM, Toulouse, France) and his collaborators.

**Growth conditions.** Phenotype characterization was performed under controlled conditions in the high-throughput PHENOPSIS phenotyping platform<sup>77</sup> to track daily growth. Seeds were kept in the dark at 4°C for at least one week before sowing. Four to six seeds per genotype were sown at the soil surface in 225 ml pots filled with a 1:1 (v:v) mixture of loamy soil and organic compost (Neuhaus N2). The soil surface was moistened with one-tenth strength Hoagland solution, and pots were kept in the dark during 48 h under controlled environmental conditions (20°C, 70% air relative humidity). Then, pots were placed in the PHENOPSIS growth chamber at 20°C, 12 h photoperiod, 70% relative humidity, 175 µmol m⁻² s⁻¹ PPFD. Pots were sprayed with deionized water three times per day until germination, and then soil water content was adjusted to 0.35 g H₂O g⁻¹ dry soil (−0.07 MPa soil water potential) to ensure optimal growth<sup>74,78</sup>. After emergence of the fourth leaf, one individual plant was left in each pot.

**Measurements of plant traits.** In order to standardize measurements for all genotypes, all traits were quantified when flower buds were macroscopically visible (i.e. bolting stage), and leaf traits were measured on the last adult leaf, fully exposed to light.

Net photosynthetic rate, relative expansion rate, lifespan, vegetative dry weight, as well as leaf area were determined for the leaf and the plant canopy. Net photosynthetic rate was measured at leaf (leaf A, nmol CO₂ s⁻¹) and whole-plant levels (plant A, nmol CO₂ s⁻¹) under growing conditions using, respectively, the leaf cuvette provided with the infrared gas analyzer system (CIRAS 2, PP systems, USA), and a whole-plant chamber prototype designed for *A. thaliana* by M. Dauzat (INRA, Montpellier, France) and K. J. Parkinson (PP System, UK) (see<sup>28</sup>). Leaf and whole-plant photosynthetic rates were both expressed on dry mass basis (leaf A<sub>mass</sub> and plant A<sub>mass</sub>, nmol CO₂ g⁻¹ s⁻¹). Due to time constraints, we measured photosynthetic rates for 319 and 348 genotypes at the leaf and whole-plant levels (306 in common), respectively. We estimated the age of maturity by the number of days from germination to the appearance of the flower bud. Then, plants were harvested, and individual fresh weight was determined. The leaf used for photosynthetic measurements was identified and processed separately, and detached rosettes were kept in deionized water at 4°C for 24 h, and water-saturated weight was determined. Individual leaves were then attached to a sheet of paper and scanned for subsequent determination of the leaf number and total leaf area using ImageJ<sup>79</sup>. Dry weight of laminas and petioles were obtained after drying for 72 h at 65°C. Rosette dry weight was expressed as the sum of lamina and petiole dry weights. Leaf mass per area was both calculated for the leaf used for photosynthetic measurements (LMA, g m⁻²) and for the whole-rosette (plant LMA, g m⁻²) as the ratio of lamina dry mass to lamina area. Relative growth rate (RGR, mm² mm⁻² d⁻¹) and leaf lifespan (LLS, d) were estimated from automated daily pictures of the rosettes. More precisely, a sigmoid curve was fitted to rosette area as a function of time in order to extract growth parameters, where RGR was calculated as the slope at the inflection point<sup>80-82</sup>. Using daily pictures, we tracked three consecutive leaves from birth (emergence) to death (full senescence). For each plant, leaf duration was calculated as the average number of days from leaf emergence to senescence.

**F<sub>ST</sub> and P<sub>ST</sub> estimates.** In order to perform population genetic analyses, genetic groups were identified by genetic clustering of 378 genotypes, using the 250 K SNPs data available from Horton *et al.*<sup>44</sup>. Clustering was performed with ADMIXTURE<sup>83</sup> after linkage disequilibrium pruning (*r*² < 0.1 in a 50 kb window with a step size of
50 SNPs) with PLINK⁸⁴, resulting in 24,562 independent SNPs used for subsequent analyses. A cross-validation for different numbers of clusters (k = 1 to k = 10) showed that the set of studied genotypes can be separated into five groups representative of different genetic lineages (cross validation error = 0.89). Following the same approach as the 1001 genomes project⁴⁵, we assigned each genotype to a group if more than 50% of its genome derived from the corresponding cluster. The 74 genotypes not matching this criterion were labelled “Admixed” and were not used for the FST and PST calculation. The groups genetically defined were also geographically distinct as shown by early studies⁴⁴. We measured the genetic population structure using Weir and Cockerham FST statistic for all the 24,562 SNPs, as well as mean FST genome-wide. To determine the neutral FST value, we calculated the median of the significant FST values of the intergenic SNPs. More specifically, we filtered the SNPs dataset keeping the intergenic SNPs only, assuming that they are rather prone to endure neutral differentiation processes. We then calculated a FST value for each intergenic SNP with the attribution of population described above. To test for the significance of the FST values, we randomized 1000 times the population attribution to the genotypes in order to provide a null distribution of FST. FST values higher than the 95th quantile of their null distribution were stated as significant. The phenotypic population structure is ideally evaluated using QST, an analogue of FST measure⁸⁵,⁸⁶. We estimated QST using a phenotype-based surrogate (PST)⁸⁶ which depends on among-population and overall heritability of the trait (c and h², respectively) computed as follow:

$$P_{ST} = \frac{c \sigma_{B}^{2}}{c \sigma_{B}^{2} + h^{2} \sigma_{W}^{2}}$$

where B and W are the between- and within-population part of variance, respectively. Note that within-population variance is not multiplied by 2 since A. thaliana is mainly autogamous and homozygote at all loci. Under the null assumption (H0: c = h² = 1), PST and QST are analogous. Sensitivity analysis consists in calculating the PST value and its 95% confidence interval for a gradient of hypothetical c and a given h². PST is a robust estimator of QST when the critical value of c/h² is low, i.e. the minimum value of c/h² for which the lower limit of PST 95% CI is higher than the neutral FST. Then, a value of QST higher than neutral FST means that the phenotypic differentiation between populations is larger than expected by demographic events alone, in particular genetic drift, and is thus indicative of diversifying selection on traits⁴¹,⁸⁷. We used parametric bootstrap method to generate 95% CI around PST values with the package MCMCglmm in R (10,000 iterations).

**Statistical analysis.** Climate variables at the sampling sites of each genotype were extracted from the CHELSA database (http://www.chelsa-climate.org/), with a 2.5 arc-minutes resolution. The effect of climatic variables on traits was tested using linear model regressions. All analyses were performed in R 3.4.1 (R Core Team, 2017). Whittaker's biomes were plotted using the BIOMEplot function provided by G. Kunstler (https://rdrr.io.github/kunstler/BIOMEplot/src/R/biomes-plot.R). All leaf and plant traits, but RGR, were log₁₀ transformed when Gaussian distribution is required for statistical analyses. Principal component analysis (PCA) was performed using the package FactoMineR. The package nlme was used to perform linear models and phylogenetic generalized least squares regressions. We performed phylogenetic regressions including a relatedness matrix as covariance matrix, obtained after running the PLINK-make-rel command across the 250 K SNPs from the RegMap data. SMA regressions between leaf and plant traits were performed using the package SMATR⁸⁸, and phylogenetic SMA regressions using the Phyl.RMA function of the Phytools package. The phylogenetic tree required for SMA regression has been produced with Tassel using the RegMap SNPs data⁴⁴.

## Data Availability

Data are deposited in the Dataverse repository: https://doi.org/10.15454/B3W0OS.

## References

1. Dammhahn, M., Dingemanse, N. J., Niemelä, P. T. & Réale, D. Pace-of-life syndromes: a framework for the adaptive integration of behaviour, physiology and life history. *Behavioral Ecology and Sociobiology* 72 (2018).

2. Reich, P. B. *et al.* Generality of Leaf Trait Relationships: A Test across Six Biomes. *Ecology* 80, 1955 (1999).

3. Wright, I. J. *et al.* The worldwide leaf economics spectrum. *Nature* 428, 821–827 (2004).

4. Reich, P. B. The world-wide 'fast-slow' plant economics spectrum: a traits manifesto. *Journal of Ecology* 102, 275–301 (2014).

5. Reich, P. B., Walters, M. B. & Ellsworth, D. S. From tropics to tundra: Global convergence in plant functioning. *Proceedings of the National Academy of Sciences* 94, 13730–13734 (1997).

6. Reich, P. B. *et al.* The Evolution of Plant Functional Variation: Traits, Spectra, and Strategies. *International Journal of Plant Sciences* 164, S143–S164 (2003).

7. Lloyd, J., Bloomfield, K., Domingues, T. F. & Farquhar, G. D. Photosynthetically relevant foliar traits correlating better on a mass vs an area basis: of ecophysiological relevance or just a case of mathematical imperatives and statistical quicksand? *New Phytologist* 199, 311–321 (2013).

8. Westoby, M., Falster, D. S., Moles, A. T., Vesk, P. A. & Wright, I. J. Plant Ecological Strategies: Some Leading Dimensions of Variation Between Species. *Annual Review of Ecology and Systematics* 33, 125–159 (2002).

9. Pierce, S., Brusa, G., Vagge, I. & Cerabolini, B. E. L. Allocating CSR plant functional types: the use of leaf economics and size traits to classify woody and herbaceous vascular plants. *Functional Ecology* 27, 1002–1010 (2013).

10. Pierce, S., Vagge, I., Brusa, G. & Cerabolini, B. E. L. The intimacy between sexual traits and Grime's CSR strategies for orchids coexisting in semi-natural calcareous grassland at the Olive Lawn. *Plant Ecology* 215, 495–505 (2014).

11. Diaz, S. *et al.* The global spectrum of plant form and function. *Nature* 529, 167–171 (2015).

12. Garnier, E., Navas, M.-L. & Grigulis, K. Plant functional diversity: organism traits, community structure, and ecosystem properties. (Oxford University Press 2016).

13. Craine, J. M., Engelbrecht, B. M. J., Lusk, C. H., McDowell, N. G. & Poorter, H. Resource limitation, tolerance, and the future of ecological plant classification. *Frontiers in Plant Science* 3 (2012).

14. Poorter, H., Lambers, H. & Evans, J. R. Trait correlation networks: a whole-plant perspective on the recently criticized leaf economic spectrum. *New Phytologist* 201, 378–382 (2014).

15. Horn, H. S. *The adaptive geometry of trees*. (Princeton Univ. Pr 1971).
16. Givnish, T. Adaptation to Sun and Shade: a Whole-Plant Perspective. *Australian Journal of Plant Physiology* **15**, 63 (1988).

17. Agrawal, A. A., Conner, J. K. & Rasmann, S. Tradeoffs and negative correlations in evolutionary ecology. In *Evolution After Darwin: the First 150 Years*. (Sinauer Associates 2010).

18. Adler, P. B. et al. Functional traits explain variation in plant life history strategies. *Proceedings of the National Academy of Sciences* **111**, 740–745 (2014).

19. Salguero-Gómez, R. Applications of the fast-slow continuum and reproductive strategy framework of plant life histories. *New Phytologist* **213**, 1618–1624 (2017).

20. Wright, I. J. et al. Modulation of leaf economic traits and trait relationships by climate. *Global Ecology and Biogeography* **14**, 411–421 (2005).

21. Heberling, J. M. & Fridley, J. D. Biogeographic constraints on the world-wide leaf economics spectrum: Biogeography and leaf economics. *Global Ecology and Biogeography* **21**, 1137–1146 (2012).

22. Donovan, L. A., Maherali, H., Caruso, C. M., Huber, H. & de Kroon, H. The evolution of the worldwide leaf economics spectrum. *Trends in Ecology & Evolution* **26**, 88–95 (2011).

23. Anderegg, L. D. L. et al. Within-species patterns challenge our understanding of the leaf economics spectrum. *Ecology Letters*. https://doi.org/10.1111/ele.12945 (2018).

24. Moles, A. T. Being John Harper: Using evolutionary ideas to improve understanding of global patterns in plant traits. *Journal of Ecology* **106**, 1–18 (2018).

25. Albert, C. H. et al. Intraspecific functional variability: extent, structure and sources of variation. *Journal of Ecology* **98**, 604–613 (2010).

26. Salguero-Gómez, R., Violle, C., Gimenez, O. & Childs, D. Delivering the promises of trait-based approaches to the needs of demographic approaches, and vice versa. *Functional Ecology* **32**, 1424–1435 (2018).

27. Albert, C. H. et al. A multi-trait approach reveals the structure and the relative importance of intra- vs. interspecific variability in plant traits: Intra- vs. interspecific variability in plant traits. *Functional Ecology* **24**, 1192–1201 (2010).

28. Vasseur, F., Violle, C., Enquist, B. J., Granier, C. & Vile, D. A common genetic basis to the origin of the leaf economics spectrum and metabolic scaling allometry. *Ecology Letters* **15**, 1149–1157 (2012).

29. Blonder, B. et al. Testing models for the leaf economics spectrum with leaf and whole-plant traits in Arabidopsis thaliana. *AoB Plants* **7**, plv049 (2015).

30. Hu, Y.-K. et al. Novel evidence for within-species leaf economics spectrum at multiple spatial scales. *Frontiers in Plant Science* **6** (2015).

31. Martin, A. R. et al. Intraspecific trait variation across multiple scales: the leaf economics spectrum in coffee. *Functional Ecology* **31**, 604–612 (2017).

32. Isaac, M. E. et al. Intraspecific Trait Variation and Coordination: Root and Leaf Economics Spectra in Coffee across Environmental Gradients. *Frontiers in Plant Science* **8** (2017).

33. Niinemets, Ü. Is there a species spectrum within the world-wide leaf economics spectrum? Major variations in leaf functional traits in the Mediterranean sclerophyll *Quercus ilex*. *New Phytologist* **205**, 79–96 (2015).

34. Blonder, B., Violle, C. & Enquist, B. J. Assessing the causes and scales of the leaf economics spectrum using venation networks in *Populus tremuloides*. *Journal of Ecology* **101**, 981–989 (2013).

35. Brouillette, L. C., Mason, C. M., Shirk, R. Y. & Donovan, L. A. Adaptive differentiation of traits related to resource use in a desert annual along a resource gradient. *New Phytologist* **201**, 1316–1327 (2014).

36. Preite, V., Stöcklin, J., Armbruster, G. F. J. & Scheepens, J. F. Adaptation of flowering phenology and fitness-related traits across environmental gradients in the widespread Campanula rotundifolia. *Evolutionary Ecology* **29**, 249–267 (2015).

37. Michalski, S. G., Malyshev, A. V. & Kreyling, J. Trait variation in response to varying winter temperatures, diversity patterns and signatures of selection along the latitudinal distribution of the widespread grassland plant *Arrhenatherum elatius*. *Ecology and Evolution* **7**, 3268–3280 (2017).

38. Dittberner, H. et al. Natural variation in stomata size contributes to the local adaptation of water-use efficiency in Arabidopsis thaliana. *Molecular Ecology* **27**, 4052–4065 (2018).

39. Ramírez-Valiente, J. A. et al. Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak *Quercus oleoides*. *Molecular Ecology* **27**, 2176–2192 (2018).

40. Merilä, J. & Crnokrak, P. Comparison of genetic differentiation at marker loci and quantitative traits: Natural selection and genetic differentiation. *Journal of Evolutionary Biology* **14**, 892–903 (2001).

41. Leinonen, T., McCairns, R. J. S., O'Hara, R. B. & Merilä, J. QST-FST comparisons: evolutionary and ecological insights from genomic heterogeneity. *Nature Reviews Genetics* **14**, 179–190 (2013).

42. Clark, R. M. et al. Common Sequence Polymorphisms Shaping Genetic Diversity in Arabidopsis thaliana. *Science* **317**, 338–342 (2007).

43. Platt, A. et al. The Scale of Population Structure in Arabidopsis thaliana. *PLoS Genetics* **6**, e1000843 (2010).

44. Horton, M. W. et al. Genome-wide patterns of genetic variation in worldwide Arabidopsis thaliana accessions from the RegMap panel. *Nature Genetics* **44**, 212–216 (2012).

45. Alonso-Blanco, C. et al. 1,135 Genomes Reveal the Global Pattern of Polymorphism in Arabidopsis thaliana. *Cell* **166**, 481–491 (2016).

46. Montesinos-Navarro, A., Wig, J., Xavier Pico, F. & Tonsor, S. J. Arabidopsis thaliana populations show clinal variation in a climatic gradient associated with altitude. *New Phytologist* **189**, 282–294 (2011).

47. Weigel, D. Natural Variation in Arabidopsis: From Molecular Genetics to Ecological Genomics. *Plant Physiology* **158**, 2–22 (2012).

48. Simpson, G. G. Arabidopsis, the Rosetta Stone of Flowering Time? *Science* **296**, 285–289 (2002).

49. Mendez-Vigo, B., Pico, F. X., Ramiro, M., Martinez-Zapater, J. M. & Alonso-Blanco, C. Altitudinal and Climatic Adaptation Is Mediated by Flowering Traits and FRI, FLC, and PHYC Genes in Arabidopsis. *Plant Physiology* **157**, 1942–1955 (2011).

50. Baron, E., Richirt, J., Villoutreix, R., Amsellem, L. & Roux, F. The genetics of intra- and interspecific competitive response and effect in a local population of an annual plant species. *Functional Ecology* **29**, 1361–1370 (2015).

51. Vasseur, F. et al. Adaptive diversification of growth allometry in the plant *Arabidopsis thaliana*. *Proceedings of the National Academy of Sciences* **115**, 3416–3421 (2018).

52. Takou, M. et al. Linking genes with ecological strategies in *A rabidopsis thaliana*. *Journal of Experimental Botany*, https://doi.org/10.1093/jxb/ery447 (2018).

53. Keddy, P. A. Assembly and response rules: two goals for predictive community ecology. *Journal of Vegetation Science* **3**, 157–164 (1992).

54. Kikuzawa, K. A Cost-Benefit Analysis of Leaf Habit and Leaf Longevity of Trees and Their Geographical Pattern. *The American Naturalist* **138**, 1250–1263 (1991).

55. Kikuzawa, K. The basis for variation in leaf longevity of plants. *Vegetatio* **121**, 89–100 (1995).

56. Reich, P. B., Walters, M. B. & Ellsworth, D. S. Leaf Life-Span in Relation to Leaf, Plant, and Stand Characteristics among Diverse Ecosystems. *Ecological Monographs* **62**, 365–392 (1992).

57. Mason, C. M. & Donovan, L. A. Evolution of the leaf economics spectrum in herbs: Evidence from environmental divergences in leaf physiology across *Helianthus* (Asteraceae): evolution of leaf in diverse sunflowers. *Evolution* **69**, 2705–2720 (2015).
58. Mason, C. M. et al. Macroevolution of leaf defenses and secondary metabolites across the genus Helianthus. *New Phytologist* **209**, 1720-1733 (2016).

59. Pilote, A. J. & Donovan, L. A. Evidence of correlated evolution and adaptive differentiation of stem and leaf functional traits in the herbaceous genus, Helianthus. *American Journal of Botany* **103**, 2096-2104 (2016).

60. Muir, C. D., Conesa, M. À., Roldán, E. J., Molins, A. & Galmés, J. Weak coordination between leaf structure and function among closely related tomato species. *New Phytologist* **213**, 1642-1653 (2017).

61. Garnier, E. et al. Consistency of species ranking based on functional leaf traits. *New Phytologist* **152**, 69-83 (2001).

62. Ang, K. Y., Lucas, P. W. & Tan, H. T. W. Novel way of measuring the fracture toughness of leaves and other thin films using a single inclined razor blade. *New Phytologist* **177**, 830-837 (2008).

63. Albert, C. H., Grassein, F., Schurr, F. M., Vieilledent, G. & Violle, C. When and how should intraspecific variability be considered in trait-based plant ecology? *Perspectives in Plant Ecology, Evolution and Systematics* **13**, 217-225 (2011).

64. Violle, C., Reich, P. B., Pacala, S. W., Enquist, B. J. & Kattge, J. The emergence and promise of functional biogeography. *Proceedings of the National Academy of Sciences* **111**, 13690-13696 (2014).

65. Poorter, H. & Remkès, C. Leaf area ratio and net assimilation rate of 24 wild species differing in relative growth rate. *Oecologia* **83**, 553-559 (1990).

66. Long, Q. et al. Massive genomic variation and strong selection in Arabidopsis thaliana lines from Sweden. *Nature Genetics* **45**, 884-890 (2013).

67. Li, P. et al. Multiple FLC haplotypes defined by independent cis -regulatory variation underpin life history diversity in *Arabidopsis thaliana*. *Genes & Development* **28**, 1635-1640 (2014).

68. Exposito-Alonso, M. et al. Genomic basis and evolutionary potential for extreme drought adaptation in Arabidopsis thaliana. *Nature Ecology & Evolution* **2**, 352–358 (2018).

69. Onoda, Y. et al. Global patterns of leaf mechanical properties: Global patterns of leaf mechanical properties. *Ecology Letters* **14**, 301-312 (2011).

70. Onoda, Y. et al. Physiological and structural tradeoffs underlying the leaf economics spectrum. *New Phytologist* **214**, 1447-1463 (2017).

71. Borgy, B. et al. Plant community structure and nitrogen inputs modulate the climate signal on leaf traits. *Global Ecology and Biogeography* **26**, 1138-1152 (2017).

72. Šimová, I. et al. Spatial patterns and climate relationships of major plant traits in the New World differ between woody and herbaceous species. *Journal of Biogeography*. https://doi.org/10.1111/jbi.13171 (2018).

73. Körner, C. Growth Controls Photosynthesis - Mostly. *Nova Acta Leopoldina* **114**, 273-283 (2013).

74. Vile, D. et al. Arabidopsis growth under prolonged high temperature and water deficit: independent or interactive effects?: Plant responses to high temperature and water deficit. *Plant, Cell & Environment* **35**, 702-718 (2012).

75. Vasseur, F., Bontpart, T., Dauzat, M., Granier, C. & Vile, D. Multivariate genetic analysis of plant responses to water deficit and high temperature revealed contrasting adaptive strategies. *Journal of Experimental Botany* **65**, 6457-6469 (2014).

76. Moles, A. T. et al. Which is a better predictor of plant traits: temperature or precipitation? *Journal of Vegetation Science* **25**, 1167-1180 (2014).

77. Granier, C. et al. Phenopsis, an automated platform for reproducible phenotyping of plant responses to soil water deficit in Arabidopsis thaliana permitted the identification of an accession with low sensitivity to soil water deficit. *New Phytol.* **169**, 623-635 (2006).

78. Aguirrezabal, L. et al. Plasticity to soil water deficit in Arabidopsis thaliana: dissection of leaf development into underlying growth dynamic and cellular variables reveals invisible phenotypes. *Plant, Cell and Environment* **29**, 2216-2227 (2006).

79. Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. *Nat. Methods* **9**, 671–675 (2012).

80. Dhondt, S. et al. High-resolution time-resolved imaging of *in vitro* Arabidopsis rosette growth. *The Plant Journal* **80**, 172-184 (2014).

81. Minervini, M., Giuffrida, M. V., Perata, P. & Tsaftaris, S. A. Phenotiki: an open software and hardware platform for affordable and easy image-based phenotyping of rosette-shaped plants. *The Plant Journal* **90**, 204-216 (2017).

82. Dobrescu, A., Scorza, L. C. T., Tsaftaris, S. A. & McCormick, A. J. A "Do-It-Yourself" phenotyping system: measuring growth and morphology throughout the diel cycle in rosette shaped plants. *Plant Methods* **13** (2017).

83. Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. *Genome Research* **19**, 1655-1664 (2009).

84. Purcell, S. et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. *The American Journal of Human Genetics* **81**, 559-575 (2007).

85. Lande, R. Neutral Theory of Quantitative Genetic Variance in an Island Model with Local Extinction and Colonization. *Evolution* **46**, 381 (1992).

86. Brommer, J. E. Whither Pst? The approximation of Qst by Pst in evolutionary and conservation biology: Whither Pst? *Journal of Evolutionary Biology* **24**, 1160-1168 (2011).

87. José, R.-V. A. et al. Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak *Quercus oleoides*. *Molecular Ecology*, https://doi.org/10.1111/mec.14566 (2018).

88. Warton, D. I., Wright, I. J., Falster, D. S. & Westoby, M. Bivariate line-fitting methods for allometry. *Biological Reviews* **81**, 259 (2006).

## Acknowledgements

We thank Myriam Dauzat, Alexis Bédiée and Gaëlle Rolland for their technical assistance during trait measurements and for the environmental control in PHENOPSIS. We thank Fabrice Roux and his group (Laboratory of Plant-Microbe Interactions, Toulouse, France) for the seeds of the natural accessions used in our study. This work was supported by INRA Dpt EA grant, the European Research Council (ERC) ('CONSTRAINTS': grant ERC-StG-2014-639706-CONSTRAINTS), and the French Agency for Research (ANR grant ANR-17-CE02-0018-01, 'AraBreed'). This publication has been written with the support of the Agreenskills fellowship program which has received funding from the EU's Seventh Framework Program under the agreement N° FP7-609398 (Agreenskills contract 3215 'AraBreed').

## Author Contributions

D.V., F.V. and C.V. designed the study, K.S., E.B., M.G., O.A.-G., A.C., L.G.D.J., D.M., E.H., M.D.R.G. and A.C. conducted the experiments. K.S. and F.V. performed statistical analyses. K.S. wrote the first draft of the manuscript, and K.S., F.V., C.V., E.B., M.G., N.R., O.A.-G., A.C., L.G.D.J., D.M., E.H., M.D.R.G., A.C., E.K. and D.V. contributed to revisions.
**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-46878-2.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019