Historians have debated whether pathways and events from the past to the present are influenced largely by contingency, the dependence of outcomes on particular prior conditions, or whether there is long-term emergent directional change. Previous arguments for predictability in evolutionary history relied on the high frequency of convergence, but the repeated evolution of widely favored adaptations need not imply long-term directionality. Using evidence from the fossil record and arguments concerning the metabolic evolution of organisms, I show here that power (total energy taken up and expended per unit time) has increased stepwise over time at ecosystem-level and global scales thanks to the ratchet-like, cumulative effects of competition and cooperation and to the disproportionate influence of powerful top competitors and opportunistic species on emergent ecosystem properties and processes. The history of life therefore exhibits emergent directionality at large ecosystem-wide scales toward greater power.

The Ordovician was a time of drastic biological and geological change. Previous work has suggested that there was a dramatic increase in global diversity during this time, but also has indicated that regional dynamics and dynamics in specific environments might have been different. Here, we contrast two paleocontinents that have different geological histories through the Ordovician, namely Laurentia and Baltica. The first was situated close to the equator throughout the whole Ordovician, while the latter has traversed tens of latitudes during the same time. We predict that Baltica, which was under long-term environmental change, would show greater average and interval-to-interval origination and extinction rates than Laurentia. In addition, we are interested in the role of the environment in which taxa originated, specifically, the patterns of onshore–offshore dynamics of diversification, where onshore and offshore areas represent high-energy and low-energy environments, respectively. Here, we predict that high-energy environments might be more conducive for originations.

Our new analyses show that the global Ordovician spike in genus richness from the Dapingian to the Darriwilian Stage resulted from a very high origination rate at the Dapingian/Darriwilian boundary, while the extinction rate remained low. We found substantial interval-to-interval variation in the origination and extinction rates in Baltica and Laurentia, but the probabilities of origination and extinction are somewhat higher in Baltica than Laurentia. Onshore and offshore areas have largely indistinguishable origination and extinction rates, in contradiction to our predictions. The global spike in origination rates at the Dapingian/Darriwilian boundary is apparent in Baltica, Laurentia, and onshore and offshore areas, and abundant variability in diversification rates is apparent over other time intervals for these paleocontinents and paleoenvironments. This observation hints at global mechanisms for the spike in origination rates at the Dapingian/Darriwilian boundary but a domination of more regional and local mechanisms over other time intervals in the Ordovician.

Farming is a behavior in which an organism promotes the growth and reproduction of other organisms in or on a substrate as a food source. A number of trace fossils have been suggested to record the occurrence of farming behavior. These include the deep-sea graphoglyptid trace fossils, proposed to be microbial farms on the seafloor, and terrestrial fossil social insect nests thought to represent fungicultural behavior. The presumed farming behavior of graphoglyptids is the basis of the ethological category agrichnia. Four criteria have been proposed as diagnostic of farming behavior, and these can be applied to both observed modern and proposed trace fossil examples of farming behavior. The evidence for farming behavior in the social insect trace record is strong but is much weaker in the case of graphoglyptids. The use of agrichnia as an ethological category should be limited to well-supported cases.

Evolutionary history studies depend on having reliable chronologies of macroevolutionary processes. Construction of such chronologies often yields discrepancies between paleontological and molecular dates, which are sometimes viewed as conflicting. Nevertheless, each macroevolutionary process is composed of two main phases: emergence of a trait or clade and success of that trait or clade, which differ in mechanisms, drivers, and types of evidence. Moreover, emergence may be observed as gene divergence (which may be trait-coding or trait-unrelated genes), trait emergence, and clade emergence; whereas success can be observed as increase in abundance, diffusion, and/or diversity or as overall persistence over geologic time. Therefore, to fully and correctly understand any macroevolutionary process, it is of paramount importance to understand what event each date refers to, and how dates of various events and their integration reveal the complexity of macroevolutionary processes. I demonstrate this through three examples: the chronological gap between oxygenic photosynthesis emergence and the Great Oxidation Event, the chronological gap between paleontological and molecular dates of angiosperm emergence, and the evolution of plant silicon accumulation.

