Plasticity in functional traits in the context of climate change: a case study of the subalpine forb Boechera stricta (Brassicaceae)

Environmental variation often induces shifts in functional traits, yet we know little about whether plasticity will reduce extinction risks under climate change. As climate change proceeds, phenotypic plasticity could enable species with limited dispersal capacity to persist in situ, and migrating populations of other species to establish in new sites at higher elevations or latitudes. Alternatively, climate change could induce maladaptive plasticity, reducing fitness, and potentially stalling adaptation and migration. Here, we quantified plasticity in life history, foliar morphology, and ecophysiology in Boechera stricta (Brassicaceae), a perennial forb native to the Rocky Mountains. In this region, warming winters are reducing snowpack and warming springs are advancing the timing of snow melt. We hypothesized that traits that were historically advantageous in hot and dry, low‐elevation locations will be favored at higher elevation sites due to climate change. To test this hypothesis, we quantified trait variation in natural populations across an elevational gradient. We then estimated plasticity and genetic variation in common gardens at two elevations. Finally, we tested whether climatic manipulations induce plasticity, with the prediction that plants exposed to early snow removal would resemble individuals from lower elevation populations. In natural populations, foliar morphology and ecophysiology varied with elevation in the predicted directions. In the common gardens, trait plasticity was generally concordant with phenotypic clines from the natural populations. Experimental snow removal advanced flowering phenology by 7 days, which is similar in magnitude to flowering time shifts over 2–3 decades of climate change. Therefore, snow manipulations in this system can be used to predict eco‐evolutionary responses to global change. Snow removal also altered foliar morphology, but in unexpected ways. Extensive plasticity could buffer against immediate fitness declines due to changing climates.     

Drought stress limits the geographic ranges of two tree species via different physiological mechanisms

Range shifts are among the most ubiquitous ecological responses to anthropogenic climate change and have large consequences for ecosystems. Unfortunately, the ecophysiological forces that constrain range boundaries are poorly understood, making it difficult to mechanistically project range shifts. To explore the physiological mechanisms by which drought stress controls dry range boundaries in trees, we quantified elevational variation in drought tolerance and in drought avoidance‐related functional traits of a widespread gymnosperm (ponderosa pine – Pinus ponderosa) and angiosperm (trembling aspen – Populus tremuloides) tree species in the southwestern USA. Specifically, we quantified tree‐to‐tree variation in growth, water stress (predawn and midday xylem tension), drought avoidance traits (branch conductivity, leaf/needle size, tree height, leaf area‐to‐sapwood area ratio), and drought tolerance traits (xylem resistance to embolism, hydraulic safety margin, wood density) at the range margins and range center of each species. Although water stress increased and growth declined strongly at lower range margins of both species, ponderosa pine and aspen showed contrasting patterns of clinal trait variation. Trembling aspen increased its drought tolerance at its dry range edge by growing stronger but more carbon dense branch and leaf tissues, implying an increased cost of growth at its range boundary. By contrast, ponderosa pine showed little elevational variation in drought‐related traits but avoided drought stress at low elevations by limiting transpiration through stomatal closure, such that its dry range boundary is associated with limited carbon assimilation even in average climatic conditions. Thus, the same climatic factor (drought) may drive range boundaries through different physiological mechanisms – a result that has important implications for process‐based modeling approaches to tree biogeography. Further, we show that comparing intraspecific patterns of trait variation across ranges, something rarely done in a range‐limit context, helps elucidate a mechanistic understanding of range constraints.     

Plant functional traits differ in adaptability and are predicted to be differentially affected by climate change

1. Climate change is testing the resilience of forests worldwide pushing physiological tolerance to climatic extremes. Plant functional traits have been shown to be adapted to climate and have evolved patterns of trait correlations (similar patterns of distribution) and coordinations (mechanistic trade‐off). We predicted that traits would differentiate between populations associated with climatic gradients, suggestive of adaptive variation, and correlated traits would adapt to future climate scenarios in similar ways.
2. We measured genetically determined trait variation and described patterns of correlation for seven traits: photochemical reflectance index (PRI), normalized difference vegetation index (NDVI), leaf size (LS), specific leaf area (SLA), δ¹³C (integrated water‐use efficiency, WUE), nitrogen concentration (N_CONC), and wood density (WD). All measures were conducted in an experimental plantation on 960 trees sourced from 12 populations of a key forest canopy species in southwestern Australia.
3. Significant differences were found between populations for all traits. Narrow‐sense heritability was significant for five traits (0.15–0.21), indicating that natural selection can drive differentiation; however, SLA (0.08) and PRI (0.11) were not significantly heritable. Generalized additive models predicted trait values across the landscape for current and future climatic conditions (>90% variance). The percent change differed markedly among traits between current and future predictions (differing as little as 1.5% (δ¹³C) or as much as 30% (PRI)). Some trait correlations were predicted to break down in the future (SLA:N_CONC, δ¹³C:PRI, and N_CONC:WD).
4. Synthesis: Our results suggest that traits have contrasting genotypic patterns and will be subjected to different climate selection pressures, which may lower the working optimum for functional traits. Further, traits are independently associated with different climate factors, indicating that some trait correlations may be disrupted in the future. Genetic constraints and trait correlations may limit the ability for functional traits to adapt to climate change.

     

