利用专性互利共生代理解决稀有物种建模悖论——以 Philoliche longirostris 与 Roscoea purpurea 为例

《Conservation Science and Practice》：Leveraging obligate mutualism to resolve the rare species modeling paradox

　　

过去二十年间，物种分布模型（species distribution models, SDMs）已成为制定生物多样性保护规划（biodiversity conservation planning）的重要工具，但建模稀有、数据缺乏物种（rare, data-deficient species）仍是一项重大挑战。为应对这一“稀有物种建模悖论”（rare species modelling paradox），研究人员提出一种新方法：当两个专性互利共生伙伴（obligate mutualistic partners）共享基础生态位（fundamental ecological niche）时，一个记录充分的互利共生伙伴（well-documented mutualistic partner）的分布数据可作为可靠代理（proxy），用于建模其稀有对应物种的分布。研究人员利用一种蝇（Philoliche longirostris）与其食源植物（Roscoea purpurea）之间的专性互利共生（obligatory mutualism）检验了该假设。研究人员首先使用 Schoener's D 和 Hellinger's I 作为代理指标，检验蝇与植物之间的生态位重叠（niche overlap）。随后，利用 R. purpurea 的出现数据，研究人员实施集合物种分布模型（ensemble SDM），为蝇构建代理分布模型（proxy distribution model），并通过在代理模型预测的40个地点对蝇出现情况进行实地验证（ground truthing），检验代理模型的准确性。最后，研究人员将通过实地验证生成的点位与当前已知分布点整合，开发 P. longirostris 的最终分布模型，并对两个模型的出现概率（probabilities of occurrence）进行 Pearson 相关分析（Pearson correlation analysis）。分析结果证实两个物种之间存在显著生态位重叠（Schoener's D = 0.552；Hellinger's I = 0.737）。实地验证证实 P. longirostris 存在于37个地点（准确率92.5%），表明代理模型具有较高预测性能。最终模型与初始代理模型表现出强相关（r = 0.836）和高空间重叠（spatial overlap）。该研究表明，互利共生代理（mutualistic proxies）可有效规避稀有物种的数据限制问题。研究人员提出的框架为增强处于紧密生态相互作用（tight ecological interactions）中的稀有/数据缺乏物种的保护规划提供了一种具有生物学依据的方法，并对其他生态相互作用具有广泛应用前景。

研究背景与问题：全球生物多样性持续下降，栖息地破坏与破碎化、土地利用集约化、生物入侵、杀虫剂和农药使用增加以及全球气候变化等共同威胁物种存续。物种分布模型（species distribution models, SDMs）已成为表征生态位（ecological niche）和预测地理栖息地适宜性的关键工具，但稀有/数据缺乏物种（rare/data-deficient species）因出现点少、地点不确定性高，模型易过拟合，预测精度随样本量下降而降低。传统种群调查虽可靠但昂贵耗时。已有群落水平模型（CLM）、联合SDMs、多个双变量模型集合及伪缺失生成等方法，但仍存在方法复杂、计算量大或未利用生物互作等问题。Philoliche longirostris 为喜马拉雅特有种，曾被认为灭绝，后重新发现为 Roscoea purpurea 的专性传粉者，目前仅约11个准确出现点，分布、幼虫生物学、扩散距离和生态需求知之甚少。R. purpurea 分布较清楚，且二者呈现专性互利共生，植物可无蝇存活，蝇却依赖植物作为食物，因此研究人员提出可用 R. purpurea 分布代理建模 P. longirostris 分布，并开展生态位重叠检验、代理建模、实地验证和最终模型比较。

关键技术与方法（≤250字）：研究人员从文献、GBIF、iNaturalist、博物馆及 KATH、TUCH 等来源获得 R. purpurea 和 P. longirostris 出现点，并空间稀释至 1-km² 分辨率；使用27个环境变量，经方差膨胀因子（VIF）筛选后保留11个变量。采用“ecospat”进行主成分分析（PCA）和核密度估计，计算 Schoener's D、Hellinger's I，并实施生态位等价性与相似性检验。用“sdm”实施八种算法，剔除 GLM 和 FDA 后，以 Boosted Regression Tree（BRT）、Multiple Adaptive Regression Splines（MARS）、Maximum Entropy（MAXENT）、Random Forest（RF）、Recursive Partitioning and Regression Trees（RPART）和 Support Vector Machine（SVM）按 True Skill Statistics（TSS）加权构建集合模型。从 R. purpurea 适宜区随机抽取40个点，于2014–2024年高峰期实地验证，计算混淆矩阵指标与 Partial Relative Bias（PRB）。最终纳入验证点建模，并用 Pearson 相关和空间叠加比较代理与最终模型。

