Analysis of postmortem intestinal microbiota successional patterns with application in postmortem interval estimation

DOI: https://doi.org/10.21203/rs.3.rs-411962/v1

Abstract

Microorganism plays a vital role in the decomposition of vertebrate remains in nature nutrient cycling, where the postmortem microbial succession patterns during decomposition remain unclear. Here, hierarchical clustering based on Manhattan distances analyzed the similarity and difference among postmortem intestine microbial succession patterns from microbial 16S rDNA sequences in a mouse decomposition model. Based on the similarity, seven different classes of succession patterns are obtained. Generally, the normal intestinal flora in the cecum gradually decreases with changes in living conditions after death, while some facultative anaerobes and obligate anaerobes grew and multiplied with oxygen consumption. Furthermore, a random forest regression model is developed to predict postmortem interval with the microbial succession trend data set. The model demonstrates a mean absolute error of 20.01 hours and a squared correlation coefficient of 0.95 during 15-day decomposition. Lactobacillus, Dubosiella, Enterococcus, and Lachnospiraceae NK4A136 group are considered as significant biomarkers for this model according to the ranked list. This study explores microbial succession patterns in terms of relative abundances and variety, aiding in predicting postmortem intervals, offering some knowledge of microbial behaviors in decomposition ecology.

Introduction

As a key aspect of the material cycle and energy flow of the ecosystem, the decomposition of vertebrate remains is one of the main processes involved in maintaining ecosystem functions. This process begins immediately after death. A lack of oxygen results in cell autolysis, and tissue breakdown, releasing cellular components and macromolecules. These nutrients are subsequently metabolized by predominately microbial communities, particularly those concentrated in the gastrointestinal tract. These microbes begin to exponentially proliferate, transmigrate, and create specialized proteins that digest host tissues during putrefaction [1]. That the intestinal microbiome plays an essential role in the remaining decomposition [2,3]. However, the number of studies of the dynamic succession patterns in these microbial communities is relatively small, but these studies may facilitate estimation of the postmortem interval (PMI) in forensic investigations and improve understanding of the mechanisms governing succession processes during decomposition to accelerate decomposition, facilitate nutrition recycling, and reduce environmental pollution.

Recently, some studies reported microbial succession trends during decomposition. For example, Liu R et al. [4] established a microbial clock providing an accurate estimate of PMI in a mouse model system and detected Enterococcus faecalis, Anaerosalibacter bizertensis, Lactobacillus reuteri, and others as the most informative species in the decomposition process for PMI estimation. Zhang J. et al [5]. also established a random forest model to estimate PMI with burial microbial sequence data set. Additionally, Debruyn et al. [6] documented that Bacteroidales (Bacteroides, Parabacteroides) significantly declined while Clostridiales (Clostridium, Anaerosphaera) and the fly-associated Gammaproteobacteria Ignatzschineria and Wohlfahrtiimonas increased in human gut bacterial communities. Adserias et al. [7] monitored the oral microbiota of decaying human bodies and found that Firmicutes and Actinobacteria were the predominant phyla in the fresh stage, while the presence of Tenericutes corresponded to the bloat stage, and Firmicutes was the predominant phylum in advanced decay. As decomposition progressed, Firmicutes replaced Proteobacteria as the dominant phylum in the late decay stages of swine models [8]. Additionally, researchers found that the relative abundances of Clostridiaceae and Enterococcaceae increased during decomposition in murine models [2,9]. These similar results regarding microbial succession between animal decomposition models and human remains showed the possibility to use animal decomposition models as alternatives for microbial succession in human corpses to mathematically establish microbial community change models during decomposition.

However, these studies only partially focused on microbial succession during decomposition. Few studies have tried to mathematically describe microbial succession patterns. Fortunately, recent studies have shown that utilizing mathematical models to describe temporal patterns of microbial diversity is a developing strategy that facilitates the inference of relative growth rates, mechanisms of interaction, and responses to perturbations [10-12]. A recent increase in the number of microbial time series studies with mathematical methods offers some insights into the stability and dynamics of microbial communities in a wide range of environments [13-16]. In a meta-analysis of temporal dynamics in microbial communities (including 76 sites representing air, aquatic, soil, brewery wastewater treatment, human- and plant-associated microbial biomes), approximately half had time-decay curves with negative slopes [15]. Besides, the temporal variability in microbial community diversity was comparable across studies within the same environment but varied across environments, being lowest in soil and brewery wastewater and highest in the human palm and the infant gut [15,16]. Some microbial communities go through a series of predictable states after colonization (primary succession), which occurs for instance during the formation of dental plaque, where oxygen-tolerant early colonizers prepare the ground for later oxygen-sensitive colonizers [17]. In some cases, such as infant gut microbiota colonization [18], communities vary considerably in the initial stages of succession, but stabilize at similar states. These studies indicated the potential of using mathematical methods in microbial dynamic analysis.

Thus, this study aimed to explore intestinal microbial succession patterns during decomposition by a mathematical algorithm and establish a PMI estimation model based on successional trend data set. The statistically significant changes in microbial relative abundance during decomposition were analyzed by the Kruskal-Wallis test (P< 0.001). Utilizing hierarchical classification analysis to classify the successional patterns based on the Manhattan distance, seven classes of successional patterns were obtained. Moreover, A random forest model combined with the postmortem intestinal microbial successional trend data set yielded a mean absolute error of 20.01 hours within 15-day decomposition. Some initial colonizers of the mammalian gastrointestinal microbiota and some obligately anaerobic microbes were regarded as important taxa in this model according to the ranked list of all taxa.

