Understanding cranial sutures and how they relieve and dissipate stress is essential to assess their role in cranial biomechanics and to develop highly accurate predictive models. This involves examining how ontogeny affects cranial sutures, as well as their morphology and function, and how these changes through time may impact essential biomechanical loadings such as chewing or direct biting. In this work, we study the cranial sutures of Crocodylus niloticus in detail using contact elements under finite element analysis. Contact elements permit the creation of a physical relationship between two bones that are in contact and even the configuration of these relationships, for example, in terms of movement or flexibility. The definition of bone contacts may require linear and/or nonlinear computational solutions to attain higher accuracy. Herein, skull geometry is tested to determine how bones may be altered by different types of contacts under various conditions. As predicted, the absence of sutures or cranial kinesis leads to a reduction in stress distribution across the skull, whereas sutures and cranial kinesis help the skull relieve stress and prevent certain bones from sustaining high stress levels. The type of contact used in individual sutures has a significant effect on model outcomes. Additionally, feeding behaviors significantly impact cranial biomechanics, reflecting the influence of other variables that may be applied to the models. As highlighted by the results, in order to obtain accurate results when analyzing fossil taxa, the nature of the cranial sutures should be taken into account. Therefore, developing predictive models based on living taxa is invaluable, because it facilitates the study of extinct taxa for which there is a lack of information on the fibrous joints due to poor or no preservation in the fossil record.

Phylogenetic inferences using combined datasets of both extant and extinct species have grown increasingly popular, in part thanks to the development of the fossilized birth–death (FBD) process. The FBD process provides a powerful model for the evolution of past and present lineages and can be used for both inference and simulation. Simulations in particular are very helpful for new users to gain better understanding of the model and its different components. In this work, we present FossilSimShiny, a visual application for simulating phylogenies, fossil samples, and fossil taxonomies under the FBD process. The app integrates a wide range of simulation models and presents the simulation results in clear, customizable figures. As a teaching tool, FossilSimShiny allows lecturers to create illustration plots and students to directly experiment with the model. For research applications, the app can help researchers save time and effort by testing and calibrating simulation setups before running them on a large scale.

Tribrachidium heraldicum is an Ediacaran body fossil characterized by triradial symmetry. Previous work has suggested that the anatomy of Tribrachidium was conducive to passive suspension feeding; however, these analyses used an inaccurate model and a relatively simple set of simulations. Using computational fluid dynamics, we explore the functional morphology of Tribrachidium in unprecedented detail by gauging how the presence or absence of distinctive anatomical features (e.g., apical pits and arms) affects flow patterns. Additionally, we map particle pathways, quantify deposition rates at proposed feeding sites, and assess gregarious feeding habits to more fully reconstruct the lifestyle of this enigmatic taxon. Our results provide strong support for interpreting Tribrachidium as a macroscopic suspension feeder, with the apical pits representing loci of particle collection (and possibly ingestion) and the triradial arms representing morphological adaptations for interrupting flow and inducing settling. More speculatively, we suggest that the radial grooves may represent ciliated pathways through which food particles accumulating in the wake of the organism were transported toward the apical pits. Finally, our results allow us to generate new functional hypotheses for other Ediacaran taxa with a triradial body plan. This work refines our understanding of the appearance of suspension feeding in shallow-water paleoenvironments, with implications for the radiation of Metazoa across the Ediacaran/Cambrian boundary.

The trend of global cooling across the Cenozoic transformed the North American landscape from closed forest to more open grasslands, resulting in dietary adaptations in herbivores in response to shifting resources. In contrast, the material properties of the predator food source (muscle, skin, and bone) have remained constant over this transition, suggesting a corresponding lack of change in predator dietary adaptations. We investigated the North American mammalian predator fossil record using a tooth-shape metric and body mass, predicting that the former would exhibit stability. Instead, we found that mean molar morphology became more blade-like, with our tooth-shape metric sharply increasing in the late Eocene and remaining high from the Oligocene onward. Subsequent tests in extant carnivorans reveal taxa with more bladelike teeth are prevalent in more open environments. Our results reveal an unexpected functional shift among North American predators in response to large-scale environmental changes across the Cenozoic.

