Bivalves and evolutionary resilience: old skills and new strategies to recover from the P/T and T/J extinction events

Diversity dynamics among bivalves during the Triassic and Early Jurassic provides the opportunity to analyse the recovery patterns after two mass extinctions: Permian/Triassic and Triassic/Jurassic (T/J). The results presented here are based on a newly compiled worldwide genus-level database and are contrasted to the main morphological characters of the different taxonomical (orders and their constituent families and genera) and ecological groups. Many of such morphological characters are innovations appearing during the time span considered. Diversity and evolutionary rates were assessed and compared between these groups. During the Early Triassic there was a slow recovery, dominated by epifaunal taxa, the order Pectinida being the most diverse. The major post-Permian radiation took place during the Anisian, with several morphological and ecological innovations appearing and/or diversifying. The Late Triassic was a time of great diversification and ecological specialisation. Although the T/J was a true mass extinction for bivalves, it was not indiscriminate as its impact was stronger on specialised orders and not all ecological categories were equally affected. Recovery during earliest Jurassic was fast, confirming the high-evolutionary resilience of bivalve molluscs, except for groups with thick shells and tropical distribution, probably because of a biocalcification crisis.     

Phylogenetic Clustering of Origination and Extinction across the Late Ordovician Mass Extinction

Mass extinctions can have dramatic effects on the trajectory of life, but in some cases the effects can be relatively small even when extinction rates are high. For example, the Late Ordovician mass extinction is the second most severe in terms of the proportion of genera eliminated, yet is noted for the lack of ecological consequences and shifts in clade dominance. By comparison, the end-Cretaceous mass extinction was less severe but eliminated several major clades while some rare surviving clades diversified in the Paleogene. This disconnect may be better understood by incorporating the phylogenetic relatedness of taxa into studies of mass extinctions, as the factors driving extinction and recovery are thought to be phylogenetically conserved and should therefore promote both origination and extinction of closely related taxa. Here, we test whether there was phylogenetic selectivity in extinction and origination using brachiopod genera from the Middle Ordovician through the Devonian. Using an index of taxonomic clustering (R_CL) as a proxy for phylogenetic clustering, we find that A) both extinctions and originations shift from taxonomically random or weakly clustered within families in the Ordovician to strongly clustered in the Silurian and Devonian, beginning with the recovery following the Late Ordovician mass extinction, and B) the Late Ordovician mass extinction was itself only weakly clustered. Both results stand in stark contrast to Cretaceous-Cenozoic bivalves, which showed significant levels of taxonomic clustering of extinctions in the Cretaceous, including strong clustering in the mass extinction, but taxonomically random extinctions in the Cenozoic. The contrasting patterns between the Late Ordovician and end-Cretaceous events suggest a complex relationship between the phylogenetic selectivity of mass extinctions and the long-term phylogenetic signal in origination and extinction patterns.     

Phanerozoic changes in the high-rank suprageneric diversity structure of brachiopods: Linear and non-linear effects

