< Back to 68k.news US front page

Evolution and dispersal of snakes across the Cretaceous-Paleogene mass extinction

Original source (on modern site)


Mass extinctions have repeatedly shaped global biodiversity. The Cretaceous-Paleogene (K-Pg) mass extinction caused the demise of numerous vertebrate groups, and its aftermath saw the rapid diversification of surviving mammals, birds, frogs, and teleost fishes. However, the effects of the K-Pg extinction on the evolution of snakes—a major clade of predators comprising over 3,700 living species—remains poorly understood. Here, we combine an extensive molecular dataset with phylogenetically and stratigraphically constrained fossil calibrations to infer an evolutionary timescale for Serpentes. We reveal a potential diversification among crown snakes associated with the K-Pg mass extinction, led by the successful colonisation of Asia by the major extant clade Afrophidia. Vertebral morphometrics suggest increasing morphological specialisation among marine snakes through the Paleogene. The dispersal patterns of snakes following the K-Pg underscore the importance of this mass extinction event in shaping Earth's extant vertebrate faunas.


The Cretaceous-Paleogene (K-Pg) transition resulted in the loss of an estimated 76% of all species1,2. High-resolution records of fossil pollen and marine microfossils show that the K-Pg extinction coincided with the Chicxulub bolide impact in Mexico3,4, which generated shockwaves, tsunamis, and a brief thermal pulse caused by the re-entry of heated ejecta in its immediate aftermath4,5. Critically, sulfates, dust, and hydrocarbon soot ejected into the atmosphere are thought to have greatly reduced insolation, resulting in global cooling over a period of years4. This would have caused a global collapse of photosynthesis4, resulting in high levels of extinction across the tree of life. Among vertebrates, iconic Mesozoic groups such as non-avian dinosaurs6, pterosaurs7, and mosasaurs8 disappeared at the K-Pg boundary. Other clades persisted but suffered severe reductions in diversity, including birds9, mammals2, and squamates10. Studies of insect feeding traces also indicate major losses among herbivorous arthropods11, whereas land plants suffered species-level extinctions of 80% or more12. Despite the magnitude of the impact and subsequent extinctions, terrestrial ecosystems recovered rapidly. In the aftermath, a "fern spike" followed by pioneer angiosperm communities documents the recovery of Earth's flora13,14, while surviving vertebrate groups, including mammals15,16, birds17,18,19, frogs20, and teleost fishes21,22, recovered and rapidly radiated in the early Cenozoic.

Historically, squamates were believed to have experienced minimal extinction at the K-Pg boundary23. However, analysis of the K-Pg transition in western North America found evidence for high rates of extinction among squamates10, although it remains unclear whether this pattern holds on a global scale. The evolutionary history of snakes across the K-Pg boundary has been particularly difficult to assess. The early fossil record of crown group snakes is fragmentary, often restricted to vertebrae and afflicted by relatively high rates of homoplasy24. As a result, phylogenetic analyses largely rely on restricted character sets, which is frequently misleading for vertebrates25,26,27,28, complicating our understanding of the affinities of many fossil snakes and possibly obscuring macroevolutionary patterns across the K-Pg boundary.

Molecular divergence dating efforts across Squamata suggest post-Cretaceous diversifications of major clades such as lacertids29 and amphisbaenians30. So far, molecular divergence time analyses of snakes recover conflicting patterns. Most studies31,32,33,34 suggest that the majority of extant snake clades diverged in the Cretaceous, although several analyses hint at a more recent diversification of the major subclade Alethinophidia30,35,36. Given this uncertainty, we attempted to improve our understanding of the timescale of crown snake diversification and the methodological factors that affect these inferences.

In this study, we investigate the effects of a plurality of modelling approaches—alternative calibration schemes (one developed here, in addition to that of Head et al.37,38, as well as variations on both), clock models (uncorrelated vs. autocorrelated), and associated priors (uniform, skew-T, and skew-Normal)—upon inferred node ages within crown snakes. Our results suggest a potential diversification of snakes near the time of the K-Pg transition. We also explore the effect of the K-Pg extinction event on vertebral morphological disparity, in order to directly incorporate an extensive sample of fossil taxa, while also benefiting from the tendency for morphological disparity indices to be less sensitive to sampling biases than diversity estimates39. We find a pattern of increasing vertebral disparity in the aftermath of the extinction, with concurrent increases in average and maximum body size, and dispersal to previously unoccupied landmasses.


Divergence time estimation

Our inferred tree topology is largely congruent with previous analyses32. Using our calibration set, we consistently recovered a diversification of crown Afrophidia (used here to refer to all alethinophidians, except Tropidophiidae + Aniliidae) near the K-Pg transition, and near-contemporaneous originations of the Amerophidia, Leptotyphlopidae, and Typhlopoidea crown groups, irrespective of the statistical framework, molecular clock model, priors, or phylogeny used (Fig. 1 and Supplementary Figs. 4-16). Implementing alternative assumptions regarding the quality of the fossil record by manipulating parameterizations of priors on fossil calibrations influenced mean estimated ages, but confidence intervals for nodes of interest—in particular for crown Afrophidia—converged near the K-Pg boundary. Choice of calibration priors influenced our estimated divergence times, reflecting the differential approaches to weighting fossil calibration ages inherent in the shapes of their probability distributions: skew-T yielded the youngest mean age estimates as it places the bulk of the effective prior distributions near the age of a fossil calibration, whereas uniform priors recovered ages averaging ~6.5 Ma older. Repeating divergence time analyses with different backbone topologies had little effect on our results (Supplementary Fig. S10, 11). The major differences between the results presented here and those of previous studies therefore appear to be a function of alternative calibration strategies (see Supplementary Information 2 for further discussion of our calibration approach), where the inclusion or exclusion of various Cretaceous fossils as internal calibrations for crown snakes influenced inferences on the timing of crown snake diversification with respect to the end-Cretaceous mass extinction. The conservative calibration approach we propose in Supplementary Information 2.2 excludes all putative records of Cretaceous crown snakes and we note that inferences on the role of the K-Pg transition in structuring the diversification of crown snakes discussed herein are dependent on this calibration scheme (see Supplementary Figs. 12-21 for exploration of alternative calibration strategies)—emphasizing the need to closely re-evaluate the phylogenetic placement of these fossils.

