The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change

Species are the unit of analysis in many global change and conservation biology studies; however, species are not uniform entities but are composed of different, sometimes locally adapted, populations differing in plasticity. We examined how intraspecific variation in thermal niches and phenotypic plasticity will affect species distributions in a warming climate. We first developed a conceptual model linking plasticity and niche breadth, providing five alternative intraspecific scenarios that are consistent with existing literature. Secondly, we used ecological niche‐modeling techniques to quantify the impact of each intraspecific scenario on the distribution of a virtual species across a geographically realistic setting. Finally, we performed an analogous modeling exercise using real data on the climatic niches of different tree provenances. We show that when population differentiation is accounted for and dispersal is restricted, forecasts of species range shifts under climate change are even more pessimistic than those using the conventional assumption of homogeneously high plasticity across a species' range. Suitable population‐level data are not available for most species so identifying general patterns of population differentiation could fill this gap. However, the literature review revealed contrasting patterns among species, urging greater levels of integration among empirical, modeling and theoretical research on intraspecific phenotypic variation.     

Evolutionary potential of Chamaecrista fasciculata in relation to climate change. I. Clinal patterns of selection along an environmental gradient in the great plains

Climate change will alter natural selection on native plant populations. Little information is available to predict how selection will change in the future and how populations will respond. Insight can be obtained by comparing selection regimes in current environments to selection regimes in environments similar to those predicted for the future. To mimic predicted temporal change in climate, three natural populations of the annual legume Chamaecrista fasciculata were sampled from a climate gradient in the Great Plains and progeny of formal crosses were reciprocally planted back into common gardens across this climate gradient. In each garden, native populations produced significantly more seed than the other populations, providing strong evidence of local adaptation. Phenotypic selection analysis conducted by site showed that plants with slower reproductive development, more leaves, and thicker leaves were favored in the most southern garden. Evidence of clinal variation in selection regimes was also found; selection coefficients were ordered according to the latitude of the common gardens. The adaptive value of native traits was indicated by selection toward the mean of local populations. Repeated clinal patterns in linear and nonlinear selection coefficients among populations and within and between sites were found. To the extent that temporal change in climate into the future will parallel the differences in selection across this spatial gradient, this study suggests that selection regimes will be displaced northward and different trait values will be favored in natural populations.     

Plant population differentiation and climate change: responses of grassland species along an elevational gradient

Mountain ecosystems are particularly susceptible to climate change. Characterizing intraspecific variation of alpine plants along elevational gradients is crucial for estimating their vulnerability to predicted changes. Environmental conditions vary with elevation, which might influence plastic responses and affect selection pressures that lead to local adaptation. Thus, local adaptation and phenotypic plasticity among low and high elevation plant populations in response to climate, soil and other factors associated with elevational gradients might underlie different responses of these populations to climate warming. Using a transplant experiment along an elevational gradient, we investigated reproductive phenology, growth and reproduction of the nutrient‐poor grassland species Ranunculus bulbosus, Trifolium montanum and Briza media. Seeds were collected from low and high elevation source populations across the Swiss Alps and grown in nine common gardens at three different elevations with two different soil depths. Despite genetic differentiation in some traits, the results revealed no indication of local adaptation to the elevation of population origin. Reproductive phenology was advanced at lower elevation in low and high elevation populations of all three species. Growth and reproduction of T. montanum and B. media were hardly affected by garden elevation and soil depth. In R. bulbosus, however, growth decreased and reproductive investment increased at higher elevation. Furthermore, soil depth influenced growth and reproduction of low elevation R. bulbosus populations. We found no evidence for local adaptation to elevation of origin and hardly any differences in the responses of low and high elevation populations. However, the consistent advanced reproductive phenology observed in all three species shows that they have the potential to plastically respond to environmental variation. We conclude that populations might not be forced to migrate to higher elevations as a consequence of climate warming, as plasticity will buffer the detrimental effects of climate change in the three investigated nutrient‐poor grassland species.     

Population variation in the thermal response to climate change reveals differing sensitivity in a benthic shark

