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.     

Body length of bony fishes was not a selective factor during the biggest mass extinction of all time

The Permo‐Triassic mass extinction devastated life on land and in the sea, but it is not clear why some species survived and others went extinct. One explanation is that lineage loss during mass extinctions is a random process in which luck determines which species survive. Alternatively, a phylogenetic signal in extinction may indicate a selection process operating on phenotypic traits. Large body size has often emerged as an extinction risk factor in studies of modern extinction risk, but this is not so commonly the case for mass extinctions in deep time. Here, we explore the evolution of non‐teleostean Actinopterygii (bony fishes) from the Devonian to the present day, and we concentrate on the Permo‐Triassic mass extinction. We apply a variety of time‐scaling metrics to date the phylogeny, and show that diversity peaked in the latest Permian and declined severely during the Early Triassic. In line with previous evidence, we find the phylogenetic signal of extinction increases across the mass extinction boundary: extinction of species in the earliest Triassic is more clustered across phylogeny compared to the more randomly distributed extinction signal in the late Permian. However, body length plays no role in differential survival or extinction of taxa across the boundary. In the case of fishes, size did not determine which species survived and which went extinct, but phylogenetic signal indicates that the mass extinction was not a random field of bullets.     

Marine extinction risk shaped by trait–environment interactions over 500 million years

Perhaps the most pressing issue in predicting biotic responses to present and future global change is understanding how environmental factors shape the relationship between ecological traits and extinction risk. The fossil record provides millions of years of insight into how extinction selectivity (i.e., differential extinction risk) is shaped by interactions between ecological traits and environmental conditions. Numerous paleontological studies have examined trait‐based extinction selectivity; however, the extent to which these patterns are shaped by environmental conditions is poorly understood due to a lack of quantitative synthesis across studies. We conducted a meta‐analysis of published studies on fossil marine bivalves and gastropods that span 458 million years to uncover how global environmental and geochemical changes covary with trait‐based extinction selectivity. We focused on geographic range size and life habit (i.e., infaunal vs. epifaunal), two of the most important and commonly examined predictors of extinction selectivity. We used geochemical proxies related to global climate, as well as indicators of ocean acidification, to infer average global environmental conditions. Life‐habit selectivity is weakly dependent on environmental conditions, with infaunal species relatively buffered from extinction during warmer climate states. In contrast, the odds of taxa with broad geographic ranges surviving an extinction (>2500 km for genera, >500 km for species) are on average three times greater than narrow‐ranging taxa (estimate of odds ratio: 2.8, 95% confidence interval = 2.3–3.5), regardless of the prevailing global environmental conditions. The environmental independence of geographic range size extinction selectivity emphasizes the critical role of geographic range size in setting conservation priorities.     

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.     

Primate extinction risk and historical patterns of speciation and extinction in relation to body mass

Body mass is thought to influence diversification rates, but previous studies have produced ambiguous results. We investigated patterns of diversification across 100 trees obtained from a new Bayesian inference of primate phylogeny that sampled trees in proportion to their posterior probabilities. First, we used simulations to assess the validity of previous studies that used linear models to investigate the links between IUCN Red List status and body mass. These analyses support the use of linear models for ordinal ranked data on threat status, and phylogenetic generalized linear models revealed a significant positive correlation between current extinction risk and body mass across our tree block. We then investigated historical patterns of speciation and extinction rates using a recently developed maximum-likelihood method. Specifically, we predicted that body mass correlates positively with extinction rate because larger bodied organisms reproduce more slowly, and body mass correlates negatively with speciation rate because smaller bodied organisms are better able to partition niche space. We failed to find evidence that extinction rates covary with body mass across primate phylogeny. Similarly, the speciation rate was generally unrelated to body mass, except in some tests that indicated an increase in the speciation rate with increasing body mass. Importantly, we discovered that our data violated a key assumption of sample randomness with respect to body mass. After correcting for this bias, we found no association between diversification rates and mass.     

Nature and timing of biotic recovery in Antarctic benthic marine ecosystems following the Cretaceous–Palaeogene mass extinction