Materials And Methods

Sample collection and experimental setup

All experiments were approved by the Institutional Animal Use and Care Committee of Xi’an Jiaotong University (approval No: 2017-288). Adult male C57BL/6J mice (18-25 g, 6-10 weeks, n = 240) were obtained from the Experimental Animal Center of Xi’an Jiaotong University. All mice were housed in micro isolator cages under standard lighting (light/dark periods of 12h) and were given autoclaved food and water ad libitum. After a sacrifice by cervical dislocation, the mice were placed on sterile plates with ambient temperature (Ta, 25.0±1.5℃) and moderate relative humidity (RH, 50.0±7%) in an artificial climate box (PRX-35B, Zuole, Shanghai, China). The artificial climate box and clean bench were treated with 75% alcohol solution and an hour of UV radiation to minimize and fragment environmental contaminant DNA. The vent of the artificial climate box was wrapped by sterile gauze at 200 mesh to resist insect access. Every surgical instrument was autoclaved and burned with the outer flame of an alcohol lamp. Ten-time points (0 hours, 8 hours, 12 hours, 1 day, 2 days, 4 days, 7 days, 10 days, 13 days, 15 days after death) were set spanning 15-day decomposition. At each time point, ceca were obtained from 20 mice under strict aseptic operation by multiple people in a short period, to avoid the influence of sampling on microbiota variability. This 20*10 sample set was used for the following model training and validation. Besides, we set another four mice to every time point as an external sample set (4*10) according to the single-blind principle. Concretely, another experimenter collected postmortem ceca without telling the PMI of the samples. Samples were stored at −80°C immediately after harvest until further processing.

16S rRNA gene sequencing, curation, and community analysis.

Following bead beating, DNA was isolated from each cecum sample using the QIAamp® DNA Mini Kit (Qiagen, Germany) according to the manufacturer's specifications on a super clean bench. The variable region 3-4 (V3, V4) of the 16S rRNA gene was then amplified from each DNA sample by the bacterial primers 341F (5'-CCTAYGGGRBGCASCAG-3') and 806R (5'-GGACTACNNGGGTATCTAAT-3'). Each amplicon was sequenced on an IonS5TMXL sequencer. DNA sequences were analyzed using the mothur (version 1.33.3) pipeline. Potentially chimeric sequences were removed using the UCHIME algorithm[19]. Sequences were assigned to OTUs based on a 3% dissimilarity cutoff using the average neighbor algorithm in the UPARSE pipeline (version 7.0.1001)[20]. Each representative sequence was matched against the SSUrRNA database of SILVA reference alignments using mothur for taxonomic annotation[21]. OTU abundance information from each sample was homogenized, with the standard of the fewest sequences in the samples. Rarefaction and rank-abundance curves were drawn using R software (version 2.15.3). Alpha diversity was measured by the Shannon, Chao1, Simpson, and ACE indexes in QIIME software[22-24]. Additionally, phylogenetic beta diversity was measured using unweighted UniFrac and weighted UniFrac distances in QIIME software [22,25,26].

Data preprocessing

Every individual sample contained the relative abundance values of 5099 different taxa. First, the data were filtered to remove the noninformative taxa because the relative abundance values of many taxa were mostly zero. Specifically, samples were grouped according to the time point of death. For the data set, if the relative abundance of a taxon was zero in more than 60% of samples in all groups, then this taxon was discarded. Thus, the number of taxa was reduced from 5099 to 1810. Second, the dataset was normalized by Z-score standardization. After normalization, the relative abundance values for each taxon were scaled to be between 0 and 1 so that the change tendencies could be explored.

Microbial community dynamic analysis

The Kruskal-Wallis test was utilized to analyze the significance of the differences in the relative abundances between adjacent time points. We counted a significant increase as +1 and a significant decrease as -1. The trend values were accumulated in sequence according to the PMI time series. Finally, we obtained a microbial community succession trend dataset. In figure 2B, the red lines indicate that the changing trend of that taxon is significant, while the blue lines indicate that the changing trend of that taxon is not significant (P < 0.001, Kruskal-Wallis test). The decomposition patterns were explored by an HCA method based on the Manhattan distances between the trend values of different taxa. HCA is a commonly used unsupervised pattern recognition method that seeks to build a hierarchy of classes based on similarity measurement. The results of HCA are presented in a dendrogram. Different arrangements of classes were obtained when choosing different threshold values. Besides, the class results were significantly affected by the choice of dissimilarity or distance measurement method. In this study, the Manhattan distance, or taxicab distance, was used, in which the distance between two points is the sum of the absolute differences of their Cartesian coordinates [27], which was defined as

Evaluation of microbial succession patterns