Many species with broad distributions are exposed to different thermal regimes which often select for varied phenotypes. This intraspecific variation is often overlooked but may be critical in dictating the vulnerability of different populations to environmental change. We reared Port Jackson shark (Heterodontus portusjacksoni) eggs from two thermally discrete populations (i.e. Jervis Bay and Adelaide) under each location's present‐day mean temperatures, predicted end‐of‐century temperatures and under reciprocal‐cross conditions to establish intraspecific thermal sensitivity. Rearing temperatures strongly influenced ?O_2Max and critical thermal limits, regardless of population, indicative of acclimation processes. However, there were significant population‐level effects, such that Jervis Bay sharks, regardless of rearing temperature, did not exhibit differences in ?O_2Rest, but under elevated temperatures exhibited reduced maximum swimming activity with step‐wise increases in temperature. In contrast, Adelaide sharks reared under elevated temperatures doubled their ?O_2Rest, relative to their present‐day temperature counterparts; however, maximum swimming activity was not influenced. With respect to reciprocal‐cross comparisons, few differences were detected between Jervis Bay and Adelaide sharks reared under ambient Jervis Bay temperatures. Similarly, juveniles (from both populations) reared under Adelaide conditions had similar thermal limits and swimming activity (maximum volitional velocity and distance) to each other, indicative of conserved acclimation capacity. However, under Adelaide temperatures, the ?O_2Rest of Jervis Bay sharks was greater than that of Adelaide sharks. This indicates that the energetics of cooler water population (Adelaide) is likely more thermally sensitive than that of the warmer population (Jervis Bay). While unique to elasmobranchs, these data provide further support that by treating species as static, homogeneous populations, we ignore the impacts of thermal history and intraspecific variation on thermal sensitivity. With climate change, intraspecific variation will manifest as populations move, demographics change or extirpations occur, starting with the most sensitive populations.     

Evidence of adaptive divergence in plasticity: density- and site-dependent selection on shade-avoidance responses in Impatiens capensis

We investigated the conditions under which plastic responses to density are adaptive in natural populations of Impatiens capensis and determined whether plasticity has evolved differently in different selective environments. Previous studies showed that a population that evolved in a sunny site exhibited greater plasticity in response to density than did a population that evolved in a woodland site. Using replicate inbred lines in a reciprocal transplant that included a density manipulation, we asked whether such population differentiation was consistent with the hypothesis of adaptive divergence. We hypothesized that plasticity would be more strongly favored in the sunny site than in the woodland site; consequently, we predicted that selection would be more strongly density dependent in the sunny site, favoring the phenotype that was expressed at each density. Selection on internode length and flowering date was consistent with the hypothesis of adaptive divergence in plasticity. Few costs or benefits of plasticity were detected independently from the expressed phenotype, so plasticity was selected primarily through selection on the phenotype. Correlations between phenotypes and their plasticity varied with the environment and would cause indirect selection on plasticity to be environment dependent. We showed that an appropriate plastic response even to a rare environment can greatly increase genotypic fitness when that environment is favorable. Selection on the measured characters contributed to local adaptation and fully accounted for fitness differences between populations in all treatments except the woodland site at natural density.     

Limitations of cross‐ and multigenerational plasticity for marine invertebrates faced with global climate change

Although cross generation (CGP) and multigenerational (MGP) plasticity have been identified as mechanisms of acclimation to global change, the weight of evidence indicates that parental conditioning over generations is not a panacea to rescue stress sensitivity in offspring. For many species, there were no benefits of parental conditioning. Even when improved performance was observed, this waned over time within a generation or across generations and fitness declined. CGP and MGP studies identified resilient species with stress tolerant genotypes in wild populations and selected family lines. Several bivalves possess favourable stress tolerance and phenotypically plastic traits potentially associated with genetic adaptation to life in habitats where they routinely experience temperature and/or acidification stress. These traits will be important to help ‘climate proof’ shellfish ventures. Species that are naturally stress tolerant and those that naturally experience a broad range of environmental conditions are good candidates to provide insights into the physiological and molecular mechanisms involved in CGP and MGP. It is challenging to conduct ecologically relevant global change experiments over the long times commensurate with the pace of changing climate. As a result, many studies present stressors in a shock‐type exposure at rates much faster than projected scenarios. With more gradual stressor introduction over longer experimental durations and in context with conditions species are currently acclimatized and/or adapted to, the outcomes for sensitive species might differ. We highlight the importance to understand primordial germ cell development and the timing of gametogenesis with respect to stressor exposure. Although multigenerational exposure to global change stressors currently appears limited as a universal tool to rescue species in the face of changing climate, natural proxies of future conditions (upwelling zones, CO₂ vents, naturally warm habitats) show that phenotypic adjustment and/or beneficial genetic selection is possible for some species, indicating complex plasticity–adaptation interactions.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