Organismal metabolic rates reflect the interaction of environmental and physiological factors. Thus, calcifying organisms that record growth history can provide insight into both the ancient environments in which they lived and their own physiology and life history. However, interpreting them requires understanding which environmental factors have the greatest influence on growth rate and the extent to which evolutionary history constrains growth rates across lineages. We integrated satellite measurements of sea-surface temperature and chlorophyll-a concentration with a database of growth coefficients, body sizes, and life spans for 692 populations of living marine bivalves in 195 species, set within the context of a new maximum-likelihood phylogeny of bivalves. We find that environmental predictors overall explain only a small proportion of variation in growth coefficient across all species; temperature is a better predictor of growth coefficient than food supply, and growth coefficient is somewhat more variable at higher summer temperatures. Growth coefficients exhibit moderate phylogenetic signal, and taxonomic membership is a stronger predictor of growth coefficient than any environmental predictor, but phylogenetic inertia cannot fully explain the disjunction between our findings and the extensive body of work demonstrating strong environmental control on growth rates within taxa. Accounting for evolutionary history is critical when considering shells as historical archives. The weak relationship between variation in food supply and variation in growth coefficient in our data set is inconsistent with the hypothesis that the increase in mean body size through the Phanerozoic was driven by increasing productivity enabling faster growth rates.

Recent excavations of Ediacaran assemblages have revealed striking bed-to-bed variation in diversity–abundance structure, offering potential insight into the ecology and taphonomy of these poorly understood early Metazoan ecosystems. Here we compare faunal variability in Ediacaran assemblages to that of younger benthic assemblages, both fossil and modern. We decompose the diversity of local assemblages into within-collection (α) and among-collection (β) components and show that β diversity in Ediacaran assemblages is unusually high relative to younger assemblages. Average between-bed ecological dissimilarities in the Phanerozoic fossil record are comparable to within-habitat dissimilarities typically observed over meter to kilometer scales in modern benthic marine habitats, but dissimilarities in Ediacaran assemblages are comparable to those typically observed over 10–100 km scales in modern habitats. We suggest that the unusually variable diversity–abundance structure of Ediacaran assemblages is due both to their preservation as near snapshots of benthic communities and to original ecological differences, in particular the paucity of motile taxa and the near lack of predation and infaunalization.

Previous studies of cranial shape have established a consistent interspecific allometric pattern relating the relative lengths of the face and braincase regions of the skull within multiple families of mammals. In this interspecific allometry, the facial region of the skull is proportionally longer than the braincase in larger species. The regularity and broad taxonomic occurrence of this allometric pattern suggests that it may have an origin near the base of crown Mammalia, or even deeper in the synapsid or amniote forerunners of mammals. To investigate the possible origins of this allometric pattern, we used geometric morphometric techniques to analyze cranial shape in 194 species of nonmammalian synapsids, which constitute a set of successive outgroups to Mammalia. We recovered a much greater diversity of allometric patterns within nonmammalian synapsids than has been observed in mammals, including several instances similar to the mammalian pattern. However, we found no evidence of the mammalian pattern within Therocephalia and nonmammalian Cynodontia, the synapsids most closely related to mammals. This suggests that the mammalian allometric pattern arose somewhere within Mammaliaformes, rather than within nonmammalian synapsids. Further investigation using an ontogenetic series of the anomodont Diictodon feliceps shows that the pattern of interspecific allometry within anomodonts parallels the ontogenetic trajectory of Diictodon. This indicates that in at least some synapsids, allometric patterns associated with ontogeny may provide a “path of least resistance” for interspecific variation, a mechanism that we suggest produces the interspecific allometric pattern observed in mammals.

