中国新生儿脑膜炎的时空趋势、病原分布与负担预测（1990–2023）：基于2023年全球疾病负担研究的描述性分析

《Brain‐X》：Temporal trends, pathogen distribution, and forecasted burden of neonatal meningitis in China (1990–2023): A descriptive analysis based on Global Burden of Disease 2023

　　

摘要专业翻译 背景 新生儿脑膜炎仍然是导致死亡和健康损失的重要原因。研究人员利用2023年全球疾病负担（Global Burden of Disease, GBD）估算数据，描述了中国新生儿脑膜炎的时间趋势、病原分布及风险归因负担。 方法 研究人员分析了GBD模型估算的早期（0–6天）和晚期（7–27天）新生儿的发病率、死亡率、伤残调整生命年（disability-adjusted life years, DALYs）、损失寿命年（years of life lost, YLLs）和伤残生存年（years lived with disability, YLDs）。采用估算年百分比变化（estimated annual percentage change, EAPC）和分段对数线性回归评估趋势。对数指数平滑（exponential smoothing, ETS）模型预测了截至2035年的死亡率，并以对数自回归积分滑动平均（auto-regressive integrated moving average, ARIMA）模型作为敏感性分析。同时评估了病原构成比和GBD估算的风险归因负担。 结果 1990年至2023年，模型估算的发病率和死亡率分别从每10万人口715.7例和98.0例下降至79.0例和4.5例（EAPC = ?7.6% [95%置信区间{CI} ?8.2至?6.9]和?9.6% [95% CI ?9.9至?9.2]）。早期新生儿死亡率仍高于晚期新生儿死亡率（2023年率比为1.8）。分段回归识别出1997–2002年为下降最陡峭的时期。2035年对数ETS预测值分别为总体、早期和晚期新生儿每10万人口4.3（95%预测区间[PI] 0.5–10.4）、5.5（95% PI 0.6–12.5）和3.9（95% PI 0.5–9.7）例。前五种病原体归因死亡的累计占比从58.9%上升至62.9%。低出生体重和短妊娠期具有最大的归因分数（47.7% [95%不确定性区间44.0%–51.1%]）。 结论 尽管负担显著下降，但持续存在的早期新生儿脆弱性和分歧性的探索性预测强调，需要开展临床监测以验证这些模型发现并为预防策略提供信息。

论文解读文章

一、研究背景与问题提出

新生儿脑膜炎是全球重大公共卫生挑战之一。世界卫生组织（World Health Organization, WHO）发布的《2030年战胜脑膜炎》全球路线图强调降低死亡率、预防残疾并改善幸存者生活质量。新生儿期在生命周期中具有独特的脆弱性：免疫系统不成熟、围产期暴露、生理储备有限以及早期症状隐匿，这些因素共同导致诊断延迟和严重不良结局。即使在死亡率已下降的地区，新生儿脑膜炎仍与神经损伤、听力丧失、发育迟缓、延续至青春期的精神科后遗症以及长期医疗需求密切相关。

中国现有的新生儿脑膜炎临床证据主要来源于医院队列、极早产儿人群和儿童细菌性脑膜炎监测，这些研究虽然提供了详细的临床微生物学和结局数据，但并非连续的全国性监测，其特定场景和人群限制了对全国数十年变化趋势的推断能力。此外，早期和晚期新生儿群体可能呈现不同的负担轨迹，但这一全国性年龄特异模式尚未得到充分描述。一项基于GBD 2021的全球研究虽涵盖0–14岁儿童脑膜炎并包含早晚期新生儿的独立分析，但其广泛的全球设计未能提供聚焦于新生儿期的中国特异性综合评估。

因此，本研究利用GBD 2023更新估算数据，系统整合了负担趋势、年龄特异性比较、病原归因、出生体重（birth weight, BW）和孕周相关归因、分段趋势敏感性分析以及截至2035年的探索性预测，填补了这一研究空白。

二、研究设计与关键方法

研究人员开展了一项基于GBD 2023模型估算的二次描述性生态学研究，遵循健康估算准确透明报告指南（Guidelines for Accurate and Transparent Health Estimates Reporting, GATHER）。数据提取自中国（GBD位置ID=6）、两性合并（性别ID=3）、脑膜炎（病因ID=332）、1990–2023年，覆盖早期新生儿（0–6天；年龄ID=2）、晚期新生儿（7–27天；年龄ID=3）和总体新生儿（年龄<28天；年龄ID=42）年龄组。研究使用公开可获取的聚合模型估算数据，无需伦理审批和知情同意。

