《One Health》：Host diversity and landscape structure drive genospecies-specific Lyme disease risk across Eurasia
编辑推荐：
莱姆病（Lyme disease）对全球健康构成日益严峻的挑战，部分原因在于气候变化和土地利用转型，而对其传播机制的机理性认识仍十分有限。研究人员分析了来自43个欧亚国家的流行率记录，以绘制病原体生物地理学图谱，并应用理论驱动的结构方程模型（structura
莱姆病（Lyme disease）对全球健康构成日益严峻的挑战，部分原因在于气候变化和土地利用转型，而对其传播机制的机理性认识仍十分有限。研究人员分析了来自43个欧亚国家的流行率记录，以绘制病原体生物地理学图谱，并应用理论驱动的结构方程模型（structural equation modeling, SEM）评估气候、土地覆盖、景观结构和宿主生物多样性对病原体流行率的直接、间接及级联影响路径。基因种（genospecies）表现出由媒介-宿主特异性塑造的独特生物地理格局，其中欧洲传播以篦子硬蜱（Ixodes ricinus）为主，亚洲则以全沟硬蜱（I. persulcatus）为主。Borrelia afzelii流行率与气候、森林破碎化和景观异质性相关，并在大陆尺度上表现出与啮齿动物物种丰富度的负相关关系。Borrelia garinii未表现出宿主多样性效应，但直接响应温度和景观结构，反映出其对高度移动性鸟类宿主的依赖。这些发现揭示了不同基因种因宿主-贮主生态学差异而产生的根本性传播异质性，对疾病控制具有重要启示意义。有效的监测必须将基因种特异性生态学与景观管理相结合，以协调生物多样性保护、气候适应和公共卫生。
论文解读
一、研究背景与科学问题
莱姆病（Lyme disease）是由伯氏疏螺旋体广义种（Borrelia burgdorferi sensu lato, B. burgdorferi s.l.）引起的蜱传疾病，是北半球最常见的蜱传疾病，也是全球公共卫生日益关注的焦点。每年美国约有47.6万人被诊断并接受治疗，西欧每年估计发生23.2万例，亚洲也正出现新的传播疫源地。这些格局与人为环境变化密切相关，后者可改变蜱、宿主与人类之间的相互作用，促进媒介种群扩张和病原体循环。
Ixodes属蜱是B. burgdorferi s.l.的主要媒介，其幼虫、若虫和成虫三个生活阶段均需在脊椎动物宿主上吸血，这一专性吸血行为使蜱-宿主相互作用成为疾病生态学的核心。气候变暖可影响蜱的季节性活动和分布边界，同时改变宿主丰度与群落组成；土地利用变化（包括森林破碎化和农业集约化）重塑景观格局，调节蜱-宿主接触和病原体传播动态。然而，这些环境效应的方向和强度可能因病原体基因种而异。B. burgdorferi s.l.包含20多个基因种，其中B. afzelii和B. garinii是欧亚大陆主要的人类病原体。B. afzelii主要依赖啮齿动物宿主，常与皮肤病相关；B. garinii利用鸟类宿主，与神经系统感染相关。二者不同的宿主关联表明其传播的生态途径不同。
尽管已有大量研究，关键不确定性仍然存在。关于B. burgdorferi s.l.生态学的研究仍较为分散，很少整合跨大陆的发现。现有病原体流行率驱动因素研究多局限于局部到区域尺度，侧重于非生物因素（如气候变化），仅少数研究纳入实地调查获得的宿主群落数据等关键生物因素，其结果难以推广到更广泛的地理范围。机制仍未阐明，现有研究很少区分不同环境因素的直接效应与间接及级联效应。宿主生物多样性增加了另一层复杂性：在北美，生物多样性的稀释效应（dilution effect）得到强有力实证支持，即大量非适宜宿主限制了白足鼠（Peromyscus leucopus）等适宜贮主中的感染；但稀释效应假说在欧亚大陆仍存争议。复杂的宿主-媒介-病原体系统阻碍了将生物多样性纳入疾病监测框架。
鉴于上述问题，本研究旨在绘制B. burgdorferi s.l.及主要致病基因种在欧亚大陆的分布图，评估不同基因种与优势蜱属的关联，并利用结构方程模型（SEM）分析基因种特异性传播路径，量化气候、土地覆盖、景观和宿主多样性对蜱中病原体流行率的直接、间接和级联效应。宿主多样性根据基因种特异性生态学进行评估，区分与B. afzelii相关的小型哺乳动物群落和与B. garinii相关的鸟类群落。
二、主要技术方法
研究人员使用了一套综合性的B. burgdorferi s.l.流行率数据集，整合了2000—2023年间来自43个欧亚国家522篇出版物的2528条流行率记录。潜在驱动因素分为气候、土地覆盖、景观和宿主多样性四个域，所有变量重采样至统一的5弧分空间网格（赤道处约10×10 km）。气候变量包括最湿季度温度和最暖季度降水量，来自WorldClim数据库（2.1版）；土地覆盖变量包括森林覆盖百分比和建成区百分比，来自EarthEnv数据库；景观指标包括森林破碎化指数（forest fragmentation index, FFI）和景观多样性（Shannon多样性指数, SHDI）；宿主多样性按分类群分层，包括鸟类宿主、小型哺乳动物宿主和中大型哺乳动物宿主，物种丰富度通过IUCN红色名录物种分布多边形空间叠加计算。结构方程模型采用最大似然法估计参数，因数据偏离多元正态性，使用Bollen-Stine自助程序评估模型拟合，从2000次重抽样（随机种子=1）估计偏差校正的95%自助置信区间。根据生态理论和已发表证据指定先验路径模型，拟合后依次移除不显著路径（P>0.05）获得简约最终模型。分别构建了B. burgdorferi s.l.、B. afzelii和B. garinii的结构方程模型，模型拟合指标均显示可接受的拟合度。
三、研究结果
3.1 大陆尺度监测与流行格局
数据集包含来自43个欧亚国家的2528条流行率记录，报告涉及从国家到第四级行政区的五个行政级别。西欧和中欧国家主要报告精细尺度数据，而亚洲国家空间分辨率较粗，反映区域蜱监测和研究能力的差异。采样强度极不均匀，中欧和西欧记录密度最高，亚洲以中国贡献记录最多，但中亚和北亚大片区域采样不足。国家平均B. burgdorferi s.l.流行率差异显著，马来西亚和摩尔多瓦最高。
3.2 关键媒介与基因种的生物地理学
媒介-病原体关联在不同蜱属间差异显著。Ixodes属为主要媒介，占阳性记录的大多数，在欧亚广泛分布；Haemaphysalis属感染模式较局限，集中于东亚尤其是中国；Dermacentor、Rhipicephalus和Hyalomma属阳性记录极少且流行率很低，在该区域Borrelia传播中作用次要。基因种表现出独特的生物地理格局：B. afzelii在中欧和西欧占优势，东亚偶发；B. garinii地理范围最广，从西欧经俄罗斯延伸至东亚，中国东北部流行率较高；B. burgdorferi s.s.主要局限于西欧，B. valaisiana、B. lusitaniae和B. spielmanii分布受限。在18个已鉴定基因种中，B. afzelii和B. garinii最为普遍和广泛。B. afzelii在篦子硬蜱（I. ricinus）中流行率最高（2–5%），B. garinii在I. pavlovskyi中达10–20%。共感染在Ixodes物种中发生，I. ricinus携带的基因种多样性最高。
3.3 气候与宿主丰富度驱动总体流行率
在大陆尺度上，B. burgdorferi s.l.流行率受气候和小型哺乳动物丰富度直接影响。最湿季度温度有正直接效应（β=0.172, P<0.001），小型哺乳动物丰富度亦然（β=0.139, P<0.001）。鸟类和中大型哺乳动物丰富度无显著直接效应。温度通过小型哺乳动物丰富度产生正间接效应（间接效应：0.071）。土地覆盖和景观变量不直接影响流行率，但通过森林破碎化或宿主丰富度产生小幅度间接效应（间接效应：0.041）。
3.4 B. afzelii：与啮齿动物丰富度负相关
B. afzelii表现出独特生态学，具有多条直接路径。小型哺乳动物丰富度对流行率有显著负直接效应（β=?0.197, P<0.001）。气候也产生正直接影响：最湿季度温度（β=0.149, P<0.01）和最暖季度降水量（β=0.132, P<0.01）。景观结构具有对比效应：森林破碎化指数（FFI）正向影响流行率（β=0.108, P<0.05），而Shannon多样性指数（SHDI）呈负相关（β=?0.132, P<0.01）。气候和土地覆盖通过双路径产生间接效应：一是通过小型哺乳动物丰富度和景观多样性的负间接效应，二是通过森林破碎化的正间接效应。与总体B. burgdorferi s.l.模型相比，B. afzelii的小型哺乳动物效应由正转负，并受景观结构直接调控。
3.5 B. garinii：气候-景观控制，宿主效应最小
B. garinii传播主要与温度和景观多样性相关，未检测到宿主多样性的直接效应。流行率随最湿季度温度（β=0.142, P<0.001）和景观多样性（β=0.250, P<0.001）增加。降水量和所有宿主丰富度均无直接效应。土地覆盖仅间接影响B. garinii：较高森林面积增加破碎化（β=0.311, P<0.001），后者降低景观多样性（β=?0.145, P<0.001），进而通过SHDI影响流行率。与B. burgdorferi s.l.和B. afzelii不同（宿主群落组成发挥核心调节作用），B. garinii在大陆尺度的动态主要由非生物和景观特征塑造，宿主多样性仅作为气候-景观级联中的中间节点，与流行率无直接关联。
四、讨论总结与结论翻译
讨论总结：本研究表明，B. burgdorferi s.l.在欧亚大陆的传播生态学在不同基因种间存在显著差异，尤其是B. afzelii和B. garinii之间。B. afzelii在中欧和西欧占优势，而B. garinii具有跨大陆传播特征，这些格局与病原体-媒介-宿主特异性密切相关。B. afzelii流行率与气候、宿主多样性和景观格局等多种因素相关；B. garinii流行率主要与气候和景观格局相关，未检测到与宿主多样性的直接关联。啮齿动物物种丰富度与B. afzelii流行率在大陆尺度呈显著负相关，可能通过多样性介导的接触稀释和种间竞争抑制优势贮主等途径实现。森林破碎化通过相反的直接和间接路径影响B. afzelii流行率，其中密度介导的扩增效应超过多样性介导的抑制效应。B. garinii流行率与景观多样性正相关，但与局部宿主多样性无显著关联，这可能与鸟类宿主的高移动性及其跨位点整合暴露有关。所有基因种流行率均与最湿季度温度正相关，B. afzelii还额外与最暖季度降水量正相关。研究支持从泛化监测转向基因种特异性管理策略，并强调将基因种特异性生态学信息与景观管理相结合以协调生物多样性保护、气候适应和公共卫生。研究局限性包括空间采样不均、未测量因素可能影响路径估计、横断面观察设计无法进行因果推断等。
结论翻译：总之，B. burgdorferi s.l.流行率的生态关联因基因种而异。B. afzelii流行率与气候、啮齿动物多样性和景观配置相关，而B. garinii流行率主要与温度和景观多样性相关，在所分析的空间尺度上未检测到与宿主丰富度的直接关系。这些差异表明，未来莱姆病监测应将基因种特异性生态学信息与景观管理相结合，以更好地协调生物多样性保护、气候适应和公共卫生优先事项。
