Oak leaf morphology may be more strongly shaped by climate than by phylogeny

Key message Despite been grown under the same climate, oak species are able to correlate with looser, but still identifiable, leaf morphological syndromes, composed by morphological traits with an ecological role in their respec‑ tive macroclimates. Context Environmental restrictions imposed by climate have been shown to modulate leaf morphology. A reduction of leaf area in hot and dry climates reduces overheating because of a thinner boundary layer. Lobed shapes enhance hydraulic conductivity and faster cooling. Elongated leaves drain more quickly under high precipitation. Trichomes may help to reduce the effects of excessive sun exposure in hot and dry environments. Leaf mass per area (LMA) increases in stressful environments. Aim To assess the influence of global climate and clade on oak leaf syndromes comprising morphological traits with ecological roles. Methods We analyzed seven morphological traits in 141 oak species grown in a botanical garden, characterized into five macroclimates, and explored the partial effects of clade and climate. Results We found significant associations between macroclimate and every morphological trait measured. Tem‑ perate species tend to have large, obovate, lobed, malacophyllous leaves. Species occurring in dry habitats usually present small, rounded, pubescent, sclerophyllous leaves. Warm and wet climates induce the development of slender, lanceolate, glabrous leaves with an acuminate apex. Conclusion The functional roles of the different morphological traits are partially confirmed in genus Quercus as a response to the different macroclimates, where different leaf syndromes can be distinguished.