研究结果：
3.1 R. purpurea 与 P. longirostris 的生态位重叠：通过核密度估计与生态位检验得出，二者 Schoener's D 为0.552，Hellinger's I 为0.737；生态位等价性检验 p 值分别为0.931和0.921，不能拒绝两个生态位统计等价的零假设；生态位相似性检验 p 值为0.018和0.016，拒绝随机重叠零假设；观察值位于模拟值95%范围内。结论是二者显著共享相似栖息地需求和气候空间，R. purpurea 出现记录可作为 P. longirostris 的建模代理。
3.2 R. purpurea 分布模型：通过六算法集合模型及 AUC、TSS、Continuous Boyce Index（CBI）评价，模型表现高相关。温度季节性（Bio4）、最干季度平均温度（Bio9）、蒸散和蒸气压为主要贡献变量；R. purpurea 最适最干季度平均温度为0–10°C，Bio4约为600。二值化集合模型显示，其潜在分布西部和中部与当前已知范围接近，东部潜在适宜区更广；西藏虽有少量记录但未显示适宜栖息地；适宜区主要位于尼泊尔中丘至高山、不丹相似生态区以及印度东北和西北。
3.3 地面验证与模型评价：在40个验证点中，37个点发现 P. longirostris，3个点未发现；混淆矩阵为36个真阳性、3个假阳性、1个假阴性、0个真阴性；敏感性0.973、总体准确率0.900、精确率0.923、F1分数0.947，特异性为0，TSS为?0.027，PRB为?0.028。结论是代理模型对蝇实际存在具有较强预测性能，但真缺失样本流行率低。
3.4 P. longirostris 最终模型：仅用11个已知点建立的模型高度碎片化且不现实；纳入实地验证点后的最终模型在六算法中表现高预测性能，Bio4 贡献最大。二值化集合预测在最终模型出现点范围内与代理模型几乎一致，潜在分布从印度 Himachal 西部经尼泊尔、不丹至印度 Arunachal 东部，尼泊尔中丘和高山、相似生态区以及印度东北和西北为适宜区。代理模型与最终模型出现概率呈显著强相关（Pearson r = 0.8357，p < 0.05），空间分布几乎完全重叠。结论是代理模型预测性能与最终模型相当，使用 R. purpurea 分布点建模 P. longirostris 合理。

讨论部分总结：研究人员认为专性互利共生可克服稀有物种建模悖论，无需为数据缺乏伙伴追加大量出现数据。其前提是互利伙伴生态位显著重叠；R. purpurea 地理生态位较宽，而 P. longirostris 生态位较窄并嵌套其中，样本量分别为 n = 161 和 n = 47，植物克隆繁殖可形成更宽实际生态位，蝇的实际生态位则受植物可用性限制。实地验证中3/40点无蝇，证明植物生态位更宽；92.5%成功率说明嵌套并未导致不可接受假阳性。由于采用存在–背景方法，这3个缺失点未作为真缺失。代理模型与最终模型 AUC 均超过0.9，Pearson r = 0.8357，敏感性、准确率、精确率、F1分数高，PRB 接近0，空间重叠高，表明代理模型与最终模型性能相当。方法上，少出现点加多预测因子不能产生可靠预测，已有创新方法仍有局限。该框架适用于真实稀有或采样不足的数据缺乏物种，前提是显著生态位重叠；可扩展至其他紧密生态相互作用，但建议满足 Schoener's D 和 Hellinger's I > 0.5、相似性检验 p < 0.05、等价性检验 p > 0.05，且代理种空间稀释后至少50个出现点。中度重叠可谨慎应用并独立验证。局限性包括地面验证限于尼泊尔喜马拉雅且主要在预测适宜区内，未设预测范围外真负对照，主要评估敏感性；未来建议纳入30–50个负对照点。保护意义上，该互利系统为喜马拉雅唯一已知长舌蝇传粉系统，研究新增37个 P. longirostris 地点，建议在印度 Himachal、尼泊尔中丘和高喜马拉雅、印度 Assam、不丹及 Arunachal 至 Brahmaputra 河边界进一步调查，为稀有物种及其栖息地保护提供可扩展框架。

研究结论翻译：本研究以喜马拉雅一个知之甚少的传粉者及其专性食物来源为案例，在物种分布模型（SDM）领域引入了一种新方法。通过整合互利共生伙伴共享相似宏观栖息地需求的生态学理论与先进建模技术，本研究为稀有/数据缺乏物种建模悖论提供了有说服力的解决方案。互利共生代理的使用不仅提高了分布预测的准确性，也强调了保护生态相互作用的重要性。研究人员得出结论，其发现为保护稀有物种及其栖息地提供了可扩展框架，并在全球气候变化背景下有助于更广泛的生物多样性保护目标。