Phanerozoic survivors: Actinopterygian evolution through the Permo‐Triassic and Triassic‐Jurassic mass extinction events

Actinopterygians (ray‐finned fishes) successfully passed through four of the big five mass extinction events of the Phanerozoic, but the effects of these crises on the group are poorly understood. Many researchers have assumed that the Permo‐Triassic mass extinction (PTME) and end‐Triassic extinction (ETE) had little impact on actinopterygians, despite devastating many other groups. Here, two morphometric techniques, geometric (body shape) and functional (jaw morphology), are used to assess the effects of these two extinction events on the group. The PTME elicits no significant shifts in functional disparity while body shape disparity increases. An expansion of body shape and functional disparity coincides with the neopterygian radiation and evolution of novel feeding adaptations in the Middle‐Late Triassic. Through the ETE, small decreases are seen in shape and functional disparity, but are unlikely to represent major changes brought about by the extinction event. In the Early Jurassic, further expansions into novel areas of ecospace indicative of durophagy occur, potentially linked to losses in the ETE. As no evidence is found for major perturbations in actinopterygian evolution through either extinction event, the group appears to have been immune to two major environmental crises that were disastrous to most other organisms.     

Recovery from the most profound mass extinction of all time

The end-Permian mass extinction, 251 million years (Myr) ago, was the most devastating ecological event of all time, and it was exacerbated by two earlier events at the beginning and end of the Guadalupian, 270 and 260 Myr ago. Ecosystems were destroyed worldwide, communities were restructured and organisms were left struggling to recover. Disaster taxa, such as Lystrosaurus, insinuated themselves into almost every corner of the sparsely populated landscape in the earliest Triassic, and a quick taxonomic recovery apparently occurred on a global scale. However, close study of ecosystem evolution shows that true ecological recovery was slower. After the end-Guadalupian event, faunas began rebuilding complex trophic structures and refilling guilds, but were hit again by the end-Permian event. Taxonomic diversity at the alpha (community) level did not recover to pre-extinction levels; it reached only a low plateau after each pulse and continued low into the Late Triassic. Our data showed that though there was an initial rise in cosmopolitanism after the extinction pulses, large drops subsequently occurred and, counter-intuitively, a surprisingly low level of cosmopolitanism was sustained through the Early and Middle Triassic.     

The early evolution of ray‐finned fishes

Ray‐finned fishes (Actinopterygii) constitute approximately half of all living vertebrate species. A stable hypothesis of relationships among major modern lineages has emerged over the past decade, supported by both anatomy and molecules. Diversity is unevenly partitioned across the actinopterygian tree, with most species concentrated within a handful of geologically young (i.e. Cretaceous) teleost clades. Extant non‐teleost groups are portrayed as ‘living fossils’, but this moniker should not be taken as evidence of especially primitive structure: each of these lineages is characterized by profound specializations. Attribution of fossils to the crowns and apical stems of Cladistia, Chondrostei and Neopterygii is uncontroversial, but placements of Palaeozoic taxa along deeper branches of actinopterygian phylogeny are less secure. Despite these limitations, some major outlines of actinopterygian diversification seem reasonably clear from the fossil record: low richness and disparity in the Devonian; elevated morphological variety, linked to increases in taxonomic dominance, in the early Carboniferous; and further gains in taxonomic dominance in the Early Triassic associated with earliest appearance of trophically diverse crown neopterygians.     

Archosauromorph extinction selectivity during the Triassic–Jurassic mass extinction

Many traits have been linked to extinction risk among modern vertebrates, including mode of life and body size. However, previous work has indicated there is little evidence that body size, or any other trait, was selective during past mass extinctions. Here, we investigate the impact of the Triassic–Jurassic mass extinction on early Archosauromorpha (basal dinosaurs, crocodylomorphs and their relatives) by focusing on body size and other life history traits. We built several new archosauromorph maximum‐likelihood supertrees, incorporating uncertainty in phylogenetic relationships. These supertrees were then employed as a framework to test whether extinction had a phylogenetic signal during the Triassic–Jurassic mass extinction, and whether species with certain traits were more or less likely to go extinct. We find evidence for phylogenetic signal in extinction, in that taxa were more likely to become extinct if a close relative also did. However, there is no correlation between extinction and body size, or any other tested trait. These conclusions add to previous findings that body size, and other traits, were not subject to selection during mass extinctions in closely‐related clades, although the phylogenetic signal in extinction indicates that selection may have acted on traits not investigated here.     

The complex effects of mass extinctions on morphological disparity

