条条大路通感染：解析野生有蹄类种群中布鲁氏菌病传播路径的机制模型

《Functional Ecology》：All roads lead to infection: A mechanistic model to assess transmission routes of brucellosis in a wild ungulate population

　　

摘要 野生动物与家养种群常共享病原体，存在跨物种传播风险，对动物健康、人类健康、经济发展及生物多样性保护构成重大威胁。量化自然种群中的疾病传播过程是理解传染病动态的核心，但传播过程难以直接观测，其机制评估极具挑战。由马尔他布鲁氏菌（Brucella melitensis）和流产布鲁氏菌（Brucella abortus）引起的布鲁氏菌病是全球最重要的人兽共患病之一，可感染家养与野生反刍动物。研究人员以阿尔卑斯羱羊（Capra ibex）种群为对象，探究马尔他布鲁氏菌的传播动态——该种群是首个已知可维持此病原体的野生宿主。研究旨在量化种群内的疾病传播，拆解流产产物、感染性分娩产物、性传播及垂直传播四条替代传播路径的贡献。研究人员构建了空间显式、分性别与年龄结构的个体为本模型（IBM），模拟宿主种群动态及易感与感染个体接触导致的病原传播，结合13年的宿主种群与疾病监测数据拟合模型，估算四条传播路径的相对重要性，量化各路径的 transmission 占比，并区分年龄与性别特异性贡献以估算有效再生数（Reff）。结果显示，流产产物是马尔他布鲁氏菌传播的主要驱动因素，占所有传播的57%–85%；其余路径均为次要作用，其中性传播的占比存在不确定性（4%–25%）。雌性个体在疾病维持中起核心作用，有效再生数远高于1（2.1），而雄性的贡献可忽略不计（<1，0.4）。研究结果为制定靶向管理措施提供了关键依据：由于多数传播通过流产产物发生，靶向清除近期感染的雌性可将有效再生数降至1以下，推动病原消亡；此外，流产高峰预计出现在晚春，旨在限制羱羊向包括家畜在内的其他物种传播的控制措施应优先在该高风险时段实施。

论文解读

研究背景与意义

布鲁氏菌病是由布鲁氏菌属（Brucella）细菌引起的全球流行人兽共患病，可感染家养与野生反刍动物，导致流产、繁殖力下降，每年约造成200万人类感染病例，且无安全有效的人用疫苗。当前，野生动物与家畜的疫病界面传播已成为生态保护与公共卫生领域的重点关注问题，但野生种群的疾病传播机制难以直接观测，且家养动物的流行病学参数无法直接外推至野生种群——野生有蹄类的社会行为呈季节性聚散的动态变化，性别与年龄结构、繁殖模式均与家养群体存在显著差异。阿尔卑斯羱羊（Capra ibex）种群是首个被证实可自我维持马尔他布鲁氏菌（B. melitensis）的野生宿主，自2012年在法国阿尔卑斯山巴尔吉地块发现以来，已开展十余年监测与管理，但针对其传播路径的定量解析仍十分匮乏。明确该种群中布鲁氏菌的传播机制，既是实现欧盟布鲁氏菌病根除目标的科学基础，也是降低野生种群向家畜传播风险的关键前提。该研究发表于《Functional Ecology》，通过整合长期监测数据与机制模型，首次量化了野生羱羊种群中四条布鲁氏菌传播路径的相对贡献。

关键技术方法

研究基于法国巴尔吉地块2012–2024年的阿尔卑斯羱羊种群监测数据，涵盖种群数量标记重捕估计、个体血清学检测、空间分布记录及历年扑杀管理操作档案。核心方法是改进已有的空间显式个体为本模型（IBM），新增流产与感染性分娩的传播概率区分、性别特异性性传播模块及经产雌性的感染后流产概率校准，模型包含宿主种群动态、病原传播、管理措施三大模块，区分个体的性别、年龄、健康状态（易感、暴露、感染、恢复）及所属地理分区。采用近似贝叶斯计算（ABC）方法拟合模型，估算流产暴露传播概率、感染性分娩暴露传播概率、主导雄性与追逐雄性的性传播概率及种群环境容纳量四个关键参数。基于拟合后的模型，分别通过年度传播事件统计与有效再生数（R_eff，即单个感染个体在其传染期内产生的二代病例数的期望值）计算，量化各传播路径及不同年龄-性别组的贡献，并通过调整垂直传播概率、流产后恢复抗体持续时间、流产暴露群体规模等参数开展敏感性分析验证结果稳健性。

