《Forensic Science International》：Machine Learning and Metabolic Signatures of Drowning: A Pathway to Uncovering Cause of Death in Aquatic Environments
编辑推荐：
本研究针对溺亡尸体解剖诊断困难、特异性体征缺乏的难题，创新性地采用死后代谢组学分析，通过构建基于股血样品的多变量模型，成功将溺亡与慢性心脏病、中毒、创伤及缢死等主要鉴别诊断区分开来，模型展现出高灵敏度（83-87%）与特异性（78-89%）。该研究为法医科学提供了潜在的客观生化诊断工具，有望提升水环境中死亡案件的死因判定准确性。
当一具尸体在水中被发现，法医病理学家面临的挑战往往超乎想象。是单纯的溺亡，还是死者因心脏病突发、药物中毒或遭遇创伤后落水？传统上，法医依赖尸检发现来推断死因，但对于溺亡，这些迹象却充满了不确定性。典型的溺亡体征，如水性肺气肿（emphysema aquosum）、气道泡沫和肺水肿，虽能提供线索，却并非特异——它们同样可能出现在其他死因中，如中毒或某些自然疾病。更棘手的是，这些体征可能在尸体浸泡数天后消失，或因死后浸泡本身而产生，使得仅凭尸检做出明确诊断变得异常困难。这种诊断的不确定性，可能导致司法误判，使真相沉入水底。因此，法医学界亟需一种更客观、更可靠的辅助诊断方法。
在此背景下，一项发表于《法医科学国际》（Forensic Science International）的研究带来了新的希望。由阿尔伯特·埃尔姆舍（Albert Elmsj?）等来自瑞典国家法医学委员会法医遗传与毒理学系的研究团队，将目光投向了代谢组学（metabolomics）——这门全面研究生物系统内小分子代谢物的学科。他们设想，溺亡这一独特的窒息过程，是否会在人体内留下特定的“代谢指纹”？这些微观的化学痕迹，或许能成为区分溺亡与其他死因的关键。为了验证这一假设，研究人员开展了一项大规模的回顾性研究，旨在探究死后代谢组学分析能否有效区分溺亡与其他竞争性死因。
研究者利用了瑞典国家法医学委员会在2017年11月至2020年10月间进行法医尸检的存档样本。研究共纳入503例溺亡案例，并设置了四个对照组，分别代表在水体中发现尸体时可能需考虑的其他死因：慢性心脏病（CHD，510例）、中毒（516例）、创伤（497例）。此外，还纳入了缢死（511例）作为第四个“阳性”对照组，以测试模型区分两种不同窒息过程的能力。所有对照组案例均与溺亡案例在性别、年龄、身体质量指数（BMI）和死后间隔（PMI）上进行了匹配。核心技术方法包括：使用标准化流程对尸检时采集的股血样本进行蛋白沉淀前处理；采用超高效液相色谱-电喷雾-四极杆飞行时间质谱（UHPLC-ESI-QToF）系统进行代谢物数据采集；利用R语言中的“XCMS”和“CAMERA”软件包进行代谢组学数据预处理和特征注释；最后，使用SIMCA软件进行多变量建模，主要采用正交偏最小二乘判别分析（OPLS-DA）来识别区分溺亡与各对照组的代谢特征变量，并通过训练集与验证集的分割来评估模型的预测性能。
3.1. 研究人群的人口统计学特征
溺亡组与四个对照组在性别和死后间隔（PMI）上无显著差异。然而，创伤组和慢性心脏病（CHD）组在年龄分布上与溺亡组存在显著差异，这主要是由于低分位数不同所致。缢死组的身体质量指数（BMI）则显著低于CHD组和创伤组。尽管存在这些统计学上的差异，但由于样本量较大（每组约500例），其对新陈代谢组的实际影响尚不确定。
3.2. 溺亡事件的情况
对尸检和警方报告的回顾显示，23%的溺亡案例有目击者，54%有明确的溺亡迹象（如落水前正在进行游泳或钓鱼等活动）。大多数溺亡发生在淡水环境（63%），其次是微咸水（19%）、自来水（11%）、海水（6%）和游泳池水（1%）。36%的案例接受了心肺复苏（CPR），其中31%未后续住院，5%后续住院。
3.3. 尸检结果
在溺亡案例中，水性肺气肿的出现率为80%，气道泡沫为35%，肺水肿为61%，胸腔积液为39%。这些发现凸显了传统尸检体征的变异性与非特异性。
3.4. 多变量建模与模型评估
质谱数据处理后，保留了2610个特征用于建模。首先，将所有五组（溺亡、缢死、CHD、中毒、创伤）纳入一个OPLS-DA模型，结果显示各组间存在重叠，模型对训练集的总体分类准确率为65%，对溺亡类的灵敏度为77%，特异性为91%。在验证集中，模型对溺亡的预测灵敏度为69%，特异性为84%。
鉴于法医实践中很少同时考虑五种死因假设，为了更贴近实际场景，研究团队构建了四个独立的二元OPLS-DA模型，分别将溺亡与一个对照组进行比较。这些二元模型的性能显著提升，R2值介于0.61-0.76之间，Q2值介于0.40-0.56之间。溺亡案例的分类灵敏度在83%至87%之间，特异性在78%至89%之间，表现出优异的区分能力。
为了探究水体浸泡本身（而非溺亡过程）对新陈代谢组的影响，研究者还建立了仅包含短死后浸泡间隔（小于24小时）溺亡案例的模型，并与包含各种浸泡间隔案例的模型进行比较。结果发现，仅使用短浸泡间隔案例的模型性能仅有微小差异。此外，基于不同“捞出后间隔”（即尸体从水中移出到尸检的时间）的模型比较也显示，短间隔与长间隔案例的模型之间没有明显差异，且用短间隔案例训练的模型能够预测长间隔案例，反之亦然。这表明死后浸泡时间和捞出后间隔对新陈代谢组分类模型的影响微乎其微。
3.5. 代谢物鉴定与通路分析
通过数据库匹配，研究鉴定出52个能区分溺亡组与四个对照组的特征代谢物。代谢物分析显示，不同死因之间存在 distinct 的模式。例如，1-(β-D-核糖呋喃基)-1,4-二氢烟酰胺、5-羟基-L-色氨酸、甲基鸟苷和二甲基鸟苷等在CHD和中毒案例中的值显著高于溺亡和缢死案例。溶血磷脂酰胆碱（LysoPCs）在CHD和中毒案例中的值普遍低于缢死、创伤和溺亡案例。肉碱相关代谢物则显示出混合趋势。
通路富集分析揭示了溺亡与各对照组之间差异显著的代谢通路。在至少一种二元比较中，以下通路显示出联合p值低于0.05：C21-类固醇激素生物合成与代谢、脂肪酸活化、糖酵解与糖异生、尿素循环/氨基代谢、辅酶A（CoA）分解代谢、维生素B5 - 从泛酸合成CoA、药物代谢 - 细胞色素P450、甘油磷脂代谢以及生物蝶呤。
该研究的结论与讨论部分深入剖析了其发现的意义与局限。首先，研究证实了死后代谢组学在区分溺亡与其他常见死因方面具有很高的准确度，其二元模型表现出的高灵敏度和特异性使其有潜力成为法医实践中有价值的辅助工具。特别是在传统尸检发现模棱两可、需要评估两种竞争性死因假设时（例如，是高浓度药物导致的中毒，还是溺亡伴有药物影响），代谢组学模型可以提供客观的化学证据支持。
其次，研究揭示了不同死因间独特的代谢特征。例如，CHD与中毒案例在代谢谱上显示出相似性，这可能与两者都可能涉及较长的濒死期或特定的代谢紊乱有关；而创伤与缢死案例的代谢谱也较为接近，可能反映了急性机械性损伤与窒息的共同应激反应。溺亡则显示出最为独特的代谢轮廓，这可能与水吸入、缺氧挣扎以及可能的电解质紊乱等特定过程相关。通路分析进一步指出，甘油磷脂代谢、类固醇激素生物合成和细胞色素P450相关药物代谢等通路在区分溺亡中扮演重要角色。
然而，研究也存在局限性。最主要的假阴性案例多发生于死亡前曾住院治疗的溺亡者，这些案例常被误分类为中毒。这可能是由于医疗干预（如药物输注、复苏努力）导致了延长的呼吸衰竭过程，使其代谢谱更接近中毒。因此，在应用此类模型时，可能需要排除或有区别地分析此类经医院治疗的案例。此外，虽然通过亚组分析表明死后浸泡时间对模型影响很小，但所有对照组均未接触水体，理论上完美的对照应是死于其他原因后再被浸入水中的尸体。由于这类案例罕见，本研究未能纳入，这是未来研究需要补足的方向。
尽管存在这些挑战，本研究清晰地展示了代谢组学作为法医病理学新型工具的广阔前景。它能够从常规毒理学筛查已获取的血液样本中挖掘出更深层的信息，通过分析内源性代谢物的模式来预测可能的死因，为法医工作者提供了一个强大的、数据驱动的决策支持手段。这不仅有助于提高溺亡诊断的准确性，其方法论和思路也可扩展到其他复杂死因的鉴别中，推动法医科学向更精准、更客观的方向发展。