The Fezouata Shale Formation has dramatically impacted our understanding of Early Ordovician marine ecosystems before the great Ordovician biodiversification event (GOBE), thanks to the abundance and quality of exceptionally preserved animals within it. Systematic work has noted that the shelly fossil subassemblages of the Fezouata Shale biota are typical of open-marine deposits from the Lower Ordovician, but no studies have tested the quantitative validity of this statement. We extracted 491 occurrences of recalcitrant fossil genera from the Paleobiology Database to reconstruct 31 subassemblages to explore the paleoecology of the Fezouata Shale and other contemporary, high-latitude (66°S–90°S) deposits from the Lower Ordovician (485.4–470 Ma) and test the interpretation that the Fezouata Shale biota is typical for an Ordovician open-marine environment. Sørensen's dissimilarity metrics and Wilcoxon tests indicate that the subassemblages of the Tremadocian-aged lower Fezouata Shale are approximately 20% more heterogenous than the Floian-aged upper Fezouata Shale. Dissimilarity metrics and visualization suggest that while the lower Fezouata and upper Fezouata share faunal components, the two sections have distinct faunas. We find that the faunal composition of the lower Fezouata Shale is comparable with other Tremadocian-aged subassemblages from high latitudes, suggesting that it is typical for an Early Ordovician open-marine environment. We also find differences in faunal composition between Tremadocian- and Floian-aged deposits. Our results corroborate previous field-based and qualitative systematic studies that concluded that the shelly assemblages of the Fezouata Shale are comparable with those of other Lower Ordovician deposits from high latitudes. This establishes the first quantitative baseline for examining the composition and variability within the assemblages of the Fezouata Shale and will be key to future studies attempting to discern the degree to which it can inform our understanding of marine ecosystems just before the start of the GOBE.

Humans have dramatically transformed ecosystems over the previous millennia and are potentially causing a mass extinction event comparable to the others that shaped the history of life. However, only a fraction of these impacts has been directly recorded, limiting conservation actions. Conservation paleobiology leverages geohistorical records to offer a long-term perspective on biodiversity change in the face of anthropogenic stressors. Nevertheless, the field's on-the-ground contributions to conservation outcomes are still developing. Here, we present an overview of directions in which paleobiological research could progress to aid conservation in the coming decades using elasmobranchs (sharks, rays, and skates)—a highly threatened group with a rich fossil record—as a model. These research directions are guided by areas of overlap between an expert-led list of current elasmobranch conservation priorities and available fossil and historical records. Four research topics emerged for which paleobiological research could address open questions in elasmobranch science and conservation: (1) baselines, (2) ecological roles, (3) threats, and (4) conservation priorities. Increasingly rich datasets and novel analytical frameworks present exciting opportunities to apply the elasmobranch fossil record to conservation practice. A similar approach could be extended to other clades. Given the synthetic nature of these research topics, we encourage collaboration across timescales and with conservation practitioners to safeguard the future of our planet's rapidly disappearing species.

It is a pleasure to be invited to contribute to the celebration of the publication of “punctuated equilibria” (“punk eek”) 50 years ago—the canonical version I did with Steve Gould (Eldredge and Gould 1972) at the behest of Tom Schopf for his visionary project to inject more thought, more interpretation and theory, into the working lives of paleontologists.

A traditional typological approach to taxonomy often does not adequately account for intraspecific variation and can result in taxonomic oversplitting. For many groups, including ammonoids of the Placenticeras genus, intraspecific variation documented in recent studies (e.g., ontogenetic changes, sexual dimorphism, polymorphism) challenges the historic proliferation of species names. Here, we used a population approach to taxonomy and quantitatively evaluated morphometric variation in a sample of Late Cretaceous (Santonian–Campanian) Placenticeras from Alabama and adjacent counties.

We used linear mixed models (LMMs) to characterize how morphological variables scale with conch size across the sample, exploiting mixed longitudinal data to evaluate individual variation in growth and inform interpretations of multivariate analyses. Extended LMMs incorporating geological formation evaluated morphological changes through time. Principal component analysis and clustering analysis were then used to evaluate the number of distinct clusters that emerged in multivariate morphospace independent of previous taxon name assignments.