Drought responses,phenotypic plasticity and survival of Mediterranean species in two different microclimatic sites

• Climate models predict a further drying of the Mediterranean summer. One way for plant species to persist during such climate changes is through acclimation. Here, we determine the extent to which trait plasticity in response to drought differs between species and between sites, and address the question whether there is a trade‐off between drought survival and phenotypic plasticity.
• Throughout the summer we measured physiological traits (photosynthesis – A_max, stomatal conductance – g_s, transpiration – E, leaf water potential – ψl) and structural traits (specific leaf area – SLA, leaf density – LD, leaf dry matter content – LDMC, leaf relative water content – LRWC) of leaves of eight woody species in two sites with slightly different microclimate (north‐ versus south‐facing slopes) in southern Spain. Plant recovery and survival was estimated after the summer drought period.
• We found high trait variability between species. In most variables, phenotypic plasticity was lower in the drier site. Phenotypic plasticity of SLA and LDMC correlated negatively with drought survival, which suggests a trade‐off between them. On the other hand, high phenotypic plasticity of SLA and LDMC was positively related to traits associated with rapid recovery and growth after the drought period.
• Although phenotypic plasticity is generally seen as favourable during stress conditions, here it seemed beneficial for favourable conditions. We propose that in environments with fluctuating drought periods there can be a trade‐off between drought survival and growth during favourable conditions. When climate become drier, species with high drought survival but low phenotypic plasticity might be selected for.

     

Characterizing the contribution of plasticity and genetic differentiation to community‐level trait responses to environmental change

The match between functional trait variation in communities and environmental gradients is maintained by three processes: phenotypic plasticity and genetic differentiation (intraspecific processes), and species turnover (interspecific). Recently, evidence has emerged suggesting that intraspecific variation might have a potentially large role in driving functional community composition and response to environmental change. However, empirical evidence quantifying the respective importance of phenotypic plasticity and genetic differentiation relative to species turnover is still lacking. We performed a reciprocal transplant experiment using a common herbaceous plant species (Oxalis montana) among low‐, mid‐, and high‐elevation sites to first quantify the contributions of plasticity and genetic differentiation in driving intraspecific variation in three traits: height, specific leaf area, and leaf area. We next compared the contributions of these intraspecific drivers of community trait–environment matching to that of species turnover, which had been previously assessed along the same elevational gradient. Plasticity was the dominant driver of intraspecific trait variation across elevation in all traits, with only a small contribution of genetic differentiation among populations. Local adaptation was not detected to a major extent along the gradient. Fitness components were greatest in O. montana plants with trait values closest to the local community‐weighted means, thus supporting the common assumption that community‐weighted mean trait values represent selective optima. Our results suggest that community‐level trait responses to ongoing climate change should be mostly mediated by species turnover, even at the small spatial scale of our study, with an especially small contribution of evolutionary adaptation within species.     

Plasticity of functional traits varies clinally along a rainfall gradient in Eucalyptus tricarpa

Widespread species often occur across a range of climatic conditions, through a combination of local genetic adaptations and phenotypic plasticity. Species with greater phenotypic plasticity are likely to be better positioned to cope with rapid anthropogenic climate changes, while those displaying strong local adaptations might benefit from translocations to assist the movement of adaptive genes as the climate changes. Eucalyptus tricarpa occurs across a climatic gradient in south‐eastern Australia, a region of increasing aridity, and we hypothesized that this species would display local adaptation to climate. We measured morphological and physiological traits reflecting climate responses in nine provenances from sites of 460 to 1040 mm annual rainfall, in their natural habitat and in common gardens near each end of the gradient. Local adaptation was evident in functional traits and differential growth rates in the common gardens. Some traits displayed complex combinations of plasticity and genetic divergence among provenances, including clinal variation in plasticity itself. Provenances from drier locations were more plastic in leaf thickness, whereas leaf size was more plastic in provenances from higher rainfall locations. Leaf density and stomatal physiology (as indicated by δ¹³C and δ¹⁸O) were highly and uniformly plastic. In addition to variation in mean trait values, genetic variation in trait plasticity may play a role in climate adaptation.     