Phanerozoic evolution of brachiopods produced many linear (established by a comparison of successive geologic time units) and non-linear (established by a comparison of non-successive geologic time units) effects, which can be examined quantitatively by using the similarity coefficients (Czekanowski's Quantified Coefficient and Gower Index) and correlation tools. The high-rank suprageneric diversity structure accounts for a number of superfamilies in each of 26 orders for every epoch of geological time. The intensity of turnovers in this structure was generally low during the entire Phanerozoic. It was slightly stronger during the Early Paleozoic, but close to zero during the Cenozoic, when the high-rank suprageneric diversity structure of brachiopods stabilized finally. Significant turnovers took place at the Middle Cambrian–Early Ordovician, the Late Ordovician–Early Silurian, the Late Silurian–Early Devonian, the Middle Devonian–Mississippian, and the Permian–Triassic transitions. Influences of mass extinctions, both major like those End Ordovician or Permian/Triassic and minor like Early Jurassic or Jurassic/Cretaceous, on the high-rank suprageneric diversity structure of brachiopods is registered. The strongest was the consequences of the Permian/Triassic catastrophe, which perhaps even reset the brachiopod evolution. No evident direct relationships are established between intensity of turnovers and eustatic fluctuations. However, the changes in the diversity structure recorded with the Gower Index provide evidence that eustatic lowstands were more favorable for intensification in these changes.     

Trophic network models explain instability of Early Triassic terrestrial communities

Studies of the end-Permian mass extinction have emphasized potential abiotic causes and their direct biotic effects. Less attention has been devoted to secondary extinctions resulting from ecological crises and the effect of community structure on such extinctions. Here we use a trophic network model that combines topological and dynamic approaches to simulate disruptions of primary productivity in palaeocommunities. We apply the model to Permian and Triassic communities of the Karoo Basin, South Africa, and show that while Permian communities bear no evidence of being especially susceptible to extinction, Early Triassic communities appear to have been inherently less stable. Much of the instability results from the faster post-extinction diversification of amphibian guilds relative to amniotes. The resulting communities differed fundamentally in structure from their Permian predecessors. Additionally, our results imply that changing community structures over time may explain long-term trends like declining rates of Phanerozoic background extinction.     

Permian–Triassic Osteichthyes (bony fishes): diversity dynamics and body size evolution

The Permian and Triassic were key time intervals in the history of life on Earth. Both periods are marked by a series of biotic crises including the most catastrophic of such events, the end‐Permian mass extinction, which eventually led to a major turnover from typical Palaeozoic faunas and floras to those that are emblematic for the Mesozoic and Cenozoic. Here we review patterns in Permian–Triassic bony fishes, a group whose evolutionary dynamics are understudied. Based on data from primary literature, we analyse changes in their taxonomic diversity and body size (as a proxy for trophic position) and explore their response to Permian–Triassic events. Diversity and body size are investigated separately for different groups of Osteichthyes (Dipnoi, Actinistia, ‘Palaeopterygii’, ‘Subholostei’, Holostei, Teleosteomorpha), within the marine and freshwater realms and on a global scale (total diversity) as well as across palaeolatitudinal belts. Diversity is also measured for different palaeogeographical provinces. Our results suggest a general trend from low osteichthyan diversity in the Permian to higher levels in the Triassic. Diversity dynamics in the Permian are marked by a decline in freshwater taxa during the Cisuralian. An extinction event during the end‐Guadalupian crisis is not evident from our data, but ‘palaeopterygians’ experienced a significant body size increase across the Guadalupian–Lopingian boundary and these fishes upheld their position as large, top predators from the Late Permian to the Late Triassic. Elevated turnover rates are documented at the Permian–Triassic boundary, and two distinct diversification events are noted in the wake of this biotic crisis, a first one during the Early Triassic (dipnoans, actinistians, ‘palaeopterygians’, ‘subholosteans’) and a second one during the Middle Triassic (‘subholosteans’, neopterygians). The origination of new, small taxa predominantly among these groups during the Middle Triassic event caused a significant reduction in osteichthyan body size. Neopterygii, the clade that encompasses the vast majority of extant fishes, underwent another diversification phase in the Late Triassic. The Triassic radiation of Osteichthyes, predominantly of Actinopterygii, which only occurred after severe extinctions among Chondrichthyes during the Middle–Late Permian, resulted in a profound change within global fish communities, from chondrichthyan‐rich faunas of the Permo‐Carboniferous to typical Mesozoic and Cenozoic associations dominated by actinopterygians. This turnover was not sudden but followed a stepwise pattern, with leaps during extinction events.     

Large extinctions of articulate brachiopods in the paleozoic and their ecological and evolutionary consequences

The largest Paleozoic extinctions of articulate brachiopods occurred at the Frasnian—Famennian boundary in the Late Devonian and at the Permian—Triassic boundary. Both extinctions affected taxa of all levels, including orders, but differed in scale, course, and ecological and evolutionary consequences. The Frasnian—Famennian extinction event was selective and evolutionary activity after the crisis varied in different orders. However, in the Early Carboniferous, the brachiopod diversity was mostly restored in comparison with the Devonian maximum. In particular groups, preadaptation played a role in changes in diversity and reconstruction of communities. The brachiopod composition at this boundary changed sharply. The extinction event at the end of Permian was global and accompanied by changes in the biota. Later, in the Meso-Cenozoic, the brachiopod diversity was not restored, and bivalves acquired primary importance in various bottom communities of different sea zones where Paleozoic brachiopods previously dominated. Extinction of brachiopods at this boundary was long and gradual. The symptoms of the ecological crisis in the development of Permian brachiopods are recognized beginning from the Roadian Age, which was probably the onset of this crisis.     

Cyclism revisited: extinction and ‘Achilles’ Heels' keep diversification in check on macroevolutionary time scales

Mass extinctions of varying magnitude prune the continuous diversification predicted by Darwinian evolutionary processes. They are caused by events that are too rare to become adaptatively accommodated. Their effects depend not only on the nature and magnitude of the triggering event but also on the state of the biosphere at the particular time. This is most clearly shown by the existence of Golden Ages preceding all Phanerozoic mass extinctions. These coincide with greenhouse periods, in which doomed clades gave rise to heteromorphs, deviating in strange ways from established bauplans. When critically examined, the seemingly ‘decadent’ morphologies of Schindewolf's ‘typolytic stages’ turn out to have been highly functional. The paradoxical link between adaptive peaks and evolutionary failure can now be explained. Specialisation tends to increase vulnerability (1) by narrowing niches and (2) by the retention of clade-specific conservative features that happen to become fatal Achilles’ Heels for entire clades in the face of a particular perturbation. Following extinctions, the availability of open niches favoured relatively rapid diversification of more innovative clades and their rise to ecological dominance (Schindewolf's ‘typogenetic stage’). Although the long-term changes can be observed only in the fossil record, Golden Biotopes in the present biosphere show that the Darwinian process may also be promoted by ecological isolation. As a result, clade histories do resemble individual biographies, but for ecological rather than orthogenetic reasons. This insight may help us to deal with the present mass extinction caused by our own species.     

