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.     

DISENTANGLING THE ROLE OF PHENOTYPIC PLASTICITY AND GENETIC DIVERGENCE IN CONTEMPORARY ECOTYPE FORMATION DURING A BIOLOGICAL INVASION

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full‐sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.     

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 in reproduction and growth among 52 range‐wide populations of a Mediterranean conifer: adaptive responses to environmental stress

A plastic response towards enhanced reproduction is expected in stressful environments, but it is assumed to trade off against vegetative growth and efficiency in the use of available resources deployed in reproduction [reproductive efficiency (RE)]. Evidence supporting this expectation is scarce for plants, particularly for long‐lived species. Forest trees such as Mediterranean pines provide ideal models to study the adaptive value of allocation to reproduction vs. vegetative growth given their among‐population differentiation for adaptive traits and their remarkable capacity to cope with dry and low‐fertility environments. We studied 52 range‐wide Pinus halepensis populations planted into two environmentally contrasting sites during their initial reproductive stage. We investigated the effect of site, population and their interaction on vegetative growth, threshold size for female reproduction, reproductive–vegetative size relationships and RE. We quantified correlations among traits and environmental variables to identify allocation trade‐offs and ecotypic trends. Genetic variation for plasticity was high for vegetative growth, whereas it was nonsignificant for reproduction. Size‐corrected reproduction was enhanced in the more stressful site supporting the expectation for adverse conditions to elicit plastic responses in reproductive allometry. However, RE was unrelated with early reproductive investment. Our results followed theoretical predictions and support that phenotypic plasticity for reproduction is adaptive under stressful environments. Considering expectations of increased drought in the Mediterranean, we hypothesize that phenotypic plasticity together with natural selection on reproductive traits will play a relevant role in the future adaptation of forest tree species.     

Host environment and local genetic adaptation determine phenotype in parasitic Rhinanthus angustifolius

Plants have a remarkable capacity to adapt to local environmental conditions, which can result in ecotypic differentiation. Patterns of differentiation can, however, also be influenced by the extensive phenotypic plasticity exhibited by many plant species. In this study, we evaluated the distinctness of two putative ecotypes of the parasitic herb Rhinanthus angustifolius. We compared population means of characters commonly used to distinguish between the putative ecotypes after growing individuals of R. angustifolius with a variety of host species in a common garden. Resulting data were also pooled over environments to study how phenotypic plasticity affects the distinctness of ecotypes and individual populations. Except for node number, most of the characters were plastic. The pattern of differentiation was consistent with the existence of two, or possibly three, habitat‐related groups of populations; however, we observed considerable overlap in flowering time and morphological characters after pooling data across host environments. The results show that the complex phenological and morphological variation in R. angustifolius is caused by a combination of genetically determined ecotypic differentiation and plastic responses to the host environment and other factors. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2016, 180 , 89–103.     

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.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Adaptive plasticity and epigenetic variation in response to warming in an Alpine plant

Environmentally induced phenotypic plasticity may be a critical component of response to changing environments. We examined local differentiation and adaptive phenotypic plasticity in response to elevated temperature in half‐sib lines collected across an elevation gradient for the alpine herb, Wahlenbergia ceracea. Using Amplified Fragment Length Polymorphism (AFLP), we found low but significant genetic differentiation between low‐ and high‐elevation seedlings, and seedlings originating from low elevations grew faster and showed stronger temperature responses (more plasticity) than those from medium and high elevations. Furthermore, plasticity was more often adaptive for plants of low‐elevation origin and maladaptive for plants of high elevation. With methylation sensitive‐AFLP (MS‐AFLP), we revealed an increase in epigenetic variation in response to temperature in low‐elevation seedlings. Although we did not find significant direct correlations between MS‐AFLP loci and phenotypes, our results demonstrate that adaptive plasticity in temperature response to warming varies over fine spatial scales and suggest the involvement of epigenetic mechanisms in this response.     

Invasive plants do not display greater phenotypic plasticity than their native or non‐invasive counterparts: a meta‐analysis