Toothed cetacean (Odontoceti) lineages in the Miocene and Pliocene evolved rostra that are proportionally more elongate than any other aquatic mammal or reptile, living or extinct. Their similarities in cranial proportions to billfish may suggest a convergent feeding style, where the rostrum is swept through the water to hit and stun prey. Here we calculated second moment of area from rostral cross sections of these fossil odontocete taxa, as well as from extant ecological analogues, to infer variation in feeding behavior. Our results show that the extremely long rostra of extinct toothed whales vary considerably in functionally relevant measures of shape and likely exhibited a diversity of feeding behaviors, ranging from those similar to modern odontocetes to those convergent with billfish. Eustatic sea-level and temperature maxima of the Miocene likely led to changes in prey characteristics or abundance that enabled the repeated evolution of this extreme morphotype, which later went extinct during late Pliocene climatic deterioration.

The estimation of origination and extinction rates and their temporal variation is central to understanding diversity patterns and the evolutionary history of clades. The fossil record provides the only direct evidence of extinction and biodiversity changes through time and has long been used to infer the dynamics of diversity changes in deep time. The software PyRate implements a Bayesian framework to analyze fossil occurrence data to estimate the rates of preservation, origination, and extinction while incorporating several sources of uncertainty. Building upon this framework, we present a suite of methodological advances including more complex and realistic models of preservation and the first likelihood-based test to compare the fit across different models. Further, we develop a new reversible jump Markov chain Monte Carlo algorithm to estimate origination and extinction rates and their temporal variation, which provides more reliable results and includes an explicit estimation of the number and temporal placement of statistically significant rate changes. Finally, we implement a new C++ library that speeds up the analyses by orders of magnitude, therefore facilitating the application of the PyRate methods to large data sets. We demonstrate the new functionalities through extensive simulations and with the analysis of a large data set of Cenozoic marine mammals. We compare our analytical framework against two widely used alternative methods to infer origination and extinction rates, revealing that PyRate decisively outperforms them across a range of simulated data sets. Our analyses indicate that explicit statistical model testing, which is often neglected in fossil-based macroevolutionary analyses, is crucial to obtain accurate and robust results.

Gastropods often show signs of unsuccessful attacks by durophagous predators in the form of healed scars in their shells. As such, fossil gastropods can be taken as providing a record of predation through geological time. However, interpreting the number of such scars has proved to be problematic—Would a low number of scars mean a low rate of attack or a high rate of success, for example? Here we develop a model of population dynamics among individuals exposed to predation, including both lethal and nonlethal attacks. Using this model, we calculate the equilibrium distributions of ages and healed scars in the population and among fossilized specimens, based on the assumption that predation is independent of age or scar number. Based on these results, we formally show that the rates of attack and success cannot be disambiguated without further information about population structure. Nevertheless, by making the assumptions that the non-durophagous predatory death rate is both constant and low, we show that it is possible to use relatively small assemblages of gastropods to produce accurate estimates of both attack and success rates, if the overall death rate can be estimated. We consider likely violations of the assumptions in our model and what sort of information would be required to solve this problem in these more general cases. However, it is not easy to extract the relevant information from the fossil record: a variety of important biases are likely to intervene to obscure the data that gastropod assemblages may yield. Nonetheless, the model provides a theoretical framework for interpreting summary data, including for comparison between different assemblages.

