Isotopic clumping in wood as a proxy for photorespiration in trees.

Photorespiration can limit gross primary productivity in terrestrial plants. The rate of photorespiration relative to carbon fixation increases with temperature and decreases with atmospheric [CO2]. However, the extent to which this rate varies in the environment is unclear. Here, we introduce a proxy for relative photorespiration rate based on the clumped isotopic composition of methoxyl groups (R-O-CH3) in wood. Most methoxyl C-H bonds are formed either during photorespiration or the Calvin cycle and thus their isotopic composition may be sensitive to the mixing ratio of these pathways. In water-replete growing conditions, we find that the abundance of the clumped isotopologue 13CH2D correlates with temperature (18-28 °C) and atmospheric [CO2] (280-1000 ppm), consistent with a common dependence on relative photorespiration rate. When applied to a global dataset of wood, we observe global trends of isotopic clumping with climate and water availability. Clumped isotopic compositions are similar across environments with temperatures below ~18 °C. Above ~18 °C, clumped isotopic compositions in water-limited and water-replete trees increasingly diverge. We propose that trees from hotter climates photorespire substantially more than trees from cooler climates. How increased photorespiration is managed depends on water availability: water-replete trees export more photorespiratory metabolites to lignin whereas water-limited trees either export fewer overall or direct more to other sinks that mitigate water stress. These disparate trends indicate contrasting responses of photorespiration rate (and thus gross primary productivity) to a future high-[CO2] world. This work enables reconstructing photorespiration rates in the geologic past using fossil wood.

The PhanSST global database of Phanerozoic sea surface temperature proxy data.

Paleotemperature proxy data form the cornerstone of paleoclimate research and are integral to understanding the evolution of the Earth system across the Phanerozoic Eon. Here, we present PhanSST, a database containing over 150,000 data points from five proxy systems that can be used to estimate past sea surface temperature. The geochemical data have a near-global spatial distribution and temporally span most of the Phanerozoic. Each proxy value is associated with consistent and queryable metadata fields, including information about the location, age, and taxonomy of the organism from which the data derive. To promote transparency and reproducibility, we include all available published data, regardless of interpreted preservation state or vital effects. However, we also provide expert-assigned diagenetic assessments, ecological and environmental flags, and other proxy-specific fields, which facilitate informed and responsible reuse of the database. The data are quality control checked and the foraminiferal taxonomy has been updated. PhanSST will serve as a valuable resource to the paleoclimate community and has myriad applications, including evolutionary, geochemical, diagenetic, and proxy calibration studies.

Diversification in the Rosales is influenced by dispersal, geographic range size, and pre-existing species richness.

PREMISE: Biodiversity results from origination and extinction, justifying interest in identifying traits that influence this balance. Traits implicated in the success or failure of lineages include dispersal, colonization ability, and geographic range size. We investigated the impact of dispersal and range size on contemporary diversity in the Rosales. METHODS: We used the multiple-state speciation and extinction (MuSSE) method to explore the effects on genus-level diversification of two genus-level traits (geographic range size and within-genus proclivity to speciate) and two species traits (seed dispersal and growth habit) and the multiple hidden-state speciation and extinction (MuHiSSE) method for species-level associations. Finally, we conducted a PGLS (phylogenetic least-squares) analysis to distinguish between speciation within genera versus origination of new genera. RESULTS: At the species level, animal dispersal enhances diversification rate in both woody and herbaceous lineages, while woody lineages without animal dispersal have higher extinction rates than speciation rates. At the genus level, herbaceous taxa have positive diversification rates regardless of other character states. Diversification rate variation is also explained by two interactions: (1) a three-way interaction between large geographic range, animal-mediated dispersal, and high within-genus species richness, whereby genera possessing all three traits have high diversification rates, and (2) a four-way interaction by which the three-way interaction is stronger in woody genera than in herbaceous genera. CONCLUSIONS: Colonization ability may underlie the relationship between dispersal type and range size and may influence past diversification rates by decreasing extinction rates during late Cenozoic climate volatility. Thus, colonization ability could be used to predict future extinction risk to aid conservation.

Fossil papilionoids of the Bowdichia clade (Leguminosae) from the Paleogene of North America.

