Paleobiology

Evolving interactions between predators and prey constitute one of the major adaptive influences on marine animals during the Paleozoic. Crinoids and fish constitute a predator–prey system that may date back to at least the Silurian, as suggested by patterns of crinoid regeneration and spinosity in concert with changes in the predatory fauna. Here we present data on the frequency of breakage and regeneration in the spines of the Middle Devonian camerate Gennaeocrinus and late Paleozoic cladids, as well as an expanded survey of the prevalence of spinosity and infestation by platyceratid gastropods on crinoid genera during the Paleozoic. Spine regeneration frequency in the measured populations is comparable to arm regeneration frequencies from Mississippian Rhodocrinites and from modern deep-water crinoid populations. The prevalence of spinosity varies by taxon, time, and anatomy among Paleozoic crinoids; notably, spinosity in camerates increased from the Silurian through the Mississippian and decreased sharply during the Pennsylvanian, whereas spines were uncommon in cladids until their Late Mississippian diversification. Among camerates, tegmen spinosity is positively correlated with the presence of infesting platyceratid gastropods. These results allow us to evaluate several hypotheses for the effects of predation on morphological differences between early, middle, and late Paleozoic crinoid faunas. Our data corroborate the hypothesis that predators targeted epibionts on camerate crinoids and anal sacs on advanced cladids and suggest that the replacement of shearing predators by crushing predators after the Hangenberg extinction affected the locations of spines in Mississippian camerates.

Despite a global fossil record, Metatheria are now largely restricted to Australasia and South America. Most metatherian paleodiversity studies to date are limited to particular subclades, time intervals, and/or regions, and few consider uneven sampling. Here, we present a comprehensive new data set on metatherian fossil occurrences (Barremian to end Pliocene). These data are analyzed using standard rarefaction and shareholder quorum subsampling (including a new protocol for handling Lagerstätte-like localities).

Global metatherian diversity was lowest during the Cretaceous, and increased sharply in the Paleocene, when the South American record begins. Global and South American diversity rose in the early Eocene then fell in the late Eocene, in contrast to the North American pattern. In the Oligocene, diversity declined in the Americas, but this was more than offset by Oligocene radiations in Australia. Diversity continued to decrease in Laurasia, with final representatives in North America (excluding the later entry of Didelphis virginiana) and Europe in the early Miocene, and Asia in the middle Miocene. Global metatherian diversity appears to have peaked in the early Miocene, especially in Australia. Following a trough in the late Miocene, the Pliocene saw another increase in global diversity. By this time, metatherian biogeographic distribution had essentially contracted to that of today.

Comparison of the raw and sampling-corrected diversity estimates, coupled with evaluation of “coverage” and number of prolific sites, demonstrates that the metatherian fossil record is spatially and temporally extremely patchy. Therefore, assessments of macroevolutionary patterns based on the raw fossil record (as in most previous studies) are inadvisable.

For mammals today, mountains are diverse ecosystems globally, yet the strong relationship between species richness and topographic complexity is not a persistent feature of the fossil record. Based on fossil-occurrence data, diversity and diversification rates in the intermontane western North America varied through time, increasing significantly during an interval of global warming and regional intensification of tectonic activity from 18 to 14 Ma. However, our ability to infer origination and extinction rates reliably from the fossil record is affected by variation in preservation history. To investigate the influence of preservation on estimates of diversification rates, I simulated fossil records under four alternative diversification hypotheses and six preservation scenarios. Diversification hypotheses included tectonically controlled speciation pulses, while preservation scenarios were based on common trends (e.g., increasing rock record toward the present) or derived from fossil occurrences and the continental rock record. For each scenario, I estimated origination, extinction, and diversification rates using three standard methods—per capita, three-timer, and capture–mark–recapture (CMR) metrics—and evaluated the ability of the simulated fossil records to accurately recover the underlying diversification dynamics. Despite variable and low preservation probabilities, simulated fossil records retained the signal of true rates in several of the scenarios. The three metrics did not exhibit similar behavior under each preservation scenario: while three-timer and CMR metrics produced more accurate rate estimates, per capita rates tended to better reproduce true shifts in origination rates. All metrics suffered from spurious peaks in origination and extinction rates when highly volatile preservation impacted the simulated record. Results from these simulations indicate that elevated diversification rates in relation to tectonic activity during the middle Miocene are likely to be evident in the fossil record, even if preservation in the North American fossil record was variable. Input from the past is necessary to evaluate the ultimate mechanisms underlying speciation and extinction dynamics.

