《Ecology and Evolution》：Roads, Soil, Snow, and Topography Influence Genetic Connectivity: A Machine Learning Approach for a Peripheral American Badger Population
编辑推荐：
本文针对加拿大不列颠哥伦比亚省濒危的美洲獾西部边缘种群（Taxidea taxus jeffersonii），面临栖息地破碎化与道路致死威胁，基因流受阻的问题，研究人员运用梯度提升机（GBM）与基于廊道的景观遗传学模型，量化了景观特征对遗传距离的影响。研究发现，碎石质土壤、陡峭坡度、主要道路和大量积雪会阻碍基因流动，而有机质与冲积土壤则促进连通性，空间遗传自相关可达约110 km。该研究为通过栖息地保护、道路地下通道建设等针对性管理措施，以促进该濒危种群的基因流动和长期存活提供了科学依据。
在加拿大不列颠哥伦比亚省（British Columbia, BC）南部的山谷与山峦之间，生活着一群特别的居民——美洲獾的一个亚种，杰弗逊獾（Taxidea taxus jeffersonii）。然而，它们的生存现状岌岌可危。作为地理上的边缘种群，它们不仅与中心种群隔离，还面临着栖息地丧失、道路致死以及气候变化等多重威胁。其中，公路是导致其死亡的最大直接原因，而适宜它们挖掘洞穴的疏松土壤（如冲积土、湖积土）又偏偏与人类聚居的谷底高度重叠。这使得本就数量稀少（西部种群估计仅150-245只）的獾群，其个体间的基因交流（即基因流动）变得异常困难。长此以往，种群会因近亲繁殖导致遗传多样性下降，适应环境变化的能力减弱，灭绝风险大大增加。那么，究竟是哪些具体的自然与人为景观因素，在阻碍着这些獾的扩散与基因交流？又该如何精准地绘制出“基因通行图”，来指导保护行动呢？为了回答这些问题，由Eric C. Palm等人领导的研究团队展开了一项深入的景观遗传学研究，其成果发表在了期刊《Ecology and Evolution》上。
为了阐明景观特征对杰弗逊獾遗传连通性的影响，研究人员运用了一套结合了现代遗传学与机器学习技术的分析方法。首先，他们收集了2008年至2023年间在BC南部四个生态区域（Cariboo, Thompson, Nicola, Okanagan）获得的116份獾的遗传样本（主要来自路杀个体或调查），并利用14个微卫星（Microsatellite）位点进行了基因分型。基于这些基因数据，他们计算了所有个体对之间的三种遗传距离（欧氏遗传距离、共享等位基因比例、Queller and Goodnight亲缘关系估计值），以量化遗传差异。在景观分析方面，研究采用了基于廊道（或称为样带）的建模方法。他们在每对样本点之间创建了直线，并将其缓冲为1000米宽的样带，提取了13个假设与獾扩散相关的景观变量（如土壤母质类型、坡度、主要道路、年降雪量、冠层覆盖度等）在样带内的平均值。随后，他们使用梯度提升机（Gradient Boosting Machine, GBM）模型，以遗传距离为响应变量，以地理距离和景观变量为解释变量进行建模。该过程还包含了空间交叉验证和变量选择，以防止过拟合并确保模型的泛化能力。通过迭代计算基于模型预测阻力面的最小成本路径（Least-Cost Path, LCP）并重新建模，他们最终确定了最能解释遗传变异的景观特征，并据此生成了景观阻力面。最后，利用阻力核（Resistant Kernel）分析方法，预测了研究区域内高遗传连通性的潜在廊道区域。
3.1 遗传汇总统计与空间遗传结构
对116份样本的遗传分析显示，所有位点观测杂合度为0.609至0.853。空间遗传自相关分析表明，遗传相似性在约110公里距离内仍然存在，研究者将此距离定义为遗传邻域半径。稀薄等位基因丰富度分析显示，遗传多样性从西北部的Cariboo地区向东南部的Okanagan地区递增，这与中心-边缘假说一致。STRUCTURE聚类分析将种群主要分为西北部（Cariboo）和东南部（其他三区）两个遗传集群。曼特尔检验（Mantel test）证实了隔离距离的存在。
3.2 景观阻力模型结果
在诸多候选模型中，基于直线样带、以欧氏遗传距离为响应变量的GBM模型表现最佳。该模型揭示了多个景观因素对基因流动的显著影响：
- •
阻碍因素：碎石质土壤母质是预测遗传距离最有影响力的变量，与较陡的坡度高度相关，共同构成了基因流动的主要障碍。此外，主要道路和较高的年降雪量也显著增加了遗传距离，即阻碍了基因流动。
- •
促进因素：有机质土壤母质（如湿地周边）和冲积土壤母质与更低的遗传距离相关，促进了基因流动。
- •
距离作用：地理距离本身也是一个重要影响因素，超过约100公里的距离会导致遗传分化显著增加。
- •
未入选因素：与预期不同，草地/灌丛/湿地等土地覆盖类型以及冠层覆盖度并未被最终模型选中，表明它们在本研究中对于解释遗传连通性的相对贡献较小。
模型结果显示，一条主要道路对基因流动造成的阻力，大约相当于100公里地理距离带来的效应。分别对西北部和东南部区域建模的结果表明，有机质和湖积土壤的促进作用主要反映在Cariboo地区及该区与东南部地区的连接上，而主要道路的负面影响在连接Cariboo与东南部地区的样带中最为强烈。
3.3 预测遗传连通性
基于顶级模型生成的景观阻力面，研究者通过阻力核分析绘制了遗传连通性预测图。地图显示，在Cariboo地区，由于广泛分布的有机质和冲积土壤以及较平缓的地形，存在大片低阻力、高连通性区域。相比之下，在Thompson、Okanagan和Nicola地区，低阻力区域主要集中在少数狭窄的谷底。许多预测的高连通性廊道与主要公路平行或交叉（例如97号公路沿线），突显了这些自然通道同时也是高死亡率风险区。
4 讨论与结论
本研究强有力地证明，景观属性显著塑造了BC省杰弗逊獾西部边缘种群的空间遗传结构。研究发现，基因流动受到碎石质土壤（与陡坡相关）、陡峭坡度、主要道路和冬季大量积雪的阻碍，而有机质和冲积土壤则有利于连通性。与之前在该物种分布核心区域的研究不同，本研究表明，对于这个边缘种群而言，景观偏好对遗传结构的影响可能超过了其强大的扩散能力，这可能与研究区陡峭的环境梯度有关。
研究结果具有重要的保护意义。首先，它确认了主要道路是基因流动的关键屏障，这支持了在该物种恢复计划中优先考虑缓解道路影响的必要性。预测的遗传连通性地图可以直接指导保护实践，例如：1）在高连通性廊道与公路交叉的“热点”区域，优先建设或改造野生动物地下通道、设置防护栏，以降低路杀风险、提高道路渗透性；2）识别并优先保护那些狭窄但至关重要的谷底栖息地廊道，特别是在东南部地区，以维持区域间的基因交流。
研究也指出了当前样本的局限性（大多来自道路附近，可能存在性别偏差和空间聚集），并建议未来应在远离道路的区域、以及东南部种群与相邻的美国华盛顿州种群交界处增加遗传采样，以完善连通性预测，并识别连接更大、遗传多样性更丰富种群的潜在廊道。
总之，这项研究通过创新的机器学习景观遗传学方法，为理解和管理濒危边缘种群的遗传连通性提供了清晰的科学蓝图。它不仅揭示了影响美洲獾基因流动的具体自然与人为因素，更重要的是，其产出的高分辨率空间预测地图，能够帮助保护工作者和管理者将有限的资源投入到最能有效促进基因流动、从而增强种群长期存活能力的关键区域和行动上。