Phenotypic plasticity is commonly considered as a trait associated with invasiveness in alien plants because it may enhance the ability of plants to occupy a wide range of environments. Although the evidence of greater phenotypic plasticity in invasive plants is considerable, it is not yet conclusive. We used a meta‐analysis approach to evaluate whether invasive plant species show greater phenotypic plasticity than their native or non‐invasive counterparts. The outcome of such interspecific comparisons may be biased when phylogenetic relatedness is not taken into account. Consequently, species pairs belonged to the same genus, tribe or family. The meta‐analysis included 93 records from 35 studies reporting plastic responses to light, nutrients, water, CO₂, herbivory and support availability. Contrary to what is often assumed, overall, phenotypic plasticity was similar between invasive plants and native or non‐invasive closely related species. The same result was found when separate analyses were conducted for trait plasticity to nutrients, light and water availability. Thus, invasive plant species and their native or non‐invasive counterparts are equally capable of displaying functional responses to environmental heterogeneity. The colonization of a wide range of environments by invasive plants could be due to their capacity to undergo adaptive ecotypic differentiation rather than to their ability to display plastic responses. Alternatively, phenotypic plasticity might play a role in plant invasion, but only during the initial phases, when tolerance of the novel environment is essential for plant survival. Afterwards, once alien plants are identified as invaders, the magnitude of phenotypic plasticity might be reduced after selection of the optimum phenotypes in each habitat. The identification of plant traits that consistently predict invasiveness might be a futile task because different traits favor invasiveness in different environments. Approaches at the local scale, focusing on the ecology of specific invasive plants, could be more fruitful than global macro‐analyses.     

Phenotypic plasticity of four Chenopodiaceae species with contrasting saline–sodic tolerance in response to increased salinity–sodicity

It is unknown whether phenotypic plasticity in fitness‐related traits is associated with salinity–sodicity tolerance. This study compared growth and allocation phenotypic plasticity in two species with low salinity–sodicity tolerance (Chenopodium acuminatum and C. stenophyllum) and two species with high salinity–sodicity tolerance (Suaeda glauca and S. salsa) in a pot experiment in the Songnen grassland, China. While the species with low tolerance had higher growth and allocation plasticity than the highly tolerant species, the highly tolerant species only adjusted their growth traits and maintained higher fitness (e.g., plant height and total biomass) in response to increased soil salinity–sodicity, with low biomass allocation plasticity. Most plasticity is “apparent” plasticity (ontogenetic change), and only a few traits, for example, plant height:stem diameter ratio and root:shoot biomass ratio, represent “real” plasticity (real change in response to the environment). Our results show that phenotypic plasticity was negatively correlated with saline–sodic tolerance and could be used as an index of species sensitivity to soil salinity–sodicity.     

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta‐analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.     

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 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.

     

Natural selection on immune defense: A field experiment

Predicting the evolution of phenotypic traits requires an understanding of natural selection on them. Despite its indispensability in the fight against parasites, selection on host immune defense has remained understudied. Theory predicts immune traits to be under stabilizing selection due to associated trade‐offs with other fitness‐related traits. Empirical studies, however, report mainly positive directional selection. This discrepancy could be caused by low phenotypic variation in the examined individuals and/or variation in host resource level that confounds trade‐offs in empirical studies. In a field experiment where we maintained Lymnaea stagnalis snails individually in cages in a lake, we investigated phenotypic selection on two immune defense traits, phenoloxidase (PO)‐like activity and antibacterial activity, in hemolymph. We used a diverse laboratory population and manipulated snail resource level by limiting their food supply. For six weeks, we followed immune activity, growth, and two fitness components, survival and fecundity of snails. We found that PO‐like activity and growth were under stabilizing selection, while antibacterial activity was under positive directional selection. Selection on immune traits was mainly driven by variation in survival. The form of selection on immune defense apparently depends on the particular trait, possibly due to its importance for countering the present parasite community.     

