Integrated machine learning for cause-of-death classification and postmortem interval prediction: Liver and kidney metabolomics from seawater-immersed rat cadavers

Jianghuan Lu, Yuzhao Xu, Siqi Chen, Zhiao Duan, Yixin Ma, Xiaoshi Qin, Yunchao Zhou, Shan Ha, Jianhua Chen, Deng Jianqiang

PURPOSE: To assess whether liver and kidney metabolomics combined with machine learning can distinguish seawater drowning from postmortem submersion after CO2 euthanasia and estimate postmortem interval (PMI) under controlled conditions. METHODS: Sixty male Sprague-Dawley rats were assigned to a seawater-drowning group or a postmortem-submersion group after CO2 euthanasia (30 per group). Liver and kidney tissues were collected at 0, 12, 24, 36, 48, and 72 h after death and analyzed by untargeted LC-MS/MS. Principal component analysis (PCA) and cross-validated orthogonal partial least squares discriminant analysis (OPLS-DA) were used to characterize global metabolic variation. Random Forest (RF), Support Vector Machine (SVM), Multi-Layer Perceptron (MLP), and Gradient Boosting Decision Tree (GBDT) were evaluated using repeated 10-fold cross-validation (five repeats; 50 folds in total). PMI-specific XGBoost models used 20 metabolites selected within each training fold by absolute Spearman correlation and were assessed by repeated 10-fold and leave-one-time-point-out cross-validation. SHAP analysis summarized feature contributions. RESULTS: PCA showed that dominant metabolic variation was mainly associated with PMI, with substantial overlap between the two groups in the PC1-PC2 space. Cross-validated OPLS-DA identified group-associated structure in liver and kidney, with Q2(cum) values of 0.847 and 0.723, respectively. Mean classification AUCs ranged from 0.973 to 0.996 in liver and from 0.931 to 1.000 in kidney. Under repeated 10-fold cross-validation, the liver and kidney PMI models achieved MAEs of 4.78 and 4.05 h and R2 values of 0.823 and 0.892, respectively. Under leave-one-time-point-out cross-validation, MAEs increased to 13.52 and 13.64 h, while R2 values were 0.587 and 0.592. SHAP analysis showed partly different metabolite contribution patterns between the two organs. CONCLUSION: In this controlled rat model, liver and kidney metabolomic profiles combined with machine learning supported discrimination between seawater drowning and postmortem submersion after CO2 euthanasia and captured PMI-related metabolic changes. The classification models, PMI-specific XGBoost regression models, and SHAP analyses characterized organ-specific metabolite patterns associated with discrimination between the two modeled conditions and PMI prediction.