Skip to main content

Patterns and driving factors of functional traits of desert species with different elevational distributions in the Tibetan Plateau and adjacent areas

Abstract

Variations in functional traits serve as measures of plants’ ability to adapt to environment. Exploring the patterns of functional traits of desert plants along elevational gradients is helpful to understand the responses and adaptation strategies of species to changing environments. However, it is unknown whether the relationship between functional traits and elevation is affected by differences in the species’ elevational distributions (elevation preference and species’ range). Importantly, most researches have concerned with differences in mean trait values and ignored intraspecific trait variation. Here, we measured functional traits of desert plants along a wide elevational gradient in the Tibetan Plateau and adjacent areas and explored functional trait patterns over elevation in species with different elevational distributions. We decomposed trait variation and further investigated characterizations of intraspecific variation. Ultimately, the main drivers of trait variation were identified using redundancy analysis. We found that species’ elevational distributions significantly influenced the relationship of functional traits such as plant height, leaf dry matter content, leaf thickness, leaf nitrogen and carbon content with elevation. Species with a lower elevational preference showed greater trait variation than species with a higher elevational preference, suggesting that species that prefer high elevation are more conservative facing environmental changes. We provide evidence that interspecific trait variation in leaf thickness and leaf carbon content decreased with increasing species’ range, indicating that increased variations in resistance traits within species make greater responsiveness to environmental changes, enabling species a wider range. Elevation, temperature and precipitation were the main drivers of trait variation in species with a low elevational preference, while the effect of precipitation on trait variation in species with a high elevational preference was not significant. This study sheds new insights on how plants with different elevational distributions regulate their ecological strategies to cope with changing environments.

Peer Review reports

Introduction

Predicting how different specie will respond to environmental changes is challenging due to the diversity of natural ecosystems [1]. Functional traits provide a method for disentangling community responses to environmental changes by linking environment with individual performance [2,3,4]. Functional traits are measurable characteristics of an individual that represent species adaptive responses to abiotic and biotic factors [5, 6]. There has been a classic problem in ecology on what cause functional trait variations along climate gradients [4, 7]. This knowledge is crucial for foreseeing how climate changes will affect species, species interactions, and ecosystem functioning [8]. The distribution of functional traits along environmental gradients is explained by mechanisms relating to physiological constraints on species [9]. Physiology-based theories presuppose that changes in the physical environment and physiological restraints control the distribution and evolution of organisms, with implications for the distribution of morphological traits [8].

Elevational gradients offer suitable environments for enhancing inference of the mechanical causes of ecosystem functioning due to their varied environmental and climatic circumstances [6]. Some functional traits, particularly those that relate to plant height, leaf size and resource acquisition, are strong predictors of plant performance, differ between species, and can be used to infer changes in ecosystem functioning at broad ecological scales [10, 11]. For example, with increasing elevation, leaves get thicker and smaller [4]. In harsher situations, plant height and leaf size tend to decline, whereas leaf nutrient contents vary with leaf morphology, elevation, and climatic conditions [12]. The variations in functional trait over elevational gradients are expected to explain plant ecological strategies [13].

Trait-based ecology has, up to this point, mostly emphasized the differences in traits between species [2, 14, 15]. However, there is mounting evidence that intraspecific variation, rather than interspecific variations, contributes more to trait variation caused by environmental factors [1, 16, 17]. Intraspecific variation accounts for about a quarter of total trait variation globally [17], but this proportion is predicted to increase in harsh environments due to the filtering effect of environment on the trait expression [18]. Large intraspecific trait variation may conceal or change the relationships among interspecific traits, limiting the usefulness of interspecific variation for ecological prediction at different scales [19]. Therefore, investigating intraspecific trait variations may provide a more comprehensive answer to community construction and ecosystem function maintenance [17].

The Tibetan Plateau region is characterized by high average elevations and wide elevational gradients, due to the topographic features and the atmospheric circulation characteristics, harboring not only unique alpine ecosystems, but also a variety of natural ecosystem types such as forests, meadows, steppes and deserts appear from southeast to northwest [20, 21]. Therefore, the Tibetan Plateau has nurtured many unique plant resources with high biological and genetic research values, which are important for biodiversity maintenance and biological resource conservation [21, 22]. The Tibetan Plateau is not only an important ecological security barrier, but also a sensitive and fragile zone to global climate changes. As the global climate change process advances, glacial retreat, permafrost melting and desertification are becoming more prominent, accelerating the degradation of vegetation in the Tibetan Plateau [23]. Exploring the status of desert plants on the Tibetan Plateau can provide a theoretical basis for desertification control, biological resource conservation and sustainable development of the ecosystem.

Plants respond differently to different elevational gradients, but there are fewer studies on the response of traits to different elevational distribution ranges. Interspecific and intraspecific trait variations are major components of plant functional trait variation, but intraspecific trait variations across a large elevational gradient merits further research [2, 24]. The objectives of the present study were: (i) to establish the relationship between functional traits and elevation in desert species with different elevational preferences and elevational distribution ranges; (ii) to explore the sources of variation of functional traits and to determine the proportion of intraspecific variation and (iii) to identify the main environmental factors that influence functional traits of species with different elevational distributions.

Materials and methods

Study area

According to comprehensive natural geographic zoning, distribution of deserts and desertification in China and land use conditions, we selected a typical desert ecosystem of the Tibetan Plateau and adjacent areas as our study area. The region reaches an average elevation of 4000 m a.s.l., and nearly a quarter of its northwestern area is alpine, with altitudes above 5000 m a.s.l. The moisture status of the study area has large differences, with annual precipitation mostly below 900 mm, decreasing from east to west and from south to north. The sampling sites were selected in the desert ecosystems of the Tibetan Plateau and adjacent areas, ranging from 813 to 5930 m a.s.l.

Sampling and trait measurements