Intraspecific variation in seedling growth responses of a relict tree species Euptelea pleiospermum to precipitation manipulation along an elevation gradient

Improving the accuracy of predictions regarding how plants respond to climate change is crucial to protecting biodiversity. However, little is known about the effects of seed source and elevation on the response of mountain plant species to reductions in precipitation. Here, we collected seeds of a tree species (Euptelea pleiospermum) from three seed sources and carried out a two-growing-season reciprocal transplant experiment with precipitation manipulation at three sites along an elevation gradient in the Shennongjia Mountains, central China. Variations in whole-plant traits, leaf traits, and root traits were investigated. We found that most plant traits of E. pleiospermum seedlings were affected by reductions in precipitation, and responses varied among different elevations and seed sources. Whole-plant traits, root biomass, and leaf traits related to photosynthesis capacity decreased under reduced precipitation treatments at mid and high elevation sites. Thus, climate change induced drought will likely have a negative influence on seedling growth at mid and high elevation regions. In addition, a home-site advantage in whole-plant traits and root traits was observed. However, the responses of leaf traits in most cases were not affected by seed source because of higher phenotypic plasticity. Our results suggested that both local adaptation and phenotypic plasticity were important in seedling growth responses to reduced precipitation. We also highlight the importance of taking intraspecific variation into account when studying the response of plants to changes in climate.

     

Geographic range size, seedling ecophysiology and phenotypic plasticity in Australian Acacia species

Aim The degree to which eco‐physiological traits critical to seedling establishment are related to differences in geographic range size among species is not well understood. Here, we first tested the idea that seedling eco‐physiological attributes associated with establishment differ between narrowly distributed and geographically widespread plant species. Secondly, we tested the notion that species occupying wide geographic ranges have greater phenotypic plasticity in response to the environment than contrasted species with more restricted distributions. Location Eastern Australia. Methods We compared five pairs of geographically restricted and widespread Acacia species grown under glasshouse conditions for differences in seedling relative growth rate and associated allocational, morphological and physiological traits. We then examined whether widespread species displayed greater phenotypic plasticity in these traits than narrowly distributed species in response to changes in soil water availability. Results Neither relative growth rate nor any measure of biomass accumulation or allocation differed significantly between seedlings of narrowly distributed and widespread species. In addition, the plasticity of biomass allocation was not greater in widespread species. However, the leaflets of widespread species had higher photosynthetic capacity and greater plasticity of water use efficiency than the leaflets of narrowly distributed species. Main conclusions We demonstrated fundamental differences in the physiology and plasticity of leaflets of widespread and narrowly distributed species. The greater plasticity of these seedling leaflet traits may allow widespread Acacia species to utilize a wider range of environmental conditions in relation to soil moisture than restricted Acacia species. However, we did not find corresponding differences in mean or plasticity of seedling growth and allocational traits. In general, we suggest that relationships between rarity and species traits are both context and taxon specific.     

Ecophysiological plasticity and local differentiation help explain the invasion success of Taraxacum officinale (dandelion) in South America

Plasticity and local adaptation have been suggested as two main mechanisms that alien species use to successfully tolerate and invade broad geographic areas. In the present study, we try answer the question if the mechanism for the broad distributional range of T. officinale is for phenotypic plasticity, ecotypic adaptation or both. For this, we used individuals of T. officinale originated from seeds collected in five localities along its latitudinal distribution range in the southern‐hemisphere. Seedlings were acclimated at 5 and 25°C for one month. After the acclimation period we evaluated ecophysiological and cytogenetic traits. Additionally, we assessed the fitness at each temperature by recording the seed output of individuals from different localities. Finally, we performed a manipulative experiment in order to assess the tolerance to herbivory and competitive ability between T. officinale from all origins and Hypochaeris scorzonerae a co‐occurring native species. Overall, individuals of T. officinale showed high plasticity and ecotypic adaptation for all traits assessed in this study. Changes both in physiology and morphology observed in T. officinale from different origins were mostly correlated, enhancing their ecophysiological performance in temperatures similar to those of their origin. Additionally, all localities showed the same chromosome number and ploidy level. On the other hand, all individuals showed an increase the seed output at 25°C, but those from northern localities increased more. T. officinale from all origins was not significantly affected by herbivory while native showed a negative effect. On the other hand, T. officinale exerted a strong negative effect on the native species, but this former not effected significantly to the invasive T. officinale. High plasticity and local adaptation in all ecophysiological traits, seed‐set and the low cytogenetic variability in T. officinale suggests that both strategies are present in this invasive plant species and are not mutually exclusive. Finally, higher tolerance to herbivory and competitive ability suggests that T. officinale could perform successfully in environments with different climatic conditions, and thus colonize and invade South‐America.     

