基因组学揭示隐蔽冠层特化斐济 Hylaeus 蜂种群崩溃与分布区收缩：H. navai 有效种群大小约99%下降及极危评估证据

《Ecology and Evolution》：Genomics Reveal Population Crash and Range Contraction in Elusive Fijian Hylaeus Bees

　　

传粉者的保护因其在生态系统功能中的关键作用而尤为受关注，但系统性的地理和方法学偏差阻碍了基线数据通常难以获得。岛屿作为理解进化过程的简化系统很有用，已知其单位面积蜜蜂物种丰富度高于大陆地区，并且受威胁物种比例不成比例地高。研究人员采用低层植被和冠层采样，结合线粒体和基因组数据，评估了新近发现的冠层专化斐济 Hylaeus 蜂分支。研究人员发现岛间遗传结构——或初期物种形成——需要进一步系统学研究。在数量最多的物种 H. navai 中，研究人员利用全基因组单核苷酸多态性（single nucleotide polymorphisms, SNPs）结合种群历史重建，推断自欧洲殖民以来有效种群大小（effective population size, N_e）下降约99%（最保守至最不保守分别为81.6%–99.9%）。在 H. navai 种群历史中三个独立的10年期内，研究人员推断种群大小下降分别为61%、76%和82%。这些下降为国际自然保护联盟（International Union for Conservation of Nature, IUCN）评估提供了证据，这种幅度的种群下降符合极危（Critically Endangered）列级标准。研究人员发现观测杂合度低和近交升高，与遗传瓶颈（genetic bottleneck）一致。研究人员使用物种分布模型（species distribution models, SDMs）显示 H. navai 的适宜栖息地占斐济陆地面积的1.4%–2.4%，且现在可能局限于高海拔森林。通过检查历史采集记录，研究人员提出对当代从低地地区迁出/消失的担忧。研究人员认为，栖息地改变（由于人类人口增加）、入侵性非洲郁金香树（African tulip tree，潜在生态陷阱）、气旋，以及来自本地和引入传粉者的竞争加剧，可能共同相互作用驱动下降。研究人员的研究结果可能不限于斐济蜜蜂，并提供了隐藏冠层传粉者下降的实证证据，这应成为太平洋及全球评估的优先事项。

论文主体内容解读

研究背景与问题：全球昆虫下降威胁营养循环、土壤形成、分解、生物防治和传粉等生态系统功能，但昆虫保护代表性不足，且存在地理和方法学偏差。传粉者尤为关键，但其下降常在研究不足地区和隐蔽类群中未被发现。蜜蜂是多数生态系统中主要动物传粉者，但IUCN评估比例低。岛屿具有高特有性和单位面积较高蜜蜂丰富度，却是蜜蜂多样性和种群动态研究最少的系统之一；传统采样还系统性漏掉森林冠层传粉者。斐济正成为利用分子技术和特有蜜蜂研究生态、进化与保护问题的模型系统，本地地巢 Lasioglossum（Halictidae）研究较多，而六种冠层专化 Hylaeus（Colletidae）直到冠层采样常规实施后才被发现。Dorey et al. 描述这些体型小、携带花粉、可能特化红色花冠层植物的斐济 Hylaeus，并推测栖息地改变、欧洲殖民后改变、气旋和入侵非洲郁金香树可能影响种群。由于过去监测不可能、现在监测不实际，全基因组溯祖有效种群大小估计成为评估种群历史的重要途径。研究人员整合系统随机冠层采样、机会性花访调查、种群基因组学和物种分布模型，试图澄清物种边界、评估种群遗传结构和连通性、重建种群历史、描述遗传多样性和种群健康指标、建模生态位和预测当代栖息地适宜性。假设为：斐济 Hylaeus 种群结构反映残余森林斑块结构且连通性低；有效种群大小因栖息地丧失和破碎化向现在下降。

主要关键技术方法：研究人员在斐济 Viti Levu 按海拔和土地利用分层，在公共道路网络随机设置主样点和辅助样点，并补充 Mt. Tomanivi、Mt. Batilamu 和 Monasavu Dam；2025年4月14日至5月11日以低层植被和冠层扫网采样。共采集101份 Hylaeus；94份用于 Diversity Arrays Technology（DArT）基因组SNP分析（85 H. navai、5 H. veli、3 H. derectus、1 H. albaeus），16份用于线粒体细胞色素氧化酶 I（cytochrome oxidase I, COI）测序。历史/比较样本包括 Dorey et al. 先前的 Taveuni H. apertus 和 H. navai 记录等。分析包括COI贝叶斯系统发育、SNP过滤、主成分分析（principal component analysis, PCA）/判别主成分分析（discriminant analysis of principal components, DAPC）、tess3r、Mantel检验、Tajima's D、NeEstimator2、Stairway Plot 2和最大熵（maximum entropy, MaxEnt）物种分布模型，并与非洲郁金香树模型比较。

研究结果：

3.1 采样结果：研究人员在 Viti Levu 共采样54个地点，包括29个主样点、5个辅助样点和20个补充样点；34个随机生成样点中20个为森林、14个为非森林。共采集101份 Hylaeus，海拔范围371–1334 m，主要花访植物为 Geissois superba 和 Metrosideros collina。26/54个样点低于600 m，但低地仅采集到1份 H. derectus，与历史低地记录形成反差。H. navai 2025年采集91份，海拔883–1334 m；H. veli 6份，海拔1334 m；H. derectus 3份，海拔371–900 m；H. albaeus 1份，海拔923 m。