Fig. 1: Time-calibrated phylogeny of crown snakes based on the 42 fossil calibration set developed in this study, under a Bayesian Inference framework.

Green curves represent concatenated posterior age estimates for ingroup nodes using a skew-T and skew-Normal prior. Major clades discussed in the main text are shown on the right; Scolecoph., Scolecophidia; Am., Amerophidia. The K-Pg boundary is indicated by a red line and the Oligocene-Miocene boundary by a yellow line.

Disparity through time

Our analyses consider overall vertebral size independently from disparity of form, although we note that shape variables may scale allometrically. Shape disparity remained relatively constant throughout the Late Cretaceous, illustrating substantial morphological diversity early in the evolutionary history of total-clade snakes (Fig. 2 and see Supplementary Figs. 22-24 for additional axes, Supplementary Fig. 25 for labelled taxa). Disparity (sums of ranges and variances) increased across the K-Pg boundary, from the Maastrichtian to the Paleocene. Rarefaction analysis suggests that these differences were significant (the disparity of the smaller, Paleocene sample lies above the 95% confidence interval of the larger, Maastrichtian sample) and not a function of sample size differences (Supplementary Fig. 26). Higher Paleocene disparity was partly a function of a small number of morphologically eccentric taxa, notably Titanoboa40, an early palaeophiid41, and an indeterminate scolecophidian42. We note that the latter two records are late Paleocene to earliest Eocene in age—much of the apparent increase in vertebral disparity therefore appears to have been achieved after the first stage of the Paleocene. Sparse fossil sampling in the Danian means that the effects of the immediate aftermath of the K-Pg event are difficult to determine (although the very sparseness of the Danian record may point towards a decline in abundance following the extinction event43) and our approach, which necessarily involves time averaging and binning, is such that a short-lived crash in disparity following the K-Pg would be undetectable even if it were present. Moreover, empirical44 and simulation studies45 demonstrate that levels of disparity can be maintained in the face of even severe diversity loss, provided extinction is not centrifugal (i.e., preferentially concentrated away from the centre of morphospace).

Fig. 2: Principal component analysis of vertebral shape across different time bins, based on nine measurements, corrected for size.

The K-Pg boundary is indicated by a red line.

Vertebral morphological innovations in the Paleogene were driven by adaptation to marine life by palaeophiids, including extreme dorsoventral vertebral elongation. By the late Eocene, palaeophiids form a cluster distinct from all other snakes (Fig. 2). Rarefaction analyses controlling for sample size differences reveal that the Oligocene fauna was significantly less disparate than that of the late Eocene using sum of range and sum of variance (intrinsically less susceptible to sample size differences) disparity indices (Supplementary Fig. 26).

Postzygapophyseal width has been used to estimate body length of Titanoboa based on regression from a sample of boine snakes40. Using prezygapophyseal width as a proxy for body size (similar to postzygapophyseal width but allowing for a greater sample size and also strongly correlated with the mean of all measurements, as seen in Supplementary Fig. 29), we document a significant overall decrease in size from the mid-Cretaceous to the Campanian (Supplementary Fig. 28 and Supplementary Table 4). There was no significant change over the K-Pg boundary, although the Paleocene contains some particularly large taxa such as Titanoboa40 and Gigantophis46. In the late Eocene, the increasing dorsoventral vertebral elongation and overall size increase among Palaeophiidae results in a significant difference in size between larger aquatic taxa and smaller terrestrial taxa (see Supplementary Fig. 28 and Supplementary Table 4).

Biogeographical reconstruction

Of six biogeographic models tested in BioGeoBEARS47, +J models (allowing founder events) consistently outperformed their simpler counterparts, with DIVALIKE +J exhibiting the best fit to our data (Supplementary Information 4.3 and Supplementary Table 6). This result suggests that crown snake biogeographic history has largely been driven by dispersal, extinction, vicariance, and especially founder events, rather than widespread or subset sympatry (DEC and BAYAREAlike models)48. We recovered a strong signal for an Asian origin of the major clade Afrophidia (Fig. 3); our preferred model also recovered a South American origin for Amerophidia and suggested African or North American origins for Leptotyphlopidae and Indian origins for Typhlopoidea.

Fig. 3: Reconstructed historical biogeography of crown snakes, applying the DIVAlike + J model from BioGeoBEARS.

Pies represent the likelihood of the presence in each geographic area for each node; squares indicate present-day geographic distributions of extant taxa. The K-Pg boundary is indicated by a red line and the Oligocene-Miocene boundary by a yellow line.


Early snake evolution

Cretaceous snake fossils are predominantly, but not exclusively, found in Gondwanan deposits, in line with previous analyses, which suggest a Gondwanan origin for crown snakes34. Diverse assemblages have been noted from as early as the Cenomanian, such as the Moroccan Kem Kem beds, which include both aquatic and terrestrial taxa49,50,51. This is also reflected in morphological disparity, with early Late Cretaceous disparity rivalling that achieved later in the Cretaceous (Fig. 2).