A total of 414 study sites were selected in the desert ecosystems of the Tibetan Plateau and adjacent areas, and vegetation surveys were conducted over 4 years (Table S1). At each study site, typical and representative plant communities were selected and a 100 m × 100 m sampling area was established. Within the sampling area, five 10 m×10 m shrub sampling plots and nine 1 m×1 m grass sampling plots were set up to investigate the species composition in shrub and grass sampling plots, respectively.

For dominant species of the community, 6 functional traits were measured based on the relevance to plant survival strategies and the feasibility of field measurements [10, 25], including plant height, specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (LT), leaf carbon content (LCC) and leaf nitrogen content (LNC). Plant height of each species was measured at the same time as the vegetation survey distancing from soil to highest leaf. We selected 5–10 individuals and at least 10 leaves from each dominant species within a site for determination of functional traits. The SLA, LDMC and LT were measured referring to standard protocols [25]. The leaves were dried and crushed for the determination of LCC and LNC by elemental analyzer (Costech, Milano, Italy).

Elevational distribution indicators

The elevational preference (EP) and species’ range (SR) can be used to explain two aspects of elevational distributions, reflecting the species’ preference for habitat elevation and the range over which the species can be distributed, respectively [1]. A species’ EP represents its median elevation in relation to all species, calculating by the following formula. EP ranges from 0 to 1, with values close to 0 for species with median elevation that is near to the lower elevation and 1 for species with median elevation that is near to the higher elevation.

$${Elevational\, preference}_{i}=1+\left[\frac{{Ele\left(Med\right)}_{i}-Ele\left(Max\right)}{Ele\left(Max\right)-Ele\left(Min\right)}\right]$$

Where ‘Elevational preferencei’ is the elevational preference of species i, ‘Ele (Med)i’ is the median elevation of species i, ‘Ele (Max)’ is the maximum elevation of all species and ‘Ele (Min)’ is the minimum elevation of all species.

We calculated each target species’ SR, which reflects its elevational distribution in comparison to all species. SR ranges between 0 and 1 with values near to 0 for species with narrower elevation ranges, and 1 for species with wider elevational ranges. We estimated the SR by the following equation.

$${Species\, range}_{i}=\left[\frac{{Ele\left(Max\right)}_{i}-{Ele\left(Min\right)}_{i}}{Ele\left(Max\right)-Ele\left(Min\right)}\right]$$

Where ‘Species rangei’ is the distribution range of species i, ‘Ele (Max)i’ is the maximum elevation of species i and ‘Ele (Min)i’ is the minimum elevation of species i.

Statistical analysis

To examine the patterns of plant trait variation in desert species with different elevational distributions over elevation, we constructed mixed-effects models for plant height, SLA, LDMC, LT, LNC and LCC using the lmer function from the lme4 package.

By fitting linear mixed effects models with a fixed intercept and random effects for region, site, functional group and species, we quantified the amount of trait variation for each species and trait at each nested scale using variance decomposition. The random effect variances in this equation stand in for variance between regions, sites, functional groups and species, whereas the residual variance represents samples within species (intraspecific trait variation).

We ran a redundancy analysis (RDA) on all trait measurements and elevation, temperature, precipitation, soil pH, soil electrical conductivity (EC), soil clay content (Clay), soil sand content (Sand), soil nitrogen content (SNC) and soil carbon content (SCC) to determine the relationships between functional traits and environmental factors. RDA was analyzed by the rda function form the vegan package. All the data analysis was carried out using R (R Development Core Team 2022).

Results

Patterns of desert plant traits along elevational gradients

At the overall level, plant height, LT, LNC and LCC gradually decreased with increasing elevation (Figure S2). In detail, height showed a decreasing overall pattern with elevation in the responses of most species, meaning higher elevations resulted in shorter plants. Although trends varied widely among species, more than half of the species showed a decreasing trend, resulting in a significant decline in LT and LCC with elevation. However, there were few species with a decreasing trend, but overall LNC decreased significantly along elevation, probably due to interspecific differences. Moreover, Differences between species trends may account for the non-significant relationship between SLA, LDMC and elevation (Figure S3-S8).

Relationships between functional traits and elevation for species with different elevational distributions

Elevation distributions, namely EP and SR, largely influences the relationship between functional traits and elevation. There were significant interactions between elevation and EP for traits such as plant height, LDMC, LT, LNC and LCC. Plant height, LT, and LNC varied less along the elevational gradients in high EP species, whereas low EP plants showed greater variation in functional trait values (Fig. 1). Among them, plant height and LT of low EP species decreased at higher elevations, while LNC increased. On the contrary, LDMC of high EP species decreased gradually along the elevation, while that of low EP species remained at a low level. In addition, LCC was highest and lowest in median elevation for high and low EP species, respectively.

Fig. 1
figure 1

Relationship between functional traits and elevation, as influenced by species elevational preference (EP). Functional traits: Plant height, specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (LT), leaf nitrogen content (LNC) and leaf carbon content (LCC). Shade areas are the 95% confidence intervals. cR2 represents conditional R2 value, and mR2 represents marginal R2 value. Trait values and elevation were standardized

The interactions of elevation and SR had significant effects on LDMC, LT, LNC and LCC. LT, LNC and LCC varied less along the elevational gradients in wide SR species, whereas narrow SR plants displayed a greater variation in functional trait values (Fig. 2). Wide and narrow species had opposite trends in LDMC. As the elevation increased, the LDMC of wide SR species decreased, and that of narrow SR species increased.

Fig. 2
figure 2

Relationship between functional traits and elevation, as influenced by species’ elevational range (SR). Functional traits: Plant height, specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (LT), leaf nitrogen content (LNC) and leaf carbon content (LCC). Shade areas are the 95% confidence intervals. cR2 represents conditional R2 value, and mR2 represents marginal R2 value. Trait values and elevation were standardized

Species’ elevational distributions