3.2 遗传结构与系统发育：全基因组SNP的PCA显示 Viti Levu 物种间清晰分异，PC1解释46%、PC2解释18%。COI贝叶斯系统发育拓扑分辨良好，主要节点后验概率支持强；首次获得 Viti Levu 的 H. veli 序列和 H. navai 雄性标本。Viti Levu 与 Taveuni 的 H. navai 和 H. apertus 分支长度几乎等于两岛 H. veli 的分化，H. albaeus 两标本间节点更深。研究人员认为存在岛间遗传结构或初期物种形成，需要进一步系统学研究。

3.3 H. navai 的种群遗传学：女性-only SNP分析显示种群结构证据少，PCA PC1和PC2分别解释3.3%和3.1%，DAPC最可能为一个种群，tess3r交叉验证未出现平台。观测杂合度 H_o=0.08显著低于期望杂合度 H_e=0.10，近交系数 F_IS=0.18，Mantel检验无显著隔离-by-distance，Tajima's D=?1.55但不显著。Stairway Plot显示有效种群大小长期变化，约1550年开始下降，1874年后多数年份下降；1933–1942下降约61%，2002–2011下降约76%，2016–2025下降约82%；自欧洲殖民约150年下降约99%，最保守至最不保守范围为81.6%–99.9%。当代N_e估计因等位基因频率阈值不同而不一致，MAF为0.05时N_e=121，研究人员认为该值最可靠但精度有损失。

3.4 物种分布建模：H. navai-apertus模型预测性能高，AUC=0.99，适宜栖息地占斐济2.4%，主要预测因子为最热月最高温、土地利用和海拔；H. navai-only模型AUC=0.99，适宜栖息地占1.4%，主要预测因子为等温性和海拔，集中于中高海拔。非洲郁金香树模型AUC=0.88，适宜栖息地占13.4%，集中于低海拔，与 Hylaeus 模型 largely 反向。

讨论总结：

4.1 采样与海拔模式：研究人员采用系统冠层和低层采样，覆盖 Viti Levu 全海拔，并广泛采样低地；但低地 Hylaeus 几乎缺失，而历史上低地有记录，这提示低地种群可能显著下降或局部灭绝。H. navai 和 H. veli 海拔范围窄，H. derectus 和 H. albaeus 较宽。

4.2 物种边界与分化：COI系统发育与 Dorey et al. 拓扑一致。Viti Levu 与 Taveuni 的 H. navai 和 H. apertus 分化类似 H. veli 岛间分化，可能为初期物种形成或物种概念过度拆分；研究人员建议将其视为进化显著单元。H. albaeus 深分歧可能因模式标本序列质量差或该标本为 H. breviflavus。未来需多证据和博物馆标本，≥50个体可稳健重建种群历史，2个SNP样本可计算 pairwise F_ST。

4.3 H. navai 的种群结构与过去种群历史：SNP分析未发现高地采集点间种群结构，与历史连通一致；但因仅采自连续高地森林，未能有意义检验第一假设。COI显示岛间分化，可能存在更大尺度结构。当代N_e估计不一致，研究人员未解释；Stairway Plot下降与斐济生态和人为变化大致吻合，最早约1556年开始下降，后欧洲定居后持续下降，并出现三个主要10年下降。

4.4 物种分布建模与栖息地适宜性：SDM预测性能高，H. navai适宜栖息地占1.4%–2.4%，集中中高海拔；非洲郁金香树适宜低海拔且面积更大。模型未独立验证，应视为假设，但与野外观察一致。非洲郁金香树随气候变暖可能扩张，构成潜在生态陷阱威胁。

4.5 保护意义：H. navai尽管最常见，仍下降99%，三个10年下降分别可支持IUCN Vulnerable、Endangered和Critically Endangered。下降未见减弱。可能驱动包括人类人口增长、栖息地破坏和破碎化、非洲郁金香树、气旋、与 Lasioglossum (Homalictus) fijiense 和欧洲蜜蜂的竞争，以及热过滤可能。研究人员建议将 H. navai 列为极危（CR A2b）。

结论翻译：研究人员使用跨学科方法和分析发现冠层专化传粉者下降的稳健证据。该研究整合系统野外调查、种群基因组学和物种分布模型，评估特有斐济 Hylaeus 蜂的保护状况和种群历史。配对低层和冠层采样支持斐济 Hylaeus 是冠层和红花专化者。然而，物种边界不确定和岛间遗传分化突出需要扩大基因组分析和系统采样。约150年间 H. navai 有效种群大小约99%下降，与欧洲殖民、栖息地改变和入侵物种建立同时发生。SDM及历史和当代采样推断表明，H. navai 及更广泛的斐济 Hylaeus 可能从低地收缩，低地人为改变和入侵物种影响可能最高。模型表明 H. navai 局限于斐济陆地面积的1.4%–2.4%，大部分在更高海拔。采样、分子和地理空间分析提供重大分布区和种群收缩的强证据。研究人员建议这构成基于IUCN标准A2b的 H. navai 极危（CR A2b）状态基础。实证识别和抵消主要下降驱动因素对所有斐济 Hylaeus 蜂具有立即保护重要性。研究人员建议斐济 Hylaeus 可作为理解隐蔽类群过去下降的模型系统，并强调系统学在隐蔽类群灭绝前识别它们的基础重要性。由于冠层专化 Hylaeus 预计广泛分布于太平洋，连接澳大利亚和法属波利尼西亚类群，研究人员建议对隐蔽 Hylaeus 物种的威胁将广泛存在。研究人员还强调全球可能存在隐蔽冠层专化传粉者。保护生物学家应考虑使用系统和靶向多分层采样方法与分子技术相结合，以对种群动态和下降作出有力推断，这是仅用标准生态采样不可能实现的。由于大多数类群不存在基线数据，特别是了解不足的无脊椎动物，该整合框架为全球检测森林冠层隐蔽下降和驱动因素提供路线图。