Our molecular clock inferences suggest that as few as six extant lineages crossed the K-Pg boundary (Fig. 1 and Supplementary Fig. 4). However, deviations from our preferred calibration scheme find different results (see Supplementary Information 2) and the fossil record documents additional extinct clades (Nigerophiidae52, Madtsoiidae53, and Coniophiidae54), which also passed through the end-Cretaceous extinction event. The effects of the K-Pg mass extinction on snake disparity are unclear, largely because our data do not permit fine temporal resolution. We recover higher disparity estimates in the Paleocene than in the latest Cretaceous, although several scenarios could explain this pattern. First, it is consistent with the persistence of a wide range of morphologies (and presumably ecological niches) across the boundary, indicating that, overall, the disparity of snake vertebral morphology may not have been dramatically affected by the end-Cretaceous mass extinction. Fossorial adaptations of 'Scolecophidia' and early branching alethinophidians may have been key to their survival, as burrowing would have provided protection against many environmental disturbances caused by the asteroid impact, including significant and global shifts in temperature5. Similarly, freshwater ecosystems may have been buffered against extinction through the high thermal inertia of water during an initial thermal pulse4,5,55 and subsequent cooling, along with food chains based more on detritus and less on photosynthesis56,57, and may have provided a refugium for surviving aquatic nigerophiids. Surviving alethinophidians, madtsoiids, and coniophiids were probably macropredators34, but snake specializations for infrequent feeding58,59 may have facilitated their survival despite the limited availability of large prey in the aftermath of the asteroid impact60. The ability of many snakes to hunt in darkness may also have aided their search for food in the low light conditions that prevailed following the impact4. The eventual recovery of global forests could have provided new opportunities for snakes, as it did for birds61 and possibly frogs20, with the earliest known putatively arboreal snake, Corallus priscus, appearing in the Early Eocene62.

Alternatively, snakes could have experienced a more severe but short-lived reduction in disparity across the K-Pg boundary, with regions of morphospace temporarily vacated and subsequently re-occupied by the end of the Danian. Unfortunately, the rarity of early Paleocene snakes makes it impossible to test this hypothesis at present. Additional fossils from this key interval will be needed to better understand patterns of morphospace change immediately above and below the K-Pg boundary.

Post-extinction diversification

Fossil calibrations and their treatment have frequently been shown to have a profound effect on divergence time estimation63,64 and have therefore been a focus for improving molecular clock analyses65. Problematically, the snake fossil record is beset by a number of issues that complicate the development of unambiguous calibration schemes, relating primarily to the patchy and fragmentary nature of the snake fossil record. This has led to a number of disagreements regarding, e.g., the phylogenetic placement of fossils across the snake tree of life33,34,37,38; the rarity of distinctive, unambiguous synapomorphies in vertebral specimens49; and the absence of a comprehensive morphological phylogenetic framework with which to evaluate the affinities of problematic fossils66,67. These effects underlie differences in calibration choices between the scheme developed in this study and previously published calibration schemes37,38, such as the decision to remove Simoliophiidae and Australophis as calibrations due to their controversial phylogenetic placement (see Supplementary Information 2.1 for full discussion and Supplementary Figs. 4-21 to see exploration of different calibration approaches).

Alternative treatment of statistical probability distributions on fossil calibrations and chosen relaxed clock models similarly affect recovered divergence time estimates. We concatenated timetree results from two analyses applying different interpretations of the quality of the snake fossil record—using skew-T and skew-Normal priors on the fossil calibrations (see Fig. 1 and Supplementary Figs. 4 and 5). Autocorrelated clock models, in which related species are assumed to evolve at similar rates, recover younger estimated divergence times than uncorrelated clock models, where related lineages are not assumed to have similar rates of evolution. We favour the use of autocorrelated clock models on the basis that evidence for rate correlation can be detected in empirical datasets68 (see Supplementary Information 2.4 for full discussion).

Our inferred timescale for the evolution of crown Serpentes (following our conservative calibration scheme; see Supplementary Information 2.2) suggests a diversification of snakes spanning the K-Pg mass extinction. However, discerning the precise timing of divergence events near geologically instantaneous events such as the Chicxulub asteroid impact is challenging16,19. On both logical and theoretical grounds69, we consider it more likely for a diversification event to have occurred post-extinction, as this scenario does not imply the survival of numerous recently evolved and closely related taxa across one of the most severe mass extinction events in Earth's history. Moreover, post-extinction fragmentation of populations and reduction in competition for resources could all facilitate a rapid diversification, as has been posited for numerous diverse vertebrate clades, such as placental mammals15,16, neoavian birds18, and acanthomorph fishes21,22. We therefore interpret our results regarding the timing of divergence events in Afrophidia to be consistent with a post-extinction burst in diversification.

The aftermath of the extinction event may have provided an ecological release for snakes. Following an extinction, survivors can either reoccupy vacated niche space or expand into hitherto unoccupied regions of ecospace and morphospace (to the extent that the two are correlated44,70) made available by the extinction of competitors and predators. A pattern of post-K-Pg morphospace expansion is supported by our disparity analysis, in which snakes realize several new and disparate vertebral morphologies in the Paleogene (Fig. 2). Our analyses also highlight the increasing specialization of marine palaeophiids, which ultimately become morphologically distinct from all contemporaneous snakes by the late Eocene. In snakes, vertebral morphology can be indicative of ecological niche and overall body size40,71. The extinction of non-avian dinosaurs6 and other terrestrial predators may have enabled snakes to exploit relatively abundant small vertebrate prey in the early Cenozoic. Similarly, the K-Pg-associated extinction of marine reptiles and large teleosts may have facilitated the exploitation of marine ecosystems by large Paleocene marine snakes (Supplementary Fig. 28).

Our results suggest that the K-Pg mass extinction also influenced biogeographic patterns among snakes. We infer an initial dispersal by snakes into Asia near the K-Pg boundary (a landmass that they had apparently not previously occupied, see Supplementary Information 4.1), coinciding with the diversification of crown Afrophidia (Fig. 3). These patterns are consistent with the hypothesis that extinction events may help drive biogeographic patterns30. Although the probability of intercontinental dispersal should be unaffected by a mass extinction event, the extinction of potential competitors and predators may increase the likelihood that dispersed lineages ultimately become established on a newly colonized landmass.

Oligocene extinction and recovery