Previous attempts to quantify the test-flattening trend in Heterostegina depressa with water depth have been rather unsuccessful. Due to its broad depth distribution, H. depressa is a perfect model species to calibrate test flattening as a bathymetric signal for fossil assemblages. This might enable us to better reconstruct paleoenvironments of fossil communities of larger foraminifera or even provide clues to the degree of transport in allochthonous deposits. In this study, we used growth-independent functions to describe the change of test thickness throughout ontogeny. Four growth-invariant characters, deriving from these functions, clearly quantify a transition of individuals with thicker to thinner central parts along the water-depth gradient. This transition is probably controlled by light intensity, because the photosymbionts of H. depressa (diatoms) are most effective at low irradiation levels. Thus, specimens at shallower depths grow thicker to reduce light penetration, whereas specimens living deeper than the light optimum increase their surface by flattening to obtain better exposure to light.

There is a well-established discrepancy between paleontological and molecular data regarding the timing of the origin and diversification of placental mammals. Molecular estimates place interordinal diversification dates in the Cretaceous, while no unambiguous crown placental fossils have been found prior to the end-Cretaceous mass extinction. Here, the completeness of the eutherian fossil record through geological time is evaluated to assess the suggestion that a poor fossil record is largely responsible for the difference in estimates of placental origins. The completeness of fossil specimens was measured using the character completeness metric, which quantifies the completeness of fossil taxa as the percentage of phylogenetic characters available to be scored for any given taxon. Our data set comprised 33 published cladistic matrices representing 445 genera, of which 333 were coded at the species level.

There was no significant difference in eutherian completeness across the Cretaceous/Paleogene (K/Pg) boundary. This suggests that the lack of placental mammal fossils in the Cretaceous is not due to a poor fossil record but more likely represents a genuine absence of placental mammals in the Cretaceous. This result supports the “explosive model” of early placental evolution, whereby placental mammals originated around the time of the K/Pg boundary and diversified soon after.

No correlation was found between the completeness pattern observed in this study and those of previous completeness studies on birds and sauropodomorph dinosaurs, suggesting that different factors affect the preservation of these groups. No correlations were found with various isotope proxy measures, but Akaike information criterion analysis found that eutherian character completeness metric scores were best explained by models involving the marine-carbonate strontium-isotope ratios (87Sr/86Sr), suggesting that tectonic activity might play a role in controlling the completeness of the eutherian fossil record.

Sauropodomorpha included the largest known terrestrial vertebrates and was the first dinosaur clade to achieve a global distribution. This success is associated with their early adoption of herbivory, and sauropod gigantism has been hypothesized to be a specialization for bulk feeding and obligate high-fiber herbivory. Here, we apply a combination of biomechanical character analysis and comparative phylogenetic methods with the aim of quantifying the evolutionary mechanics of the sauropodomorph feeding apparatus. We test for the role of convergence to common feeding function and divergence toward functional optima across sauropodomorph evolution, quantify the rate of evolution for functional characters, and test for coincident evolutionary rate shifts in craniodental functional characters and body mass. Results identify a functional shift toward increased cranial robustness, increased bite force, and the onset of static occlusion at the base of the Sauropoda, consistent with a shift toward bulk feeding. Trends toward similarity in functional characters are observed in Diplodocoidea and Titanosauriformes. However, diplodocids and titanosaurs retain significant craniodental functional differences, and evidence for convergent adoption of a common “adaptive zone” between them is weak. Modeling of craniodental character and body-mass evolution demonstrates that these functional shifts were not correlated with evolutionary rate shifts. Instead, a significant correlation between body mass and characters related to bite force and cranial robustness suggests a correlated-progression evolutionary mode, with positive-feedback loops between body mass and dietary specializations fueling sauropod gigantism.