PREMISE: Understanding the evolutionary history of flowering plants has been enriched by the integration of molecular phylogenies and evidence from the fossil record. Fossil fruits and leaves from the late Paleocene and Eocene of Wyoming and Eocene of Kentucky and Tennessee are described as extinct genera in the tropical American Bowdichia clade of the legume subfamily Papilionoideae. Recent phylogenetic study and taxonomic revision of the Bowdichia clade have facilitated understanding of relationships of the fossil taxa and their evolutionary implications and paleoenvironmental significance. METHODS: The fossils were studied using standard methods of specimen preparation and light microscopy and compared to fruits and leaves from extant legume taxa using herbarium collections. Phylogenetic relationships of the fossil taxa were assessed using morphology and DNA sequence data. RESULTS: Two new fossil genera are described and their phylogenetic relationships are established. Paleobowdichia lamarensis is placed as sister to the extant genus Bowdichia and Tobya claibornensis is placed with the extant genera Guianodendron and Staminodianthus. CONCLUSIONS: These fossils demonstrate that the tropical American Bowdichia clade was present in North America during a period when tropical or subtropical conditions prevailed in the northern Rocky Mountains during the late Paleocene and the Mississippi Embayment during the middle Eocene. These fossils also document that the Bowdichia clade had diversified by the late Paleocene when the fossil record of the family is relatively sparse. This result suggests that future work on early fossil legumes should focus on tropical and subtropical climatic zones, wherever they may occur latitudinally.

An image dataset of cleared, x-rayed, and fossil leaves vetted to plant family for human and machine learning.

Leaves are the most abundant and visible plant organ, both in the modern world and the fossil record. Identifying foliage to the correct plant family based on leaf architecture is a fundamental botanical skill that is also critical for isolated fossil leaves, which often, especially in the Cenozoic, represent extinct genera and species from extant families. Resources focused on leaf identification are remarkably scarce; however, the situation has improved due to the recent proliferation of digitized herbarium material, live-plant identification applications, and online collections of cleared and fossil leaf images. Nevertheless, the need remains for a specialized image dataset for comparative leaf architecture. We address this gap by assembling an open-access database of 30,252 images of vouchered leaf specimens vetted to family level, primarily of angiosperms, including 26,176 images of cleared and x-rayed leaves representing 354 families and 4,076 of fossil leaves from 48 families. The images maintain original resolution, have user-friendly filenames, and are vetted using APG and modern paleobotanical standards. The cleared and x-rayed leaves include the Jack A. Wolfe and Leo J. Hickey contributions to the National Cleared Leaf Collection and a collection of high-resolution scanned x-ray negatives, housed in the Division of Paleobotany, Department of Paleobiology, Smithsonian National Museum of Natural History, Washington D.C.; and the Daniel I. Axelrod Cleared Leaf Collection, housed at the University of California Museum of Paleontology, Berkeley. The fossil images include a sampling of Late Cretaceous to Eocene paleobotanical sites from the Western Hemisphere held at numerous institutions, especially from Florissant Fossil Beds National Monument (late Eocene, Colorado), as well as several other localities from the Late Cretaceous to Eocene of the Western USA and the early Paleogene of Colombia and southern Argentina. The dataset facilitates new research and education opportunities in paleobotany, comparative leaf architecture, systematics, and machine learning.

Fossil Atmospheres: a case study of citizen science in question-driven palaeontological research.

Palaeontologists increasingly use large datasets of observations collected from museum specimens to address broad-scale questions about evolution and ecology on geological timescales. One such question is whether information from fossil organisms can be used as a robust proxy for atmospheric carbon dioxide through time. Here, we present the citizen science branch of 'Fossil Atmospheres', a project designed to refine stomatal index of Ginkgo leaves as a palaeo-CO2 proxy by involving citizen scientists in data collection through the Zooniverse website. Citizen science helped to overcome a barrier presented by the time taken to count cells in Ginkgo samples; however, a new set of challenges arose as a result. A beta-testing phase with Zooniverse volunteers provided an opportunity to improve instructions to ensure high fidelity data. Exploration of citizen scientists' estimates shows that volunteer counts of stomata are accurate with respect to counts made by the project's lead scientist. However, counts of epidermal cells have a wide range, and mean values tend to underestimate expert counts. We demonstrate a variety of approaches to reducing the inaccuracy in the calculated stomatal index that this variation causes. Zooniverse serves as an ideal tool for collection of palaeontological data where the distribution of fossils would be impossible, but where specimens can be easily imaged. Such an approach facilitates the collection of a large palaeontological dataset, as well as providing an opportunity for citizens to engage with climate research.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Binary-state speciation and extinction method is conditionally robust to realistic violations of its assumptions.