Significantly, our results indicate a second major snake extinction and recovery event in the Oligocene. We infer multiple dispersal and diversification events within caenophidian snakes beginning near the Oligocene-Miocene boundary and extending into the early Miocene (Figs. 1 and 3). We also infer a significant drop in morphospace occupation from the late Eocene to the Oligocene, supported further by the disappearance of early Cenozoic marine snakes such as Palaeophiidae and Gigantophis from the fossil record72 (Fig. 2). The Eocene-Oligocene event was associated with severe global cooling73 and may therefore have particularly affected poikilothermic taxa such as snakes. We hypothesize that caenophidians may have acquired their predominantly diurnal habits during this cooler interval as a result of being driven towards elevated daytime activity levels to take advantage of warmer daytime temperatures. Niche differentiation from their nocturnal antecedents may therefore have played a primary role underpinning caenophidian radiation and post-dispersal colonization during this interval34,74. Warming in the late Oligocene would have permitted the occupation of higher latitudes by poikilothermic snakes, facilitating their dispersal to the Americas via Beringia, as has been posited for birds75 (Fig. 3). Concurrent fragmentation of forests and the spread of newly appearing grasslands76 may also have provided opportunities for speciation and adaptation to changing environments (Fig. 1).

It is clear that integrating molecular sequence data and the fossil record can elucidate the influence of major events in Earth's history on the evolutionary history of extant clades19. Our analyses suggest that the K-Pg mass extinction influenced the evolutionary history of snakes, potentially facilitating the diversification of Afrophidia and the origination of several other snake subclades in the extinction's aftermath. This appears to have enabled the successful establishment of snakes in Asia and provided the opportunity for Palaeophiidae to specialize as marine predators in the Eocene. Along with the contemporaneous diversification and disparification of numerous other vertebrate clades16,17,18,19,20,21,22, our results help corroborate the fundamental role of the K-Pg mass extinction in shaping the vertebrate biodiversity occupying our planet today.


Phylogenetic analyses

Analyses included 115 extant snake taxa. Forty-four non-snake squamates across all major squamate clades were sampled as outgroups and to provide additional calibration points. Ten non-squamate amniote outgroups were included to calibrate deeper nodes (see Supplementary Fig. 3). Taxa were sampled from the Zheng and Wiens32,77 supermatrix. Criteria for inclusion were maximizing representation of extant families, minimizing missing data within the phylogenetic data matrix, and prioritizing slower evolving taxa. All extant snake families, except Anomochilidae and Xenophidiidae, were included, as their sparse gene sampling reduced resolution.

We extracted molecular sequence alignments for our chosen taxa from the Zheng and Wiens32 supermatrix, which includes 40 nuclear and 12 mitochondrial loci. Overall matrix completeness was 49.7%. We ran phylogenetic analyses in PhyloBayes v.4.1c78 for ~30,000 generations, under the CAT + GTR + G parameters and with a birth-death prior, as PhyloBayes does not co-estimate tree topology and divergence time, and is MPI (Message Passing Interface) enabled79. We built a consensus tree, with a burn-in of 7500 generations and retaining every tenth tree. For analyses where different tree topologies were enforced, we used Mesquite v. 3.1080 to manipulate relationships in the consensus tree, to set 'Scolecophidia' as monophyletic, and to change immediate outgroup to be Iguania + Anguimorpha, as in Reeder et al.81.


We compiled a set of 42 calibrations, comprising previously published and novel calibrations (see Supplementary Information 2.2). Multiple variations of this calibration set were also tested, including the addition of Simoliophiidae as a minimum calibration for Alethinophidia32,33,34,38, as well as testing parviraptorids as a calibration for total-group snakes (see Supplementary Table 3 for summary). In addition, the calibration schemes proposed by Head38 and Head et al.37 were also tested, as well as a revision of this set with two amendments: the removal of Simoliophiidae to assess the effect of assigning them to crown Alethinophidia and the removal of Australophis (see Supplementary Information 2.5).

Estimating divergence times

To investigate the effects of different probabilistic methods, different relaxed clock models and different prior distributions of our fossil calibrations on divergence time estimates, we performed analyses in the PAML package MCMCTree82. For these analyses, data were partitioned by gene83,84. Prior distribution shapes and scales for calibrations were calculated using the MCMCTreeR package85. We used soft maximum bounds, with a tail of 0.05, and a root prior for Amniota of 318-332.9 Ma86. Divergence time estimates were calculated under Skew-T, Skew-Normal, and Uniform distributions, to simulate progressively more liberal interpretations of the fossil record83. Further analyses were run using the uniform and/or skew-T prior distribution with both independent rates and correlated rates, alternative topologies, and with alternative calibration schemes (see Supplementary Figs. 4-21). Analyses performed under an agnostic distribution (the uniform prior) are important. These make the fewest assumptions regarding the quality of the fossil record and are modelled with a hard minimum age. If our preferred set of assumptions were to be found invalid, uniform priors help to circumscribe the range of inferences that might reasonably be made from our data. However, given the fossil record of snakes significantly improved in the Cenozoic, where most of our calibrations are found, there is no reason to prefer an agnostic distribution over other priors; for Fig. 1, results from skew-T and skew-Normal priors were concatenated. Divergence time analyses were run for 100,000 generations, with a burn-in of 25,000 and sampling every 25th remaining tree.

Biogeographical reconstruction

Historical biogeography was reconstructed using BioGeoBEARS47. Geographical data were collected for ingroup taxa from the Reptile Database87. Full methods can be found in Supplementary Information 1 and input files in Supplementary Information 4.1-3. We assessed model fit using log likelihood, AIC (Akaike Information Criterion), and AICc (AIC with correction for small data samples) values.

Disparity analyses