研究结果

参数估计

后验分布显示，流产暴露的传播概率为0.008（95%可信区间：0.002–0.019），感染性分娩暴露的传播概率为0.007（0–0.019），两者无显著差异但感染性分娩的缩放系数α的后验分布偏零，提示其传播效率可能更低；主导雄性与雌性交配的性传播概率为0.154（0.006–0.599），追逐雄性与雌性交配的性传播概率为0.069（0.003–0.270）；种群环境容纳量为554（451–684）。各参数间相关性较弱，仅流产与感染性分娩传播概率呈负相关（-0.69），提示两者存在补偿效应。

模型性能

后验预测检验表明，模型可准确复现2012–2024年的种群数量动态与血清阳性率变化：成功模拟了2013、2015年大规模扑杀导致的种群骤降，2016年后的逐步恢复，以及血清阳性率从2013年的45%持续下降至2018年后趋于平稳的趋势，与管理记录吻合。

传播路径与个体贡献的重要性

各路径的平均传播占比为：流产66%、性传播17%、感染性分娩12%、垂直传播5%，其中流产的贡献在模态估计下可达85%。有效再生数层面，流产的R_eff为0.8（0–6），是维持病原的核心路径；感染性分娩、性传播、垂直传播的R_eff分别为0.13（0–2）、0.31（0–2）、0.05（0–1），均低于1。个体贡献呈显著性别差异：雌性贡献92%，R_eff为2.1（0–8）；雄性贡献仅8%，R_eff为0.4（0–3）。年龄差异相对较小：年轻雌性（2–5岁）R_eff为2.4（0–9），老年雌性（≥6岁）为1.9（0–8）；年轻雄性R_eff为0.3（0–3），老年雄性为0.5（0–3）。

结果稳健性

调整垂直传播概率（从10%升至30%）、流产后恢复抗体持续时间（从终身阳性改为平均6年阳性）、流产暴露群体规模（从全分区个体限制为60只）均未改变流产为首要传播路径的结论；移除流产与感染性分娩路径会显著降低模型拟合度，仅移除流产路径则会导致感染性分娩传播概率升高8倍以维持拟合，进一步验证了水平传播的核心作用。

讨论与研究结论

讨论部分指出，野生羱羊的布鲁氏菌传播模式与家养反刍类一致，流产产物的环境高载量与晚春种群聚集共同驱动了传播优势；感染性分娩的传播效率更低，与分娩后雌性的短期隔离行为有关。性传播虽占比不确定且不足以单独维持病原，但雄性的高移动性可能是病原扩散至邻近种群的关键机制。垂直传播的贡献可忽略，与野牛的研究结果一致。个体贡献的性别差异与野牛、马鹿的研究吻合，但年龄偏向老年个体的特征不同于野牛的幼龄高发，源于羱羊的流产暴露无年龄限制。研究支持布鲁氏菌的水平传播优势是其在群居野生有蹄类中建立的地方性流行的一般规律，但宿主特异性免疫差异可能导致传播效率变异——阿尔卑斯羱羊因历史瓶颈导致遗传多样性低，免疫调控能力弱，可能比家养山羊更易感染与排毒。

研究最终结论为：在目标阿尔卑斯羱羊种群中，流产相关的水平传播是马尔他布鲁氏菌的首要传播路径，感染性分娩、性传播与垂直传播为次要路径。尽管次要路径的贡献存在一定不确定性，但不影响管理建议的制定：针对晚春妊娠雌性的靶向清除可将R_eff降至1以下，推动病原消亡；针对家畜的跨界传播防控应聚焦于晚春的放牧初期，此后风险显著降低。