Testing for adaptation to climate in Arabidopsis thaliana: a calibrated common garden approach

BACKGROUND AND AIMS: A recent method used to test for local adaptation is a common garden experiment where analyses are calibrated to the environmental conditions of the garden. In this study the calibrated common garden approach is used to test for patterns of adaptation to climate in accessions of Arabidopsis thaliana. METHODS: Seedlings from 21 accessions of A. thaliana were planted outdoors in College Park, MD, USA, and development was monitored during the course of a growing season. ANOVA and multiple regression analysis were used to determine if development traits were significant predictors of plant success. Previously published data relating to accessional differences in genetic and physiological characters were also examined. Historical records of climate were used to evaluate whether properties of the site of origin of an accession affected the fitness of plants in a novel environment. KEY RESULTS: By calibrating the analysis to the climatic conditions of the common garden site, performance differences were detected among the accessions consistent with a pattern of adaptation to latitude and climatic conditions. Relatively higher accession fitness was predicted by a latitude and climatic history similar to that of College Park in April and May during the main growth period of this experiment. The climatic histories of the accessions were better predictors of performance than many of the life-history and growth measures taken during the experiment. CONCLUSIONS: It is concluded that the calibrated common garden experiment can detect local adaptation and guide subsequent reciprocal transplant experiments.     

The evolution of phenotypic plasticity in spatially structured environments: implications of intraspecific competition, plasticity costs and environmental characteristics

We model the evolution of reaction norms focusing on three aspects: frequency-dependent selection arising from resource competition, maintenance and production costs of phenotypic plasticity, and three characteristics of environmental heterogeneity (frequency of environments, their intrinsic carrying capacity and the sensitivity to phenotypic maladaptation in these environments). We show that (i) reaction norms evolve so as to trade adaptation for acquiring resources against cost avoidance; (ii) maintenance costs cause reaction norms to better adapt to frequent rather than to infrequent environments, whereas production costs do not; and (iii) evolved reaction norms confer better adaptation to environments with low rather than with high intrinsic carrying capacity. The two previous findings contradict earlier theoretical results and originate from two previously unexplored features that are included in our model. First, production costs of phenotypic plasticity are only incurred when a given phenotype is actually produced. Therefore, they are proportional to the frequency of environments, and these frequencies thus affect the selection pressure to avoid costs just as much as the selection pressure to improve adaptation. This prevents the frequency of environments from affecting the evolving reaction norm. Secondly, our model describes the evolution of plasticity for a phenotype determining an individual's capability to acquire resources, and thus its realized carrying capacity. When individuals are distributed randomly across environments, they cannot avoid experiencing environments with intrinsically low carrying capacity. As selection pressures arising from the need to improve adaptation are stronger under such extreme conditions than under mild ones, better adaptation to environments with low rather than with high intrinsic carrying capacity results.     

Evolution in stressful environments II: adaptive value and costs of plasticity in response to low light in Sinapis arvensis

Plants possess a remarkable capacity to alter their phenotype in response to the highly heterogeneous light conditions they commonly encounter in natural environments. In the present study with the weedy annual plant Sinapis arvensis, we (a) tested for the adaptive value of phenotypic plasticity in morphological and life history traits in response to low light and (b) explored possible fitness costs of plasticity. Replicates of 31 half-sib families were grown individually in the greenhouse under full light and under low light (40% of ambient) imposed by neutral shade cloth. Low light resulted in a large increase in hypocotyl length and specific leaf area (SLA), a reduction in juvenile biomass and a delayed onset of flowering. Phenotypic selection analysis within each light environment revealed that selection favoured large SLA under low light, but not under high light, suggesting that the observed increase in SLA was adaptive. In contrast, plasticity in the other traits measured was maladaptive (i.e. in the opposite direction to that favoured by selection in the low light environment). We detected significant additive genetic variance in plasticity in most phenotypic traits and in fitness (number of seeds). Using genotypic selection gradient analysis, we found that families with high plasticity in SLA had a lower fitness than families with low plasticity, when the effect of SLA on fitness was statistically kept constant. This indicates that plasticity in SLA incurred a direct fitness cost. However, a cost of plasticity was only expressed under low light, but not under high light. Thus, models on the evolution of phenotypic plasticity will need to incorporate plasticity costs that vary in magnitude depending on environmental conditions.     

