HIV-猴痘共感染动力学中的疫苗依赖性分岔：协同效应与地方性稳定

《PLOS Complex Systems》：Vaccination-dependent bifurcations in HIV-mpox co-infection dynamics: Synergistic effects and endemic stability

　　

2022年全球猴痘（Mpox）暴发揭示了性传播作为主导传播模式，对男男性行为者（MSM）群体造成不成比例的影响，并在HIV地方性流行人群中带来了新的挑战。研究人员首次提出了一个捕捉HIV与猴痘双向相互作用机制的模型，该模型纳入了三项关键创新：（1）HIV诱导的猴痘进展免疫学调控，（2）抗逆转录病毒治疗（ART）依赖性传播率，以及（3）按HIV状态分层的疫苗接种策略。分析表明，HIV共感染产生后向分岔（临界阈值_Rc），使得即使在基本再生数（R₀）<1时猴痘仍能持续存在，这一现象在单病原体模型中并不存在。该模型揭示，未治疗的HIV使猴痘易感性增加2.3倍（95%置信区间：2.1-2.6），而当前疫苗策略在晚期HIV病例（CD4⁺ < 200 cells/mm³）中的效力降低38%。关键的是，研究人员确定了一个最优干预窗口，其中60%的ART覆盖率与靶向疫苗接种相结合，与单独干预方法相比可将共感染患病率降低5.7倍（95%置信区间：5.2-6.3）。这些发现解决了2022年应对工作中的三个关键空白：（i）缺乏共感染特异性传播指标，（ii）未量化的ART-疫苗接种协同效应，以及（iii）HIV分层疫苗效力估计。研究结果提供了HIV-猴痘综合控制的框架，证明协调的检测和预防活动在各种流行病学情景下比序贯干预措施表现优出21%–34%。

该研究聚焦HIV与猴痘在男男性行为者群体中的共感染动力学问题，构建了首个双向耦合的数学模型框架，发表于《PLOS Complex Systems》。研究背景源于2022年全球猴痘暴发的特殊性：性接触成为主要传播途径，而MSM群体中HIV感染率超过25%的地区面临独特的公共卫生挑战。传统模型分别处理HIV传播与猴痘传播，但无法捕捉两种病原体之间的免疫学相互作用，也无法解释为何在标准流行病学指标提示应消除时猴痘仍持续流行的现象。研究人员指出，现有研究存在三方面关键空白：缺乏共感染特异性传播指标、ART与疫苗接种的协同效应未能量化、以及缺少按HIV状态分层的疫苗效力估计，这些不足制约了精准防控策略的制定。

研究人员开展的核心工作包括：建立包含15个互斥状态的房室模型，追踪HIV与猴痘在MSM人群中的共动态；通过下一代矩阵分析推导基本再生数（R₀）并确定稳定性判据；利用2022–2023年暴发数据进行贝叶斯马尔可夫链蒙特卡洛（MCMC）校准；开展全局敏感性与不确定性分析以识别关键驱动参数。研究得出以下主要结论：HIV共感染诱导后向分岔现象，使猴痘在R₀<1时仍能持续；未治疗HIV使猴痘易感性增加2.3倍，ART可恢复接近正常的免疫保护；60% ART覆盖率联合靶向疫苗接种效果优于单独干预5.7倍；协调检测预防活动比序贯干预表现优出21%–34%。

研究所采用的关键技术方法如下：构建15状态房室模型，设置HIV阴性、HIV阳性未治疗、HIV阳性ART治疗三类亚群，每层再细分未接种/已接种状态，纳入免疫抑制易感性增强系数（σ₁、σ₂）和疫苗效力衰减因子；通过谱半径分析计算下一代矩阵得到R₀，利用中心流形理论证明后向分岔存在性，确定临界值_Rc；采用贝叶斯MCMC方法（MATLAB的fmincon和lsqnonlin实现）将模型输出校准至2022–2023年共感染患病率时间序列数据；运用偏秩相关系数（PRCC）分析和基于拉丁超立方采样（LHS）的Sobol全局敏感性分析，计算一阶指数（S_i）和总阶指数（S_Ti）以量化参数独立效应与交互效应；进行数值模拟比较理论情景（双病原体同时引入）与2022年现实情景（猴痘引入HIV地方性流行人群）的动态差异。数据来源为美国公开周报确诊病例（按HIV状态分层）及全球病例系列中的共感染比例数据（约38%–52%猴痘病例为HIV感染者）。