Nine linear measurements were taken between geometrically homologous landmarks on fossil snake vertebrae using the software ImageJ88. As some vertebrae were incomplete, missing values were inferred using multiple imputation. For this we used the missRanger89 package in R, which implements the MissForest chaining random forests algorithm of Stekhoven and Bühlmann90 with predictive mean matching and a maximum of 50 chaining iterations. To control for differences in specimen size, all measurements were scaled to the mean of all measurements for said specimen. The data were then subjected to R-mode principal component analyses using a correlation matrix. All resultant PC axes were used to calculate sums of ranges and sums of variances as indices of disparity. In order to account for the effects of sample size variation, we jackknifed both indices of disparity at all possible sample sizes using Rare91 to n = 9 (smallest time bin). Rarefied indices of disparity were then plotted against logged sample size. Changes in size through time were independently assessed in R92, based on the raw prezygapophyseal width data. Kruskal-Wallis tests were used to test for significant changes in size between time bins and between inferred ecologies in the same time bin. For specimens with age uncertainty, entries were duplicated so as to be present in both relevant time bins. Overall, resultant sample sizes were n= 14 for the early Late Cretaceous, n= 15 for the Campanian, n= 19 for the Maastrichtian, n= 9 for the Paleocene, n= 48 for the Early Eocene, n= 33 for the late Eocene, and n= 30 for the Oligocene.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

