It is frequently important to identify the prognosis of fulminant hepatic failure (FHF) patients as this will influence patient management and candidacy for liver transplantation. Therefore, a novel scoring system based on metabonomics combining with multivariate logistic regression was developed to predict the prognosis of FHF mouse model.
BALB/c mice were used to construct FHF model. Parts of plasma were collected at 4, 5, and 6-h time points after treatment, respectively, and detected using gas chromatography/time-of-flight mass spectrometry (GC/TOFMS). The acquired data were processed using partial least square discriminant analysis (PLS-DA). The metabolic markers identified were used to construct a scoring system by multivariate regression analysis.
28 mice of survival group and 28 of dead group were randomly selected and analyzed. PLS regression analysis showed that both the PLS models of 5 h and 6 h after d-galactosamine/lipopolysaccharide treatment demonstrated good performances. Loadings plot suggested that phosphate, beta-hydroxybutyrate (HB), urea, glucose and lactate concentrations in plasma had the highest weightings on the clustering differences at the three time points. By the multivariate logistic regression analysis, the death/survival index (DSI) was constructed based on relative concentrations of HB, urea and phosphate. It provided general accurate rate of prediction of 93.3% in the independent samples.
The novel scoring system based on metabonomics combining with multivariate logistic regression is accurate in predicting the prognosis of FHF mouse model and may be referred in clinical practice as a more useful prognostic tool with other available information.
Fulminant hepatic failure (FHF) is a syndrome characterized by an abrupt onset, hepatic encephalopathy and coagulation abnormality in the absence of preexistent liver disease . The syndrome carries high morbidity and mortality. Orthotopic liver transplantation (OLT) has been increasingly the most important way for FHF, which decreased fatality rate with FHF from 85% to 35% . Unfortunately, the application of OLT has remained low for its higher expense and less donor organs, and a high percentage of patients die before an organ is available . Alternatively, with the earlier use of artificial liver supporting system (ALSS), the nonsurgical medicinal therapy was increasingly improved. Early diagnosis and prognosis identification of FHF patients is likely to result in an intervention at an earlier stage of the disease and, thus, an increased likelihood of treatment success.
In previous studies, individual compounds including lactate, phosphorus among others were proposed to predict patient outcome in FHF [4,5]. However, their prognostic value is debatable . Application of multivariate models has been improved prediction accuracy of FHF outcome. The Clichy criteria developed in patients with hepatitis B consisted of factor V level, patient age, HBsAg and serum alfa-fetoprotein level. The King's College Hospital (KCH) criteria, which was performed based on a retrospective analysis of 585 FHF patients, have been the most widely applied . However, a retrospective study about both criteria in FHF patients reported lower predictive values with the positive predictive values (PPV) at 73% to 79% and the negative predictive values (NPV) at < 50%. The Model for End-Stage Liver Disease (MELD) scoring system is a reliable measure of short-term mortality risk in patients with end-stage chronic liver disease. However, as a predictor of death from FHF, MELD score did not provide more information than the King's College Hospital criteria or international normalized ratio (INR) alone . Therefore, better prognostic model to identify FHF patients who are likely to survive is needed.
After proposed by Nicholson in 1999, metabonomics has increasingly attained attention in recent years, applied to study various human diseases such as coronary artery disease and cancers, and provided new insights into the regulatory investigations [9-11]. One of the most important objectives of metabonomic studies is to identify novel metabolic biomarkers from biofluids in various physiological and pathological conditions. In the previous studies, we reported the dynamics of metabolic profiles in d-galactosamine/lipopolysaccharide (GalN/LPS)-treated BALB/c mouse model, and identified 5 biomarkers including phosphate, β-hydroxybutyrate (HB), urea, glucose and lactate which may constitute a set of markers for the early diagnosis and prognosis of FHF [12,13]. But how to predict the prognosis based on plasma levels of these metabolic biomarkers? In this study, using metabonomic methods based on gas chromatography/time-of-flight mass spectrometry (GC/TOFMS) combined with multivariate logistic regression analysis, we aimed to derive a regression model consisted of several metabolic biomarkers for monitoring FHF, and make a basis of related investigations in clinical practice.
Male BALB/c mice (18-22 g) were purchased from the Academy of Military Medical Sciences (Beijing, China), housed in a standard animal laboratory with constant conditions of temperature, humidity, and 12 h light-dark cycles, and allowed free access to standard laboratory chow and water. Animal studies were carried out in accordance with the Chinese National Research Council guidelines and approved by the Subcommittee on Research Animal Care and Laboratory Animal Resources of the Peking University People's Hospital. Absolute ethanol and chloroform were obtained from Sigma (Sigma-Aldrich Chemical Co. Steinheim, Germany) and used in the extraction of polar metabolites. N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA), methoxyamine hydrochloride and pyridine were obtained from Sigma and used to derivatize pure metabolites. Ribitol was used as an internal reference. Doubly distilled water was used in the preparation of standard and sample solutions.
Construction of experimental FHF model
The FHF model was established as described previously . The model group was intraperitoneally treated with 0.6 g/kg GalN followed by 8 μg/kg LPS, and GalN/LPS treated mice were classed into a survival group and a dead group. Five mice of each group were randomly selected for collection of retro-orbital blood at 4, 5 and 6 h after treatment. Additional samples obtained
Assessment of liver injury
Plasma levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST) and T-bilirubin (TBIL), measured with a 7170A automatic analyzer (Hitachi, Japan), were used to assess the extent of liver injury.
Plasma samples preprocessing and GC/TOFMS analysis
In this study, processing of plasma samples followed a protocol described in ref. . Briefly, the thawed plasma samples were quenched by pure methanol and extracted by chloroform. The upper phase was separated and then evaporated to dryness under a stream of N2 gas. At last, the precipitation was derivatized by methoximation and trimethylsilylation.
The GC/TOFMS system consists of an HP 6890 gas chromatograph and a time-of-flight mass spectrometer (Waters Co., Milford, Massachusetts, USA). As described in previous study , the pretreated plasma samples were analyzed. Identity of GC/MS detected peaks was established by comparing their mass spectra and retention index with those available in the NIST02 library and those of commercially available reference compounds. As for some overlapping peaks, one solution to identify all the peaks in the chromatogram is to use a deconvolution software like AMDIS prior to library searching .
After exported from MarkerLynx, All the collected data from the plasma samples were processed with SIMCA-P plus (version10.0, Umetrics AB, Sweden) as described previously [12,13], in which a range of multivariate statistical analysis were conducted. A data matrix was constructed with the sample IDs as observations and the peaks/retention times as the response variables. It was represented in a K-dimensional space (where K stands for the number of variables), and then projected and reduced to a few principal components that can maximize the separation between classes . Here, partial least square discriminant analysis (PLS-DA) was used to process the required GC/TOFMS data. For PLS analysis, data were visualized by plotting principal component (PC) scores. In the scores, each point represents an individual sample. The plots allow the recognition of clusters of samples with similar scores. Each score plot has a loading plot associated with it, which allows in identifying the spectral regions (metabolites) responsible for the sample clustering observed .
To construct a regression model, multivariate regression analysis was performed based on 56 plasma samples at 5 h after GalN/LPS treatment using the Statistical Package for Social Scientists. The analysis was carried out with the binary response being dead or survived. This model was performed with a model group of 13 randomly selected survival samples and 13 dead samples. A validation group consisted of the remained 30 samples. The linear combination of those biomarkers identified in the metabonomic analysis was then applied to the logistic regression formula to predict the probability of FHF prognosis based on a formula score.
Assessment of liver injury at 5 h and 6 h after treatment
Compared with plasma levels of ALT, AST, and TBIL in the survival group, they were significantly increased in the dead group at 6 h after GalN/LPS treatment. However, plasma ALT, AST, and TBIL levels showed no significant differences between the two groups at 5 h time point (Table 1).
Table 1. Assessment of liver injury at 5 h and 6 h after GalN/LPS treatment
Absence of ribital in plasma samples of BALB/c mice
To confirm absence of ribitol in plasma of BALA/c mice, two samples from the same mouse were analyzed. Ribitol was added to one of them during plasma pretreatment. Figure 1 showed the peak of ribitol in the total ion current (TIC) of the sample with ribitol (upper), and there was no corresponding peak in the same retention time of TIC of another sample without ribitol (lower). Therefore, ribitol can be used as an internal standard to correct for minor variations during sample preparation and analysis.
Figure 1. Absence of ribital in plasma samples of BALB/c mice. The peak of ribitol was showed in TIC of the sample with ribitol (upper), and there was no corresponding peak in the same retention time of TIC of another sample without ribitol (lower).
Partial least square-discriminant analysis
PLS is a supervised statistical method and additional knowledge (Y-variables) must be included to "supervise" the independent data (X-variables). In this study, discriminant analysis was performed using PLS. The PLS diagnostic is tabulated in Table 2 after modeling. Here R2X and R2Y are the fractions of X-variation and Y-variation modeled for specific component, respectively. Q2Y is the estimate of how well the model predicts the Y-variables through a default 7-fold cross-validation. The cumulative values of R2Y and Q2Y of the selected components close to one indicate a perfect fitting. So the diagnostic results indicated both the PLS models of 5 h and 6 h after GalN/LPS treatment demonstrated good performances based on the acquired metabonomics data.
Table 2. Modeling diagnostic of the metabolic data derived from plasma samples of various time points after GalN/LPS treatment
Identification of metabolic markers for distinguishing between the survival and dead groups
In previous studies, clustering analysis of data collected at 4 h, 5 h and 6 h was performed using PLS-DA. Score plots showed distinct clustering differences between the survival and dead groups at 5 h and 6 h after GalN/LPS treatment. Loading plots reveal the impact of metabolites on the clustering results: the closer to zero, the smaller the influence of a linear metabolite combination. All of the three loading plots in Figure 2 suggested that phosphate, HB, urea, glucose and lactate concentrations in plasma had the highest weightings on the clustering differences although there was no clear cluster at 4 h.
Figure 2. Identification of metabolic markers for distinguishing between survival and dead groups. Phosphate, HB, urea, glucose and lactate concentrations in plasma had the highest weightings on the clustering differences at 4 h, 5 h and 6 h after GalN/LPS treatment, although there was no clear cluster at 4 h.
Regression analysis of data at 5 h point
Backward stepwise, which is one of three main approaches to automating the selection of explanatory variables for inclusion in a regression model, was used in the multivariate logistic regression analysis. In backward selection, all variables are included in the initial model, and then removed sequentially until a final model is produced. After 3 steps, lactate and glucose were removed, and other variables (HB, urea and phosphate) were used to calculate a nomogram. The derived nomogram, the death/survival index (DSI) was:
Probability of survival as compared to death was calibrated from the logistic regression analysis using the equation as follows:
The resulting value was defined as the DSI, with DSI < 0 incrementally favored a prognosis of "death", while DSI > 0 corresponded to a higher likelihood of a prognosis of "survival".
Validation of the death/survival index (DSI)
We further validated and calibrated the DSI by comparing the predicted probability of FHF prognosis based on DSI with the validation group. Based on the cut-off value of DSI as above, one of 15 survival samples and one of 15 dead samples were removed, and general accurate rate of prediction was 93.3%. Figure 3 shows predicted probability of survival based on DSI. DSI of >0.65 or < - 0.65 provided a relatively clear diagnosis, corresponding to a 93.3% or 6.7% probability of survival.
Figure 3. Predicted probability of survival based on the death/survival index (DSI). DSI of >0.65 or < - 0.65 provided a relatively clear diagnosis, corresponding to a 93.3% or 6.7% probability of survival. Squares represent dead mice, and diamonds represent surviving mice.