Throughout the late Quaternary, the Sahul (Pleistocene Australia–New Guinea) vertebrate fauna was dominated by a diversity of large mammals, birds, and reptiles, commonly referred to as megafauna. Since ca. 450–400Ka, approximately 88 species disappeared in Sahul, including kangaroos exceeding 200kg in size, wombat-like animals the size of hippopotamuses, flightless birds, and giant monitor lizards that were likely venomous. Ongoing debates over the primary cause of these extinctions have typically favored climate change or human activities. Improving our understanding of the population biology of extinct megafauna as more refined paleoenvironmental data sets become available will assist in identifying their potential vulnerabilities. Here, we apply a multiproxy approach to analyze fossil teeth from deposits dated to the middle and late Pleistocene at Cuddie Springs in southeastern Australia, assessing relative aridity via oxygen isotopes as well as vegetation and megafaunal diets using both carbon isotopes and dental microwear texture analyses. We report that the Cuddie Springs middle Pleistocene fauna was largely dominated by browsers, including consumers of C4 shrubs, but that by late Pleistocene times the C4 dietary component was markedly reduced. Our results suggest dietary restriction in more arid conditions. These dietary shifts are consistent with other independently derived isotopic data from eggshells and wombat teeth that also suggest a reduction in C4 vegetation after ~45 Ka in southeastern Australia, coincident with increasing aridification through the middle to late Pleistocene. Understanding the ecology of extinct species is important in clarifying the primary drivers of faunal extinction in Sahul. The results presented here highlight the potential impacts of aridification on marsupial megafauna. The trend to increasingly arid conditions through the middle to late Pleistocene (as identified in other paleoenvironmental records and now also observed, in part, in the Cuddie Springs sequence) may have stressed the most vulnerable animals, perhaps accelerating the decline of late Pleistocene megafauna in Australia.

Morphological responses of nonmammalian herbivores to external ecological drivers have not been quantified over extended timescales. Herbivorous nonavian dinosaurs are an ideal group to test for such responses, because they dominated terrestrial ecosystems for more than 155 Myr and included the largest herbivores that ever existed. The radiation of dinosaurs was punctuated by several ecologically important events, including extinctions at the Triassic/Jurassic (Tr/J) and Jurassic/Cretaceous (J/K) boundaries, the decline of cycadophytes, and the origin of angiosperms, all of which may have had profound consequences for herbivore communities. Here we present the first analysis of morphological and biomechanical disparity for sauropodomorph and ornithischian dinosaurs in order to investigate patterns of jaw shape and function through time. We find that morphological and biomechanical mandibular disparity are decoupled: mandibular shape disparity follows taxonomic diversity, with a steady increase through the Mesozoic. By contrast, biomechanical disparity builds to a peak in the Late Jurassic that corresponds to increased functional variation among sauropods. The reduction in biomechanical disparity following this peak coincides with the J/K extinction, the associated loss of sauropod and stegosaur diversity, and the decline of cycadophytes. We find no specific correspondence between biomechanical disparity and the proliferation of angiosperms. Continual ecological and functional replacement of pre-existing taxa accounts for disparity patterns through much of the Cretaceous, with the exception of several unique groups, such as psittacosaurids that are never replaced in their biomechanical or morphological profiles.

Extant species of the genus Equus (e.g., horses, asses, and zebras) have a widespread distribution today on all continents except Antarctica. Extinct species of Equus represented by fossils were likewise widely distributed in the Pliocene and even more so during the Pleistocene. In order to understand the efficacy of “big data” for (paleo)biogeographic analyses, location records (latitude, longitude) and fossil occurrences for the genus Equus were mined and further explored from six databases, including iDigBio, Paleobiology Database, VertNet, BISON, Neotoma, and GBIF. These were chosen from a priori knowledge of where relevant data might be aggregated. We also realized that these databases have different objectives and data sources and therefore would provide a useful comparative study of the widespread taxon Equus in space and time.

