Present-day correlations between leaf physiognomic traits (shape and size) and climate are widely used to estimate paleoclimate using fossil floras. For example, leaf-margin analysis estimates paleotemperature using the modern relation of mean annual temperature (MAT) and the site-proportion of untoothed-leaf species (NT). This uniformitarian approach should provide accurate paleoclimate reconstructions under the core assumption that leaf-trait variation principally results from adaptive environmental convergence, and because variation is thus largely independent of phylogeny it should be constant through geologic time. Although much research acknowledges and investigates possible pitfalls in paleoclimate estimation based on leaf physiognomy, the core assumption has never been explicitly tested in a phylogenetic comparative framework. Combining an extant dataset of 21 leaf traits and temperature with a phylogenetic hypothesis for 569 species-site pairs at 17 sites, we found varying amounts of non-random phylogenetic signal in all traits. Phylogenetic vs. standard regressions generally support prevailing ideas that leaf-traits are adaptively responding to temperature, but wider confidence intervals, and shifts in slope and intercept, indicate an overall reduced ability to predict climate precisely due to the non-random phylogenetic signal. Notably, the modern-day relation of proportion of untoothed taxa with mean annual temperature (NT-MAT), central in paleotemperature inference, was greatly modified and reduced, indicating that the modern correlation primarily results from biogeographic history. Importantly, some tooth traits, such as number of teeth, had similar or steeper slopes after taking phylogeny into account, suggesting that leaf teeth display a pattern of exaptive evolution in higher latitudes. This study shows that the assumption of convergence required for precise, quantitative temperature estimates using present-day leaf traits is not supported by empirical evidence, and thus we have very low confidence in previously published, numerical paleotemperature estimates. However, interpreting qualitative changes in paleotemperature remains warranted, given certain conditions such as stratigraphically closely-spaced samples with floristic continuity.
Citation: Little SA, Kembel SW, Wilf P (2010) Paleotemperature Proxies from Leaf Fossils Reinterpreted in Light of Evolutionary History. PLoS ONE 5(12): e15161. doi:10.1371/journal.pone.0015161
Editor: Andrew Allen Farke, Raymond M. Alf Museum of Paleontology, United States of America
Received: August 16, 2010; Accepted: October 27, 2010; Published: December 22, 2010
Copyright: © 2010 Little et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the David and Lucile Packard Foundation to Peter Wilf, which funded Stefan Little's position at Pennsylvania State University 2006-2009 (http://www.packard.org/home.aspx, http://www.packard.org/fellowsDetails.aspx?RootCatID=3&CategoryID=70&ItemID=2951). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
In a seminal 1915 paper, Bailey and Sinnott proposed “a botanical index of Cretaceous and Tertiary climates” : in extant mesic floras, the proportion of woody “dicot” species that have untoothed leaf margins (NT) is positively related to mean annual temperature (MAT), and thus quantifying untoothed taxa in fossil floras is informative about past temperatures. Because this relationship occurs across different continents and biomes containing various plant lineages, the authors suggested that environmental convergence is the most important explanatory factor, rather than phylogenetic causes , . Therefore, the pattern should hold through time, and paleotemperature can be inferred from fossil floras without precise systematic information, an especially useful benefit because isolated fossil leaves, the most common type of plant fossil, are notoriously difficult to identify , .
Paleobotanists have since continued to make extensive use of Bailey and Sinnott's “index,” eventually developing a quantitative method known as leaf-margin analysis, based on linear regressions of extant proportions of untoothed species (NT) and mean annual temperature (MAT) –. More recently, multivariate methods to quantify temperature and other climate variables have been proposed that also use NT and additional site-based means of numerous leaf-physiognomic (size and shape) variables; crucially, these methods still derive most of their predictive power for temperature from the NT-MAT correlation , –. To date, hundreds of papers have been produced regarding leaf physiognomy, climate, and paleoclimate, following the ideas of Bailey and Sinnott in a “taxon-free” approach, i.e., without phylogenetic considerations ; many of these are compiled in an extended topical bibliography in File S1. The general, often-stated consensus (but see below) is that the climatic distribution of leaf physiognomic traits should be similar in the past and that phylogeny is a negligible component.
In particular, precise, quantitative paleotemperature estimates from taxon-free approaches operate under a still-untested core assumption that leaf-trait variation principally results from adaptive environmental convergence, and because it is thus largely independent of phylogeny it should be constant through geologic time. If the assumption is valid, then the current uniformitarian applications are warranted, and modern trait-climate relations should estimate past climate in a quantitatively precise manner. However, certain observations cast doubt on the core assumption , –. First, leaf-climate correlations, including that of NT-MAT, vary sometimes considerably across biogeographic regions, suggesting an historical influence on trait-climate relationships , , , –. Second, many clades have obligate traits independent of the environment, including being typically toothed or untoothed (e.g., Betulaceae, Lauraceae, Myrtaceae, many Nothofagus, and Rosaceae), indicating phylogenetic signal , , . Third, the relative richness of species with particular traits is affected by factors other than temperature, including differential origination or extinction among clades , , , , suggesting possible inconstancy of trait-climate relationships through time. Fourth, adaptations to climate change can involve many aspects of plant biology, including anatomy, physiology and biochemistry, and thus leaf traits are not necessarily expected to respond strongly to climate, depending on the nature of correlations between leaf traits and other aspects of the phenotype –. Despite these considerations raised in the literature, the issues above are usually considered to have insignificant effects when reconstructing paleoclimates.
An additional, more theoretical consideration is that even under an ideal model of adaptive convergence without phylogenetic signal, directional selection is expected to affect the assumption of constancy through time. Given a hypothetically unchanging climate, with constant species composition, directional selection should shift trait values over time; producing an increase in number of teeth over geological time under a constant low MAT for example. Thus, lineage persistence may contribute to some observed changes in leaf-traits through geologic time, unrelated to climate change. Further, given that lineages are expected to have differing rates of evolution, extinction and species radiation, trait-climate relations should vary through deep time, even if adaptive leaf-trait responses were ideal , .
The issues above are important because if present-day trait-climate relations do not have the expected adaptive explanation, or were not constant over evolutionary time, then the core assumption of the current uniformitarian approach would not be valid (i.e., leaf physiognomy-climate relations may have differed in the past). In this case, one would require additional information in order to precisely quantify a past trait-climate relation with confidence, such as data regarding phylogenetic placement of difficult-to-identify fossil leaves, or independent, highly accurate and well-correlated climate proxies. We note that several examples of temperature estimates from isotopic data are considered broadly concordant with associated leaf-physiognomy estimates –; however, these confirm qualitative temperature changes (presence of warming and cooling) much more robustly than they support precise quantitative estimates from leaf physiognomic data. The development of phylogenetic comparative methods provides a novel approach to investigate calibration floras used in paleoclimate estimates. Here, we use phylogenetic comparative methods to provide a more explicit investigation into the influence of phylogenetic history in modern-day leaf-trait variation with temperature.
We test the assumption of adaptive convergence of leaf traits with temperature by quantifying phylogenetic signal in a dataset comprising 21 leaf-physiognomic traits and MAT among 569 species-site pairs at 17 sites in the eastern USA and Barro Colorado Island, Republic of Panamá, published by Huff et al.  and Royer et al.. The dataset shows strong correlations between leaf-physiognomic traits and MAT , and is ideal for a conservative test of phylogenetic effects: it is restricted to a region widely used to calibrate “leaf paleothermometers,” and it thus minimizes biogeographic effects on leaf parameters to a scale where they usually are considered unimportant , , . In addition, these tests are conservative with regard to the effect of phylogenetic signal because the phylogenetic hypothesis used is resolved only to the family level (see Materials and Methods). Further, to determine the effects of evolutionary history on adaptive interpretations of leaf physiognomy–temperature relationships, we compared non-phylogenetic and phylogenetic model parameters and fits of trait-temperature relationships.
We note that the issues addressed by our analysis are separate from the numerous additional sources of uncertainty in leaf paleothermometry that have been noted, including environmental, taphonomic/preservational, sampling, and scoring biases , , , , –.
The Huff et al. and Royer et al. dataset ,  comprises 21 leaf traits of woody “dicot” angiosperms representing 569 species-site pairs at 17 sites in the eastern USA and Barro Colorado Island, Republic of Panamá, over a range in MAT of 5.6–25.8°C. Of the 21 leaf-traits, 10 are perimeter- and area-derived measurements and 11 are tooth measurements, or ratios of tooth measurements with perimeter- and/or area-derived measurements , . The Panamá site has the highest mean annual temperature and species richness, which includes lineages not found at the other sites, so its effects are also addressed (see below). We compared model parameters and fits of trait-MAT relationships between non-phylogenetic and phylogenetic generalized least squares (GLS and pGLS; ) to determine the effects of evolutionary history on the prevailing adaptive interpretations of leaf physiognomy–MAT relationships. Models were compared using the Akaike Information Criterion (AIC), which evaluates the fit of a model taking into account differences in the number of parameters included in the model . Models with lower AIC scores are considered more parsimonious and more strongly supported by the data.
We created a phylogenetic hypothesis for the species included in this study by grafting them onto a family-level phylogenetic supertree of the angiosperms  using the Phylomatic toolkit . Genera were placed as polytomies within families, and species found at several sites were placed as polytomies at the species level; each species at each site was placed as a tip of the phylogeny. Branch lengths in the resulting phylogeny represent estimates of clade age based on spacing undated nodes evenly between dated nodes in the original supertree to produce an ultrametric tree . The resulting phylogeny maintained the resolution of phylogenetic relationships at the family/genus level and below, providing a conservative data set to test for the influence of phylogenetic history on traits and relationships among traits.
Phylogenetic signal is a tendency for closely related taxa to possess similar trait values due to descent from a common ancestor. Phylogenetic signal in all traits, including MAT, was measured using the K statistic, which compares the observed phylogenetic signal in a trait to the signal expected under a Brownian motion model of evolution , . Higher values of K indicate stronger phylogenetic signal (greater similarity of closely related species), with K = 1 expected for traits evolving under a Brownian motion model of evolution . The statistical significance (P-values) of the K values for each trait was assessed by comparing the variance of independent contrasts for each trait to the expected values under a tip shuffling algorithm . This P-value provides a test of whether the phylogenetic signal in each trait is greater than the null expectation of no signal, while the K statistic provides an estimate of the magnitude of phylogenetic signal. All measures of K were calculated using the picante R package . The Panamá site's effects are addressed by measuring K with this site removed (Table 1).
Non-phylogenetic generalized least squares (GLS) and phylogenetic GLS (pGLS) regressions between climate and leaf traits were computed using the APE  and nlme  R packages. All leaf and climate traits except ‘Margin untoothed’ (see below) were log10-transformed prior to analysis in order to meet model assumptions of normality. For the phylogenetic GLS analyses, branch lengths were first scaled using the optimal value of Pagel's λ parameter  as calculated using the GEIGER R package . The explanatory power of each model was evaluated using the Akaike Information Criterion (AIC) , which allows comparison of models with different numbers of estimated parameters. For each trait, the AIC scores of phylogenetic and non-phylogenetic GLS models were compared to determine whether the inclusion of phylogenetic information improved model fits, with lower AIC score indicating a better fit. We note that phylogenetic comparative methods cannot accommodate leaf-physiognomic means by site, as typically used in leaf-paleoclimate estimates, and thus we examined MAT correlations at the species-level.
We performed analyses treating the ‘Margin untoothed’ trait as a binary and as a ternary variable to perform comparisons among the phylogenetic vs. non phylogenetic trait-temperature GLS models. Current physiognomic methods use the proportion of untoothed taxa at a site as a continuous variable in a least squares regression with site-climate variables. Because our analyses were based on species-level data and not site means, ‘Margin untoothed’ was defined as a ternary variable for each species at a site as follows: 1 = all leaves untoothed, 0.5 = both toothed and untoothed, 0 = toothed. Leaf margin analysis and multivariate physiognomic methods treat the presence-absence of leaf teeth in this way to calculate the site means that are used as continuous values , , . Because this variable is not truly continuous, we also performed GLS regressions treating ‘Margin untoothed’ as a binomial variable using the glmmPQL function in the MASS package  for the R statistical language and computing environment . For this analysis, toothedness for each species at a site was coded as 1 (all specimens untoothed) or 0 (any specimens toothed). The binomial GLS regression found similar patterns of differences in slope and parameter uncertainty in phylogenetic versus nonphylogenetic GLS models, concordant with the non-binomial GLS analyses presented here (Table 1), where slope is greatly reduced after accounting for phylogeny.
There was non-random phylogenetic signal in MAT and all measured leaf traits (Table 1: all K, P-values ≤0.001; Figure 1). Although all traits exhibited non-random phylogenetic signal, the amount of signal varied among traits. Tooth traits (e.g., number of teeth, average tooth area) exhibited the strongest phylogenetic signal (highest K values, e.g. average tooth area; K = 1.9). There was also phylogenetic signal in the climatic associations of taxa (K-value for MAT = 0.6, P = 0.001). Phylogenetic signal in traits and climate associations was not driven by the inclusion of the single tropical site (Barro Colorado Island, Republic of Panamá); when the tropical site was excluded (413 species-site pairs), signal remained significantly non-random for all traits and MAT (Table 1).
Mean tooth-area character mapped at tips of the phylogeny of all species in the community samples  to illustrate non-random variation across the phylogeny of both presence of teeth and tooth area.
All trait-climate regression models were improved by incorporation of phylogenetic relationships (Table 1; Figure 2), as demonstrated by the much lower AIC scores for phylogenetic GLS versus non-phylogenetic GLS models for all traits. The standard errors of the intercepts in pGLS regressions were generally higher than in non-phylogenetic GLS regressions (Table 1), leading to greater uncertainty in predictions of climate after phylogenetic relationships are taken into account.
Mean annual temperature (MAT) versus the six leaf traits featured in Figure 2 of Royer et al. . For each trait, the best-fit lines for nonphylogenetic GLS (dashed line) and phylogenetic GLS (solid line) are displayed. All traits except Margin untoothed and MAT were log10-transformed. Phylogenetic models used branch lengths scaled by the best-fit estimate of Pagel's λ parameter  (Table 1). A 95% confidence interval is displayed for each regression model (dashed rose lines = nonphylogenetic GLS, solid blue lines = phylogenetic GLS) to illustrate the increased uncertainty in predictions of climate from leaf traits when phylogeny is taken into account. Points were lightly jittered at each site to better visualize density of trait values.
Traits with the strongest phylogenetic signal (highest K values; Margin untoothed, Feret diameter ratio, Tooth area : Blade area, and Average tooth area) also showed the largest decreases in slope versus MAT after taking phylogeny into account (GLS vs. pGLS slopes; Figure 2, Table 1). Several trait-MAT relations with negative slopes (Perimeter, Internal Perimeter, Minor axis length, and Number of secondary teeth) or weakly positive slopes (Blade Area, Major Axis Length, and Feret Diameter) in nonphylogenetic models had steeper negative slope in pGLS. In general, tooth-related traits showed both strong phylogenetic signal and relatively large shifts in slopes between GLS and pGLS models. Conversely, non-tooth related traits such as Perimeter ratio, Compactness, and Shape Factor (a modified area:perimeter ratio), had the least altered regressions after accounting for phylogenetic relationships, in accord with their relatively low phylogenetic signal (Table 1).
The presence of non-random phylogenetic signal in all traits (Table 1) invalidates the core assumption that convergence dominates, and that phylogenetic history is an insignificant component of modern-day leaf-trait variation with climate. This assumption is required for “taxon-free”, quantitative leaf paleothermometry as currently used. The finding that tooth traits have the greatest amount of phylogenetic signal (highest K values) reinforces earlier observations that characters of leaf teeth are taxonomically informative . The presence of moderate phylogenetic signal in MAT is compatible with indications of biome conservatism in plant clades through deep time  and supports some of the ideas behind taxon-based paleoclimate methods , –. When the tropical site was excluded, signal remained significant (Table 1), showing that phylogenetic signal among leaf traits and climate associations exists even within Eastern North America, a biogeographic region with relatively minor floristic variation that is often used for leaf-physiognomic calibrations , , .
The greatly improved model fits (reduced AIC values) for all pGLS regressions indicate that incorporation of phylogenetic information is important for understanding trait-temperature relationships (Table 1; Figure 2). Incorporating phylogenetic information into trait-climate models changed both the estimates of a relationship itself (slopes and intercepts) as well as estimates of the certainty of predictions (standard errors of slopes and intercepts). Uncertainty in prediction (standard error of intercept estimates) was increased after accounting for phylogenetic relationships, indicating that previous studies have overestimated the ability to precisely predict paleoclimate via modern-day plant traits.
The nonzero slopes of phylogenetic regressions (pGLS models; Table 1; Figure 2) indicate that many leaf-physiognomic traits show varying degrees of adaptive responses to temperature (as well as significant and varying phylogenetic signal as previously discussed). Importantly, this includes tooth traits other than tooth presence (e.g., number of teeth, mean tooth area), consistent with species-level observations wherein tooth traits covary with temperature across the geographic distribution of species , . Our results are also consistent with some leaf trait responses to temperature in a single-generation experiment from seed  where Acer rubrum from a Florida population had more leaf teeth when planted in Rhode Island but still had significantly fewer teeth than the Rhode Island population, supporting the presence of both phylogenetic signal and phenotypic plasticity, as well as probable adaptive trait response in the species. Non-tooth related traits, especially Shape factor (a modified area:perimeter ratio), had the least altered correlations, in accord with their relatively low phylogenetic signal (K; Table 1). Recently developed tools to measure perimeter-based leaf traits show utility in detecting subtle phenotype variation among mutant plant lines –, which suggests that other, more sensitive, perimeter- derived measurements ,  could show differing values for phylogenetic signal than shown in this study.
Overall, our results support the prevailing idea that leaves are adaptively responding to climate, but that phylogenetic signal in leaf traits is responsible for a portion of variation in leaf-climate relationships, and phylogenetic information modifies our understanding of adaptive relationships between leaf physiognomic variables and climate. Several insights for improved understanding of adaptive relationships between leaf traits and temperature were revealed by our analyses; for example, Feret diameter did not appear to be strongly responsive to temperature in the dataset previously (GLS slope) , but after accounting for phylogenetic relationships there is a stronger relationship with temperature (i.e. steeper slope in pGLS vs. GLS; Table 1; Figure 2). This suggests that the influence of phylogenetic history can possibly mask leaf-trait climate relations that are not discernible in standard correlations that do not incorporate phylogenetic information.
Although we expected to find some influence of phylogenetic signal in the leaf traits and their relation to temperature, the most surprising results involved the high signal in presence of leaf teeth (Margin Untoothed). This trait is the basis of leaf-margin analysis (NT-MAT relation) and a key component of all multivariate approaches that estimate paleotemperature , –, . After accounting for phylogenetic relationships among species, there is an extremely weak relationship between presence of leaf teeth and temperature (NT-MAT pGLS; Figure 2), including an order of magnitude (in log space) decrease in slope, approaching a flat line, and substantially wider confidence intervals (Figure 2; Table 1). This pattern was observed whether presence of leaf teeth was treated as a binomial or ternary variable (Table 1). Because the NT-MAT correlation is fundamental for both univariate and multivariate paleotemperature inference as widely practiced, this result impacts most paleotemperature proxies using leaf physiognomy, as well as prevailing adaptive interpretations of leaf-teeth with temperature [7,8,10,12, S1–S351].
The relations between MAT and the traits Feret diameter ratio, Tooth area : blade area, and Average tooth area also had phylogenetic regressions with highly flattened slopes, similar to that of the NT-MAT relation (Figure 2; Table 1). Importantly, from this same data-set, Feret diameter ratio, and Tooth area : blade area were considered among the most useful adaptive traits in multivariate approaches to temperature inference based on standard trait-climate correlations .
The weak relationship between presence of leaf teeth and temperature after accounting for phylogeny indicates that the prevailing adaptive scenario since 1915 , whereby temperature is the primary force explaining evolutionary gain or loss of leaf teeth, is inaccurate. The small, but non-zero, slope of the phylogenetic regression suggests that temperature is at most only weakly related to the evolutionary gain and loss of teeth, and thus the proportion of toothed species in a flora at a given temperature would not be expected to be constant through time. The strong phylogenetic signal in leaf-tooth presence we observed is consistent with our observations of Southern Hemisphere data that clearly show phylogenetic conservatism in tooth presence or absence, including regional studies of Chilean and Australian forests ,  and among the lineages compiled in a large whole-Hemisphere dataset . Interestingly, as stated above, temperature does appear to have a large influence on the evolution of tooth traits other than presence of teeth (Number of teeth, Tooth area, etc; Figures 1, 2; Table 1).
Because our results indicate that historical events unrelated to temperature contributed to the majority of the present-day distribution of toothed lineages, the modern NT-MAT relationship, which is variable at a global scale, now deserves renewed investigation , . Because both temperature and inheritance influence traits of leaf teeth (number of teeth, tooth area, etc.; Figure 2; Table 1), we suggest that a component of the modern-day distribution of untoothed species may be explained as an exaptive scenario where cool-temperature selection acted on preexisting toothed lineages (Table 1; Figure 2) . Temperature selection on already-toothed lineages may have affected diversification and thereby increased the proportional richness of toothed taxa relative to warmer places and times. Thus, under this scenario, the modern-day NT-MAT correlation is primarily a result of biogeographic history, and secondarily altered by effects from exaptive temperature selection on toothed lineages.
The cool-temperature selection scenario, suggested here, is also provisionally consistent with the fossil record, in which many characteristically, or commonly, toothed clades (e.g. Betulaceae, temperate Nothofagus lineages, Ulmaceae, Rosaceae, Vitaceae) radiated at middle and high latitudes under warmer climates than today , –. The groups remained and further speciated at these latitudes through Cenozoic global cooling beginning in the Eocene. In this scenario, the incumbent, already toothed lineages would have been exapted and shown subsequent adaptive changes in traits of teeth (not their presence-absence) during cooling, hypothetically influencing speciation rates over time.
Observations of changing relative presence of leaf teeth through deep time that are qualitatively validated by correlation to independent temperature proxies, are often considered to be related to the migrations of clades along temperature gradients (e.g., latitude, altitude) –, –. However, the same studies typically invoke convergent adaptive response to justify the application of quantitative paleotemperature estimates from the same fossil floras. Our results indicate that clade migrations along temperature gradients were probably the principal causes of the observed changes in the proportion of species with leaf teeth in these fossil floras. Inference of paleoclimate from relative proportions of toothed species will only be accurate to the extent that the distribution of toothed species and the patterns of biogeographic migration along paleotemperature gradients were the same in the past as today. As we have shown, this assumption is not warranted because of the weak adaptive relationship between relative presence of teeth and temperature, after accounting for phylogeny, which greatly decreases confidence in estimates of paleoclimate.
Qualitative analyses, i.e., detection of relative warming and cooling, remain justified using physiognomic data from well-understood regional floras that have supporting data on taxonomy, paleogeography, and distribution of traits along independently inferred paleotemperature gradients (i.e. using floras from several latitudinally adjacent basins). In practical terms, these conditions are met for several heavily studied assemblages (i.e., latest Cretaceous and Paleogene floras of the Western USA and Germany). However, we reiterate that although physiognomic trait-latitude gradients clearly existed at many times and places in the past with varying similarity to the modern day , , , the explicit trait-temperature gradients are usually unknown in the past. Thus, although past trait-latitude gradients allow for qualitative climate inference, confidence is very low for quantitative models of past trait-temperature gradients based directly from modern gradients as currently practiced. This point is well-demonstrated by discrepancies between taxonomic and physiognomic paleoclimate estimates, even when derived by the same workers from the same fossil floras , . Explicit reconstruction of past gradients could potentially surmount this problem but would require major advances in independent paleotemperature proxies and their precise correlation to fossil floras.
Despite our results, it may be tempting to continue relying on current, taxon-free leaf physiognomy to generate quantitatively inferred paleotemperature estimates, by using traits that show clear adaptive responses to temperature (i.e. Number of teeth), using traits that display only slightly altered phylogenetic regressions with temperature (i.e., Shape factor), using standard leaf-margin analysis as a convenient proxy for the presence of adaptive tooth-trait response, or relying on multivariate approaches , ,  in the hope that using many traits, or site averages of traits, will reduce prediction error. In light of our results, we advise against the above strategies because all leaf traits show non-random phylogenetic signal (Figure 1; Table 1). The use of site-averages, single-variable proxies, or multiple traits will not remove this underlying signal or reconcile the inherent uncertainties from non-independent data, and would mask the true uncertainty in prediction error. If phylogenetic signal in leaf traits is a general global phenomenon, as we strongly expect, then broader sampling of present-day “calibration” floras  would have the counterintuitive effect of increasing the influence of phylogenetic signal in leaf physiognomic data as more regions of the angiosperm tree of life would be sampled.
Regional differences in the relationship between proportions of toothed species and temperature (NT-MAT), wherein temperature estimates for a given value of NT differ by >5°C, may well be due to differences in phylogenetic history among biogeographic regions , , , –, , , and may provide a rough empirical approximation of uncertainty due to phylogenetic history. However, the amount of extant variation may not apply to the past. For example, a large proportion of early-diverging angiosperm lineages are toothed  and tropical , suggesting possible shifting correlations through time, consistent with observations of the earliest angiosperm leaf assemblages having atypical latitudinal distributions of toothed leaves .
In summary, we have demonstrated that there is evidence for an adaptive response to temperature in many leaf traits. However, the presence of non-random phylogenetic signal throughout leaf physiognomic data leads to leaf trait-climate relationships that are driven both by adaptive evolution and phylogenetic history, and the adaptive signal is especially weak for the most widely used variable, presence of teeth. Non-independence of species data due to phylogenetic relatedness results in conventional, non-phylogenetic models of leaf trait-climate relationships underestimating the true uncertainty in estimates of paleotemperature from leaf traits. An approach that should permit reliable qualitative estimation of change in paleoclimate from leaf traits would be to use leaf physiognomic variables that show the strongest evolutionary correlations with climate based on pGLS models. These variables should be able to detect relative temperature changes through time, especially over geologically short intervals with floristic continuity, given sufficient evidence that the trait-climate relationships were not greatly altered in the past.
We thank D. Royer for significant assistance as well as M. Carvalho, G. Jordan, and D. Royer for helpful critiques of earlier drafts.
Conceived and designed the experiments: SAL SWK PW. Analyzed the data: SWK. Contributed reagents/materials/analysis tools: SAL. Wrote the paper: SAL SWK PW. Compiled File S1: PW.
- 1. Bailey IW, Sinnott EW (1915) A botanical index of Cretaceous and Tertiary climates. Science 41: 831–834.
- 2. Bailey IW, Sinnott EW (1916) The climatic distribution of certain types of angiosperm leaves. Am J Bot 3: 24–39.
- 3. Dilcher DL (1973) A paleoclimatic interpretation of the Eocene floras of southeastern North America. In: Graham A, editor. Vegetation and vegetational history of Northern Latin America. Amsterdam: Elsevier. pp. 39–59.
- 4. Ellis B, Daly DC, Hickey LJ, Johnson KR, Mitchell JD, et al. (2009) Manual of leaf architecture. Ithaca: Cornell University Press.
- 5. Wolfe JA (1971) Tertiary climatic fluctuations and methods of analysis of Tertiary floras. Palaeogeogr Palaeoclimatol Palaeoecol 9: 27–57.
- 6. Wolfe JA (1979) Temperature parameters of humid to mesic forests of Eastern Asia and relation to forests of other regions of the Northern Hemisphere and Australasia. U S Geol Surv Prof Pap 1106: 1–37.
- 7. Wing SL, Greenwood DR (1993) Fossils and fossil climate: the case for equable continental interiors in the Eocene. Philos Trans R Soc London Ser B 341: 243–252.
- 8. Wilf P (1997) When are leaves good thermometers? A new case for Leaf Margin Analysis. Paleobiology 23: 373–390.
- 9. Greenwood DR, Wilf P, Wing SL, Christophel DC (2004) Paleotemperature estimation using leaf-margin analysis: is Australia different? Palaios 19: 129–142.
- 10. Wolfe JA (1993) A method of obtaining climatic parameters from leaf assemblages. U S Geol Surv Bull 2040: 1–71.
- 11. Wolfe JA (1995) Paleoclimatic estimates from Tertiary leaf assemblages. Annu Rev Earth Plan Sci 23: 119–142.
- 12. Royer DL, Wilf P, Janesko DA, Kowalski EA, Dilcher DL (2005) Correlations of climate and plant ecology to leaf size and shape: potential proxies for the fossil record. Am J Bot 92: 1141–1151.
- 13. Green WA (2006) Loosening the CLAMP: an exploratory graphical approach to the Climate Leaf Analysis Multivariate Program. Palaeontol Electron 9: art. 9A.
- 14. Spicer RA, Valdes PJ, Spicer TEV, Craggs HJ, Srivastava G, et al. (2009) New developments in CLAMP: Calibration using global gridded meteorological data. Palaeogeogr Palaeoclimatol Palaeoecol 283: 91–98.
- 15. Greenwood DR (2007) Fossil angiosperm leaves and climate: from Wolfe and Dilcher to Burnham and Wilf. Cour Forschungsinst Senckenb 258: 95–108.
- 16. Jordan GJ (1997) Uncertainty in paleoclimatic reconstructions based on leaf physiognomy. Aust J Bot 45: 527–547.
- 17. Mosbrugger V, Utescher T (1997) The coexistence approach -— a method for quantitative reconstructions of Tertiary terrestrial palaeoclimate data using plant fossils. Palaeogeogr Palaeoclimatol Palaeoecol 134: 61–86.
- 18. Utescher T, Mosbrugger V, Ashraf AR (2000) Terrestrial climate evolution in northwest Germany over the last 25 million years. Palaios 15: 430–449.
- 19. Hinojosa LF, Pérez F, Gaxiola A, Sandoval I (2010) Historical and phylogenetic constraints on the incidence of entire leaf margins: insights from a new South American model. Global Ecol Biogeogr. DOI:10.1111/j.1466-8238.2010.00595.x.
- 20. Wolfe JA (1987) Late Cretaceous-Cenozoic history of deciduousness and the terminal Cretaceous event. Paleobiology 13: 215–226.
- 21. Burnham RJ (1997) Stand characteristics and leaf litter composition of a dry forest hectare in Santa Rosa National Park, Costa Rica. Biotropica 29: 384–395.
- 22. Stranks L, England P (1997) The use of a resemblance function in the measurement of climatic parameters from the physiognomy of woody dicotyledons. Palaeogeogr Palaeoclimatol Palaeoecol 131: 15–28.
- 23. Kowalski EA (2002) Mean annual temperature estimation based on leaf morphology: a test from tropical South America. Palaeogeogr Palaeoclimatol Palaeoecol 188: 141–165.
- 24. Hinojosa LF, Armesto JJ, Villagrán C (2006) Are Chilean coastal forests pre-Pleistocene relicts? Evidence from foliar physiognomy, palaeoclimate, and phytogeography. J Biogeogr 33: 331–341.
- 25. Aizen MA, Ezcurra C (2008) Do leaf margins of the temperate forest flora of southern South America reflect a warmer past? Global Ecol Biogeogr 17: 164–174.
- 26. Wolfe JA (1966) Tertiary plants from the Cook Inlet Region, Alaska. U S Geol Surv Prof Pap 398-B: 1–32.
- 27. Manchester SR, Dilcher DL, Wing SL (1998) Attached leaves and fruits of myrtaceous affinity from the middle Eocene of Colorado. Rev Palaeobot Palynol 102: 153–163.
- 28. Wolfe JA, Upchurch GR Jr (1986) Vegetation, climatic and floral changes at the Cretaceous-Tertiary boundary. Nature 324: 148–152.
- 29. Wright IJ, Ackerly DD, Bongers F, Harms KE, Ibarra-Manríquez G, et al. (2007) Relationships among ecologically important dimensions of plant trait variation in seven Neotropical forests. Annals of Botany 99: 1003–1015.
- 30. Poorter L, Wright SJ, Paz H, Ackerly DD, Condit R, et al. (2008) Are functional traits good predictors of demographic rates? Evidence from five Neotropical forests. Ecology 89: 1908–1920.
- 31. Chave J, Coomes D, Jansen S, Lewis SL, Swenson NG, et al. (2009) Towards a worldwide wood economics spectrum. Ecol Lett 12: 351–366.
- 32. Albarán-Lara AL, Mendoza-Cuenca L, Valencia-Avalos S, González-Rodríguez A, Oyama K (2010) Leaf fluctuating asymmetry increases with hybridization and introgression between Quercus magnoliifolia and Quercus resinosa (Fagaceae) through an altitudinal gradient in Mexico. Int J Plant Sci 171: 310–322.
- 33. Wolfe JA, Upchurch GR Jr (1987) North American nonmarine climates and vegetation during the Late Cretaceous. Palaeogeogr Palaeoclimatol Palaeoecol 61: 33–77.
- 34. Wolfe JA, Poore RZ (1982) Tertiary marine and nonmarine climatic trends. In: Berger W, Crowell JC, editors. Climate in Earth history. Washington, DC: National Academy of Sciences. pp. 154–15.
- 35. Wing SL, Bao H, Koch PL (2000) An early Eocene cool period? Evidence for continental cooling during the warmest part of the Cenozoic. In: Huber BT, MacLeod K, Wing SL, editors. Warm climates in earth history. Cambridge: Cambridge University Press. pp. 197–237.
- 36. Wilf P, Johnson KR, Huber BT (2003) Correlated terrestrial and marine evidence for global climate changes before mass extinction at the Cretaceous-Paleogene boundary. Proc Natl Acad Sci U S A 100: 599–604.
- 37. Fricke HC, Wing SL (2004) Oxygen isotope and paleobotanical estimates of temperature and δ18O-latitude gradients over North America during the early Eocene. Am J Sci 304: 612–635.
- 38. Huff PM, Wilf P, Azumah EJ (2003) Digital future for paleoclimate estimation from fossil leaves? Preliminary results. Palaios 18: 266–274.
- 39. Upchurch GR Jr, Wolfe JA (1993) Cretaceous vegetation of the Western Interior and adjacent regions of North America. Geol Assoc Can Spec Pap 39: 243–281.
- 40. Burnham RJ, Pitman NCA, Johnson KR, Wilf P (2001) Habitat-related error in estimating temperatures from leaf margins in a humid tropical forest. Am J Bot 88: 1096–1102.
- 41. Peppe DJ, Royer DL, Wilf P, Kowalski EA (2010) Quantification of large uncertainties in fossil leaf paleoaltimetry. Tectonics 29: TC3015. doi:10.1029/2009TC002549.
- 42. Grafen A (1989) The phylogenetic regression. Philos Trans R Soc London Ser B 326: 119–157.
- 43. Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information-theoretic approach. New York, NY, USA: Springer-verlag. 488 p.
- 44. Davies TJ, Barraclough TG, Chase MW, Soltis PS, Soltis DE, et al. (2004) Darwin's abominable mystery: insights from a supertree of the angiosperms. Proc Natl Acad Sci U S A 101: 1904–1909.
- 45. Webb CO, Donoghue MJ (2005) Phylomatic: tree assembly for applied phylogenetics. Mol Ecol Notes 5: 181–183.
- 46. Blomberg SP, Garland T Jr (2002) Tempo and mode in evolution: phylogenetic inertia, adaptation and comparative methods. J Evol Biol 15: 899–910.
- 47. Blomberg SP, Garland T Jr, Ives AR (2003) Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57: 717–745.
- 48. Kembel SW, Cowan PD, Helmus M, Cornwell WK, Morlon H, et al. (2010) Picante: R tools for integrating phylogenies and ecology. Bioinformatics 26: 1463–1464.
- 49. Paradis E, Claude J, Strimmer K (2004) APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20: 289–290.
- 50. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team (2009) nlme: linear and nonlinear mixed effects models. R package version 3: 1–92.
- 51. Pagel M (1999) Inferring the historical patterns of biological evolution. Nature 401: 877–884.
- 52. Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W (2008) GEIGER: investigating evolutionary radiations. Bioinformatics 24: 129–131.
- 53. Venables WN, Ripley BD (2002) Modern applied statistics with S. New York: Springer-Verlag.
- 54. R Development Core Team (2009) R: a language and environment for statistical computing. Software version 2.10.0. http://www.R-project.org: R Foundation for Statistical Computing.
- 55. Hickey LJ, Wolfe JA (1975) The bases of angiosperm phylogeny: vegetative morphology. Ann Mo Bot Gard 62: 538–589.
- 56. Crisp MD, Arroyo MTK, Cook LG, Gandolfo MA, Jordan GJ, et al. (2009) Phylogenetic biome conservatism on a global scale. Nature 458: 754–756.
- 57. Chaloner WG, Creeber GT (1990) Do fossil plants give a climatic signal? J Geol Soc London 147: 343–350.
- 58. Greenwood DR, Wing SL (1995) Eocene continental climates and latitudinal temperature gradients. Geology 23: 1044–1048.
- 59. McIver EE, Basinger JF (1999) Early Tertiary floral evolution in the Canadian high arctic. Ann Mo Bot Gard 86: 523–545.
- 60. Poole I, Cantrill D, Utescher T (2005) A multi-proxy approach to determine Antarctic terrestrial palaeoclimate during the Late Cretaceous and Early Tertiary. Palaeogeogr Palaeoclimatol Palaeoecol 222: 95–121.
- 61. Adams JM, Green WA, Zhang Y (2008) Leaf margins and temperature in the North American flora: recalibrating the paleoclimatic thermometer. Global Planet Change 60: 523–534.
- 62. Royer DL, McElwain JC, Adams JM, Wilf P (2008) Sensitivity of leaf size and shape to climate within Acer rubrum and Quercus kelloggii. New Phytol 179: 808–817.
- 63. Royer DL, Meyerson LA, Robertson KM, Adams JM (2009) Phenotypic plasticity of leaf shape along a temperature gradient in Acer rubrum. PLoS One 4: e7653.
- 64. Bensmihen S, Hanna AI, Langlade NB, Micol JL, Bangham A, et al. (2008) Mutational spaces for leaf shape and size. Hfsp Journal 2: 110–120.
- 65. Feng XZ, Wilson Y, Bowers J, Kennaway R, Bangham A, et al. (2009) Evolution of allometry in Antirrhinum. Plant Cell 21: 2999–3007.
- 66. Casci T (2010) Population genetics breaking down hybrids. Nature Reviews Genetics 11: 594.
- 67. Backhaus A, Kuwabara A, Bauch M, Monk N, Sanguinetti G, et al. (2010) LEAFPROCESSOR: a new leaf phenotyping tool using contour bending energy and shape cluster analysis. New Phytol 187: 251–261.
- 68. Klingenberg CP (2010) Evolution and development of shape: integrating quantitative approaches. Nature Reviews Genetics 11: 623–635.
- 69. Krieger JD, Guralnick RP, Smith DM (2007) Generating empirically determined, continuous measures of leaf shape for paleoclimate reconstruction. Palaios 22: 212–219.
- 70. Bylesjö M, Segura V, Soolanayakanahally RY, Rae AM, Trygg J, et al. (2008) LAMINA: a tool for rapid quantification of leaf size and shape parameters. BMC Plant Biol 8: doi:10.1186/1471-2229-8-82.
- 71. Greenwood DR (2005) Leaf form and the reconstruction of past climates. New Phytol 166: 355–357.
- 72. Royer DL, Kooyman RM, Little SA, Wilf P (2009) Ecology of leaf teeth: a multi-site analysis from eastern Australian subtropical rainforest. Am J Bot 96: 738–750.
- 73. Feild TS, Sage TL, Czerniak C, Iles WJD (2005) Hydathodal leaf teeth of Chloranthus japonicus (Chloranthaceae) prevent guttation-induced flooding of the mesophyll. Plant Cell Environ 28: 1179–1190.
- 74. Royer DL, Wilf P (2006) Why do toothed leaves correlate with cold climates? Gas exchange at leaf margins provides new insights into a classic paleotemperature proxy. Int J Plant Sci 167: 11–18.
- 75. Manchester SR (1999) Biogeographical relationships of North American Tertiary floras. Ann Mo Bot Gard 86: 472–522.
- 76. Hunt RJ, Poole I (2003) Paleogene West Antarctic climate and vegetation history in light of new data from King George Island. Geol Soc Am Spec Pap 369: 395–412.
- 77. Wing SL (1987) Eocene and Oligocene floras and vegetation of the Rocky Mountains. Ann Mo Bot Gard 74: 748–784.
- 78. Wilf P (2000) Late Paleocene-early Eocene climate changes in southwestern Wyoming: paleobotanical analysis. Geol Soc Am Bull 112: 292–307.
- 79. Currano ED, Labandeira CC, Wilf P (2010) Fossil insect folivory tracks paleotemperature for six million years. Ecol Monogr 80: 547–567.
- 80. Wing SL, Harrington GJ, Smith FA, Bloch JI, Boyer DM, et al. (2005) Transient floral change and rapid global warming at the Paleocene-Eocene boundary. Science 310: 993–996.
- 81. Miller IM, Brandon MT, Hickey LJ (2006) Using leaf margin analysis to estimate the mid-Cretaceous (Albian) paleolatitude of the Baja BC block. Earth Planet Sci Lett 245: 95–114.
- 82. Uhl D, Mosbrugger V, Bruch A, Utescher T (2003) Reconstructing palaeotemperatures using leaf floras - case studies for a comparison of leaf margin analysis and the coexistence approach. Rev Palaeobot Palynol 126: 49–64.
- 83. Liang M-M, Bruch A, Collinson M, Mosbrugger V, Li C-S, et al. (2003) Testing the climatic estimates from different palaeobotanical methods: an example from the middle Miocene Shanwang flora of China. Palaeogeogr Palaeoclimatol Palaeoecol 198: 279–301.
- 84. Peppe DJ, Royer DL, Cariglino B, Oliver SY, Newman S, et al. (2010) Sensitivity of leaf size and shape to climate: global patterns and paleoclimatic applications. New Phytol in press.
- 85. Doyle JA (2007) Systematic value and evolution of leaf architecture across the angiosperms in light of molecular phylogenetic analyses. Cour Forschungsinst Senckenb 258: 21–37.
- 86. Feild TS, Arens NC (2005) Form, function and environments of the early angiosperms: merging extant phylogeny and ecophysiology with fossils. New Phytol 166: 383–408.