Using EP as the horizontal axis and SR as the vertical axis, we plotted elevational distributions of desert species with observations greater than 10 times. The clustering of plant species with similar elevational distribution characteristics could be found in the Figure S9. LN Group had 8 species and was characterized by low EP and narrow SR (0< EP<0.3, 0< SR<0.3), revealing that these species prefer lower elevations and only distributed at lower elevations. All species in this group were shrubs. We divided the 4 species into LW Group with low EP and wide SR (0< EP<0.3, 0.7< SR<1.0), indicating that these species prefer lower elevations but have a broader elevational distribution. This group consisted mainly of shrubs and forbs. Four species were classified into the HN group with high EP and narrow SR (0.7< EP<1.0, 0< SR<0.3), considering that these species were exclusive to high-elevation habitats. In this group, there were only two functional types (forbs and shrubs), and the species with the largest proportion were forbs. Ultimately, We classified 7 species into HW group with high EP and wide SR (0.7< EP<1.0, 0.7< SR<1.0), showing that these species prefer higher elevations but have a broader elevational distribution. In HW group, the most represented species were forbs and graminoids (Figure S9; Table S2).

Sources of variation in functional traits

We found that differences within species explained trait variation of plant height in low EP desert plants (Fig. 3a). Both total and intraspecific variation in plant height decreased significantly with increasing EP, and the total variation decreased with increasing SR (Figure S10 a; Figure S11 a). The intraspecific variation of SLA and LDMC showed a high proportion of the total trait variation accounting for an average of 63.68% and 52.96% of total variation, respectively, and the total variation in SLA increased with increasing SR (Fig. 3b c; Figure S11 b). LT had large interspecific variation, and total variation in LT decreased significantly with increasing EP and SR (Fig. 3d; Figure S10 d; Figure S11 d). LNC of wide SR species had a large intraspecific variation, and interspecific variation in LNC increased with increasing EP (Fig. 3e; Figure S10 e). Moreover, LCC in HW group exhibited large proportion of intraspecific variation, and interspecific variation in LCC decreased with increasing SR (Fig. 3f; Figure S11 f).

Fig. 3
figure 3

Variance decomposition of height, specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (LT), leaf nitrogen content (LNC) and leaf carbon content (LCC) measured across species with different elevational distributions. LN: species with low elevational preference and narrow species’ range; LW: species with low elevational preference and wide species’ range; HN: species with high elevational preference and narrow species’ range; HW: species with high elevational preference and wide species’ range. Colored bars show proportion of total trait variance (‘% trait Var’) while black bar shows absolute amount of variance (‘Tot trait Var’)

For species presenting these four distributions, we took into account patterns of intraspecific variation over elevation. We found that the significant negative feedback of plant height with elevation was reflected at the level of individual species for most species in all groups except the HW species group (Figure S12 a). SLA and LCC had significant intraspecific trends only for species in the LN and HW groups with the mostly decrease trend (Figure S12 b f). Most species with significant intraspecific variation trends in LDMC showed negative responses to elevation (Figure S12 c). LT of species in the LN group showed more constraint, with half of the species showing an increase trend and a quarter of the species a decrease trend (Figure S12 d). LNC showed the most consistent variance constraint with elevation in HW group, with half of the species showing higher trait values at higher elevation (Figure S12 e).

Relationships between functional traits and environmental factors

The relationship between functional traits of desert plants and environmental factors was further analyzed by RDA analysis, and the results showed that the cumulative explanation rates of the first two axes were 89.23%, 78.24%, 83.31% and 95.34%, respectively (Fig. 4). The first two axes could reflect the relationship well, and were mainly determined by axis I. In detail, temperature, precipitation and elevation provided a better explanation for the variation in functional traits in the LN group (Table S3). Environmental factors explain more about plant height, LCC and LT in the LN group. Temperature, elevation, precipitation, SNC, sand content, SCC and EC significantly affected (p < 0.05) the differences in functional traits of desert plants in the LW group. LCC, LNC and height of the LW group were better explained by environmental factors (Fig. 4; Table S3). In the HN group, temperature, elevation, sand content, SNC and pH were determined to be significant environmental factors (p < 0.05) affecting SLA, LNC and LCC of desert plants (Table S3). Temperature, elevation and sand content had significant effects on the variation of functional traits in the HW group, especially for SLA, LDMC and LNC.

Fig. 4
figure 4

Redundancy analysis of functional traits and environmental factors across species with different elevational distributions. Functional traits: plant height, specific leaf area (SLA), leaf dry matter content (LDMC), leaf thickness (LT), leaf nitrogen content (LNC) and leaf carbon content (LCC). Environmental factors: elevation, annual mean temperature (Temp), annual mean precipitation (Prec), soil pH (pH), soil electrical conductivity (EC), soil clay content (Clay), soil sand content (Sand), soil nitrogen content (SNC) and soil carbon content (SCC). LN: species with low elevational preference and narrow species’ range; LW: species with low elevational preference and wide species’ range; HN: species with high elevational preference and narrow species’ range; HW: species with high elevational preference and wide species’ range. Red lines indicate functional traits, and blue lines indicate environmental factors

Discussion

Our results demonstrated that elevational distributions affected trait shifts. These shifts were manifested in both trait values and variance portioning. The significant value changes in plant height and LT over elevation were mainly driven by species with lower EP, while LDMC was mainly driven by species with higher EP. Higher EP species exhibit lower trait variation compared to lower EP species, which may have a reduced potential to respond to environmental gradients [1]. Moreover, intraspecific trait variations of species with different elevational distributions demonstrated different adaptation strategies corresponding to elevation changes. We provide preliminary evidence that elevation, temperature and precipitation were the main factors influencing functional traits in lower EP species, while higher EP species were not influenced by precipitation.

Variation in functional traits across a wide elevational gradient