The mining of Equus data from these six sources yielded a combined total of 123.8 K location records, including 116.2K fossil specimens. These include individual points that are unique, that is, only occurring in one of these databases, and those that are duplicated in multiple databases. Of the six databases, three (iDigBio, Paleobiology Database, and GBIF) were judged to be the most useful in the Equus use case. Most of the databases are biased toward North American records, thus limiting the reconstruction of the actual distribution of the genus Equus in space and time outside of this continent. Although Equus has a large number of digitally accessible records, fundamentally interesting questions pertaining to evolutionary dynamics and extinction geography are still a challenge for these kinds of biodiversity databases due primarily to the lack of sufficiently dense and precise temporal data.

Mesozoic marine ecosystems were dominated by several clades of reptiles, including sauropterygians, ichthyosaurs, crocodylomorphs, turtles, and mosasaurs, that repeatedly invaded ocean ecosystems. Previous research has shown that marine reptiles achieved great taxonomic diversity in the Middle Triassic, as they broadly diversified into many feeding modes in the aftermath of the Permo-Triassic mass extinction, but it is not known whether this initial phase of evolution was exceptional in the context of the entire Mesozoic. Here, we use a broad array of disparity, morphospace, and comparative phylogenetic analyses to test this. Metrics of ecomorphology, including functional disparity in the jaws and dentition and skull-size diversity, show that the Middle to early Late Triassic represented a time of pronounced phenotypic diversification in marine reptile evolution. Following the Late Triassic extinctions, diversity recovered, but disparity did not, and it took over 100 Myr for comparable variation to recover in the Campanian and Maastrichtian. Jurassic marine reptiles generally failed to radiate into vacated functional roles. The signatures of adaptive radiation are not seen in all marine reptile groups. Clades that diversified during the Triassic biotic recovery, the sauropterygians and ichthyosauromorphs, do show early diversifications, early high disparity, and early burst, while less support for these models is found in thalattosuchian crocodylomorphs and mosasaurs. Overall, the Triassic represented a special interval in marine reptile evolution, as a number of groups radiated into new adaptive zones.

