《International Journal of General Medicine》：Development and Internal Validation of a Nomogram for Predicting Severe Complications After Snakebite: A Retrospective Cohort Study
编辑推荐：
背景 蛇咬伤中毒（snakebite envenomation）是全球公共卫生问题，死亡率和发病率高。鉴于临床进展高度异质性以及重症病例存在快速恶化可能，早期风险分层对改善结局至关重要。本研究旨在开发基于简单实验室指标的列线图（nomogram），协助基层医疗
背景 蛇咬伤中毒（snakebite envenomation）是全球公共卫生问题，死亡率和发病率高。鉴于临床进展高度异质性以及重症病例存在快速恶化可能，早期风险分层对改善结局至关重要。本研究旨在开发基于简单实验室指标的列线图（nomogram），协助基层医疗临床医生快速识别高风险患者，从而优化干预并降低死亡率和医疗负担。
方法 研究人员对2014年3月至2025年10月陕西省人民医院医学检验中心收治的10,505例蛇咬伤中毒病例进行回顾性分析。应用严格纳排标准后，2762例患者按3:1比例随机分配至训练队列（n=2072）和验证队列（n=690）。采用LASSO回归（least absolute shrinkage and selection operator regression, LASSO）和多因素logistic回归（multivariable logistic regression）识别严重并发症的独立危险因素，并将其纳入列线图（nomogram）。通过受试者工作特征曲线下面积（area under the receiver operating characteristic curve, AUC）、校准曲线（calibration plots）和决策曲线分析（decision curve analysis, DCA）评估模型性能。
结果 2762例患者中，64例（2.32%）发生严重并发症，最常见为急性心肌损伤、急性肝损伤、凝血病（coagulopathy）和多器官功能障碍；多数表现为两种或以上。咬伤至就诊时间（OR=1.017，95% CI 1.005–1.027）、白细胞计数（white blood cell count, WBC；OR=1.116，95% CI 1.043–1.190）和凝血酶原时间（prothrombin time, PT；OR=1.103，95% CI 1.011–1.203）被识别为独立预测因子。列线图显示良好区分度，训练队列和验证队列AUC分别为0.812（95% CI 0.741–0.882）和0.839（95% CI 0.734–0.945）。校准良好（Hosmer–Lemeshow P>0.05），决策曲线分析证实临床净获益。
结论 该易于获取的列线图基于就诊时间、WBC和PT，显示出较强预测性能、高可操作性和临床可及性，有助于早期识别高风险患者并优化治疗决策。
关键词：蛇咬伤（snakebite）；白细胞计数（white blood cell count, WBC）；凝血酶原时间（prothrombin time, PT）；列线图（nomogram）；预测模型（predictive model）
研究背景与问题：蛇咬伤中毒（snakebite envenomation）是全球重要但常被忽视的公共卫生问题，具有高死亡率和发病率。全球每年约发生540万例蛇咬伤事件，其中中毒病例约180万至270万例，年死亡人数估计为81,000至138,000例；世界卫生组织（WHO）于2019年提出到2030年将蛇咬伤相关死亡和残疾减半的全球战略目标。蛇毒为复杂生化混合物，蛋白约占90%，主要成分包括蛇毒磷脂酶A2（snake venom phospholipase A2, SVPLA2）、三指毒素（three-finger toxins, 3FTx）、蛇毒金属蛋白酶（snake venom metalloproteinases, SVMP）、蛇毒丝氨酸蛋白酶（snake venom serine proteinases, SVSP）、L-氨基酸氧化酶（L-amino acid oxidases, LAAO）和富含半胱氨酸分泌蛋白（cysteine-rich secretory proteins, CRiSP）等。临床预后高度异质，多数患者仅有局部红肿、疼痛或瘀斑，但部分重症可快速出现过敏性休克、严重凝血病、血小板减少和急性肾损伤（acute kidney injury, AKI），甚至多器官衰竭和死亡。尽管已有若干蛇咬伤风险预测工具，但其适用性有限，难以直接转移至所有基层机构。例如，Gopalakrishnan et al和Paul et al均以20分钟全血凝血试验（20-minute whole blood clotting test, 20WBCT）为基础预测死亡和AKI；Cha et al依赖专用血液分析仪计算的未成熟粒细胞参数预测弥散性血管内凝血；Hamman et al需要蛇种、咬伤至就诊时间和抗蛇毒血清剂量等完整记录构建机器学习模型。中国西北基层机构常缺乏蛇种鉴定、专门凝血检测和部分专科实验室检测，且具体并发症在就诊后难以立即确定，因此基层需要整体严重风险评估而非单一结局预测，以便早期分诊、监测和转诊。
研究设计与总体结论：研究人员基于中国西北地区陕西省人民医院医学检验中心2014年3月至2025年10月收治的10,505例蛇咬伤患者资料开展回顾性队列分析，严格纳排后纳入2762例，按3:1随机分为训练队列（n=2072）和验证队列（n=690）。采用LASSO回归（least absolute shrinkage and selection operator regression, LASSO）和多因素logistic回归（multivariable logistic regression）筛选独立危险因素，并构建列线图（nomogram）。研究得出就诊时间、WBC和PT为独立预测因子，列线图在训练和验证队列中具有良好区分度、校准和临床净获益；其重要意义在于仅使用常规急诊实验室变量，具有较高临床可及性，可支持早期风险分层、优化资源配置并辅助基层患者管理，但模型在基层的实际效果仍需前瞻性外部验证。
关键技术方法（不超过250字）：研究为单中心回顾性队列研究，样本来源于陕西省人民医院医学检验中心2014年3月至2025年10月收治的中国西北地区蛇咬伤患者，初始10,505例，严格纳排后2762例；按3:1随机分为训练队列和验证队列（固定随机种子seed=123）。关键技术包括LASSO回归筛选候选指标，多因素logistic回归识别独立危险因素并建模，基于模型系数构建列线图；采用受试者工作特征曲线（receiver operating characteristic curve, ROC）及曲线下面积（area under the curve, AUC）评价区分度，校准曲线和Hosmer–Lemeshow检验评价校准，决策曲线分析（decision curve analysis, DCA）评价临床效用；缺失值>15%的变量排除，其余采用MICE多重插补（5套、20次迭代、seed=123），统计分析使用R 4.2.3，模型开发和报告遵循TRIPOD声明。
研究结果：
（一）人口学与临床特征（Demographic and Clinical Characteristics）。通过纳入2762例患者并比较训练与验证队列基线资料得出，患者中位年龄为53.00（39.00, 61.00）岁，男性1602例（58.00%）；64例（2.32%）发生严重并发症，主要表现为急性心肌损伤、急性肝损伤、凝血病和多器官功能障碍，多数患者发生两种或以上严重并发症；训练队列45例发生严重并发症（2.17%），验证队列19例（2.75%）；两组基线数据无显著差异（P>0.05）。
（二）独立危险因素与模型预测因子筛选（Screening of Independent Risk Factors and Model Predictors）。通过LASSO回归得出就诊时间（Time）、WBC、中性粒细胞/淋巴细胞比值（neutrophil/lymphocyte ratio, NLR）、血尿素氮（blood urea nitrogen, BUN）、胱抑素C（cystatin C, Cyst-C）、钠离子（Na）、PT和凝血酶时间（thrombin time, TT）与严重并发症相关；进一步通过多因素logistic回归得出仅Time（OR=1.017，95% CI 1.005–1.027）、WBC（OR=1.116，95% CI 1.043–1.190）和PT（OR=1.103，95% CI 1.011–1.203）为独立危险因素；共线性诊断显示三者无显著共线性，均作为候选预测因子纳入模型。
（三）列线图预测模型构建（Construction of Nomogram Prediction Model）。通过基于三个独立预测因子及其系数构建可视化列线图得出，严重并发症风险按三步计算：从每个预测因子轴向“Points”轴作垂线获得各预测因子得分；将各预测因子得分相加获得总分；从“Total Score”轴向“Risk”轴作垂线确定严重并发症概率。原文示例：咬伤后24 h就诊、WBC为22.22、PT为15.4的患者总分为33.2，对应概率为15.2%。
（四）列线图验证（Validation of Nomogram）。通过校准曲线、Hosmer–Lemeshow检验、ROC/AUC和DCA在训练及验证队列中评价模型得出，校准曲线接近对角线，Hosmer–Lemeshow检验P值训练队列为0.633、验证队列为0.376（均>0.05）；AUC训练队列为0.812（0.741–0.882）、验证队列为0.839（0.734–0.945）；敏感度和特异度在训练队列分别为68.89%和83.52%，在验证队列分别为78.95%和78.09%；DCA显示列线图提供显著临床净获益。
讨论部分总结：作者讨论确认就诊时间、WBC和PT为蛇咬伤严重并发症的独立危险因素，模型区分度和校准良好，DCA证实临床净获益；与依赖特定生物标志物或复杂影像的模型不同，本模型纳入因素均为常规急诊检测项目，且预测复合终点更符合基层整体风险评估需求，具有较好可操作性和资源有限地区适用性。机制上，就诊时间与毒液作用及药代动力学相关：微循环损伤随中毒时间延长而加重，SVMP可特异性降解基底膜IV型胶原等关键成分，导致毛细血管壁结构完整性丧失、微血管渗漏和出血；凝血功能紊乱随时间累积，毒素持续酶解纤维蛋白原和凝血因子，削弱纤维蛋白凝块稳定性；微循环和凝血系统破坏常先于临床症状，延迟治疗可使潜在损伤演变为不可逆器官功能障碍，淋巴管滞留毒素可持续进入循环，6 h内干预可显著改善凝血功能恢复。WBC反映蛇毒诱导的病理炎症反应，SVPLA2和SVMP为核心毒性成分，可绕过经典免疫识别直接水解细胞膜磷脂，产生高水平花生四烯酸及其代谢物（如前列腺素和白三烯），快速激活免疫细胞；固有免疫细胞释放IL-6、TNF-α和趋化因子，形成正反馈，进一步招募和激活白细胞迁移至损伤组织；WBC升高常提示广泛组织坏死和细胞毒性，Naja annulifera研究显示SVPLA2介导的细胞毒性可导致严重肌肉坏死，并伴中性粒细胞和单核细胞显著增加。PT较昂贵血栓弹力图和缺乏定量标准的20WBCT具有标准化检测和临床可及性优势，PT延长是毒液诱导消耗性凝血病（venom-induced consumptive coagulopathy, VICC）的直接后果，涉及凝血因子过度消耗和功能抑制双机制：抗凝毒素可直接抑制凝血因子活性或产生缺乏XIII因子交联的不稳定纤维蛋白原，导致纤维蛋白聚合失败和抗凝效应；促凝毒素激活凝血因子X和凝血酶原，虽生成纤维蛋白但消耗凝血因子，导致多种凝血因子缺乏和自发性出血高风险，因此PT延长是早期识别出血风险和严重并发症的关键指标。BUN和Cyst-C在单因素分析中与严重并发症相关，但最终多因素模型无独立预测价值，这不否定AKI在预后中的意义，而提示肾功能指标变化存在滞后；蛇咬伤相关AKI涉及SVPLA2和SVMP直接水解肾小球基底膜、损伤肾小管上皮细胞的直接肾毒性，以及通过横纹肌溶解、血栓性微血管病或血流动力学改变导致肾低灌注和缺血性坏死，作为继发性病理改变，AKI指标升高具有明显时间窗。作者团队此前基于WBC、全身炎症反应指数（systemic inflammation response index, SIRI）和BUN在607例Hymenoptera（膜翅目）蜇伤队列开发列线图；本研究将建模方法扩展至蛇咬伤，后者病因不同，以消耗性凝血病为主而非过敏和肾损伤，样本更大（n=2762）且预测因子不同，反映蛇毒时间依赖、以凝血为中心的病理生理机制；两项研究支持列线图风险分层跨不同毒动物损伤的一般性，并强调病因特异性预测因子选择。训练和验证队列阳性预测值分别为8.49%和9.26%，阴性预测值分别为99.18%和99.24%；低阳性预测值主要反映严重并发症低发生率，列线图应主要用于识别需加强监测或转诊的高风险患者，而非作为严重并发症的确定性诊断试验。局限性包括：单中心回顾性分析，样本限于中国西北，未按不同蛇种进行亚组分析，可能引入选择偏倚并限制模型在其他地区的适用性；病历完整性限制导致院前急救措施和某些实验室指标未纳入分析，可能影响模型准确性；临床实施需确保三个预测因子按标准方案测量，任一预测因子不可用或质量差时应谨慎应用，怀疑质量问题时建议重复实验室检测而非插补缺失值，因为该模型的插补策略未在实时临床环境中验证；PT参与凝血病定义并作为预测因子，PT与复合结局的关联可能涉及部分结局泄漏，但单纯PT延长不足以定义凝血病，且预测因子和结局基于不同时间窗，降低了循环性。
研究结论翻译：本研究开发并内部验证了预测蛇咬伤后严重并发症的列线图。模型在训练和验证队列中显示良好区分度、校准和临床净获益。由于仅使用三个常规急诊实验室变量，模型具有较高临床可及性，并具有支持早期风险分层、优化资源配置和辅助基层医疗环境中患者管理的潜力。研究结论及模型在基层医疗环境中的实际效果需在未来研究中开展前瞻性外部验证。