FOSSILS AND A LARGE MOLECULAR PHYLOGENY SHOW THAT THE EVOLUTION OF SPECIES RICHNESS,GENERIC DIVERSITY,AND TURNOVER RATES ARE DISCONNECTED

The magnitude and extent of global change during the Cenozoic is remarkable, yet the impacts of these global changes on the biodiversity and evolutionary dynamics of species diversification remain poorly understood. To investigate this question, we combine paleontological and neontological data for the angiosperm order Fagales, an ecologically important clade of about 1370 species of trees with an exceptional fossil record. We show differences in patterns of accumulation of generic diversity, species richness, and turnover rates for Fagales. Generic diversity evolved rapidly since the Late Cretaceous and peaked during the Eocene or Oligocene. Turnover rates were high during periods of extreme global climate change, but relatively low when the climate remained stable. Species richness accumulated gradually throughout the Cenozoic, possibly at an accelerated pace after the Middle Miocene. Species diversification occurred in new environments: Quercoids radiating in Oligocene subtropical seasonally arid habitats, Casuarinaceae in Australian pyrophytic biomes, and Betula in Late Neogene holarctic habitats. These radiations were counterbalanced by regional extinctions in Late Neogene mesic warm‐temperate forests. Thus, the overall diversification at species level is linked to regional radiations of clades with appropriate ecologies exploiting newly available habitats.     

HOW DID LIFE BECOME SO DIVERSE? THE DYNAMICS OF DIVERSIFICATION ACCORDING TO THE FOSSIL RECORD AND MOLECULAR PHYLOGENETICS