Introduction
Leaves are the primary photosynthetic organs of plants.Their morphology should be adapted to environmental conditions (Traiser et al. 2005), as seemingly fine differences in leaf form may strongly affect the capacity of plants to grow and survive under different environments (Givnish 1987;Desmond et al. 2021).It has long been recognized that natural selection by climatic restrictions may influence leaf morphology (Stahl 1880(Stahl , 1883;;Schimper 1903).In the past two decades, numerous studies have found correlations between leaf morphological traits (e.g.leaf size and shape) and environmental factors (mainly annual mean temperature and precipitation) from global to local scales (Royer et al. 2005;Traiser et al. 2005;Adams et al. 2008;Peppe et al. 2011;Li and Wang 2021).Leaf area, shape, pubescence, and dry mass per area in particular have proven to be of high functional significance (Alonso-Forn et al. 2020;Cavender-Bares and Ramírez-Valiente 2017;Du et al. 2022;Sancho-Knapik et al. 2021).
Leaf area varies and evolves in response to water availability, temperature, and other aspects of climate (Webb 1968;Dolph and Dilcher 1980;Givnish 1984;Peppe et al. 2011;Li et al. 2020).Small leaves have been linked to harsh conditions such as hot, dry and lightsaturating environments (Peguero-Pina et al. 2020;Nicotra et al. 2011) due to the thinner boundary layer of small leaves (Yates et al. 2010), which facilitates heat loss (Roth-Nebelsick 2001).Leaf temperature tends to be higher in large leaves due to their thicker boundary layers (Vogel 1970;Parkhurst and Loucks 1972;Givnish 1987;Vogel 2009).
Leaf shape-e.g.length-width ratio, toothing and lobation, apex structure-interacts with leaf area in adapting plants to climate and habitat.Leaf shape has various functional roles, including thermoregulation, hydraulic constraints, and leaf expansion patterns in deciduous species (Nicotra et al. 2011).Plants in cool regions tend to have broader leaves with lower length-width ratios to increase the thickness of their boundary layers and prevent wind cooling, maintaining higher leaf temperatures for photosynthesis (Wright et al. 2017;Li et al. 2020).In contrast, plants in warm and humid regions often have elliptic leaves with higher length-width ratios to avoid overheating and enhance water transport (Wright et al. 2017;Hu and Hu 2015), analogous to the effect of having smaller leaves.Mean annual temperature has long been shown to influence leaf toothing and lobation, with narrow and deeply lobated leaves better adapted to and more frequent in warmer climates due to thinner boundary layers and enhanced hydraulic efficiency (Sisó et al. 2001;Sack and Tyree 2005;Leigh et al. 2017).However, toothed leaves are more common in cool climates, where they are considered an adaptation that enhances sap flow and carbon gain early in the growing season (Bailey and Sinnott 1916;Baker-Brosh and Peet 1997;Royer and Wilf 2006;Nicotra et al. 2011;Peppe et al. 2011;MacKee et al. 2019;Li and Wang 2021).
Leaf dry mass per unit area (LMA) is a leaf morphological trait related to shape widely used as a proxy for sclerophylly (see Alonso-Forn et al. 2020 and references therein).LMA commonly varies along environmental gradients (Niinemets 2001;Poorter et al. 2009).It is widely recognized that an increase in LMA is a general response to environmental stress (Onoda et al. 2011;Alonso-Forn et al. 2020), such as drought (Turner 1994;Reich et al. 1999;Jordan et al. 2005) or low winter temperatures (Ogaya and Peñuelas 2007;González-Zurdo et al. 2016;Niinemets 2016). Recently, Sancho-Knapik et al. (2021) showed that LMA in deciduous oak species may be modulated by aridity, whereas diversification of LMA within evergreen oaks seems modulated by both aridity and winter cold.
The functions of leaf pubescence are more ambiguous.Functions that have been demonstrated include providing resistance against biotic or abiotic damages, absorbing nutrients and water, and protecting against excessive solar radiation (Ehleringer et al. 1976;Sakai et al. 1980;Ehleringer 1981;Werker 2000;Morales et al. 2002;Agrawal et al. 2009;Fernández et al. 2014).Adaxial trichomes may increase leaf reflectance, which could be advantageous in dry and hot climates by reducing heat load and leaf temperature without increasing water loss through transpiration (Ehleringer and Björkman 1978;Ehleringer and Mooney 1978;Pérez-Estrada et al. 2000;Abdulrahaman and Oladele 2011;Bickford 2016).Additionally, trichomes have been explored for their role in repelling water from the leaf surface to keep stomata dry and ensure free gas exchange (Brewer and Smith 1997;Pierce et al. 2001).However, some studies have found conflicting results regarding the role of trichomes in regulating leaf transpiration and coping with aridity, indicating that further research is needed to understand better their functional significance and possible contextdependence of different roles (Aronne and DeMicco 2001;Haworth and McElwain 2008).
There are few studies within Fagaceae genera that analyze the relationships between morphological traits and their response to climatic variables in a broad set of species and climatic conditions.For example, Chen and Kohyama (2022) studied in 91 Lithocarpus species how fruit and vegetative traits correlate with temperature and wetness, however, Lithocarpus is a genus with little morphological and climatic variation compared to other Fagaceae genera.Most of the research focuses on other proxies, such as how a reduced set of species, or even a single species, responds to a climatic gradient or to climate change, typically using the most common species (e.g.Fagus sylvatica or Castanea sativa).Quercus (oaks) is the genus with more research in this matter since it is an excellent system for investigating how co-variation in leaf traits reflects complex adaptations to varying environments (Ramírez-Valiente et al. 2020;Cavender-Bares 2019).Oaks show considerable variation amongspecies in leaf morphological traits, climatic niches and geographical range (Gil-Pelegrín et al. 2017;Cavender-Bares 2019;Jiang et al. 2019;Sancho-Knapik et al. 2021).Although some studies have explored the relationships between climate and morphology within oaks (Kaproth et al. 2023;Ramírez-Valiente et al. 2020) they are limited by the number of species or the biogeographical range.To our knowledge, there is not any study taking into account a great number of species, representing almost every biome occupied by oak species and involving species from each infrageneric clade.
The study aims to assess in a wide range of oak species ( 141) to what extent leaf morphological traits (leaf area, leaf mass per area, length-width ratio, index of lobation, leaf shape, leaf apex, leaf margin, pubescence and leaf habit) are associated with the prevailing macroclimate within the species distribution.The study also includes clade as a covariate to investigate the partial effects of clade and climate in the evolution of oak leaf syndromes.These analyses help illuminate the relative importance of phylogenetic niche conservatism and adaptation in the broad evolution of oak leaf syndromes.

Plant material
One hundred and forty-one oak species (Quercus spp.) (Martín-Sánchez et al. 2024, Table S1) were sampled during midsummer from living collections in the Jardín Botánico de Iturrarán (43°130N, 02°010W, 70 m asl, Gipuzkoa, Spain, www.iturr aran.org).Using of the botanical garden as a quasi-common garden in our study enables interspecific comparisons of traits that result from long-term adaptation to the species' native climate without conflating that adaptive history with the short-term plastic response of each tree to the climate where it occurs naturally.Note that some variation we observe is plastic in response to site conditions within the botanic garden, but this variation is essentially irrelevant to our question about among-biome adaptation.Ten fully developed leaves per species (3-4 leaves per tree × 3 different trees) were collected from mature trees (15-25 years old) in midsummer.Leaves were sampled from south-exposed branches during early morning (07:00-09:00 h, solar time) and they were stored in sealed plastic bags and carried to the laboratory in portable coolers to measure both quantitative and qualitative morphological traits.

Morphological traits
In this study, we considered five qualitative variables: overall shape, apex shape, margin, pubescence and leaf habit.We also measured four quantitative variables: leaf area (LA), leaf mass per area (LMA), length-width ratio (LW) and index of lobation (ILB).Leaf shape, apex and margin were classified according to Ellis et al. (2009) as one of six shapes (obovate, elliptic, oblong, ovate, lanceolate or circular), one of four apex shapes (rounded, acute, acuminate or straight), and one of six margin types (entire, serrate, dentate, crenate, undulate or spinose).Although quantitative variables are always preferred, we decided to include qualitative traits to support the continuous traits.The variability in nominal margin is gathered by the continuous ILB.In the same way, LW strongly overlaps with the shape.However, nominal traits may explain nuances that are not achieved with these continuous traits.Regarding pubescence, we followed Gil-Pelegrín et al. (2017), classifying species into three categories: glabrous (or with scattered trichomes), pubescent, and densely pubescent, according to the botanical description of the species in their respective floras (Flora of North America, Nixon 1993; Flora of China, Huang et al. 1999;and Flora Europaea, Schwarz 1964).The same floras were taken into account to classify species as deciduous or evergreen.

Climatic and geographic data
Climate analysis for each studied species were carried out following a similar procedure used by Gil-Pelegrín et al. (2017).We first obtained geographical distribution coordinates for each species individually from GBIF (https:// www.gbif.org/) taking into account their distribution ranges (when available, see Euforgen and US Tree Atlas databases) or the geographical descriptions in the beforementioned floras to crop the areas from where points were downloaded.Then, the downloaded datasheets were thinned using SDMtune R package so that redundant points within a square kilometre were removed (Vignali et al. 2022).The resultant information (see Martín-Sánchez et al. 2024, Table S2) was used to obtain the values of each climatic variable from the WorldClim version 2.1 database (Fick and Hijmans 2017.WorldClim 2. https:// www.world clim.org/) for each presence point using Raster R package.We calculated the annual precipitation range (bio20) as BIO16 minus BIO17 (see Martín-Sánchez et al. 2024, Table S3 for the meaning of each bioclimatic variable).Points were also overlaid on a climatic Köppen map (Köppen 1936) to extract their Köppen climate (Martín-Sánchez et al. 2024, Table S2).S4).Species were assigned to each macroclimate according to where most of their occurrences fall compared with the rest of the other macroclimates (Appendix Fig. 8).

Statistical analysis
ANOVAs were performed to quantify the variance explained by macroclimate, subgenus, section-nested into subgenus-on each quantitative variable (LA, LMA, LW and ILB) (Table 1).We also performed ANOVAs to compare quantitative traits among macroclimates.Multiple comparisons were carried out among macroclimates for LA, LMA, LW and ILB using the post-hoc Tukey's honestly significant difference test.Heteroscedasticity and normality of the data were checked throughout analyses of the residuals, and Kolmogorov-Smirnov tests were carried out using the DHARMa R package (Hartig 2022).In addition, contingency tables based on pairwise χ 2 tests-corrected by Holm-Bonferroni method for multiple comparisons-were used to test the correlation between macroclimate and the qualitative traits: margin, apex, shape, pubescence and leaf habit.
Factor analysis of mixed data (FAMD) was carried out to summarize the multivariate relationships among First dimension scores (PC1 Morphology) were treated as a generalized leaf morphological response variable influenced by all four quantitative traits together.An additional ANOVA was performed to quantify the variance explained by macroclimate, subgenus, and section when all traits are taken into account altogether (Table 1).Principal components analysis (PCA) was performed with climatic variables extracted from Worldclim database, with a previous selection of those that accounted for more variance in the first two principal components via the FactoExtra R package (Kassambara and Mundt 2020).Every data included in both FAMD and climatic PCA represent the mean value for each variable and species.
Finally, we did cross-correlations between the first and second dimensions of both FADM and climatic PCA to check if there is any significance between both analyses and to explore the existence of relationships between morphological traits and climatic variables.We first tested if there was a phylogenetic signal performing phylogenetic generalized least squares (PGLS) regressions using the most complete published phylogenetic tree for Quercus (Hipp et al. 2020).Section and subgenus names correspond to named clades in Hipp et al. (2020).Ninetynine species intersected between phylogenetic and trait datasets.We first performed a PGLS with all species.The results do not allow us to refuse the lack of phylogenetic signal.However, since the slopes of both regressions: ordinary least squares (OLS) and PGLS, are minimally different (Appendix Fig. 9) and we want to include all taxa, we decided to present the results of the OLS regression.We also performed additional PGLS within the different clades: both subgenera and the two sections with more species, i.e.Quercus and Lobatae (Martín-Sánchez et al. 2024, Table S5).

Results
ANOVAs reveal that subgenus and macroclimate significantly (P < 0.001) account for 6.9% and 35.4% of the morphological variance when traits are analyzed altogether (Table 1).The taxonomic section is a non-significant predictor of 3.2% of the variance (P = 0.18).54.5% of variance was unexplained.By contrast, when traits are analysed separately, we find that macroclimate explains more variance for LA (29%), LW (20.3%) and ILB (14.4%) than any other factor.However, in the case of LMA, not only macroclimate is significant, but also subgenus and section, with this last one explaining the highest percentage of variance (16.4%) (Table 1).
Contingency tables show significant associations between the five qualitative morphological traits and macroclimate.In terms of shape, Obovate leaves relate to Cold Temperate (P = 0.003) and Warm Temperate (P < 0.001) species; besides, elliptic leaves to Table 1 Summary of percentage of variance explained by ANOVAs performed for each variable separately: Leaf area (LA), lengthwidth ratio (LW), index of lobation (ILB), leaf mass per area (LMA) and altogether (PC1 scores from FAMD).Each row represents a separate ANOVA.All models include all factors: macroclimate, subgenus and section nested in subgenus.Number of asterisks represents the level of significance (P < 0.001 ***, P < 0.01 **, P < 0.05 *).The proportion of variance is calculated dividing the sum squares for each factor by the total sum of squares Winter-dry temperate to Tropical species (P < 0.001) (Fig. 3a).For apex, acute (P < 0.001) and acuminate (P = 0.01) apices are linked to Winter-dry temperate to Tropical species (Fig. 3b).Lobated margins are associated with Cold Temperate (P = 0.014) and Warm Temperate (P < 0.001) species, entire margins with Warm Temperate (P = 0.014) and Winter-dry temperate to Tropical species (P < 0.014) and serrate margins with Winter-dry temperate to Tropical species (P < 0.001) (Fig. 3c).Glabrous leaves are associated with Winterdry temperate to Tropical species (P < 0.001) (Fig. 3d).Finally, deciduousness is significantly correlated with Cold Temperate (P = 0.008) and Warm Temperate species (P < 0.001) (Fig. 3).FAMD on the combined morphological groups species together by macroclimate, but with strong overlap (Fig. 5a).Traits cluster loosely by climate into three groups: (1) obovate, lobated leaves with large areas corresponding mainly to species from Cold Temperate climates, with some Warm Temperate and deciduous Mediterranean species; (2) elongate leaves with acuminate apices, lacking abaxial trichomes tend to be associated with Winter-dry temperate to Tropical and Warm Temperate species; and (3) oblong and undulate leaves with dentate margins, rounded apices, a dense tomentum abaxially and a high LMA are concentrated in Arid and evergreen Mediterranean species (Fig. 5a).
Climatic PCA clusters species into well-defined groups, with the most overlap between Arid and Mediterranean species (Fig. 5b).Cold Temperate species, are defined by a high temperature annual range (BIO3, BIO7), with very cold winters (BIO6).By contrast, Winter-dry temperate to Tropical species are characterized by warmer winters (BIO6) and more temperature-stable climates (BIO3 and BIO4) but pronounced seasonality in precipitation regimes (BIO15, annual precipitation range), with high rainfall during the wet season (BIO16, BIO18, annual precipitation range).Warm temperate species do not suffer a great reduction in precipitation during their driest month.Arid and Mediterranean species occupy climates with a very dry season during the warmest period of the year (BIO18, annual precipitation range) and low annual precipitation (BIO12).This relationship between climate and morphology is supported as well by the OLS regression between morphology PCA axis 1 (PC1-Morph) and climate PCA axis 1 (PC1-Clim), which is significant (r = 0.543, P < 0.001) (Fig. 6).The same correlation performing a PGLS shows a similar slope, significance (P < 0.001) and a low phylogenetic signal (λ = 0.28) (Appendix Fig. 9).The additional PGLS within the different clades: subgenus Cerris, subgenus Quercus, section Lobatae and section Quercus (clades named as in Hipp et al. 2020) revealed that we cannot claim that both subgenera are  S5).

Individual leaf traits correlate with climate
Leaf area exhibits high evolutionary lability in some clades, with high potential to evolve in response to climate variation (Jones et al. 2009;Peppe et al. 2011).Our work shows a strong negative correlation between leaf size, growing season temperature, and aridity.Thus, the smallest leaves are found in Arid and Mediterranean climates, the largest ones in Cold Temperate climates (Fig. 2a).This may be a result of oaks reducing LA to decrease their boundary layer, enhancing heat loss (Vogel 1970;2009;Givnish 1987).A reduced LA also contributes to a higher leaf-specific conductivity (LSC) of the shoot which is advantageous in dry habitats (Peguero-Pina et al. 2014, 2016).Conversely, larger leaves confer to species a greater photosynthetic rate when compared with smaller leaves.However, many small-leaf species usually compensate for this reduction of photosynthetic tissue with a higher ability to photosynthesize throughout a higher carboxylation rate, mesophyll and stomatal conductance (Peguero-Pina et al. 2016).Of leaf shape traits, lobation, length-width ratio, and leaf apex exhibit the strongest correlation with climate.Lobing assessed qualitatively and higher ILB, both  were associated significantly with temperate climates (Fig. 2d, Fig. 3c).Lobation is also exhibited in sub-Mediterranean species, which occupy a climate transitional between Mediterranean and temperate (Abadía et al. 1996;Himrane et al. 2004;Gil-Pelegrín et al. 2017).
Cold Temperate, sub-Mediterranean and some of the Warm Temperate (primarily American) species gather under the same group in the FAMD (Fig. 5a), which is defined by four traits: high ILB, high LA and obovate and lobated shapes.Leaf lobation may enhance hydraulic efficiency by reducing the ratio of mesophyll to conductive tissue (Sisó et al. 2001;Sack and Tyree 2005) and decreasing the boundary layer to speed leaf cooling (Vogel 1970(Vogel , 2009;;Givnish 1987).Both potential benefits of lobation are especially relevant in sub-Mediterranean species, which experience very high, dry temperatures during summer compared to temperate species.The sub-Mediterranean species in our study also compensate for the reduction of photosynthetic tissue in their relatively small and / or lobated leaves by increasing net CO 2 assimilation per unit leaf area, as Peguero-Pina et al. (2016) showed when comparing the sub-Mediterranean Q. faginea with the cold temperate Q. robur.In temperate environments, by contrast, leaf lobation may play a role in bud packing, which is crucial for competitive growth in winter-freezing environments (Edwards et al. 2016).
By contrast, Winter-dry temperate to Tropical and some Warm Temperate species exhibit a high LW (Fig. 2c), besides a significant percentage of lanceolate shapes (Fig. 3a) and acuminate apices (Fig. 3b).These three traits associate in the FAMD with Winterdry temperate to Tropical and Warm Temperate species (Fig. 5a).Together with a tilted leaf position, these traits have been hypothesized to limit rainfall damage (Wang et al. 2020).High rains in the tropics and Warm Temperate regions can slow photosynthesis (Ishibashi and Terashima 1995), leach minerals from the leaves (Minorsky 2019), damage leaves mechanically (Wang et al. 2020), and promote pathogens (Burd 2007).The morphological mechanisms we document, which evacuate water from the leaf surface, are particularly needed in the Warm Temperate and Winter-dry temperate to Tropical climates where these traits dominate.These climates, in the Quercus distribution, are most characterized by high precipitation (Fig. 5b).
LMA follows a completely different pattern in our study.Although it is significantly modulated by climate (Table 1), the major clade in which a species occurs plays a more significant role in explaining differences in LMA among species (Table 1) (cf.Sancho-Knapik et al. 2021).The effect of clade at section level is likely explained by the LMA-homogeneous sections, i.e., Ilex, Protobalanus and Virentes.After accounting for the clade, Arid and Mediterranean species tend to have a higher LMA than the rest of the species, although not significant (Fig. 2b).This relatively weak effect of climate is particularly remarkable given the plasticity that LMA presents in response to several environmental factors, such as drought (Niinemets 2001), exposure to light (Poorter et al. 2006) or herbivory (Mediavilla et al. 2018;Onoda et al. 2011).For example, Sancho-Knapik et al. (2021) demonstrated that oak species with long cold periods and/or humid climates have lower LMA.
Despite the poor correlation between climate and LMA, we find a significant correlation between sclerophylly-which is strongly correlated with LMA-and continuous climatic variables (Alonso-Forn et al. 2020).Species with low LMA are usually deciduous and malacophyllous.In contrast, species with a high LMA tend to be evergreen and sclerophyllous (Mediavilla et al. 2008) (Appendix Fig. 10).Thus, Arid, evergreen Mediterranean, most of Winter-dry temperate to Tropical and some Warm Temperate species are located at low Dim1 values in the FAMD (Fig. 5a).By contrast, deciduous Mediterranean, Cold Temperate and the rest of Warm Temperate species are located at high Dim1 values.
The increment in LMA, both via increasing leaf density or leaf thickness, is a major resource investment for the plant.Species which make a higher investment per leaf must keep their photosynthetic tissues active for longer to recover the resources expended according to the worldwide leaf economics spectrum (Wright et al. 2004).However, there are several factors within an individual to take into account in the construction costs and payback time, such as light exposure of the leaves, chemical compounds and leaf habit.Poorter et al. (2006) showed how shade leaves present longer payback times than light leaves, in spite of the fact that these long times can be partially compensated by decreasing the LMA.Differences in payback time can be found between deciduous and evergreen species too (Villar and Merino 2001).Thus, species with a higher LMA-correlated with sclerophylly and longer payback times-tend to be evergreen.Winter-dry temperate to Tropical and some Warm Temperate species can develop sclerophyllous leaves increasing their LMA more than temperate deciduous species but less than Arid and Mediterranean (Alonso-Forn et al. 2020), as the high contribution of LMA to Dim2 of FAMD reflects (Appendix Fig. 11b).This might be related to the two kinds of sclerophyllous leaves: pycnophylls and pachyphylls (Grubb 1986).Differences between them are found in the mesophyll since sclerophylly is mainly defined by an investment to reinforce the epidermis, which is shared by both types.Pycnophyll leaves present a very compact mesophyll, usually small cells with little space among them, i.e., an increase in LMA via investment in leaf density.By contrast, pachyphyll leaves are characterized by a larger, laxer and spongier mesophyll, i.e., an increase in LMA via leaf thickness (Grubb 1986).Thereby, Winter-dry temperate to Tropical and Warm Temperate species might be mainly pachyphylls, exhibiting thick leaves with low density; meanwhile, Arid and evergreen Mediterranean might be pycnophylls.
Leaf water relations relations are also strongly influenced by pubescence, which appears in our study to vary independently of macroclimate, except glabrous leaves associated with Winter-dry temperate to tropical climates (P < 0.001).Previous studies have suggested that abaxial pubescence in oaks could be a xeromorphic trait in response to drier climates (Fig. 3d) (Hardin 1979).This trend is observed in a major presence of pubescent species in Arid and Mediterranean species, although not significant (P = 0.41).(Fig. 3d).The functional roles of trichomes in the genus may be broader, including protection against a range of biotic and abiotic factors (Zvereva et al. 1998;Sato and Kudoh 2015;Ehleringer et al. 1976;Coltharp et al. 2021).A trichome layer may act as a protection against some environmental factors, especially excessive solar radiation, by modifying the leaf 's optical properties, lowering heat load, and reducing leaf temperature and transpiration (Ehleringer et al. 1976;Ehleringer 1981;Sandquist and Ehleringer 1997;Pérez-Estrada et al. 2000;Abdulrahaman and Oladele 2011;Bickforf 2016).Nonetheless, the hydrophobic nature of some trichomes may also be advantageous under wet climates, since they can keep the stomata dry and ensure the free exchange of CO 2 (Brewer and Smith 1997;Pierce et al. 2001).However, trichomes do not necessarily have to be hydrophobic, as Fernández et al. (2014) demonstrated in Quercus ilex subsp.rotundifolia, whose adaxial trichomes seem to help leaves to rehydrate during summer rainfall events.This double behaviour of trichomes and the lack of significance in our results indicate that their role as xeromorphic trait cannot be taken for granted, and further research to clarify the role of pubescence is needed.

Multivariate leaf syndromes associate with macroclimate
Our study demonstrates that the world's oak species assemble into well-defined climatic groups that correlate with loose but still identifiable leaf morphological syndromes.The Arid and Mediterranean species are an exception, as they mainly overlap in climatic space (Fig. 5b).Because species were classified into their macroclimate according to their presence within a specific Köppen climate, the climatic separation is a result of the Köppen system itself: PC1 and PC2 together account for 82% of the explained variance, which suggests that Köppen climates are well defined by WorldClim climatic variables.
A more interesting finding is the significant correlation of climatic PC1 with morphological PC1 (Fig. 6).Morphological PC1 is explained most strongly by lobation (ILB and lobated shapes), LA and Obovate shapes (Appendix Fig. 11a).In contrast, climatic variables with the highest scores are the difference in precipitation between the wettest month and the driest one (annual precipitation range), temperature annual range (BIO7), temperature seasonality (BIO4) and precipitation of the wettest quarter (BIO16) (Appendix Fig. 12a).Climatic PC1 thus separates most strongly between Winter-dry temperate to Tropical (with negative values) and the rest of the species (with positive ones) (Fig. 5b).The segregation of species along this climatic gradient represents both an ecological split and a phylogenetic split replicated twice in the Americas (Hipp et al. 2018) and between the primarily Mediterranean to temperate sections Cerris and Ilex of Eurasia versus section Cyclobalanopsis of the southeast Asia's broad-leaved evergreen forests (Deng et al. 2018).
Mediterranean, Arid, and temperate species live in climates with a marked thermal seasonality.The stressful dry or cold periods that mark these climates are much reduced in Winter-dry temperate to Tropical climates.Winter-dry temperate to Tropical, Arid and Mediterranean species-three macroclimates with higher temperatures, on average, than temperate climates-are represented with negative values in PC1 of FAMD (Fig. 4a).Simultaneously, temperate and sub-Mediterranean species appear together and are characterized by a high lobation, and large, obovate leaves (Fig. 5a), which are the four variables with highest loadings in PC1 (Appendix Fig. 11a).Thus, we see once again a correlation between area and lobation with temperature.
The significant regression between climate and morphology PC2 (Appendix Fig. 13) shows the correlation of glabrous, slender leaves with drip-tips with higherprecipitation climates and pubescent rounded leaves with dry climates (Appendix Fig. 11b and 12b).As morphological and climatic PC2 represents a water-availability gradient, Mediterranean, Arid, and Winter-dry temperate to Tropical with monsoon-dry season cluster in the FAMD separate from most Warm Temperate and the rest of the Winter-dry temperate to Tropical species (Fig. 5a).
Cold Temperate species are generally characterized by large, obovate, deciduous leaves with complex-generally lobated-margins and low LMA (Fig. 7).Winter-dry temperate to Tropical species are defined by long, lanceolate and acuminate leaves, generally glabrous, evergreen or tardily-deciduous and sclerophyllous, likely pachypylltype (anatomical studies in further researches are needed) (Fig. 7).Warm Temperate species tend to share or present intermediary traits between Cold Temperate and Winter-dry temperate to Tropical species (Fig. 7).Arid and Mediterranean species exhibit the smallest leaves, short, generally with a LW ratio near 1, and entire or relatively simple margins.They are usually evergreen or semi-evergreen, with the highest values of LMA probably due to their hypothetical sclerophyll pycnophyll-type condition (also further researches are needed) (Fig. 7).Finally, sub-Mediterranean species (e.g.Q. pyrenaica, Q. frainetto or Q. cerris) are characterized by larger leaves than their genuine Mediterranean relatives, with deeply lobated margins, usually obovate, deciduous and pubescent (Fig. 7).Overall, our work shows that we can relate leaf syndromes in oaks to broad-scale patterns in global climate.
Finally, although climate explains most of the variance found in oak leaf morphology, clade seems to have a minor but still significative effect since we cannot reject the presence of a phylogenetic signal in these morphological traits, especially in LMA.This phylogenetic signal could be appearing due to some morphologically homogeneous clades, e.g.Sect.Protobalanus, Virentes or Ilex.

Conclusion
In this study we analysed 141 Quercus species across a broad range of northern hemisphere biomes, ca.25% of the species in the genus.Although we cannot refuse the presence of a phylogenetic signal, our results strongly suggest that leaf morphological traits of oaks are adapted to their macroclimates, as previous work has demonstrated for leaf anatomy (Sancho-Knapik et al. 2021;Alonso-Forn et al. 2023).Our study demonstrates that even in a common garden, growing under the same climatic conditions, a wide sample of the world's oak species develop a set of morphological traits with different ecological roles that conclude into recognisable leaf syndromes according to the climate that these species occupy.Except for Mediterranean species that, on the one hand, tend to develop a common leaf syndrome with Arid species.On the other hand, another group composed by sub-Mediterranean species gathers into a temperate syndrome.
Even though climate explains most of the variance found in leaf morphology, the partial effects of clade cannot be disregarded, as we see within the LMA, which resulted to be an outlier trait in this study.LMA is still explained significantly by climate; however, section explains twice variance as much as the climate for this trait.
Finally, pubescence did not show significant enough results to be considered a xeromorphic trait, as many authors have proposed.
Köppen climates inhabited by Quercus were grouped into five macroclimates: Arid, Mediterranean, Winter-dry temperate to Tropical, Warm Temperate and Cold Temperate, similar to the schema used by Gil-Pelegrín et al. (2017) (Fig. 1).Arid macroclimate (ARID) is defined by climates where annual precipitation is lower than evapotranspiration (BSk and BSh primarily, with some peripheral populations in BWk or BWh).Mediterranean macroclimate (MED) comprises climates with the dry period in summer (Csa, Csb, Csc, Dsa and Dsb).Winter-dry temperate to Tropical macroclimate (WdT-TROP) includes warm and temperate climates with a dry season other than summer (Aw, Am, Cwa and Cwb).Warm Temperate macroclimate (WARM_T) comprises Cfa Köppen climates, since it represents a climate widely colonized by Quercus species, especially in southeastern regions of the United States and China.The Cold Temperate (COLD_T) macroclimate comprises Dfa, Dfb, Dwa, Dwb, Cfb, and Cfc primarily, with some peripheral populations in Dsc, Dsd, Dwc, Dwd, Dfc or Dfd.Several species of oaks, especially some Cyclobalanopsis species, live in an equatorial climate (Af ), which was not treated in this study as these last species are not considered in this study (for a more detailed definition of the Köppen climate categories see Martín-Sánchez et al. 2024, Table

Fig. 1
Fig. 1 Map of the five macroclimates considered in this study: Arid in yellow (BSk, BSh, according to Köppen Classification), Mediterranean in red (Csa, Csb, Csc, Dsa, Dsb), Winter-dry to Tropical in cyan (Aw, Am, Cwa, Cwb), Warm Temperate in light green (Cfa) and Cold Temperate in dark blue (Cfb, Cfc, Dfa, Dfb, Dwa, Dwb).Colored areas represent the regions where Quercus species are present.The Equatorial climate (not considered in this study) is represented by Tropical Rainforest Climate (Af ) and coloured with dark green

Fig. 13
Fig. 13 Correlations between PC2 scores resulting from FAMD (PC2-Morph) and climatic PCA (PC2-Clim).Individual species are represented by dots and their colour represent the macroclimate which they have been classified: Arid (yellow), Mediterranean (red), Winter-dry temperate to Tropical (cyan), Warm Temperate (light green) and Cold Temperate (dark blue).A linear model was adjusted and showed a significant p-value (P < 0.001)