静脉注射吸毒人群中治疗依从性对丙型肝炎传播动力学影响的建模

《Mathematical and Computational Applications》：Modeling the Impact of Treatment Adherence on Hepatitis C Transmission Dynamics Among Intravenous Drug Users

　　

丙型肝炎病毒（HCV）感染在静脉注射吸毒人群（PWID）中高度流行，治疗效果受依从性和行为复发的强烈影响。研究人员开发并分析了一个确定性分区模型，以考察治疗依从性、复发和风险降低行为对HCV传播动力学的影响。该模型纳入了关键行为路径，包括治疗启动、部分依从性以及向高风险行为的复发。分析结果表明解的正性与有界性，并利用下一代矩阵法推导了基本再生数（R0）。研究人员证明当再生数低于1时疾病无病平衡全局渐近稳定，而当传播超过该阈值时出现一致持续性与全局稳定的地方病平衡。数值模拟与敏感性分析表明，每次共用针头事件的传播概率与共用伴侣数是流行持续的主导驱动因素；相比之下，治疗接受与持续依从性显著减少传播。向药物滥用的复发通过扩大感染池而削弱治疗收益。这些发现凸显了将抗病毒疗法与行为支持整合以实现PWID中HCV长期控制的重要性。

研究背景与动机

丙型肝炎病毒（Hepatitis C virus, HCV）感染在全球约5800万慢性感染者中仍是重大公共卫生负担，静脉注射吸毒人群（people who inject drugs, PWID）估计52%有既往或现症感染，是WHO消除目标的关键群体。尽管直接抗病毒药物（direct-acting antiviral, DAA）使病毒血症人数下降，但无预防性疫苗，且治疗收益受依从性不足与戒毒后复发削弱。既往模型多关注治疗或 abstinence 单一维度，Mushayabasa与Bhunu曾考虑抗病毒治疗、戒断与复发，但未系统量化依从性异质性与行为降级（reduced-risk）通道。为回答“治疗参与及依从相关行为改变如何调控R₀与持久性”，研究人员构建五室确定性分区模型，发表于《Mathematical and Computational Applications》。

主要技术方法

研究人员采用常微分方程（ODE）分区建模，设S（易感）、I（未治感染）、T（治疗后戒断吸毒）、H（治疗后仍注射但降低共用伴侣数）、R（康复）五室，假设均匀混合与常数招募Λ、自然死亡率μ。利用下一代矩阵（next-generation matrix）于无病平衡推导R₀=βc/(μ+γ+α)·(1+... )形式（含治疗分流比例与复发结构）；以Lyapunov函数与LaSalle不变原理证无病与地方病平衡全局稳定性；用Python 3.10的SciPy 1.11中solve_ivp(LSODA)做30年数值积分，初始条件(S,I,T,H,R)=(…)；以偏秩相关系数（PRCC）与超立方采样做参数敏感性，所有参数取文献基线值，无外部队列拟合。

研究结果

1. 2.
   模型构建（Model Formulation）
   研究人员将PWID分为五室，力感染项λ=βc(I+c_H H/N)，其中c_H为H类伴侣获得率折减、κ为H类传播折扣、γ为自发清除率。I以率α求治，比例ρ进T、(1?ρ)进H；T、H分别以率σ_T、σ_H康复；T以率ω复发，比例r回I、(1?r)进H；H以同率ω退行至I。总种群N有界，可行域Ω正不变且吸引。
2. 3.
   模型分析（Model Analysis）
   3.1 无病平衡与R₀：DFE为(S^*,0,0,0,0)=(Λ/μ,0,0,0,0)。F、V矩阵给出R₀=βc/(μ+γ+α)·Φ（Φ含ρ,σ_T,σ_H,ω,r代数组合）。定理1证R₀<1时DFE局部渐近稳定，R₀>1不稳；定理2用比较原理与Gronwall不等式证R₀<1时DFE全局渐近稳定；定理3据边界半流理论证R₀>1时系统一致持久（uniform persistence），存在ξ使liminf I(t)>ξ。
   3.2 地方病平衡：R₀>1时存在唯一内点EE，以R₀显式表出各室。定理4构造Volterra型Lyapunov函数V=Σa_i(X_i?X_i^*?X_i^*ln(X_i/X_i^*))，配EE恒等式与AM-GM不等式得V?≤0，等号仅于EE，由LaSalle原理证EE全局渐近稳定。
3. 4.
   数值模拟（Numerical Simulations）
   PRCC（图2）显示β与c对R₀正向最强（PRCC≈0.8），α负向最强（≈?0.7），μ、γ中等负向，ω、σ弱效。图3等值线：α↑或ω↓压低R₀，σ_T与σ_H协同α可推R₀<1。图4超立方：β、c单增使R₀陡升，α、σ_H抑制最强，ω回升抵消治疗收益。图5变ω：ω↑使12年后I略降但长期感染池增大，T/H池峰值降低，R终值略减，证复发扩充感染池。图6变γ：γ↑微降I峰值，对各室影响弱于α。图7变ρ（依从戒断比例）：ρ↑使H↓、T↑，I与R几乎不动，表明依从性改善治疗内部分流而非直接压低I。图8变r（戒断者回高危比例）：仅12年后H微动，整体影响有限。图9变c：c是最强传播杠杆，c↑使I峰与长期池倍增。图10变α：α↑显著降低I峰与长期水平，放大T/H/R。

讨论与结论翻译

HCV在PWID中持续流行虽DAA有效但仍具挑战。研究人员建模型纳入治疗依从、行为修正与复发，证R₀<1时DFE全局稳、R₀>1时一致持久收束唯一EE。敏感性指认每次共用传播概率β与平均伴侣数c为持久主导正驱，求治率α为最强负驱；治疗收益高度依赖依从与持续行为改变，未戒断或复发者仍以折扣率传毒，异质性削弱人群层收益。c与α为两最大杠杆——c经力感染直接放大，α经缩短未治传染期起效。局限在均匀混合忽略网络聚类与再感染。结论：PWID中HCV控制需超越治疗可及，须依依从、防复发、降风险并行；将抗病毒治疗与危害减少（针具项目、同伴教育）及行为支持整合，方能将R₀压至1下、趋近消除。