Taxonomic and ecological recovery from the Cretaceous–Palaeogene (K–Pg) mass extinction 66 million years ago shaped the composition and structure of modern ecosystems. The timing and nature of recovery has been linked to many factors including palaeolatitude, geographical range, the ecology of survivors, incumbency and palaeoenvironmental setting. Using a temporally constrained fossil dataset from one of the most expanded K–Pg successions in the world, integrated with palaeoenvironmental information, we provide the most detailed examination of the patterns and timing of recovery from the K–Pg mass extinction event in the high southern latitudes of Antarctica. The timing of biotic recovery was influenced by global stabilization of the wider Earth system following severe environmental perturbations, apparently regardless of latitude or local environment. Extinction intensity and ecological change were decoupled, with community scale ecological change less distinct compared to other locations, even if the taxonomic severity of the extinction was the same as at lower latitudes. This is consistent with a degree of geographical heterogeneity in the recovery from the K–Pg mass extinction. Recovery in Antarctica was influenced by local factors (such as water depth changes, local volcanism, and possibly incumbency and pre‐adaptation to seasonality of the local benthic molluscan population), and also showed global signals, for example the radiation of the Neogastropoda within the first million years of the Danian, and a shift in dominance between bivalves and gastropods.     

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.     

Geographical distribution and extinction risk: lessons from Triassic–Jurassic marine benthic organisms

Aim To evaluate the influence of geographical distribution on the extinction risk of benthic marine invertebrates using data from the fossil record, both during times of background extinction and across a mass‐extinction episode. Total geographical range is contrasted with proxies of global abundance to assess the relationships between the two essential components of geographical distribution and extinction risk. Location A global occurrence data base of fossil benthic macro‐organisms from the Triassic and Jurassic periods was used for this study. Methods Geographical distributions and biodiversity dynamics were assessed for each genus (all taxa) or species (bivalves) based on a sample‐standardized data set and palaeogeographical reconstructions. Geographical ranges were measured by the maximum great circle distance of a taxon within a stratigraphic interval. Global abundance was assessed by the number of localities at which a taxon was recorded. Widespread and rare taxa were separated using median and percentile values of the frequency distributions of occurrences. Results The frequency distribution of geographical ranges is very similar to that for modern taxa. Although no significant correlation could be established between local abundance and geographical range, proxies of global abundance are strongly correlated with geographical range. Taxon longevities are correlated with both mean geographical range and mean global abundance, but range size appears to be more critical than abundance in determining extinction risk. These results are valid when geographical distribution is treated as a trait of taxa and when assessed for individual geological stages. Main conclusions Geographical distribution is a key predictor of extinction risk of Triassic and Jurassic benthic marine invertebrates. An important exception is in the end‐Triassic mass extinction, which equally affected geographically restricted and widespread genera, as well as common and rare genera. This suggests that global diversity crises may curtail the role of geographical distribution in determining extinction risk.     

Scaling the extinction vortex: Body size as a predictor of population dynamics close to extinction events

1. Mutual reinforcement between abiotic and biotic factors can drive small populations into a catastrophic downward spiral to extinction—a process known as the “extinction vortex.” However, empirical studies investigating extinction dynamics in relation to species'' traits have been lacking.
2. We assembled a database of 35 vertebrate populations monitored to extirpation over a period of at least ten years, represented by 32 different species, including 25 birds, five mammals, and two reptiles. We supplemented these population time series with species‐specific mean adult body size to investigate whether this key intrinsic trait affects the dynamics of populations declining toward extinction.
3. We performed three analyses to quantify the effects of adult body size on three characteristics of population dynamics: time to extinction, population growth rate, and residual variability in population growth rate.
4. Our results provide support for the existence of extinction vortex dynamics in extirpated populations. We show that populations typically decline nonlinearly to extinction, while both the rate of population decline and variability in population growth rate increase as extinction is approached. Our results also suggest that smaller‐bodied species are particularly prone to the extinction vortex, with larger increases in rates of population decline and population growth rate variability when compared to larger‐bodied species.
5. Our results reaffirm and extend our understanding of extinction dynamics in real‐life extirpated populations. In particular, we suggest that smaller‐bodied species may be at greater risk of rapid collapse to extinction than larger‐bodied species, and thus, management of smaller‐bodied species should focus on maintaining higher population abundances as a priority.

     

Length–mass allometry in snakes

