Yin Wen
Qing Ye
Cristian Román-Palacios
Hui Liu
Guilin Wu- 1Key Laboratory of Vegetation Restoration and Management of Degraded Ecosystems, Guangdong Provincial Key Laboratory of Applied Botany, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- 2College of Life Sciences, Gannan Normal University, Ganzhou, China
- 3School of Information, University of Arizona, Tucson, AZ, United States
Understanding how plants respond to thermal stress is central to predicting plant responses and community dynamics in natural ecosystems under projected scenarios of climate change. Although physiological tolerance is suggested to evolve slower than climatic niches, this comparison remains to be addressed in plants using a phylogenetic comparative approach. In this study, we compared i) the evolutionary rates of physiological tolerance to extreme temperatures with ii) the corresponding rates of climatic niche across three major vascular plant groups. We further accounted for the potential effects of hardening when examining the association between physiological and climatic niche rates. We found that physiological cold tolerance evolves faster than heat tolerance in all three groups. The coldest climatic-niche temperatures evolve faster than the warmest climatic-niche temperatures. Importantly, evolutionary rates of physiological cold tolerance were faster than rates of change in climatic niches. However, an inverse association between physiological cold tolerance and responding climatic niche for plants without hardening was detected. Our results indicated that plants may be sensitive to changes in warmer temperatures due to the slower evolutionary rates of heat tolerance. This pattern has deep implications for the framework that is being used to estimate climate-related extinctions over the upcoming century.
1 Introduction
Climate change (average 0.015°C/year in 1959–2015) (Glikson, 2016) is occurring faster than it did in any other geological time in the past (the maximum of 0.0015°C/year in the Palaeocene–Eocene thermal maximum [PETM]) (Zeebe et al., 2009). In principle, species exposed to stressful conditions could potentially avoid climate-related extinctions over the upcoming decades by either 1) dispersing to areas with more suitable climatic conditions or 2) shifting their niche through plastic/adaptive responses. However, previous studies have suggested that dispersal rates (latitudinal 16.9 km/decade, elevational 11 m/decade) are much slower than the ongoing rates of climate change (Chen et al., 2011). Therefore, it is reasonable to suggest that many plant species will still need to shift their niches, instead of dispersing, in order to avoid experiencing climate-related extinctions (Corlett and Westcott, 2013). Thus, understanding the specific mechanisms that underlie niche shifts bears special importance in predicting plant species responses under projected climate change.
Plant thermal tolerance, generally understood from physiological tolerance measurements or approached from the description of climatic niches, can be used to understand spatial and temporal patterns of species’ survival and growth. Physiological tolerance, including plant photosynthetic cold/heat tolerance, is derived from the physiological response curves to temperature changes. These physiological descriptors are usually calculated using relatively short-term experimental conditions (Geange et al., 2021). Conversely, climatic niches are calculated from climate variables across species’ geographical ranges and tend to long-term adaptation to climatic conditions (Colwell and Rangel, 2009). Climatic niche reflects the complex tolerance to thermal and water stress, as well as to a single physiological process (Bush et al., 2018). Therefore, shifts in thermal tolerance can be achieved through either physiological or climatic niche shifts. For example, congeneric species distributed in contrasting climatic zones are expected to evolve divergent physiological tolerance, such as the temperature at 50% loss of photosynthesis rate under cold stress, to adapt to local climates by enhancing the stress resistance directly (Armstrong et al., 2020). In contrast, the niche expansion of angiosperms into freezing areas is associated with deciduous phenology, herbaceous growth form, and narrow vessel conduit (Zanne et al., 2014). These growth form or phenology changes avoid (e.g., deciduous phenology and short life cycle take advantage of seasons with optimum temperatures) or indirectly strengthen (e.g., conduit narrowing decreases the potential for freezing/thawing-induced embolism) the freezing tolerance of plants, which, in turn, expand their climatic niche.
To our knowledge, only a few studies have estimated evolutionary rates of thermal tolerance in plant lineages, especially based on physiological tolerance measurements (Araujo et al., 2013; Jezkova and Wiens, 2016; Lancaster and Humphreys, 2020; Liu et al., 2020b; Perez and Feeley, 2020). However, no study has explicitly compared evolutionary rates of physiological tolerance and corresponding rates of climatic niches in plants. For instance, a recent study focused on both plants and animals suggested that the evolutionary rate of the coldest temperatures in the climatic niche was higher than that of the hottest temperatures (Liu et al., 2020b). Congruently, another study using a large database of physiological tolerance for 1,028 plant species found that evolutionary rates of cold tolerance are much faster than heat tolerance, although paired comparison of the two tolerances was not performed in the study (Lancaster and Humphreys, 2020). Nevertheless, despite the obvious constraints that physiological tolerance is expected to impose on climatic niches, previous studies have not compared the evolutionary rates between physiological tolerance and climatic niche in plants. In fact, although studies on the association between rates of climatic niche and physiological tolerance have been conducted in animals (Qu and Wiens, 2020), the relationship between climatic niche and physiological tolerance remains largely unexplored in plants. Given that physiological tolerance contributes to species’ climatic niches and that the evolution of climatic niches should lag behind in the evolution of physiological tolerance, we suggest that evolutionary rates of physiological tolerance are faster than the responding climatic niche in plant species.
The association between rates of climatic and physiological evolution in plants is likely mediated by additional factors. For instance, hardening is a form of plasticity that physiological tolerance enhances after a pre-experience of stress. Exposure to stress (e.g., extreme temperature) is expected to have consequences on plant structure and function, which, in turn, would alter the physiological tolerance (Buchner et al., 2017; Janmohammadi et al., 2018; Armstrong et al., 2020). Therefore, hardening should increase or decrease the physiological tolerance value, which consequently is expected to affect the ancestral state reconstructions and evolutionary rate estimates of physiological tolerance. For instance, Lancaster and Humphreys (2020) found that the estimates of evolutionary rates are higher in hardened than non-hardened plants for both cold and heat physiological tolerance (Lancaster and Humphreys, 2020). Despite the potential hardening plasticity for physiological tolerance, the climatic niche of a species would keep stable during a relatively long time period of decades of years (Liu et al., 2020a). The uncertain change of rate in hardening status may obscure the rate comparison between climatic niche and physiological tolerance. Therefore, we suggest that the hardening condition should influence the association of the comparison between climatic niche and physiological tolerance.
We compiled a dataset of variables related to plant physiological tolerance to test 1) whether evolutionary rates of physiological tolerance are faster than the responding climatic niche in plant species and 2) if hardening affects the association between climatic niche and physiological tolerance. We calculated the climatic niche of each species based on their native geographic distribution. We performed phylogenetic comparative analyses to calculate the overall evolutionary rates of physiological tolerance and climatic niches in both heat and cold tolerance to test our hypothesis that 1) evolutionary rates of physiological tolerance are higher than the responding climatic niche in plants and 2) hardening status would influence the results of the comparison between climatic niche and physiological tolerance.
2 Materials and methods
2.1 Data collection
Physiological tolerance data of plants were compiled from different published articles and books (Table S1), but primarily from the GlobTherm database (Bennett et al., 2018) and Lancaster and Humphreys (2020). Note that Lancaster and Humphreys (2020) focused on the global variation of plant physiological tolerance and provided insights into the sources of variation by focusing on local environments, phylogeny, and biogeographic histories. However, our study expands on this original analytical framework by conducting paired comparisons between the evolutionary rate of physiological tolerance and the corresponding climatic niche.
Study sites of the physiological measurements covered warm-to-cold areas across the globe for both heat tolerance and cold tolerance (Figure 1). To maximize the taxonomic breadth of our dataset, multiple traits for cold and heat physiological tolerance were included in our dataset, such as freezing resistance (FR), freezing temperature (FT), temperature at 50% leak (LT50), and temperature at 100% leak (LT100) for cold tolerance and temperature at 50% leak (LT50), temperature at 100% leak (LT100), maximum temperature at photosynthetic and respiratory machinery can function (Tmax), and critical temperature (CT) for heat tolerance. Data for physiological cold tolerance were summarized in our dataset as Tmin and physiological heat tolerance as Tmax. Although different approaches reflect different mechanisms to estimate species’ physiological tolerance values, alternative measurements tend to be strongly correlated (O'sullivan et al., 2013; O'sullivan et al., 2016). The use of physiological tolerance values derived from different approaches is not expected to significantly affect the big picture of rate comparison (Table S2).
Figure 1 Global distribution of the physiological heat tolerance (Tmax, red circle, 1,271 records) and cold tolerance (Tmin, green circle, 882 records) records analyzed in this study. Each point indicates one site on the map.
Following Lancaster and Humphreys (2020), our dataset was divided by considering the hardening status of plants (a binary feature; plants with and without hardening). However, there are still some plants that lack the necessary data to determine whether they have undergone hardening. Plants that have been identified as “hardened” undergo pre-acclimation through temperature treatments conducted under laboratory or greenhouse conditions. Alternatively, they may have experienced seasonal acclimation, which is documented in the original studies and highlights the seasonal variations in physiological tolerance. For species where physiological tolerance was available from multiple datasets, the average value across measurements was used. In total, our dataset summarized 882 cold tolerance records and 1,271 heat tolerance records, in 1,145 angiosperms, 99 gymnosperms, 111 ferns, and 42 bryophytes.
For each species, we retrieved occurrence data from the Global Biodiversity Information Facility (GBIF; https://www.gbif.org). We analyzed only records within the species’ native distribution areas (sensu Plants of the World Online database; https://powo.science.kew.org/). This means that any records outside the species’ original distribution areas, referred to as “outliers”, were manually excluded from the dataset. For each occurrence, we used the extract function in the raster package (version 3.5-2) to extract climatic variables of each species from 30 seconds resolution climatic data in the WorldClim 2.1 database (Fick and Hijmans, 2017). We focused on the maximum temperature of the warmest month (MTWM; Bio5) and the minimum annual temperature of the coldest month (MTCM; Bio6) for each occurrence. These two variables, MTWM and MTCM are, in fact, expected to represent the warmest and coldest temperatures in certain coordinates, potentially informing both the higher and lower bounds of climatic niche across the species range. Therefore, we described cold and warm climatic niches using 1) the average value of the MTWM and MTCM, 2) the top 90% values of MTWM (upper MTWM) and the bottom 10% of MTCM (lower MTCM) values across species’ ranges, and 3) the local MTCM or MTWM of the experimental site. As the comparison results using the upper MTWM or lower MTCM or local climate (MTCM or MTWM) were similar to those using average values (Tables S3, S4), we only reported the results derived from average values.
2.2 Time-calibrated phylogeny reconstruction
We used the taxonstand R package to standardize the species names in the database (Cayuela et al., 2012). The phylogenetic tree analyzed in this study was generated by subsampling the GBOTB.extend megatree (Smith and Brown, 2018; Jin and Qian, 2019) using the drop.tip function implemented in the ape package (Revell, 2012). The GBOTB.extend tree is the updated version of the fossil-calibrated time-tree GBOTB (Smith and Brown, 2018; Jin and Qian, 2019). This phylogeny samples over 70,000 vascular plant taxa and is widely used in macroecology and macroevolution (Landis et al., 2018; Igea and Tanentzap, 2020; Carta et al., 2022). This tree provides access to a high-resolution phylogeny of thousands of vascular plant species. From the 1,397 species with trait-level data, 42 species were not matched with the sampling in the tree. Therefore, we analyzed a total of 1,355 species with trait-level and phylogenetic information.
2.3 Data analyses
We compared the overall evolutionary rates of physiological tolerance to extreme temperatures and corresponding climatic-niche temperatures. Note that Lancaster and Humphreys (2020) also calculated the evolutionary rates of physiological cold and heat tolerance. However, they did not make paired comparisons between physiological tolerance and climatic niche. Here, we compared the rates between i) physiological cold tolerance (Tmin) and minimum temperatures in the coldest season (MTCM) and ii) physiological heat tolerance (Tmax) and maximum temperatures in the hottest season (MTWM). Additionally, we compared rates between iii) Tmin and Tmax or iv) MTCM and MTWM. We estimated the overall evolutionary rate (σ2) based on the Brownian motion model under maximum likelihood (Adams, 2013). We performed likelihood-ratio tests (LRTs) to test for rate differences in each group comparison. LRT values provide support to identify whether the rates are statistically different. We note that given the limited taxonomic sampling, the estimated overall rates (σ2) may be biased, but we only used the values in paired comparison for each comparison group. The R code associated with these analyses is presented in the supplement to this article and follows the structure of two previous studies (Adams, 2013; Liu et al., 2020b).
We also estimated absolute evolutionary rates of physiological tolerance (Tmin and Tmax) and climatic niche (MTWM and MTCM). We estimated these rates for the i–iv sections presented above to obtain the rate comparison results between physiological tolerance and climatic niche. We calculated absolute evolutionary rates as the difference between the estimated ancestor’s state and the descendant species’ state divided by the absolute divergence time (see below for ancestral state reconstructions; Quintero and Wiens, 2013). The resulting estimate represents how much the trait changes per unit of evolutionary time for a certain species (Quintero and Wiens, 2013; Jezkova and Wiens, 2016). We used four alternative models of continuous trait evolution to reconstruct the ancestor’s state of physiological tolerance and corresponding climatic niche: Brownian motion (BM; lambda = 1) (Felsenstein, 1973), Ornstein–Uhlenbeck (OU; with one optimum value) (Butler and King, 2004), white noise (WN), and Lambda (LA; 0 < lambda < 1) (Prokoph et al., 2008). We fit the models using the fitContinuous function in the geiger R package (version 2.0.10) (Pennell et al., 2014). For each trait, we based ancestral states on the best-fitting model, which was selected for each trait as the one with the lowest Akaike information criterion (AIC). We transformed the branch lengths based on the best-fitting model, estimated the ancestor’s value, and then calculated the absolute evolutionary rates for each terminal tip. We performed paired t-tests to compare the absolute evolutionary rate (ln-transformed) at the species level in each group (i to iv). The results of the absolute rate comparison were similar to those of overall rates (Tables S5, S6). Finally, to test the relationship between physiological tolerance and climatic niche, we used phylogenetic generalized least squares regression (PGLS) in R package caper (version 1.0.1) (Orme, 2018). We performed all analyses, and we generated graphs in R 4.0.2 (R Core Team, 2020).
2.4 Methodological caveats
First, taxon sampling and physiological data were limited in each plant group. Because the limited sampling ratio would overestimate the branch length and thus influence the values of the calculated evolutionary rates, instead of utilizing the absolute rate value (which is with biological significance) for quantitative description, we used the overall rates and only used pairwise comparison in the same group for qualitative comparison. In this case, incomplete sampling contributed the same effect for the rate bias and did not distort the results of our comparison. Previous studies also suggested that incomplete sampling has no effect on rate comparison (Liu et al., 2020b; Qu and Wiens, 2020). Second, physiological tolerance was not discriminated against on the basis of the measuring methods as suggested by Lancaster and Humphreys (Lancaster and Humphreys, 2020; Lancaster and Humphreys, 2021). Different methods used for physiological tolerance measurements have different biological significance and values and consequently influence the rate estimation. However, because the results of the present study depended on the rate comparison, and not on rate values, the different physiological methods may have no effect on the overall tendency of the results (Table S2).
3 Results
The overall evolutionary rate (σ2) of physiological cold tolerance (Tmin) was significantly higher than the corresponding climatic niche (MTCM) in all three major plant groups (all p < 0.001; Figure 2). Plants with hardening showed similar results, except the ferns that did not show significantly higher evolutionary rates in Tmin relative to those in MTCM (Table 1). However, for plants without hardening, the evolutionary rates of Tmin were lower than MTCM in all three major groups (Table 1).
Figure 2 Comparison of estimated overall evolutionary rates (σ2) between (A) cold tolerance (Tmin and minimum temperature of coldest month (MTCM)) and (B) heat tolerance (Tmax and maximum temperature of warmest month (MTWM)). LRTs were conducted to obtain the p-value. N, number of species; LRT, likelihood-ratio test.
Table 1 Comparison of estimated evolutionary rates (σ2) between physiological tolerance [cold tolerance (Tmin) and (Tmax)] and corresponding climatic niche [minimum temperature of coldest month (MTCM) and maximum temperature of warmest month (MTWM)].
In terms of heat tolerance for gymnosperms and angiosperms, there were no significant differences between the rates of physiological tolerance (Tmax) and corresponding climatic niche (MTWM). In ferns, Tmax evolved slower than MTWM (Figure 2). Plants with hardening showed no significant difference in the evolutionary rates of Tmax and MTWM (Table 1).
We found cold tolerance, in both physiological and climatic niches, to evolve faster than corresponding heat tolerance in all three major groups (p < 0.001) (Figure 3; Table S7). Tmin was significantly correlated with MTCM (R2 = 0.20, p < 0.001), whereas the correlation of hardened plants was stronger (R2 = 0.58, p < 0.001) than that of non-hardened plants (R2 = 0.10, p < 0.001). For Tmax, a significant but weak correlation between Tmax and MTWM was found in all species (R2 = 0.04, p < 0.001) and hardened plants (R2 = 0.02, p < 0.001), but not in non-hardened plants (R2 = 0.01, p = 0.12) (Figure 4).
Figure 3 Comparison of estimated evolutionary rates (σ2) between physiological tolerance and responding climatic niche in both (A) cold tolerance [minimum temperature of coldest month (MTCM) versus physiological tolerance (Tmin)] and (B) heat tolerance [maximum temperature of warmest month (MTWM) versus physiological tolerance (Tmax)] among different plant groups. LRTs were conducted to obtain the p-value. N, number of species; LRT, likelihood-ratio test.
Figure 4 Relationships between physiological cold tolerance and corresponding climatic niche. (A–C) Physiological cold tolerance (Tmin) versus minimum temperature of coldest month (MTCM). (D–F) Physiological heat tolerance (Tmax) versus maximum temperature of warmest month (MTWM). Different hardening statuses were presented in different plots: all data (A, D), hardened (B, E), and not hardened (C, F). Regression lines, R2, and p-values are based on phylogenetic generalized least squares (PGLS) regression analyses.
4 Discussion
We leveraged a large species-level dataset on species climatic responses based on physiological and climatic niches, combined with phylogenetic information. We found that the evolutionary rates of physiological cold tolerance were higher than those of the corresponding climatic niche, while the cold tolerance of plants without hardening showed the opposite trend. However, the evolutionary rates of physiological heat tolerance were not different from those of the corresponding climatic niche, in both hardening and non-hardening plants. Our results suggest that physiological responses are likely the basis of ecological adaptation under cold stress. Consistent with previous studies (Lancaster and Humphreys, 2020; Liu et al., 2020b), we show that cold tolerance evolves faster than heat tolerance in both physiological tolerance and climatic niche. Thus, cold temperatures may be one of the key drivers of plant evolution. However, the limited evolutionary rate of heat tolerance may limit the future performance of plants in the coming warm world.
4.1 Physiological tolerances underlie climatic niches
Our results show that physiological cold tolerance evolved faster than the responding climatic niche for hardened plants, whereas the non-hardened plants showed the opposite trend. First, local microclimate, rather than macroclimate, could explain the plant responses in communities because it is directly the temperature that plant individuals experience and is more diverse compared to sole and more stable macroclimate (Zellweger et al., 2020). Physiological tolerance faces a higher selective pressure for the sensitivity of the physiological process to microclimate change, therefore showing higher evolutionary rates than macro climatic niches. Second, hardened plants in physiological cold tolerance and physiological heat tolerance were not consistent with the higher rates in climatic niches reported for animals (Qu and Wiens, 2020). This may be attributed to the ability of animals to move in order to avoid extreme temperatures (behavior avoidance) or to keep body temperature largely constant (homeothermy) (Hetem et al., 2012; Pereira et al., 2020). In contrast, plant migration is achieved only in next-generation individuals via seed dispersal, which is a slow process that lags behind climate warming (Corlett and Westcott, 2013). Therefore, physiological changes or acclimation are essential for the survival of plants.
Our results highlight the importance of hardening to physiological tolerance, especially to cold tolerance. Previous studies have addressed the physiological basis of the hardening effect on plant tolerance; for example, the changes in carbohydrate metabolism (Knaupp et al., 2011) and downregulation of photosynthesis (Frechette et al., 2015) enhanced cold tolerance in a relatively short time. This hardening process also accelerates the evolutionary rate of physiological tolerance at higher rates than the responding climatic niche. However, without hardening, plants with limited tolerance to both cold and heat contributed to lower or similar evolutionary rates than the corresponding climatic niche. In this case, the results were partly consistent with the observations in animals (Qu and Wiens, 2020).
It has been suggested that the physiological tolerance to extreme temperatures is correlated to the species’ climatic niche. However, our study revealed enhanced correlations in hardened plants and weakened relationships in non-hardened plants, suggesting that these correlations may partly be due to the physiological acclimation to natural climate-induced hardening. Thus, our results indicate that physiological tolerance, especially cold tolerance, constituted the basis of the climatic niche.
4.2 Why does cold tolerance evolve faster than heat tolerance?
Several potential explanations have been postulated regarding the higher rates of cold tolerance. First, the spatial distribution of the Earth’s surface temperature shapes both physiological tolerance and climatic niche. The lower variance in annual maximum temperature relative to minimum temperature across Earth biomes contributes to lower variance of heat tolerance and consequently lowers the rates of evolution for heat tolerance. Second, the historical thermal dynamics in geological time could have influenced plant evolution. For instance, at ~65 Ma, the surface temperature of Earth became cooler, despite a few rapid warming events (e.g., PETM and the last glacial termination) (Prokoph et al., 2008; Veizer and Prokoph, 2015). Cold climates were thus the main drivers of plant adaptation to the environment because high temperatures were relatively absent. These also contributed to the strong correlations between physiological cold tolerance and MTCM and the weak correlations between heat tolerance and MTWM. Third, the ancestors of most major clades of land plants (e.g., angiosperms) could have originated in warm-climate areas, but the descendants migrated and radiated into cooler areas; therefore, cold tolerance evolution is needed. One classic example is the poleward expansion of angiosperms during the Cenozoic (Zanne et al., 2014). Fourth, the evolution of climatic adaptation might also be associated with the macroevolution and diversification of plants. The younger age and higher speciation rates of plant taxa in cooler climates and higher latitudes indicate their potential higher selective pressure in cool areas (Lu et al., 2018; Igea and Tanentzap, 2020), further suggesting the cold climates may be the main driver in the macroevolution of plants. Some of the ideas are not exclusive of each other; however, they suggest alternative but potentially overlapping paths explaining faster cold tolerance evolution in plants.
4.3 Plants may be more sensitive to projected warming but still affected by cold hazards
Lower evolutionary rates of heat tolerance, as well as its narrow variance (Araujo et al., 2013; Lancaster and Humphreys, 2020), indicate that plants may have difficulty adapting to warmer climates. Our results also revealed that even for hardened plants, the rates of heat tolerance were not significantly higher than MTWM, indicating the limited plasticity of heat tolerance and higher risk of plants to climate warming. Particularly for tropical plants, existence in a uniformly warm climate after the last glacial maximum is conservative to heat tolerance and with increased risk of climate warming (Costion et al., 2015; Perez and Feeley, 2020; Tierney et al., 2020).
The rates of current warming are considerably higher than they had been at any geological time in the past, whereas the evolutionary rates of climatic niches are lower than those of climate change (Quintero and Wiens, 2013; Jezkova and Wiens, 2016). As the current global mean temperature is near the maximum value estimated for the past 1.2 million years (Foster et al., 2017; Steffen et al., 2018), such higher temperatures and warming are unnatural and may impact the global ecosystem similar to the past rapid warming events in geological history (Jaramillo et al., 2010; Bellard et al., 2012; Wing and Currano, 2013; Smith et al., 2020). In addition to heat, the cold hazards may affect plants, especially populations of recently latitudinally or altitudinally migrated plants, which would experience more frequent cold hazards than native habitats. Cold events are projected to occur with high anomaly (e.g., occurrence of frost events in spring after unusual warming) (Ipcc, 2014; Vitasse et al., 2018). Cold hazard affects the physiological performance of plants and slows down their migration to cooler climate areas (Wen et al., 2018). Therefore, it is necessary to assess the physiological tolerance to heat and cold tolerance of plants for their future conservation and potential risk to climate change.
5 Conclusions
The evolutionary rates of physiological cold tolerance were higher than those of the corresponding climatic niche, while the cold tolerance without hardening showed the opposite trend, suggesting that physiological tolerance constitutes the basis of the climatic niche. Physiological heat tolerance was not evolving faster than the corresponding climatic niche in both hardening and non-hardening plants. Cold tolerance evolved faster than heat tolerance in both physiological tolerance and climatic niche. Plants may be sensitive and with a high risk of future warming, and assessing the physiological tolerance of plants is a prerequisite for their successful conservation.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://figshare.com/s/06708b691f153febca98.
Author contributions
YW: Conceptualization, Writing – review & editing, Formal Analysis, Visualization, Writing – original draft. QY: Conceptualization, Writing – review & editing, Supervision. CR: Writing – review & editing. HL: Writing – review & editing. GW: Writing – review & editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by the National Natural Science Foundation of China (32201286).
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.2023.1257499/full#supplementary-material
References
Adams, D. C. (2013). Comparing evolutionary rates for different phenotypic traits on a phylogeny using likelihood. Syst. Biol. 62, 181–192. doi: 10.1093/sysbio/sys083
Araujo, M. B., Ferri-Yanez, F., Bozinovic, F., Marquet, P. A., Valladares, F., Chown, S. L. (2013). Heat freezes niche evolution. Ecol. Lett. 16, 1206–1219. doi: 10.1111/ele.12155
Armstrong, J. J., Takebayashi, N., Wolf, D. E. (2020). Cold tolerance in the genus Arabidopsis. Am. J. Bot. 107, 489–497. doi: 10.1002/ajb2.1442
Bellard, C., Bertelsmeier, C., Leadley, P., Thuiller, W., Courchamp, F. (2012). Impacts of climate change on the future of biodiversity. Ecol. Lett. 15, 365–377. doi: 10.1111/j.1461-0248.2011.01736.x
Bennett, J. M., Calosi, P., Clusella-Trullas, S., Martinez, B., Sunday, J., Algar, A. C., et al. (2018). GlobTherm, a global database on thermal tolerances for aquatic and terrestrial organisms. Sci. Data 5, 180022. doi: 10.1038/sdata.2018.22
Buchner, O., Roach, T., Gertzen, J., Schenk, S., Karadar, M., Stöggl, W., et al. (2017). Drought affects the heat-hardening capacity of alpine plants as indicated by changes in xanthophyll cycle pigments, singlet oxygen scavenging, α-tocopherol and plant hormones. Environ. Exp. Bot. 133, 159–175. doi: 10.1016/j.envexpbot.2016.10.010
Bush, A., Catullo, R. A., Mokany, K., Thornhill, A. H., Miller, J. T., Ferrier, S. (2018). Truncation of thermal tolerance niches among Australian plants. Global Ecol. Biogeogr. 27, 22–31. doi: 10.1111/geb.12637
Butler, M. A., King, A. A. (2004). Phylogenetic comparative analysis: A modeling approach for adaptive evolution. Am. Nat. 164, 683–695. doi: 10.1086/426002
Carta, A., Peruzzi, L., Ramirez-Barahona, S. (2022). A global phylogenetic regionalization of vascular plants reveals a deep split between Gondwanan and Laurasian biotas. New Phytol. 233, 1494–1504. doi: 10.1111/nph.17844
Cayuela, L., Granzow-De La Cerda, Í., Albuquerque, F. S., Golicher, D. J. (2012). taxonstand: An r package for species names standardisation in vegetation databases. Methods Ecol. Evol. 3, 1078–1083. doi: 10.1111/j.2041-210X.2012.00232.x
Chen, I. C., Hill, J. K., Ohlemuller, R., Roy, D. B., Thomas, C. D. (2011). Rapid range shifts of species associated with high levels of climate warming. Science 333, 1024–1026. doi: 10.1126/science.1206432
Colwell, R. K., Rangel, T. F. (2009). Hutchinson's duality: The once and future niche. Proc. Natl. Acad. Sci. U.S.A. 106, 19651–19658. doi: 10.1073/pnas.0901650106
Corlett, R. T., Westcott, D. A. (2013). Will plant movements keep up with climate change? Trends Ecol. Evol. 28, 482–488. doi: 10.1016/j.tree.2013.04.003
Costion, C. M., Simpson, L., Pert, P. L., Carlsen, M. M., John Kress, W., Crayn, D. (2015). Will tropical mountaintop plant species survive climate change? Identifying key knowledge gaps using species distribution modelling in Australia. Biol. Conserv. 191, 322–330. doi: 10.1016/j.biocon.2015.07.022
Felsenstein, J. (1973). Maximum-likelihood estimation of evolutionary trees from continuous characters. Am. J. Hum. Genet. 25, 471–492.
Fick, S. E., Hijmans, R. J. (2017). WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315. doi: 10.1002/joc.5086
Foster, G. L., Royer, D. L., Lunt, D. J. (2017). Future climate forcing potentially without precedent in the last 420 million years. Nat. Commun. 8, 14845. doi: 10.1038/ncomms14845
Frechette, E., Wong, C. Y., Junker, L. V., Chang, C. Y., Ensminger, I. (2015). Zeaxanthin-independent energy quenching and alternative electron sinks cause a decoupling of the relationship between the photochemical reflectance index (PRI) and photosynthesis in an evergreen conifer during spring. J. Exp. Bot. 66, 7309–7323. doi: 10.1093/jxb/erv427
Geange, S. R., Arnold, P. A., Catling, A. A., Coast, O., Cook, A. M., Gowland, K. M., et al. (2021). The thermal tolerance of photosynthetic tissues: a global systematic review and agenda for future research. New Phytol. 229, 2497–2513. doi: 10.1111/nph.17052
Glikson, A. (2016). Cenozoic mean greenhouse gases and temperature changes with reference to the Anthropocene. Glob Chang Biol. 22, 3843–3858. doi: 10.1111/gcb.13342
Hetem, R. S., Strauss, W. M., Fick, L. G., Maloney, S. K., Meyer, L. C., Shobrak, M., et al. (2012). Activity re-assignment and microclimate selection of free-living Arabian oryx: responses that could minimise the effects of climate change on homeostasis? Zoology 115, 411–416. doi: 10.1016/j.zool.2012.04.005
Igea, J., Tanentzap, A. J. (2020). Angiosperm speciation cools down in the tropics. Ecol. Lett. 23, 692–700. doi: 10.1111/ele.13476
Ipcc (2014) Climate Change 2014: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Available at: https://www.ipcc.ch/report/ar5/syr.
Janmohammadi, M., Sabaghnia, N., Mahfoozi, S. (2018). Frost tolerance and metabolite changes of rye (Secale cereale) during the cold hardening and overwintering. Acta Physiol. Plant 40, 42. doi: 10.1007/s11738-018-2620-0
Jaramillo, C., Ochoa, D., Contreras, L., Pagani, M., Carvajal-Ortiz, H., Pratt, L. M., et al. (2010). Effects of rapid global warming at the Paleocene-Eocene boundary on neotropical vegetation. Science 330, 957–961. doi: 10.1126/science.1193833
Jezkova, T., Wiens, J. J. (2016). Rates of change in climatic niches in plant and animal populations are much slower than projected climate change. Proc. Biol. Sci. 283, 20162104. doi: 10.1098/rspb.2016.2104
Jin, Y., Qian, H. (2019). V.PhyloMaker: an R package that can generate very large phylogenies for vascular plants. Ecography 42, 1353–1359. doi: 10.1111/ecog.04434
Knaupp, M., Mishra, K. B., Nedbal, L., Heyer, A. G. (2011). Evidence for a role of raffinose in stabilizing photosystem II during freeze-thaw cycles. Planta 234, 477–486. doi: 10.1007/s00425-011-1413-0
Lancaster, L. T., Humphreys, A. M. (2020). Global variation in the thermal tolerances of plants. Proc. Natl. Acad. Sci. U.S.A. 117, 13580–13587. doi: 10.1073/pnas.1918162117
Lancaster, L. T., Humphreys, A. M. (2021). Reply to Perez et al.: Experimental duration unlikely to bias global variation in plant thermal tolerances. Proc. Natl. Acad. Sci. U.S.A. 118, e2102037118. doi: 10.1073/pnas.2102037118
Landis, J. B., Soltis, D. E., Li, Z., Marx, H. E., Barker, M. S., Tank, D. C., et al. (2018). Impact of whole-genome duplication events on diversification rates in angiosperms. Am. J. Bot. 105, 348–363. doi: 10.1002/ajb2.1060
Liu, C., Wolter, C., Xian, W., Jeschke, J. M. (2020a). Most invasive species largely conserve their climatic niche. Proc. Natl. Acad. Sci. U.S.A. 117, 23643–23651. doi: 10.1073/pnas.2004289117
Liu, H., Ye, Q., Wiens, J. J. (2020b). Climatic-niche evolution follows similar rules in plants and animals. Nat. Ecol. Evol. 4, 753–763. doi: 10.1038/s41559-020-1158-x
Lu, L. M., Mao, L. F., Yang, T., Ye, J. F., Liu, B., Li, H. L., et al. (2018). Evolutionary history of the angiosperm flora of China. Nature 554, 234–238. doi: 10.1038/nature25485
O'sullivan, O. S., Heskel, M. A., Reich, P. B., Tjoelker, M. G., Weerasinghe, L. K., Penillard, A., et al. (2016). Thermal limits of leaf metabolism across biomes. Global Change Biol. 23, 209–223. doi: 10.1111/gcb.13477
O'sullivan, O. S., Weerasinghe, K. W. L. K., Evans, J. R., Egerton, J. J. G., Tjoelker, M. G., Atkin, O. K. (2013). High-resolution temperature responses of leaf respiration in snow gum (Eucalyptus pauciflora) reveal high-temperature limits to respiratory function. Plant Cell Environ. 36, 1268–1284. doi: 10.1111/pce.12057
Orme, D. (2018) The Caper Package: Comparative Analysis of Phylogenetics and Evolution in R. Available at: https://cran.r-project.org/web/packages/caper/vignettes/caper.pdf.
Pennell, M. W., Eastman, J. M., Slater, G. J., Brown, J. W., Uyeda, J. C., Fitzjohn, R. G., et al. (2014). geiger v2.0: an expanded suite of methods for fitting macroevolutionary models to phylogenetic trees. Bioinformatics 30, 2216–2218. doi: 10.1093/bioinformatics/btu181
Pereira, A. M. F., Vilela, R. A., Titto, C. G., Leme-Dos-Santos, T. M. C., Geraldo, A. C. M., Balieiro, J. C. C., et al. (2020). Thermoregulatory responses of heat acclimatized buffaloes to simulated heat waves. Animals 10, 756. doi: 10.3390/ani10050756
Perez, T. M., Feeley, K. J. (2020). Weak phylogenetic and climatic signals in plant heat tolerance. J. Biogeogr. 48, 91–100. doi: 10.1111/jbi.13984
Prokoph, A., Shields, G. A., Veizer, J. (2008). Compilation and time-series analysis of a marine carbonate δ18O, δ13C, 87Sr/86Sr and δ34S database through Earth history. Earth Sci. Rev. 87, 113–133. doi: 10.1016/j.earscirev.2007.12.003
Qu, Y. F., Wiens, J. J. (2020). Higher temperatures lower rates of physiological and niche evolution. Proc. Biol. Sci. 287, 20200823. doi: 10.1098/rspb.2020.0823
Quintero, I., Wiens, J. J. (2013). Rates of projected climate change dramatically exceed past rates of climatic niche evolution among vertebrate species. Ecol. Lett. 16, 1095–1103. doi: 10.1111/ele.12144
R Core Team (2020) R: A Language and Environment for Statistical Computing. Available at: https://www.R-project.org.
Revell, L. J. (2012). phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223. doi: 10.1111/j.2041-210X.2011.00169.x
Smith, S. A., Brown, J. W. (2018). Constructing a broadly inclusive seed plant phylogeny. Am. J. Bot. 105, 302–314. doi: 10.1002/ajb2.1019
Smith, M. N., Taylor, T. C., Van Haren, J., Rosolem, R., Restrepo-Coupe, N., Adams, J., et al. (2020). Empirical evidence for resilience of tropical forest photosynthesis in a warmer world. Nat. Plants 6, 1225–1230. doi: 10.1038/s41477-020-00780-2
Steffen, W., Rockstrom, J., Richardson, K., Lenton, T. M., Folke, C., Liverman, D., et al. (2018). Trajectories of the earth system in the anthropocene. Proc. Natl. Acad. Sci. U.S.A. 115, 8252–8259. doi: 10.1073/pnas.1810141115
Tierney, J. E., Poulsen, C. J., Montanez, I. P., Bhattacharya, T., Feng, R., Ford, H. L., et al. (2020). Past climates inform our future. Science 370. doi: 10.1126/science.aay3701
Veizer, J., Prokoph, A. (2015). Temperatures and oxygen isotopic composition of Phanerozoic oceans. Earth Sci. Rev. 146, 92–104. doi: 10.1016/j.earscirev.2015.03.008
Vitasse, Y., Schneider, L., Rixen, C., Christen, D., Rebetez, M. (2018). Increase in the risk of exposure of forest and fruit trees to spring frosts at higher elevations in Switzerland over the last four decades. Agric. For. Meteorol. 248, 60–69. doi: 10.1016/j.agrformet.2017.09.005
Wen, Y., Qin, D. W., Leng, B., Zhu, Y. F., Cao, K. F. (2018). The physiological cold tolerance of warm-climate plants is correlated with their latitudinal range limit. Biol. Lett. 14, 20180277. doi: 10.1098/rsbl.2018.0277
Wing, S. L., Currano, E. D. (2013). Plant response to a global greenhouse event 56 million years ago. Am. J. Bot. 100, 1234–1254. doi: 10.3732/ajb.1200554
Zanne, A. E., Tank, D. C., Cornwell, W. K., Eastman, J. M., Smith, S. A., Fitzjohn, R. G., et al. (2014). Three keys to the radiation of angiosperms into freezing environments. Nature 506, 89–92. doi: 10.1038/nature13842
Zeebe, R. E., Zachos, J. C., Dickens, G. R. (2009). Carbon dioxide forcing alone insufficient to explain Palaeocene-Eocene Thermal Maximum warming. Nat. Geosci. 2, 576–580. doi: 10.1038/Ngeo578
Keywords: physiological tolerance, climatic niche, climate change, evolution, heat tolerance
Citation: Wen Y, Ye Q, Román-Palacios C, Liu H and Wu G (2023) Physiological cold tolerance evolves faster than climatic niches in plants. Front. Plant Sci. 14:1257499. doi: 10.3389/fpls.2023.1257499
Received: 12 July 2023; Accepted: 21 August 2023;
Published: 08 September 2023.
Edited by:
Runguo Zang, Chinese Academy of Forestry, ChinaReviewed by:
Peili Mao, Shandong Agricultural University, ChinaZhicheng Chen, Chinese Academy of Forestry, China
Copyright © 2023 Wen, Ye, Román-Palacios, Liu and Wu. 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: Qing Ye, qye@scbg.ac.cn