关键技术方法包括：其一，采用对数线性回归模型计算估算年百分比变化（EAPC）及其95%置信区间（CI），量化长期负担变化的方向和幅度；其二，运用Python实现的连续分段对数线性回归识别时间转折点，主分析允许最多三个断点、每段至少4年，以贝叶斯信息准则选择最优模型，并对最大断点设置（2、3、4）和最小段长设置（3、4、5年）进行27种组合的敏感性分析；其三，死亡率预测采用自然对数变换后的阻尼趋势指数平滑模型ETS(A, Ad, N)，以自动选择的ARIMA模型作为敏感性分析，通过滚动起点验证（2014–2022年预测起点）评估平均绝对误差、均方根误差、对称平均绝对百分比误差（sMAPE）、95%预测区间（PI）覆盖率和平均PI宽度，主分析95% PI基于5000条残差自助路径生成；其四，病原结构采用排名、个体占比和前3、前5、前10累计占比进行描述性评估。数据分析使用Python 3.14.3，分段回归使用NumPy 2.4.4，预测模型使用statsmodels 0.14.6。

三、研究结果

3.1 中国新生儿脑膜炎的发病率、死亡率和残疾负担

1990年至2023年，GBD估算的新生儿（<28天）脑膜炎发病例数从12,300例（95% UI 8558–17,346）降至558例（95% UI 427–707），发病率从每10万人口715.7例（95% UI 497.9–1009.2）降至79.0例（95% UI 60.5–100.1），降幅达89.0%，EAPC为?7.6%（95% CI ?8.2至?6.9）。死亡率从每10万人口98.0例（95% UI 62.9–144.4）降至4.5例（95% UI 3.0–7.2），降幅达95.4%，EAPC为?9.6%（95% CI ?9.9至?9.2）。DALY和YLL率呈现几乎相同的降幅，YLD率绝对值较小，从每10万人口1.44例（95% UI 0.89–2.23）降至0.23例（95% UI 0.14–0.36），EAPC为?6.5%（95% CI ?7.2至?5.8）。DALY仍以YLL为主导。总体而言，死亡率下降速度快于发病率或YLDs，表明疾病发生、致命负担和模型化非致命健康损失之间的改善并不均衡。

3.2 新生儿脑膜炎负担的年龄分段差异

两个新生儿年龄组的负担均有所下降，但模型化发病率与死亡率的趋同程度不同。早期新生儿发病率从每10万人口1173.2例（95% UI 813.3–1661.0）降至81.6例（95% UI 62.1–104.0），晚期新生儿从562.1例（95% UI 392.0–790.5）降至78.2例（95% UI 59.9–98.8）。早期新生儿死亡率从161.1降至每10万人口6.9例，晚期从76.8降至3.8例。早期与晚期发病率比从1990年的约2.1收窄至2023年的1.0，而死亡率比从约2.1降至1.8。这表明到2023年两个年龄组的发病率已接近趋同，但模型化死亡率仍存在早期新生儿梯度。

3.3 新生儿脑膜炎死亡率的时间模式

分段分析显示各新生儿年龄组的死亡率趋势呈非线性但一致的模式。总体新生儿的死亡率APC在1990–1997年为?5.3%（95% CI ?6.4至?4.2），1997–2002年加速至?17.9%（95% CI ?19.0至?16.7），2002–2011年放缓至?6.4%（95% CI ?7.0至?5.7），2011–2023年再次加速至?10.1%（95% CI ?10.6至?9.6）。在最陡峭的1997–2002年段，早期和晚期新生儿的APC分别为?18.3%（95% CI ?19.3至?17.3）和?17.6%（95% CI ?18.8至?16.4）。在27种年龄组/参数组合中，所有敏感性设置均识别出重叠于1990年代末至2000年代初的陡峭下降区间，表明这一加速是模型化死亡率序列的共同特征。

3.4 截至2035年新生儿脑膜炎死亡率的探索性预测

从2024年到2035年，点预测显示总体新生儿死亡率从4.5降至每10万人口4.3例，早期新生儿从6.7降至5.5例，而晚期新生儿仅从3.8微调至3.9例。2035年预测值分别为总体、早期和晚期新生儿每10万人口4.3（95% PI 0.5–10.4）、5.5（95% PI 0.6–12.5）和3.9（95% PI 0.5–9.7）例。对数ARIMA敏感性预测在总体和早期新生儿组方向相似，但两个模型家族均显示晚期新生儿点轨迹接近平坦或略有上升，伴随长期预测区间增宽。滚动起点sMAPE在模型家族间差异仅为0.028、0.399和0.209个百分点。Jarque–Bera检验拒绝了所有拟合模型的残差正态性，因此预测被视为探索性而非确定性预测。

3.5 归因发病率和归因死亡率的病原谱