Abstract: The long‐term diversification of life probably cannot be modelled as a simple equilibrial process: the time scales are too long, the potential for exploring new ecospace is too large and it is unlikely that ecological controls can act at global scales. The sum of many clade expansions and reductions, each of which happens according to its own dynamic, probably approximates more a damped exponential curve when translated into a global‐scale species diversification curve. Unfortunately, it is not possible to plot such a meaningful global‐scale species diversification curve through time, but curves at higher taxonomic levels have been produced. These curves are subject to the vagaries of the fossil record, but it is unlikely that the sources of error entirely overwhelm the biological signal. Clades radiate when the external and internal conditions are right: a new territory or ecospace becomes available, and the lineage has acquired a number of characters that open up a new diet or mode of life. Modern high levels of diversity in certain speciose clades may depend on such ancient opportunities taken. Dramatic climatic changes through the Quaternary must have driven extinctions and originations, but many species responded simply by moving to more favourable locations. Ecological communities appear to be no more than merely chance associations of species, but there may be real interactions among species. Ironically, high species diversity may lead to more speciation, not, as had been assumed, less: more species create more opportunities and selective pressures for other species to respond to, rather than capping diversity at a fixed equilibrium level. Studies from the scale of modern ecosystems to global long‐term patterns in the fossil record support a model for the exponential diversification of life, and one explanation for a pattern of exponential diversification is that as diversity increases, new forms become ever more refinements of existing forms. In a sense the world becomes increasingly divided into finer niche space. Organisms have a propensity to speciate freely, species richness within ecosystems appears to generate opportunities for more speciation, clades show all kinds of patterns from sluggish speciation rates and constant diversity through time to apparently explosive speciation, and there is no evidence that rapidly speciating clades have reached a limit, nor that they are driving other clades to extinction. A corollary of this view is that current biodiversity must be higher than it has ever been. Limits to infinite growth are clearly local, regional, and global turnover and extinction events, when climate change and physical catastrophes knock out species and whole clades, and push the rising exponential curve down a notch or two.     

The best sections method of studying mass extinctions

Phanerozoic mass extinctions have been studied primarily by analysing global diversity patterns compiled from the published literature. However, such compilations are beset by problems of incorrect correlation, imprecise age assignments and changing taxonomy. An alternative approach is to analyse mass extinctions by the ‘best sections’ method. This method identifies abundantly fossiliferous, well‐studied, stratigraphically dense and temporally extensive fossil records in strata that contain geochemical and other relevant non‐palaeontological data from a single depositional basin or geographically restricted outcrop area as the ‘best sections’ by which to analyse extinctions. A strength of the best sections method is that it allows the extinctions identified to be compared directly to changes in facies and other factors recorded in the best section. And, the hypothesis of a widespread extinction based on an extinction seen in a best section can be tested by its presence or absence in temporally equivalent sections. What we need are more field‐based studies of the best sections that encompass mass extinctions (real and hypothetical) and less of a reliance on literature‐based diversity compilations to produce a more reliable and comprehensive understanding of the history of extinctions.     

Were Phanerozoic mass extinctions among brachiopod superfamilies selective by taxa longevity?

All mass extinctions are characterized by certain kind of selectivity. An analysis of stratigraphic ranges of 112 brachiopod superfamilies implies that some Phanerozoic mass extinctions (Late Ordovician, Frasnian/Famennian and Devonian/Carboniferous, Early Jurassic, and Cretaceous/Paleogene) were selective by taxa longevity. They preferentially affected relatively old superfamilies and favoured a survival of relatively young superfamilies. No explanation of this selectivity as an apparent phenomenon is fully satisfactory. The Permian/Triassic mass extinction did not favour a survival of “young” superfamilies because of abnormally low rate of origination established since the Pennsylvanian and the absence of these “young” taxa. This study confirms tentatively a difference between Paleozoic and post-Paleozoic times by the importance of post-extinction recovery intervals for taxa longevity.     

Developmental stages of Permian bivalves of northeast Asia

In the development of Permian bivalves of northeastern Asia, the following five large stages have been recognized: Asselian-Middle Artinskian, Late Artinskian-Kungurian, Roadian-Wordian, Capitanian-Early Wuchiapingian, and Late Wuchiapingian-Changhsingian. The boundaries between stages correspond to great biotic events and frequently display a sharp change in species diversity of bivalves. During the Permian, brachiopods were gradually replaced by bivalves; this was promoted by repeated sharp changes in the environment.     

Shell mineralogical trends in epifaunal Mesozoic bivalves and their relationship to seawater chemistry and atmospheric carbon dioxide concentration