BACKGROUND: Phylogenetic comparative methods allow us to test evolutionary hypotheses without the benefit of an extensive fossil record. These methods, however, make simplifying assumptions, among them that clades are always increasing or stable in diversity, an assumption we know to be false. This study simulates hypothetical clades to test whether the Binary State Speciation and Extinction (BiSSE) method can be used to correctly detect relative differences in diversification rate between ancestral and derived character states even as net diversification rates are declining overall. We simulate clades with declining but positive diversification rates, as well those in which speciation rates decline below extinction rates so that they are losing richness for part of their history. We run these analyses both with simulated symmetric and asymmetric speciation rates to test whether BiSSE can be used to detect them correctly. RESULTS: For simulations with a neutral character, the fit for a BiSSE model with a neutral character is better than alternative models so long as net diversification rates remain positive. Once net diversification rates become negative, the BiSSE model with the greatest likelihood often has a non-neutral character, even though there is no such character in the simulation. BiSSE's usefulness in detecting real asymmetry in speciation rates improves with clade age, even well after net diversification rates have become negative. CONCLUSIONS: BiSSE is most useful in analyzing clades of intermediate age, before they have reached peak diversity and gone into decline. After this point, users of BiSSE risk incorrectly inferring differential evolutionary rates when none exist. Fortunately, most studies using BiSSE and similar models focus on rapid, recent diversifications, and are less likely to encounter the biases BiSSE models are subject to for older clades. For extant groups that were once more diverse than now, however, caution should be taken in inferring past diversification patterns without fossil data.

Consequences of elevated temperature and pCO2 on insect folivory at the ecosystem level: perspectives from the fossil record.

Paleoecological studies document the net effects of atmospheric and climate change in a natural laboratory over timescales not accessible to laboratory or ecological studies. Insect feeding damage is visible on well-preserved fossil leaves, and changes in leaf damage through time can be compared to environmental changes. We measured percent leaf area damaged on four fossil leaf assemblages from the Bighorn Basin, Wyoming, that range in age from 56.1 to 52.65 million years (Ma). We also include similar published data from three US sites 49.4 to ~45 Ma in our analyses. Regional climate was subtropical or warmer throughout this period, and the second oldest assemblage (56 Ma) was deposited during the Paleocene-Eocene Thermal Maximum (PETM), a geologically abrupt global warming event caused by massive release of carbon into the atmosphere. Total and leaf-chewing damage are highest during the PETM, whether considering percent area damaged on the bulk flora, the average of individual host plants, or a single plant host that occurs at multiple sites. Another fossil assemblage in our study, the 52.65 Ma Fifteenmile Creek paleoflora, also lived during a period of globally high temperature and pCO 2, but does not have elevated herbivory. Comparison of these two sites, as well as regression analyses conducted on the entire dataset, demonstrates that, over long timescales, temperature and pCO 2 are uncorrelated with total insect consumption at the ecosystem level. Rather, the most important factor affecting herbivory is the relative abundance of plants with nitrogen-fixing symbionts. Legumes dominate the PETM site; their prevalence would have decreased nitrogen limitation across the ecosystem, buffering generalist herbivore populations against decreased leaf nutritional quality that commonly occurs at high pCO 2. We hypothesize that nitrogen concentration regulates the opposing effects of elevated temperature and CO 2 on insect abundance and thereby total insect consumption, which has important implications for agricultural practices in today's world of steadily increasing pCO 2.

Computer vision cracks the leaf code.

Understanding the extremely variable, complex shape and venation characters of angiosperm leaves is one of the most challenging problems in botany. Machine learning offers opportunities to analyze large numbers of specimens, to discover novel leaf features of angiosperm clades that may have phylogenetic significance, and to use those characters to classify unknowns. Previous computer vision approaches have primarily focused on leaf identification at the species level. It remains an open question whether learning and classification are possible among major evolutionary groups such as families and orders, which usually contain hundreds to thousands of species each and exhibit many times the foliar variation of individual species. Here, we tested whether a computer vision algorithm could use a database of 7,597 leaf images from 2,001 genera to learn features of botanical families and orders, then classify novel images. The images are of cleared leaves, specimens that are chemically bleached, then stained to reveal venation. Machine learning was used to learn a codebook of visual elements representing leaf shape and venation patterns. The resulting automated system learned to classify images into families and orders with a success rate many times greater than chance. Of direct botanical interest, the responses of diagnostic features can be visualized on leaf images as heat maps, which are likely to prompt recognition and evolutionary interpretation of a wealth of novel morphological characters. With assistance from computer vision, leaves are poised to make numerous new contributions to systematic and paleobotanical studies.

