基于大鼠骨骼肌微生物组的死后浸没时间高精度预测模型：一种独立于死因的法医学新方法

《Frontiers in Microbiology》：Postmortem submersion interval prediction model based on the rat muscle microbiome

本研究系统揭示了大鼠骨骼肌微生物群在死后水体浸没（PMSI）过程中呈现规律演替（从变形菌门到厚壁菌门优势转变），并成功构建了不受死因（溺死与死后浸没）影响的二阶段机器学习预测模型（随机森林分类准确率90.9%，早期/晚期回归MAE分别为0.303天与1.293天），为法医学长期PMSI估算提供了基于内部组织的稳定"微生物时钟"新策略。

1 引言

准确估算死后浸没时间（PMSI）是法医学实践的关键环节，对刑事侦查、嫌疑人排查和案件时间线重建具有重要意义。传统方法（如体温、尸僵、尸斑评估）主要适用于早期死后阶段（通常72小时内），且易受环境温度等干扰。法医昆虫学对水下尸体的PMSI估算信息有限。近年来，高通量测序技术推动了法医微生物学发展，"微生物时钟"概念逐渐形成。当前研究多聚焦于易受环境污染的肠道、口腔和皮肤微生物组，而内部组织的潜力尚未充分探索。骨骼肌作为生物量大、相对封闭的内部组织，可能为PMSI估算提供更稳定的信息源。本研究假设骨骼肌微生物群在水浸没后会发生规律演替，且该模式不受死因（溺死vs.死后浸没）影响，旨在通过大鼠模型系统追踪14天内骨骼肌微生物组的动态变化，并评估其构建高精度PMSI预测模型的潜力。

2 材料与方法

2.1 动物模型与实验设计

研究使用72只健康成年雄性SD大鼠，随机分为溺死组（D组）和死后浸没组（PS组）。D组大鼠通过重复浸没-暴露模拟溺死过程；PS组大鼠采用CO₂吸入法处死。两组尸体均置于自然水体环境中，分别于0、1、3、5、7、14天六个时间点采集骨骼肌样本（每个时间点n=6）。

2.2 样本采集与微生物DNA测序

无菌条件下采集后肢深层肌肉组织，采用细菌基因组DNA提取试剂盒提取总微生物DNA，扩增16S rRNA基因V3-V4区（引物341F/805R），在Illumina NovaSeq 6000平台上进行PE250测序。

2.3 生物信息学与统计分析

使用QIIME2流程处理原始数据，DADA2算法生成ASV特征表。采用SILVA数据库（Release 138）进行物种注释。在R语言环境中计算α多样性（Chao1、Shannon指数）和β多样性（Bray-Curtis距离），使用PERMANOVA分析组间差异。以微生物属水平相对丰度为特征，分别构建随机森林（RF）分类模型（区分早、晚期）和回归模型（预测PMSI），通过AUC和MAE评估模型性能。

3 结果

3.1 微生物群落组成动态变化

骨骼肌微生物群落呈现显著时间依赖性演替。在门水平，随着PMSI延长，变形菌门（Proteobacteria）相对丰度从初期（<20%）持续上升至14天时超过60%，而厚壁菌门（Firmicutes）丰度逐步下降。属水平上，早期阶段（0-3天）以假单胞菌属（Pseudomonas）、不动杆菌属（Acinetobacter）和沃格菌属（Vogesella）为主；晚期阶段（5-14天）则转变为梭菌属（Clostridium）和类梭菌属（Paeniclostridium）优势。Sankey图进一步显示晚期厚壁菌门的增加几乎完全由梭菌属和类梭菌属贡献。

3.2 死因对微生物多样性的影响有限

PCoA分析显示D组与PS组样本点在排序空间中高度混杂，95%置信椭圆大量重叠。PERMANOVA证实两组间群落结构无显著差异（R²=0.008, p=0.251），死因仅能解释0.8%的群落变异。α多样性分析（Chao1、Shannon指数）在各时间点均未显示组间显著差异（p>0.05），表明死因并非塑造骨骼肌微生物演替模式的主要驱动因素。

3.3 微生物群落演替的阶段划分

基于Bray-Curtis距离的PCoA和NMDS分析均显示样本点沿清晰时间梯度分布，可明确划分为早期（EP, 0-3天）和晚期（LP, 5-14天）两个阶段。PERMANOVA证实两阶段群落结构差异极显著（p<0.001），95%置信椭圆在排序空间中几乎无重叠。

3.4 基于阶段演替模型的PMSI预测

比较RF、LR、KNN、SVM四种分类模型，RF模型在独立测试集上表现最佳，区分早、晚期的准确率达90.9%，AUC为0.9504。随后构建的阶段特异性RF回归模型显示：早期模型预测极精准，测试集MAE仅0.303天；晚期模型测试集MAE为1.293天。这种"先分类后回归"的二阶段建模策略有效提升了PMSI预测精度。

4 讨论

本研究首次系统揭示了大鼠骨骼肌微生物群在水体浸没后呈现规律的二阶段演替模式，并成功构建了高精度PMSI预测模型。骨骼肌作为相对封闭的内部组织，其微生物演替受外部环境干扰小，死因影响不显著，这为法医学提供了比体表微生物更稳定的"微生物时钟"。演替模式（从变形菌门到厚壁菌门/梭菌属的转变）反映了内部微环境从缺氧到严重分解的生态转折，与陆地分解模型有一定共性。水-肌肉微生物共享分析表明，尽管溺死过程引入了大量水体微生物，但骨骼肌内部的厌氧微环境可能通过竞争排斥作用筛选了优势菌群，维持了核心演替轨迹的稳定性。研究的局限性包括物种差异、未采用积温（ADD）更精确建模以及未深入探讨功能代谢变化。

5 结论

大鼠骨骼肌微生物组在死后水体浸没过程中发生规律且可预测的演替，与死因无关。基于此构建的二阶段机器学习模型可实现高精度PMSI估算，为法医学提供了新型稳定的研究靶点和建模范式。未来需在更接近人类的模型和复杂环境中验证，并结合宏功能基因组学等方法，推动该"微生物时钟"转化为法医学实践的有力工具。