A number of factors, including temperature, precipitation, solar radiation, and atmospheric pressure, can influence the patterns of functional traits along elevational gradients [26]. Plant height, LT, LNC and LCC significantly decreased with increasing elevation. These traits of plant growth and resource utilization indicate a survival strategy of desert plants along a wide elevational gradient. Variations in plant height are supported by previous studies [12, 21]. Long documented, the negative correlation between plant height and elevation has been considered to be linked to community assembly and plastic variation of plants [12, 27]. For example, plant communities at higher elevations have lower height compared to those at lower elevations [28]. Individuals of the same species at higher elevations tend to be shorter than those at lower elevations, according to homogeneous garden trials, indicating intraspecific adaptation of plant size to elevation [12].

Variations in leaf size along elevational gradients are determined by different climatic factors and soil conditions [29]. It is generally accepted that the leaves become smaller and thicker with increasing elevation resulting in low SLA, high LDMC and high LT [4, 25, 30]. This can be explained in terms of both water and heat. Firstly, because leaf size is a key factor in controlling evapotranspiration and is used as a proxy for energy and water balances [29, 31], variations in leaf size caused by rising elevation may be a result of feedback from soil moisture [32]. In this experiment, both LT and SLA were strongly correlated with the precipitation. Secondly, variations in leaf size with elevation probably reflect the divergence in temperature between the day and night [4], with accumulation of matter in leaves [33]. Our results showed that the SLA in the LW group and the LDMC in the HW group were strongly related to temperature, which confirmed this opinion. Consequently, the variations in these traits may reflect the widespread and pervasive role of water-heat exchange in influencing plant plasticity [4]. However, our results suggested that the overall patterns of SLA and LDMC along elevation were not significant. This may be due to the complexity of trait variation within different species. Results from a study that measured SLA for 21 species at different elevational ranges found that SLA increased across elevational gradients [34]. This pattern may indicate that plants in higher elevations have lager photosynthetic capacity to ensure rapid carbon uptake [34, 35].

In line with previous studies, LNC and LCC decreased with increasing elevation mainly due to temperature and precipitation decreased [13, 36]. This result is consistent with the plant temperature-physiological theory, which postulates that plants’ metabolic activity slows down in the cold [12]. Transpiration, along with feedback from the soil and atmosphere, is the primary factor influencing changes in leaf nutrient concentrations [37]. These directional shifts in traits along elevational gradients collectively suggest that plants adjust their height, leaf size and nutrient concentration to make up for decreased photosynthetic capacity in response to lower temperatures, reduced precipitation, rising solar radiation and increased atmospheric pressure [21].

Species occupying different elevational distributions

Different plants have different elevational preferences, which is partially reflected in their spatial distribution [1]. The EP and SR of desert species were considered as X and Y axes, and we identified four categories of desert plants. Species in LN group characterized by low EP and narrow SR, contain Alhagi camelorum, Anabasis brevifolia, Halostachys caspica, Haloxylon ammodendron, Kalidium foliatum, Nitraria tangutorum, Tamarix chinensis and Zygophyllum xanthoxylum. These species, mostly drought and salinity tolerant plants, usually have small and tough leaves to conserve water and prevent transpiration [38], thus tending to colonize in specific areas, such as the saline and arid lands. LW group contains drought-tolerant species, and it is also characterized by a high resistance to harsh environments and high reproductive capacity, thus having a relatively wide distribution. Species in HN group with high EP and narrow SR, contain Asteraceae wellbyi, Christolea crassifolia, Oxytropis microphylla and Stellera chamaejasme. These species are cold tolerant and grow at higher elevations. Especially, Asteraceae wellbyi is endemic to Tibet. Species in the HW group are highly adaptable and barrenness tolerant. The wide distribution may be due to their high vigor and diverse reproduction modes, which enable them to occupy the ground quickly, as well as a well-developed fibrous root system and some physiological and ecological characteristics essential for adaptation to environmental stress [39].

Trait variation in species with different elevational distributions

For species that inhabit various elevational distributions, we would anticipate different trait values and variation portioning. First of all, we considered functional trait values over elevation for species with different elevational distributions. Our study suggested that the patterns of trait change with elevation depended on the EP and SR. In particular, the main reason for the decrease in plant height and LT with elevation is the species with a low EP. In contrast, elevation-induced LDMC reduction was primarily caused by high EP species. These results suggest that species with different elevational preferences may have different strategies for functional trait variation in response to environmental changes. Desert plants with a high EP may be subject to more abiotic stressors and less interspecific competition than species with a low EP, which may help them stick to their conservative growth strategy of staying small [1]. Desert plants with a low EP tend to be shrubs with succulent leaves that generally have lower LDMC to withstand drought [25, 40]. Most importantly, the inconsistency in the relationship between LDMC and elevation may stem from the species with different SR. In a narrow SR, the changes in LDMC were consistent with most studies [4, 25, 30], suggesting that hydrothermal conditions play an important role in the trait response process [4, 32]. Over a wide SR, however, variation in LDMC may be caused primarily by characteristics of different species. As mentioned before, species of LN group are mainly drought-tolerant shrubs and salt plants, and species of HN group are mainly cold-tolerant graminoids.

Trait variation across elevational gradients may also be a means by which desert species convey their varying preferences for habitat. We offered preliminary proof that desert species with various elevational distributions have diverse patterns of trait variation portioning. According to our results, structural traits of high EP plants show relatively little variation with elevation, which may point to a higher capacity for adapting to environmental changes [1, 34, 41]. While species with a high EP show greater interspecific variation in nutrient trait values with elevation. Given the potential effects of climate change, plant species with relatively high trait variability may be more adaptable to different environmental situations than those with relatively low trait variability [1]. Moreover, intraspecific variations in LNC and LCC were higher in species with a wide SR than in species with a narrow SR, highlighting the high trait plasticity of plant carbon and nitrogen content in widely distributed species. The global mean value of intraspecific variation was 25% [42], and the contribution of intraspecific variation to total trait variation was either equal to or greater than this value for the different groups in this study, despite tough environmental conditions and wide species ranges in Tibetan Plateau. However, intraspecific trait variation varied considerably for different traits, e.g., intraspecific trait variation in SLA and LDMC accounted for about 50% of the total variation, which was consistent with previous studies [30]. Our results provide evidence that the distribution of species along environmental gradients is constrained by intraspecific trait variation [34]. Taken together, the results of this study revealed that, plant establishment and adaption success under varying environmental conditions can be attributed to differences in functional traits [6, 29].