Reading the leaves: A comparison of leaf rank and automated areole measurement for quantifying aspects of leaf venation.

The reticulate venation that is characteristic of a dicot leaf has excited interest from systematists for more than a century, and from physiological and developmental botanists for decades. The tools of digital image acquisition and computer image analysis, however, are only now approaching the sophistication needed to quantify aspects of the venation network found in real leaves quickly, easily, accurately, and reliably enough to produce biologically meaningful data. In this paper, we examine 120 leaves distributed across vascular plants (representing 118 genera and 80 families) using two approaches: a semiquantitative scoring system called "leaf ranking," devised by the late Leo Hickey, and an automated image-analysis protocol. In the process of comparing these approaches, we review some methodological issues that arise in trying to quantify a vein network, and discuss the strengths and weaknesses of automatic data collection and human pattern recognition. We conclude that subjective leaf rank provides a relatively consistent, semiquantitative measure of areole size among other variables; that modal areole size is generally consistent across large sections of a leaf lamina; and that both approaches-semiquantitative, subjective scoring; and fully quantitative, automated measurement-have appropriate places in the study of leaf venation.

Biomechanical and leaf-climate relationships: a comparison of ferns and seed plants.

PREMISE OF THE STUDY: Relationships of leaf size and shape (physiognomy) with climate have been well characterized for woody non-monocotyledonous angiosperms (dicots), allowing the development of models for estimating paleoclimate from fossil leaves. More recently, petiole width of seed plants has been shown to scale closely with leaf mass. By measuring petiole width and leaf area in fossils, leaf mass per area (MA) can be estimated and an approximate leaf life span inferred. However, little is known about these relationships in ferns, a clade with a deep fossil record and with the potential to greatly expand the applicability of these proxies. METHODS: We measured the petiole width, MA, and leaf physiognomic characters of 179 fern species from 188 locations across six continents. We applied biomechanical models and assessed the relationship between leaf physiognomy and climate using correlational approaches. KEY RESULTS: The scaling relationship between area-normalized petiole width and MA differs between fern fronds and pinnae. The scaling relationship is best modeled as an end-loaded cantilevered beam, which is different from the best-fit biomechanical model for seed plants. Fern leaf physiognomy is not influenced by climatic conditions. CONCLUSIONS: The cantilever beam model can be applied to fossil ferns. The lack of sensitivity of leaf physiognomy to climate in ferns argues against their use to reconstruct paleoclimate. Differences in climate sensitivity and biomechanical relationships between ferns and seed plants may be driven by differences in their hydraulic conductivity and/or their differing evolutionary histories of vein architecture and leaf morphology.

Plant response to a global greenhouse event 56 million years ago.

PREMISE OF THE STUDY: The fossil record provides information about the long-term response of plants to CO2-induced climate change. The Paleocene-Eocene Thermal Maximum (PETM), a 200000-yr-long period of rapid carbon release and warming that occurred ∼56 million years ago, is analogous to future anthropogenic global warming. METHODS: We collected plant macrofossils in the Bighorn Basin, Wyoming, United States, from a period spanning the PETM and studied changes in floristic composition. We also compiled and summarized published records of floristic change during the PETM. KEY RESULTS: There was radical floristic change in the Bighorn Basin during the PETM reflecting local or regional extirpation of mesophytic plants, notably conifers, and colonization of the area by thermophilic and dry-tolerant species, especially Fabaceae. This floristic change largely reversed itself as the PETM ended, though some immigrant species persisted and some Paleocene species never returned. Less detailed records from other parts of the world show regional variation in floristic response, but are mostly consistent with the Bighorn Basin trends. CONCLUSIONS: Despite geologically rapid extirpation, colonization, and recolonization, we detected little extinction during the PETM, suggesting the rate of climate change did not exceed the dispersal capacity of terrestrial plants. Extrapolating the response of plants from the PETM to future anthropogenic climate change likely underestimates risk because rates of climate change during the PETM may have been an order of magnitude slower than current rates of change and because the abundant, widespread species common as fossils are likely resistant to extinction.

Evolution of the earliest horses driven by climate change in the Paleocene-Eocene Thermal Maximum.