Body size and body shape are tightly related to an animal's physiology, ecology and life history, and, as such, play a major role in understanding ecological and evolutionary phenomena. Because organisms have different shapes, only a uniform proxy of size, such as mass, may be suitable for comparisons between taxa. Unfortunately, snake masses are rarely reported in the literature. On the basis of 423 species of snakes in 10 families, we developed clade‐specific equations for the estimation of snake masses from snout–vent lengths and total lengths. We found that snout–vent lengths predict masses better than total lengths. By examining the effects of phylogeny, as well as ecological and life history traits on the relationship between mass and length, we found that viviparous species are heavier than oviparous species, and diurnal species are heavier than nocturnal species. Furthermore, microhabitat preferences profoundly influence body shape: arboreal snakes are lighter than terrestrial snakes, whereas aquatic snakes are heavier than terrestrial snakes of a similar length. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●● , ●●–●●.     

Escalation and extinction selectivity: morphology versus isotopic reconstruction of bivalve metabolism

Studies that have tested and failed to support the hypothesis that escalated species (e.g., those with predation-resistant adaptations) are more susceptible to elimination during mass extinctions have concentrated on the distribution and degree of morphological defenses in molluscan species. This morphological approach to determining level of escalation in bivalves may be oversimplified because it does not account for metabolic rate, which is an important measure of escalation that is less readily accessible for fossils. Shell growth rates in living bivalves are positively correlated with metabolic rate and thus are potential indicators of level of escalation. To evaluate this approach, we used oxygen isotopes to reconstruct shell growth rates for two bivalve species (Macrocallista marylandica and Glossus markoei) from Miocene-aged sediments of Maryland. Although both species are classified as non-escalated based on morphology, the isotopic data indicate that M. marylandica was a faster-growing species with a higher metabolic rate and G. markoei was a slower-growing species with a lower metabolic rate. Based on these results, we predict that some morphologically non-escalated species in previous tests of extinction selectivity should be reclassified as escalated because of their fast shell growth rates (i.e., high metabolic rates). Studies that evaluate the level of escalation of a fauna should take into account the energetic physiology of taxa to avoid misleading results.     

Questions of mass extinction

Earth's biodiversity is being overtaken by a mass extinction which, if allowed to proceed unchecked, could well eliminate between one quarter and one half of all species. Our conservation responses must be science-based if we are to address the problem in its full scope and with most productive use of conservation resources. Yet our scientific understanding of the impending mass extinction is inadequate in many salient respects. We have only a rudimentary grasp of the number of species at risk, of biodiversity depletion processes, of island biogeography in practice, and of evolutionary consequences, to cite but a few leading questions. The same applies to the issue of the most efficient strategies to confront the conservation challenge. Worse, there is scant evidence (due in part to gross lack of funding) of a comprehensive and coordinated campaign to mount a research effort of scope to match the problem. The paper broaches ten key questions that warrant urgent attention.     

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.     

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.     

Body mass and extinction risk in Australian marsupials: The ‘Critical Weight Range’ revisited

Australian mammals have suffered an exceptionally high rate of decline and extinction over the last two hundred years. Body mass is linked to extinction risk in Australian mammals, but the nature of this association is controversial. A widely held view is that species of intermediate body mass (between 35 and 5500 g, the ‘critical weight range’, CWR) have declined most severely. However, the existence of the CWR has been disputed. In this paper we clarify the relationship of decline status and body mass in Australian marsupials. We show that the form of this relationship differs for ground‐living versus arboreal species, and for species from low versus high rainfall areas. Among ground‐living species and those from low‐rainfall areas, declines were strongly size‐selective and concentrated on species within the CWR. For the remaining species, decline was only weakly related to body mass with no evidence of heightened risk for species of intermediate size. We conclude that for terrestrial species in low rainfall areas, species within the CWR are most at risk of decline and extinction.     

Severe extinction and rapid recovery of mammals across the Cretaceous–Palaeogene boundary,and the effects of rarity on patterns of extinction and recovery