Response of functional traits to environmental factors

In addition to the influence of the species’ elevational distributions of the plant itself, external environmental factors are important for variation in plant functional traits. Temperature and precipitation are important determinants of the regional climate type and significantly affect plant growth and development. We discovered that species with a high EP faced constraints from elevation and temperature, those with a low EP mostly derived their functional traits from elevation temperature and precipitation. Species with high EP are mostly located in the alpine desert of the Tibetan plateau, where the presence of cold climatic conditions and permafrost prevent plants from efficiently utilizing water [43]. Therefore, changes in the functional traits of desert plants are not significantly influenced by precipitation in this area. Our results suggest that LNC was negatively correlated with temperature, which is consistent with previous studies in field surveys and simulated controlled experiments [44, 45]. This may because that high nitrogen content at low temperatures is needed to compensate for the reduced biochemical efficiency caused by the reduction of high-nitrogen enzyme activity [46]. Alternatively, warmer climate accelerates the plant growth process, thus diluting LNC [45]. Precipitation negatively correlates with SLA for low EP species, suggesting the plant adaptation strategies to maximize carbon income and minimize water consumption under drought stress [4, 47]. In this study, soil properties such as sand content and SNC have important effects on the formation of functional traits in desert plants, as soil is a material and energy source for plant growth and development.

Conclusions

We discovered that desert plant species displayed different trait trends over elevation, and that these associations relied on the elevational distributions (elevational preferences and species’ ranges) of the individual species. In particular, species with lower elevational preferences expressed higher trait variation in structure trait than those with higher elevational preferences. It was suggested by the increased intraspecific variation of SLA and LDMC that these species may be better adapted to biotic and abiotic changes. Plant species with lower elevational ranges have trait-elevation connections that are widely applicable globally, but LDMC at wider elevational ranges show opposite trends, suggesting that interspecific variation plays an important role in size-related traits at large scales. Most importantly, the main controlling factors of functional traits differed among species with different elevational distributions. Our experiments provide preliminary evidence that desert species with different elevational distributions have different trait distribution patterns and adaptation mechanisms.

Data availability

Data is provided within the manuscript or supplementary information files.

Abbreviations

EC:

Electrical conductivity

EP:

Elevational preference

LCC:

Leaf carbon content

LDMC:

Leaf dry matter content

LNC:

Leaf nitrogen content

LT:

Leaf thickness

RDA:

Redundancy analysis

SCC:

Soil carbon content

SLA:

Specific leaf area

SNC:

Soil nitrogen content

SR:

Species’ range