Input and output files for molecular clock and disparity analyses are deposited at https://doi.org/10.5061/dryad.tv055. Information on extant taxon ranges was gathered from the Reptile Database https://reptile-database.reptarium.cz/. Gene sequences from Zheng and Wiens32,77 can be found at https://doi.org/10.5061/dryad.tv055. All other data can be found in the Supplementary Information. Different formats can be requested from corresponding authors.


  1. 1.

    Barnosky, A. D. et al. Has the Earth's sixth mass extinction already arrived? Nature 471, 51-57 (2011).

    ADS  CAS  PubMed  Article  Google Scholar 

  2. 2.

    Longrich, N. R., Scriberas, J. & Wills, M. A. Severe extinction and rapid recovery of mammals across the Cretaceous-Palaeogene boundary, and the effects of rarity on patterns of extinction and recovery. J. Evol. Biol. 29, 1495-1512 (2016).

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Alvarez, L. W., Alvarez, W., Asaro, F. & Michel, H. V. Extraterrestrial cause for the Cretaceous-Tertiary extinction. Science 208, 1095-1108 (1980).

    ADS  CAS  PubMed  Article  Google Scholar 

  4. 4.

    Schulte, P. et al. The Chicxulub asteroid impact and mass extinction at the Cretaceous-Paleogene boundary. Science 327, 1214-1218 (2010).

    ADS  CAS  PubMed  Article  Google Scholar 

  5. 5.

    Robertson, D. S., McKenna, M. C., Toon, O. B., Hope, S. & Lillegraven, J. A. Survival in the first hours of the Cenozoic. GSA Bull. 116, 760-768 (2004).

    Article  Google Scholar 

  6. 6.

    Brusatte, S. L. et al. The extinction of the dinosaurs. Biol. Rev. 90, 628-642 (2015).

    PubMed  Article  Google Scholar 

  7. 7.

    Longrich, N. R., Martill, D. M. & Andres, B. Late Maastrichtian pterosaurs from North Africa and mass extinction of Pterosauria at the Cretaceous-Paleogene boundary. PLoS Biol. 16, e2001663 (2018).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  8. 8.

    Polcyn, M. J., Jacobs, L. L., Araújo, R., Schulp, A. S. & Mateus, O. Physical drivers of mosasaur evolution. Palaeogeogr. Palaeoclimatol. Palaeoecol. 400, 17-27 (2014).

    Article  Google Scholar 

  9. 9.

    Longrich, N. R., Tokaryk, T. & Field, D. J. Mass extinction of birds at the Cretaceous - Paleogene (K-Pg) boundary. Proc. Natl Acad. Sci. USA 108, 15253-15257 (2011).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Longrich, N. R., Bhullar, B.-A. S. & Gauthier, J. A. Mass extinction of lizards and snakes at the Cretaceous-Paleogene boundary. Proc. Natl Acad. Sci. USA 109, 21396-21401 (2012).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Labandeira, C. C., Johnson, K. R. & Wilf, P. Impact of the terminal Cretaceous event on plant - insect associations. Proc. Natl Acad. Sci. USA 99, 2061-2066 (2001).

    ADS  Article  CAS  Google Scholar 

  12. 12.

    Wilf, P. & Johnson, K. R. Land plant extinction at the end of the Cretaceous: a quantitative analysis of the North Dakota megafloral record. Paleobiology 30, 347-368 (2004).

    Article  Google Scholar 

  13. 13.

    Wolfe, J. A. & Upchurch, G. R. Jr Vegetation, climatic and floral changes at the Cretaceous-Tertiary boundary. Nature 324, 148-152 (1986).

    ADS  Article  Google Scholar 

  14. 14.

    Nichols, D. J. & Johnson, K. R. Plants and the K-T Boundary (Cambridge Univ. Press, 2008).

  15. 15.

    Alroy, J. The fossil record of North American mammals: evidence for a Paleocene evolutionary radiation. Syst. Biol. 48, 107-118 (1999).

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    dos Reis, M. et al. Phylogenomic datasets provide both precision and accuracy in estimating the timescale of placental mammal phylogeny. Proc. R. Soc. B Biol. Sci. 279, 3491-3500 (2012).

    Article  Google Scholar 

  17. 17.

    Ksepka, D. T., Stidham, T. A. & Williamson, T. E. Early Paleocene landbird supports rapid phylogenetic and morphological diversification of crown birds after the K-Pg mass extinction. Proc. Natl Acad. Sci. USA 114, 8047-8052 (2017).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Prum, R. O. et al. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature 526, 569-573 (2015).

    ADS  CAS  PubMed  Article  Google Scholar 

  19. 19.

    Berv, J. S. & Field, D. J. Genomic signature of an avian Lilliput Effect across the K-Pg extinction. Syst. Biol. 67, 1-13 (2018).

    PubMed  Article  Google Scholar 

  20. 20.

    Feng, Y. et al. Phylogenomics reveals rapid, simultaneous diversification of three major clades of Gondwanan frogs at the Cretaceous - Paleogene boundary. Proc. Natl Acad. Sci. USA 114, E5864-E5870 (2017).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Friedman, M. Explosive morphological diversification of spiny-finned teleost fishes in the aftermath of the end-Cretaceous extinction. Proc. Biol. Sci. 277, 1675-1683 (2010).

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Alfaro, M. E. et al. Explosive diversification of marine fishes at the Cretaceous-Palaeogene boundary. Nat. Ecol. Evol. 2, 688-696 (2018).

    PubMed  Article  Google Scholar 

  23. 23.

    Evans, S. E. At the feet of the dinosaurs: the early history and radiation of lizards. Biol. Rev. Camb. Philos. Soc. 78, 513-551 (2003).

    PubMed  Article  Google Scholar 

  24. 24.

    Rage, J.-C. & Escuillié, F. The Cenomanian: Stage of Hindlimbed Snakes. Carnets de Géologie 1-11 (2003).

  25. 25.

    Mounce, R. C. P., Sansom, R. & Wills, M. A. Sampling diverse characters improves phylogenies: Craniodental and postcranial characters of vertebrates often imply different trees. Evolution 70, 666-668 (2016).

    PubMed  Article  Google Scholar 

  26. 26.

    Sansom, R. S., Wills, M. A. & Williams, T. Dental data perform relatively poorly in reconstructing mammal phylogenies: morphological partitions evaluated with molecular benchmarks. Syst. Biol. 66, 813-822 (2017).

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Li, Y., Ruta, M. & Wills, M. A. Craniodental and postcranial characters of non-avian Dinosauria often imply different trees. Syst. Biol. 69, 638-659 (2020).

  28. 28.

    Sansom, R. S. & Wills, M. A. Differences between hard and soft phylogenetic data. Proc. R. Soc. B 284, 20172150 (2017).

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Hipsley, C. A., Himmelmann, L., Metzler, D. & Müller, J. Integration of Bayesian molecular clock methods and fossil-based soft bounds reveals early Cenozoic origin of African lacertid lizards. BMC Evol. Biol. 9, 151 (2009).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Longrich, N. R., Vinther, J., Pyron, R. A., Pisani, D. & Gauthier, J. A. Biogeography of worm lizards (Amphisbaenia) driven by end-Cretaceous mass extinction. Proc. R. Soc. Lond. B Biol. Sci. 282, 20143034 (2015).

    Google Scholar 

  31. 31.

    Vidal, N. & Hedges, S. B. The molecular evolutionary tree of lizards, snakes, and amphisbaenians. C. R. Biol. 332, 129-139 (2009).

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Zheng, Y. & Wiens, J. J. Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species. Mol. Phylogenet. Evol. 94, 537-547 (2016).

    PubMed  Article  Google Scholar 

  33. 33.

    Harrington, S. M. & Reeder, T. W. Phylogenetic inference and divergence dating of snakes using molecules, morphology and fossils: new insights into convergent evolution of feeding morphology and limb reduction. Biol. J. Linn. Soc. 121, 379-394 (2017).

    Article  Google Scholar 

  34. 34.

    Hsiang, A. Y. et al. The origin of snakes: revealing the ecology, behavior, and evolutionary history of early snakes using genomics, phenomics, and the fossil record. BMC Evol. Biol. 15, 87 (2015).

    PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Jones, M. E. H. et al. Integration of molecules and new fossils supports a Triassic origin for Lepidosauria (lizards, snakes, and tuatara). BMC Evol. Biol. 13, 208 (2013).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Pyron, R. A. Novel approaches for phylogenetic inference from morphological data and total-evidence dating in squamate reptiles (lizards, snakes, and amphisbaenians). Syst. Biol. 66, 38-56 (2017).

    PubMed  Google Scholar 

  37. 37.

    Head, J. J., Mahlow, K. & Müller, J. Fossil calibration dates for molecular phylogenetic analysis of snakes 2: Caenophidia, Colubroidea, Elapoidea, Colubridae. Palaeontol. Electron. 19, 1-21 (2016).

    Article  Google Scholar 

  38. 38.

    Head, J. J. Fossil calibration dates for molecular phylogenetic analysis of snakes 1: Serpentes, Alethinophidia, Boidae, Pythonidae. Palaeontol. Electron. 18, 1-17 (2015).

    Google Scholar 

  39. 39.

    Butler, R. J., Brusatte, S. L., Andres, B. & Benson, R. B. J. How do geological sampling biases affect studies of morphological evolution in deep time? A case study of pterosaur (Reptilia: Archosauria) disparity. Evolution 66, 147-162 (2011).

    PubMed  Article  Google Scholar 

  40. 40.

    Head, J. J. et al. Giant boid snake from the Palaeocene neotropics reveals hotter past equatorial temperatures. Nature 457, 715-717 (2009).

    ADS  CAS  PubMed  Article  Google Scholar 

  41. 41.

    Kristensen, H. V., Cuny, G., Rasmussen, A. R. & Madsen, H. Earliest record of the fossil snake Palaeophis from the Paleocene/Eocene boundary in Denmark. Bull. Soc. Géol. Fr. 183, 621-625 (2012).

    Article  Google Scholar 

  42. 42.

    Augé, M. & Rage, J.-C. Herpetofaunas from the Upper Paleocene and Lower Eocene of Morocco. Ann. Paléontologie 92, 235-253 (2006).

    Article  Google Scholar 

  43. 43.

    Field, D. J. et al. Timing the extant avian radiation: The rise of modern birds, and the importance of modeling molecular rate variation. Bull. Am. Mus. Nat. Hist. 440, 1590191 (2020).

  44. 44.

    Korn, D., Hopkins, M. J. & Walton, S. A. Extinction space - a method for the quantification and classification of changes in morphospace across extinction boundaries. Evolution 67, 2795-2810 (2013).

    PubMed  Google Scholar 

  45. 45.

    Puttick, M. N., Guillerme, T. & Wills, M. A. The complex effects of mass extinctions on morphological disparity. Evolution 74, 2207-2220 (2020).

    PubMed  Article  Google Scholar 

  46. 46.

    Rio, J. P. & Mannion, P. D. The osteology of the giant snake Gigantophis garstini from the upper Eocene of North Africa and its bearing on the phylogenetic relationships and biogeography of Madtsoiidae. J. Vertebr. Paleontol. 37, e1347179 (2017).

    Article  Google Scholar 

  47. 47.

    Matzke, N. Founder-event speciation in BioGeoBEARS package dramatically improves likelihoods and alters parameter inference in Dispersal-Extinction-Cladogenesis (DEC) analyses. Front. Biogeogr. 4, 210 (2012).

    Google Scholar 

  48. 48.

    Matzke, N. Probabilistic historical biogeography: new models for founder-event speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Front. Biogeogr. 5, 242-248 (2013).

    Article  Google Scholar 

  49. 49.

    Klein, C. G., Longrich, N. R., Ibrahim, N., Zouhri, S. & Martill, D. M. A new basal snake from the mid-Cretaceous of Morocco. Cretac. Res. 72, 134-141 (2017).

    Article  Google Scholar 

  50. 50.

    Rage, J.-C. & Dutheil, D. B. Amphibians and squamates from the Cretaceous (Cenomanian) of Morocco. A preliminary study, with description of a new genus of pipid frog. Palaeontogr. Palaozool. Stratigraph. 285, 1-22 (2008).

    Google Scholar 

  51. 51.

    Vullo, R. A new species of Lapparentophis from the mid-Cretaceous Kem Kem beds, Morocco, with remarks on the distribution of lapparentophiid snakes. C. R. Palevol. 1132, 6-11 (2019).

    Google Scholar 

  52. 52.

    Rage, J.-C. in Handbuch der Paläoherpetologie (ed. Wellnhofer, P.) 1-80 (Gustav Fischer, 1984).

  53. 53.

    Rage, J.-C. Fossil snakes from the Palaeocene of São José de Itaboraí, Brazil. Part I. Madtsoiidae, Aniliidae. Palaeovertebrata 27, 109-1144 (1998).

    Google Scholar 

  54. 54.

    Smith, K. T. New constraints on the evolution of the snake clades Ungaliophiinae, Loxocemidae and Colubridae (Serpentes), with comments on the fossil history of erycine boids in North America. Zool. Anz. 252, 157-182 (2013).

    Article  Google Scholar 

  55. 55.

    Goldin, T. J. & Melosh, H. J. Self-shielding of thermal radiation by Chicxulub impact ejecta: firestorm or fizzle? Geology 37, 1135-1138 (2009).

    ADS  Article  Google Scholar 

  56. 56.

    Robertson, D. S., Lewis, W. M., Sheehan, P. M. & Toon, O. B. K-Pg extinction patterns in marine and freshwater environments: the impact winter model. J. Geophys. Res. Biogeosci. 118, 1006-1014 (2013).

    Article  Google Scholar 

  57. 57.

    Sheehan, P. M. & Fastovsky, D. E. Major extinctions of land-dwelling vertebrates at the Cretaceous-Tertiary boundary, eastern Montana. Geology 20, 556-560 (1992).

    ADS  Article  Google Scholar 

  58. 58.

    Greene, H. W. Snakes: The Evolution of Mystery in Nature (Univ. California Press, 1997).

  59. 59.

    Secor, S. M. & Diamond, J. M. Evolution of regulatory responses to feeding in snakes. Physiol. Biochem. Zool. 73, 123-141 (2000).

    CAS  PubMed  Article  Google Scholar 

  60. 60.

    Wilson, G. P. Mammals across the K/Pg boundary in northeastern Montana, U.S.A.: dental morphology and body-size patterns reveal extinction selectivity and immigrant-fueled ecospace filling. Paleobiology 39, 429-469 (2013).

    Article  Google Scholar 

  61. 61.

    Field, D. J. et al. Early evolution of modern birds structured by global forest collapse at the End-Cretaceous mass extinction. Curr. Biol. 28, 1825-1831.e2 (2018).

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Rage, J.-C. Fossil snakes from the Palaeocene of São José de Itaboraí, Brazil. Part II. Boidae. Palaeovertebrata (Montp.) 30, 111-150 (2001).

    Google Scholar 

  63. 63.

    Duchêne, S., Lanfear, R. & Ho, S. Y. W. The impact of calibration and clock-model choice on molecular estimates of divergence times. Mol. Phylogenet. Evol. 78, 277-289 (2014).

    PubMed  Article  Google Scholar 

  64. 64.

    Warnock, R. C. M., Parham, J. F., Joyce, W. G., Lyson, T. R. & Donoghue, P. C. J. Calibration uncertainty in molecular dating analyses: there is no substitute for the prior evaluation of time priors. Proc. R. Soc. B 282, 20141013 (2015).

    PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Parham, J. F. et al. Best practices for justifying fossil calibrations. Syst. Biol. 61, 346-359 (2012).

    PubMed  Article  Google Scholar 

  66. 66.

    Gauthier, J. A., Kearney, M., Maisano, J. A., Rieppel, O. & Behlke, A. D. B. Assembling the squamate tree of life: perspectives from the phenotype and the fossil record. Bull. Peabody Mus. Nat. Hist. 53, 3-308 (2012).

    Article  Google Scholar 

  67. 67.

    Conrad, J. L. Phylogeny and systematics of Squamata (Reptilia) based on morphology. Bull. Am. Mus. Nat. Hist. 310, 1-182 (2008).

    Article  Google Scholar 

  68. 68.

    Tao, Q., Tamura, K., Battistuzzi, F. U. & Kumar, S. A machine learning method for detecting autocorrelation of evolutionary rates in large phylogenies. Mol. Biol. Evol. 36, 811-824 (2019).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  69. 69.

    Budd, G. E. & Mann, R. P. The dynamics of stem and crown groups. Sci. Adv. 6, eaaz1626 (2020).

    ADS  PubMed  PubMed Central  Article  Google Scholar 

  70. 70.

    Hughes, M., Gerber, S. & Wills, M. A. Clades reach highest morphological disparity early in their evolution. Proc. Natl Acad. Sci. USA 110, 13875-13879 (2013).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Lawing, A. M., Head, J. J. & Polly, P. D. In Paleontology in Ecology and Conservation (ed., Louys, J.) 117-146 (Springer Earth System Sciences, 2012). https://doi.org/10.1007/978-3-642-25038-5

  72. 72.

    Rage, J.-C. & Augé, M. Squamates from the Cainozoic of the Western part of Europe. A review. Rev. Paléobiol. Spec. 7, 199-216 (1993).

    Google Scholar 

  73. 73.

    Liu, Z. et al. Eocene-Oligocene climate transition. Science 323, 1187-1190 (2008).

    ADS  Article  CAS  Google Scholar 

  74. 74.

    Simões, B. F. et al. Visual system evolution and the nature of the ancestral snake. J. Evol. Biol. 28, 1309-1320 (2015).

    PubMed  Article  CAS  Google Scholar 

  75. 75.

    Oliveros, C. H. et al. Earth history and the passerine superradiation. Proc. Natl Acad. Sci. 116, 7916-7925 (2019).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Strömberg, C. A. E. Evolution of grasses and grassland ecosystems. Annu. Rev. Earth Planet. Sci. USA 39, 517-544 (2011).

    ADS  Article  CAS  Google Scholar 

  77. 77.

    Zheng, Y. & Wiens, J. J. Data from: Combining phylogenomic and supermatrix approaches, and a time-calibrated phylogeny for squamate reptiles (lizards and snakes) based on 52 genes and 4162 species. DataDryad. https://doi.org/10.5061/dryad.tv055 (2016).

  78. 78.

    Lartillot, N., Lepage, T. & Blanquart, S. PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics 25, 2286-2288 (2009).

  79. 79.

    dos Reis, M., Donoghue, P. C. J. & Yang, Z. Bayesian molecular clock dating of species divergences in the genomics era. Nat. Rev. 17, 71-80 (2015).

    Article  CAS  Google Scholar 

  80. 80.

    Maddison, W. P. & Maddison, D. R. Mesquite: a modular system for evolutionary analysis. http://www.mesquiteproject.org (2016).

  81. 81.

    Reeder, T. W. et al. Integrated analyses resolve conflicts over squamate reptile phylogeny and reveal unexpected placements for fossil taxa. PLoS ONE 10, 1-22 (2015).

    Google Scholar 

  82. 82.

    Yang, Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    Betts, H. C. et al. Integrated genomic and fossil evidence illuminates life's early evolution and eukaryote origin. Nat. Ecol. Evol. 2, 1556-1562 (2018).

  84. 84.

    dos Reis, M. et al. Uncertainty in the timing of origin of animals and the limits of precision in molecular timescales. Curr. Biol. 25, 2939-2950 (2015).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  85. 85.

    Morris, J. L. et al. The timescale of early land plant evolution. Proc. Natl Acad. Sci. USA 115, E2274-E2283 (2018).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  86. 86.

    Benton, M. J. et al. Constraints on the timescale of animal evolutionary history. Palaeontol. Electron. 18, 1-106 (2015).

    Google Scholar 

  87. 87.

    Uetz, P., Freed, P. & Hosek, J. The Reptile Database http://www.reptile-database.org/ (2020).

  88. 88.

    Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671-675 (2012).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Mayer, M. missRanger: fast imputation of missing values. R package version 2.1.0 https://github.com/mayer79/missRanger (2019).

  90. 90.

    Stekhoven, D. J. & Bühlmann, P. MissForest — non-parametric missing value imputation for mixed-type data. Bioinformatics 25, 303-317 (2011).

    Google Scholar 

  91. 91.

    Wills, M. A. Crustacean disparity through the Phanerozoic: comparing morphological and stratigraphic data. Biol. J. Linn. Soc. 65, 455-500 (1998).

    Article  Google Scholar 

  92. 92.

    R CoreTeam. R: A Language and Environment for Statistical Computing (R CoreTeam, 2017).