Studies of biodiversity through deep time have been a staple for biologists and paleontologists for over 60 years. Investigations of species richness (diversity) revealed that at least five mass extinctions punctuated the last half billion years, each seeing the rapid demise of a large proportion of contemporary taxa. In contrast to diversity, the response of morphological diversity (disparity) to mass extinctions is unclear. Generally, diversity and disparity are decoupled, such that diversity may decline as morphological disparity increases, and vice versa. Here, we develop simulations to model disparity changes across mass extinctions using continuous traits and birth-death trees. We find no simple null for disparity change following a mass extinction but do observe general patterns. The range of trait values decreases following either random or trait-selective mass extinctions, whereas variance and the density of morphospace occupation only decline following trait-selective events. General trends may differentiate random and trait-selective mass extinctions, but methods struggle to identify trait selectivity. Long-term effects of mass extinction trait selectivity change support for phylogenetic comparative methods away from the simulated Brownian motion toward Ornstein-Uhlenbeck and Early Burst models. We find that morphological change over mass extinction is best studied by quantifying multiple aspects of morphospace occupation.     

Reduction in animal abundance and oxygen availability during and after the end-Triassic mass extinction

The end-Triassic biodiversity crisis was one of the most severe mass extinctions in the history of animal life. However, the extent to which the loss of taxonomic diversity was coupled with a reduction in organismal abundance remains to be quantified. Further, the temporal relationship between organismal abundance and local marine redox conditions is lacking in carbonate sections. To address these questions, we measured skeletal grain abundance in shallow-marine limestones by point counting 293 thin sections from four stratigraphic sections across the Triassic/Jurassic boundary in the Lombardy Basin and Apennine Platform of western Tethys. Skeletal abundance decreased abruptly across the Triassic/Jurassic boundary in all stratigraphic sections. The abundance of skeletal organisms remained low throughout the lower-middle Hettangian strata and began to rebound during the late Hettangian and early Sinemurian. A two-way ANOVA indicates that sample age (p < .01, η² = 0.30) explains more of the variation in skeletal abundance than the depositional environment or paleobathymetry (p < .01, η² = 0.15). Measured I/Ca ratios, a proxy for local shallow-marine redox conditions, show this same pattern with the lowest I/Ca ratios occurring in the early Hettangian. The close correspondence between oceanic water column oxygen levels and skeletal abundance indicates a connection between redox conditions and benthic organismal abundance across the Triassic/Jurassic boundary. These findings indicate that the end-Triassic mass extinction reduced not only the biodiversity but also the carrying capacity for skeletal organisms in early Hettangian ecosystems, adding to evidence that mass extinction of species generally leads to mass rarity among survivors.     

Body size changes in bivalves of the family Limidae in the aftermath of the end‐Triassic mass extinction: the Brobdingnag effect