References

  1. Rixen C, Wipf S, Rumpf SB, Giejsztowt J, Millen J, Morgan JW, Nicotra AB, Venn S, Zong S, Dickinson KJM, Freschet GT, Kurzböck C, Li J, Pan H, Pfund B, Quaglia E, Su X, Wang W, Wang X, Yin H, Deslippe JR. Intraspecific trait variation in alpine plants relates to their elevational distribution. J Ecol. 2022;110:860–75. https://doi.org/10.1111/1365-2745.13848

    Article  Google Scholar 

  2. Enquist BJ, Norberg J, Bonser SP, Violle C, Webb CT, Henderson A, Sloat LL, Savage VM. Scaling from traits to ecosystems, in: advances in Ecological Research. Elsevier; 2015. pp. 249–318. https://doi.org/10.1016/bs.aecr.2015.02.001

  3. Mcgill B, Enquist B, Weiher E, Westoby M. Rebuilding community ecology from functional traits. Trends Ecol Evol. 2006;21:178–85. https://doi.org/10.1016/j.tree.2006.02.002

    Article  PubMed  Google Scholar 

  4. Wright IJ, Dong N, Maire V, Prentice IC, Westoby M, Díaz S, Gallagher RV, Jacobs BF, Kooyman R, Law EA, Leishman MR, Niinemets Ü, Reich PB, Sack L, Villar R, Wang H, Wilf P. Global climatic drivers of leaf size. Science. 2017;357:917–21. https://doi.org/10.1126/science.aal4760

    Article  CAS  PubMed  Google Scholar 

  5. Niu S, Classen AT, Luo Y. Functional traits along a transect. Funct Ecol. 2018;32:4–9. https://doi.org/10.1111/1365-2435.13023

    Article  Google Scholar 

  6. Wang J, Hu A, Meng F, Zhao W, Yang Y, Soininen J, Shen J, Zhou J. Embracing mountain microbiome and ecosystem functions under global change. New Phytol. 2022;234:1987–2002. https://doi.org/10.1111/nph.18051

    Article  PubMed  Google Scholar 

  7. Reich PB. The world-wide ‘fast-slow’ plant economics spectrum: a traits manifesto. J Ecol. 2014;102:275–301. https://doi.org/10.1111/1365-2745.12211

    Article  Google Scholar 

  8. Classen A, Steffan-Dewenter I, Kindeketa WJ, Peters MK. Integrating intraspecific variation in community ecology unifies theories on body size shifts along climatic gradients. Funct Ecol. 2017;31:768–77. https://doi.org/10.1111/1365-2435.12786

    Article  Google Scholar 

  9. Allen JA. The influence of physical conditions in the genesis of species. Radical Reviews. 1877;1:108–40. https://doi.org/10.1038/scientificamerican05181907-26226supp

    Article  Google Scholar 

  10. Cornelissen JHC, van Bodegom PM, Aerts R, Callaghan TV, van Logtestijn RSP, Alatalo J, Stuart Chapin F, Gerdol R, Gudmundsson J, Gwynn-Jones D, Hartley AE, Hik DS, Hofgaard A, Jónsdóttir IS, Karlsson S, Klein JA, Laundre J, Magnusson B, Michelsen A, Molau U, Onipchenko VG, Quested HM, Sandvik SM, Schmidt IK, Shaver GR, Solheim B, Soudzilovskaia NA, Stenström A, Tolvanen A, Totland Ø, Wada N, Welker JM, Zhao X, Team† MOL. Global negative vegetation feedback to climate warming responses of leaf litter decomposition rates in cold biomes. Ecol Lett. 2007;10:619–27. https://doi.org/10.1111/j.1461-0248.2007.01051.x

    Article  PubMed  Google Scholar 

  11. Freschet GT, Roumet C, Comas LH, Weemstra M, Bengough AG, Rewald B, Bardgett RD, De Deyn GB, Johnson D, Klimešová J, Lukac M, McCormack ML, Meier IC, Pagès L, Poorter H, Prieto I, Wurzburger N, Zadworny M, Bagniewska-Zadworna A, Blancaflor EB, Brunner I, Gessler A, Hobbie SE, Iversen CM, Mommer L, Picon‐Cochard C, Postma JA, Rose L, Ryser P, Scherer‐Lorenzen M, Soudzilovskaia NA, Sun T, Valverde‐Barrantes OJ, Weigelt A, York LM, Stokes A. Root traits as drivers of plant and ecosystem functioning: current understanding, pitfalls and future research needs. New Phytol. 2021;232:1123–58. https://doi.org/10.1111/nph.17072

    Article  PubMed  Google Scholar 

  12. Halbritter AH, Fior S, Keller I, Billeter R, Edwards PJ, Holderegger R, Karrenberg S, Pluess AR, Widmer A, Alexander JM. Trait differentiation and adaptation of plants along elevation gradients. J Evol Biol. 2018;31:784–800. https://doi.org/10.1111/jeb.13262

    Article  PubMed  Google Scholar 

  13. Midolo G, De Frenne P, Hölzel N, Wellstein C. Global patterns of intraspecific leaf trait responses to elevation. Glob Change Biol. 2019;25:2485–98. https://doi.org/10.1111/gcb.14646

    Article  Google Scholar 

  14. MacLean SA, Beissinger SR. Species’ traits as predictors of range shifts under contemporary climate change: a review and meta-analysis. Glob Change Biol. 2017;23:4094–105. https://doi.org/10.1111/gcb.13736

    Article  Google Scholar 

  15. Shipley B, De Bello F, Cornelissen JHC, Laliberté E, Laughlin DC, Reich PB. Reinforcing loose foundation stones in trait-based plant ecology. Oecologia. 2016;180:923–31. https://doi.org/10.1007/s00442-016-3549-x

    Article  PubMed  Google Scholar 

  16. Giejsztowt J, Classen AT, Deslippe JR. Climate change and invasion may synergistically affect native plant reproduction. Ecology. 2020;101. https://doi.org/10.1002/ecy.2913

  17. Siefert A, Violle C, Chalmandrier L, Albert CH, Taudiere A, Fajardo A, Aarssen LW, Baraloto C, Carlucci MB, Cianciaruso MV, Dantas L, Bello V, Duarte F, Fonseca LDS, Freschet CR, Gaucherand GT, Gross S, Hikosaka N, Jackson K, Jung B, Kamiyama V, Katabuchi C, Kembel M, Kichenin SW, Kraft E, Lagerström NJB, Bagousse-Pinguet A, Li YL, Mason Y, Messier N, Nakashizuka J, Overton T, Peltzer JMC, Pérez‐Ramos DA, Pillar IM, Prentice VD, Richardson HC, Sasaki S, Schamp T, Schöb BS, Shipley C, Sundqvist B, Sykes M, Vandewalle MT, Wardle M, D.A. A global meta‐analysis of the relative extent of intraspecific trait variation in plant communities. Ecol Lett. 2015;18:1406–19. https://doi.org/10.1111/ele.12508

    Article  PubMed  Google Scholar 

  18. Fajardo A, Piper FI. Intraspecific trait variation and covariation in a widespread tree species (Nothofagus Pumilio) in southern Chile. New Phytol. 2011;189:259–71. https://doi.org/10.1111/j.1469-8137.2010.03468.x

    Article  PubMed  Google Scholar 

  19. Anderegg LDL, Berner LT, Badgley G, Sethi ML, Law BE, HilleRisLambers J. Within-species patterns challenge our understanding of the leaf economics spectrum. Ecol Lett. 2018;21:734–44. https://doi.org/10.1111/ele.12945

    Article  PubMed  Google Scholar 

  20. Pei ZY, Ouyang H, Zhou CP, Xu XL. Carbon Balance in an Alpine Steppe in the Tibetan Plateau. J Integr Plant Biol. 2009;51:521–6. https://doi.org/10.1111/j.1744-7909.2009.00813.x

    Article  CAS  PubMed  Google Scholar 

  21. Sigdel SR, Liang E, Rokaya MB, Rai S, Dyola N, Sun J, Zhang L, Zhu H, Chettri N, Chaudhary RP, Camarero JJ, Peñuelas J. Functional traits of a plant species fingerprint ecosystem productivity along broad elevational gradients in the Himalayas. Funct Ecol. 2023;37:383–94. https://doi.org/10.1111/1365-2435.14226

    Article  CAS  Google Scholar 

  22. Qiu J. China: the third Pole. Nature. 2008;454:393–6. https://doi.org/10.1038/454393a

    Article  CAS  PubMed  Google Scholar 

  23. Cheng G, Wu T. Responses of permafrost to climate change and their environmental significance. Tibetan Plateau J Geophys Res. 2007;112:F02S03. https://doi.org/10.1029/2006JF000631

  24. Maharjan SK, Sterck FJ, Dhakal BP, Makri M, Poorter L. Functional traits shape tree species distribution in the Himalayas. J Ecol. 2021;109:3818–34. https://doi.org/10.1111/1365-2745.13759

    Article  Google Scholar 

  25. Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P, Bret-Harte MS, Cornwell WK, Craine JM, Gurvich DE, Urcelay C, Veneklaas EJ, Reich PB, Poorter L, Wright IJ, Ray P, Enrico L, Pausas JG, De Vos AC, Buchmann N, Funes G, Quétier F, Hodgson JG, Thompson K, Morgan HD, Ter Steege H, Sack L, Blonder B, Poschlod P, Vaieretti MV, Conti G, Staver AC, Aquino S, Cornelissen JHC. New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61:167. https://doi.org/10.1071/BT12225

    Article  Google Scholar 

  26. Körner C. The use of ‘altitude’ in ecological research. Trends Ecol Evol. 2007;22:569–74. https://doi.org/10.1016/j.tree.2007.09.006

    Article  PubMed  Google Scholar 

  27. Read QD, Moorhead LC, Swenson NG, Bailey JK, Sanders NJ. Convergent effects of elevation on functional leaf traits within and among species. Funct Ecol. 2014;28:37–45. https://doi.org/10.1111/1365-2435.12162

    Article  Google Scholar 

  28. Kergunteuil A, Descombes P, Glauser G, Pellissier L, Rasmann S. Plant physical and chemical defence variation along elevation gradients: a functional trait-based approach. Oecologia. 2018;187:561–71. https://doi.org/10.1007/s00442-018-4162-y

    Article  PubMed  Google Scholar 

  29. Díaz S, Kattge J, Cornelissen JHC, Wright IJ, Lavorel S, Dray S, Reu B, Kleyer M, Wirth C, Colin Prentice I, Garnier E, Bönisch G, Westoby M, Poorter H, Reich PB, Moles AT, Dickie J, Gillison AN, Zanne AE, Chave J, Joseph Wright S, Sheremet’ev SN, Jactel H, Baraloto C, Cerabolini B, Pierce S, Shipley B, Kirkup D, Casanoves F, Joswig JS, Günther A, Falczuk V, Rüger N, Mahecha MD, Gorné LD. The global spectrum of plant form and function. Nature. 2016;529:167–71. https://doi.org/10.1038/nature16489

    Article  CAS  PubMed  Google Scholar 

  30. Thomas HJD, Bjorkman AD, Myers-Smith IH, Elmendorf SC, Kattge J, Diaz S, Vellend M, Blok D, Cornelissen JHC, Forbes BC, Henry GHR, Hollister RD, Normand S, Prevéy JS, Rixen C, Schaepman-Strub G, Wilmking M, Wipf S, Cornwell WK, Beck PSA, Georges D, Goetz SJ, Guay KC, Rüger N, Soudzilovskaia NA, Spasojevic MJ, Alatalo JM, Alexander HD, Anadon-Rosell A, Angers-Blondin S, Te Beest M, Berner LT, Björk RG, Buchwal A, Buras A, Carbognani M, Christie KS, Collier LS, Cooper EJ, Elberling B, Eskelinen A, Frei ER, Grau O, Grogan P, Hallinger M, Heijmans MMPD, Hermanutz L, Hudson JMG, Johnstone JF, Hülber K, Iturrate-Garcia M, Iversen CM, Jaroszynska F, Kaarlejarvi E, Kulonen A, Lamarque LJ, Lantz TC, Lévesque E, Little CJ, Michelsen A, Milbau A, Nabe-Nielsen J, Nielsen SS, Ninot JM, Oberbauer SF, Olofsson J, Onipchenko VG, Petraglia A, Rumpf SB, Shetti R, Speed JDM, Suding KN, Tape KD, Tomaselli M, Trant AJ, Treier UA, Tremblay M, Venn SE, Vowles T, Weijers S, Wookey PA, Zamin TJ, Bahn M, Blonder B, Van Bodegom PM, Bond-Lamberty B, Campetella G, Cerabolini BEL, Chapin FS, Craine JM, Dainese M, Green WA, Jansen S, Kleyer M, Manning P, Niinemets Ü, Onoda Y, Ozinga WA, Peñuelas J, Poschlod P, Reich PB, Sandel B, Schamp BS, Sheremetiev, S.N., De Vries. F.T., 2020. Global plant trait relationships extend to the climatic extremes of the tundra biome. Nat Commun 11, 1351. https://doi.org/10.1038/s41467-020-15014-4

  31. Van Der Plas F, Schröder-Georgi T, Weigelt A, Barry K, Meyer S, Alzate A, Barnard RL, Buchmann N, De Kroon H, Ebeling A, Eisenhauer N, Engels C, Fischer M, Gleixner G, Hildebrandt A, Koller-France E, Leimer S, Milcu A, Mommer L, Niklaus PA, Oelmann Y, Roscher C, Scherber C, Scherer-Lorenzen M, Scheu S, Schmid B, Schulze E-D, Temperton V, Tscharntke T, Voigt W, Weisser W, Wilcke W, Wirth C. Plant traits alone are poor predictors of ecosystem properties and long-term ecosystem functioning. Nat Ecol Evol. 2020;4:1602–11. https://doi.org/10.1038/s41559-020-01316-9

    Article  PubMed  Google Scholar 

  32. Anderegg WRL, Trugman AT, Bowling DR, Salvucci G, Tuttle SE. Plant functional traits and climate influence drought intensification and land–atmosphere feedbacks. Proc Natl Acad Sci U S A. 2019;116:14071–6. https://doi.org/10.1073/pnas.1904747116

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Poorter L, Wright SJ, Paz H, Ackerly DD, Condit R, Ibarra-Manríquez G, Harms KE, Licona JC, Martínez-Ramos M, Mazer SJ, Muller-Landau HC, Peña-Claros M, Webb CO, Wright IJ. ARE FUNCTIONAL TRAITS GOOD PREDICTORS OF DEMOGRAPHIC RATES? EVIDENCE FROM FIVE NEOTROPICAL FORESTS. Ecology. 2008;89:1908–20. https://doi.org/10.1890/07-0207.1

    Article  CAS  PubMed  Google Scholar 

  34. Sides CB, Enquist BJ, Ebersole JJ, Smith MN, Henderson AN, Sloat LL. Revisiting Darwin’s hypothesis: does greater intraspecific variability increase species’ ecological breadth? Am J Bot. 2014;101:56–62. https://doi.org/10.3732/ajb.1300284

    Article  PubMed  Google Scholar 

  35. Kudo G. Intraspecific variation of leaf traits in several deciduous species in relation to length of growing season. Écoscience. 1996;3:483–9. https://doi.org/10.1080/11956860.1996.11682367

    Article  Google Scholar 

  36. Drollinger S, Müller M, Kobl T, Schwab N, Böhner J, Schickhoff U, Scholten T. Decreasing nutrient concentrations in soils and trees with increasing elevation across a treeline ecotone in Rolwaling Himal. Nepal J Mt Sci. 2017;14:843–58. https://doi.org/10.1007/s11629-016-4228-4

    Article  Google Scholar 

  37. Westerband AC, Knight TM, Barton KE. Intraspecific trait variation and reversals of trait strategies across key climate gradients in native hawaiian plants and non-native invaders. Ann Botany. 2021;127:553–64. https://doi.org/10.1093/aob/mcaa050

    Article  CAS  Google Scholar 

  38. Prieto I, Kikvidze Z, Pugnaire FI. Hydraulic lift: soil processes and transpiration in the Mediterranean leguminous shrub Retama sphaerocarpa (L.) Boiss. Plant Soil. 2010;329:447–56. https://doi.org/10.1007/s11104-009-0170-3

    Article  CAS  Google Scholar 

  39. Wellstein C, Poschlod P, Gohlke A, Chelli S, Campetella G, Rosbakh S, Canullo R, Kreyling J, Jentsch A, Beierkuhnlein C. Effects of extreme drought on specific leaf area of grassland species: a meta-analysis of experimental studies in temperate and sub-mediterranean systems. Glob Change Biol. 2017;23:2473–81. https://doi.org/10.1111/gcb.13662

    Article  Google Scholar 

  40. Wilson PJ, Thompson K, Hodgson JG. Specific leaf area and leaf dry matter content as alternative predictors of plant strategies. New Phytol. 1999;143:155–62. https://doi.org/10.1046/j.1469-8137.1999.00427.x

    Article  Google Scholar 

  41. Rumpf SB, Hülber K, Klonner G, Moser D, Schütz M, Wessely J, Willner W, Zimmermann NE, Dullinger S. 2018. Range dynamics of mountain plants decrease with elevation. Proc. Natl. Acad. Sci. U.S.A. 115, 1848–1853. https://doi.org/10.1073/pnas.1713936115

  42. Lavorel S, Garnier E. Predicting changes in community composition and ecosystem functioning from plant traits: revisiting the Holy Grail: plant response and effect groups. Funct Ecol. 2002;16:545–56. https://doi.org/10.1046/j.1365-2435.2002.00664.x

    Article  Google Scholar 

  43. Tian H, Chen G, Liu M, Zhang C, Sun G, Lu C, Xu X, Ren W, Pan S, Chappelka A. Model estimates of net primary productivity, evapotranspiration, and water use efficiency in the terrestrial ecosystems of the southern United States during 1895–2007. For Ecol Manag. 2010;259:1311–27. https://doi.org/10.1016/j.foreco.2009.10.009

    Article  Google Scholar 

  44. Navarro T, Oualidi EI, Taleb J, Pascual MS, Cabezudo V, Milla B, R. Leaf patterns, leaf size and ecologically related traits in high Mediterranean mountain on the Moroccan high Atlas. Plant Ecol. 2010;210:275–90. https://doi.org/10.1007/s11258-010-9756-3

    Article  Google Scholar 

  45. Day TA, Ruhland CT, Xiong F. Warming increases aboveground plant biomass and C stocks in vascular-plant‐dominated Antarctic tundra. Glob Change Biol. 2008;14:1827–43. https://doi.org/10.1111/j.1365-2486.2008.01623.x

    Article  Google Scholar 

  46. Wang G. Leaf trait co-variation, response and effect in a chronosequence. J Veg Sci. 2007;18:563–70. https://doi.org/10.1111/j.1654-1103.2007.tb02570.x

    Article  Google Scholar 

  47. Westoby M, Falster DS, Moles AT, Vesk PA, Wright IJ. Plant ecological strategies: some leading dimensions of variation between species. Annu Rev Ecol Syst. 2002;33:125–59. https://doi.org/10.1146/annurev.ecolsys.33.010802.150452

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was financially supported by the Second Tibetan Plateau Scientific Expedition and Research program (No. 2019QZKK0305), National Natural Science Foundation of China (No. 42071140) and CAS “Light of West China” Program (No. E129050301).

Author information

Authors and Affiliations

Authors

Contributions

X.Z. designed this experiment. X.Li., S.W., P.L., P.Y. and M.C. contributed significantly to analysis and manuscript preparation. Y.H. performed the data analyses and wrote the manuscript. All authors reviewed the manuscript.

Corresponding author

Correspondence to Xiaoan Zuo.

Ethics declarations

Ethics approval and consent to participate

This article does not contain any studies with human participants or animals and did not involve any endangered or protected species. The plant materials sampled and experiments performed in this research were conducted in accordance with local legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, Y., Li, X., Wang, S. et al. Patterns and driving factors of functional traits of desert species with different elevational distributions in the Tibetan Plateau and adjacent areas. BMC Plant Biol 24, 371 (2024). https://doi.org/10.1186/s12870-024-05080-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05080-x

Keywords