The end‐Cretaceous mass extinction ranks among the most severe extinctions of all time; however, patterns of extinction and recovery remain incompletely understood. In particular, it is unclear how severe the extinction was, how rapid the recovery was and how sampling biases might affect our understanding of these processes. To better understand terrestrial extinction and recovery and how sampling influences these patterns, we collected data on the occurrence and abundance of fossil mammals to examine mammalian diversity across the K‐Pg boundary in North America. Our data show that the extinction was more severe and the recovery more rapid than previously thought. Extinction rates are markedly higher than previously estimated: of 59 species, four survived (93% species extinction, 86% of genera). Survival is correlated with geographic range size and abundance, with widespread, common species tending to survive. This creates a sampling artefact in which rare species are both more vulnerable to extinction and less likely to be recovered, such that the fossil record is inherently biased towards the survivors. The recovery was remarkably rapid. Within 300 000 years, local diversity recovered and regional diversity rose to twice Cretaceous levels, driven by increased endemicity; morphological disparity increased above levels observed in the Cretaceous. The speed of the recovery tends to be obscured by sampling effects; faunas show increased endemicity, such that a rapid, regional increase in diversity and disparity is not seen in geographically restricted studies. Sampling biases that operate against rare taxa appear to obscure the severity of extinction and the pace of recovery across the K‐Pg boundary, and similar biases may operate during other extinction events.     

Biogeography of worm lizards (Amphisbaenia) driven by end-Cretaceous mass extinction

Worm lizards (Amphisbaenia) are burrowing squamates that live as subterranean predators. Their underground existence should limit dispersal, yet they are widespread throughout the Americas, Europe and Africa. This pattern was traditionally explained by continental drift, but molecular clocks suggest a Cenozoic diversification, long after the break-up of Pangaea, implying dispersal. Here, we describe primitive amphisbaenians from the North American Palaeocene, including the oldest known amphisbaenian, and provide new and older molecular divergence estimates for the clade, showing that worm lizards originated in North America, then radiated and dispersed in the Palaeogene following the Cretaceous-Palaeogene (K-Pg) extinction. This scenario implies at least three trans-oceanic dispersals: from North America to Europe, from North America to Africa and from Africa to South America. Amphisbaenians provide a striking case study in biogeography, suggesting that the role of continental drift in biogeography may be overstated. Instead, these patterns support Darwin and Wallace''s hypothesis that the geographical ranges of modern clades result from dispersal, including oceanic rafting. Mass extinctions may facilitate dispersal events by eliminating competitors and predators that would otherwise hinder establishment of dispersing populations, removing biotic barriers to dispersal.     

Evolutionary and biogeographical shifts in response to the Late Ordovician mass extinction

The Late Ordovician mass extinction was an interval of high extinction with inferred low ecological selectivity, resulting in little change in community structure after the event. In contrast, the mass extinction may have fundamentally changed evolutionary dynamics in the surviving groups. We investigated the phylogenetic relationships among strophomenoid brachiopods, a diverse brachiopod superfamily that was a primary component of Ordovician ecosystems. Four Ordovician families/subfamilies sampled in the analysis (Rafinesquinidae, Strophomeninae, Glyptomenidae and Furcitellinae) were reconstructed as monophyletic groups, and the base of the strophomenoid clade that dominated the Silurian recovery was reconstructed as diversifying alongside these families during the Middle Ordovician. We time‐calibrated the phylogeny and used geographical occurrences to investigate biogeographical changes in the strophomenoids through time with the R package BiogeoBEARS . Our results indicate that extinction was higher in taxa whose ranges were constrained to tropical or subtropical regions. Furthermore, our results suggest important shifts in the diversification patterns of these brachiopods after the mass extinction. While most of the strophomenoid families survived the Late Ordovician event, ecologically abundant taxonomic groups during the Ordovician were either driven to extinction, reduced in diversity, or slowly died off during the Silurian. The new abundant strophomenoid taxa derived from one clade (consisting of Silurian–Devonian groups such as Douvillinidae, Strophodontidae and Amphistrophiidae) that diversified during the post‐extinction radiation. Our results suggest the selective diversification during the Silurian radiation, rather than selective extinction in the Late Ordovician, had a greater impact on the evolutionary history of strophomenoid brachiopods.     

Drivers of extinction: the case of Azorean beetles

Oceanic islands host a disproportionately high fraction of endangered or recently extinct endemic species. We report on species extinctions among endemic Azorean beetles following 97% habitat loss since AD 1440. We infer extinctions from historical and contemporary records and examine the influence of three predictors: geographical range, habitat specialization and body size. Of 55 endemic beetle species investigated (out of 63), seven can be considered extinct. Single-island endemics (SIEs) were more prone to extinction than multi-island endemics. Within SIEs restricted to native habitat, larger species were more extinction-prone. We thus show a hierarchical path to extinction in Azorean beetles: species with small geographical range face extinction first, with the larger bodied ones being the most threatened. Our study provides a clear warning of the impact of habitat loss on island endemic biotas.