Body size plays a critical role in mammalian ecology and physiology. Previous research has shown that many mammals became smaller during the Paleocene-Eocene Thermal Maximum (PETM), but the timing and magnitude of that change relative to climate change have been unclear. A high-resolution record of continental climate and equid body size change shows a directional size decrease of ~30% over the first ~130,000 years of the PETM, followed by a ~76% increase in the recovery phase of the PETM. These size changes are negatively correlated with temperature inferred from oxygen isotopes in mammal teeth and were probably driven by shifts in temperature and possibly high atmospheric CO(2) concentrations. These findings could be important for understanding mammalian evolutionary responses to future global warming.

Paleocene Malvaceae from northern South America and their biogeographical implications.

PREMISE OF THE STUDY: The clade Bombacoideae + Malvoideae ('Malvatheca group' sensu Baum et al.) in Malvaceae comprises a mostly tropical lineage with derived taxa that now thrive in higher latitudes. The sparse fossil record, especially for Malvoideae, obscures interpretations of past distributions. We describe fossil leaves of Malvoideae from the middle-late Paleocene Cerrejón Formation in Colombia, which contains evidence for the earliest known neotropical rainforest. METHODS: Fifty-six leaf compressions belonging to Malvaceae were collected from the Cerrejón Formation in northern Colombia. Leaf architectural characters were scored and optimized for 81 genera of Malvaceae. Synapomorphic characters and unique character combinations support natural affinities for the fossil leaves. Fossil pollen from the same formation was also assessed. KEY RESULTS: Despite convergence of overall leaf architecture among many Malvaceae, Malvaciphyllum macondicus sp. nov. can be assigned to the clade Eumalvoideae because of distal and proximal bifurcations of the costal secondary and agrophic veins, a synapomorphy for this clade. CONCLUSIONS: The leaf compressions, the oldest fossils for Eumalvoideae, indicate a minimum divergence time of 58-60 Ma, older than existing estimates from molecular analyses of living species. The abundance of eumalvoid leaves and of bombacoid pollen in the midlate Paleocene of Colombia suggests that the Malvatheca group (Malvoideae + Bombacoideae) was already a common element in neotropical forests and does not support an Australasian origin for Eumalvoideae.

Sensitivity of leaf size and shape to climate: global patterns and paleoclimatic applications.

• Paleobotanists have long used models based on leaf size and shape to reconstruct paleoclimate. However, most models incorporate a single variable or use traits that are not physiologically or functionally linked to climate, limiting their predictive power. Further, they often underestimate paleotemperature relative to other proxies. • Here we quantify leaf-climate correlations from 92 globally distributed, climatically diverse sites, and explore potential confounding factors. Multiple linear regression models for mean annual temperature (MAT) and mean annual precipitation (MAP) are developed and applied to nine well-studied fossil floras. • We find that leaves in cold climates typically have larger, more numerous teeth, and are more highly dissected. Leaf habit (deciduous vs evergreen), local water availability, and phylogenetic history all affect these relationships. Leaves in wet climates are larger and have fewer, smaller teeth. Our multivariate MAT and MAP models offer moderate improvements in precision over univariate approaches (± 4.0 vs 4.8°C for MAT) and strong improvements in accuracy. For example, our provisional MAT estimates for most North American fossil floras are considerably warmer and in better agreement with independent paleoclimate evidence. • Our study demonstrates that the inclusion of additional leaf traits that are functionally linked to climate improves paleoclimate reconstructions. This work also illustrates the need for better understanding of the impact of phylogeny and leaf habit on leaf-climate relationships.

Global patterns in leaf 13C discrimination and implications for studies of past and future climate.