The Late Triassic-Early Jurassic change from aragonite- to calcite-facilitating conditions in the oceans, which was caused by a decrease of the Mg²⁺/Ca²⁺ ratio of seawater in combination with an increase of the partial pressure of carbon dioxide, also affected the shell mineralogy of epifaunal bivalves. In the “calcite sea” of the Jurassic and Cretaceous, the most diverse and abundant families of epifaunal bivalves had largely calcitic shells. Some of them, such as the Inoceramidae, acquired this shell mineralogy earlier in Earth's history but did not significantly diversify until the onset of “calcite sea” conditions. Others, however, replaced aragonite by calcite in their shell at the beginning of the Jurassic, as shown for the Ostreidae, Gryphaeidae, Pectinidae, Plicatulidae, and Buchiidae. In these families, replacement of aragonite by calcite took place in the middle and inner layer of the shell and was not associated with changes in morphology and life habit. It is therefore proposed that lower metabolic costs rather than higher resistance against dissolution or advantageous physical properties triggered the calcite expansion in their shells. This explanation fits well the observation that clades of thin-shelled bivalves were less affected by the change of seawater chemistry. Thick-shelled clades, by contrast, may suffer a severe decline in diversity until they adapt their shell mineralogy, as demonstrated by the Hippuritoida: The diversity of the Megalodontoidea, which failed to adapt their shell mineralogy to “calcite sea” conditions, dramatically decreased at the end of the Triassic, whereas their descendents became dominant carbonate producers during the Late Mesozoic after they acquired a calcitic outer shell layer in the Late Jurassic. These examples indicate that changes in the seawater chemistry and in the partial pressure of carbon dioxide are factors that influence the diversity of carbonate-secreting animals, and, as in the case of the decline of the Megalodontoidea, may contribute to mass extinctions.     

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.     

Evolution of Lycopodiaceae (Lycopsida): estimating divergence times from rbcL gene sequences by use of nonparametric rate smoothing

By use of nonparametric rate smoothing and nucleotide sequences of the rbcL gene, divergence times in Lycopodiaceae are estimated. The results show that much extant species diversity in Lycopodiaceae stems from relatively recent cladogenic events. These results corroborate previous ideas based on paleobotanical and biogeographical data. Previous molecular phylogenetic analyses recognized a split into neotropical and paleotropical clades in Huperzia, which contains 85-90% of all living species. Connecting this biogeographical pattern with continent movements, the diversification of this epiphytic group was suggested to coincide with that of angiosperms in the mid to Late Cretaceous. Results presented here are consistent with this idea, and the diversification of the two clades is resolved as Late Cretacous (78 and 95 Myr). In the related genera Lycopodium and Lycopodiella, the patterns are somewhat different. Here species diversity is scattered among different subgeneric groups. Most of the high-diversity subgeneric groups seem to have diversified very recently (Late Tertiary), whereas the cladogenic events leading to these groups are much older (Early to Late Cretaceous). Our analysis shows that, although much living species diversity stems from relatively recent cladogenesis, the origins of the family (Early Carboniferous) and generic crown groups (Early Permian to Early Jurassic) are much more ancient events.     

The radiation of cynodonts and the ground plan of mammalian morphological diversity

Cynodont therapsids diversified extensively after the Permo-Triassic mass extinction event, and gave rise to mammals in the Jurassic. We use an enlarged and revised dataset of discrete skeletal characters to build a new phylogeny for all main cynodont clades from the Late Permian to the Early Jurassic, and we analyse models of morphological diversification in the group. Basal taxa and epicynodonts are paraphyletic relative to eucynodonts, and the latter are divided into cynognathians and probainognathians, with tritylodonts and mammals forming sister groups. Disparity analyses reveal a heterogeneous distribution of cynodonts in a morphospace derived from cladistic characters. Pairwise morphological distances are weakly correlated with phylogenetic distances. Comparisons of disparity by groups and through time are non-significant, especially after the data are rarefied. A disparity peak occurs in the Early/Middle Triassic, after which period the mean disparity fluctuates little. Cynognathians were characterized by high evolutionary rates and high diversity early in their history, whereas probainognathian rates were low. Community structure may have been instrumental in imposing different rates on the two clades.     

Signatures of random and selective mass extinctions in phylogenetic tree balance

Current models of diversification with evolving speciation rates have trouble mimicking the extreme imbalance seen in estimated phylogenies. However, these models have not incorporated extinction. Here, we report on a simple simulation model that includes heritable and evolving speciation rates coupled with mass extinctions, Random (but not selective) mass extinctions, coupled with evolving among-lineage variation in speciation rates, increase imbalance of postrecovery clades. Thus, random mass extinctions are plausible contributors to the imbalance of modern clades. Paleontological evidence suggests that mass extinctions are often random with respect to ecological and morphological traits, consistent with our simulations. In contrast, evidence that the current anthropogenic mass extinction is phylogenetically selective suggests that the current extinction episode may be qualitatively different from past ones in the way it reshapes future biotas.     