Patterns and drivers of intraspecific variation in avian life history along elevational gradients: a meta‐analysis

Elevational gradients provide powerful natural systems for testing hypotheses regarding the role of environmental variation in the evolution of life‐history strategies. Case studies have revealed shifts towards slower life histories in organisms living at high elevations yet no synthetic analyses exist of elevational variation in life‐history traits for major vertebrate clades. We examined (i) how life‐history traits change with elevation in paired populations of bird species worldwide, and (ii) which biotic and abiotic factors drive elevational shifts in life history. Using three analytical methods, we found that fecundity declined at higher elevations due to smaller clutches and fewer reproductive attempts per year. By contrast, elevational differences in traits associated with parental investment or survival varied among studies. High‐elevation populations had shorter and later breeding seasons, but longer developmental periods implying that temporal constraints contribute to reduced fecundity. Analyses of clutch size data, the trait for which we had the largest number of population comparisons, indicated no evidence that phylogenetic history constrained species‐level plasticity in trait variation associated with elevational gradients. The magnitude of elevational shifts in life‐history traits were largely unrelated to geographic (altitude, latitude), intrinsic (body mass, migratory status), or habitat covariates. Meta‐population structure, methodological issues associated with estimating survival, or processes shaping range boundaries could potentially explain the nature of elevational shifts in life‐history traits evident in this data set. We identify a new risk factor for montane populations in changing climates: low fecundity will result in lower reproductive potential to recover from perturbations, especially as fewer than half of the species experienced higher survival at higher elevations.     

Global patterns of intraspecific leaf trait responses to elevation

Elevational gradients are often used to quantify how traits of plant species respond to abiotic and biotic environmental variations. Yet, such analyses are frequently restricted spatially and applied along single slopes or mountain ranges. Since we know little on the response of intraspecific leaf traits to elevation across the globe, we here perform a global meta‐analysis of leaf traits in 109 plant species located in 4 continents and reported in 71 studies published between 1983 and 2018. We quantified the intraspecific change in seven morpho‐ecophysiological leaf traits along global elevational gradients: specific leaf area (SLA), leaf mass per area (LMA), leaf area (LA), nitrogen concentration per unit of area (Narea), nitrogen concentration per unit mass (Nmass), phosphorous concentration per unit mass (Pmass) and carbon isotope composition (δ¹³C). We found LMA, Narea, Nmass and δ¹³C to significantly increase and SLA to decrease with increasing elevation. Conversely, LA and Pmass showed no significant pattern with elevation worldwide. We found significantly larger increase in Narea, Nmass, Pmass and δ¹³C with elevation in warmer regions. Larger responses to increasing elevation were apparent for SLA of herbaceous compared to woody species, but not for the other traits. Finally, we also detected evidences of covariation across morphological and physiological traits within the same elevational gradient. In sum, we demonstrate that there are common cross‐species patterns of intraspecific leaf trait variation across elevational gradients worldwide. Irrespective of whether such variation is genetically determined via local adaptation or attributed to phenotypic plasticity, the leaf trait patterns quantified here suggest that plant species are adapted to live on a range of temperature conditions. Since the distribution of mountain biota is predominantly shifting upslope in response to changes in environmental conditions, our results are important to further our understanding of how plants species of mountain ecosystems adapt to global environmental change.