Discontinuous scaling relationships and distinct clusters in multivariate space suggest sexual dimorphism characterized by differences in adult size and, secondarily, shape. Previous Stantonoceras and Placenticeras assignments broadly overlap in our morphospace, failing to justify this historic distinction (as sexual dimorphs or as genera or subgenera). Placenticeras conch morphology and ornament placement changed through time, suggesting a potential utility for coarse (stage-level) biostratigraphy. However, temporal changes were not associated with distinct clusters in morphospace, and our data fail to support the plethora of reported species names. As few as one or two (successive) species may be present in our sample (representing 130 years of collection effort). In addition to highlighting the need for a significant taxonomic revision of the Placenticeras genus, this study demonstrates the utility of LMMs for distinguishing between different sources of morphological variation, improving interpretations of morphospace under a population approach to taxonomy, and maximizing the amount of ontogenetic information that can be obtained nondestructively.

Theropods are obligate bipedal dinosaurs that appeared 230 Ma and are still extant as birds. Their history is characterized by extreme variations in body mass, with gigantism evolving convergently between many lineages. However, no quantification of hindlimb functional morphology has shown whether these body mass increases led to similar specializations between distinct lineages. Here we studied femoral shape variation across 41 species of theropods (n = 68 specimens) using a high-density 3D geometric morphometric approach. We demonstrated that the heaviest theropods evolved wider epiphyses and a more distally located fourth trochanter, as previously demonstrated in early archosaurs, along with an upturned femoral head and a mediodistal crest that extended proximally along the shaft. Phylogenetically informed analyses highlighted that these traits evolved convergently within six major theropod lineages, regardless of their maximum body mass. Conversely, the most gracile femora were distinct from the rest of the dataset, which we interpret as a femoral specialization to “miniaturization” evolving close to Avialae (bird lineage). Our results support a gradual evolution of known “avian” features, such as the fusion between lesser and greater trochanters and a reduction of the epiphyseal offset, independent from body mass variations, which may relate to a more “avian” type of locomotion (more knee than hip driven). The distinction between body mass variations and a more “avian” locomotion is represented by a decoupling in the mediodistal crest morphology, whose biomechanical nature should be studied to better understand the importance of its functional role in gigantism, miniaturization, and higher parasagittal abilities.

Scaling fluctuation analyses of marine animal diversity and extinction and origination rates based on the Paleobiology Database occurrence data have opened new perspectives on macroevolution, supporting the hypothesis that the environment (climate proxies) and life (extinction and origination rates) are scaling over the “megaclimate” biogeological regime (from ≈1 Myr to at least 400 Myr). In the emerging picture, biodiversity is a scaling “crossover” phenomenon being dominated by the environment at short timescales and by life at long timescales with a crossover at ≈40 Myr. These findings provide the empirical basis for constructing the Fractional MacroEvolution Model (FMEM), a simple stochastic model combining destabilizing and stabilizing tendencies in macroevolutionary dynamics, driven by two scaling processes: temperature and turnover rates.

Macroevolution models are typically deterministic (albeit sometimes perturbed by random noises) and are based on integer-ordered differential equations. In contrast, the FMEM is stochastic and based on fractional-ordered equations. Stochastic models are natural for systems with large numbers of degrees of freedom, and fractional equations naturally give rise to scaling processes.

The basic FMEM drivers are fractional Brownian motions (temperature, T) and fractional Gaussian noises (turnover rates, E+) and the responses (solutions), are fractionally integrated fractional relaxation noises (diversity [D], extinction [E], origination [O], and E− = O − E). We discuss the impulse response (itself representing the model response to a bolide impact) and derive the model's full statistical properties. By numerically solving the model, we verified the mathematical analysis and compared both uniformly and irregularly sampled model outputs with paleobiology series.

Apicobasal ridges are longitudinal ridges of enamel that are particularly common in several clades of aquatic-feeding predatory amniotes, including Plesiosauria, Ichthyosauria, Mosasauridae, Crocodylia, and Spinosauridae, as well as some early members of Cetacea. Although the repeated evolution of these dental ridges in unrelated clades suggests an adaptive benefit, their primary function in feeding is debated. Hypothesized functions range from increasing tooth strength to improving prey puncture or removal efficiency, but these have never been quantitatively tested. This study utilizes finite element analysis (FEA) to assess the impact of apicobasal ridges upon tooth crown strength in aquatic-feeding amniotes. Drawing on morphometric data from fossilized tooth crowns, a set of digital models was constructed to calculate the performance of smooth and ridged tooth variants under simulated bite force loadings. The similarities in overall stress distribution patterns across models of the same tooth shape, regardless of the presence or morphology of ridges, indicate that apicobasal ridges have little impact on stress reduction within the tooth crown. Ultimately, these findings suggest that apicobasal ridges have a minimal role in improving crown strength and form a framework for future research into the remaining hypotheses.

