评估热生物学知情的基本再生数(R0)解释欧洲西尼罗热发病模式的能力

《One Health》：Evaluating the ability of a thermal biology-informed reproduction number to explain patterns of West Nile incidence in Europe

　　

摘要：西尼罗病毒（West Nile virus, WNV）是欧洲公共卫生与兽医健康日益增长的威胁，近几十年来随着气候快速变暖，欧洲WNV暴发日趋频繁。基于实验室蚊媒及病原体实验参数化的热生物学模型（thermal biology model），在因果理解温度对WNV传播的影响方面发挥核心作用。研究人员评估了由尖音库蚊（Culex pipiens）传播的WNV相对基本再生数（relative basic reproduction number, R0rel）的热生物学模型解释2010—2023年欧洲人类西尼罗神经侵袭性疾病（West Nile neuroinvasive disease, WNND）真实世界传播风险模式的能力。研究人员评估了R0rel估算值与WNND观测值的时空吻合度，基于不同分辨率温度数据计算R0rel并评估其相比单独使用温度在WNND风险排序中的价值。此外，研究人员采用广义加性模型（generalized additive model, GAM）检验欧洲WNND发病率对温度的温度响应是否镜像R0的温度响应（含传播最适温度）。研究发现R0rel能准确捕捉WNND的季节动态、峰值风险对应温度，且在WNND风险排序上较单独使用温度略有改善；但R0rel对WNND不规则的年际变化及报告病例较局限的地理分布解释力较差。此外，基于高时空分辨率温度数据经非线性平均计算R0rel在高温条件下略改善风险排序，而基于平均温度计算的估算值在包含低温时表现更优。敏感性分析表明最适传播温度的验证主要由希腊观测值驱动，其他国家在此高温区间信息贡献有限。研究结果表明热生物学模型可捕捉欧洲WNND风险的重要方面，但需纳入额外生态学与流行病学驱动因子以完善风险预测。

论文解读——评估热生物学知情的基本再生数(R₀)解释欧洲西尼罗热发病模式的能力（发表于《One Health》）

一、研究背景与意义

西尼罗病毒（West Nile virus, WNV）是目前欧洲导致人类蚊媒感染的首要病原，通过尖音库蚊（Culex pipiens，欧洲主要媒介）在鸟类储存宿主与蚊虫间维持传播循环，偶发溢出感染人类等终末宿主，可致西尼罗神经侵袭性疾病（West Nile neuroinvasive disease, WNND），病死率约10%。近10—20年欧洲WNV流行范围扩大，部分归因于气候变暖使温度升高更适宜病毒传播。温度通过影响蚊媒生活史参数、媒介传毒能力（vector competence）及外潜伏期（extrinsic incubation period, EIP）等非线性的蚊—病原体性状调控WNV传播，这类关系可通过实验室实验参数化为热驱动的机制模型，并简化为相对基本再生数（relative basic reproduction number, R₀^rel）作为热适宜性（thermal suitability）指标。已有实验室研究给出尖音库蚊传播WNV的R₀^rel最适温度约24.5℃。然而，此类实验室推导的热生物学指数在欧洲人群中经独立现场数据验证仍较少，且欧洲与美国在媒介/宿主种类、病毒株、气候上存在差异，最适传播温度是否吻合尚不清楚；此外，用高时空分辨率还是聚合后温度计算R₀^rel更有价值亦缺乏评估。为此，Heidecke J等研究人员利用2010—2023年欧洲人类WNND月度监测数据与ERA5-Land再分析温度数据，系统验证R₀^rel模型对欧洲WNND时空模式、风险排序及单峰温度响应的解释能力，并比较不同温度聚合方式下R₀^rel计算方法的优劣。

二、主要技术方法概述

研究人员获取欧洲疾控中心（ECDC）TESSy系统中2010—2023年欧盟/欧洲经济区确诊WNND病例，按NUTS3（欧盟三级行政区划）聚合为月度病例数与发病率（Eurostat人口数据校正），排除无症状及非特异性临床表现病例。温度数据取自ERA5-Land逐日2米气温（0.1°分辨率，2010—2023），经人口加权空间平均及1—4个月滑动平均（最终选用3个月滑动平均）匹配NUTS3尺度。R₀^rel依据已发表并微调的尖音库蚊—WNV温度响应函数计算，取值归一化至[0,1]。为比较计算方式影响，分别在时空聚合前后计算R₀^rel（逐日/月网格→聚合 vs 聚合温度→R₀^rel）。采用Kendall's τ_B秩相关评估温度与R₀^rel对WNND发病率的风险排序能力，置换检验比较差异显著性。采用含负二项分布的广义加性模型（generalized additive model, GAM），log(人口)为偏移量，惩罚立方回归样条拟合温度/R₀^rel与WNND的非线性关联（共享平滑+国家特异随机平滑），估计经验最适温度并与实验室R₀^rel最适温比对。部分GAM引入CORINE土地覆被（城市/耕地比例）、3个月滑动累积降水、人均GDP及≥65岁人口占比作校正。开展剔除希腊数据的敏感性分析以评估高温端信息驱动来源。

