《BMC Ecology and Evolution》：Ecological niche divergence and specialization in Dianthus pseudocrinitus, a neo-endemic species: ecological evidence challenges sister-taxon delimitation
编辑推荐：
这篇研究性论文（非综述）通过整合多维生态位分析，揭示了伊朗东北部石竹属三个近缘类群间的显著生态分化。研究表明，土壤特性是驱动其生态位分离的主要因素，为Dianthus pseudocrinitus的独立物种地位提供了有力证据，并支持了同域分化假说，强调了微生境特化在山区生物多样性形成中的关键作用。
本研究题为“生态位分化与特化：生态证据挑战石竹属新特有类群的姐妹种划分”，聚焦于分类学上存在争议的石竹属新特有类群Dianthus pseudocrinitus及其近缘类群D. crinitus subsp. turkomanicus和D. orientalis subsp. stenocalyx。文章采用了空间显式、整合性的生态学方法，旨在通过探究其生态位动态，为解决该形态和遗传复杂类群的物种边界问题提供证据。
引言
理解生态位如何分化并驱动物种形成与分布，是进化生态学的核心课题，尤其是在环境梯度复杂和历史过程交织的山区。石竹属是一个分类学上公认具有挑战性的类群，具有高度的形态多样性、快速的物种辐射以及频繁出现的细微形态和低分子水平分化。这通常与多倍体、杂交和基因组重复等进化过程有关，模糊了系统发育关系，为传统分类方案带来了困扰。因此，解决石竹属的物种边界问题，需要整合形态、细胞、解剖、分子以及生态信息的综合性分类框架。生态位理论为此提供了一个强大的概念透镜，特别是当遗传或形态分化较新或较微弱时，生态位保守或分化模式可以为谱系独立性、适应策略和进化轨迹提供关键见解。本研究关注的Khorassan–Kopet Dagh植物区是伊朗-图兰地区的主要植物多样性中心之一。该地区拥有半干旱到温带山地气候、显著的海拔梯度和异质性的岩石-钙质基质，形成了支持高水平特有现象的精细环境镶嵌。在此背景下，本研究旨在：评估D. pseudocrinitus与其近缘类群的生态位相似性，以阐明其分类学亲缘关系和物种边界；评估这些类群间的生态位分化程度，以检验物种水平的独特性；识别驱动其生态位分化的关键环境因子；并考察生态位保守与生态转移的信号，以阐明近期分化和山区物种形成的生态机制。
材料与方法
研究数据收集了三个类群在伊朗东北部及土库曼斯坦南部的自然分布范围内的发生记录，并通过文献、标本馆、GBIF及SpeciesLink等来源进行补充。为减少空间聚集，使用R包spThin进行了空间稀释处理。最终数据集分别包含66、61和20个独立的发生记录用于模型构建和分析。
为建立生态学上稳健的物种分布模型，研究综合了四类环境预测变量：气候、植被生产力、地形和土壤性质。所有变量均重采样至30弧秒（约1公里）分辨率。为解决多重共线性并降低维度，对每一类环境变量分别进行了主成分分析，保留解释大部分方差的前三个主成分，最终得到12个正交的预测变量用于后续分析。
考虑到校准区域选择的重要性，本研究使用R包grinnell，基于各物种的发生数据和所有12个环境预测变量，通过扩散模拟为每个类群定义了特定的可达区域，作为模型校准的范围。
生态位建模使用基于最大熵原理的Maxent算法，通过R包kuenm实现。建模过程分为两阶段，针对样本量较小的D. pseudocrinitus进行了专门的参数优化。模型评估采用了折刀交叉验证，并根据部分受试者工作特征曲线统计显著性和AICc（赤池信息量准则）等指标筛选最优模型。最终模型通过10次自助法重采样构建，并基于最低训练存在阈值将连续预测转换为二元适宜性分布图。
为全面评估生态位特征和环境关系，研究采用了四种互补的分析方法。首先，对12个PCA衍生变量进行标准化后，进行多变量方差分析和线性判别分析，以检验和可视化类群间的整体环境差异及关键区分变量。其次，对每个环境变量进行单变量密度重叠分析，量化其贡献。最后，基于显著环境变量，使用hypervolume R包构建每个类群的高维生态位超体积，以量化生态位体积、交集、并集以及物种特异性独特部分，并计算Jaccard和S?rensen相似性指数来衡量生态位相似性与分化程度。
结果
可达区域模拟结果显示，D. crinitus subsp. turkomanicus的可达区覆盖了Khorassan–Kopet Dagh的大部分地区，而D. orientalis subsp. stenocalyx和D. pseudocrinitus的可达区则主要位于该区域的北部。
生态位模型显示，环境变量对模型的贡献在不同类群间差异显著。在D. crinitus subsp. turkomanicus中，地形是主要预测因子；在D. orientalis subsp. stenocalyx中，土壤变量贡献最大；而D. pseudocrinitus的生态位则受土壤和地形的共同影响。最终的最优模型参数各不相同，但均表现出较好的预测性能。高栖息地适宜性区域集中分布于山区。D. crinitus subsp. turkomanicus的高适宜区范围最广，覆盖多个山系；D. orientalis subsp. stenocalyx和D. pseudocrinitus的高适宜区则更为集中。模型不确定性在核心高适宜区较低，在边缘和生态过渡带较高。基于二元分布图的重叠分析表明，D. crinitus subsp. turkomanicus和D. orientalis subsp. stenocalyx共享的适宜栖息地面积最大，而与D. pseudocrinitus相关的重叠则更加破碎和空间受限，表明其栖息地偏好更为特化。
多变量方差分析揭示了一个高度显著的物种整体效应对组合环境变量的影响，证实了三个石竹类群间存在不同的环境生态位。土壤变量是所有差异中最显著的驱动因素，三个土壤主成分均表现出高度显著性。NDVI和部分生物气候变量也显示出显著分化。线性判别分析进一步支持了这些发现。LDA1解释了99.13%的方差，其载荷主要由土壤变量主导，成功将D. crinitus subsp. turkomanicus与其他两个类群区分开来。LDA2则主要区分了D. pseudocrinitus。密度重叠分析显示，部分生物气候变量重叠度较高，表明存在共享的气候条件；而NDVI和所有土壤主成分则存在显著分化，其中pc1_soil的重叠度最低，表明土壤特性是生态位分化的主要轴线。
生态位超体积分析量化了这些模式。生态位体积差异巨大：D. orientalis subsp. stenocalyx拥有最大的生态位超体积，其次是D. crinitus subsp. turkomanicus，而D. pseudocrinitus的生态位体积则要小得多，表明其生态广度非常狭窄。最大的成对生态位重叠发生在D. crinitus subsp. turkomanicus和D. orientalis subsp. stenocalyx之间。相比之下，D. pseudocrinitus与另外两个类群的重叠都非常小，Jaccard指数均≤0.162。独特的生态位部分呈现出强烈的非对称性，D. crinitus subsp. turkomanicus和D. orientalis subsp. stenocalyx保留了相对于D. pseudocrinitus大得多的独特生态位空间。体积比也凸显了D. pseudocrinitus极为受限的生态位宽度。
讨论
本研究通过整合生态位建模、多元统计分析和超体积量化等多种方法，一致且显著地证明了D. pseudocrinitus与其形态和遗传上相似的近缘类群之间存在生态位分化。空间分析精确地描绘了三个类群在伊朗东北部的潜在分布，并揭示了栖息地适宜性在环境梯度上的显著变异。高适宜性区域与已知分布区基本吻合，但模型在边缘地区存在一定不确定性。栖息地适宜性预测，特别是基于D. pseudocrinitus生态位校准的预测，揭示了类群间潜在的生态重叠与分化。紫色区域所指示的空间重叠暗示了生态位保守性的存在，这可能反映了共享的环境偏好，或指示了潜在的生态接触区，其中可能发生种间相互作用。这些发现凸显了将空间生态学与进化框架整合，以更全面理解山区环境中塑造生态多样化的复杂过程的重要性。
多变量方差分析和线性判别分析的结果共同表明，土壤特性是驱动生态位分离的主要环境因素，而生物气候和NDVI变量则在塑造更精细尺度的差异中起次要作用。密度重叠分析进一步阐明了生态位划分的模式：尽管某些生物气候变量重叠度较高，但NDVI和所有土壤主成分均存在显著分化。这为D. pseudocrinitus表现出功能分化的观点提供了支持，尤其是在受干扰的退化生境中。观察到的生物气候高度重叠与土壤和植被相关因子的强烈分化之间的二分法，为理解这些山区生境中近期分化和物种形成的生态机制提供了关键见解。这种模式对于理解物种如何在同域中共存至关重要：在资源利用或非生物耐受性（如土壤pH或养分有效性）方面的细微差异可以减少种间竞争，并允许近缘物种稳定共存。
生态位超体积的定量分析为三个石竹类群间的生态位分化提供了强有力的支持。D. pseudocrinitus表现出极为狭窄的生态位体积，与其近缘类群的重叠度极低。这种显著的生态隔离和狭窄的生态位，连同其与地理上相邻近缘类群的比较，为同域物种形成提供了令人信服的证据。该物种占据了高度受限的、在环境空间中明显可区分的生态位，这为其在存在形态相似性和地理邻近性的情况下，仍具有独立的物种地位提供了强有力的生态学指标。这一强烈的、精细尺度的生态分化模式，连同一些更广泛的生物气候生态位保守性，突显了塑造物种边界的因素的动态相互作用。研究结果强调了微生境特化和局部适应在驱动复杂山地环境多样化方面的潜在重要性，并倡导采用整合性的分类学方法来解决复杂的物种划分挑战。