The warm, equable, and ice-free early Eocene Epoch permits investigation of ecosystem function and macro-ecological patterns during a very different climate regime than exists today. It also provides insight into what the future may entail, as anthropogenic CO2 release drives Earth toward a comparable hothouse condition. Studying plant–insect herbivore food webs during hothouse intervals is warranted, because these account for the majority of nonmicrobial terrestrial biodiversity. Here, we report new plant and insect herbivore damage census data from two floodplain sites in the Wind River Basin of central Wyoming, one in the Aycross Formation (50–48.25 Ma) at the basin edge (WRE) and the second in the Wind River Formation in the interior of the basin (WRI). The WRI site is in stratigraphic proximity to a volcanic ash that is newly dated to 52.416 ± 0.016/0.028/0.063 (2σ). We compare the Wind River Basin assemblages to published data from a 52.65 Ma floodplain flora in the neighboring Bighorn (BH) Basin and find that only 5.6% of plant taxa occur at all three sites and approximately 10% occur in both basins. The dissimilar floras support distinct suites of insect herbivores, as recorded by leaf damage. The relatively low-diversity BH flora has the highest diversity of insect damage, contrary to hypotheses that insect herbivore diversity tracks floral diversity. The distinctiveness of the WRE flora is likely due to its younger age and cooler reconstructed paleotemperature, but these factors are nearly identical for the WRI and BH floras. Site-specific microenvironmental factors that cannot be measured easily in deep time may account for these differences. Alternatively, the Owl Creek Mountains between the two basins may have provided a formidable barrier to the thermophilic organisms that inhabited the basin interiors, supporting Janzen's hypothesis that mountain passes appear higher in tropical environments.

To assess evolutionary processes in deep time, it is essential to understand the roles of development and environment, both recorded through the morphological variability of fossil assemblages. Thanks to their great abundance and the high temporal resolution of their fossil record, conodont elements are ideal to address this issue. In this paper, we present the first quantitative study of a Carnian–Norian (Late Triassic) assemblage of closely related P1 conodont elements. Using geometric morphometrics (landmarks, sliding landmarks, and elliptic Fourier analysis), we explore the main axes of phenotypic variation and relate them to classically used taxonomic characters. We show that some important taxonomic features follow laws of covariation, hence highlighting developmental constraints. Furthermore, the intraspecific variation within all considered species, either Carnian or Norian forms, is similarly restricted, emphasizing, for the first time in conodont P1 elements, a common line of least resistance to evolution, which means that similar intrinsic (developmental) factors were acting on these taxa and likely biased the evolutionary trajectories of all these taxa in a similar way. Because the evolution between Carnian and Norian forms is known to have followed a trajectory that is significantly different from the line of least resistance, strong extrinsic pressures, such as environmental disturbances, were probably at play around the Carnian/Norian boundary to counteract the effects of these intrinsic, developmental constraints.

The time separating the first appearances of species from their divergences from related taxa affects assessments of macroevolutionary hypotheses about rates of anatomical or ecological change. Branch durations necessarily posit stratigraphic gaps in sampling within a clade over which we have failed to sample predecessors (ancestors) and over which there are no divergences leading to sampled relatives (sister taxa). The former reflects only sampling rates, whereas the latter reflects sampling, origination, and extinction rates. Because all three rates vary over time, the probability of a branch duration of any particular length will differ depending on when in the Phanerozoic that branch duration spans. Here, I present a birth–death-sampling model allowing interval-to-interval variation in diversification and sampling rates. Increasing either origination or sampling rates increases the probability of finding sister taxa that diverge both during and before intervals of high sampling/origination. Conversely, elevated extinction reduces the probability of divergences from sampled sister taxa before and during intervals of elevated extinction. In the case of total extinction, a Signor-Lipps will reduce expected sister taxa leading up to the extinction, with the possible effect stretching back many millions of years when sampling is low. Simulations indicate that this approach provides reasonable estimates of branch duration probabilities under a variety of circumstances. Because current probability models for describing morphological evolution are less advanced than methods for inferring diversification and sampling rates, branch duration priors allowing for time-varying diversification could be a potent tool for phylogenetic inference with fossil data.

