Abstract
Lichen associations are overwhelmingly supported by carbon produced by photosynthetic algal symbionts. These algae have diversified to occupy nearly all climates and continents; however, we have a limited understanding of how their climatic niches have evolved through time. Here we extend previous work and ask whether phylogenetic signal in, and the evolution of, climatic niche, varies across climatic variables, phylogenetic scales, and among algal lineages in Trebouxia—the most common genus of lichen-forming algae. Our analyses reveal heterogeneous levels of phylogenetic signal across variables, and that contrasting models of evolution underlie the evolution of climatic niche divergence. Together these analyses demonstrate the variable processes responsible for shaping climatic tolerance in Trebouxia, and provide a framework within which to better understand potential responses to climate change-associated perturbations. Such predictions reveal a disturbing trend in which the pace at which modern climate change is proceeding will vastly exceed the rate at which Trebouxia climatic niches have previously evolved.
Introduction
Closely related species have a tendency to resemble one another as a consequence of their shared evolutionary history (Darwin, 1859; Blomberg and Garland, 2002). Thus, related lineages may occupy similar ecological niches and overlap in their climatic tolerance. However, the extent to which lineages resemble one another may vary both by trait and across clades as a result of variable evolutionary processes. As climate change has, and is predicted to, disrupt ecosystems and force ranges to shift, it is important to understand which lineages may be especially vulnerable. This vulnerability, is in part, a product of climatic preference. Lineages with slow rates of climatic niche evolution may have a more limited ability to adapt to and persist through climatic perturbations than those with a legacy of evolutionary lability (Holt, 1990; Wiens et al., 2010). Indeed, estimates of predicted climate change often vastly outpace past rates of climatic niche evolution in diverse lineages, suggesting these species may not be able to adapt to the rapid climate change that is underway (Quintero and Wiens, 2013; Jezkova and Wiens, 2016). Thus, it is important to understand how the climatic niche of a lineage has evolved through time in order to identify its potential sensitivity or ability to adapt to shifting climates.
Symbioses are expected to be especially vulnerable to climate change (Dunn et al., 2009; Colwell et al., 2012; Kikuchi et al., 2016; Sampayo et al., 2016; Baker et al., 2018; Rolshausen et al., 2018). Lichens—symbiotic associations between heterotrophic fungi and photosynthetic microbes (eukaryotic microalgae and/or cyanobacteria)—make important and diverse ecosystem contributions (Seaward, 2008; Elbert et al., 2012; Porada et al., 2014, 2018; Rutherford et al., 2017), and are well-known bioindicators (Nimis et al., 2002; Giordani and Incerti, 2008; Will-Wolf et al., 2015; Loppi, 2019; Will-Wolf and Jovan, 2019) with certain lichens predicted to exhibit strong responses to changing climates (Allen and Lendemer, 2016; Ellis, 2019; Nelsen and Lumbsch, 2020). Thus, understanding the evolutionary legacy of climatic niche evolution in these symbionts is paramount to predicting their extinction vulnerability and potential to adapt to changing climates.
The eukaryotic algal genus Trebouxia is the most common lichen photobiont (Tschermak-Woess, 1988; Rambold et al., 1998; Voytsekhovich et al., 2011; Miadlikowska et al., 2014; Nelsen et al., 2021). It occurs in nearly all terrestrial ecosystems, but is especially frequent in more temperate-boreal/alpine ecosystems (Honegger, 1998, 2009; Sanders, 2001; Pérez-Ortega et al., 2012; Ruprecht et al., 2012; Wagner et al., 2020; Nelsen et al., 2021). Global efforts to characterize the environmental preferences of Asterochloris, a Trebouxia relative, have revealed the differential climatic tolerance of individual species-level groups, while suggesting closely related accessions shared similar environmental preferences and exhibited strong phylogenetic conservatism (Peksa and Škaloud, 2011; Vančurová et al., 2018; Kosecka et al., 2021). Previous multivariate efforts in Trebouxia have revealed coarse-scaled trends such that members of a specific clade (Trebouxia clade “C”) tended to exhibit a more limited climatic preference than others, and a joint analysis of two climatic variables suggested rates of climatic niche diversification across Trebouxia have remained constant through time (Nelsen et al., 2021). However, these two variables (mean annual temperature and precipitation) exhibited similar responses in a phylogenetic principal components analysis, and it remains unclear whether preference toward other climatic variables exhibit similar trends and were shaped by similar evolutionary processes. Here we extend previous work (Nelsen et al., 2021) and ask whether climatic niche evolution varies across climatic variables. These past estimates are then used to make predictions about projected responses to climate change.
Materials and Methods
General
Here we focused on characterizing phylogenetic signal in, and the evolution of, climatic niche, in the green algal genus Trebouxia. While comparable analyses of the over 7,000 fungal species associated with Trebouxia algae (Nelsen et al., 2021) would certainly be of interest, this was regarded as beyond the scope of the present study. We extended previous work (Nelsen et al., 2021) by characterizing the phylogenetic signal and evolution of climatic preferences for Trebouxia operational taxonomic units (OTUs) across a set of phylogenies, thereby accounting for temporal and topological uncertainty. Briefly, we utilized a previous dataset (Nelsen et al., 2021) that subjected over 6,000 publicly available Trebouxia ITS sequences to quality control filters, binned accessions into OTUs, and inferred time-scaled phylogenies using a previously published secondary calibration (Nelsen et al., 2020). Here we utilized 1,001 time-scaled phylogenies (Maximum Clade Credibility [MCC] tree, and a random set of 1,000 trees derived from the posterior) of Trebouxia OTUs delimited at a sequence similarity of 97.5% (Nelsen et al., 2021). Geographic coordinates for individual accessions passing quality control filters were previously retrieved from NCBI or obtained from the primary literature and estimates for each of 19 bioclimatic variables (Vega et al., 2017) were extracted for each of the 2,359 accessions representing 81 OTUs (min = 1 accession/OTU, mean = 29.1 accessions/OTU, max = 454 accessions/OTU) (Nelsen et al., 2021). Mean estimates of each bioclimatic variable for individual OTUs were used in analyses below (Nelsen et al., 2021), and sampling error estimated from raw values (described below). We note that these are not experimentally derived climatic tolerances, and the observed distribution of Trebouxia accessions may be due to other factors including but not limited to biotic interactions such as competition, and sampling bias. However, in the absence of experimental studies characterizing tolerance of individual OTUs to thermal and precipitation extremes and variability—both in and out of the lichenized state—we regard these as provisional proxies for their climatic preference.
Previous work identified contrasting correlations among climatic variables and principal component (PC) axes (Nelsen et al., 2021). These variables can be binned into six groups (Supplementary Figure 1). Three variable groups exhibited strongly negative loadings along PC1: Group 1 represented annual mean temperature and specific humidity (BIO1, BIO12); Group 2 consisted of variables representing mean temperature and specific humidity during the warmest and most humid periods (BIO5, BIO8, BIO10, BIO13, BIO16, BIO18); Group 3 was comprised of variables representing mean temperature and specific humidity during the coldest and least humid periods (BIO6, BIO9, BIO11, BIO14, BIO17, BIO19). Two sets of variables were characterized by strong negative loadings along PC2: Group 4 was composed of variables representing temperature seasonality and range (BIO4, BIO7) and Group 5 contained one variable related to specific humidity seasonality (BIO15). Group 6 included one variable linked to isothermality (BIO3) and was most strongly associated with PC3. Finally, mean diurnal range of temperature (BIO2) exhibited an equally negative loading along PC2 and PC3. Here we assess whether patterns of phylogenetic signal and trait evolution are consistent among individual variables within each of these groupings, as well as among groupings. Analyses were performed in the R programming environment (R Core Team, 2014) at the Grainger Bioinformatics Center (Field Museum).
Phylogenetic Signal in Climatic Niche
We first asked to what extent closely related Trebouxia OTUs resembled one another in their climatic niches by quantifying phylogenetic signal individually for 19 climatic variables. Phylogenetic signal represents the tendency (pattern) for related species to resemble one another as expected under a model of random drift or Brownian motion (BM) (Blomberg and Garland, 2002). The degree of phylogenetic signal in climatic preference was estimated for each climatic variable by calculating the K statistic (Blomberg et al., 2003) for mean climate values in phytools (Revell, 2012). K-values of 0 indicate trait evolution proceeded independently of the phylogeny, while K-values of 1 reflect the expectation under pure BM/random drift along the phylogeny. Values greater than 1 indicate that trait values of closely related taxa are more similar than expected under a BM model. Significance tests were then conducted to determine whether K-values were significantly different from K = 0 through comparison with 1,000 randomizations. We included within-OTU variance estimates, but as sample sizes for many OTUs were equal to one, our OTU-specific sampling error was instead derived by calculating the pooled variance and dividing by the square root of the sample size for each OTU (Garamszegi, 2014; Slater and Friscia, 2019). Analyses were conducted for all bioclimatic variables using the MCC tree and 1,000 trees randomly sampled from the posterior, thereby integrating over topological and temporal uncertainty.
Climatic Niche Disparity Through Time
We conducted genus-level disparity through time (DTT) analyses on each climatic variable in geiger (Pennell et al., 2014), which allowed us to ascertain whether relative climatic niche disparity deviated substantially from that expected under BM. Disparity was measured using the average squared Euclidean distance, and the morphological disparity index (MDI) estimated, which reflects the difference in relative disparity of the observed clade with that expected under BM. Negative MDI values indicate that subclade disparity is lower than expected under BM, and is a commonly observed pattern in adaptive radiations, such that trait diversification is concentrated disproportionately early in the phylogeny (Harmon et al., 2003). Significance was assessed by estimating the probability of obtaining a negative MDI value, and was achieved by calculating the proportion of instances in which the relative disparity of the observed clade was equal to or exceeded that observed in each of 1,000 datasets simulated under Brownian motion (BM) (Slater et al., 2010). Analyses were again conducted on the MCC tree and 1,000 topologies derived from the posterior distribution.
Climatic Niche Evolution
While the estimates of phylogenetic signal measured above quantify an observed pattern, they do not provide insight into the underlying evolutionary and non-evolutionary processes giving rise to the observed signal (Blomberg and Garland, 2002; Revell et al., 2008; Münkemüller et al., 2012). For instance, phylogenetic signal may be a consequence of multiple evolutionary processes (Revell et al., 2008). Consequently, we fit three explicit models of trait evolution to narrow the potential underlying evolutionary processes structuring climatic niche evolution. These three models were used to determine whether the rate of trait evolution has remained constant, increased or decreased through time, and included: (1) a BM or random walk model in which trait variance between pairs of taxa is determined by the amount of shared ancestry, and the rate of trait evolution remains constant through time (Felsenstein, 1973); (2) an Ornstein-Uhlenbeck (OU) model, which is essentially a BM model in which traits are pulled toward an adaptive optimum (Butler and King, 2004). However, on an ultrametric tree, this model is also identical to one in which trait evolution exponentially increases through time (Slater et al., 2012; Uyeda et al., 2015); consequently, an OU model should be favored if rates of evolution are increasing. Finally, we fit (3) an early burst (EB) model in which trait evolution decreases exponentially through time (Harmon et al., 2010). Analyses were conducted at the genus-level, and we included within-OTU variance estimates as described above. Models were fit in geiger (Pennell et al., 2014) for each variable on the MCC tree and a random set of 1,000 trees derived from the posterior, and fits assessed using relative weights of Akaike Information Criterion corrected for small samples (AICcw).
In addition, we compared diffusion rate (σ2) estimates from two pairs of thermal and moisture variables to determine whether tolerance to minima and maxima evolve at different rates (Araújo et al., 2013; Liu et al., 2020). As in previous work (Liu et al., 2020), thermal tolerance was assessed by comparing rates from BIO5 (Max Temperature of Warmest Month) and BIO6 (Min Temperature of Coldest Month), while moisture tolerance was assessed by comparing rates from BIO16 (Specific Humidity Mean of Most Humid Quarter) and BIO17 (Specific Humidity Mean of Least Humid Quarter). Rate estimates were derived from the best-fit models (above) for the MCC tree and 1,000 trees derived from the posterior. Rates for the MCC tree are reported, together with 95% highest posterior density (HPD) values from across the distribution of trees.
Finally, we estimated absolute rates of climatic niche change for three thermal variables (BIO1: Mean Annual Temperature; BIO5: Max Temperature of Warmest Month; BIO6: Min Temperature of Coldest Month). We then compared projected thermal changes over the next 80 years due to global heating (IPCC, 2014) with historic rates of thermal niche evolution in Trebouxia algae to assess whether these algae are likely to adapt in place to increasing temperatures, or if they will be required to migrate to suitable climates. Absolute rates of change were calculated using an approach similar to that applied in previous work (Quintero and Wiens, 2013; Jezkova and Wiens, 2016) for the MCC tree and a random set of 1,000 trees derived from the posterior. For each trait, we transformed individual trees using the best-fit model of evolution and parameter estimates (above) in geiger (Pennell et al., 2014), and subsequently reconstructed ancestral states in phytools. We then calculated the absolute difference between the tip state and its immediate ancestor, and divided this by the time separating the two. In addition, minimum and maximum variable values for each OTU were extracted from previous work (Nelsen et al., 2021) and used to calculate absolute rates of change. Rates for the MCC tree are reported, together with 95% highest posterior density (HPD) values from across the distribution of trees. The MCC tree was plotted together with rate estimates in phytools.
Results
We used the K-statistic to evaluate the extent to which climatic preferences of related OTUs resembled one another, relative to those drawn at random. The K-values obtained typically ranged from 0.3 to 1 (Figure 1) and were significant in all cases (except for some BIO3 analyses), indicating that climatic preferences were not independent of phylogeny, and those of closely related taxa tended to be more similar to one another than to those of distantly related taxa. Within climatic variable groups, values were broadly consistent. K-values for traits from climatic variable Groups 1 (BIO1, BIO12), 2 (BIO5, BIO8, BIO10, BIO13, BIO16, BIO18), and 3 (BIO6, BIO9, BIO11, BIO14, BIO17, BIO19) were all consistent with expectations under BM (K = 1) or when trait values for closely related lineages are more similar than expected under a BM model (K > 1). In contrast, K-values for climatic variable Groups 4 (BIO4, BIO7), 5 (BIO15), 6 (BIO3), and BIO2 were lower than expected under pure BM (K = 1).
FIGURE 1
Our genus-wide disparity through time analyses for individual climatic variables further demonstrated that disparity in these traits was not disproportionately concentrated early in the phylogeny. Mean disparity index (MDI) values for the MCC (Figure 2 and Supplementary Figures 2, 3) and the set of trees derived from the posterior (Supplementary Table 1) were not significantly lower than expected under BM, and therefore, inconsistent with an EB model. In contrast, climatic variable Groups 5 (BIO15), 6 (BIO3), BIO2, and a member of Group 4 (BIO4) regularly yielded MDI values significantly more positive (1—p-value, since test is one-sided) than expected under BM, and instead consistent with an OU model of evolution. The MDI values and significance levels achieved for the other variable in climatic variable Group 4 (BIO7) were also consistent with an OU model of evolution, but significance and MDI values were not as extreme as those for BIO4.
FIGURE 2
FIGURE 3
Our subsequent model-fitting procedure confirmed that a BM model of evolution best characterized the evolution of traits in climatic variable Groups 1 (BIO1, BIO12), 2 (BIO5, BIO8, BIO10, BIO13, BIO16, BIO18), and 3 (BIO6, BIO9, BIO11, BIO14, BIO17, BIO19) (Figure 3). In contrast The OU model was favored for climatic variable Groups 4 (BIO4, BIO7), 5 (BIO15), and 6 (BIO3), while support for an OU model was equal to, or slightly higher than that of BM for climatic variable BIO2. Diffusion rate estimates (σ2) for targeted thermal (BIO5, BIO6) and moisture (BIO16, BIO17) variables derived from the best-fit models of evolution (Figure 3) generally suggest that thermal maxima (BIO5) evolved at a slower rate than thermal minima (BIO6), and moisture maxima (BIO16) largely evolved at a faster rate than moisture minima (BIO17) (Figure 4). While there is some overlap among rate estimates within each variable class, distributions appear largely segregated. Finally, our absolute rate estimates of climatic niche evolution varied among OTUs and across variables, but were generally less than 1°C/Myr, regardless of thermal values and approaches employed (Figure 5 and Supplementary Tables 2–4).
FIGURE 4
FIGURE 5
Discussion
Significant phylogenetic signal at the genus level was detected across nearly all variables; however, the magnitude of this signal varied substantially across climatic variables. All variables linked to PC1 (BIO1, BIO12, BIO5, BIO8, BIO10, BIO13, BIO16, BIO18, BIO6, BIO9, BIO11, BIO14, BIO17, BIO19) in previous work (Nelsen et al., 2021) yielded K-values consistent with a BM expectation (K = 1), DTT plots and MDI values consistent with BM, and had a BM model of evolution overwhelmingly favored in our model-fitting analyses. These climatic variables reflect averages and extreme thermal and specific humidity values, and suggest that at this phylogenetic scale, Trebouxia algae have evolved thermal and moisture tolerances under a neutral evolutionary process and at a relatively constant rate.
In contrast, lower (K < 1) or non-significant levels of signal were detected at the genus level for a limited number of climatic variables—all (BIO4, BIO7, BIO15, BIO2, BIO3) of which were linked to PC2 and PC3 in previous work (Nelsen et al., 2021), and reflected thermal or moisture variability (seasonally, diurnally). These variables also yielded DTT plots inconsistent with BM, and MDI values consistent with a non-BM process. Moreover, our model-fitting analyses suggested an OU model was almost always overwhelmingly favored for these variables. Together, this may indicate that, at this phylogenetic scale, the evolution of thermal and moisture variability tolerance may be under selection and evolve toward an adaptive optimum. However, it may instead suggest the rate of trait evolution is exponentially increasing through time (AC model, Blomberg et al., 2003), as these two models are equivalent on an ultrametric tree lacking fossil taxa (Slater et al., 2012; Uyeda et al., 2015); thus, we are unable to distinguish among these two patterns. Regardless, it does reveal a departure from Brownian motion and, together with our findings for other climatic variables, indicates variable evolutionary processes underlie climatic niche evolution in this clade of keystone symbionts.
Our DTT plots, MDI values and model-fitting analyses also demonstrate that at this phylogenetic scale, climatic preference under an exponentially decreasing rate of evolution (EB) consistent with that of an adaptive radiation is highly unlikely. Thus, clades have not evolved exceptionally distinct climatic preferences at deep phylogenetic scales, and instead exhibit overlapping climatic tolerances. This is consistent with previous analyses demonstrating the presence of multiple clades in the same geographic locality (Doering and Piercey-Normore, 2009; Werth, 2012; Leavitt et al., 2013, 2015), and observations that climatic niche divergence of annual mean temperature (BIO1) and specific humidity (BIO12) did not jointly evolve under an EB process (Nelsen et al., 2021). Extending our analyses to include numerous individual climatic variables further demonstrates that at the genus level, early bursts in climatic niche tolerance appear rare.
Our analyses of the evolution of thermal and moisture climatic niche limits suggest an asymmetry such that upper thermal limits of Trebouxia algae evolve at a slower rate relative to lower thermal limits, and that upper moisture limits evolve at a more rapid rate than lower thermal limits. Both of these findings are consistent with previous work. For instance, a greater range of lower thermal limits relative to upper limits was detected in animals (Araújo et al., 2013), and subsequent work integrating evolutionary the age and relationships of a diverse range of animals and plants has further revealed phylogenetic constraints on thermal limits, with lower limits evolving faster than upper limits (Muñoz et al., 2014; Culumber and Tobler, 2018; Liu et al., 2020; Bennett et al., 2021). We emphasize that these limits may represent hard physiological boundaries that may restrict the potential for lineages to adapt in situ to changing climates. Together, these findings suggest the limited potential for several lineages to evolve elevated upper thermal limits in response to increasing global temperatures (Colwell et al., 2008; Liu et al., 2020; Bennett et al., 2021). Similarly, previous work on a range of plants and animals has demonstrated that lower moisture limits evolve at a slower rate than upper limits, suggesting that it may be more difficult to adapt to increased aridity than to increased moisture (Liu et al., 2020). Consistent with other systems, our analyses of the evolutionary legacies of Trebouxia algae suggest a limited capacity to adapt in situ to increased temperature and decreased moisture (relative to decreased temperature and increased moisture). However, we note that lower thermal limits frequently decrease with increasing latitude, while upper thermal limits remain constant across latitudes (Addo-Bediako et al., 2000); this trend could in turn underlie the observed asymmetric rates of thermal niche evolution, as thermal minima may exhibit greater global variation relative to thermal maxima (Liu et al., 2020). We thus encourage further work to disentangle the roles intrinsic and extrinsic variables have played in structuring these discordant rates of evolution.
Our estimates of absolute rates of climatic niche evolution suggested Trebouxia OTUs have evolved at a relatively slow pace, on order with some of the slower rates previously estimated from a diverse range of animals and plants (Quintero and Wiens, 2013; Jezkova and Wiens, 2016). Given that the global heating we are currently experiencing is expected to increase temperatures by 1–4°C (relative to 1986–2005) by the year 2100 (IPCC, 2014)—less than 80 years from now, these algae with their comparatively slow rates of thermal niche evolution [mean rate estimates of 0.13–0.62°C/Myr (BIO1), 0.07–0.71°C/Myr (BIO5), 0.1–1.08°C/Myr (BIO6)] will likely undergo range modification. This suggests that in certain portions of an OTU’s range, individuals are unlikely to adapt in situ to increased temperatures experienced at such a rapid rate, and will instead be required to disperse to more habitable locations, thereby altering their distribution. Previous work has emphasized species migration as a means to persist in the face of changing climates; for instance, the Quaternary record illustrates how species ranges have changed in response to shifting climates (Ackerly, 2003). However, this migration is likely to also affect lichen fungi that depend on these algae for nutrition. Such disruptions of biotic interactions are expected to be one of the major drivers of climate change-mediated local extinction (Cahill et al., 2013), with symbioses in particular anticipated to be especially vulnerable to the effects of climate change (Dunn et al., 2009; Colwell et al., 2012; Kikuchi et al., 2016; Sampayo et al., 2016; Baker et al., 2018; Rolshausen et al., 2018). Thus, algal range shifts may further disrupt interaction webs and force fungi to also disperse or acquire new algal symbionts. We emphasize that our climatic niche variables are not direct estimates of the physiological limits of species, but are instead regarded as proxies. Similarly, our characterization of their niches only encompassed their realized niche, rather than their fundamental niche. We therefore stress the need for detailed studies to identify the physiological limits of these algae, both in culture and in the lichenized state. Despite this, our analyses of past rates of climatic niche evolution paint a troubling picture by suggesting a limited evolutionary capacity for Trebouxia algae to adapt in situ to rapidly increasing temperatures.
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.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. Data and code used here are provided at: https://github.com/mpnelsen/Trebouxia_Niche_Divergence.
Author contributions
MN, SL, KH, LM, and HL designed research and wrote the manuscript. MN performed research. MN and KH collected, analyzed, and interpreted data. All authors contributed to the article and approved the submitted version.
Funding
The Field Museum of Natural History, Negaunee Foundation (Field Museum), Grainger Bioinformatics Center (Field Museum), and University of Chicago Odyssey Scholars program are kindly acknowledged for support.
Acknowledgments
MN thanks Graham Slater and Dean Adams for discussion of phylogenetic comparative methods. Three reviewers are gratefully acknowledged for their comments, which strengthened the manuscript.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.791546/full#supplementary-material
Supplementary Figure 1The first two axes of the phylogenetic principal component analysis (pPCA) of 19 bioclimatic variables from Nelsen et al. (2021) with OTUs shaded by clade. Axes are labeled with the proportion of variance explained by each axis, and correlations with bioclimatic variables (1–19) are indicated with gray lines. Climatic variables can be binned into six groups with similar correlations to individual axes. Group 1 includes: Annual Mean Temperature (BIO1), Annual Mean Specific Humidity (BIO12). Climatic variable Group 2 includes: Max Temperature of Warmest Month (BIO5), Mean Temperature of Most Humid Quarter (BIO8), Mean Temperature of Warmest Quarter (BIO10), Specific Humidity of Most Humid Month (BIO13), Specific Humidity Mean of Most Humid Quarter (BIO16), Specific Humidity Mean of Warmest Quarter (BIO18). Climatic variable Group 3 includes: Min Temperature of Coldest Month (BIO6), Mean Temperature of Least Humid Quarter (BIO9), Mean Temperature of Coldest Quarter (BIO11), Specific Humidity of Least Humid Month (BIO14), Specific Humidity Mean of Least Humid Quarter (BIO17), Specific Humidity Mean of Coldest Quarter (BIO19). Climatic variable Group 4 includes: Temperature Seasonality (Standard Deviation *100) (BIO4), Temperature Annual Range (BIO5–BIO6) (BIO7). Climatic variable Group 5 is limited to Specific Humidity Seasonality (Coefficient of Variation) (BIO15), while Group 6 is restricted to Isothermality (BIO2/BIO7) (* 100) (BIO3). It is unclear which group Mean Diurnal Range Temperature (BIO2) belongs to. Figure originally published in Nelsen et al. (2021) and reproduced here with permission from Oxford University Press.
Supplementary Figure 2Disparity through time (DTT) plot of BIO7–BIO12 variables. The solid line indicates the observed subclade disparity. The dashed line indicates median estimates from simulated data, and the gray cloud indicates the 95% confidence interval derived from simulated data. The MDI value and its associated p-value are provided. Vertical shading indicates geologic periods, which are abbreviated at the bottom.
Supplementary Figure 3Disparity through time (DTT) plot of BIO13–BIO19 variables. The solid line indicates the observed subclade disparity. The dashed line indicates median estimates from simulated data, and the gray cloud indicates the 95% confidence interval derived from simulated data. The MDI value and its associated p-value are provided. Vertical shading indicates geologic periods, which are abbreviated at the bottom.
References
1
AckerlyD. D. (2003). Community assembly, niche conservatism, and adaptive evolution in changing environments.Int. J. Plant Sci.164S165–S184. 10.1086/368401
2
Addo-BediakoA.ChownS. L.GastonK. J. (2000). Thermal tolerance, climatic variability and latitude.Proc. R. Soc. Lond. Ser. B Biol. Sci.267739–745. 10.1098/rspb.2000.1065
3
AllenJ. L.LendemerJ. C. (2016). Climate change impacts on endemic, high-elevation lichens in a biodiversity hotspot.Biodivers. Conserv.25555–568. 10.1007/s10531-016-1071-4
4
AraújoM. B.Ferri-YáñezF.BozinovicF.MarquetP. A.ValladaresF.ChownS. L. (2013). Heat freezes niche evolution.Ecol. Lett.161206–1219. 10.1111/ele.12155
5
BakerD. M.FreemanC. J.WongJ. C. Y.FogelM. L.KnowltonN. (2018). Climate change promotes parasitism in a coral symbiosis.ISME J.12921–930. 10.1038/s41396-018-0046-8
6
BennettJ. M.SundayJ.CalosiP.VillalobosF.MartínezB.Molina-VenegasR.et al (2021). The evolution of critical thermal limits of life on Earth.Nat. Commun.12:1198. 10.1038/s41467-021-21263-8
7
BlombergS. P.GarlandT. (2002). Tempo and mode in evolution: phylogenetic inertia, adaptation and comparative methods.J. Evol. Biol.15899–910. 10.1046/j.1420-9101.2002.00472.x
8
BlombergS. P.GarlandT.IvesA. R. (2003). Testing for phylogenetic signal in comparative data: behavioral traits are more labile.Evolution57717–745. 10.1111/j.0014-3820.2003.tb00285.x
9
ButlerM. A.KingA. A. (2004). Phylogenetic comparative analysis: a modeling approach for adaptive evolution.Am. Natural.164683–695. 10.1086/426002
10
CahillA. E.Aiello-LammensM. E.Fisher-ReidM. C.HuaX.KaranewskyC. J.Yeong RyuH.et al (2013). How does climate change cause extinction?Proc. R. Soc. B Biol. Sci.280:20121890. 10.1098/rspb.2012.1890
11
ColwellR. K.BrehmG.CardelúsC. L.GilmanA. C.LonginoJ. T. (2008). Global warming, elevational range shifts, and lowland biotic attrition in the wet tropics.Science322258–261. 10.1126/science.1162547
12
ColwellR. K.DunnR. R.HarrisN. C. (2012). Coextinction and persistence of dependent species in a changing world.Annu. Rev. Ecol. Evol. Syst.43183–203. 10.1146/annurev-ecolsys-110411-160304
13
CulumberZ. W.ToblerM. (2018). Correlated evolution of thermal niches and functional physiology in tropical freshwater fishes.J. Evol. Biol.31722–734. 10.1111/jeb.13260
14
DarwinC. (1859). On the Origin of Species.London: J. Murray.
15
DoeringM.Piercey-NormoreM. D. (2009). Genetically divergent algae shape an epiphytic lichen community on Jack Pine in Manitoba.Lichenologist4169–80. 10.1017/S0024282909008111
16
DunnR. R.HarrisN. C.ColwellR. K.KohL. P.SodhiN. S. (2009). The sixth mass coextinction: are most endangered species parasites and mutualists?Proc. R. Soc. B Biol. Sci.2763037–3045. 10.1098/rspb.2009.0413
17
ElbertW.WeberB.BurrowsS.SteinkampJ.BüdelB.AndreaeM. O.et al (2012). Contribution of cryptogamic covers to the global cycles of carbon and nitrogen.Nat. Geosci.5459–462. 10.1038/ngeo1486
18
EllisC. J. (2019). Climate change, bioclimatic models and the risk to lichen diversity.Diversity11:54. 10.3390/d11040054
19
FelsensteinJ. (1973). Maximum-likelihood estimation of evolutionary trees from continuous characters.Am. J. Hum. Genet.25471–492.
20
GaramszegiL. Z. (2014). “Uncertainties due to within-species variation in comparative studies: measurement errors and statistical weights,” in Modern phylogenetic Comparative Methods and Their Application in Evolutionary Biology: Concepts and Practice, ed.GaramszegiL. Z. (Berlin: Springer), 157–199. 10.1007/978-3-662-43550-2_7
21
GiordaniP.IncertiG. (2008). The influence of climate on the distribution of lichens: a case study in a borderline area (Liguria, NW Italy).Plant Ecol.195257–272. 10.1007/s11258-007-9324-7
22
HarmonL. J.LososJ. B.DaviesT. J.GillespieR. G.GittlemanJ. L.JenningsW. B.et al (2010). Early bursts of body size and shape evolution are rare in comparative data.Evolution642385–2396. 10.1111/j.1558-5646.2010.01025.x
23
HarmonL. J.SchulteJ. A.LarsonA.LososJ. B. (2003). Tempo and mode of evolutionary radiation in iguanian lizards.Science301961–964. 10.1126/science.1084786
24
HoltR. D. (1990). The microevolutionary consequences of climate change.Trends Ecol. Evol.5311–315. 10.1016/0169-5347(90)90088-U
25
HoneggerR. (1998). The lichen symbiosis—what is so spectacular about it?Lichenologist30193–212. 10.1006/lich.1998.0140
26
HoneggerR. (2009). “Lichen-forming fungi and their photobionts,” in Plant Relationships the Mycota, ed.DeisingH. B. (Berlin: Springer), 307–333. 10.1007/978-3-540-87407-2_16
27
IPCC (2014). “Climate change 2014: synthesis report,” in Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edsCore Writing TeamPachauriR. K.MeyerL. A. (Geneva: IPCC).
28
JezkovaT.WiensJ. J. (2016). Rates of change in climatic niches in plant and animal populations are much slower than projected climate change.Proc. R. Soc. B Biol. Sci.283:20162104. 10.1098/rspb.2016.2104
29
KikuchiY.TadaA.MusolinD. L.HariN.HosokawaT.FujisakiK.et al (2016). Collapse of insect gut symbiosis under simulated climate change.mBio7:e01578-16. 10.1128/mBio.01578-16
30
KoseckaM.Guzow-KrzemińskaB.ČernajováI.ŠkaloudP.JabłońskaA.KukwaM. (2021). New lineages of photobionts in Bolivian lichens expand our knowledge on habitat preferences and distribution of Asterochloris algae.Sci. Rep.11:8701. 10.1038/s41598-021-88110-0
31
LeavittS. D.KraichakE.NelsenM. P.AltermannS.DivakarP. K.AlorsD.et al (2015). Fungal specificity and selectivity for algae play a major role in determining lichen partnerships across diverse ecogeographic regions in the lichen-forming family Parmeliaceae (Ascomycota).Mol. Ecol.243779–3797. 10.1111/mec.13271
32
LeavittS. D.NelsenM. P.LumbschH. T.JohnsonL. A.ClairL. L. S. (2013). Symbiont flexibility in subalpine rock shield lichen communities in the Southwestern USA.Bryologist116149–161. 10.1639/0007-2745-116.2.149
33
LiuH.YeQ.WiensJ. J. (2020). Climatic-niche evolution follows similar rules in plants and animals.Nat. Ecol. Evol.4753–763. 10.1038/s41559-020-1158-x
34
LoppiS. (2019). May the diversity of epiphytic lichens be used in environmental forensics?Diversity11:36. 10.3390/d11030036
35
MiadlikowskaJ.KauffF.HögnabbaF.OliverJ. C.MolnárK.FrakerE.et al (2014). A multigene phylogenetic synthesis for the class Lecanoromycetes (Ascomycota): 1307 fungi representing 1139 infrageneric taxa, 317 genera and 66 families.Mol. Phylogenet. Evol.79132–168. 10.1016/j.ympev.2014.04.003
36
MünkemüllerT.LavergneS.BzeznikB.DrayS.JombartT.SchiffersK.et al (2012). How to measure and test phylogenetic signal.Methods Ecol. Evol.3743–756. 10.1111/j.2041-210X.2012.00196.x
37
MuñozM. M.StimolaM. A.AlgarA. C.ConoverA.RodriguezA. J.LandestoyM. A.et al (2014). Evolutionary stasis and lability in thermal physiology in a group of tropical lizards.Proc. R. Soc. B Biol. Sci.281:20132433. 10.1098/rspb.2013.2433
38
NelsenM. P.LumbschH. T. (2020). A data-driven evaluation of lichen climate change indicators in Central Europe.Biodivers. Conserv.293959–3971. 10.1007/s10531-020-02057-8
39
NelsenM. P.LeavittS. D.HellerK.MuggiaL.LumbschH. T. (2021). Macroecological diversification and convergence in a clade of keystone symbionts.FEMS Microbiol. Ecol.97:fiab072. 10.1093/femsec/fiab072
40
NelsenM. P.LückingR.BoyceC. K.LumbschH. T.ReeR. H. (2020). No support for the emergence of lichens prior to the evolution of vascular plants.Geobiology183–13. 10.1111/gbi.12369
41
NimisP. L.ScheideggerC.WolseleyP. (eds) (2002). Monitoring with Lichens - Monitoring Lichens.Dordrecht: Springer, 10.1007/978-94-010-0423-7
42
PeksaO.ŠkaloudP. (2011). Do photobionts influence the ecology of lichens? A case study of environmental preferences in symbiotic green alga Asterochloris (Trebouxiophyceae).Mol. Ecol.203936–3948. 10.1111/j.1365-294X.2011.05168.x
43
PennellM. W.EastmanJ. M.SlaterG. J.BrownJ. W.UyedaJ. C.FitzJohnR. G.et al (2014). geiger v2.0: an expanded suite of methods for fitting macroevolutionary models to phylogenetic trees.Bioinformatics302216–2218. 10.1093/bioinformatics/btu181
44
Pérez-OrtegaS.Ortiz-ÁlvarezR.GreenT. G. A.de los RíosA. (2012). Lichen myco- and photobiont diversity and their relationships at the edge of life (McMurdo Dry Valleys, Antarctica).FEMS Microbiol. Ecol.82429–448. 10.1111/j.1574-6941.2012.01422.x
45
PoradaP.StanJ. T. V.KleidonA. (2018). Significant contribution of non-vascular vegetation to global rainfall interception.Nat. Geosci.11563–567. 10.1038/s41561-018-0176-7
46
PoradaP.WeberB.ElbertW.PöschlU.KleidonA. (2014). Estimating impacts of lichens and bryophytes on global biogeochemical cycles.Glob. Biogeochem. Cycl.2871–85. 10.1002/2013GB004705
47
QuinteroI.WiensJ. J. (2013). Rates of projected climate change dramatically exceed past rates of climatic niche evolution among vertebrate species.Ecol. Lett.161095–1103. 10.1111/ele.12144
48
R Core Team. (2014). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
49
RamboldG.FriedlT.BeckA. (1998). Photobionts in lichens: possible indicators of phylogenetic relationships?The Bryologist101392–397. 10.2307/3244177
50
RevellL. J. (2012). phytools: an R package for phylogenetic comparative biology (and other things).Methods Ecol. Evol.3217–223. 10.1111/j.2041-210X.2011.00169.x
51
RevellL. J.HarmonL. J.CollarD. C. (2008). Phylogenetic signal, evolutionary process, and rate.Syst. Biol.57591–601. 10.1080/10635150802302427
52
RolshausenG.GrandeF. D.Sadowska-DeśA. D.OtteJ.SchmittI. (2018). Quantifying the climatic niche of symbiont partners in a lichen symbiosis indicates mutualist-mediated niche expansions.Ecography411380–1392. 10.1111/ecog.03457
53
RuprechtU.BrunauerG.PrintzenC. (2012). Genetic diversity of photobionts in Antarctic lecideoid lichens from an ecological view point.Lichenologist44661–678. 10.1017/S0024282912000291
54
RutherfordW. A.PainterT. H.FerrenbergS.BelnapJ.OkinG. S.FlaggC.et al (2017). Albedo feedbacks to future climate via climate change impacts on dryland biocrusts.Sci. Rep.7:44188. 10.1038/srep44188
55
SampayoE. M.RidgwayT.FranceschinisL.RoffG.Hoegh-GuldbergO.DoveS. (2016). Coral symbioses under prolonged environmental change: living near tolerance range limits.Sci. Rep.6:36271. 10.1038/srep36271
56
SandersW. B. (2001). Lichens: the interface between mycology and plant morphology.Bioscience511025–1035. 10.1641/0006-3568(2001)051[1025:ltibma]2.0.co;2
57
SeawardM. R. D. (2008). “Environmental role of lichens,” in Lichen Biology, ed.NashT. H. I. (Cambridge: Cambridge University Press), 274–298. 10.1017/cbo9780511790478.015
58
SlaterG. J.FrisciaA. R. (2019). Hierarchy in adaptive radiation: a case study using the Carnivora (Mammalia).Evolution73524–539. 10.1111/evo.13689
59
SlaterG. J.HarmonL. J.AlfaroM. E. (2012). Integrating fossils with molecular phylogenies improves inference of trait evolution.Evolution663931–3944. 10.1111/j.1558-5646.2012.01723.x
60
SlaterG. J.PriceS. A.SantiniF.AlfaroM. E. (2010). Diversity versus disparity and the radiation of modern cetaceans.Proc. R. Soc. B Biol. Sci.2773097–3104. 10.1098/rspb.2010.0408
61
Tschermak-WoessE. (1988). “The algal partner,” in CRC Handbook of Lichenology, ed.GalunM. (Boca Raton, FL: CRC Press), 39–92.
62
UyedaJ. C.CaetanoD. S.PennellM. W. (2015). Comparative analysis of principal components can be misleading.Syst. Biol.64677–689. 10.1093/sysbio/syv019
63
VančurováL.MuggiaL.PeksaO.ŘídkáT.ŠkaloudP. (2018). The complexity of symbiotic interactions influences the ecological amplitude of the host: a case study in Stereocaulon (lichenized Ascomycota).Mol. Ecol.273016–3033. 10.1111/mec.14764
64
VegaG. C.PertierraL. R.Olalla-TárragaM. Á (2017). MERRAclim, a high-resolution global dataset of remotely sensed bioclimatic variables for ecological modelling.Sci. Data4:170078. 10.1038/sdata.2017.78
65
VoytsekhovichA. A.MikhailyukT. I.DarienkoT. M. (2011). Lichen photobionts. 2: origin and correlation with mycobiont.Algologia21151–177.
66
WagnerM.BathkeA. C.CaryS. C.GreenT. G. A.JunkerR. R.TrutschnigW.et al (2020). Myco- and photobiont associations in crustose lichens in the McMurdo Dry Valleys (Antarctica) reveal high differentiation along an elevational gradient.Polar Biol.431967–1983. 10.1007/s00300-020-02754-8
67
WerthS. (2012). Fungal-algal interactions in Ramalina menziesii and its associated epiphytic lichen community.Lichenologist44543–560. 10.1017/S0024282912000138
68
WiensJ. J.AckerlyD. D.AllenA. P.AnackerB. L.BuckleyL. B.CornellH. V.et al (2010). Niche conservatism as an emerging principle in ecology and conservation biology.Ecol. Lett.131310–1324. 10.1111/j.1461-0248.2010.01515.x
69
Will-WolfS.JovanS. (2019). Lichen species as element bioindicators for air pollution in the eastern United States of America.Plant Fungal Syst.64137–147. 10.2478/pfs-2019-0015
70
Will-WolfS.JovanS.NeitlichP.PeckJ. E.RosentreterR. (2015). Lichen-based indices to quantify responses to climate and air pollution across northeastern U.S.A.Bryologist11859–82.
Summary
Keywords
climate, niche, diversification, lichen, photobiont, Trebouxia
Citation
Nelsen MP, Leavitt SD, Heller K, Muggia L and Lumbsch HT (2022) Contrasting Patterns of Climatic Niche Divergence in Trebouxia—A Clade of Lichen-Forming Algae. Front. Microbiol. 13:791546. doi: 10.3389/fmicb.2022.791546
Received
08 October 2021
Accepted
11 January 2022
Published
15 February 2022
Volume
13 - 2022
Edited by
Bernardo González, Adolfo Ibáñez University, Chile
Reviewed by
Victor Manuel Bandala, Instituto de Ecología (INECOL), Mexico; Ulrike Ruprecht, University of Salzburg, Austria; Rebecca Yahr, Royal Botanic Garden Edinburgh, United Kingdom
Updates
Copyright
© 2022 Nelsen, Leavitt, Heller, Muggia and Lumbsch.
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: Matthew P. Nelsen, mpnelsen@gmail.com
†ORCID: Matthew P. Nelsen, orcid.org/0000-0002-6866-815X; Steven D. Leavitt, orcid.org/0000-0002-5034-9724; Kathleen Heller, orcid.org/0000-0002-1030-1136; Lucia Muggia, orcid.org/0000-0003-0390-6169; H. Thorsten Lumbsch, orcid.org/0000-0003-1512-835X
This article was submitted to Microbe and Virus Interactions with Plants, a section of the journal Frontiers in Microbiology
Disclaimer
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.