Random forest, proposed by Breiman [28], is a popular approach in applied statistics due to its easy applicability to classification and regression problems. It creates random decision trees based on subsets of features (e.g., taxa) in the data and then chooses the subsets of features that are best able to classify samples into predefined groups using only part of the data. This technique shows high predictive accuracy and applicability even in high-dimensional problems with highly correlated variables, which often occurs in bioinformatics [29]. Recently, some researchers used RF in PMI prediction studies [8,30] and other microbial ecology studies [31], providing high classifier accuracy. In this study, a PMI prediction model was established based on the microbial community succession trend dataset mentioned above with the RF algorithm to evaluate the microbial succession patterns. The X-variable corresponded to the matrix of microbial community succession trends and the responding Y-variable was associated with PMI values. The data set was randomly divided into two parts: 80% as the training set, 20% as the validation set. Training refers to fitting or building the model while validation is equivalent to predicting. For murine body decomposition, each subset was partitioned based on an individual for cross-validation so that the samples from the same individuals are either in the training set or validation set, but not both. In the current study, the best RF regression model after hyperparameter tuning through cross-validation was selected to represent the final model. To evaluate the stability and predictive ability of the RF model, the external set was preprocessed and analyzed following the above method and then formed the testing data set. This external testing set was performed using constructed models to predict the PMI. The MAE, RMSE, and , as the three main parameters of the model calibrated and predicted results, were used to evaluate the regression model reliabilities. The MAE was calculated as the deviation of the predicted from observed values and representing the average prediction error in the same unit of the original data. The RMSE represents the quadratic mean of these differences between predicted values and observed values. The can evaluate the goodness of fitting between actual and predicted values. High values of and low values of RMSE and MAE demonstrate a well-established RF regression model. In particular, the model was run and evaluated fifteen times in the validation and test sets, and the mean value of the MAE and as shown in Table 1. The modeling was done with Matlab software version R2019b (MathWorks, Natick, MA, USA) equipped with Windows-Precompiled-RF_Mexstandalone-v0.02toolbox (https://randomforest-matlab.googlecode.com/files/Windows-Precompiled RF_Mexstandalone-v0.02-.zip), and the number of trees and predictors sampled for splitting at each node were both 500.

Additionally, mean decrease accuracy (MDA) [32] values were used to evaluate variable importance in this PMI prediction model. The MDA takes the importance of out-of-bag error and the impurity of the variable into consideration. The high value of the MDA showed the great importance of the variables to this model.

Statistics and reproducibility

Graphical illustrations are plotted using Prism software (GraphPad, San Diego, CA). The P-value cut-off used in this study is 0.05. *, **, and *** in the figures refer to P-values ≤0.05, ≤0.01, and ≤0.001, respectively. For all comparisons between two groups, the unpaired Student's t-test or the Kruskal-Wallis test was performed depending on the data normality and homogeneity of variance. We considered biological replicates the sequencing results of biologically distinct samples.

Results

Overview of microbial composition during decomposition

The average raw library size after paired-end assembly was 78,545 ± 10,670 (SD). After sequence trimming, quality filtering, and removal of chimeras, 17,696,903 high-quality sequences remained, with an average of 73,771 ± 9,688 (SD) reads per sample. Samples with library sizes close to, or larger than, 3000 reads showed convergence in their sampling effort curves. These high-quality sequences were assigned to a total of 5,099 operational taxonomic units (OTUs) by the UPARSE pipeline using a threshold of 97% identity. The OTUs belonged to 29 phyla, 67 classes,160 orders, 250 families, 473 genera, and 835 species. Among the 29 phyla, Firmicutes, Bacteroidetes, and Actinobacteria were dominant in all samples. (Fig. 1a). Lactobacillus was dominant in all groups at the genus level (Fig. 1b). The Lachnospiraceae NK4A136 group accounted for more than 15% of the total sequences before day 1. At the species level, Lactobacillus reuteri, Lactobacillus johnsonii, Enterococcus faecalis, Firmicutes bacterium M10-2, and Clostridium tetani E88 were the most abundant species (Fig. 1c). In general, the microbial community was dominated by OTUs with the Firmicutes (Fig. 1d) affiliated within the Lactobacillaceae.

Dynamics of microbial communities during decomposition

Although there appeared to be a significant amount of stochastic variation across the duration of this experiment, the overall reproducibility of these patterns between animals suggested that deterministic microbial community succession decreased with time (Fig. 1). As a result, only the following genera remained at day 15 postmortem: Gordonibacter, Bifidobacterium, Enterorhabdus, Lactococcus, Clostridium sensu stricto 18, Clostridium sensu stricto 15, Anaerosalibacter, Enterococcus, Dubosiella, and Lactobacillus (Fig. 1b). At the species level (Fig. 1c), the relative abundance of Lactobacillus reuteri, Clostridium tetani E88, Lactobacillus johnsonii, and Enterococcus faecalis showed notable variation during decomposition, which may aid in PMI estimation. Lactobacillus reuteri increased and exhibited peaks in abundances at both day 7 and day 15. Clostridium tetani E88 appeared on day 7 and decreased until day 15. Lactobacillus johnsonii showed a higher relative abundance one week after death than in the advanced stage of decomposition. Besides, Firmicutes bacterium M10-2 appeared on day 2 and showed an immediate increase from day 2 to day 4. Enterococcus faecalis appeared on day 2 and increased until day 10. Clostridium tetani E88 appeared on day 4 and increased and then decreased from then on. At the OTU level, OTU 1613 gradually increased in relative abundance during the 15-day decomposition and finally became the dominant OUT at day 15 (Fig. 1d). Among the Bacteroidetes and Actinobacteria, the OTUs had a similar change trend, increasing up to their peak abundances at day 2 and decreasing until day 4. However, OTU 2202, annotated to Helicobacter, gradually decreased during the 15-day decomposition (Fig. 1d). A heat map also showed discrepancies in microbial community abundance change trends during decomposition, with a color gradient from deep blue (low abundance) to deep red (high abundance) (Fig. 1e).

Classification of microbial succession patterns during decomposition

To analyze the microbial succession patterns, we applied hierarchical cluster analysis (HCA) based on the Manhattan distances between different taxa (Fig. 2a). Seven classes of taxon change patterns were obtained (Fig. 2b). Classes 1, 2, and 5showed a similar decreasing trend, which appeared with different PMIs. The taxa in lass 1 and class 5exhibited a fast decline before four days after death and a slow decline from then on. However, the declination range was greater in class 5 than in class 1. The taxa in Class 2 exhibited a continuous decrease during 15-day decomposition. In contrast, classes 3 and 4 showed an increase during the early PMI points. Class 3 showed a small increase two days after death, while class 4 showed a fast increase, a slow increase during 7-day decomposition, and a decrease from day 7 to day 13. Class 7 and class 6 assembled in clustering tree as both contained a lot of taxa with non-significant change trend during decomposition. Class 6 showed a relatively continuous profile, while class 7 showed a slow decline during the whole 15-day decomposition. After species annotation, class 1 mainly contained Lachnospiraceae, Ruminococcaceae, and Muribaculaceae at the family level and Lachnospiraceae NK4A136 group, Lactobacillus, Bacteroides, Lachnoclostridium, and Prevotellaceae UCG-001 at the genus level. Class 2 mainly contained Muribaculaceae, Lachnospiraceae, Ruminococcaceae, and Lachnospiraceae at the family level, Lactobacillus, Lachnospiraceae NK4A136 group, Ruminococcaceae UCG-014, Ruminiclostridium 5, Helicobacter, Ruminiclostridium 9, Alistipes, and Lachnoclostridium at the genus level. The microbial composition in class 1 and class 2 were similar. Class 5 also mainly contained Muribaculaceae, Lachnospiraceae, and Ruminococcaceae at the family level, and Lachnospiraceae NK4A136 group, Ruminococcaceae UCG-014, Lachnoclostridium, Bacteroides, Prevotellaceae UCG-001, Lactobacillus, Helicobacter, and Ruminiclostridium 9 at the genus level. Class 7 contained Enterorhabdus and Dubosiella. Dubosiella, Lactobacillus, Lactococcus, and Enterococcus were dominant genera in class 3 and class 4, increasing during early decomposition. We also counted the number of taxa in each class (Fig. 2b). Classes 6 and 7 contained 58.95% of the microbial taxa in all patterns. However, class 4 contained the fewest microbial taxa among all classes and was unimodal with an increasing profile. The blue lines in the plot represent relatively invariable taxa based on the Kruskal-Wallis test, which were mostly included in class 6. The other classes are mainly represented by the red lines, with significant variation during decomposition (P < 0.001, Kruskal-Wallis test). To display the changes in the patterns more obviously, the relative abundance changes of microbes in every class are shown in Fig. 3 (at the phylum level) and Supplementary Figure 1 (at the genus level). Specifically, in class 1, the relative abundance of Firmicutes gradually decreased, while Proteobacteria increased with PMI. The relative abundance of Bacteroidetes decreased during decomposition in class 2. Firmicutes in class 3 gradually increased. Latescibacteria in class 4 increased during 15-day decomposition. Bacteroidetes in class 5 gradually decreased during decomposition. The phyla in class 6 and class 7 did not show obvious variation during decomposition. In conclusion, the Lachnospiraceae NK4A136 group, Lactobacillus, Muribaculaceae, and Lachnospiraceae showed a declining profile, while Dubosiella, Enterococcus, Lactococcus, Clostridium_sensu_stricto_15, and Enterorhabdus were mainly increased during the 15-day decomposition. Then, we calculated the percentage of OTU annotation for the top 20 genera and every phylum. Firmicutes were the most dominant phylum in all classes (Fig. 4a). Class 6 contained the greatest number of phyla among the classes. Class 4 contained the fewest phyla, including Firmicutes, Actinobacteria, Bacteroidetes, Deferribacteres, Epsilonbacteraeota, and Proteobacteria. At the genus level (Fig. 4b), Lachnospiraceae NK4A136 group was the dominant genus in class 1. Lactobacillus was dominant in class 2. Dubosiella, Lactobacillus, and Enterococcus were dominant genera in class 3. Lactobacillus was the most dominant genus in class 4. Class 6 and class 7 contained many genera, while the distribution of the taxa in class 6 was discrete. Class 4 contained the fewest genera of the classes.

A PMI estimation model based on the microbial community succession trend data set

A random forests (RF) regression model was established to predict the PMI based on the microbial succession trend data set. The root mean square error (RMSE), the mean absolute error (MAE), and the squared correlation coefficient (R-squared, ) values in the validation set were 30.03 hours, 20.42 hours, 0.945, respectively. Those values in the testing set were 28.38 hours, 20.01 hours, and 0.947, respectively. Besides, the importance of variables was evaluated using the mean decrease accuracy (MDA) in the prediction illustrated in Fig. 5. In this model, the taxa with MDA values above 17.93, as 1.5 times the mean value of all MDA values, were regarded as significant strains, which was represented by the orange line in Fig. 5. The microbial taxa with MDA values over 17.93 were from all classes. According to the species annotation, significant taxa were mainly related to Lactobacillus, Dubosiella, Enterococcus, Lachnospiraceae NK4A136 group, Prevotellaceae UCG-001 and Clostridium sensu stricto 15 at the genus level.

Discussion

Decay and decomposition in mammalian and other vertebrate taxa are a key process in biological nutrient cycling. Cooperation between vertebrate and invertebrate scavengers (including bacteria, archaea, fungi, and protists) plays an important role in the chemical decomposition of animal waste. In particular, microbial decomposers are expected to lead to the conservation of key biochemical metabolic pathways and cross-kingdom ecological interactions for the efficient recycling of nutrient reserves. When a mammalian body is decomposing, microbial and biochemical activities result in a series of decomposition stages[9] that are associated with reproducible microbial succession across mice[4], swine[8], and human cadavers. Thus, a deep investigation of the behavior of this group during decomposition could be very helpful.

In the current study, we utilized hierarchical classification analysis with Manhattan distances to classify the relative abundance changes in the microbiota during decomposition. Seven classes of succession patterns were detected in the current study. The seven changing patterns could be further divided into unimodal, continuous decline, and a single peak followed by an ascending phase. The declining pattern was the most dominant pattern of microbial succession. We focused on intestinal microbiome changes during decomposition. At the very beginning of death, the microbes in the cecum were almost all normal intestinal tract bacterial flora. These microbes play a vital role before death in functions such as the decomposition of food, degradation of drugs, promotion of fat accumulation, and synthesis of vitamins and amino acids. Once the body dies, the living environment of the microbial organism changes, inducing a decline in microbial diversity and richness. The Chao1 estimator showed a decline in microbial community richness within 4 days of death and a subsequent gradual increase until 13 days (Supplementary Figure 2a). The Shannon index demonstrated that the diversity of the microbial community also declined during 15-day decomposition (Supplementary Figure 2b). These results coincided with the dominant microbial succession pattern profiles (Fig. 2b).

However, there remained discrepancies among 7 classes of postmortem microbial successional patterns. Class 1, 2, and 5 showed similar decrease patterns. Specifically, class 1 showed a fast decline four days after death and a slow decline from then on. The microbiome in classes 2 and 5 exhibited a continuous decrease during 15-day decomposition, while the range of decrease was higher in class 5 than in class 2. The reasons for these changes may be that the microbial compositions in these classes differ. After species annotation, the microbial compositions in classes 1 and 2 were similar. A hierarchical cluster plot (Fig. 2a) also showed the similarity between class 1 and class 2. However, the relative abundances of Muribaculaceae and Lactobacillus were higher in class 2, while the relative abundance of Lachnospiraceae NK4A136 group was higher in class 1. Lachnospiraceae NK4A136 group, a genus of the gut microbiota, is involved in the pathogenesis of age-related cognitive impairment in mice [33]. Classes 1, 2, and 5 all contained Lachnospiraceae, Ruminococcaceae, and Muribaculaceae at the family level. The bacteria in Ruminococcaceae are obligate anaerobes, while Lachnospiraceae is a family of anaerobic, spore-forming bacteria in the order Clostridiales that ferment diverse plant polysaccharides to short-chain fatty acids (butyrate, acetate) and alcohols (ethanol). Based on NCBI taxonomy, Muribaculaceae, Lachnospiraceae, and Deferribacteraceae were the predominant families identified in the healthy mouse gut microbiome. Researchers [34] also found that the genes for carbohydrate metabolism were upregulated in Muribaculaceae, while genes for cofactors and vitamin metabolism and amino acid metabolism were upregulated in Deferribacteraceae. These bacteria were reported in the mammalian gut microbiota [35,36]. Class 7 showed a slow decline during 15-day decomposition.The microbial composition of class 7 was slightly different from those of classes 1, 2, and 5, which was also demonstrated in the hierarchical cluster results. Class 7 included Enterorhabdus and Dubosiella, while classes 1, 2, and 5 did not. Some studies reported that the relative abundance of the genus Enterorhabdus increased in subjects with prediabetes compared to healthy subjects [37]. In conclusion, classes 1, 2, and 5 showed similar decreasing profiles, and they were mainly annotated to Ruminococcaceae, Lachnospiraceae,Deferribacteraceae, and Muribaculaceae at the family level. The transformations listed above might be caused by changes in the internal conditions of the mouse after death; in particular, the oxygen content decreased, which was not suitable for obligate anaerobes. Conversely, classes 3 and 4 showed an increased profile during the early PMIs. Class 6 showed a relatively invariable profile during the whole 15-day decomposition. Dubosiella, Lactobacillus, Lactococcus, Clostridium sensu stricto 15, and Enterococcus were dominant genera in class 3 and class 4, with increasing patterns during early decomposition. Dubosiella, a gram-positive genus from the family Erysipelotrichidae, has been isolated from the intestinal content of mice. Lactococcus species are capable of proteolysis, a process that converts amino acids into flavor complexes, which may aid in decomposing the protein of corpses. Lactobacillus is a genus of gram-positive, facultative anaerobic or microaerophilic, rod-shaped, non-spore-forming bacteria. The relative abundances of the genus Enterococcus, as facultative anaerobic and gram-positive bacteria, suggested that as time continued after death, the proportion of normal intestinal flora became maladjusted, and opportunistic pathogens rose immediately. Classes 3 and 4 also contained Enterorhabdus at the genus level, which are gram-positive, facultative anaerobe, non-spore-forming rods that exist in mammalian intestines [38]. The possible reason for this facultative anaerobe increase is that body rupture permits oxygen entry. At the species level, Lactobacillus reuteri and Lactobacillus animals were dominant in class 4. These species were reported as commensal probiotic bacteria that produce antimicrobials to fight off harmful infections and mediate the immune system [39,40]. They both showed the highest relative abundances on day 7 after death. This finding provides the possibility of extracting useful microbes from the decomposing microbiota. Thus, local microbes are most dominant in the intestinal tract at the early stage. With corpse bloating, the oxygen inside the body is consumed by the decomposition of nutrient substances. Meanwhile, the bacteria of the natural gut microbiome spread as tissue liquefies. When the body enters advanced decay, there is a proliferation of obligate anaerobes. After body rupture, facultative anaerobes proliferated inside the gut track. Moreover, we observed the relative abundance change patterns of each OTU in every class. The original curves of taxa were aligned with the fitting curves of every class. Thus, it is obvious that the classifying and fitting methods used in this study are suitable and reliable. Interestingly, the succession patterns showed similar trends in general, and the change rates for 4-day decomposition were higher than those for the subsequent processes. From this point, the microbial decomposers may have similar behavior during the degradation of bodies, which improves understanding of the ecological cycling rules in the microbial kingdom.

A RF model was established with high at 0.95 and low MAE at 20.01 hours during 15-day decomposition based on postmortem microbial successional trend data (Table 1). To evaluate the PMI prediction ability of each class, the MDA value was calculated for every microbial taxon. Classes 2 and 7 contained the largest number of important variables for predicting PMI based on the MDA value of each taxon (Fig. 5). However, class 6 contained the fewest among the 7 classes. This result is explained in Fig. 2b that class 6 contained many invariable taxa, while the abundances of members of classes 2 and 7 were constantly changing during the 15-day decomposition. The important taxa in class 2 were mainly initial colonizers of the mammalian gastrointestinal microbiota, most of which are anaerobic or aerotolerant anaerobic bacteria, such as Helicobacter, Bacteroides, Alistipes, Muribaculum, Lactobacillus, Enterorhabdus, Adlercreutzia, Anaerotruncus, and Mucispirillum. When the bodies died, the living conditions of the intestinal microbiota gradually changed, resulting in a decline in the relative abundances of these initial microbes. The significant taxa in class 7 were mainly obligately anaerobic microbes, for example, Desulfovibrio (found in a variety of habitats, including soil and the intestines, and feces of animals), Roseburia (gram-positive anaerobic bacteria that inhabit the human colon), Anaerovorax (strictly anaerobic), and Allobaculum (strictly anaerobic). The relative abundances of these anaerobic microbes decreased immediately after body rupture, which may facilitate the prediction of the early PMI. As a result, Lactobacillus, Dubosiella, Enterococcus, Lachnospiraceae NK4A136 group, and Clostridium sensu stricto 15 were the most important genera for the PMI estimation model, which also coincides with results of previous studies.

This study provides basic knowledge for future studies on microbial decomposers. However, there are still knowledge gaps existing in the natural decomposing kingdom. Future work needs to focus on the succession patterns of fungi, heterotrophic organisms (scavengers), and their relationship with PMI estimation. Our group is also dedicated to expanding the application of this work to human remains in forensic practice.

Conclusion

The present study contributes to our knowledge about postmortem microbiology as the thorough characterization of microbial succession patterns. We demonstrated that postmortem microbial communities differ over time. However, their changing trends have some similarities, as the successional rate changes faster at early decomposition and slower later. A random forest model was established to predict the PMI with an MAE of 20.01 hours during 15-day decomposition based on the trend data set. Taken together, these succession patterns are confirmed as critical information to consider when developing the use of microorganisms as physical evidence for PMI estimation.

Declarations

Funding

This work has been supported by the funders of the National Natural Science Foundation of China (No. 81730056).

Conflicts of interest

The authors declare that there is no conflict of interest.

Availability of data and material

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Code availability

The code used in this work is available on request from the corresponding author.

Authors' contributions

Conceptualization: Ruina Liu, Zhenyuan Wang; Methodology: Ruina Liu, Qi Wang, Kai Zhang; Formal analysis and investigation: Ruina Liu, Hao wu, Wumin Cai, Gongji Wang; Writing - original draft preparation: Ruina Liu; Writing - review and editing: Zhenyuan Wang, Ruina Liu, Kai yu; Funding acquisition: Zhenyuan Wang; Resources: Qinru Sun; Supervision: Shuanliang Fan.

Ethics approval

This work used only dead animals that had been sacrificed as a byproduct of other research and were deemed not to be animal research by the Institutional Animal Use and Care Committee of Xi’an Jiaotong University. Mice were euthanized humanely under approved approval No: 2017-288.

Consent to participate

Not applicable

Consent for publication

Not applicable

Acknowledgments

This work has been supported by the funders of the National Natural Science Foundation of China (No. 81730056). The authors also thank Novogene Co., Ltd. for insightful assistance with the DNA sequencing and metagenomics analysis. The supporters had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

References

  1. Carter DO, Yellowlees D, Tibbett M (2007) Cadaver decomposition in terrestrial ecosystems. Naturwissenschaften 94 (1):12-24. doi:10.1007/s00114-006-0159-1
  2. Liu R, Gu Y, Shen M, Li H, Zhang K, Wang Q, Wei X, Zhang H, Wu D, Yu K, Cai W, Wang G, Zhang S, Sun Q, Huang P, Wang Z Predicting postmortem interval based on microbial community sequences and machine learning algorithms. n/a (n/a). doi:10.1111/1462-2920.15000
  3. Burcham ZM, Hood JA, Pechal JL, Krausz KL, Bose JL, Schmidt CJ, Benbow ME, Jordan HR (2016) Fluorescently labeled bacteria provide insight on post-mortem microbial transmigration. Forensic Science International 264:63-69. doi:https://doi.org/10.1016/j.forsciint.2016.03.019
  4. Liu R, Gu Y, Shen M, Li H, Zhang K, Wang Q, Wei X, Zhang H, Wu D, Yu K, Cai W, Wang G, Zhang S, Sun Q, Huang P, Wang Z (2020) Predicting postmortem interval based on microbial community sequences and machine learning algorithms. Environmental Microbiology 22 (6):2273-2291. doi:10.1111/1462-2920.15000
  5. Zhang J, Wang M, Qi X, Shi L, Zhang J, Zhang X, Yang T, Ren J, Liu F, Zhang G, Yan J (2021) Predicting the postmortem interval of burial cadavers based on microbial community succession. Forensic Science International: Genetics 52:102488. doi:https://doi.org/10.1016/j.fsigen.2021.102488
  6. Debruyn JM, Hauther KAJP (2017) Postmortem succession of gut microbial communities in deceased human subjects. 5 (6):e3437
  7. Adserias-Garriga J, Quijada NM, Hernandez M, Rodríguez Lázaro D, Steadman D, Garcia-Gil LJ (2017) Dynamics of the oral microbiota as a tool to estimate time since death. Molecular oral microbiology 32 (6):511-516. doi:10.1111/omi.12191
  8. Pechal JL, Crippen TL, Benbow ME, Tarone AM, Dowd S, Tomberlin JKJIJoLM (2014) The potential use of bacterial community succession in forensics as described by high throughput metagenomic sequencing. 128 (1):193-205
  9. Metcalf JL, Xu ZZ, Weiss S, Lax S, Van Treuren W, Hyde ER, Song SJ, Amir A, Larsen P, Sangwan N, Haarmann D, Humphrey GC, Ackermann G, Thompson LR, Lauber C, Bibat A, Nicholas C, Gebert MJ, Petrosino JF, Reed SC, Gilbert JA, Lynne AM, Bucheli SR, Carter DO, Knight R (2016) Microbial community assembly and metabolic function during mammalian corpse decomposition. Science (New York, NY) 351 (6269):158-162. doi:10.1126/science.aad2646
  10. Larsen P, Hamada Y, Gilbert J (2012) Modeling microbial communities: Current, developing, and future technologies for predicting microbial community interaction. Journal of Biotechnology 160 (1):17-24. doi:https://doi.org/10.1016/j.jbiotec.2012.03.009
  11. O'Dwyer JP, Lake JK, Ostling A, Savage VM, Green JL (2009) An integrative framework for stochastic, size-structured community assembly. 106 (15):6170-6175. doi:10.1073/pnas.0813041106 %J Proceedings of the National Academy of Sciences
  12. Gerber GK, Onderdonk AB, Bry L (2012) Inferring dynamic signatures of microbes in complex host ecosystems. PLoS computational biology 8 (8):e1002624-e1002624. doi:10.1371/journal.pcbi.1002624
  13. David LA, Materna AC, Friedman J, Campos-Baptista MI, Blackburn MC, Perrotta A, Erdman SE, Alm EJ (2014) Host lifestyle affects human microbiota on daily timescales. Genome biology 15 (7):R89. doi:10.1186/gb-2014-15-7-r89
  14. Chow C-ET, Kim DY, Sachdeva R, Caron DA, Fuhrman JA (2014) Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists. The ISME Journal 8 (4):816-829. doi:10.1038/ismej.2013.199
  15. Shade A, Gregory Caporaso J, Handelsman J, Knight R, Fierer N (2013) A meta-analysis of changes in bacterial and archaeal communities with time. The ISME Journal 7 (8):1493-1506. doi:10.1038/ismej.2013.54
  16. Xia LC, Ai D, Cram J, Fuhrman JA, Sun F (2012) Efficient statistical significance approximation for local similarity analysis of high-throughput time series data. Bioinformatics 29 (2):230-237. doi:10.1093/bioinformatics/bts668
  17. Kolenbrander PE, Andersen RN, Blehert DS, Egland PG, Foster JS, Palmer RJ (2002) Communication among Oral Bacteria. Microbiology and Molecular Biology Reviews 66 (3):486. doi:10.1128/MMBR.66.3.486-505.2002
  18. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO (2007) Development of the human infant intestinal microbiota. PLoS biology 5 (7):e177. doi:10.1371/journal.pbio.0050177
  19. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) UCHIME improves sensitivity and speed of chimera detection. 27 (16):2194
  20. Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, Ciulla D, Tabbaa D, Highlander SK, Sodergren E, Methe B, DeSantis TZ, Petrosino JF, Knight R, Birren BW (2011) Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome research 21 (3):494-504. doi:10.1101/gr.112730.110
  21. Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JIJN The Human Microbiome Project. 449 (7164):804-810
  22. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7 (5):335-336. doi:10.1038/nmeth.f.303
  23. Tai N, Peng J, Liu F, Gulden E, Hu Y, Zhang X, Chen L, Wong FS, Wen L (2016) Microbial antigen mimics activate diabetogenic CD8 T cells in NOD mice. Journal of Experimental Medicine 213 (10):2129-2146
  24. Africa JKJJotLSoS (2004) Measuring Biological Diversity. 29 (2):285-286
  25. Graham CH, Letters PVAFJE Phylogenetic beta diversity: linking ecological and evolutionary processes across space in time. 11 (12):1265-1277
  26. Wang J, Shen J, Wu Y, Tu C, Soininen J, Stegen JC, He J, Liu X, Zhang L, Zhang EJIJ Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes. 7 (7):1310-1321
  27. Krause EFJMG (1988) Taxicab Geometry: Adventure in Non-Euclidean Geometry. 72 (461):255
  28. Breiman LJML (2001) Random Forests. 45 (1):5-32. doi:10.1023/a:1010933404324
  29. Strobl C, Boulesteix A-L, Kneib T, Augustin T, Zeileis A (2008) Conditional variable importance for random forests. BMC Bioinformatics 9 (1):307. doi:10.1186/1471-2105-9-307
  30. Metcalf JL, Wegener Parfrey L, Gonzalez A, Lauber CL, Knights D, Ackermann G, Humphrey GC, Gebert MJ, Van Treuren W, Berg-Lyons DJE A microbial clock provides an accurate estimate of the postmortem interval in a mouse model system. 2:e01104---
  31. Knights D, Parfrey LW, Zaneveld J, Lozupone C, Knight RJCh, microbe (2011) Human-associated microbial signatures: examining their predictive value. 10 (4):292-296
  32. Hong H, Xiaoling G, Hua Y Variable selection using Mean Decrease Accuracy and Mean Decrease Gini based on Random Forest. In: 2016 7th IEEE International Conference on Software Engineering and Service Science (ICSESS), 26-28 Aug. 2016 2016. pp 219-224. doi:10.1109/ICSESS.2016.7883053
  33. Zhan G, Yang N, Li S, Huang N, Fang X, Zhang J, Zhu B, Yang L, Yang C, Luo A (2018) Abnormal gut microbiota composition contributes to cognitive dysfunction in SAMP8 mice. Aging 10 (6):1257-1267. doi:10.18632/aging.101464
  34. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto JM, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P, Sicheritz-Ponten T, Turner K, Zhu H, Yu C, Li S, Jian M, Zhou Y, Li Y, Zhang X, Li S, Qin N, Yang H, Wang J, Brunak S, Doré J, Guarner F, Kristiansen K, Pedersen O, Parkhill J, Weissenbach J, Bork P, Ehrlich SD, Wang J (2010) A human gut microbial gene catalogue established by metagenomic sequencing. Nature 464 (7285):59-65. doi:10.1038/nature08821
  35. Rajilić-Stojanović M, de Vos WM (2014) The first 1000 cultured species of the human gastrointestinal microbiota. FEMS Microbiology Reviews 38 (5):996-1047. doi:10.1111/1574-6976.12075
  36. Kanki PJ (2012) Infectious Diseases: Selected Entries from the Encyclopedia of Sustainability Science and Technology. Springer,
  37. Li Z, Yang J, Carslon P, Henning S, Hsu M, Tseng C-H, Thames G, Finegold S, Heber D (2015) Xylooligosaccharide Induced Changes in Gut Microbiota in Healthy and Prediabetic Adults. The FASEB Journal 29 (S1):924.928. doi:10.1096/fasebj.29.1_supplement.924.28
  38. Clavel T, Charrier C, Braune A, Wenning M, Blaut M, Haller D (2009) Isolation of bacteria from the ileal mucosa of TNFdeltaARE mice and description of Enterorhabdus mucosicola gen. nov., sp. nov. International journal of systematic and evolutionary microbiology 59 (Pt 7):1805-1812. doi:10.1099/ijs.0.003087-0
  39. Jones SE, Versalovic J (2009) Probiotic Lactobacillus reuteribiofilms produce antimicrobial and anti-inflammatory factors. BMC Microbiology 9 (1):35. doi:10.1186/1471-2180-9-35
  40. Biagi G, Cipollini I, Pompei A, Zaghini G, Matteuzzi D (2007) Effect of a Lactobacillus animalis strain on composition and metabolism of the intestinal microflora in adult dogs. Veterinary microbiology 124 (1-2):160-165. doi:10.1016/j.vetmic.2007.03.013

Tables

Table 1. The performances of the RF model in validation and testing sets.

Model

Validation set

Testing set

 

RMSE

MAE

R2

RMSE

MAE

R2

PMI prediction model

30.03

20.42

0.945

28.38

20.01

0.947