Institut de Biologie de l'École Normale Supérieure (IBENS), École Normale Supérieure, CNRS, INSERM, PSL Research University, Paris, France.
ED 227, Sorbonne Universités, Paris, France.
Syst Biol. 2018 Nov 1;67(6):1025-1040. doi: 10.1093/sysbio/syy030.
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 ($\beta$, which controls for tree balance) with two additional parameters: a new parameter $\alpha$ quantifying the relation between age and richness of subclades, and a parameter $\eta$ 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, $\beta$, $\alpha,$ and $\eta$. 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" ($\beta<-1$ or $\alpha<0$; $\eta>1$). Besides, when high extinction risks are clustered within large clades, the loss of PD can be higher in trees that are more balanced ($\beta>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 data set comprising 120 bird clades (class Aves) with known range sizes, we show that parameter estimates precisely fall close to the danger zone: the combination of their ranking tree shape and nonrandom extinctions risks makes them prone to a sudden collapse of PD.
系统发育多样性 (PD) 是衡量一组物种进化遗产的指标,可用于确定保护优先级。研究表明,物种多样性的重要损失有时会导致 PD 的损失小得多,这取决于物种树的拓扑结构及其分支长度的分布。然而,PD 的减少速度强烈取决于树中节点的相对深度以及物种灭绝的顺序。我们引入了一种新的、采样一致的、三参数模型,该模型生成具有共变拓扑、分支相对深度和分支相对灭绝风险的随机树。该模型可以看作是 Aldous 的一个参数分裂模型($\beta$,控制树平衡)的扩展,增加了两个额外的参数:一个新参数 $\alpha$,用于量化子分支年龄和丰富度之间的关系,以及一个参数 $\eta$,用于量化子分支相对丰度和丰富度之间的关系,在此用作总体灭绝风险的代理。我们在模拟系统发育树上表明,PD 的损失取决于所有三个参数,即 $\beta$、$\alpha$和$\eta$的综合效应。特别是,当高灭绝风险集中在小的、古老的分支中时,PD 可能会像物种多样性一样迅速下降,这对应于我们称为“危险区”的参数范围($\beta<-1$或$\alpha<0$;$\eta>1$)。此外,当高灭绝风险集中在大分支中时,在更平衡的树中($\beta>0$),PD 的损失可能更高,这与基于更简单模型的早期研究的预测相反。我们提出了一种蒙特卡罗算法,在模拟数据上进行了测试,以推断所有三个参数。将其应用于包含已知范围大小的 120 个鸟类分支(类鸟类)的真实数据集,我们表明参数估计值恰好接近危险区:它们的排序树形状和非随机灭绝风险的组合使它们容易出现 PD 的突然崩溃。