三、研究结果

3.1. Summary of data（数据概况）

研究期共报告3844例人类WNND，>96%发生于7—9月，>92%来自希腊、意大利、罗马尼亚和匈牙利；热点区域（连续≥6年有病例的相邻NUTS3）见于希腊北部、意大利北部、罗马尼亚南部及匈牙利东部。有病例年份中7—9月平均人口加权3个月滑动均温为22.3℃（范围14.1—28.7℃），仅17.5%超过实验室R₀^rel最适温24.4℃。

3.2. Spatial-temporal comparison（时空模式比较）

空间上，7—9月多年平均R₀^rel高值区与WNND累计高发区大致吻合，但相当多具相似热适宜性的地区未报告病例，表明R₀^rel高估了实际发病地理范围。时间上，热点区域R₀^rel与标准化月度WNND在时间窗口及季节峰值一致，但R₀^rel无法解释WNND强烈的年际波动（无病例年份与有病例年份R₀^rel可相近）。

3.3. Rank correlation and aggregation methods（秩相关与聚合方法比较）

全数据集：R₀^rel（先温度时空聚合再算R₀^rel，标记为R_{0,NUTS3,window}^rel）与月度WNND发病率Kendall τ_B=0.440，略高于温度的0.419（p<0.05）；先算R₀^rel再空间聚合（网格级先算）反而降低秩相关。限定温度>15℃子集：先逐日算R₀^rel再聚合（R_{0,NUTS3,daily}^rel）τ_B=0.413为最优但不显著优于温度（0.408）。说明整体而言经聚合温度计算的R₀^rel略优，高温子集高时分辨率先算R₀^rel略好，符合热限值附近平均温度输入更适、最适温附近非线性高分辨先算更适的理论预期。

3.4. GAM analysis（广义加性模型分析）

GAM显示WNND对3个月滑动均温呈单峰（unimodal）响应，经验估计最适温度为24.9℃（95% CI: 24.5—25.4℃），与实验室R₀^rel最适温24.4℃（95% CI: 22.9—25.9℃）高度接近；此单峰形态主要由希腊数据支撑（希腊最适温25.1℃, 95% CI: 24.7—25.4℃），他国高温观测少致不确定。加入土地覆被、降水、GDP及老龄人口比校正后最适温微升至25.1℃（95% CI: 24.5—26.1℃），希腊调至25.9℃（95% CI: 25.2—27.0℃），高温段置信区间变宽。WNND对R₀^rel呈单调但非线性递增关系。温度作自变量的GAM的AIC略优于R₀^rel作自变量者。

3.5. Sensitivity analysis（敏感性分析）

剔除希腊数据后，温度—WNND单峰形态不确定性大增难以区分单峰/单调；R₀^rel—WNND仍呈单调非线性关系，R₀^rel秩相关仍优于温度；高温子集最优算法变为先网格级算R₀^rel再空间聚合略优，提示希腊复杂海岸地形致ERA5-Land网格精度差异影响原先结果。

四、讨论与结论总结（翻译浓缩）

机制性热生物学模型可为媒介传播疾病热适宜性评估提供工具。本研究表明实验室推导的R₀^rel略改善欧洲WNND风险排序能力，且其预测的最适传播温度与欧洲WNND经验峰值温度高度吻合；3个月滑动平均R₀^rel较好复现WNND季节性与峰值时间。但R₀^rel未能充分解释WNND年际变异及较窄的实际发病地理分布，反映当前模型未考虑宿主免疫、非温气候因子、土地利用及候鸟储存宿主动态等。先时空聚合温度再算R₀^rel整体更优，高温段高时分辨率先算R₀^rel略好，符合疟疾蚊相关研究结论；网格先算R₀^rel再空间平均效果略差可能与再分析温度沿海精度有关。GAM揭示欧洲WNND—温度呈单峰响应且最适温近实验室值，此前欧洲多项研究因设计差异仅报道正向关联，未来暖化使高温超限抑制传播效应将更易检测，需在预测中纳入单峰温度效应。敏感性分析提示最适温验证强烈依赖希腊主动监测数据，他国高温观测不足限制跨国异质性探讨。局限性含各国监测差异、用人病例代表鸟—蚊放大循环存在偏差、统计关联非因果推断。

结论（Conclusion）： 实验室推导的热生物学驱动模型相较单独使用温度可改善欧洲WNND风险排序，并能捕捉观测动态的若干特征（包括季节性及近似峰值风险温度）。然而仅靠热适宜性无法解释WNND大量变异，支持将热生物学模型用作WNV传播潜力温度相关变化的机制指示指标，同时强调需整合或扩展模型以纳入额外驱动因子并修正简化假设，从而实现更全面的风险评估。