Interpreting changes in ecosystem structure from the fossil record can be challenging. In a prominent example, the traditional view that brachiopods were ecologically dominant over bivalves in the Paleozoic has been disputed on both taphonomic and metabolic grounds. Aragonitic bivalves may be underrepresented in many fossil assemblages due to preferential dissolution. Abundance counts may further understate the ecological importance of bivalves, which tend to have more biomass and higher metabolic rates than brachiopods. We evaluate the relative importance of the two clades in exceptionally preserved, bulk-sampled fossil assemblages from the Pennsylvanian Breathitt Formation of Kentucky, where aragonitic bivalves are preserved as shells, not molds. At the regional scale, brachiopods were twice as abundant as bivalves and were collectively equivalent in biomass and energy use. Analyses of samples from the Paleobiology Database that contain abundance counts are consistent with these results and show no clear trend in the relative ecological importance of bivalves during the middle and late Paleozoic. Bivalves were probably more important in Paleozoic ecosystems than is apparent in many fossil assemblages, but they were not clearly dominant over brachiopods until after the Permian–Triassic extinction, which caused the shelly benthos to shift from bivalve and brachiopod dominated to merely bivalve dominated.

Preferential dissolution of the biogenic carbonate polymorph aragonite promotes preservational bias in shelly marine faunas. While field studies have documented the impact of preferential aragonite dissolution on fossil molluscan diversity, its impact on regional and global biodiversity metrics is debated. Epicontinental seas are especially prone to conditions that both promote and inhibit preferential dissolution, which may result in spatially extensive zones with variable preservation. Here we present a multifaceted evaluation of aragonite dissolution within the Late Cretaceous Western Interior Seaway of North America. Occurrence data of mollusks from two time intervals (Cenomanian/Turonian boundary, early Campanian) are plotted on new high-resolution paleogeographies to assess aragonite preservation within the seaway. Fossil occurrences, diversity estimates, and sampling probabilities for calcitic and aragonitic fauna were compared in zones defined by depth and distance from the seaway margins. Apparent range sizes, which could be influenced by differential preservation potential of aragonite between separate localities, were also compared. Our results are consistent with exacerbated aragonite dissolution within specific depth zones for both time slices, with aragonitic bivalves additionally showing a statistically significant decrease in range size compared with calcitic fauna within carbonate-dominated Cenomanian–Turonian strata. However, we are unable to conclusively show that aragonite dissolution impacted diversity estimates. Therefore, while aragonite dissolution is likely to have affected the preservation of fauna in specific localities, time averaging and instantaneous preservation events preserve regional biodiversity. Our results suggest that the spatial expression of taphonomic biases should be an important consideration for paleontologists working on paleobiogeographic problems.

Oceanic anoxic events (OAEs) are contemporaneous with 11 of the 18 largest Phanerozoic extinction events, but the magnitude and selectivity of their paleoecological impact remains disputed. OAEs are associated with abrupt, rapid warming and increased CO2 flux to the atmosphere; thus, insights from this study may clarify the impact of current anthropogenic climate change on the biosphere. We investigated the influence of the Late Cretaceous Bonarelli event (OAE2; Cenomanian/Turonian stage boundary; ~94 Ma) on generic- and species-level molluscan diversity, extinction rates, and ecological turnover. Cenomanian/Turonian results were compared with changes across all Cretaceous stage boundaries, some of which are coincident with less severe OAEs. We found increased generic turnover, but not species-level turnover, associated with several Cretaceous OAEs. The absence of a species-level pattern may reflect species occurrence data that are too temporally coarse to robustly detect patterns. Five hypotheses of ecological selectivity relating anoxia to survivorship were tested across stage boundaries with respect to faunality, mobility, and diet using generalized linear models. Interestingly, benthic taxa were consistently selected against throughout the Cretaceous regardless of the presence or absence of OAEs. These results suggest that: (1) the Cenomanian/Turonian boundary (OAE2) was associated with a decline in molluscan diversity and increase in extinction rate that were significantly more severe than Cretaceous background levels; and (2) no differential ecological selectivity was associated with OAE-related diversity declines among the variables tested here.