Intraspecific trait variation across elevation predicts a widespread tree species' climate niche and range limits

Global change is widely altering environmental conditions which makes accurately predicting species range limits across natural landscapes critical for conservation and management decisions. If climate pressures along elevation gradients influence the distribution of phenotypic and genetic variation of plant functional traits, then such trait variation may be informative of the selective mechanisms and adaptations that help define climatic niche limits. Using extensive field surveys along 16 elevation transects and a large common garden experiment, we tested whether functional trait variation could predict the climatic niche of a widespread tree species (Populus angustifolia) with a double quantile regression approach. We show that intraspecific variation in plant size, growth, and leaf morphology corresponds with the species' total climate range and certain climatic limits related to temperature and moisture extremes. Moreover, we find evidence of genetic clines and phenotypic plasticity at environmental boundaries, which we use to create geographic predictions of trait variation and maximum values due to climatic constraints across the western US. Overall, our findings show the utility of double quantile regressions for connecting species distributions and climate gradients through trait‐based mechanisms. We highlight how new approaches like ours that incorporate genetic variation in functional traits and their response to climate gradients will lead to a better understanding of plant distributions as well as identifying populations anticipated to be maladapted to future environments.     

Environmental Heterogeneity and Population Differentiation in Plasticity to Drought in <Emphasis Type="Italic">Convolvulus Chilensis</Emphasis> (Convolvulaceae)

Plant populations may show differentiation in phenotypic plasticity, and theory predicts that greater levels of environmental heterogeneity should select for higher magnitudes of phenotypic plasticity. We evaluated phenotypic responses to reduced soil moisture in plants of Convolvulus chilensis grown in a greenhouse from seeds collected in three natural populations that differ in environmental heterogeneity (precipitation regime). Among several morphological and ecophysiological traits evaluated, only four traits showed differentiation among populations in plasticity to soil moisture: leaf area, leaf shape, leaf area ratio (LAR), and foliar trichome density. In all of these traits plasticity to drought was greatest in plants from the population with the highest interannual variation in precipitation. We further tested the adaptive nature of these plastic responses by evaluating the relationship between phenotypic traits and total biomass, as a proxy for plant fitness, in the low water environment. Foliar trichome density appears to be the only trait that shows adaptive patterns of plasticity to drought. Plants from populations showing plasticity had higher trichome density when growing in soils with reduced moisture, and foliar trichome density was positively associated with total biomass. Co-ordinating editor: F. Stuefer     

Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak Quercus oleoides

The impacts of drought are expanding worldwide as a consequence of climate change. However, there is still little knowledge of how species respond to long‐term selection in seasonally dry ecosystems. In this study, we used Q_ST‐F_ST comparisons to investigate (i) the role of natural selection on population genetic differentiation for a set of functional traits related to drought resistance in the seasonally dry tropical oak Quercus oleoides and (ii) the influence of water availability at the site of population origin and in experimental treatments on patterns of trait divergence. We conducted a thorough phenotypic characterization of 1912 seedlings from ten populations growing in field and greenhouse common gardens under replicated watering treatments. We also genotyped 218 individuals from the same set of populations using eleven nuclear microsatellites. Q_ST distributions for leaf lamina area, specific leaf area, leaf thickness and stomatal pore index were higher than F_ST distribution. Results were consistent across growth environments. Genetic differentiation among populations for these functional traits was associated with the index of moisture at the origin of the populations. Together, our results suggest that drought is an important selective agent for Q. oleoides and that differences in length and severity of the dry season have driven the evolution of genetic differences in functional traits.     

Tolerance to multiple climate stressors: a case study of Douglas‐fir drought and cold hardiness

相似文献   

Drought affects the coordination of belowground and aboveground resource‐related traits in Solidago canadensis in China