The Paleobiology Database (PBDB; https://paleobiodb.org) consists of geographically and temporally explicit, taxonomically identified fossil occurrence data. The taxonomy utilized by the PBDB is not static, but is instead dynamically generated using an algorithm applied to separately managed taxonomic authority and opinion data. The PBDB owes its existence to many individuals, some of whom have entered more than 1.26 million fossil occurrences and over 570,000 taxonomic opinions, and some of whom have developed and maintained supporting infrastructure and analysis tools. Here, we provide an overview of the data model currently used by the PBDB and then briefly describe how this model is exposed via an Application Programming Interface (API). Our objective is to outline how PBDB data can now be accessed within individual scientific workflows, used to develop independently managed educational and scientific applications, and accessed to forge dynamic, near real-time connections to other data resources.

Detailed quantitative data has previously been collected from plant megafossil assemblages from a Middle Jurassic (Aalenian) plant bed from Hasty Bank, North Yorkshire, UK. We conducted a similar analysis of palynological dispersed sporomorph (spore and pollen) assemblages collected from the same section using the same sampling regime: 67 sporomorph taxa were recorded from 50 samples taken at 10 cm intervals through the plant bed. Basic palynofacies analysis was also undertaken on each sample. Both dispersed sporomorph and plant megafossil assemblages display consistent changes in composition, diversity (richness), and abundance through time. However, the dispersed sporomorph and plant megafossil records provide conflicting evidence for the nature of parent vegetation. Specifically, conifers and ferns are underrepresented in plant megafossil assemblages, bryophytes and lycopsids are represented only in sporomorph assemblages, and sphenophytes, pteridosperms, Caytoniales, Cycadales, Ginkgoales and Bennettitales are comparatively underrepresented in sporomorph assemblages. Combined multivariate analysis (correspondence analysis and nonmetric multidimensional scaling) of sporomorph occurrence/abundance data demonstrates that temporal variation in sporomorph assemblages is the result of depositional change through the plant bed. The reproductive strategies of parent plants are considered to be a principal factor in shaping many of the major abundance and diversity irregularities between dispersed sporomorph and plant megafossil data sets that seemingly reflects different parent vegetation. Preferential occurrence/preservation of sporomorphs and equivalent parent plants is a consequence of a complex array of biological, ecological, geographical, taphonomic, and depositional factors that act inconsistently between and within fossil assemblages, which results in notable discrepancies between data sets.

Proterozoic strata host evidence of global “Snowball Earth” glaciations, large perturbations to the carbon cycle, proposed changes in the redox state of oceans, the diversification of microscopic eukaryotes, and the rise of metazoans. Over the past half century, the number of fossils described from Proterozoic rocks has increased exponentially. These discoveries have occurred alongside an increased understanding of the Proterozoic Earth system and the geological context of fossil occurrences, including improved age constraints. However, the evaluation of relationships between Proterozoic environmental change and fossil diversity has been hampered by several factors, particularly lithological and taphonomic biases. Here we compile and analyze the current record of eukaryotic fossils in Proterozoic strata to assess the effect of biases and better constrain diversity through time. Our results show that mean within assemblage diversity increases through the Proterozoic Eon due to an increase in high diversity assemblages, and that this trend is robust to various external factors including lithology and paleogeographic location. In addition, assemblage composition changes dramatically through time. Most notably, robust recalcitrant taxa appear in the early Neoproterozoic Era, only to disappear by the beginning of the Ediacaran Period. Within assemblage diversity is significantly lower in the Cryogenian Period than in the preceding and following intervals, but the short duration of the nonglacial interlude and unusual depositional conditions may present additional biases. In general, large scale patterns of diversity are robust while smaller scale patterns are difficult to discern through the lens of lithological, taphonomic, and geographic variability.

The extinct shark Carcharocles megalodon is one of the largest marine apex predators ever to exist. Nonetheless, little is known about its body-size variations through time and space. Here, we studied the body-size trends of C. megalodon through its temporal and geographic range to better understand its ecology and evolution. Given that this species was the last of the megatooth lineage, a group of species that shows a purported size increase through time, we hypothesized that C. megalodon also displayed this trend, increasing in size over time and reaching its largest size prior to extinction. We found that C. megalodon body-size distribution was left-skewed (suggesting a long-term selective pressure favoring larger individuals), and presented significant geographic variation (possibly as a result of the heterogeneous ecological constraints of this cosmopolitan species) over geologic time. Finally, we found that stasis was the general mode of size evolution of C. megalodon (i.e., no net changes over time), contrasting with the trends of the megatooth lineage and our hypothesis. Given that C. megalodon is a relatively long-lived species with a widely distributed fossil record, we further used this study system to provide a deep-time perspective to the understanding of the body-size trends of marine apex predators. For instance, our results suggest that (1) a selective pressure in predatory sharks for consuming a broader range of prey may favor larger individuals and produce left-skewed distributions on a geologic time scale; (2) body-size variations in cosmopolitan apex marine predators may depend on their interactions with geographically discrete communities; and (3) the inherent characteristics of shark species can produce stable sizes over geologic time, regardless of the size trends of their lineages.

Although Phanerozoic increases in the global richness, local richness, and evenness of marine invertebrates are well documented, a common explanation for these patterns has been difficult to identify. Evidence is presented here from marine invertebrate communities that there is a Phanerozoic increase in the fundamental biodiversity number (θ), which describes diversity and relative abundance distributions in neutral ecological theory. If marine ecosystems behave according to the rules of Hubbell’s Neutral Theory of Biodiversity and Biogeography, the Phanerozoic increase in θ suggests three possible mechanisms for the parallel increases in global richness, local richness, and evenness: (1) an increase in the per-individual probability of speciation, (2) an increase in the area occupied by marine metacommunities, and (3) an increase in the density (per-area abundance) of marine organisms. Because speciation rates have declined over time and because there is no clear evidence for an increase in metacommunity area through the Phanerozoic, the most likely of these is an increase in the spatial density of marine invertebrates over the Phanerozoic, an interpretation supported by previous studies of fossil abundance. This, coupled with a Phanerozoic rise in body size, suggests that an increase in primary productivity through time is the primary cause of Phanerozoic increases in θ, global richness, local richness, local evenness, abundance, and body size.