Data from: Ranked tree shapes, non-random extinctions and the loss of phylogenetic diversity

Phylogenetic diversity (PD) is a measure of the evolutionary legacy of a group of species, which can be used to define conservation priorities. It has been shown that an important loss of species diversity can sometimes lead to a much less important loss of PD, depending on the topology of the species tree and on the distribution of its branch lengths. However, the rate of decrease of PD strongly depends on the relative depths of the nodes in the tree and on the order in which species become extinct. We introduce a new, sampling-consistent, three-parameter model generating random trees with covarying topology, clades relative depths and clades relative extinction risks. This model can be seen as an extension to Aldous' one parameter splitting model (β, which controls for tree balance) with two additional parameters: a new parameter α quantifying the relation between age and richness of subclades, and a parameter η quantifying the relation between relative abundance and richness of subclades, taken herein as a proxy for overall extinction risk. We show on simulated phylogenies that loss of PD depends on the combined effect of all three parameters, β, α and η. In particular, PD may decrease as fast as species diversity when high extinction risks are clustered within small, old clades, corresponding to a parameter range that we term the `danger zone' (β<-1 or α<0; η>1). Besides, when high extinction risks are clustered within large clades, the loss of PD can be higher in trees that are more balanced (β>0), in contrast to the predictions of earlier studies based on simpler models. We propose a Monte-Carlo algorithm, tested on simulated data, to infer all three parameters. Applying it to a real dataset comprising 120 bird clades (class Aves) with known range sizes, we show that parameter estimates precisely fall close to a danger zone: the combination of their ranking tree shape and non-random extinctions risks makes them prone to a sudden collapse of PD.系统发育多样性（Phylogenetic Diversity, PD）是衡量一类物种演化遗留特征的指标，可用于划定保护优先级。已有研究表明，物种多样性的显著丧失有时仅会导致系统发育多样性的小幅下降，这一结果取决于物种树的拓扑结构及其分支长度的分布。然而，系统发育多样性的下降速率强烈依赖于物种树中节点的相对深度，以及物种灭绝的先后顺序。 本文提出一种全新的、采样一致（sampling-consistent）的三参数模型，可生成兼具协同变化拓扑结构、支系相对深度与支系相对灭绝风险的随机树。该模型可视为对奥尔达斯（Aldous）单参数分裂模型的扩展：原模型仅通过参数β控制树的平衡度，本模型新增两个参数——参数α用于量化亚支系的年龄与其丰富度之间的关联，参数η用于量化亚支系相对多度与丰富度之间的关联（本文中将其作为整体灭绝风险的替代指标）。通过模拟系统发育树的实验表明，系统发育多样性的丧失取决于β、α、η三个参数的共同作用。 具体而言，当高灭绝风险聚集于小型古老支系时，系统发育多样性的下降速率可与物种多样性的下降速率相当，对应的参数范围我们称之为“危险区”（β<-1 或 α<0；η>1）。此外，当高灭绝风险聚集于大型支系时，树结构越平衡（β>0），系统发育多样性的损失反而越高，这与早期基于简化模型得出的研究结论相悖。我们提出一种经模拟数据验证的蒙特卡洛（Monte-Carlo）算法，用于推断上述三个参数。将该算法应用于包含120个已知分布范围的鸟类支系（鸟纲Aves）的真实数据集后发现，参数估计值恰好落在危险区附近：鸟类支系的分级树结构与非随机的灭绝风险组合，使得其极易出现系统发育多样性的骤降。