Local adaptation and cogradient selection in the alpine plant, Poa hiemata, along a narrow altitudinal gradient

Alpine environments are particularly susceptible to environmental changes associated with global warming but there is potential for alpine plants to adapt to warming if local adaptation occurs and gene flow allows genotypes adapted to low altitudes to colonize higher altitude sites. Here we examine the adaptive potential of a common alpine grass, Poa hiemata, within the restricted alpine habitat of Australian mountains, across a narrow altitudinal gradient replicated in three areas. Grasses at high altitude sites had shorter leaf lengths and larger circumferences than those at lower sites. Transplant experiments with clonal material and plants grown from seed indicated that these differences were partly genetic, with environmental and genetic factors both contributing to the differences between altitudes. Differences in altitudinal forms were also evident in a common garden experiment. Plants showed a home-site advantage in terms of survival. A fitness analysis indicated that at high altitude sites, selection favored plants with short leaves and larger circumferences, whereas these traits were selected in the opposite direction at the low altitude sites. These findings indicate cogradient selection and potential for both plastic and genotypic shifts in response to climate change in P. hiemata.     

Transgenerational plasticity as an important mechanism affecting response of clonal species to changing climate

In spite of the increasing number of studies on the importance of transgenerational plasticity for species response to novel environments, its effects on species ability to respond to climate change are still largely unexplored. We study the importance of transgenerational plasticity for response of a clonal species Festuca rubra. Individuals from four natural populations representing two levels of temperature and two levels of precipitation were cultivated in four growth chambers that simulate the temperature and precipitation of origin of the populations (maternal phase). Each population was represented in each growth chamber. After 6 months, single young ramets of these plants were reshuffled among the growth chambers and let to grow for additional 2 months (offspring phase). The results show that transgenerational effects (i.e., maternal phase conditions) significantly modify species response to novel climates, and the direction and intensity of the response depend on the climate of origin of the plants. For traits related to recourse acquisition, the conditions of maternal phase, either alone or in interaction mainly with climate of origin, had stronger effect than the conditions of cultivation. Overall, the maternal climate interacted more intensively with the climate of origin than with the offspring climate. The direction of the effect of the maternal climate was of different directions and intensities depending on plant origin and trait studied. The data demonstrated strong significant effects of conditions during maternal phase on species response to novel climates. These transgenerational affects were, however, not adaptive. Still, transgenerational plasticity may be an important driver of species response to novel conditions across clonal generations. These effects thus need to be carefully considered in future studies exploring species response to novel climates. This will also have strong effects on species performance under increasingly variable climates expected to occur with the climate change.     

Leaf‐trait plasticity and species vulnerability to climate change in a Mongolian steppe

Climate change is expected to modify plant assemblages in ways that will have major consequences for ecosystem functions. How climate change will affect community composition will depend on how individual species respond, which is likely related to interspecific differences in functional traits. The extraordinary plasticity of some plant traits is typically neglected in assessing how climate change will affect different species. In the Mongolian steppe, we examined whether leaf functional traits under ambient conditions and whether plasticity in these traits under altered climate could explain climate‐induced biomass responses in 12 co‐occurring plant species. We experimentally created three probable climate change scenarios and used a model selection procedure to determine the set of baseline traits or plasticity values that best explained biomass response. Under all climate change scenarios, plasticity for at least one leaf trait correlated with change in species performance, while functional leaf‐trait values in ambient conditions did not. We demonstrate that trait plasticity could play a critical role in vulnerability of species to a rapidly changing environment. Plasticity should be considered when examining how climate change will affect plant performance, species' niche spaces, and ecological processes that depend on plant community composition.     

Contrasting plant physiological adaptation to climate in the native and introduced range of Hypericum perforatum

