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.

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.

Holocene shifts in the assembly of plant and animal communities implicate human impacts.

Understanding how ecological communities are organized and how they change through time is critical to predicting the effects of climate change. Recent work documenting the co-occurrence structure of modern communities found that most significant species pairs co-occur less frequently than would be expected by chance. However, little is known about how co-occurrence structure changes through time. Here we evaluate changes in plant and animal community organization over geological time by quantifying the co-occurrence structure of 359,896 unique taxon pairs in 80 assemblages spanning the past 300 million years. Co-occurrences of most taxon pairs were statistically random, but a significant fraction were spatially aggregated or segregated. Aggregated pairs dominated from the Carboniferous period (307 million years ago) to the early Holocene epoch (11,700 years before present), when there was a pronounced shift to more segregated pairs, a trend that continues in modern assemblages. The shift began during the Holocene and coincided with increasing human population size and the spread of agriculture in North America. Before the shift, an average of 64% of significant pairs were aggregated; after the shift, the average dropped to 37%. The organization of modern and late Holocene plant and animal assemblages differs fundamentally from that of assemblages over the past 300 million years that predate the large-scale impacts of humans. Our results suggest that the rules governing the assembly of communities have recently been changed by human activity.

StomataCounter: a neural network for automatic stomata identification and counting.

Stomata regulate important physiological processes in plants and are often phenotyped by researchers in diverse fields of plant biology. Currently, there are no user-friendly, fully automated methods to perform the task of identifying and counting stomata, and stomata density is generally estimated by manually counting stomata. We introduce StomataCounter, an automated stomata counting system using a deep convolutional neural network to identify stomata in a variety of different microscopic images. We use a human-in-the-loop approach to train and refine a neural network on a taxonomically diverse collection of microscopic images. Our network achieves 98.1% identification accuracy on Ginkgo scanning electron microscropy micrographs, and 94.2% transfer accuracy when tested on untrained species. To facilitate adoption of the method, we provide the method in a publicly available website at http://www.stomata.science/.

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.

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.

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.

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.

Scaling and structure of dicotyledonous leaf venation networks.

There have been numerous attempts to derive general models for the structure and function of resource delivery networks in biology. Such theories typically predict the quantitative structure of vascular networks across scales. For example, fractal branching models of plant structure predict that the network dimensions within plant stems or leaves should be scale-free. However, very few empirical examples of such networks are available with which to evaluate such hypotheses. Here, we apply recently developed leaf network extraction software to a global leaf dataset. We find that leaf networks are neither entirely scale-free nor governed entirely by a characteristic scale. Indeed, we find many network properties, such as vein length distributions, which are governed by characteristic scales, and other network properties, notably vein diameter distributions, which are typified by power-law behaviour. Our findings suggest that theories of network structure will remain incomplete until they address the multiple constraints on network architecture.

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.

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.

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.

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.

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.

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.