Decoupling of morphological disparity and taxic diversity during the adaptive radiation of anomodont therapsids

Adaptive radiations are central to macroevolutionary theory. Whether triggered by acquisition of new traits or ecological opportunities arising from mass extinctions, it is debated whether adaptive radiations are marked by initial expansion of taxic diversity or of morphological disparity (the range of anatomical form). If a group rediversifies following a mass extinction, it is said to have passed through a macroevolutionary bottleneck, and the loss of taxic or phylogenetic diversity may limit the amount of morphological novelty that it can subsequently generate. Anomodont therapsids, a diverse clade of Permian and Triassic herbivorous tetrapods, passed through a bottleneck during the end-Permian mass extinction. Their taxic diversity increased during the Permian, declined significantly at the Permo–Triassic boundary and rebounded during the Middle Triassic before the clade''s final extinction at the end of the Triassic. By sharp contrast, disparity declined steadily during most of anomodont history. Our results highlight three main aspects of adaptive radiations: (i) diversity and disparity are generally decoupled; (ii) models of radiations following mass extinctions may differ from those triggered by other causes (e.g. trait acquisition); and (iii) the bottleneck caused by a mass extinction means that a clade can emerge lacking its original potential for generating morphological variety.     

Thecamoebians (Testate Amoebae) Straddling the Permian-Triassic Boundary in the Guryul Ravine Section,India: Evolutionary and Palaeoecological Implications

Exceptionally well-preserved organic remains of thecamoebians (testate amoebae) were preserved in marine sediments that straddle the greatest extinction event in the Phanerozoic: the Permian-Triassic Boundary. Outcrops from the Late Permian Zewan Formation and the Early Triassic Khunamuh Formation are represented by a complete sedimentary sequence at the Guryul Ravine Section in Kashmir, India, which is an archetypal Permian-Triassic boundary sequence [1]. Previous biostratigraphic analysis provides chronological control for the section, and a perspective of faunal turnover in the brachiopods, ammonoids, bivalves, conodonts, gastropods and foraminifera. Thecamoebians were concentrated from bulk sediments using palynological procedures, which isolated the organic constituents of preserved thecamoebian tests. The recovered individuals demonstrate exceptional similarity to the modern thecamoebian families Centropyxidae, Arcellidae, Hyalospheniidae and Trigonopyxidae, however, the vast majority belong to the Centropyxidae. This study further confirms the morphologic stability of the thecamoebian lineages through the Phanerozoic, and most importantly, their apparent little response to an infamous biological crisis in Earth’s history.     

Global Taxonomic Diversity of Anomodonts (Tetrapoda,Therapsida) and the Terrestrial Rock Record Across the Permian-Triassic Boundary

The end-Permian biotic crisis (∼252.5 Ma) represents the most severe extinction event in Earth''s history. This paper investigates diversity patterns in Anomodontia, an extinct group of therapsid synapsids (‘mammal-like reptiles’), through time and in particular across this event. As herbivores and the dominant terrestrial tetrapods of their time, anomodonts play a central role in assessing the impact of the end-Permian extinction on terrestrial ecosystems. Taxonomic diversity analysis reveals that anomodonts experienced three distinct phases of diversification interrupted by the same number of extinctions, i.e. an end-Guadalupian, an end-Permian, and a mid-Triassic extinction. A positive correlation between the number of taxa and the number of formations per time interval shows that anomodont diversity is biased by the Permian-Triassic terrestrial rock record. Normalized diversity curves indicate that anomodont richness continuously declines from the Middle Permian to the Late Triassic, but also reveals all three extinction events. Taxonomic rates (origination and extinction) indicate that the end-Guadalupian and end-Permian extinctions were driven by increased rates of extinction as well as low origination rates. However, this pattern is not evident at the final decline of anomodont diversity during the Middle Triassic. Therefore, it remains unclear whether the Middle Triassic extinction represents a gradual or abrupt event that is unique to anomodonts or more common among terrestrial tetrapods. The end-Permian extinction represents the most distinct event in terms of decline in anomodont richness and turnover rates.