**安第斯汉坦病毒在长尾稻鼠中栖息地介导的溢出风险：一种机制性生态流行病学模型**

《Viruses》：Habitat-Mediated Spillover Risk of Andes Hantavirus in Oligoryzomys longicaudatus: A Mechanistic Eco-Epidemiological Model

　　

**摘要**

安第斯汉坦病毒（Andes hantavirus, ANDV）是一种啮齿动物传播的正汉坦病毒（orthohantavirus），与南美洲南部的汉坦病毒心肺综合征相关。其在自然宿主中的维持及溢出风险取决于其主要宿主动物——长尾稻鼠（*Oligoryzomys longicaudatus*）的生态学，以及改变栖息地可用性、宿主丰度和人–鼠接触界面的环境变化。研究人员构建了一个机制性生态流行病学模型，将有效栖息地覆盖度与易感–感染–移除（SIR）框架耦合，用于模拟ANDV在长尾稻鼠中的传播。栖息地退化与补偿性恢复可改变啮齿动物的环境容纳量、出生率及感染压力，其中感染压力取决于感染宿主负荷相对于瞬时生态容纳量的比值。研究人员推导了栖息地平衡态、基本再生数（basic reproduction number, R₀）和时变有效再生数（effective reproduction number, R_t），并在退化–恢复情景下评估了感染负担和阈值指标。分析表明，R₀独立于平衡栖息地覆盖度，因为在无病平衡点处易感者丰度随环境容纳量成比例缩放；相反，R_t、累积发病率和感染负荷则依赖于栖息地介导的瞬时拥挤效应。恢复会增加宿主动物丰度和绝对感染负担，而退化虽可降低丰度，却会增强拥挤驱动的传播压力，并延长超临界窗口期。上述结果界定了栖息地变化可能加剧宿主动物感染压力的生态条件，并为人类–啮齿动物界面的人畜共患监测提供指导。

**论文解读文章**

**1. 研究背景与问题**

人畜共患病的出现常与人类活动导致的栖息地扰动密切相关。砍伐森林、栖息地破碎化等过程可改变野生动物储主、媒介、家养动物与人类之间的接触频率，从而重塑病原体传播风险。啮齿动物是多种人畜共患病原体的重要储主，其传播方式包括直接接触、环境污染或向相关宿主种群的溢出。正汉坦病毒属（*Orthohantavirus*）是典型代表，其主要储主属于鼠科（Muridae），感染通常呈持续无症状状态，受染个体通过尿液、粪便和唾液排毒。人类主要通过吸入气溶胶化的排泄物感染，在美洲主要表现为汉坦病毒心肺综合征，具有较高病死率。安第斯汉坦病毒（ANDV）尤其特殊，已被证实具备持续人际传播能力，对感染控制和医院隔离具有直接意义。

尽管人为改变的栖息地可能促进传播和病原体持续存在，但其机制难以分离。干扰可能增加接触率、减少栖息地可用性、改变资源分布或重塑宿主群落。在啮齿动物传播系统中，这些过程可改变宿主丰度、局部拥挤程度、扩散行为以及与人接触的频率。既有经验研究支持栖息地减少与传播增强之间的关联：强制聚集可提高病原体流行率，破碎化可增加辛诺柏病毒（Sin Nombre virus）感染率，储主密度可调节汉坦病毒传播。然而，多数汉坦病毒模型仅强调储主种群动态，未将栖息地作为动态状态变量；多物种兽疫模型虽已提出，但多数未严格区分栖息地丧失、密度调节和传播拥挤。这一空白限制了模型对恢复与退化情景作为明确管理干预措施的评估能力。

**2. 研究目标与模型框架**

为弥补上述空白，研究者开发、分析并模拟了一个生态流行病学模型，用于描述长尾稻鼠（*Oligoryzomys longicaudatus*）种内ANDV传播。长尾稻鼠是阿根廷和智利的主要汉坦病毒储主，属机会主义泛化种，具有广泛生态可塑性，从湿润温带森林到草原交错带均有分布，在具有灌木覆盖的茂密林下植被中可达到高丰度。其丰度可在竹子开花、作物脉冲或景观干扰后增加，因此栖息地介导的种群统计学和接触结构对ANDV传播动态至关重要。

模型将有效栖息地覆盖度视为动态生态状态变量，它调节储主环境容纳量、出生率和传播压力。栖息地变化通过自然更新、人为退化和补偿性恢复实现。自然更新遵循逻辑斯蒂增长，退化以速率q作用于可用栖息地，恢复以速率p作用于栖息地赤字。啮齿动物种群分为易感（S）、感染（I）和移除（R）三类。出生率遵循依赖密度的调节，感染压力与感染宿主负荷相对于瞬时生态容纳量的比值成正比，即拥挤负荷项。该线性逆形式隔离了栖息地拥挤机制并保持阈值分析可处理。研究者将β视为有效传播系数，聚合了直接接触与通过环境沉积的间接传播路径；未引入独立环境病毒库，以保证模型可辨识性。