Quantifying patterns of variation and coordination of plant functional traits can help to understand the mechanisms underlying both invasiveness and adaptation of plants. Little is known about the coordinated variations of performance and functional traits of different organs in invasive plants, especially in response to their adaptation to environmental stressors. To identify the responses of the invasive species Solidago canadensis to drought, 180 individuals were randomly collected from 15 populations and 212 ramets were replanted in a greenhouse to investigate both the response and coordination between root and leaf functional traits. Drought significantly decreased plant growth and most of the root and leaf functional traits, that is, root length, surface area, volume and leaf size, number, and mass fraction, except for the root length ratio and root mass fraction. Phenotypic plasticity was higher in root traits than in leaf traits in response to drought, and populations did not differ significantly. The plasticity of most root functional traits, that is, root length (RL), root surface area (RSA), root volume (RV), and root mass fraction (RMF), were significantly positively correlated with biomass between control and drought. However, the opposite was found for leaf functional traits, that is, specific leaf area (SLA), leaf area ratio (LAR), and leaf mass fraction (LMF). Drought enhanced the relationship between root and leaf, that is, 26 pairwise root–leaf traits were significantly correlated under drought, while only 15 pairwise root–leaf traits were significantly correlated under control conditions. Significant correlations were found between biomass and all measured functional traits except for leaf size. RV, root length ratio, RMF, total area of leaves, and LMF responded differently to water availability. These responses enable S. canadensis to cope with drought conditions and may help to explain the reason of the vast ecological amplitude of this species.     

Tree hydraulic traits are coordinated and strongly linked to climate‐of‐origin across a rainfall gradient

Plant hydraulic traits capture the impacts of drought stress on plant function, yet vegetation models lack sufficient information regarding trait coordination and variation with climate‐of‐origin across species. Here, we investigated key hydraulic and carbon economy traits of 12 woody species in Australia from a broad climatic gradient, with the aim of identifying the coordination among these traits and the role of climate in shaping cross‐species trait variation. The influence of environmental variation was minimized by a common garden approach, allowing us to factor out the influence of environment on phenotypic variation across species. We found that hydraulic traits (leaf turgor loss point, stomatal sensitivity to drought [P_gs], xylem vulnerability to cavitation [P_x], and branch capacitance [C_branch]) were highly coordinated across species and strongly related to rainfall and aridity in the species native distributional range. In addition, trade‐offs between drought tolerance and plant growth rate were observed across species. Collectively, these results provide critical insight into the coordination among hydraulic traits in modulating drought adaptation and will significantly advance our ability to predict drought vulnerability in these dominant trees species.     

Variation in and adaptive plasticity of flower size and drought-coping traits

Although phenotypic plasticity of morphological and physiological traits in response to drought could be adaptive, there have been relatively few tests of plasticity variation or of adaptive plasticity in drought-coping traits across populations with different moisture availabilities. We measured floral size, vegetative size, and physiological traits in four field populations of Leptosiphon androsaceus (Polemoniaceae) that were distributed across a rainfall gradient in California, USA. Measurements were made over 5 years that varied in precipitation. We also conducted a growth chamber experiment in which half-sibs from three populations were divided equally among a well-watered and a drought treatment. We tested for selection on traits in each of the watering treatments, and evaluated whether traits exhibited plasticity. In the field, plant traits exhibited substantial variation across populations and years. Flower size, leaf size, and water-use efficiency (WUE) were generally higher for populations that received greater average rainfall. However, in dry years, we observed a decrease in flower and leaf size, but an increase in WUE across the populations. In the growth chamber experiment, leaf and physiological traits exhibited plasticity, with smaller leaves and higher WUE found in the drought, as compared to the well-watered treatment. Only specific leaf area exhibited differentiation in plasticity among populations. Although there was no observed plasticity in floral size, selection favored smaller flowers in the drought treatment and larger flowers in the well-watered treatment. Our results suggest that moisture availability has led to trait variation in L. androsaceus via a combination of selection and phenotypic plasticity.     

Experimental migration upward in elevation is associated with strong selection on life history traits

One of the strongest biological impacts of climate change has been the movement of species poleward and upward in elevation. Yet, what is not clear is the extent to which the spatial distribution of locally adapted lineages and ecologically important traits may also shift with continued climate change. Here, we take advantage of a transplant experiment mimicking up‐slope seed dispersal for a suite of ecologically diverse populations of yellow monkeyflower (Mimulus guttatus sensu lato) into a high‐elevation common garden during an extreme drought period in the Sierra Nevada mountains, California, USA. We use a demographic approach to quantify fitness and test for selection on life history traits in local versus lower‐elevation populations and in normal versus drought years to test the potential for up‐slope migration and phenotypic selection to alter the distribution of key life history traits in montane environments. We find that lower‐elevation populations tend to outperform local populations, confirming the potential for up‐slope migration. Although selection generally favored some local montane traits, including larger flowers and larger stem size at flowering, drought conditions tended to select for earlier flowering typical of lower‐elevation genotypes. Taken together, this suggests that monkeyflower lineages moving upward in elevation could experience selection for novel trait combinations, particularly under warmer and drier conditions that are predicted to occur with continued climate change.     