Download references


This project was funded by the Fonds National de la Recherche Luxembourg (AFR grant 10142968 to C.G.K.) and the National Environment Research Council (NE/L002434/1 to R.J.L. and NE/K014951/1 to M.A.W.). D.J.F. is supported by a UK Research and Innovation Future Leaders Fellowship (MR/S032177/1). M.A.W. is supported by John Templeton Foundation grant 61408. We thank Jakob Vinther and Jason Head for helpful discussions during the conceptualization of this project. We also thank Mark Puttick and Aubrey Roberts for support with analyses.

Author information

Author notes

  1. Catherine G. Klein

    Present address: GeoZentrum Nordbayern, Friedrich-Alexander University Erlangen-Nürnberg, Loewenichstr. 28, Erlangen, Germany


  1. The Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Bath, UK

    Catherine G. Klein, Daniel J. Field, Rebecca Lakin, Matthew A. Wills & Nicholas R. Longrich

  2. School of Earth Sciences, School of Biological Sciences, University of Bristol, Bristol, UK

    Davide Pisani

  3. Department of Earth Sciences, University of Cambridge, Cambridge, Cambridgeshire, UK

    Daniel J. Field


N.R.L. conceived and designed the study, with input from D.P., D.J.F., M.A.W., and C.G.K. C.G.K. realized the study, collated and analysed molecular and morphological data, and implemented all phylogenetic and clock analyses with advice from D.P. N.R.L. designed the fossil calibration scheme. R.J.L. collected additional morphological data. M.A.W. performed morphological and disparity analyses. C.G.K. and N.R.L. wrote the paper, with contributions from D.J.F., D.P., and M.A.W. All authors approved the final manuscript.

Corresponding authors

Correspondence to Catherine G. Klein or Nicholas R. Longrich.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Tod Reeder and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Cite this article

Klein, C.G., Pisani, D., Field, D.J. et al. Evolution and dispersal of snakes across the Cretaceous-Paleogene mass extinction. Nat Commun 12, 5335 (2021). https://doi.org/10.1038/s41467-021-25136-y

Download citation


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

< Back to 68k.news US front page