《Perspectives in Ecology and Conservation》：Geographic mismatch between venomous snake and treatment center distributions in Brazil: a warning in the face of climate change
编辑推荐：
气候变化正在显著改变全球物种分布，扩大部分物种的分布范围，同时减少另一些物种的适生栖息地。就蛇类而言，这些变化具有直接的公共卫生意义，因为毒蛇物种可能侵入新的区域，包括人口稠密的城市地带。这提高了对适当抗蛇毒血清分配策略的需求，尤其是在偏远地区或医疗基础设施有
气候变化正在显著改变全球物种分布，扩大部分物种的分布范围，同时减少另一些物种的适生栖息地。就蛇类而言，这些变化具有直接的公共卫生意义，因为毒蛇物种可能侵入新的区域，包括人口稠密的城市地带。这提高了对适当抗蛇毒血清分配策略的需求，尤其是在偏远地区或医疗基础设施有限的区域。尽管已有若干研究关注气候变化对蛇类栖息地的影响，但较少有研究探讨此类变化如何影响蛇咬伤治疗中心的物流配置与组织结构。在此背景下，研究人员旨在评估：由气候变化驱动的蛇类分布变化是否需要巴西抗蛇毒血清治疗中心网络（ATCN）的重构。为此，研究整合了爬虫学馆藏中的出现记录、生物气候变量以及土地利用变量，采用5种不同算法构建物种分布模型，并在经缓解情景（SSP 2?4.5）和高排放失控情景（SSP 5–8.5）下考虑2061–2080年的气候变化预测。所建立的模型表明，蛇类适生区域及物种丰富度均将下降，且在SSP 5–8.5情景下尤为明显。研究识别出明显的空间错配：毒蛇适生区的抗蛇毒血清治疗中心网络（ATCN）覆盖缺口最高可达63%，尤其集中于极偏远地区，并以矛头蝮属（Bothrops）和真珊瑚蛇属（Micrurus）最为突出。相反，配备抗蛇毒血清的抗蛇毒血清治疗中心网络（ATCN），尤其是针对巨蝮属（Lachesis）和响尾蛇属（Crotalus）的治疗中心，在最高达92%的毒蛇不适生区域中仍有分布。上述发现表明，有必要结合未来气候变化情景，持续重新评估抗蛇毒血清治疗中心网络（ATCN）的区位布局，以确保对蛇咬伤事件作出更有效的响应；同时也强调，应采用整合式健康路径，将人类、动物与环境健康之间的相互作用纳入考虑，以提升巴西蛇咬伤风险管理水平。
Introduction
文章引言从气候变化与生物多样性压力之间的关系切入，指出多数生物的生存与其气候生态位密切相关，因此温度、降水格局及极端天气事件的变化，不仅威胁生态系统与物种，也对人类生存构成风险。蛇类作为外温动物，对环境温度依赖显著，故对气候变化尤为敏感，既可能直接受到热胁迫影响，也可能因疾病暴发、食物资源减少等间接机制而受冲击。文章进一步指出，全球及巴西范围内已有研究显示，气候变化会促使蛇类适生区发生扩张、收缩或位移，这种变化可能将毒蛇带入新的区域，包括城郊交错地带和城市绿色空间，从而增加人与蛇接触及蛇咬伤事件风险。
在公共卫生层面，文章强调蛇咬伤仍是重要的被忽视热带病问题。世界卫生组织（WHO）估计全球每年约有540万例蛇咬伤事件，其中180万–270万例发生中毒，死亡人数约81,000–138,000。在巴西，每年约发生29,000例蛇咬伤，平均死亡139例。巴西毒蛇主要涉及4个属，即矛头蝮属（Bothrops lato sensu）、响尾蛇属（Crotalus）、巨蝮属（Lachesis）和真珊瑚蛇属（Micrurus），其中矛头蝮属导致的病例占绝大多数。由于不同属毒液所致临床表现差异明显，而治疗效果高度依赖于正确抗蛇毒血清的及时使用，因此抗蛇毒血清治疗中心网络（ATCN）的空间布局具有关键意义。文章还指出，自2014年以来，巴西抗蛇毒血清生产能力受监管与生产设施调整影响而下降，这进一步削弱了医疗服务响应能力。
在研究空白方面，作者指出，以往研究多聚焦于气候变化对毒蛇适生栖息地的扩张或收缩影响，或关注治疗中心可达性如何影响蛇咬伤救治时效，但尚缺乏研究系统分析当前抗蛇毒血清治疗中心网络（ATCN）与毒蛇现实分布之间的空间组织关系，也缺乏对未来气候变化情景下二者空间匹配关系的评估。基于此，本文提出研究目标，即评估气候变化导致的潜在毒蛇分布转移是否需要对巴西抗蛇毒血清治疗中心网络（ATCN）进行重构，并提出当前ATCN在当下及未来均可能无法充分覆盖潜在蛇咬伤风险区的假设。
Material and methods
在数据来源方面，研究整合了多个爬虫学馆藏的数字化记录，包括巴西及阿根廷的重要标本馆资料，并纳入既有的全国蛇类分布数据库。所有记录经过严格筛选，剔除了无采集地点、圈养个体、坐标重复、缺乏地点信息及鉴定错误的样本。对于形态相似且鉴定较难的矛头蝮属（Bothrops）与真珊瑚蛇属（Micrurus）标本，研究特别考虑了鉴定者可靠性与采样地点信息。对于巨蝮属（Lachesis），由于记录缺少采集时间，研究采用保守性数据清理原则，依据历史生物地理学与既有分布证据，排除了若干生态学上不合理或可能代表历史性误差的记录，以降低模型过预测风险。
环境变量方面，研究使用了WorldClim数据库提供的5角分生物气候变量，以及土地利用/土地覆被（LULC）变量，并通过最近邻方法统一空间分辨率。考虑到环境预测变量分辨率约为10 km，研究对出现记录实施10 km空间稀疏化（spatial thinning），以减少空间自相关和局部过采样造成的伪重复。物种建模方面，矛头蝮属（Bothrops）和真珊瑚蛇属（Micrurus）按种分别建模，随后叠加生成属级分布图；响尾蛇属（Crotalus）因巴西仅有1种而直接建模；巨蝮属（Lachesis）则基于合并记录进行建模。对于初始记录数少于80条的Bothrops和Micrurus物种，研究未纳入模型，以保证稀疏化后仍具备足够样本量。
变量筛选采用标准化框架下的物种特异性方案。研究基于既有生态学知识预选变量池，再通过Pearson相关分析与方差膨胀因子（VIF）筛除高相关与高共线性变量，其中阈值分别设为|r| ≥ 0.7和VIF ≥ 3。尽管土地利用/土地覆被（LULC）变量在生态位模型中存在类别属性、时空不匹配和模型复杂度提高等方法学挑战，研究仍将其纳入，以反映森林覆盖、栖息地结构与人为景观转化对毒蛇分布的重要影响，从而避免仅使用气候变量导致对公共卫生风险的刻画不完整。
Modeling protocol
模型构建采用RStudio平台中BIOMOD2程序包的5种算法：广义线性模型（GLM）、广义加性模型（GAM）、随机森林（RF）、梯度提升模型（GBM）和最大熵模型（MaxEnt）。各物种的模型校准区域（M）依据世界自然保护联盟（IUCN）分布多边形确定，并额外设置1°缓冲区，以限制模型校准在生态地理上合理的范围内，同时减少过度外推。对于MaxEnt模型，研究在校准区内随机生成10,000个背景点；对于GLM、GAM、RF和GBM，则随机生成500个伪缺失点，并确保其距已知出现点至少30 km，以降低空间自相关和误将潜在适生点作为缺失点的风险。
模型评估采用重复5折交叉验证，每种算法重复10次，患病率（prevalence）固定为0.5。研究以受试者工作特征曲线下面积（AUC）与真实技能统计量（TSS）评估模型表现，仅保留AUC ≥ 0.8且TSS ≥ 0.7的模型用于集成模型（ensemble models）构建。随后生成连续适宜度图与二值化预测图，二值化阈值采用最大化敏感性与特异性之和的方法。未来投影则基于跨学科气候研究模型（MIROC6），覆盖2061–2080年，在SSP 2?4.5和SSP 5–8.5两种情景下进行。
Post-processing
后处理部分主要通过栅格代数分析当前与未来的二值化适生图，量化栖息地增益、损失及属级物种丰富度变化。研究进一步将预测得到的毒蛇适生区与巴西卫生部公开的抗蛇毒血清治疗中心网络（ATCN）地理信息进行叠加，从而量化两类核心指标：其一，当前及未来情景下无ATCN覆盖的毒蛇适生区面积；其二，位于毒蛇不适生区内的ATCN数量。
为估算2 h内可达的ATCN服务范围，研究依据巴西地理统计局（IBGE）的地理可达性指数，将市镇分为极偏远、偏远、可达和高可达4类，并通过Google Maps随机抽取不同类别市镇之间的线路，估计平均出行距离与时间。在此基础上，为不同可达性等级分别设定40 km、80 km、100 km和120 km的缓冲区，用于近似表示2 h内可达范围，进而评估各可达性水平下毒蛇适生区的ATCN覆盖缺口与不适生区中的过剩覆盖。
Results
研究共汇编17,727条毒蛇出现记录，最终对4个属共33个物种进行了建模，包括18种矛头蝮属（Bothrops）、12种真珊瑚蛇属（Micrurus）、2种巨蝮属（Lachesis）和1种响尾蛇属（Crotalus）。所有集成模型均显示出令人满意的预测性能。总体而言，未来气候情景下4个属均表现出适生区缩减趋势，尽管少数物种出现有限范围扩张，但整体物种丰富度与适生面积均下降，且在SSP 5–8.5情景下更为显著。
结果表明，当前巴西ATCN分布显著偏向人口密集区，导致北部、东北部与中西部多个毒蛇适生区存在明显覆盖空缺。对矛头蝮属（Bothrops）而言，配备抗Bothrops血清（SAB）的ATCN在适生区中的覆盖不足最为突出，在高可达地区未覆盖比例约为6%，而在极偏远地区可达63%，且未来情景下该格局基本持续。真珊瑚蛇属（Micrurus）同样显示较高缺口，当前及未来情景下，在极偏远地区未覆盖比例超过55%。相比之下，巨蝮属（Lachesis）适生区的未覆盖比例相对较低，但配备抗Lachesis血清（SAL）的ATCN大量位于不适生区，平均可达84%，且在更严峻的气候情景下趋于加剧。响尾蛇属（Crotalus）适生区的覆盖状况最佳，平均仅约1%的适生区未被配备抗Crotalus血清（SAC）的ATCN覆盖，但不适生区中的ATCN覆盖比例却较高，当前约31%，在最极端气候情景下可升至37%。
此外，结果显示，未来气候变化不仅会延续现有错配，还会进一步强化这种错位格局，即毒蛇适生区缩减的同时，越来越多ATCN将位于气候上不再适宜相关毒蛇分布的区域。变量重要性分析进一步表明，气候因子对适生性格局的解释力总体高于土地利用变量，其中等温性和降水相关变量最为关键，而人为土地利用变量贡献中等，提示若干毒蛇物种对改变环境的耐受性有限，而非普遍适应人类改造景观。
Discussion
讨论部分强调，本研究的核心贡献不止于揭示气候变化将导致毒蛇适生区与物种丰富度下降，更在于证明这些生态变化会重组蛇类风险空间，并加剧其与巴西抗蛇毒血清治疗中心网络（ATCN）之间原已存在的空间错配。通过整合生态适宜性、地理可达性与医疗基础设施配置，研究揭示了抗蛇毒血清供给与生态风险之间的结构性失衡，而这一失衡在未来气候情景下预计将进一步扩大。
文章指出，在当前与未来情景下，北部、东北部和中西部等以亚马孙和塞拉多（Cerrado）生物群系为主的地区，适生蛇类区域长期缺乏足够的ATCN覆盖，而相当数量ATCN又持续分布于未来将变得不适宜相关毒蛇生存的区域，反映出公共卫生资源配置效率不足。尽管响尾蛇属（Crotalus）表面上显示较高的覆盖匹配度，但其不适生区内仍存在大量治疗中心，说明即便看似较优的情形之下，结构性配置问题仍然存在。
从公共卫生后果看，这种空间错配对偏远和社会经济脆弱人群影响尤为严重。蛇咬伤发生后，患者至医疗机构的旅行时间延长，可能增加重症并发症及死亡风险，也会加重地方医疗系统负担。研究据此呼吁，巴西应尽快重新评估抗蛇毒血清分配策略，将生态风险较高且社会脆弱性较强的区域列为优先对象，并将气候变化预测纳入卫生规划，以避免服务缺口扩大。
文章同时指出，在无医学重要毒蛇分布的区域持续配置ATCN及抗蛇毒血清，不仅意味着公共资源浪费，也会挤占真正高风险地区的供给能力。因此，抗蛇毒血清配置应建立在稳健的流行病学证据基础之上，综合考虑毒蛇地理分布、蛇咬伤发生格局和医疗机构响应能力，并通过以数据驱动的公共政策与治理机制提升资源配置透明度与效率。
讨论还从生态系统与产业层面扩展了研究意义。蛇类在调控小型脊椎动物、特别是啮齿动物种群方面具有关键作用，其衰退可能造成营养级失衡并增加人畜共患病风险。同时，蛇毒也是抗蛇毒血清生产和生物医学研究的重要资源，因此蛇类种群减少还可能影响制药与生物技术领域。文章进一步提醒，蛇咬伤风险并不与自然分布完全重合，人口稀疏区漏报、人工饲养环境中的意外事件等因素，都会使生态适生性与观测发病率之间的关系更加复杂，因此需要综合生态、社会与基础设施维度的整合研究框架。
最后，作者审慎说明了研究局限，包括部分Bothrops和Micrurus物种因记录不足未被纳入建模、极端气候事件未被显式模拟，以及出现记录中存在一定时间信息缺失等。但作者认为，这些限制并不改变研究揭示的广域空间趋势与结构性脆弱性。文章最终提出，应以“同一健康”（One Health）理念统筹人类、动物与环境健康，在持续监测、适应性政策制定和跨部门协作基础上，建立更具韧性、更公平且面向未来的蛇咬伤防控与救治体系。