Reduction in body size of organisms following mass extinctions is well‐known and often ascribed to the Lilliput effect. This phenomenon is expressed as a temporary body size reduction within surviving species. Despite its wide usage the term is often loosely applied to any small post‐extinction taxa. Here we assess the size of bivalves of the family Limidae (Rafineque) prior to, and in the aftermath of, the end‐Triassic mass extinction event. Of the species studied only one occurs prior to the extinction event, though is too scarce to test for the Lilliput effect. Instead, newly evolved species originate at small body sizes and undergo a within‐species size increase, most dramatically demonstrated by Plagiostoma giganteum (Sowerby) which, over two million years, increases in size by 179%. This trend is seen in both field and museum collections. We term this within‐species size increase of newly originated species in the aftermath of mass extinction, the Brobdingnag effect, after the giants that were contemporary with the Lilliputians in Swift's Gulliver's Travels. The size increase results from greater longevity and faster growth rates. The cause of the effect is unclear, although it probably relates to improved environmental conditions. Oxygen‐poor conditions in the Early Jurassic are associated with populations of smaller body size caused by elevated juvenile mortality but these are local/regional effects that do not alter the long‐term, size increase. Although temperature‐size relationships exist for many organisms (Temperature‐Size Rule and Bergmann's Rule), the importance of this is unclear here because of a poorly known Early Jurassic temperature record.     

Oceanic stable isotope composition and a scenario for the Permo‐Triassic crisis

Stable carbon isotope data from brachiopod shells from the Upper Permian Kapp Starostin Formation (West Spitsbergen) indicate that the oceanic carbon isotopic ratio, which had already been very high in the late Permian, rapidly increased by almost 4 per mil and then dramatically declined by more than 10 per mil in the very latest Permian. This pattern is essentially repeated by the oxygen isotope curve. These data show that a geologically rapid switch between two fundamentally different states of the Earth's exosystem occurred near the Permo‐Triassic transition. The late Permian state of the global system was profoundly different from the modern one in that vast amounts of organic carbon were stored, presumably in the form of easy‐to‐mobilize sapropel‐like deposits, below the oceanic redoxcline. Under such conditions—which we propose to call overfed ocean—nutrients were intensely recycled to seawater, thus allowing the ocean to sustain a huge standing crop of the biosphere. Deposition of large amounts of organic carbon in the ocean liberated corresponding amounts of oxygen, thus leading to high oxygen levels in the atmosphere. In the latest Permian, the organic matter decaying in the ocean was subject to rapid oxidation due to appearance of the modern type of ocean which is characterized by vigorous bottom circulation and net heterotrophy. The appearance of these conditions—which we propose to call hungry ocean—led to removal of nutrients from seawater and to a substantial drop in atmospheric oxygen contents. The resulting nutrient deficiency in the ocean, oxygen depletion in the atmosphere, and other effects of this paleoceanographic change must have caused major extinctions.     

Hunting to extinction: biology and regional economy influence extinction risk and the impact of hunting in artiodactyls

Half of all artiodactyls (even-toed hoofed mammals) are threatened with extinction, around double the mammalian average. Here, using a complete species-level phylogeny, we construct a multivariate model to assess for the first time which intrinsic (biological) and extrinsic (anthropogenic and environmental) factors influence variation in extinction risk in artiodactyls. Globally artiodactyls at greatest risk live in economically less developed areas, have older weaning ages and smaller geographical ranges. Our findings suggest that identifying predictors of threat is complicated by interactions between both biological and anthropogenic factors, resulting in differential responses to threatening processes. Artiodactyl species that experience unregulated hunting live in significantly less economically developed areas than those that are not hunted; however, hunted species are more susceptible to extinction if they have slower reproductive rates (older weaning ages). In contrast, risk in non-hunted artiodactyls is unrelated to reproductive rate and more closely associated with the economic development of the region in which they live.     

Bias in phylogenetic measurements of extinction and a case study of end‐Permian tetrapods

Extinction risk in the modern world and extinction in the geological past are often linked to aspects of life history or other facets of biology that are phylogenetically conserved within clades. These links can result in phylogenetic clustering of extinction, a measurement comparable across different clades and time periods that can be made in the absence of detailed trait data. This phylogenetic approach is particularly suitable for vertebrate taxa, which often have fragmentary fossil records, but robust, cladistically‐inferred trees. Here we use simulations to investigate the adequacy of measures of phylogenetic clustering of extinction when applied to phylogenies of fossil taxa while assuming a Brownian motion model of trait evolution. We characterize expected biases under a variety of evolutionary and analytical scenarios. Recovery of accurate estimates of extinction clustering depends heavily on the sampling rate, and results can be highly variable across topologies. Clustering is often underestimated at low sampling rates, whereas at high sampling rates it is always overestimated. Sampling rate dictates which cladogram timescaling method will produce the most accurate results, as well as how much of a bias ancestor–descendant pairs introduce. We illustrate this approach by applying two phylogenetic metrics of extinction clustering (Fritz and Purvis's D and Moran's I) to three tetrapod clades across an interval including the Permo‐Triassic mass extinction event. These groups consistently show phylogenetic clustering of extinction, unrelated to change in other quantitative metrics such as taxonomic diversity or extinction intensity.     

Impact of the Late Triassic mass extinction on functional diversity and composition of marine ecosystems

Mass extinctions have profoundly influenced the history of life, not only through the death of species but also through changes in ecosystem function and structure. Importantly, these events allow us the opportunity to study ecological dynamics under levels of environmental stress for which there are no recent analogues. Here, we examine the impact and selectivity of the Late Triassic mass extinction event on the functional diversity and functional composition of the global marine ecosystem, and test whether post‐extinction communities in the Early Jurassic represent a regime shift away from pre‐extinction communities in the Late Triassic. Our analyses show that, despite severe taxonomic losses, there is no unequivocal loss of global functional diversity associated with the extinction. Even though no functional groups were lost, the extinction event was, however, highly selective against some modes of life, in particular sessile suspension feeders. Although taxa with heavily calcified skeletons suffered higher extinction than other taxa, lightly calcified taxa also appear to have been selected against. The extinction appears to have invigorated the already ongoing faunal turnover associated with the Mesozoic Marine Revolution. The ecological effects of the Late Triassic mass extinction were preferentially felt in the tropical latitudes, especially amongst reefs, and it took until the Middle Jurassic for reef ecosystems to fully recover to pre‐extinction levels.     

Discrete and morphometric traits reveal contrasting patterns and processes in the macroevolutionary history of a clade of scorpions

Many palaeontological studies have investigated the evolution of entire body plans, generally relying on discrete character‐taxon matrices. In contrast, macroevolutionary studies performed by neontologists have mostly focused on morphometric traits. Although these data types are very different, some studies have suggested that they capture common patterns. Nonetheless, the tests employed to support this claim have not explicitly incorporated a phylogenetic framework and may therefore be susceptible to confounding effects due to the presence of common phylogenetic structure. We address this question using the scorpion genus Brachistosternus Pocock 1893 as case study. We make use of a time‐calibrated multilocus molecular phylogeny, and compile discrete and traditional morphometric data sets, both capturing the overall morphology of the organisms. We find that morphospaces derived from these matrices are significantly different, and that the degree of discordance cannot be replicated by simulations of random character evolution. Moreover, we find strong support for contrasting modes of evolution, with discrete characters being congruent with an ‘early burst’ scenario whereas morphometric traits suggest species‐specific adaptations to have driven morphological evolution. The inferred macroevolutionary dynamics are therefore contingent on the choice of character type. Finally, we confirm that metrics of correlation fail to detect these profound differences given common phylogenetic structure in both data sets, and that methods incorporating a phylogenetic framework and accounting for expected covariance should be favoured.     

The extinction of the dinosaurs

Non‐avian dinosaurs went extinct 66 million years ago, geologically coincident with the impact of a large bolide (comet or asteroid) during an interval of massive volcanic eruptions and changes in temperature and sea level. There has long been fervent debate about how these events affected dinosaurs. We review a wealth of new data accumulated over the past two decades, provide updated and novel analyses of long‐term dinosaur diversity trends during the latest Cretaceous, and discuss an emerging consensus on the extinction's tempo and causes. Little support exists for a global, long‐term decline across non‐avian dinosaur diversity prior to their extinction at the end of the Cretaceous. However, restructuring of latest Cretaceous dinosaur faunas in North America led to reduced diversity of large‐bodied herbivores, perhaps making communities more susceptible to cascading extinctions. The abruptness of the dinosaur extinction suggests a key role for the bolide impact, although the coarseness of the fossil record makes testing the effects of Deccan volcanism difficult.     

The influence of denitrifying seawater on graptolite extinction and diversification during the Hirnantian (latest Ordovician) mass extinction event

A continuous trench exposure within the uppermost type Vinini Formation at Vinini Creek, Roberts Mountains, Nevada, provides an unparalleled opportunity to examine the fate of graptolites, prominent Paleozoic zooplankton, during most of the Hirnantian mass extinction event. On the basis of a detailed biostratigraphic and sedimentological dataset, the relatively complete extinction record is examined in the context of ecological constraints, and it is found to reflect an ecological collapse driven by glacio-eustatic sea-level fall and associated changes in oceanic circulation. Diverse graptolite populations of the Dicranograptidae-Diplograptidae-Orthograptidae (DDO) fauna, which flourished in denitrifying waters within the oceanic oxygen-minimum zone (OMZ) during sea-level highstand, largely vanished with the loss of these conditions during glacio-eustatic sea-level fall. However, populations of one clade, the normalograptids, which inhabited the oxygenated waters of the photic zone, not only survived but diversified. These survivors gave rise to rapid recolonization and diversification with re-establishment of the oxygen-minimum and denitrifying conditions during post-Hirnantian sea-level rise. This ecological model also applies globally to other well-documented coeval stratigraphic intervals, representing both oceanic and platform sea settings.     

A multiple peak adaptive landscape based on feeding strategies and roosting ecology shaped the evolution of cranial covariance structure and morphological differentiation in phyllostomid bats

We explored the evolution of morphological integration in the most noteworthy example of adaptive radiation in mammals, the New World leaf‐nosed bats, using a massive dataset and by combining phylogenetic comparative methods and quantitative genetic approaches. We demonstrated that the phenotypic covariance structure remained conserved on a broader phylogenetic scale but also showed a substantial divergence between interclade comparisons. Most of the phylogenetic structure in the integration space can be explained by splits at the beginning of the diversification of major clades. Our results provide evidence for a multiple peak adaptive landscape in the evolution of cranial covariance structure and morphological differentiation, based upon diet and roosting ecology. In this scenario, the successful radiation of phyllostomid bats was triggered by the diversification of dietary and roosting strategies, and the invasion of these new adaptive zones lead to changes in phenotypic covariance structure and average morphology. Our results suggest that intense natural selection preceded the invasion of these new adaptive zones and played a fundamental role in shaping cranial covariance structure and morphological differentiation in this hyperdiverse clade of mammals. Finally, our study demonstrates the power of combining comparative methods and quantitative genetic approaches when investigating the evolution of complex morphologies.     

Historical and contingent factors affect re-evolution of a complex feature lost during mass extinction in communities of digital organisms

Re-evolution of complex biological features following the extinction of taxa bearing them remains one of evolution's most interesting phenomena, but is not amenable to study in fossil taxa. We used communities of digital organisms (computer programs that self-replicate, mutate and evolve), subjected to periods of low resource availability, to study the evolution, loss and re-evolution of a complex computational trait, the function EQU (bit-wise logical equals). We focused our analysis on cases where the pre-extinction EQU clade had surviving descendents at the end of the extinction episode. To see if these clades retained the capacity to re-evolve EQU, we seeded one set of multiple subreplicate 'replay' populations using the most abundant survivor of the pre-extinction EQU clade, and another set with the actual end-extinction ancestor of the organism in which EQU re-evolved following the extinction episode. Our results demonstrate that stochastic, historical, genomic and ecological factors can lead to constraints on further adaptation, and facilitate or hinder re-evolution of a complex feature.     

Body mass and foraging ecology predict evolutionary patterns of skeletal pneumaticity in the diverse "waterbird" clade

Extensive skeletal pneumaticity (air-filled bone) is a distinguishing feature of birds. The proportion of the skeleton that is pneumatized varies considerably among the >10,000 living species, with notable patterns including increases in larger bodied forms, and reductions in birds employing underwater pursuit diving as a foraging strategy. I assess the relationship between skeletal pneumaticity and body mass and foraging ecology, using a dataset of the diverse "waterbird" clade that encompasses a broad range of trait variation. Inferred changes in pneumaticity and body mass are congruent across different estimates of phylogeny, whereas pursuit diving has evolved independently between two and five times. Phylogenetic regressions detected positive relationships between body mass and pneumaticity, and negative relationships between pursuit diving and pneumaticity, whether independent variables are considered in isolation or jointly. Results are generally consistent across different estimates of topology and branch lengths. "Predictive" analyses reveal that several pursuit divers (loons, penguins, cormorants, darters) are significantly apneumatic compared to their relatives, and provide an example of how phylogenetic information can increase the statistical power to detect taxa that depart from established trait correlations. These findings provide the strongest quantitative comparative support yet for classical hypotheses regarding the evolution of avian skeletal pneumaticity.     

Exploring the power of Bayesian birth‐death skyline models to detect mass extinction events from phylogenies with only extant taxa

Mass extinction events (MEEs), defined as significant losses of species diversity in significantly short time periods, have attracted the attention of biologists because of their link to major environmental change. MEEs have traditionally been studied through the fossil record, but the development of birth‐death models has made it possible to detect their signature based on extant‐taxa phylogenies. Most birth‐death models consider MEEs as instantaneous events where a high proportion of species are simultaneously removed from the tree (“single pulse” approach), in contrast to the paleontological record, where MEEs have a time duration. Here, we explore the power of a Bayesian Birth‐Death Skyline (BDSKY) model to detect the signature of MEEs through changes in extinction rates under a “time‐slice” approach. In this approach, MEEs are time intervals where the extinction rate is greater than the speciation rate. Results showed BDSKY can detect and locate MEEs but that precision and accuracy depend on the phylogeny's size and MEE intensity. Comparisons of BDSKY with the single‐pulse Bayesian model, CoMET, showed a similar frequency of Type II error and neither model exhibited Type I error. However, while CoMET performed better in detecting and locating MEEs for smaller phylogenies, BDSKY showed higher accuracy in estimating extinction and speciation rates.     

Divergence of ovipositor length and egg shape in a brood parasitic bitterling fish through the use of different mussel hosts

Bitterling fishes deposit their eggs on the gills of living mussels using a long ovipositor. We examined whether ovipositor length (OL) and egg shape correlated with differences in host mussel species in the family Unionidae among populations of the tabira bitterling (Acheilognathus tabira) in Japan. Bitterling populations that use mussels in the sub-family Anodontinae possessed longer ovipositors and more elongated eggs than those using mussels of Unioninae, as expected from the difference in host size between the sub-families (anodontine mussels are larger than unionine mussels). Based on a robust phylogeny of A. tabira populations, we demonstrated that the evolution of both OL and egg shape were correlated with host differences, but not with each other, suggesting that these traits have been selected for independently. Our study demonstrates how adaptive traits for brood parasitism may diverge with host shift due to different host availability and/or interspecific competition for hosts.