Crocodylomorpha is a large and diverse clade with a long evolutionary history now restricted to modern crocodilians. Tethysuchia is a less-inclusive clade of semi-amphibious taxa that crossed two biological crises: the second Oceanic Anoxic Event (OAE 2) and the Cretaceous/Paleogene (K/Pg) crisis. Numerous studies have sought to find the driving factors explaining crocodylomorph evolution, producing contradictory conclusions. Studies of included groups may be useful. Here, we study factors driving tethysuchian evolution using phylogenetically informed statistical analyses. First, we tested the phylogenetic structure of tethysuchian extinction at the OAE 2 and K/Pg crises. We then used phylogenetic comparative methods to test the influence of intrinsic (body size, snout proportion) and extrinsic (temperature, paleolatitude) factors on the evolution of tethysuchian diversity at the OAE 2 and the K/Pg crises. Finally, we tested whether temperature influenced the evolution of body size. We conclude that (1) extinction was not random in regard to phylogeny for Tethysuchia at the OAE 2 and K/Pg crises; (2) while an important tethysuchian turnover follows OAE 2, the K/Pg crisis was followed by an explosion in diversity of tethysuchians, probably linked to the colonization of emptied ecological niches; (3) tethysuchians lived in warmer environments after the OAE 2 crisis, possibly because of both global warming and latitudinal distribution shifts; (4) there is a significant change of snout proportion after the OAE 2 and the K/Pg crises, likely caused by niche partitioning; and (5) there is a positive correlation between body size and temperature, possibly because of a longer growth season.

The fitness of groups is often considered to be the average fitness among constituent members. This assumption has been useful for developing models of multilevel selection, but its uncritical adoption has held back our understanding of how multilevel selection actually works in nature. If group fitness is only equal to mean member fitness, then it is a simple task to erode the importance of group-level selection in all multilevel scenarios—species selection could then be reduced to organismal selection as easily as group selection can. Because selection from different levels can act on a single trait, body size, for example, there must be a way to translate one level of fitness to another. This allows the calculation of the contributions of selection at each level. If high-level fitness is not a simple function of member fitness, then how do they interlace? Here we reintroduce Leigh Van Valen’s argument for the inclusion of expansion as a component of fitness. We show that expansion is an integral part of fitness even if one does not subscribe to the energetic view of fitness from which Van Valen originally derived it. From a hierarchical perspective, expansion is the projection of demographic fitness from one level to the next level up; differential births and deaths at one level produce differential expansion one level above. Including expansion in our conceptual tool kit helps allay concerns about our ability to identify the level of selection using a number of methods as well as allowing for the various forms of multilevel selection to be seen as manifestations of the same basic process.

The Paleozoic represents a key time interval in the origins and early diversification of chondrichthyans (cartilaginous fishes), but their diversity and macroevolution are largely obscured by heterogenous spatial and temporal sampling. The predominantly cartilaginous skeletons of chondrichthyans pose an additional limitation on their preservation potential and hence on the quality of their fossil record. Here, we use a newly compiled genus-level dataset and the application of sampling standardization methods to analyze global total-chondrichthyan diversity dynamics through time from their first appearance in the Ordovician through to the end of the Permian. Subsampled estimates of chondrichthyan genus richness were initially low in the Ordovician and Silurian but increased substantially in the Early Devonian. Richness reached its maximum in the middle Carboniferous before dropping across the Carboniferous/Permian boundary and gradually decreasing throughout the Permian. Sampling is higher in both the Devonian and Carboniferous compared with the Silurian and most of the Permian stages. Shark-like scales from the Ordovician are too limited to allow for some of the subsampling techniques. Our results detect two Paleozoic radiations in chondrichthyan diversity: the first in the earliest Devonian, led by acanthodians (stem-group chondrichthyans), which then decline rapidly by the Late Devonian, and the second in the earliest Carboniferous, led by holocephalans, which increase greatly in richness across the Devonian/Carboniferous boundary. Dispersal of chondrichthyans, specifically holocephalans, into deeper-water environments may reflect a niche expansion following the faunal displacement in the aftermath of the Hangenberg extinction event at the end of the Devonian.