How introduced plants, which may be locally adapted to specific climatic conditions in their native range, cope with the new abiotic conditions that they encounter as exotics is not well understood. In particular, it is unclear what role plasticity versus adaptive evolution plays in enabling exotics to persist under new environmental circumstances in the introduced range. We determined the extent to which native and introduced populations of St. John's Wort (Hypericum perforatum) are genetically differentiated with respect to leaf-level morphological and physiological traits that allow plants to tolerate different climatic conditions. In common gardens in Washington and Spain, and in a greenhouse, we examined clinal variation in percent leaf nitrogen and carbon, leaf delta(13)C values (as an integrative measure of water use efficiency), specific leaf area (SLA), root and shoot biomass, root/shoot ratio, total leaf area, and leaf area ratio (LAR). As well, we determined whether native European H. perforatum experienced directional selection on leaf-level traits in the introduced range and we compared, across gardens, levels of plasticity in these traits. In field gardens in both Washington and Spain, native populations formed latitudinal clines in percent leaf N. In the greenhouse, native populations formed latitudinal clines in root and shoot biomass and total leaf area, and in the Washington garden only, native populations also exhibited latitudinal clines in percent leaf C and leaf delta(13)C. Traits that failed to show consistent latitudinal clines instead exhibited significant phenotypic plasticity. Introduced St. John's Wort populations also formed significant or marginally significant latitudinal clines in percent leaf N in Washington and Spain, percent leaf C in Washington, and in root biomass and total leaf area in the greenhouse. In the Washington common garden, there was strong directional selection among European populations for higher percent leaf N and leaf delta(13)C, but no selection on any other measured trait. The presence of convergent, genetically based latitudinal clines between native and introduced H. perforatum, together with previously published molecular data, suggest that native and exotic genotypes have independently adapted to a broad-scale variation in climate that varies with latitude.     

Phenotypic distribution models corroborate species distribution models: A shift in the role and prevalence of a dominant prairie grass in response to climate change

Phenotypic distribution within species can vary widely across environmental gradients but forecasts of species’ responses to environmental change often assume species respond homogenously across their ranges. We compared predictions from species and phenotype distribution models under future climate scenarios for Andropogon gerardii, a widely distributed, dominant grass found throughout the central United States. Phenotype data on aboveground biomass, height, leaf width, and chlorophyll content were obtained from 33 populations spanning a ~1000 km gradient that encompassed the majority of the species’ environmental range. Species and phenotype distribution models were trained using current climate conditions and projected to future climate scenarios. We used permutation procedures to infer the most important variable for each model. The species‐level response to climate was most sensitive to maximum temperature of the hottest month, but phenotypic variables were most sensitive to mean annual precipitation. The phenotype distribution models predict that A. gerardii could be largely functionally eliminated from where this species currently dominates, with biomass and height declining by up to ~60% and leaf width by ~20%. By the 2070s, the core area of highest suitability for A. gerardii is projected to shift up to ~700 km northeastward. Further, short‐statured phenotypes found in the present‐day short grass prairies on the western periphery of the species’ range will become favored in the current core ~800 km eastward of their current location. Combined, species and phenotype models predict this currently dominant prairie grass will decline in prevalence and stature. Thus, sourcing plant material for grassland restoration and forage should consider changes in the phenotype that will be favored under future climate conditions. Phenotype distribution models account for the role of intraspecific variation in determining responses to anticipated climate change and thereby complement predictions from species distributions models in guiding climate adaptation strategies.     

Factors influencing plasticity in the arrival‐breeding interval in a migratory species reacting to climate change