归因发病率和归因死亡率的病原排名呈现不同时间模式。归因发病例数中，1990年非脊髓灰质炎肠道病毒排名第一（34.0%），B族链球菌（Group B Streptococcus, GBS）第二（20.2%）；到2023年位置反转，GBS第一（25.7%），非脊髓灰质炎肠道病毒第二（22.7%），随后为其他病毒性脑膜炎病原（6.8%）、肺炎链球菌（Streptococcus pneumoniae）（5.9%）和大肠埃希菌（Escherichia coli）（5.0%）。归因发病率前五累计占比从69.5%降至66.2%。归因死亡方面，GBS在1990年和2023年均排名第一（32.7%和35.7%），肺炎链球菌保持第二（从8.1%升至11.3%），2023年前五还包括金黄色葡萄球菌（Staphylococcus aureus）（6.3%）、其他病毒性脑膜炎病原（5.0%）和念珠菌属（Candida spp.）（4.6%）。归因死亡率前五累计占比从58.9%小幅升至62.9%，前10占比保持稳定（82.4%至81.8%），表明时间变化主要发生在广泛稳定的前10集中度内的重新排序。

3.6 归因于出生体重和孕周类别的新生儿脑膜炎死亡率

2023年，低出生体重和短妊娠期类别的GBD估算归因负担最大：15.3例死亡（95% UI 10.6–23.3），每10万人口2.17例（95% UI 1.50–3.30），人群归因分数（population-attributable fraction, PAF）为47.7%（95% UI 44.0%–51.1%）。低出生体重伴正常孕期为13.8例死亡（95% UI 9.3–20.9），率1.95（95% UI 1.31–2.96），PAF为42.9%（95% UI 39.8%–45.8%）。短妊娠期伴正常出生体重仅0.6例死亡（95% UI 0.3–1.2），率0.08（95% UI 0.04–0.17），PAF为1.8%（95% UI 0.9%–3.1%）。早期和晚期新生儿组均观察到相同的点估计排序，低出生体重和短妊娠期类别点估计最高。这些分类估算未被解释为独立、可加或直接因果的临床效应。

四、讨论与结论

讨论总结

本研究基于GBD 2023模型估算，识别出中国1990–2023年新生儿脑膜炎发病率、死亡率、DALYs、YLLs和YLDs的显著下降，早期新生儿死亡率仍高于晚期，1997–2002年为最陡峭的死亡率下降期且在27种分段回归设置中稳健。病原排名在归因发病率和归因死亡率之间不同，但死亡率前五占比变化有限。BW和孕周类别占GBD估算归因负担的很大比例。截至2035年的预测因长期不确定性宽泛而仍属探索性。

与既往研究相比，本研究提供了聚焦新生儿期的中国特异性设计、GBD 2023更新序列、归因发病率和死亡率分别的病原谱、BW和孕周风险归因、27设置分段趋势敏感性分析以及经诊断验证的探索性预测。YLD率下降慢于死亡率但绝对值低，DALYs仍以YLLs为主，支持幸存者结局需在纵向临床队列中直接研究的假设。早期与晚期死亡率比持续存在可能反映生物学脆弱性、围产期暴露、早产、早期全身不稳定和诊断困难等综合因素，但本研究未检验这些机制。GBD数据无法区分医疗相关感染与社区获得性感染，这一区分对预防策略至关重要。分段回归识别了变化时间但未确定原因，1990年代末至2000年代初的加速与中国妇幼健康和新生儿照护扩展时期重叠，但归因于具体政策需独立因果分析。全球研究中GBS（22.70%）、肺炎克雷伯菌（16.97%）和病毒性病原（15.35%）为2021年新生儿归因死亡率的主要类别，而本研究中国特异性GBD 2023估算中GBS（35.7%）和肺炎链球菌（11.3%）排名前两位。大肠埃希菌在2023年归因死亡率中排名第10，而中国NICU研究和抗菌药物耐药监测常识别出临床重要的革兰阴性感染，这一差异与地理范围、GBD版本、模型输入和归因框架的不同相容，不应被解读为微生物学 prevalence 的直接差异。前五累计占比4.0个百分点的小幅增加和前10占比稳定表明集中度和排名变化有限，金黄色葡萄球菌在2023年模型化位置的上升可能识别监测优先事项，但本研究数据无法确定其来源。

研究结论翻译

GBD 2023模型估算表明，1990年至2023年中国新生儿脑膜炎负担显著下降，同时早期新生儿死亡率仍高于晚期新生儿。1997–2002年最陡峭的下降在分段模型设置中稳健。探索性预测还显示不同年龄段的未来轨迹存在分歧，晚期新生儿死亡率保持近平坦。病原占比变化有限，风险结果代表GBD估算的归因负担，截至2035年的预测仍不确定且属探索性。需要中国直接临床注册、微生物学监测和纵向结局研究来验证模型化模式并指导具体干预措施。