Fractionation of carbon isotopes by plants during CO(2) uptake and fixation (Delta(leaf)) varies with environmental conditions, but quantitative patterns of Delta(leaf) across environmental gradients at the global scale are lacking. This impedes interpretation of variability in ancient terrestrial organic matter, which encodes climatic and ecological signals. To address this problem, we converted 3,310 published leaf delta(13)C values into mean Delta(leaf) values for 334 woody plant species at 105 locations (yielding 570 species-site combinations) representing a wide range of environmental conditions. Our analyses reveal a strong positive correlation between Delta(leaf) and mean annual precipitation (MAP; R(2) = 0.55), mirroring global trends in gross primary production and indicating stomatal constraints on leaf gas-exchange, mediated by water supply, are the dominant control of Delta(leaf) at large spatial scales. Independent of MAP, we show a lesser, negative effect of altitude on Delta(leaf) and minor effects of temperature and latitude. After accounting for these factors, mean Delta(leaf) of evergreen gymnosperms is lower (by 1-2.7 per thousand) than for other woody plant functional types (PFT), likely due to greater leaf-level water-use efficiency. Together, environmental and PFT effects contribute to differences in mean Delta(leaf) of up to 6 per thousand between biomes. Coupling geologic indicators of ancient precipitation and PFT (or biome) with modern Delta(leaf) patterns has potential to yield more robust reconstructions of atmospheric delta(13)C values, leading to better constraints on past greenhouse-gas perturbations. Accordingly, we estimate a 4.6 per thousand decline in the delta(13)C of atmospheric CO(2) at the onset of the Paleocene-Eocene Thermal Maximum, an abrupt global warming event approximately 55.8 Ma.

Late Paleocene fossils from the Cerrejon Formation, Colombia, are the earliest record of Neotropical rainforest.

Neotropical rainforests have a very poor fossil record, making hypotheses concerning their origins difficult to evaluate. Nevertheless, some of their most important characteristics can be preserved in the fossil record: high plant diversity, dominance by a distinctive combination of angiosperm families, a preponderance of plant species with large, smooth-margined leaves, and evidence for a high diversity of herbivorous insects. Here, we report on an approximately 58-my-old flora from the Cerrejón Formation of Colombia (paleolatitude approximately 5 degrees N) that is the earliest megafossil record of Neotropical rainforest. The flora has abundant, diverse palms and legumes and similar family composition to extant Neotropical rainforest. Three-quarters of the leaf types are large and entire-margined, indicating rainfall >2,500 mm/year and mean annual temperature >25 degrees C. Despite modern family composition and tropical paleoclimate, the diversity of fossil pollen and leaf samples is 60-80% that of comparable samples from extant and Quaternary Neotropical rainforest from similar climates. Insect feeding damage on Cerrejón fossil leaves, representing primary consumers, is abundant, but also of low diversity, and overwhelmingly made by generalist feeders rather than specialized herbivores. Cerrejón megafossils provide strong evidence that the same Neotropical rainforest families have characterized the biome since the Paleocene, maintaining their importance through climatic phases warmer and cooler than present. The low diversity of both plants and herbivorous insects in this Paleocene Neotropical rainforest may reflect an early stage in the diversification of the lineages that inhabit this biome, and/or a long recovery period from the terminal Cretaceous extinction.

Palms (Arecaceae) from a Paleocene rainforest of northern Colombia.

Palms are a monophyletic group with a dominantly tropical distribution; however, their fossil record in low latitudes is strikingly scarce. In this paper, we describe fossil leaves, inflorescences, and fruits of palms from the middle to late Paleocene Cerrejón Formation, outcropping in the Ranchería River Valley, northern Colombia. The fossils demonstrate the presence of at least five palm morphospecies in the basin ca. 60 Ma. We compare the morphology of the fossils with extant palms and conclude that they belong to at least three palm lineages: the pantropical Cocoseae of the subfamily Arecoideae, the monotypic genus Nypa, and either Calamoideae or Coryphoideae. The fossil fruits and inflorescences are among the oldest megafossil records of these groups and demonstrate that the divergence of the Cocoseae was more than 60 Ma, earlier than has previously been thought. These fossils are useful in tracing the range expansion or contraction of historical or current neotropical elements and also have profound implications for the understanding of the evolution of neotropical rainforests.

Sharply increased insect herbivory during the Paleocene-Eocene Thermal Maximum.

The Paleocene-Eocene Thermal Maximum (PETM, 55.8 Ma), an abrupt global warming event linked to a transient increase in pCO2, was comparable in rate and magnitude to modern anthropogenic climate change. Here we use plant fossils from the Bighorn Basin of Wyoming to document the combined effects of temperature and pCO2 on insect herbivory. We examined 5,062 fossil leaves from five sites positioned before, during, and after the PETM (59-55.2 Ma). The amount and diversity of insect damage on angiosperm leaves, as well as the relative abundance of specialized damage, correlate with rising and falling temperature. All reach distinct maxima during the PETM, and every PETM plant species is extensively damaged and colonized by specialized herbivores. Our study suggests that increased insect herbivory is likely to be a net long-term effect of anthropogenic pCO2 increase and warming temperatures.