Differentiating hemorrhagic shock and organophosphate poisoning through integrated skin microbiome–metabolome signatures

Accurate determination of cause of death and estimation of postmortem interval (PMI) are critical yet challenging tasks in forensic science, particularly in cases with rapid demise and absence of obvious morphological abnormalities. We employed an integrative multi-omics approach to characterize postmortem microbial succession and metabolic alterations on facial skin in mouse models of hemorrhagic shock (HS) and organophosphorus poisoning (OP) across three decomposition stages: bloating (2 days), active decay (8 days), and advanced decay (16 days). Metagenomic profiling revealed significantly reduced α-diversity in HS compared with OP throughout all stages ( p < 0.001), accompanied by stage-dependent compositional shifts, including early enrichment of Firmicutes in HS and Proteobacteria in OP. A total of 237 differential taxa were identified, with Providencia and Morganella predominating in OP, whereas Staphylococcus and Corynebacterium dominated bloating stage of HS. Untargeted metabolomics uncovered distinct cause-of-death–linked metabolites, notably elevated 2′-deoxycytidine-5′-diphosphate in early OP and persistent cholic acid/cholate accumulation in HS at later PMI. Functional analysis highlighted histidine and phosphate/phosphonate metabolism as key discriminatory pathways, exhibiting stage-specific oscillations and strong correlations with characteristic taxa. These findings demonstrate that skin-based metagenomic–metabolomic integration provides robust, mechanistically informed biomarkers for both PMI estimation and cause-of-death differentiation, offering a minimally invasive and temporally dynamic tool for forensic investigations.