Climate change is profoundly affecting the phenology of many species. In migratory birds, there is evidence for advances in their arrival time at the breeding ground and their timing of breeding, yet empirical studies examining the interdependence between arrival and breeding time are lacking. Hence, evidence is scarce regarding how breeding time may be adjusted via the arrival‐breeding interval to help local populations adapt to local conditions or climate change. We used long‐term data from an intensively monitored population of the northern wheatear (Oenanthe oenanthe) to examine the factors related to the length of 734 separate arrival‐to‐breeding events from 549 individual females. From 1993 to 2017, the mean arrival and egg‐laying dates advanced by approximately the same amount (~5–6 days), with considerable between‐individual variation in the arrival‐breeding interval. The arrival‐breeding interval was shorter for: (a) individuals that arrived later in the season compared to early‐arriving birds, (b) for experienced females compared to first‐year breeders, (c) as spring progressed, and (d) in later years compared to earlier ones. The influence of these factors was much larger for birds arriving earlier in the season compared to later arriving birds, with most effects on variation in the arrival‐breeding interval being absent in late‐arriving birds. Thus, in this population it appears that the timing of breeding is not constrained by arrival for early‐ to midarriving birds, but instead is dependent on local conditions after arrival. For late‐arriving birds, however, the timing of breeding appears to be influenced by arrival constraints. Hence, impacts of climate change on arrival dates and local conditions are expected to vary for different parts of the population, with potential negative impacts associated with these factors likely to differ for early‐ versus late‐arriving birds.     

Forest tree species adaptation to climate across biomes: Building on the legacy of ecological genetics to anticipate responses to climate change

Intraspecific variation plays a critical role in extant and future forest responses to climate change. Forest tree species with wide climatic niches rely on the intraspecific variation resulting from genetic adaptation and phenotypic plasticity to accommodate spatial and temporal climate variability. A centuries-old legacy of forest ecological genetics and provenance trials has provided a strong foundation upon which to continue building on this knowledge, which is critical to maintain climate-adapted forests. Our overall objective is to understand forest trees intraspecific responses to climate across species and biomes, while our specific objectives are to describe ecological genetics models used to build our foundational knowledge, summarize modeling approaches that have expanded the traditional toolset, and extensively review the literature from 1994 to 2021 to highlight the main contributions of this legacy and the new analyzes of provenance trials. We reviewed 103 studies comprising at least three common gardens, which covered 58 forest tree species, 28 of them with range-wide studies. Although studies using provenance trial data cover mostly commercially important forest tree species from temperate and boreal biomes, this synthesis provides a global overview of forest tree species adaptation to climate. We found that evidence for genetic adaptation to local climate is commonly present in the species studied (79%), being more common in conifers (87.5%) than in broadleaf species (67%). In 57% of the species, clines in fitness-related traits were associated with temperature variables, in 14% of the species with precipitation, and in 25% of the species with both. Evidence of adaptation lags was found in 50% of the species with range-wide studies. We conclude that ecological genetics models and analysis of provenance trial data provide excellent insights on intraspecific genetic variation, whereas the role and limits of phenotypic plasticity, which will likely determine the fate of extant forests, is vastly understudied.     

Plant responses to simulated warming and drought: a comparative study of functional plasticity between congeneric mid and high elevation species

Aims Effects of climate change, especially changes in temperatures and precipitation patterns, are particularly pronounced in alpine regions. In response, plants may exhibit phenotypic plasticity in key functional traits allowing short-term adjustment to novel conditions. However, little is known about the degree of phenotypic plasticity of high elevation species relative to mid elevation congeners.Methods We transplanted 14 herbaceous perennial species from high elevation into two common gardens (1050 and 2000 m.a.s.l.) in the Swiss Alps, and we examined plastic responses in key functional traits to changes in temperature and soil water availability. This design was replicated with 14 congeneric species from mid elevation to assess if the degree of phenotypic plasticity differs between mid and high elevation species. Survival was assessed across two growing seasons, while aboveground biomass and specific leaf area (SLA) were measured after the first growing season, and biomass allocation to belowground and reproductive structures after the second. Moreover, a phenotypic plasticity index was calculated for the functional traits to compare the degree of plasticity between mid and high elevation species.Important findings Aboveground biomass was higher in mid elevation species relative to high elevation congeners in all treatments, yet decreased for both with elevation and drought. Similarly, SLA decreased with elevation and drought. Root mass fraction (RMF) was generally higher in high elevation species, and decreased with drought at the lower site. Drought increased the allocation to reproductive structures, especially when plants were grown at their elevation of origin. Interestingly, no difference was found in the degree of phenotypic plasticity averaged across mid and high elevation species for any of the studied functional traits. These results indicate that phenotypic plasticity in the focal traits did not depend on the elevation of origin of the species. Plasticity was not related to environmental heterogeneity, nor constrained by selective pressures at high elevation. However, both species groups showed a remarkable capacity for short-term acclimation to a prospective climate through rapid adjustments in key functional traits.