Yicheng He1,2‡
Tingting Li1‡Ruiyang Zhang1
Jinsong Wang1†
Juntao Zhu1†Yang Li1Xinli Chen3†Junxiao Pan1Ying Shen1,4†Furong Wang1
Jingwen Li2*
Dashuan Tian1,5*†- 1Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences (CAS), Beijing, China
- 2School of Ecology and Nature Conservation, Beijing Forestry University, Beijing, China
- 3Faculty of Natural Resources Management, Lakehead University, Thunder Bay, ON, Canada
- 4Key Laboratory of Animal Ecology and Conservation Biology, China Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China
- 5College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, China
Aims: Leaf chlorophyll (Chl) is a fundamental component and good proxy for plant photosynthesis. However, we know little about the large-scale patterns of leaf Chl and the relative roles of current environment changes vs. plant evolution in driving leaf Chl variations.
Locations: The east to west grassland transect of the Tibetan Plateau.
Methods: We performed a grassland transect over 1,600 km across the Tibetan Plateau, measuring leaf Chl among 677 site-species.
Results: Leaf Chl showed a significantly spatial pattern across the grasslands in the Tibetan Plateau, decreasing with latitude but increasing with longitude. Along with environmental gradient, leaf Chl decreased with photosynthetically active radiation (PAR), but increased with water availability and soil nitrogen availability. Furthermore, leaf Chl also showed significant differences among functional groups (C4 > C3 species; legumes < non-legume species), but no difference between annual and perennial species. However, we surprisingly found that plant evolution played a dominant role in shaping leaf Chl variations when comparing the sum and individual effects of all the environmental factors above. Moreover, we revealed that leaf Chl non-linearly decreased with plant evolutionary divergence time. This well-matches the non-linearly increasing trend in PAR or decreasing trend in temperature during the geological time-scale uplift of the Tibetan Plateau.
Main Conclusion: This study highlights the dominant role of plant evolution in determining leaf Chl variations across the Tibetan Plateau. Given the fundamental role of Chl for photosynthesis, these results provide new insights into reconsidering photosynthesis capacity in alpine plants and the carbon cycle in an evolutionary view.
Introduction
Leaf chlorophyll (Chl) is fundamental for harvesting light for plant photosynthesis (Croft et al., 2017). The importance of leaf Chl for photosynthesis has been comprehensively evaluated at both ecosystem and species levels. Across terrestrial ecosystems, numerous studies, using remote sensing techniques, reveal the close relationships between leaf Chl fluorescence and photosynthetic processes, such as rubisco concentration, maximum rate of carboxylation (Vcmax), and electron transport (Jmax) (Carswell et al., 2000; Kattge et al., 2009; Luo et al., 2019). At the species level, leaf Chl has been mostly studied at the local scale (Croft et al., 2017). However, across contrasting environments, species-level Chl might show great variations due to the difference in plant evolutionary histories, climate, soil, and functional groups (Li et al., 2018a,b; Zhang et al., 2020). To date, we still know little about the large-scale patterns and key drivers of species-level leaf Chl variations (Ryu et al., 2019).
Plant evolutionary history may be crucial to drive leaf Chl variations (Schmerler et al., 2012). Previously, leaf Chl is found to show evolutionary divergence among 823 woody species across forest ecosystems (Li et al., 2018a), denoting the potential role of evolutionary history in affecting leaf Chl. However, we do not know the direction of evolutionary effects on leaf Chl (Li et al., 2018a). Coincidently, the geological time-scale uplift of the Tibetan Plateau keeps pace with plant evolutionary history, providing an ideal opportunity to test how plant evolution regulates the direction of leaf Chl variations. Specifically, following the collision of the India and Eurasia continent (c. 50 million years ago), the Tibetan Plateau begins to uplift rapidly (Royden et al., 2008; Wang et al., 2008). This should cause directional changes in environmental conditions, especially for photosynthetically active radiation (PAR) and temperature. Based on the well-known correlations of PAR or temperature with altitude, PAR should gradually increase but temperature decreases during the uplift (Blumthaler et al., 1997). These environmental changes during plant evolutionary history may decrease the leaf Chl in newly evolved species across the Tibetan Plateau, due to the evolutionary adaptation to extremely high PAR (avoiding radiation damage) and cold environment (more resource allocation to plant survival) (Knapp and Smith, 2001).
In addition to plant evolution, current climate and soil nutrients might affect leaf Chl as well (Anjum et al., 2011; Genesio et al., 2020). Generally, plants tend to enhance leaf Chl in light-limited environments, helping to compete for the limiting light resource and sustain plant growth (Hallik et al., 2009; Li et al., 2018a; Genesio et al., 2020). For example, leaf Chl is found to be higher for trees' understory than in the canopy (Niinemets, 2010). Low precipitation normally causes soil water stress and air dryness, negatively influencing the elasticity of chloroplast membrane and the synthesis of leaf Chl (Rastogi et al., 2020; Xiong and Nadal, 2020). Soil nutrient availability may influence leaf Chl via altering plant nutrient uptake and use. For instance, leaf Chl is generally stimulated by more nitrogen (N) input (Cai et al., 2011; Guo et al., 2014; Rhein and Silva, 2016), because N is the most essential but limiting element to synthesize leaf Chl (Cai et al., 2011; Croft et al., 2017). Phosphorus (P) is also the basic element for plant photosynthesis, such as adenosine triphosphate (ATP), phosphodiester linkage (Lambers et al., 2011), and 5-aminolevulinic acid (ALA, a basic compound for Chl synthesis) (Masuda et al., 1996). Based on the importance of P for energy biosynthesis (ATP) and ALA, leaf Chl may vary with soil P availability.
Plant functional group (PFG) has the potential to determine leaf Chl variations (Zhang et al., 2020; Kelly et al., 2021). For example, C4 plants have the ability to increase intercellular CO2 concentration, which leads to a stronger photosynthetic capacity than C3 plants (Wang et al., 2011; Osborne et al., 2014). It is thus expected that C4 species might have a higher leaf Chl to utilize the higher CO2 concentration for maximizing photosynthesis. Among PFG with different life spans, annual plants tend to allocate more resources to aboveground growth, while perennial plants invest more in belowground growth to maintain the advantage in competing for soil resources in the long term (Onoda et al., 2017; Monroe et al., 2019). This suggests that annual species adopt a resource-use strategy, but perennial species with a resource-conservative strategy (Wright et al., 2004). This strategy difference likely results in a higher leaf Chl in annual than perennial species, which enables to increase the photosynthetic capacity of annual species for its fast growth. Legume species generally have advantages in obtaining soil N, possibly leading to higher leaf N and Chl than non-legume species (Guo et al., 2017; Zhang et al., 2020). Consequently, leaf Chl may be different among PFGs, due to the difference in the photosynthetic pathway, life span, and nitrogen-use strategy. In spite of these potential mechanisms above, we still do not know the relative importance of plant evolutionary history vs. current environment changes in shaping the patterns of leaf Chl across a range of contrasting environments.
Based on the reasoning above, the uplift of the Tibetan Plateau provides an ideal opportunity to explore the role of plant evolution in shaping leaf Chl variations (Xu et al., 2010; Chen et al., 2013). Furthermore, recent global changes (e.g., warming, precipitation change, and N deposition) have substantially impacted leaf Chl and photosynthetic processes across this area (Fu et al., 2014; Fu and Shen, 2016). Thus, we conducted a grassland transect survey (97 sites) of leaf Chl among common species (93 species from 27 families) and related variables across the Tibetan Plateau. Specifically, we addressed the following questions: (1) What are the regional patterns and key environmental control of leaf Chl concentration among alpine herbaceous species? (2) What is the role of plant evolution in affecting leaf Chl concentration relative to current environmental changes?
Materials and Methods
Site Description
We conducted a grassland transect survey of leaf Chl from the east to west of the Tibetan Plateau, covering a large gradient (1,600 km) of climate, soil nutrients, and ecosystem type. In total, we investigated 97 sites (Supplementary Figure 1), spanning a latitude range from 30.24 to 33.67°N, longitude from 80.36 to 92.13°E, and elevation from 4,181 to 5,016 m. The average annual PAR (1981–2010) is between 2461.8 and 2860.4 mol m−2, annual temperature between −5.37 and 0.99°C, and annual precipitation between 61 and 441 mm. Ecosystem types included alpine steppe, alpine meadow, and alpine desert grassland. Soil available N varied from 6.48 to 443.71 mg kg−1 and available P from 0.94 to 7.11 mg kg−1. Detailed geographical information of the sampling sites is presented in Supplementary Table 1.
Field Survey and Measurement
The field survey was carried out from July to August during the growing season in 2019 and 2020. According to local vegetation, in each site, three 5 m*5 m plots were randomly established in a 10 m*200 m area. In each sampling site, we randomly selected at least three plant individuals (3~6 replications) for each of the common species (above 80% of the local species) across three plots. Specifically, healthy leaf Chl was measured by the SPAD-502 Chl meter (Minolta, Osaka, Japan) (Uddling et al., 2007; Yuan et al., 2016). We measured the value of SPAD in the middle of the leaf, and leaf veins were excluded in each measurement. To more accurately estimate leaf Chl, the value was obtained by the mean of at least three readings in different positions of each individual. Finally, the mean value of leaf Chl in each site species was estimated at least by three plant individuals. In total, we measured 2,031 individuals among 93 herbaceous species from 27 families (many species occurring in multiple sites). All the species are grouped into legumes (10) vs. non-legume species (83), C4 (12) vs. C3 species (81), or annual (20) vs. perennial species (68).
Accordingly, three soil samples were randomly collected by a 7.5 cm-diameter soil auger at 0–10 cm depth in each plot. Then, we mixed three soil samples into one composite sample in each plot. The composite soil samples were screened by a 2-mm mesh to remove roots and stones (Pan et al., 2021). Totally, there were three composite soil samples at each site. A part of the soil was in cold storage by the portable refrigerator, and the other part was air-dried in labs. We measured soil available N (AN) by the alkali solution diffusion absorption method. Soil available P (AP) concentration was firstly extracted by the NaHCO3 solution (0.5 M) and then analyzed using the UV-visible spectrophotometer (Olsen, 1954).
Statistical Analyses
Before analyses, we obtained the data for mean annual precipitation (MAP) and temperature (MAT) from the Worldclim Database (http://www.worldclim.org), based on our location information. We then calculated humidity index (HI) as follows: HI = MAP / (MAT +10) (Martonne, 1926). PAR data were acquired from the dataset of “Global radiation, photosynthetically active radiation, and the diffuse components in China, 1981–2010” (Ren et al., 2018).
To obtain the divergence time of our sampling species, we first confirmed whether all the species (93 species from 25 families) have been included by The Plant List (http://www.theplantlist.org). Then, we established the phylogenetic tree using the “phytools” and “S.PhyloMaker” packages (Qian and Jin, 2016). The divergence time of each species was quantified by the branch length of each species in the phylogenetic tree (Qian and Jin, 2016). At last, we tested the strength of the phylogenetic signal for leaf Chl, using Blomberg's K and Pagel'λ K analyses with the “Picante” package (Blomberg et al., 2003).
Based on the data above, we first used the t-test analysis to compare the difference in leaf Chl between C3 and C4 species, annuals and perennial species, or legumes vs. non-legume species. Then, the linear or non-linear regression analyses were applied to analyze the pairwise relationships of leaf Chl with plant evolutionary divergence time (family and genus levels), climate (PAR and HI), soil nutrients (AN and AP), and location information (latitude, longitude, and altitude). All the regression analyses were performed by the “stats” package (Flynt and Dean, 2016).
Moreover, we performed the Structural Equation Model (SEM) to distinguish the direct and indirect effects of environmental factors and plant evolution on leaf Chl. Before the SEM analysis, we firstly established a conceptual model of potential casual pathways in driving leaf Chl, based on the traditional knowledge and the correlations between environmental factors (Supplementary Figure 2). Then, the model performance was assessed by the chi-square test (χ2 = 10.178, p = 0.12) with the “Lavaan” package (Rosseel, 2012). Finally, we analyzed and compared the standardized pathway coefficients of plant evolution and environmental factors on leaf Chl by summing all the direct and indirect effects. All statistical analysis was carried out in R software (version 3.6.3, R Foundation for Statistical Computing) and all the figures were made using the “ggplot2” package (Ginestet, 2011).
Results
The Regional Patterns and Environmental Controls of Leaf Chl
Across grasslands on the Tibetan Plateau, a significant spatial pattern of leaf Chl was detected among 677 site species. Leaf Chl was lower with increasing latitude or decreasing longitude (Figures 1A,B), but showed no significant trend with altitude (Figure 1C). Along with environmental gradients, leaf Chl reduced with PAR significantly (Figure 2A), but increased with HI (Figure 2B). Furthermore, leaf Chl was enhanced by higher soil N availability (Figure 2C), but had no response to soil P availability (Figure 2D).
Figure 1. The latitudinal (A), longitudinal (B), or elevation (C) patterns of leaf chlorophyll (Chl) across grasslands in the Tibetan Plateau. Red dots represent species-specific observations in each site. The blue lines and shades represent the regression lines with a 95% confidence band.
Figure 2. Relationships of leaf chlorophyll (Chl) with PAR [(A) photosynthetically active radiation], humidity index (B), AN [(C) soil available nitrogen], and AP [(D) soil available phosphorus]. Red dots represent species-specific observations in each site. The blue lines and shades represent the regression lines with a 95% confidence band.
Among different PFGs, we found a greater leaf Chl in C4 species (37.61 ± 0.91) than in C3 species (30.04 ± 0.85) (Figure 3). Leaf Chl in legume species (22.22 ± 0.59) was lower than in non-legume species (33.03 ± 0.79). However, leaf Chl showed no significant difference between annual and perennial species.
Figure 3. Leaf chlorophyll (Chl) in different functional groups. Different colors represent different comparisons of functional groups, such as C3 vs. C4 species, annual vs. perennial species, or legume vs. non-legume species. Bar height indicates the mean of leaf Chl with its standard error. Different letters represent a significant difference among different functional groups. The numbers in the bar represent the observations of leaf Chl in different functional groups.
The Dominant Role of Plant Evolution in Affecting Leaf Chl
Blomberg's K value indicated that the phylogenetic signal was not significant for leaf Chl among 93 species. Additionally, Pagel'λ values showed that variations in leaf Chl were mainly due to external environmental changes (not genetic variations) during plant evolutionary history (Supplementary Figure 3). Furthermore, we found that plant evolutionary divergence time was a good predictor for leaf Chl variations. A decreasingly evolutionary trend was detected for leaf Chl in a non-linear manner at both family (Figure 4A, R2 = 0.421, p < 0.001) and genus levels (Figure 4B, R2 = 0.163, p < 0.01). That is, modern species (lately evolved species) had a lower leaf Chl than ancestral species (early evolved species).
Figure 4. The relationships between leaf chlorophyll (Chl) and plant evolutionary time (million years) at both family (A) and genus levels (B). Different colors represent different plant families or genera. The size of dots represents the number of site species at family or genus levels. The blue lines represent the regression lines.
Based on the SEM analyses, we further quantified that plant evolution was the dominant driver for leaf Chl variations across the Tibetan Plateau (Figure 5). Specifically, the path coefficient of the plant evolution effect was 0.347. However, the sum of coefficients of other influencing pathways was lower than that of the plant evolution effect. The coefficient of the HI direct effect was 0.159. Those of the indirect effects of HI were 0.090 (HI→PAR → Chl), 0.071 (HI→soil N → Chl), and −0.011 (HI→soil P→Chl). The individual coefficients of PAR, soil available N and P were −0.106, 0.090, and −0.093, respectively.
Figure 5. The direct and indirect impact pathways of plant evolutionary time, climate [humidity index, photosynthetically active radiation (PAR)], and soil nutrient [soil available nitrogen (N) and phosphorus (P)] in determining leaf Chl variations. The solid lines indicate significant pathways, while dashed ones represent marginally significant pathways (0.05 < p < 0.1). Blue lines are positive pathways, while red ones mean negative pathways. The numbers near the arrows represent standardized path coefficients, and the arrow width is proportional to the size of path coefficients. Significance levels are indicated by **p < 0.01; ***p < 0.001; ∧p < 0.1.
Discussion
As opposed to most previous studies conducted at the local scale (Croft et al., 2015, 2017; Lu et al., 2020), this study presented a large-scale pattern of leaf Chl among 677 site-species across grasslands on the Tibetan Plateau. Surprisingly, we found that plant evolutionary history was more dominant than current environmental changes in impacting the regional patterns of leaf Chl. This is different from the traditional view that leaf Chl is mainly affected by environmental factors, such as PAR (Genesio et al., 2020), soil water, and nutrient availability (Li et al., 2018b; Zhang et al., 2020). Moreover, we further detected the direction of evolutionary effects on leaf Chl, with the decreasing trend of leaf Chl along with evolution across the alpine grasslands. This Chl evolution response well-matches the increase in PAR during the uplift of the Tibetan Plateau. These indicate that the uplift of the Tibetan Plateau provides an ideal opportunity to study the evolutionary effects on plant traits. As leaf Chl is a fundamental component for photosynthesis, these results offer new insights into reconsidering plant photosynthesis capacity and the carbon cycle in an evolutionary view.
The Spatial Patterns and Environmental Control of Leaf Chl
Across grasslands on the Tibetan Plateau, the spatial variations of leaf Chl concentration are mainly impacted by climate (PAR, HI), soil nutrients, and PFG. Along climate gradient, we revealed a decrease in leaf Chl with increasing PAR, which is consistent with our expectation and supported by previous studies (Niinemets, 2010; Genesio et al., 2020; Zhang et al., 2020). It is likely because plants are usually exposed to strong radiation across the Tibetan Plateau (Valladares and Niinemets, 2008). To avoid excessive radiation damage, alpine plants usually reduce leaf Chl and chloroplast number to offset the harmful effect of high radiation (photo-protection response), as well as to protect chloroplast interior structure (Klem et al., 2012; Genesio et al., 2020). Furthermore, alpine plants usually invest extra resources to protect them from suffering high radiation, such as cell membrane (thylakoid and mesophyll) thickening (Hallik et al., 2009). This perhaps leads to the reduction in resources for Chl biosynthesis and leaf Chl as well. Along the water gradient, leaf Chl was suppressed in drier environments (decreasing HI). This is caused by the fact that long-term drought stress could damage the thylakoid membrane and mesophyll cells, further reducing leaf Chl concentration (Alberte and Thornber, 1977). Besides, plants tend to invest more resources for survival under drought stress (Alberte and Thornber, 1977; Mihailović et al., 1997; Lei et al., 2006), indirectly decreasing resource allocation to leaf Chl (Bokhari, 1976; Estill et al., 1991; Mihailović et al., 1997).
Soil nutrient availability was also found to determine leaf Chl variations. The increase in leaf Chl with soil available N revealed in this study is in line with previous studies. For example, a positive relationship is normally detected between leaf N and Chl concentration (Croft et al., 2017; Zhang et al., 2020), and leaf N concentration is found to increase with soil available N (Zhao et al., 2003, 2005; Boussadia et al., 2011). This results from the fundamental role of the N element in constituting leaf Chl (Croft et al., 2017). In contrast, we found little effect of soil P availability on leaf Chl, likely caused by stronger N limitation than P limitation in alpine grasslands (Gao et al., 2018; Du et al., 2020).
A significant difference was detected among PFGs for leaf Chl as well. The greater leaf Chl in C4 than C3 species is partly due to the fact that C4 species have a stronger capacity to harvest CO2 than C3 species (Wang et al., 2011; Christin and Osborne, 2014). In contrast with our expectations and previous findings (Zhang et al., 2020), we revealed no difference in leaf Chl between annual and perennial species. Annual species are usually expected to invest more resources in fast aboveground growth and reproduction (high leaf Chl with resource-use strategy), but perennial species allocate more to belowground growth for the advantages in long-term persistence (low Chl with resource-conservative strategy) (Wright et al., 2004). However, across our grassland transect in the Tibetan Plateau, plants face multiple stresses of extremely low temperature, high radiation, and drought (Chen et al., 2013). Thus, the greatest pressure faced by alpine plants might be survival, not fast growth for resource competition. Furthermore, annual plants tend to invest more resources in reproduction. These together cause little difference in leaf Chl between annual and perennial species. Contrary to our expectation, we found a lower leaf Chl of legume than non-legume species. Generally, legumes have the advantage of obtaining soil N via rhizobium, possibly promoting leaf N and Chl (Zahran, 1999; Graham and Vance, 2003; Wei et al., 2019). Nevertheless, our results indicate that more N uptake in legumes in alpine may not be mainly used for leaf Chl biosynthesis. The symbiont of rhizobium costs a lot of resources and energy (Zahran, 1999). Besides, more N needs to be invested in resisting the stress of coldness, high radiation, and drought (Wang et al., 2018). These indicate that legume tends to increase investment in belowground resources and abiotic stress resistance, which may reduce aboveground allocation and further decrease leaf Chl. In addition, we found that leaf Chl in the grass species was lower than in the forb species (Supplementary Figure 5), which is supported by recent studies in alpine and temperate grasslands (Zhang et al., 2020).
The Role of Plant Evolution in Affecting Leaf Chl
Across grasslands in the Tibetan Plateau, we revealed a decreasing trend of leaf Chl evolution among 677 site-species (93 species from 25 families). Our result is among the first to detect the direction of leaf Chl evolution, complementing previous studies that found the evolutionary divergence of leaf Chl among grassland species (Li et al., 2018a). This evolution pattern well-matches the gradual increase in PAR during the geological uplift of the Tibetan Plateau. Specifically, about 50 million years ago, the India-Eurasia continental collision rapidly uplifted the Tibetan Plateau until now (Royden et al., 2008). This geological uplift causes a gradual enhancement in PAR during plant evolutionary history (Blumthaler et al., 1997). Facing high radiation, plants do not necessarily invest more resources in leaf Chl for light capture (Qiu et al., 2011; Genesio et al., 2020; Zhang et al., 2020). To avoid strong radiation damage, alpine plants tend to evolve toward the strategy with lower leaf Chl. The additional resource from decreasing leaf Chl might be allocated to resist high radiation, such as densely villous aggregated globose, thicken bracts, and other morphological adaptions (Wen et al., 2014; Ibanez et al., 2017). Besides increased PAR, the temperature gradually reduces with the uplift (Zhisheng et al., 2001; Liu et al., 2006; Ge et al., 2014). On one hand, low temperature suppresses plant physiological activities (Bokhari, 1976; Fu et al., 2014; Li et al., 2018b), further reducing the biosynthesis of leaf Chl (Bokhari, 1976). On the other hand, low temperature forces alpine plants to invest more resources in cold resistance (Li et al., 2018b), likely cutting down allocation to leaf Chl during plant evolutionary history.
Moreover, we found that leaf Chl non-linearly declined with plant evolutionary divergence time, with a faster decrease recently. This non-linear evolution pattern also matches well with the different uplift stages. Before the continental collision (c. 50 million years ago), the crust movement is relatively inactive (Royden et al., 2008), and the evolutionary rate of leaf Chl was slow as revealed in our results. The Tibetan Plateau began to uplift after the collision (Rowley and Currie, 2006), while the uplift has been again accelerated about 20 million years ago (Royden et al., 2008; Wang et al., 2008). The different rates of the Tibetan Plateau uplift among multiple stages again support our results that the decreasing rates of leaf Chl evolution were quicker. These suggest that the uplift of the Tibetan Plateau provides an excellent opportunity for future studies to test the evolutionary effects on plant traits and functions. The greatest decrease in leaf Chl evolution recently calls our attention to the evolutionary constraint on plant photosynthesis and carbon cycle, especially for the regions where climate is quickly becoming warmer and wetter over the past decades (Chen et al., 2013).
Interestingly, the evolution of leaf Chl showed great variations in the lately evolved stage. This implies that some species evolve quickly toward the strategy with lower leaf Chl (Kita et al., 2004; Wang et al., 2004; Liu et al., 2013), but other species tend to inherit from ancestral traits when facing great environmental changes during evolutionary history (Wen et al., 2014; Volis et al., 2018). Furthermore, we found that the quickly evolved species were mainly coming from Poaceae, while the conservatively evolved species were from Asteraceae and C4 functional group. These suggest that it is interesting and valuable to study the mechanisms underlying different Chl evolutionary responses among different species with the uplift of the Tibetan Plateau, as well as the consequence for plant photosynthesis and the carbon cycle in alpine grasslands.
Overall, the SEM analysis concluded that plant evolution was the dominant driver for leaf Chl variations when compared with the individual or sum effects of current environmental changes (i.e., HI, PAR, and soil nutrient availability). Contrary to the traditional view that leaf Chl is mainly affected by climate, soil nutrients, and PFG (Zhang et al., 2020), our result highlights the dominant role of plant evolution in driving leaf Chl of alpine plants across the Tibetan Plateau. Interestingly, future efforts are deserved to quantify the other sources of great variations in leaf Chl across contrasting environments, such as soil Mg and Fe availability. As leaf Chl is fundamental to harvesting light for plant photosynthesis, these results have important implications as follows. (1) The geological uplift of the Tibetan Plateau offers an ideal opportunity to explore the impacts of plant evolution on leaf photosynthetic traits and the underlying mechanisms. (2) These results suggest the potential role of plant evolutionary history in photosynthetic capacity, which is largely ignored in previous studies. (3) Plant evolution might function on the ecosystem-level carbon cycle, arousing our awareness to reconsider the carbon budget across the Tibetan Plateau in an evolutionary view.
Conclusion
As opposed to previously local-scale studies, our study presented a large-scale pattern of leaf Chl among 677 site-species across grasslands in the Tibetan Plateau. More importantly, we found that it was plant evolution, rather than current environment changes, dominantly shaping leaf Chl variations in alpine plants across contrasting environments. This challenges the traditional view that leaf Chl is mainly regulated by climate, soil resource, and PFG. Interestingly, we also detected the direction of evolutionary effects on leaf Chl, with a non-linear decreasing trend during plant evolution history. This evolutionary response well matches the non-linear enhancement in PAR during different uplift stages of the Tibetan Plateau. These suggest that the geological uplift of the Tibetan Plateau provides a great chance to study the impact of plant evolution on alpine plant traits and functions. These new findings arouse our attention to reconsidering ecosystem carbon balance across the Tibetan Plateau from an evolutionary perspective.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
YH and DT designed the experiments. YH, RZ, and YL performed the field investigation and collected the data. YH, TL, and DT analyzed the data and wrote the paper. TL, DT, XC, and JL offered helpful comments. All authors contributed to the draft and revision of the manuscript.
Funding
This study was financially supported by the Second Tibetan Plateau Scientific Expedition and Research (STEP) Program (Grant No. 2019QZKK0302), the National Key R&D Program of China (2017YFA0604802), the Kezhen-Bingwei Young Talents (2020RC003), and the Youth Innovation Promotion Association CAS (2021050).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.941983/full#supplementary-material
References
Alberte, R. S., and Thornber, J. P. (1977). Water stress effects on the content and organization of chlorophyll in mesophyll and bundle sheath chloroplasts of maize. Plant Physiol. 59, 351–353. doi: 10.1104/pp.59.3.351
Anjum, S. A., Farooq, M., Wang, L. C., Xue, L. L., Wang, S. G., Wang, L., et al. (2011). Gas exchange and chlorophyll synthesis of maize cultivars are enhanced by exogenously-applied glycinebetaine under drought conditions. Plant Soil Environ. 57, 326–331. doi: 10.17221/41/2011-PSE
Blomberg, S. P., Garland, T Jr., and Ives, A. R. (2003). Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 57, 717–745. doi: 10.1111/j.0014-3820.2003.tb00285.x
Blumthaler, M., Ambach, W., and Ellinger, R. (1997). Increase in solar uv radiation with altitude. J. Photochem. Photobiol. B. Biol. 39, 130–134. doi: 10.1016/S1011-1344(96)00018-8
Bokhari, U. G. (1976). Influence of temperatures, water stress, and nitrogen treatments on chlorophyll and dry matter of western wheatgrass. J. Range Manag. 29, 127. doi: 10.2307/3897408
Boussadia, O., Steppe, K., Zgallai, H., El Hadj, S. B., Braham, M., Lemeur, R., et al. (2011). Nondestructive determination of nitrogen and chlorophyll content in olive tree leaves and the relation with photosynthesis and fluorescence parameters. Photosynthetica. 49, 149–153. doi: 10.1007/s11099-011-0019-x
Cai, H., Chu, Q., Yuan, L., Liu, J., Chen, X., Chen, F., et al. (2011). Identification of quantitative trait loci for leaf area and chlorophyll content in maize (zea mays) under low nitrogen and low phosphorus supply. Mol. Breed. 30, 251–266. doi: 10.1007/s11032-011-9615-5
Carswell, F. E., Meir, P., Wandelli, E. V., Bonates, L. C., Kruijt, B., Barbosa, E. M., et al. (2000). Photosynthetic capacity in a central amazonian rain forest. Tree Physiol. 20, 179–186. doi: 10.1093/treephys/20.3.179
Chen, H., Zhu, Q., Peng, C., Wu, N., Wang, Y., Fang, X., et al. (2013). The impacts of climate change and human activities on biogeochemical cycles on the qinghai-tibetan plateau. Glob Chang Biol. 19, 2940–2955. doi: 10.1111/gcb.12277
Christin, P. A., and Osborne, C. P. (2014). The evolutionary ecology of c4 plants. New Phytol. 204, 765–781. doi: 10.1111/nph.13033
Croft, H., Chen, J. M., Froelich, N. J., Chen, B., and Staebler, R. M. (2015). Seasonal controls of canopy chlorophyll content on forest carbon uptake: Implications for gpp modeling. J. Geophys. Res. Biogeosci. 120, 1576–1586. doi: 10.1002/2015JG002980
Croft, H., Chen, J. M., Luo, X., Bartlett, P., Chen, B., and Staebler, R. M. (2017). Leaf chlorophyll content as a proxy for leaf photosynthetic capacity. Glob. Chang. Biol. 23, 3513–3524. doi: 10.1111/gcb.13599
Du, E., Terrer, C., Pellegrini, A. F. A., Ahlström, A., van Lissa, C. J., Zhao, X., et al. (2020). Global patterns of terrestrial nitrogen and phosphorus limitation. Nat. Geosci. 13, 221–226. doi: 10.1038/s41561-019-0530-4
Estill, K., Delaney, R. H., Smith, W. K., and Ditterline, R. L. (1991). Water relations and productivity of alfalfa leaf chlorophyll variants. Crop Sci. 31, 1229–1233. doi: 10.2135/cropsci1991.0011183X003100050030x
Flynt, A., and Dean, N. (2016). A survey of popular r packages for cluster analysis. J. Educ. Behav. Stat. 41, 205–225. doi: 10.3102/1076998616631743
Fu, G., and Shen, Z.-X. (2016). Response of alpine plants to nitrogen addition on the tibetan plateau: a meta-analysis. J. Plant Growth Regul. 35, 974–979. doi: 10.1007/s00344-016-9595-0
Fu, G., Shen, Z.-X., Sun, W., Zhong, Z.-M., Zhang, X.-Z., and Zhou, Y.-T. (2014). A meta-analysis of the effects of experimental warming on plant physiology and growth on the tibetan plateau. J. Plant Growth Regul. 34, 57–65. doi: 10.1007/s00344-014-9442-0
Gao, Y., Cooper, D. J., and Zeng, X. (2018). Nitrogen, not phosphorus, enrichment controls biomass production in alpine wetlands on the tibetan plateau, China. Ecol. Eng. 116, 31–34. doi: 10.1016/j.ecoleng.2018.02.016
Ge, X., Liu, J., Ren, S., and Yuan, S. (2014). Tectonic uplift of tibetan plateau: impacts on the formation of landforms, climate changes and ancient human migration in china. Geol. China. 41, 698–714.
Genesio, L., Bassi, R., and Miglietta, F. (2020). Plants with less chlorophyll: a global change perspective. Glob. Chang. Biol. 27, 959–967. doi: 10.1111/gcb.15470
Ginestet, C. (2011). Ggplot2: Elegant graphics for data analysis. J. R. Stat. Soc. Ser. A. 174, 245–246. doi: 10.1111/j.1467-985X.2010.00676_9.x
Graham, P. H., and Vance, C. P. (2003). Legumes: importance and constraints to greater use. Plant Physiol. 131, 872–877. doi: 10.1104/pp.017004
Guo, X., Wang, R., Chang, R., Liang, X., Wang, C., Luo, Y., et al. (2014). Effects of nitrogen addition on growth and photosynthetic characteristics of acer truncatum seedlings. Dendrobiology. 72, 151–161. doi: 10.12657/denbio.072.013
Guo, Y., Yang, X., Schob, C., Jiang, Y., and Tang, Z. (2017). Legume shrubs are more nitrogen-homeostatic than non-legume shrubs. Front. Plant Sci. 8, 1662. doi: 10.3389/fpls.2017.01662
Hallik, L., Kull, O., Niinemets, Ü., and Aan, A. (2009). Contrasting correlation networks between leaf structure, nitrogen and chlorophyll in herbaceous and woody canopies. Basic Appl. Ecol. 10, 309–318. doi: 10.1016/j.baae.2008.08.001
Ibanez, V. N., Berli, F. J., Masuelli, R. W., Bottini, R. A., and Marfil, C. F. (2017). Influence of altitude and enhanced ultraviolet-b radiation on tuber production, seed viability, leaf pigments and morphology in the wild potato species solanum kurtzianum bitter and wittm collected from an elevational gradient. Plant Sci. 261, 60–68. doi: 10.1016/j.plantsci.2017.04.014
Kattge, J., Knorr, W., Raddatz, T., and Wirth, C. (2009). Quantifying photosynthetic capacity and its relationship to leaf nitrogen content for global-scale terrestrial biosphere models. Glob. Chang. Biol. 15, 976–991. doi: 10.1111/j.1365-2486.2008.01744.x
Kelly, R., Healy, K., Anand, M., Baudraz, M. E. A., Bahn, M., Cerabolini, B. E. L., et al. (2021). Climatic and evolutionary contexts are required to infer plant life history strategies from functional traits at a global scale. Ecol Lett. 24, 970–983. doi: 10.1111/ele.13704
Kita, Y., Fujikawa, K., Ito, M., Ohba, H., and Kato, M. (2004). Molecular phylogenetic analyses and systematics of the genus saussurea and related genera (asteraceae, cardueae). Taxon. 53, 679–690. doi: 10.2307/4135443
Klem, K., A.č, A., Holub, P., Kov,áč, D., Špunda, V, Robson, T. M., et al. (2012). Interactive effects of par and uv radiation on the physiology, morphology and leaf optical properties of two barley varieties. Environ. Exp. Bot. 75, 52–64. doi: 10.1016/j.envexpbot.2011.08.008
Knapp, A. K., and Smith, M. D. (2001). Variation among biomes in temporal dynamics of aboveground primary production. Science. 291, 481–484. doi: 10.1126/science.291.5503.481
Lambers, H., Finnegan, P. M., Laliberte, E., Pearse, S. J., Ryan, M. H., Shane, M. W., et al. (2011). Update on phosphorus nutrition in proteaceae. Phosphorus nutrition of proteaceae in severely phosphorus-impoverished soils: Are there lessons to be learned for future crops? Plant Physiol. 156, 1058–1066. doi: 10.1104/pp.111.174318
Lei, Y., Yin, C., and Li, C. (2006). Differences in some morphological, physiological, and biochemical responses to drought stress in two contrasting populations of populus przewalskii. Physiologia Plantarum. 127, 182–191. doi: 10.1111/j.1399-3054.2006.00638.x
Li, Y., He, N., Hou, J., Xu, L., Liu, C., Zhang, J., et al. (2018a). Factors influencing leaf chlorophyll content in natural forests at the biome scale. Front. Ecol. Evol. 6, 10. doi: 10.3389/fevo.2018.00064
Li, Y., Liu, C., Zhang, J., Yang, H., Xu, L., Wang, Q., et al. (2018b). Variation in leaf chlorophyll concentration from tropical to cold-temperate forests: association with gross primary productivity. Ecol. Indic. 85, 383–389. doi: 10.1016/j.ecolind.2017.10.025
Liu, B., Opgenoorth, L., Miehe, G., Zhang, D., Wan, D., Zhao, C., et al. (2013). Molecular bases for parallel evolution of translucent bracts in an alpine “glasshouse” plantrheum alexandrae(polygonaceae). J. Syst. Evol. 51, 134–141. doi: 10.1111/j.1759-6831.2012.00225.x
Liu, J. Q., Wang, Y. J., Wang, A. L., Hideaki, O., and Abbott, R. J. (2006). Radiation and diversification within the ligularia-cremanthodium-parasenecio complex (asteraceae) triggered by uplift of the qinghai-tibetan plateau. Mol. Phylogenet. Evol. 38, 31–49. doi: 10.1016/j.ympev.2005.09.010
Lu, X., Ju, W., Li, J., Croft, H., Chen, J. M., Luo, Y., et al. (2020). Maximum carboxylation rate estimation with chlorophyll content as a proxy of rubisco content. J. Geophys. Res. Biogeosci. 125:e2020JG005748. doi: 10.1029/2020JG005748
Luo, X., Croft, H., Chen, J. M., He, L., and Keenan, T. F. (2019). Improved estimates of global terrestrial photosynthesis using information on leaf chlorophyll content. Glob. Chang. Biol. 25, 2499–2514. doi: 10.1111/gcb.14624
Martonne, D. (1926). Une nouvelle fonction climatologique: L'indice d'aridité. La Meteorologie. 2, 449–458. doi: 10.3406/geo.1926.8506
Masuda, T., Takabe, K., Ohta, H., Shioi, Y., and Takamiya, K. (1996). Enzymatic activities for the synthesis of chlorophyll in pigment-deficient variegated leaves of euonymus japonicus. Plant Cell Physiol. 37, 481–487. doi: 10.1093/oxfordjournals.pcp.a028970
Mihailović, N., Lazarević, M., DŽeletović, Z., Vučković, M., and 90Aurdević, M. (1997). Chlorophyllase activity in wheat, triticum aestivum l. Leaves during drought and its dependence on the nitrogen ion form applied. Plant Sci. 129, 141–146. doi: 10.1016/S0168-9452(97)00189-1
Monroe, J. G., Gill, B., Turner, K. G., and McKay, J. K. (2019). Drought regimens predict life history strategies in heliophila. New Phytol. 223, 2054–2062. doi: 10.1111/nph.15919
Niinemets, Ü. (2010). A review of light interception in plant stands from leaf to canopy in different plant functional types and in species with varying shade tolerance. Ecol. Res. 25, 693–714. doi: 10.1007/s11284-010-0712-4
Olsen, S. R. (1954). Estimation of Available Phosphorus in Soils by Extraction With Sodium Bicarbonate. Estimation of available phosphorus in soils by extraction with sodium bicarbonate. Washington, DC: United States Department of Agriculture.
Onoda, Y., Wright, I. J., Evans, J. R., Hikosaka, K., Kitajima, K., Niinemets, U., et al. (2017). Physiological and structural tradeoffs underlying the leaf economics spectrum. New Phytol. 214, 1447–1463. doi: 10.1111/nph.14496
Osborne, C. P., Salomaa, A., Kluyver, T. A., Visser, V., Kellogg, E. A., Morrone, O., et al. (2014). A global database of c4 photosynthesis in grasses. New Phytol. 204, 441–446. doi: 10.1111/nph.12942
Pan, J., Wang, J., Zhang, R., Tian, D., Cheng, X., Wang, S., et al. (2021). Microaggregates regulated by edaphic properties determine the soil carbon stock in tibetan alpine grasslands. Catena. 206, 105570. doi: 10.1016/j.catena.2021.105570
Qian, H., and Jin, Y. (2016). An updated megaphylogeny of plants, a tool for generating plant phylogenies and an analysis of phylogenetic community structure. J. Plant Ecol. 9, 233–239. doi: 10.1093/jpe/rtv047
Qiu, Y. X., Fu, C. X., and Comes, H. P. (2011). Plant molecular phylogeography in china and adjacent regions: Tracing the genetic imprints of quaternary climate and environmental change in the world's most diverse temperate flora. Mol. Phylogenet. Evol. 59, 225–244. doi: 10.1016/j.ympev.2011.01.012
Rastogi, A., Antala, M., Gabka, M., Rosadzinski, S., Strozecki, M., Brestic, M., et al. (2020). Impact of warming and reduced precipitation on morphology and chlorophyll concentration in peat mosses (sphagnum angustifolium and s. Fallax). Sci. Rep. 10, 8592. doi: 10.1038/s41598-020-65032-x
Ren, X., He, H., Zhang, L., and Yu, G. (2018). Global radiation, photosynthetically active radiation, and the diffuse component dataset of china, 1981–2010. Earth Syst. Sci. Data. 10, 1217–1226. doi: 10.5194/essd-10-1217-2018
Rhein, A., and Silva, M. A. (2016). Nitrogen doses on physiological attributes and yield of sugarcane grown under subsurface drip fertigation. J. Plant. Nutr. 40, 227–238. doi: 10.1080/01904167.2016.1237646
Rosseel, Y. (2012). Lavaan: an r package for structural equation modeling. J. Stat. Softw. 48, 1–36. doi: 10.18637/jss.v048.i02
Rowley, D. B., and Currie, B. S. (2006). Palaeo-altimetry of the late eocene to miocene lunpola basin, central tibet. Nature. 439, 677–681. doi: 10.1038/nature04506
Royden, L. H., Burchfiel, B. C., and van der Hilst, R. D. (2008). The geological evolution of the tibetan plateau. Science. 321, 1054–1058. doi: 10.1126/science.1155371
Ryu, Y., Berry, J. A., and Baldocchi, D. D. (2019). What is global photosynthesis? History, uncertainties and opportunities. Remote Sens. Environ. 223, 95–114. doi: 10.1016/j.rse.2019.01.016
Schmerler, S. B., Clement, W. L., Beaulieu, J. M., Chatelet, D. S., Sack, L., Donoghue, M. J., et al. (2012). Evolution of leaf form correlates with tropical-temperate transitions in viburnum (adoxaceae). Proc. Biol. Sci. 279, 3905–3913. doi: 10.1098/rspb.2012.1110
Uddling, J., Gelang-Alfredsson, J., Piikki, K., and Pleijel, H. (2007). Evaluating the relationship between leaf chlorophyll concentration and spad-502 chlorophyll meter readings. Photosynth. Res. 91, 37–46. doi: 10.1007/s11120-006-9077-5
Valladares, F., and Niinemets, Ü. (2008). Shade tolerance, a key plant feature of complex nature and consequences. Annu. Rev. Ecol. Evol. Syst. 39, 237–257. doi: 10.1146/annurev.ecolsys.39.110707.173506
Volis, S., Fogel, K., Tu, T., Sun, H., and Zaretsky, M. (2018). Evolutionary history and biogeography of mandragora l. (solanaceae). Mol. Phylogenet. Evol. 129, 85–95. doi: 10.1016/j.ympev.2018.08.015
Wang, C., Zhao, X., Liu, Z., Lippert, P. C., Graham, S. A., Coe, R. S., et al. (2008). Constraints on the early uplift history of the tibetan plateau. Proc. Natl. Acad. Sci. U. S. A. 105, 4987–4992. doi: 10.1073/pnas.0703595105
Wang, D., Ling, T., Wang, P., Jing, P., Fan, J., Wang, H., et al. (2018). Effects of 8-year nitrogen and phosphorus treatments on the ecophysiological traits of two key species on tibetan plateau. Front. Plant Sci. 9, 1290. doi: 10.3389/fpls.2018.01290
Wang, L., Peterson, R. B., and Brutnell, T. P. (2011). Regulatory mechanisms underlying c4 photosynthesis. New Phytol. 190, 9–20. doi: 10.1111/j.1469-8137.2011.03649.x
Wang, Y. J., Li, X. J., Hao, G., and Liu, J. Q. (2004). Molecular phylogeny and biogeography of androsace (primulaceae) and the convergent evolution of cushion morphology. Acta Phytotaxonomica Sinica. 42, 481–499.
Wei, X., Reich, P. B., and Hobbie, S. E. (2019). Legumes regulate grassland soil n cycling and its response to variation in species diversity and n supply but not co2. Glob. Chang. Biol. 25, 2396–2409. doi: 10.1111/gcb.14636
Wen, J., Zhang, J. Q., Nie, Z. L., Zhong, Y., and Sun, H. (2014). Evolutionary diversifications of plants on the qinghai-tibetan plateau. Front Genet. 5, 4. doi: 10.3389/fgene.2014.00004
Wright, I. J., Reich, P. B., Westoby, M., Ackerly, D. D., Baruch, Z., Bongers, F., et al. (2004). The worldwide leaf economics spectrum. Nature. 428, 821–827. doi: 10.1038/nature02403
Xiong, D., and Nadal, M. (2020). Linking water relations and hydraulics with photosynthesis. Plant J. 101, 800–815. doi: 10.1111/tpj.14595
Xu, Z., Wan, C., Xiong, P., Tang, Z., Hu, R., Cao, G., et al. (2010). Initial responses of soil co2 efflux and c, n pools to experimental warming in two contrasting forest ecosystems, eastern tibetan plateau, china. Plant Soil. 336, 183–195. doi: 10.1007/s11104-010-0461-8
Yuan, Z., Cao, Q., Zhang, K., Ata-Ul-Karim, S. T., Tian, Y., Zhu, Y., et al. (2016). Optimal leaf positions for SPAD meter measurement in rice. Front. Plant. Sci. 7, 719. doi: 10.3389/fpls.2016.00719
Zahran, H. H. (1999). Rhizobium-legume symbiosis and nitrogen fixation under severe conditions and in an arid climate. Microbiol. Mol. Biol. Rev. 63, 968–989, table of contents. doi: 10.1128/MMBR.63.4.968-989.1999
Zhang, Y., Li, Y., Wang, R., Xu, L., Li, M., Liu, Z., et al. (2020). Spatial variation of leaf chlorophyll in northern hemisphere grasslands. Front. Plant Sci. 11, 1244. doi: 10.3389/fpls.2020.01244
Zhao, D., Raja Reddy, K., Kakani, V. G., Read, J. J., and Carter, G. A. (2003). Corn (zea mays l.) growth, leaf pigment concentration, photosynthesis and leaf hyperspectral reflectance properties as affected by nitrogen supply. Plant Soil. 257, 205–218. doi: 10.1023/A:1026233732507
Zhao, D., Reddy, K. R., Kakani, V. G., and Reddy, V. R. (2005). Nitrogen deficiency effects on plant growth, leaf photosynthesis, and hyperspectral reflectance properties of sorghum. Eur. J. Agron. 22, 391–403. doi: 10.1016/j.eja.2004.06.005
Keywords: alpine grassland, evolutionary history, leaf chlorophyll, photosynthesis, photosynthetically active radiation, plant functional group
Citation: He Y, Li T, Zhang R, Wang J, Zhu J, Li Y, Chen X, Pan J, Shen Y, Wang F, Li J and Tian D (2022) Plant Evolution History Overwhelms Current Environment Gradients in Affecting Leaf Chlorophyll Across the Tibetan Plateau. Front. Plant Sci. 13:941983. doi: 10.3389/fpls.2022.941983
Received: 12 May 2022; Accepted: 13 June 2022;
Published: 11 July 2022.
Edited by:
Xiaoming Kang, Chinese Academy of Forestry, ChinaReviewed by:
Xingyun Liang, South China Botanical Garden (CAS), ChinaKailiang Yu, Princeton University, United States
Copyright © 2022 He, Li, Zhang, Wang, Zhu, Li, Chen, Pan, Shen, Wang, Li and Tian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingwen Li, lijingwenhy@bjfu.edu.cn; Dashuan Tian, tiands@igsnrr.ac.cn
†ORCID: Jinsong Wang orcid.org/0000-0002-3425-7387
Juntao Zhu orcid.org/0000-0003-0379-3213
Xinli Chen orcid.org/0000-0003-0542-5959
Ying Shen orcid.org/0000-0002-6874-3500
Dashuan Tian orcid.org/0000-0001-8023-1180
‡ These authors have contributed equally to this work