Data
For our analysis of jawed vertebrate evolutionary history, we included eight major monophyletic clades—amphibians; birds; crocodylians; lizards, snakes and the Tuatara (Lepidosauria); mammals; ray-finned fish (Actinopterygii); sharks, rays and chimaeras (Chondrichthyes); and turtles and tortoises (Testudines)—which have strong phylogenetic and extinction risk data coverage and together account for >99% of jawed vertebrate species. For amphibian taxonomy, we used Frost’s Amphibian Species of the World 76, which included 8,024 species. We matched this taxonomy to the phylogeny of Jetz and Pyron22, for which we retained 7,002 species. For birds we used BirdLife International’s Handbook of the Birds of the World v5.0 77, identifying 10,988 valid and extant species, and matched the taxonomy to the phylogeny of Jetz et al.19, retaining 9,645 species. For mammals we used the mammal taxonomy and phylogenetic trees from Gumbs et al.35, which adopted the Mammal Diversity Database v1.1 78 of the American Society of Mammologists, identifying 6,253 extant and valid mammals species, and matched the taxonomy to the phylogeny of Upham et al.79, retaining 5,853 species. For chondrichthyans we used Fishbase80, which included 1,290 species, and matched the taxonomy to the phylogeny of Stein et al.24, retaining 1,165 species. For ray-finned fish we also used Fishbase80, which comprised 32,760 species, and matched this taxonomy to the phylogeny of Chang et al.81, retaining 31,506 species. For crocodylians, testudines and lepidosaurs, we used the Reptile Database82, which included 25, 351, and 10,735 species, respectively. We matched the reptilian taxonomy to the phylogeny of Colston et al.23 for crocodylians and testudines, retaining 25 crocodylian and 325 testudine species, and to the phylogeny of Tonini et al.21 for lepidosaurs, retaining 9,599 species.
For all clades except mammals we randomly sampled 1,000 trees from the available distributions to adequately capture uncertainty in phylogenetic relationships and node ages 19,83. To maximise the inclusion of species, and enable effective comparison between clades, we imputed species missing from the phylogenies, following taxonomy matching, to generate 1,000 phylogenetic trees for each clade that comprised all valid species present in each taxonomic treatment (1,022 amphibians, 1,343 birds, 1,691 ray-finned fish, 135 chondrichthyans, 26 testudines and tortoises, and 1,136 lepidosaurs). For imputation we followed earlier approaches to insert missing species into their genus along the existing phylogenetic branches25,35,84, using the ‘congeneric.impute’ function in the R package ‘pez’ 85. This approach has been used to provide estimates of threatened phylogenetic diversity 84 and was used to generate the mammal trees used here, which have been shown to produce robust species prioritisations25,35. For mammals, we used the imputed trees of Gumbs et al.35, which included the 400 missing species and utilised the same imputation approach applied here. This produced 1,000 sets of phylogenies for each of the eight clades, together comprising 70,426 jawed vertebrate species. We used extinction risk data for all species available from the IUCN Red List of Threatened Species version 2021.159, for which there were 44,530 data sufficient assessments (Least Concern [LC], Near Threatened [NT], Vulnerable [VU], Endangered [EN], Critically Endangered [CR] and Extinct in the Wild [EW]) and 3,273 Data Deficient assessments.
Global status of jawed vertebrate evolutionary history
To estimate the global status of jawed vertebrate evolutionary history, we calculated the amount of evolutionary history we currently stand to lose, given the extinction risk of species within each clade50,86. This approach calculates the amount of evolutionary history from each phylogenetic branch we expect to lose, given the extinction probability of all descendant species of that branch. For example, if we consider a branch descended by two species, both of which have a probability of extinction of 0.5, there is a 0.5*0.5 (0.25) probability of losing the branch’s evolutionary history, and thus the expected loss is 25% of the branch’s evolutionary history. When this is calculated and summed for all branches in the phylogenetic tree, we have the total amount of evolutionary history that we stand to lose for the entire clade, and from this we can derive the proportion of total evolutionary history expected to be lost. The proportion of the total evolutionary history of a clade we expect to lose underpins the Phylogenetic Diversity indicator adopted by the CBD’s GBF17,33.
We calculated the amount of evolutionary history at risk of being lost for all 1,000 trees of each clade under three existing extinction risk weightings onto which we mapped the IUCN Red List categories of species: the ‘IUCN50’ and ‘IUCN500’ extinction risk weightings34, which represent conversions of IUCN Red List categories to probability of extinctions of today’s threatened species at 50 years and 500 years into the future, respectively; and the ‘EDGE2’ weighting, which is used to identify priority EDGE species for conservation35 and underpins the PD and EDGE indicators included in the GBF17,33.
We incorporated uncertainty around extinction risk in two ways, following the EDGE2 prioritisation framework35: 1) we generated a distribution of extinction risk values associated with each Red List category derived from a fitted curve, for which the median of each set of values associated with a given Red List category aligned with its original IUCN50, IUCN500, or EDGE2 value; and 2) for each iteration, Data Deficient (DD) and Not Evaluated (NE) species had their extinction risk weighting selected from the entire distribution of extinction risk values at random, generating a highly uncertain distribution of scores with a median equivalent to the Vulnerable Red List category. Thus, during the calculation of threatened evolutionary history on any given phylogenetic tree, DD and NE species have a 60% chance of being assigned an extinction risk weighting associated with one of the three threatened Red List categories (VU, EN, CR) and a 40% chance of being assigned an extinction risk weighting associated with the two non-threatened Red List categories (LC, NT). This proportion is comparable to some predictions that 59% of jawed vertebrates may be at risk of extinction51.
An alternative approach to incorporating DD and NE species would be to draw their extinction risk weighting species in proportion with the observed distribution of Red List categories for a given clade. However, the extent to which extinction risk is dispersed—or clumped—on the Tree of Life varies between and within clades87–89. Using clade-specific parameters to draw the extinction risk values of DD/NE species would lead to significant inter-clade variation in the impact of DD/NE species on the weighting of internal branches for calculations of threatened evolutionary history, further compounded by the different proportion of DD/NE between clades. Further, DD/NE terrestrial vertebrate species are likely to be similarly range-restricted and human-impacted to EN and CR species2, and more than half of DD species are likely threatened51. Thus, using the observed distribution of extinction risk for all jawed vertebrates, whether clade-specific or pooled for the entire group, to parameterise the extinction risk values of DD/NE species would significantly underestimate their extinction risk (e.g. the observed proportion of threatened amphibians is ~40%, whereas 85% of DD species are predicted to be threatened, and the observed proportion of threatened jawed vertebrates is ~19% whereas 59% of DD species are predicted to be threatened)51.
We used terminal branch lengths as a measure of minimum evolutionary distinctiveness (ED) that is insensitive to the choice of extinction risk weighting on internal branches. We calculated the median terminal branch length for each species across the 1,000 phylogenies of their respective clade. To compare terminal branch length between clades, and for all subsequent cross-clade analyses, we used ANOVA with Tukey’s Honest Significant Difference test for pairwise comparison. We used Welch’s t-test to compare the terminal branch lengths of threatened and non-threatened species within each clade.
Priority vertebrates for conservation
We calculated the amount of evolutionary history currently at risk of being lost for each species using the ‘EDGE2’ protocol, which sums the threatened evolutionary history associated with each phylogenetic branch ancestral to a given species and can be derived from the extinction risk-transformed trees generated to calculate the clade-level expected loss of evolutionary history outlined above33,35,90. For each species, we calculated the median EDGE score from the distribution of 1,000 extinction risk-transformed phylogenies. To estimate the unique evolutionary history of each species, we calculated its median terminal branch length from the 1,000 phylogenetic trees. For each species we also calculated the amount of evolutionary history the species is expected to be responsible for in the future, given the extinction risk of all other species in the tree, previously referred to as ‘Heightened ED’90 or ‘ED2’ scores35. To ensure these findings are immediately applicable for conservation action, we used the ‘EDGE2’ extinction risk weightings to derive EDGE scores for all species, as this approach underpins both the EDGE Lists produced for all clades by the Zoological Society of London’s EDGE of Existence programme, and the PD-related indicators included in the Kunming-Montreal GBF 17,33,35.
To identify priority jawed vertebrates whose conservation would capture large amounts of threatened evolutionary history, we ranked all species by their EDGE score in descending order. We arbitrarily selected the top 1% (704 spp.) to represent our priority set of species for further exploration91. To explore the proportional representation of clades in the top 1%, we calculated: 1) the proportion of the top 1% species comprised by each clade; 2) the proportion of each clade’s total richness present in the top 1%; and 3) the proportion of threatened richness of each clade present in the top 1%. We also identified all species that are listed as ‘possibly extinct’ on the IUCN Red List for all clades, and also identified the ‘lost’ taxa identified by Martin et al.37 for tetrapods to estimate the number of priority EDGE species that are possibly extinct or likely to be lost, by calculating the evolutionary history that is likely to have been lost with the extinction of all ‘lost’ or possibly extinct species. For this, we used the conservative estimate of summing the terminal branch lengths for all lost or possibly extinct species. We also applied the EDGE2 criteria to identify sets of priority species for conservation in each clade (EDGE species): above median EDGE for the clade in 95% of calculations and in a threatened IUCN Red List category to guide conservation action (full list of species and their EDGE scores in Supplementary Data 1).
Evolutionarily distinct lineages
Following IUCN’s 2012 resolution to halt the loss of evolutionarily distinct lineages18, there has been, to our knowledge, no subsequent direct assessment of the status of the lineages considered of high conservation importance. Here, as per the IUCN resolution, we focus our assessment on threatened monotypic families and families in which all species are threatened as of utmost conservation importance to avoid large losses of irreplaceable biodiversity.
To calculate a family-level EDGE score, we calculated the mean EDGE score for all species within the family and assigned that to the family. We hereby define ‘EDGE Lineages’ as the families for which: (1) all data sufficient species are threatened; (2) the family-level EDGE score is above median for all families; and (3) at least 50% of species in the family have been assessed in a data sufficient category on the IUCN Red List. We calculated the mean, rather than the median, EDGE score for each family as the mean EDGE score in this case would represent the average amount of threatened evolutionary history we could expect to conserve with conservation action on a random species in the family. In practice, there was minimal difference between the use of median versus mean EDGE scores to characterise EDGE at the family level (>98% similarity in families selected as EDGE lineages).
We also calculated the overall age of each family as the median stem age (i.e., the distance from the tip of the tree to the most recent internal node shared by the family and its sister) across all phylogenetic trees. We took this measurement, rather than the crown age of families, for consistent inclusion of monotypic families, which lack an analogous crown age. For each clade and all jawed vertebrates combined, we then calculated the number and proportion of monotypic families, the number and proportion of fully-threatened families, and the number and proportion of families in which all data sufficient species are threatened.
For monotypic families we calculated the proportion of species listed as having decreasing, increasing, stable or unknown population trends on the IUCN Red List and compared this with a null expectation. For the null expectation we selected random sets of species that matched those in monotypic families (both in clade affinity and Red List category) from the pool of all species with available population trend data. In other words, the sets of species comprising the null distribution are equal to the set of monotypic families in terms of proportion of species belonging to each clade and to each Red List category. We repeated this to generate a distribution of 1,000 sets of trends representing the background expectation for jawed vertebrates. We compared the numbers of species listed under each trend category for monotypic families and all jawed vertebrates using one-sampled t-tests, with the observed number of monotypic species in each category as the ‘true’ value of the mean against which to compare the null distribution.
Using taxonomy to identify candidate EDGE species
Currently, the identification of priority EDGE species is limited to groups for which extensive phylogenetic and IUCN Red List data are available. However, this is currently the case for a small fraction of the entire Tree of Life. We used taxonomic and Red List information with the view of determining criteria for identifying candidate EDGE species for conservation action from clades lacking the phylogenetic data to conduct comprehensive EDGE assessments (i.e. ~99% of the Tree of Life65). To do this we calculated the proportion of species in each clade that: (1) are in monotypic or fully-threatened families; and (2) are Vulnerable, Endangered, Critically Endangered or Extinct in the Wild on the IUCN Red List (i.e. threatened); and also meet the EDGE species criteria (above median EDGE for the clade in 95% of calculations and in a threatened IUCN Red List category) 20,35. This would then tell us whether threatened species from either monotypic families or families where all species are threatened consistently meet the EDGE criteria.