《Human Vaccines & Immunotherapeutics》：Global burden of varicella in children aged 5–9?y, relationship with vaccination policy indicators, and projected trends to 2030
编辑推荐：
水痘（varicella）仍是一种全球常见的儿童感染。厘清长期趋势并识别疾病负担的关键决定因素，对制定疫苗接种策略与公共卫生规划至关重要。研究人员从全球疾病负担（GBD）2021研究中获取数据，并结合全球水痘疫苗接种政策信息，训练极端梯度提升（Extreme
水痘（varicella）仍是一种全球常见的儿童感染。厘清长期趋势并识别疾病负担的关键决定因素，对制定疫苗接种策略与公共卫生规划至关重要。研究人员从全球疾病负担（GBD）2021研究中获取数据，并结合全球水痘疫苗接种政策信息，训练极端梯度提升（Extreme Gradient Boosting, XGBoost）模型预测2030年疾病负担，估算全球5–9岁儿童水痘的发病率、患病率及伤残调整生命年（disability-adjusted life years, DALYs）。采用Shapley Additive exPlanations（SHAP）识别关键因素以指导公共卫生策略。1990至2021年，全球发病率由1028.46/10万升至1048.24/10万，患病率由30.81/10万升至31.07/10万，DALYs由21.09/10万降至12.36/10万；至2030年预计分别降至897.01/10万、28.96/10万和11.71/10万。1990至2030年间，东南亚与南亚发病率最高，亚洲及太平洋岛国患病率最高，撒哈拉以南非洲与南亚DALYs最高。SHAP显示人口是最强预测因子（发病率平均绝对SHAP值0.167，患病率0.103），较高疫苗接种项目层级（vaccination program tier）可能与较低DALYs相关。5–9岁儿童水痘全球负担总体下降，但模型预测区域差异将持续存在。缺乏疫苗接种覆盖率与队列免疫数据限制了疫苗效应解读；未来研究应纳入覆盖率与队列暴露信息以强化政策评估。
研究背景与立项依据
水痘由水痘-带状疱疹病毒（varicella-zoster virus, VZV）引起，主要通过呼吸道飞沫与直接接触传播，在儿童中虽多自限但造成显著全球健康损失。世界卫生组织（World Health Organization, WHO）估计全球每年约1.4亿例发病，其中约420万例出现严重并发症。既往年龄-时期-队列（age-period-cohort, APC）、贝叶斯年龄-时期-队列（Bayesian APC, BAPC）与自回归积分滑动平均（autoregressive integrated moving average, ARIMA）等传统统计模型在应对疫苗接种等非线性外生干预时存在结构断点响应差、难以融合多源数据与复杂交互等局限。因此，研究人员基于GBD 2021数据，采用XGBoost联合SHAP框架，分析1990–2021年5–9岁儿童水痘负担并预测至2030年，以期为靶向公共卫生策略提供证据。
主要技术方法概述
研究人员以GBD 2021结果工具获取1990–2021年204个国家与地区的5–9岁儿童水痘发病率、患病率与DALYs（伤残调整生命年）率，人口数据来自健康指标与评估研究所（Institute for Health Metrics and Evaluation, IHME）；疫苗接种政策来自Lee等综述与WHO免疫数据，构建四级分类变量vaccination_program_tier（0无程序、1有限接种、2全国一剂次、3全国两剂次）。建模以性别、年份、接种层级与5–9岁分性别人口对数作特征，对比弹性网络正则化线性回归与XGBoost，通过滚动起源验证（三扩展窗口+2020–2021独立测试集）选优，XGBoost测试集RMSE（发病率354.06/10万、患病率7.38/10万、DALYs 10.21/10万）与Pearson相关（0.5253、0.5399、0.5589）均优于线性回归，故用于预测。预测期假设各国层级不变，用IHME人口投影与500次bootstrap给出95%不确定区间（uncertainty interval, UI）。SHAP用于全局特征重要性与依赖方向解释。敏感性分析在五类社会人口学指数（socio-demographic index, SDI）层与基线/扩张/保守三种政策情景下重投影。
研究结果
模型性能：XGBoost在独立测试集上三项指标RMSE与相关均优于正则化线性回归，被选为最优预测模型。
全球与SDI级趋势（1990–2021）：全球发病率由1028.47/10万微升至1048.25/10万，患病率由30.81/10万微升至31.07/10万，DALYs由21.09/10万降至12.35/10万。多数SDI区呈倒U型，低SDI区DALY最高且降幅最小（仅降约35%），高SDI区患病率峰值明显。
至2030年全球预测：发病率率降至897.01/10万（95%UI 683.17–1170.62），患病率28.96/10万（22.31–37.58），DALY率11.71/10万（8.63–15.82）；发病数降至592万，患病数19万，DALY数继续微降。
国家层面异质性（1990–2030）：1990与2021年菲律宾发病率最高（约1962–1966/10万），马尔代夫、泰国、塞舌尔患病率居前，中非共和国DALY最高（53.90→56.01/10万）；2030年圣文森特和格林纳丁斯发病率最高（1355.71/10万），马尔代夫患病率最高（38.95/10万），尼日利亚DALY最高（23.21/10万）。低SDI的撒哈拉以南非洲与南亚DALY持续偏高。
SHAP变量解析：对发病率与患病率，log(5–9岁人口)平均绝对SHAP最高（0.167与0.103），男性风险权重更高，年份对发病率影响渐增、对DALY渐减；对DALYs，vaccination_program_tier为第二重要（0.121），较高层级对应较低预测DALY，但因未纳入覆盖率，层级仅为情境代理而非因果。
讨论部分总结
研究人员指出，男性在模型中获更高发病/患病权重，与部分流行病学报道一致，可能受流动人群性别比、职业暴露与就医延迟影响。区域差异方面，低SDI持续高DALY与其他传染病（如结核、狂犬病、脑型疟疾）分布一致。SHAP中人口为主驱动，接种层级与低DALY相关但属关联非因果——高SDI国同时具备高层级与高监测灵敏度，使发病/患病记录偏高，形成表面“反直觉”。研究局限包括GBD估算依赖异质源、接种变量仅为政策采纳层级而非实现覆盖率、SHAP仅反映关联。结论部分表明：1990–2021年全球5–9岁儿童水痘发病/患病略升、DALY持续下降；2030年预测略降；人口为最强预测因子，高接种层级与低DALY预测值相关但无覆盖数据不能立因果；XGBoost+SHAP组合对传染病预测与解释有价值，未来需补覆盖率与准实验设计以评政策真实效应。论文发表于《Human Vaccines》。
