编辑推荐：
本研究针对撒哈拉以南非洲（SSA）疟疾控制的瓶颈问题——为何在干季仍存在低度传播热点并引发后续雨季大流行——开展了深入探讨。研究人员以布基纳法索首都瓦加杜古为案例，通过构建融合多类型人口迁移（常规、季节性、永久性）的数据驱动传染病模型，并结合统计推断方法分析。结果揭示，季节性劳工迁移是维持干季低度传播、并促成雨季疫情爆发的关键驱动力。这项研究的意义在于，强调了将人类流动模式（特别是季节性迁移）纳入疟疾传播模型和防控策略的重要性，为SSA地区制定针对移动人群的精准干预措施、最终迈向消除疟疾提供了关键科学依据。
疟疾，这一由受感染雌性按蚊叮咬传播的传染病，依然是全球尤其是撒哈拉以南非洲（SSA）地区最严峻的公共卫生挑战之一。该地区承载了全球超过90%的疟疾病例和死亡负担。尽管全球在防控和治疗上投入巨大，但SSA地区的疟疾传播呈现显著的时空异质性。一个令人费解且棘手的问题是：在炎热干燥的季节，西非SSA的许多城郊和城乡结合部会持续出现低度传播“热点”（每1万人周发病率在3至150例之间）。这些“小火星”常在雨季来临前出现，并往往演变成大规模的疫情爆发。传统的控制手段在高传播区成效显著，但这些在干季持续存在的低度传播热点究竟缘何“顽固不化”，其背后的机制尚不明确。这就像在一片防火重点区域，明火已得到控制，但总有几处暗火在看似不可能的条件下阴燃，并在环境适宜时迅速引燃整片森林。破解这些“暗火”的维持机制，对于最终扑灭疟疾这团“大火”至关重要。
为了解决这个谜题，一项发表在《Journal of Epidemiology and Global Health》上的研究，将目光投向了人类自身的行为——迁移。研究人员猜测，由农业和经济周期驱动的季节性农村-城市劳动力迁移，可能在维持干季传播、并为后续雨季流行病“播下火种”方面扮演着关键角色。为了验证这一假设，他们以布基纳法索首都瓦加杜古为研究焦点，该地的疟疾发病率展现出鲜明的时空分布模式。
为开展这项研究，研究人员整合了多种数据并运用了数理建模与统计推断相结合的方法。首先，他们收集并分析了瓦加杜古市从2011年1月到2023年12月的月度疟疾报告病例数据，这些数据来源于布基纳法索卫生部的健康数据平台。其次，研究利用了先前关于该地区人口流动的横断面调查数据，该数据量化了不同月份的季节性迁移模式。此外，来自布基纳法索国家气象局的月度降雨和温度数据也被纳入分析。基于这些数据，研究团队构建了一个“起源地-目的地”框架下的疟疾传播动力学模型，该模型基于经典的易感-感染-恢复（SIR）仓室结构，并明确纳入了人类迁移（分为常规、季节性、永久性迁移）和蚊媒种群动态。他们采用了小波分析来揭示病例报告与气候因素（温度、降雨）在时频域上的关联性。更重要的是，他们运用最大似然估计（MLE）等统计推断方法，将不同迁移模式（常规、季节性、永久性）的模型与实际的疟疾发病率数据进行拟合和比较，以确定哪种迁移模式最能解释观察到的流行病学模式。
研究结果从多个层面揭示了季节性迁移的核心作用：
3.1 数据分析：通过对瓦加杜古疟疾病例时间序列的小波分析，研究人员发现病例报告存在显著的12个月周期，并且与温度和降雨在1年频带上有强相关性。变化点分析进一步确认，疟疾发病率在干热月份（1-5月）处于低水平，而在湿润的季风月份（7-11月）则显著升高。
3.2 模型分析：通过对构建的数学模型进行定性分析和随机模拟，研究人员深入探究了迁移对疾病负担的影响。
3.2.1 平衡点与基本再生数：他们计算了模型在不同迁移情景下的平衡态和基本再生数（R0）。研究发现，目的地的R0d由两部分组成：一部分由本地因素驱动，另一部分则与迁入移民的规模及其相对于城市居民更高的蚊虫叮咬暴露风险有关。这表明，迁入移民会增加目的地的有效再生数。全局基本再生数R0的计算显示，随着迁入目的地的人口增加，R0最初会下降（因为个体离开了高流行度的起源地），但当迁入量足够大时，目的地的传播风险会急剧上升，导致整个系统超过流行阈值（R0> 1）。然而，有趣的是，当迁入率（η）和迁出率（ν）数值相当时，即使迁移强度很高，最终的R0也可能降至1以下，意味着疾病会全局消亡。0在(η, ν)参数空间上的变化，显示了高、低传播风险区域。">
3.2.2 随机模拟：模拟不同迁移类型对目的地疾病负担影响的数值结果显示，常规迁移会导致病例迅速增加；而季节性和永久性迁移导致的病例增长则更为渐进和渐近。季节性迁移者通常在目的地停留较长时间（2-4个月），并在干热月份大量涌入，这显著提升了传播潜力。模拟还证实，即使目的地蚊虫出生率极低、死亡率很高且初始无感染病例，极低的迁移率也足以维持疾病在目的地的持续存在。
3.3 蚊虫生长的季节性：研究将降雨和温度数据通过插值转化为连续函数，并将其融入蚊虫种群增长模型，以捕捉气候对媒介数量的季节性驱动作用，从而构建了一个季节性强迫模型用于后续的参数估计。
3.4 参数估计与模型选择：研究的关键步骤是比较不同迁移模型（无迁移、常规迁移、季节性迁移、永久性迁移）对实际数据的拟合效果。通过最大似然估计和计算赤池信息准则（AIC），研究人员发现，季节性迁移模型最能准确捕捉瓦加杜古疟疾发病率的年度模式，包括干热月份的低水平持续和雨季的大规模爆发。该模型成功再现了干季约25例/千人的低水平流行，以及雨季约115例/千人的爆发高峰。相比之下，其他模型无法复现这一特征。
进一步的敏感性分析显示，减少季节性迁移的强度不仅能降低干季的稳态发病率，还能减小雨季的爆发规模。计算表明，将迁移量从较高水平降至中等水平，可使年总病例数减少约20%。
在结论与讨论部分，研究明确指出了季节性农村-城市迁移是解释SSA地区，特别是布基纳法索瓦加杜古，疟疾在干季维持低度传播并在雨季引发大规模爆发的关键机制。这一发现填补了现有知识空白，因为多数模型忽视了SSA地区普遍存在的、与农业周期同步的循环性城乡迁移对维持低传播地方性的潜在作用。从政策角度看，这项研究的意义重大。它强有力地表明，疟疾防控策略必须超越静态的、基于固定地点的干预措施。应将人类流动性，特别是季节性劳工迁移的数据和动态，整合到国家疟疾规划中。这包括针对流动人群的监测、在迁移路线上和高风险时段进行精准预防（如发放驱蚊剂、在迁徙前后提供预防性治疗），以及加强区域间的协调。通过理解并应对由迁移驱动的传播机制，公共卫生决策者可以更有效地预测疫情风险，将资源精准投向脆弱社区，从而巩固在低传播区取得的成果，并向消除疟疾的长期目标迈出更坚实的一步。尽管研究存在一些局限性，如未考虑无症状携带者、年龄结构以及气候与迁移的交互影响等，但它提供了一个强有力的机制框架，为量化评估SSA地区针对人口流动的疟疾控制政策和策略效果奠定了重要基础。