研究结果部分按以下结构展开：

模型构建与假设：建立封闭MSM人群的15状态房室模型，假设均匀混合、疫苗部分保护且免疫随时间衰减、ART恢复免疫功能但可能存在残余损伤、共感染者不对同种病原体重复感染。模型排除垂直传播、年龄结构及非性传播途径，专注于直接传播动力学。

基本再生数与稳定性分析：通过下一代矩阵得到R₀ = max(R₀^Mpox, R₀^HIV)，证明疾病自由平衡点在R₀<1时局部渐近稳定。当HIV诱导的易感性增强系数σ₁超过临界值σ_crit = f(ν, ε)时出现后向分岔，其中ν为疫苗接种率，ε为疫苗效力，这导致R_0c < 1 < R₀ 区域存在双稳态（稳定疾病自由平衡点与稳定地方病平衡点共存）。

参数敏感性分析：PRCC分析显示猴痘传播率β_m（0.42）为R₀最强正向驱动因素，HIV传播率β_h（-0.38）为最强抑制因素。疫苗效力ε因与HIV共感染动态的交互作用呈意外正相关（0.35）。Sobol分析表明β_m和γ_m（进展率）为R₀主导独立驱动因素，而σ₁的总阶指数显著大于其一阶指数，提示存在重要参数交互作用。峰值共感染患病率分析中，HIV诱导易感性因子σ₁ emerges为最关键参数（S_Ti = 0.47），而假设参数（ART加速因子）影响可忽略，证实模型结论对该假设不确定性具有稳健性。

疫苗与ART协同效应：确定关键疫苗接种率阈值ν_c，超过该值可实现猴痘消除。等高线分析揭示R₀=1曲线上疫苗效力与ART覆盖率的非线性交互：当ε ≥ 0.75时，达到R₀<1所需ART覆盖率比无疫苗情景低40%。时序动态显示疫苗接种使峰值共感染患病率降低68%并延迟流行高峰14周。

模型校准验证：利用美国2022年5–10月26周累积发病率数据，通过负二项分布似然函数和自适应Metropolis-Hastings采样进行贝叶斯校准。后验分布显示σ₁均值为2.3（95% CrI: 2.1-2.6），Watanabe-Akaike信息准则（WAIC）表明优于简化嵌套模型，后验预测检查确认无系统性偏差。

数值模拟与情景比较：理论情景与2022年现实情景均呈现相似分岔模式，但后者显示临界疫苗接种阈值ν_c降低15%，提示预先存在的HIV地方性流行降低了消除猴痘所需的疫苗接种努力。时变再生数（_Rt）分析显示，联合干预（ART扩大+疫苗接种）比单一干预提前约67天实现_Rt < 1。

讨论部分重点阐述了研究的公共卫生意义。研究人员指出，后向分岔现象解释了2022年暴发现场观察到的核心矛盾：为何在标准消除阈值下猴痘仍在HIV地方性流行社区持续传播。这一发现颠覆了单病原体模型的预测，要求重新修订考虑共感染率的暴发阈值。ART的免疫恢复效应量化为降低猴痘传播性38%，为"治疗即预防"策略在共感染情境中的应用提供了量化依据。研究识别的最优干预窗口——60% ART覆盖率联合靶向疫苗接种——为资源有限情景下的优先级设定提供了决策支持。

研究结论部分指出：HIV-猴痘共感染产生超线性传播效应，峰值患病率超出单病原体模型预测2.3倍；在接触率β_c = 0.42处识别的分岔揭示共感染者主导暴发轨迹的临界阈值，此现象在现有HIV-only框架中不存在；ART通过免疫恢复降低传播性38%，解释了ART依从性个体猴痘获得率较预测低60%的观察差异；联合干预实现流行控制的速度比相加预期快67天；80% ART覆盖率联合50%疫苗接种可产生五倍于单一干预的暴发削减效果；研究建议实施HIV/猴痘综合检测、CD4⁺ < 350 cells/mm³的HIV感染者疫苗优先接种、以及考虑共感染率的修订暴发阈值。未来研究方向包括年龄结构化接触模式、长期免疫衰减效应和随机超级传播者动力学。该模型的模块化设计使其能够快速适应进化中的病原体，为公共卫生机构提供实时暴发响应优化工具。