Stratigraphic paleobiology uses a modern understanding of the construction of the stratigraphic record—from beds to depositional sequences to sedimentary basins—to interpret patterns and guide sampling strategies in the fossil record. Over the past 25 years, its principles have been established primarily through forward numerical modeling, originally in shallow-marine systems and more recently in nonmarine systems. Predictions of these models have been tested through outcrop-scale and basin-scale field studies, which have also revealed new insights. At multi-basin and global scales, understanding the joint development of the biotic and sedimentary records has come largely from macrostratigraphy, the analysis of gap-bound packages of sedimentary rock. Here, we present recent advances in six major areas of stratigraphic paleobiology, including critical tests in the Po Plain of Italy, mass extinctions and recoveries, contrasts of shallow-marine and nonmarine systems, the interrelationships of habitats and stratigraphic architecture, large-scale stratigraphic architecture, and the assembly of regional ecosystems. We highlight the potential for future research that applies stratigraphic paleobiological concepts to studies of climate change, geochemistry, phylogenetics, and the large-scale structure of the fossil record. We conclude with the need for more stratigraphic thinking in paleobiology.

The fossil record is spatiotemporally heterogeneous: taxon occurrence data have patchy spatial distributions, and this patchiness varies through time. Large-scale quantitative paleobiology studies that fail to account for heterogeneous sampling coverage will generate uninformative inferences at best and confidently draw wrong conclusions at worst. Explicitly spatial methods of standardization are necessary for analyses of large-scale fossil datasets, because nonspatial sample standardization, such as diversity rarefaction, is insufficient to reduce the signal of varying spatial coverage through time or between environments and clades. Spatial standardization should control both geographic area and dispersion (spread) of fossil localities. In addition to standardizing the spatial distribution of data, other factors may be standardized, including environmental heterogeneity or the number of publications or field collecting units that report taxon occurrences. Using a case study of published global Paleobiology Database occurrences, we demonstrate strong signals of sampling; without spatial standardization, these sampling signatures could be misattributed to biological processes. We discuss practical issues of implementing spatial standardization via subsampling and present the new R package divvy to improve the accessibility of spatial analysis. The software provides three spatial subsampling approaches, as well as related tools to quantify spatial coverage. After reviewing the theory, practice, and history of equalizing spatial coverage between data comparison groups, we outline priority areas to improve related data collection, analysis, and reporting practices in paleobiology.

The process of evolution and the structures it produces are best understood in the light of hierarchy theory. The biota traditionally is described by either the genealogical Linnaean hierarchy or economic hierarchies of communities or ecosystems. Here we describe the Bretskyan hierarchy—a hybrid eco-genealogical hierarchy that consists of nested sets of different-sized, usually polyphyletic communities of interacting individuals separated from other such communities in space and time at multiple scales. The Bretskyan hierarchy consists of elements that have both genealogical and economic properties and functions—situated between, and connecting the elements of, the economic hierarchies (Vernadskyan) and the genealogical (Linnaean) hierarchy. The described hierarchy at lower tiers is populated by holobionts, individuals composed of multiple polyphyletic lineages integrated by functional interactions or biotically fabricated structures, such as membranes. At larger spatial tiers and longer time scales, the members of the Bretskyan hierarchy are of a more diffuse nature, partially due to the small size and relatively short duration of us as observers of larger and longer-lasting structures, here described as geobiomes. Their individuality is externally forced and directly tied to the spatial and temporal physical structures of our planet. These are sub-bioprovinces and bioprovinces—large and effectively isolated spatiotemporal structures of biota integrated internally by coevolution and individuated externally by a hierarchy of barriers. Gaia is here understood as the largest eco-genealogical individual compartmentalized by the outer space of the Earth and integrated at long time scales by biotic interactions and plate tectonic mixing of biota. The existence of a hierarchy of barriers and multilevel allopatry suggests that geographic isolation takes part not only in individuating species lineages, but also in producing coherent complexes of separate lineages forming bioprovinces at multiple space and time scales. The sizes, configurations, and durations of Bretskyan units are directly tied to geodynamics, demonstrating the central role of the physical planet in the processes of individuation and merging of geobiomes and the control of coevolution, and all its ramifications, at multiple space and time scales. The Bretskyan hierarchy also allows the integration of previously unconnected themes—“egalitarian” major transitions in individuality (e.g., eukaryogenesis) and some of the megatrajectories in the history of life—into a single theoretical framework of spatial and temporal scaling of eco-genealogy. The pervasive scaling of geodynamical processes and the direct connection of geodynamics to the dynamics of Bretskyan units allows us to formulate conjectures on the scales and limits of spatial and temporal contingency and competitiveness of biotas in evolution.