Conodonts are among the first vertebrates to have evolved mineralized tooth-like structures. Among these, the so-called P1 elements are known to have been used to break down food and display a wide variety of morphologies. In particular, the repeated independent evolution of platform-like P1 elements have been suggested to correspond to similar functional constraints linked to diet. To test this hypothesis of convergence, we measured size (as element length) for various conodont taxa and used it as a proxy for trophic level. We then tested the correlation between size and platform presence/absence, both on raw data and in a phylogenetic context. Retaining or excluding the platform traits from the character matrix has limited impact on the resulting phylogeny. Contrary to platform presence/absence, size shows no phylogenetic signal. Using the raw data, size and platform presence appear positively correlated. That correlation, however, is no longer significant if one corrects for the phylogeny. We conclude that platform presence cannot be explained by an enlargement of the conodont element, be it via a trophic-level change or developmental constraints. This suggests that conodonts as a whole, and in particular platform-bearing conodonts, were an ecologically diverse group and that the various known platform types are likely to reflect different, rather than convergent, ecological niches.

It was once thought that the endemic carnivorous mammals of South America, the metatherian sparassodonts, were driven extinct by North American carnivorans through competitive exclusion. However, sparassodonts went extinct before most groups of carnivorans entered South America; only the endemic Cyonasua-group procyonids (Cyonasua and Chapalmalania), which immigrated to South America nearly 4 million years earlier than other carnivorans, significantly overlapped with sparassodonts in time. In this study, we examine the functional morphology of the dentition of Cyonasua and Chapalmalania through quantitative analysis to determine the dietary habits of these taxa and the degree to which they may have ecologically overlapped sparassodonts and large predatory Neogene didelphimorphians. We find Cyonasua and Chapalmalania to be more carnivorous than extant procyonids, other than Bassariscus, in agreement with previous studies, but more omnivorous than most other carnivorans and all meat-eating South American metatherians, including sparassodonts. The extreme ecological dissimilarity between Cyonasua-group procyonids and members of the endemic South American predator guild may explain why procyonids were able to successfully establish themselves in South America several million years earlier than most other northern mammals (including all other carnivorans): they moved into a previously unoccupied ecological niche (large omnivore) and avoided direct competition with incumbent native species, a situation similar to that documented in historical cases of biological invasion. The omnivorous diets and climbing/swimming abilities of procyonids may have increased their chances for a successful over-water dispersal relative to other carnivorans, further favoring their successful establishment in South America.

Nonmammalian cynodonts represent a speciose and ecologically diverse group with a fossil record stretching from the late Permian until the Cretaceous. Because of their role as major components of Triassic terrestrial ecosystems and as the direct ancestors of mammals, cynodonts are an important group for understanding Mesozoic tetrapod diversity. We examine patterns of nonmammalian cynodont species richness and the quality of their fossil record. A supertree of cynodonts is constructed from recently published trees and time calibrated using a Bayesian approach. While this approach pushes the root of Cynodontia back to the earliest Guadalupian, the origins of Cynognathia and Probainognathia are close to their first appearance in the fossil record. Taxic, subsampled, and phylogenetic diversity estimates support a major cynodont radiation following the end-Permian mass extinction, but conflicting signals are observed at the end of the Triassic. The taxic diversity estimate shows high diversity in the Rhaetian and a drop across the Triassic/Jurassic boundary, while the phylogenetic diversity indicates an earlier extinction between the Norian and Rhaetian. The difference is attributed to the prevalence of taxa based solely on teeth in the Rhaetian, which are not included in the phylogenetic diversity estimate. Examining the completeness of cynodont specimens through geological time does not support a decrease in preservation potential; although the median completeness score decreases in the Late Triassic, the range of values remains consistent. Instead, the poor completeness scores are attributed to a shift in sampling and taxonomic practices: an increased prevalence in microvertebrate sampling and the naming of fragmentary material.