Leaf size variations in a dominant desert shrub,Reaumuria soongarica,adapted to heterogeneous environments

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species‐wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.     

Keeping your options open: Maintenance of thermal plasticity during adaptation to a stable environment

Phenotypic plasticity may allow species to cope with environmental variation. The study of thermal plasticity and its evolution helps understanding how populations respond to variation in temperature. In the context of climate change, it is essential to realize the impact of historical differences in the ability of populations to exhibit a plastic response to thermal variation and how it evolves during colonization of new environments. We have analyzed the real‐time evolution of thermal reaction norms of adult and juvenile traits in Drosophila subobscura populations from three locations of Europe in the laboratory. These populations were kept at a constant temperature of 18ºC, and were periodically assayed at three experimental temperatures (13ºC, 18ºC, and 23ºC). We found initial differentiation between populations in thermal plasticity as well as evolutionary convergence in the shape of reaction norms for some adult traits, but not for any of the juvenile traits. Contrary to theoretical expectations, an overall better performance of high latitude populations across temperatures in early generations was observed. Our study shows that the evolution of thermal plasticity is trait specific, and that a new stable environment did not limit the ability of populations to cope with environmental challenges.     

Earlier snow melt and reduced summer precipitation alter floral traits important to pollination

Climate change can cause changes in expression of organismal traits that influence fitness. In flowering plants, floral traits can respond to drought, and that phenotypic plasticity has the potential to affect pollination and plant reproductive success. Global climate change is leading to earlier snow melt in snow-dominated ecosystems as well as affecting precipitation during the growing season, but the effects of snow melt timing on floral morphology and rewards remain unknown. We conducted crossed manipulations of spring snow melt timing (early vs. control) and summer monsoon precipitation (addition, control, and reduction) that mimicked recent natural variation, and examined plastic responses in floral traits of Ipomopsis aggregata over 3 years in the Rocky Mountains. We tested whether increased summer precipitation compensated for earlier snow melt, and if plasticity was associated with changes in soil moisture and/or leaf gas exchange. Lower summer precipitation decreased corolla length, style length, corolla width, sepal width, and nectar production, and increased nectar concentration. Earlier snow melt (taking into account natural and experimental variation) had the same effects on those traits and decreased inflorescence height. The effect of reduced summer precipitation was stronger in earlier snow melt years for corolla length and sepal width. Trait reductions were explained by drier soil during the flowering period, but this effect was only partially explained by how drier soils affected plant water stress, as measured by leaf gas exchange. We predicted the effects of plastic trait changes on pollinator visitation rates, pollination success, and seed production using prior studies on I. aggregata. The largest predicted effect of drier soil on relative fitness components via plasticity was a decrease in male fitness caused by reduced pollinator rewards (nectar production). Early snow melt and reduced precipitation are strong drivers of phenotypic plasticity, and both should be considered when predicting effects of climate change on plant traits in snow-dominated ecosystems.     

Will phenotypic plasticity affecting flowering phenology keep pace with climate change?

Rising temperatures have begun to shift flowering time, but it is unclear whether phenotypic plasticity can accommodate projected temperature change for this century. Evaluating clines in phenological traits and the extent and variation in plasticity can provide key information on assessing risk of maladaptation and developing strategies to mitigate climate change. In this study, flower phenology was examined in 52 populations of big sagebrush (Artemisia tridentata) growing in three common gardens. Flowering date (anthesis) varied 91 days from late July to late November among gardens. Mixed‐effects modeling explained 79% of variation in flowering date, of which 46% could be assigned to plasticity and genetic variation in plasticity and 33% to genetics (conditional R² = 0.79, marginal R² = 0.33). Two environmental variables that explained the genetic variation were photoperiod and the onset of spring, the Julian date of accumulating degree‐days >5 °C reaching 100. The genetic variation was mapped for contemporary and future climates (decades 2060 and 2090), showing flower date change varies considerably across the landscape. Plasticity was estimated to accommodate, on average, a ±13‐day change in flowering date. However, the examination of genetic variation in plasticity suggests that the magnitude of plasticity could be affected by variation in the sensitivity to photoperiod and temperature. In a warmer common garden, lower‐latitude populations have greater plasticity (+16 days) compared to higher‐latitude populations (+10 days). Mapped climatypes of flowering date for contemporary and future climates illustrate the wide breadth of plasticity and large geographic overlap. Our research highlights the importance of integrating information on genetic variation, phenotypic plasticity and climatic niche modeling to evaluate plant responses and elucidate vulnerabilities to climate change.     

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east‐west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.