**3. 主要结果与分析**

模型解析结果呈现标准流行病学阈值结构：当R₀<1时无病平衡点局部渐近稳定；当R₀>1时地方性平衡点唯一存在且局部稳定。基本再生数R₀的表达式为β·(1/(d+γ))·((b/(d))?1)·(K0/K0)…实际上，其关键性质是R₀不依赖于平衡态栖息地覆盖度F*，因为在无病平衡点处易感者丰度与携带能力成比例缩放，栖息地相关项相互抵消。相反，时变有效再生数R_t、累积发病率和感染负荷依赖瞬时拥挤比I/(S+I+R)与K(F)的关系。数值实验在退化（q）–恢复（p）平面上揭示了绝对感染负担与拥挤介导风险之间的权衡。

研究首先分析流行病学指标对栖息地管理的依赖性。累积发病率（C_inc）和平均发病率（I_avg）在高恢复、低退化的情景下达到最大，例如在p=0.9、q=0.05时，累积发病率最大；此时储主种群接近更高环境容纳量，绝对感染负担增加。相反，平均感染负荷相对生态容纳量（L_avg）在高度退化且无恢复的情景下最大，例如q=0.8且p=0时，该值达到2.73。这一数值并非流行率，而是表示感染丰度瞬时超过环境容纳量。这表明退化虽降低了宿主动物丰度，但通过收缩有效栖息地而使剩余栖息地中的感染个体浓度升高。

其次，有效再生数及其时间持续性分析显示，R_t和峰值有效再生数R_peak在高退化、无恢复条件下达到最大，例如在q=0.8、p=0时，R_t≈1.35，R_peak≈1.37，表明中等超临界传播可持续较长时间。处于超临界状态的时间比例T_sup在q=0.5附近达到约0.99，即使退化程度较低也存在持续的传播潜力。时变分析表明，栖息地退化使R_t在瞬态过程中升高，而恢复则降低拥挤压力。所有数值均固定R₀=1.20，差异完全由时间依赖的I/K(F)比值驱动。

敏感性分析显示，R₀对传播系数β的弹性为+1，对出生率b的敏感性随基础值增大而增强，死亡率和移除率则降低R₀。

**4. 讨论与结论**

机制上，模型将病毒持续存在的宿主丰度–负荷权衡分开。恢复提升生态容纳量，支持更大宿主种群，从而增加绝对感染负担；但恢复也降低拥挤压力，使R_t在多数窗口保持亚临界。退化降低宿主丰度，但通过生态拥挤升高相对感染负荷和R_t，延长超临界窗口。该机制与智利长尾稻鼠长期监测中观察到的持续无症状感染、血清阳性率随宿主导度变化、农业改造区(Maule地区)血清阳性率最高以及多年度波动产生瞬时环境污染峰值等定性事实一致。

气候强迫（如竹子/南美山毛榉结实的资源脉冲、ENSO降水异常）可自然进入模型：通过出生率r₀和栖息地再生率ρ变化，或通过温度湿度对病毒存活的影响体现于有效传播系数β。研究人员指出，R₀在无病平衡点处独立于F*，因此季节或年际强迫不会改变入侵阈值的代数形式，而仅重塑时变量R_t、L等瞬态风险窗口的方向和幅度。

关于溢出解释，虽然模型未包含人类环节，但储主水平指标提供了机制基础。持续超临界R_t、高感染负荷比和长超临界窗口等状态标记了储主维持高感染压力的生态条件，提示应加强啮齿动物监测、环境污染检测和人类接触防护。恢复可减少拥挤但增加生物量，退化则降低丰度却提高相对感染负荷，因此将恢复单纯视为保护性干预是片面的。监测应同时追踪丰度、血清学/分子感染指标、栖息地状况和人类–啮齿动物接触，并在地方性校准后再用于预测。

研究结论明确：模型区分了绝对感染负担与拥挤介导的感染压力。恢复增加生态容量并可能通过支持更大的储主种群提高累积感染负担；退化减少宿主丰度但可通过生态拥挤提高感染负荷比和有效再生数。R₀不依赖平衡栖息地覆盖度，因为无病平衡点易感者丰度随携带能力缩放，而瞬态风险指标依赖F。因此，低宿主丰度并不必然意味着低传播压力，高储主丰度也不必然意味着高拥挤性风险。监测项目应整合丰度、感染指标、栖息地条件及人类暴露数据，在人类–啮齿动物界面同时评估栖息地恢复、退化控制、啮齿动物监测和暴露减少措施。