We know the fossil record is incomplete, but just how much biodiversity does it miss? We produce the first geographically controlled estimate by comparing the geographic ranges of 34,266 modern tetrapods with a map of the world's sedimentary basins. By modeling which tetrapods live within sedimentary basins, we produce a first-order estimate of what might be found in the fossil record of the future. In this record, nearly 30% of tetrapod species have almost no chance of fossilizing, and more stringent criteria for fossilization exclude far more diversity. This geographically structured fossil record preserves disparate patterns of taxonomic and phylogenetic diversity in different tetrapod groups and underpreserves projected extinctions. For the globally threatened amphibians, the magnitude of the extinction of all endangered species would be underestimated by 66–98% in our future record. These results raise profound questions about the structure of the fossil record. Is it capable of recording major origination and extinction events on land? Have swaths of terrestrial diversity gone unrecorded based on geography alone? There are chapters of Earth history that paleontologists can never hope to know, but what is missing, and why?

During the late Miocene and early Pliocene (latest Hemphillian, Hh4 interval, 5.7 to 4.75 Ma) a distinctive suite of four species of extinct horses (Family Equidae) were widespread in North America. This includes Nannippus aztecus, Neohipparion eurystyle, Astrohippus stocki, and Dinohippus mexicanus. In Florida, two additional equid species, Pseudhipparion simpsoni and Cormohipparion emsliei, are also typically found at Hh4 localities. Here we compare horses from four Hh4 Florida fossil sites, including three from the Bone Valley mines and a fourth from the recently discovered Montbrook site. Two of these sites have all six expected species, one has five species, and one has only four species present. To explain these differences, we used species counts from research databases and rarefaction simulation to clarify the relative abundances, species richness, and occurrences of these horses from these four sites. The Palmetto Mine (Agrico) site, with five equid species, appears to lack the sixth species owing to ecological reasons. This is different from Montbrook, the site with only four of the six species. Results indicate that Montbrook is likely lacking two missing equid species for a couple of reasons: sampling bias and biological/ecological causes. Our results demonstrate that sampling biases can account for observed equid species richness when the overall abundance of certain equid species is low. Nevertheless, other factors, including ecology and with sufficient resolution, perhaps also time, may also explain the distribution and occurrences of individual species at these and other fossil sites. In a broader perspective, analyses such as this example provide an opportunity to address a persistent challenge in paleontology, that is, how do we explain absences of extinct taxa from the fossil record?

Data compilations expand the scope of research; however, data citation practice lags behind advances in data use. It remains uncommon for data users to credit data producers in professionally meaningful ways. In paleontology, databases like the Paleobiology Database (PBDB) enable assessment of patterns and processes spanning millions of years, up to global scale. The status quo for data citation creates an imbalance wherein publications drawing data from the PBDB receive significantly more citations (median: 4.3 ± 3.5 citations/year) than the publications producing the data (1.4 ± 1.3 citations/year). By accounting for data reuse where citations were neglected, the projected citation rate for data-provisioning publications approached parity (4.2 ± 2.2 citations/year) and the impact factor of paleontological journals (n = 55) increased by an average of 13.4% (maximum increase = 57.8%) in 2019. Without rebalancing the distribution of scientific credit, emerging “big data” research in paleontology—and science in general—is at risk of undercutting itself through a systematic devaluation of the work that is foundational to the discipline.