Mass Spectrometry Based Metabolomics Comparison of Liver Grafts from Donors after Circulatory Death (DCD) and Donors after Brain Death (DBD) Used in Human Orthotopic Liver Transplantation

Abstract

Use of marginal liver grafts, especially those from donors after circulatory death (DCD), has been considered as a solution to organ shortage. Inferior outcomes have been attributed to donor warm ischaemic damage in these DCD organs. Here we sought to profile the metabolic mechanisms underpinning donor warm ischaemia. Non-targeted Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry metabolomics was applied to biopsies of liver grafts from donors after brain death (DBD; n = 27) and DCD (n = 10), both during static cold storage (T₁) as well as post-reperfusion (T₂). Furthermore 6 biopsies from DBD donors prior to the organ donation (T₀) were also profiled. Considering DBD and DCD together, significant metabolic differences were discovered between T₁ and T₂ (688 peaks) that were primarily related to amino acid metabolism, meanwhile T₀ biopsies grouped together with T₂, denoting the distinctively different metabolic activity of the perfused state. Major metabolic differences were discovered between DCD and DBD during cold-phase (T₁) primarily related to glucose, tryptophan and kynurenine metabolism, and in the post-reperfusion phase (T₂) related to amino acid and glutathione metabolism. We propose tryptophan/kynurenine and S-adenosylmethionine as possible biomarkers for the previously established higher graft failure of DCD livers, and conclude that the associated pathways should be targeted in more exhaustive and quantitative investigations.

About the Authors:

Olga Hrydziuszko

Contributed equally to this work with: Olga Hrydziuszko, M. Thamara P. R. Perera

Affiliation: School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom

M. Thamara P. R. Perera

Contributed equally to this work with: Olga Hrydziuszko, M. Thamara P. R. Perera

Affiliation: The Liver Unit, Queen Elizabeth Hospital Birmingham, Edgbaston, Birmingham, B15 2TH, United Kingdom

ORCID http://orcid.org/0000-0002-5417-3850

Richard Laing

Affiliation: The Liver Unit, Queen Elizabeth Hospital Birmingham, Edgbaston, Birmingham, B15 2TH, United Kingdom

Jennifer Kirwan

Affiliation: School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom

Michael A. Silva

Affiliation: The Liver Unit, Queen Elizabeth Hospital Birmingham, Edgbaston, Birmingham, B15 2TH, United Kingdom

Douglas A. Richards

Affiliation: The Department of Pharmacology, School of Medical and Dental Sciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom

Nick Murphy

Affiliation: Department of Critical Care and Anaesthesia, Queen Elizabeth Hospital Birmingham, Edgbaston, Birmingham, B15 2TH, United Kingdom

Darius F. Mirza

Affiliation: School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom

Mark R. Viant

* E-mail: [email protected]

Affiliation: School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, United Kingdom

Introduction

Small (<1,000 Da) metabolites regulate cell signalling, cell-to-cell communication and energy transfer amongst other cellular processes, and are the first biochemicals to respond to internal or external stimuli. This responsiveness makes the metabolome an informative measure of the cell’s dynamic state, a property that has led to a considerable and growing interest in the application of metabolomics in the health sciences. In parallel, metabolomics has begun to be used in clinical solid organ transplantation [1,2] where it holds considerable promise for the discovery of biomarkers to predict poor graft function or patient survival, as well as to elucidate the molecular mechanisms underlying pathophysiological processes such as during graft dysfunction, injury or rejection [3]. Beginning with early studies using ¹H-NMR spectroscopy, the applications of metabolomics in transplantation have expanded in more recent publications to using mass spectrometry based profiling of poor allograft function [4]. Case studies using NMR spectroscopy identified six potential metabolic biomarkers that were distinctive of a non-functional liver[5,6]. Meanwhile a metabolomics assessment in cirrhotic patients by ¹H-NMR spectroscopy linked elevated levels of high-density lipoproteins and phosphocholine with mild chronic liver failure whereas elevated levels of lactate, pyruvate, glucose and creatinine were associated with severe chronic liver failure (CLF). This study provided new insights to hepatic functional impairment in cirrhosis as well as showed an alternative approach to evaluate the severity of CLF, a key aspect in therapeutic decision making [7]. Furthermore, analyses of human bile during liver transplantation by capillary electrophoresis showed distinct metabolic fingerprints in donors and recipients [8]. Whereas more recent publication has been focussed on much larger group of liver grafts belonging to two distinct organ donor sources, namely the cadaveric donation after circulatory death (DCD), and donation after brain death (DBD) and applied lipidomics studies identify markers of early allograft dysfunction [9].

While liver transplantation is well established as the treatment of choice for many indications, the ever growing number of patients listed for transplantation has outweighed the supply of cadaveric organs, leading to greater disparity between supply and demand. As a result, most transplant programmes turned to using marginal organs, in particular liver grafts obtained by donation after circulatory death (DCD), to supplement more traditional donation after brain death (DBD) [10]. Despite rigorous donor selection criteria and best efforts to match a DCD graft with an optimal recipient, the results of the DCD transplantations remain inferior to the standard DBD transplantations and include increased occurrence of primary non-function (PNF) and biliary complications in the immediate post-operative period as well as inferior long-term graft survival rates [11–13]. These complications have been initially associated with donor warm ischaemia (Greek; isch–restriction, hema—blood), typical of DCD grafts during the surgical withdrawal phase. This phase is of variable and unpredictable duration which each DCD donor undergoes once the life support therapy is halted to allow natural passage to circulatory death [14]. Furthermore, grafts obtained from deceased donors are typically stored on ice (at a temperature between 0–4°C) before these organs are transplanted to the recipient, and this phase of cold storage further aggravates graft injury due to ischaemia [15]. Once the blood supply is reconnected within the recipient (both DCD and DBD) to the transplanted organs (termed reperfusion), further organ injury occurs through a process known as ischaemia-reperfusion or preservation-reperfusion injury. In this multifactorial process the reactive oxygen species generated in the organ during warm ischaemia initiate a cellular cascade leading to inflammation, and in severe cases to organ failure termed as PNF [16,17].

Identifying the metabolic differences between DCD and DBD liver grafts could significantly improve current clinical practise by defining biomarkers that are predictive of poor graft function prior to transplantation. Selection of grafts from donors that exhibit such metabolic biomarkers could assist in clinical decision making and the exclusion of those organs from transplantation, thereby preventing the adverse clinical sequelae after transplantation. This would also help those transplant programs that are reliant on cadaveric donor organs for transplantation by preventing the necessity to perform re-transplant operations on those who had failed liver grafts, minimising the burden and demand for organs. Furthermore, identifying the metabolic differences between DCD and DBD liver grafts could help to identify the metabolic modifications of livers prior to and after the organ procurement from the donor that would improve the organ quality, an approach called metabolic therapy. This approach is justified by the nature of current organ donation practice. In the United Kingdom alone there has been a steep rise in DCD donations in recent years, however only 27% of liver grafts from these donors are used in clinical transplantation [18] with approximately 8% of these being excluded due to the high risk of PNF. In addition, the majority of grafts are not even procured due to increased time elapsed between the treatment withdrawal phase and circulatory death, which is beyond the currently accepted criteria of donor warm ischaemia time (30 minutes). Relating the metabolic profiles of DCD and DBD liver grafts to the outcome of the transplantation could supplement and expand the traditional methods to predict organ function, in particular early on during the transplantation procedure.

Previously we reported a pilot study that demonstrated the potential of Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry to detect a few thousand metabolic features (or peaks) in biopsies obtained from liver grafts in the cold and post-reperfusion phases of orthotopic liver transplantation (OLT)[19]. We observed and characterized changes in multiple metabolic pathways showing a rapid resumption of biochemical function within the grafts following reperfusion [19]. Here, we expand considerably upon this initial investigation, specifically with the aim to investigate and characterize the metabolic differences between DCD and DBD liver grafts at two key phases of the liver transplantation, the cold storage phase (T₁) and post-reperfusion phase (T₂). We seek to reveal the underlying metabolic pathways associated with the clinical observation of reduced success of the DCD grafts.

Materials and Methods

This study was approved by the South Birmingham Regional and National Research Ethics committee. The study was aimed at discovery of biomarker differences in the liver grafts used in the clinical orthotopic liver (OLT) setting with possibility of identifying those that may predict poor graft function and primary non-function. Considering the inferior clinical outcomes reported from DCD liver transplantation, we specifically focussed on identifying whether metabolite features differed between DCD and DBD liver grafts. Biopsies were obtained while the grafts were in an ice bath after a variable period of cold ischaemia had elapsed, and again in the post reperfusion phase. In a limited number of patients biopsied were also obtained in the pre-donation setting. All recipients involved were adequately informed and signed a consent form.

Clinical data

Overall, a total of 37 (DBD; n = 27 and DCD; n = 10) liver grafts were studied and the clinical courses of the recipients were followed. In DBD donors, the brain stem death criteria had been confirmed prior to the donor operation, which involved dissection and isolation of graft blood vessels while the circulatory function was intact. In DCD donors (n = 10), death was confirmed according to the Institute of Medicine guidelines after obligatory 5 minutes standoff time from the circulatory arrest, after which the donor operation was performed as rapidly as possible [14]. The technical aspects of the organ procurement were otherwise similar. Both in the DBD and DCD procedures, grafts were perfused with Hyper-Osmolar Citrate (HOS, via cannula inserted to the aorta) and University of Wisconsin (UW, via portal vein) preservation fluids. Following donor operations, the liver grafts were packed on ice and transported to the location of the recipient operation, where the grafts were prepared for implantation while immersed in an ice bath. A wide spectrum of procedure related parameters as well as patient demographics and outcome data were recorded for comparison. These included cold ischaemia time (CIT) elapsed prior to the bench biopsy, overall CIT, donor warm ischaemia time (dWIT) in DCD, immediate post operative outcomes, acute physiological status of the recipient while in the intensive care unit, basic liver graft functions, episodes of graft rejection, graft failures and survival outcomes.

Patients underwent OLT for a variety of indications including alcoholic liver disease, primary biliary cirrhosis, primary sclerosing cholangitis, viral hepatitis (HBV/HCV), non-alcoholic steatohepatitis, polycystic disease, Cryptogenic/Autoimmune and Wilson’s disease. The median age of recipient was 53 and 56 years for DCD and DBD respectively. The median (range) model for end-stage liver disease scores of DCD and DBD graft recipients were 12 (8–22) and 16 (6–22) respectively (Table 1).

[Figure omitted. See PDF.]

Table 1. Donor and recipient demographics and surgical data.

https://doi.org/10.1371/journal.pone.0165884.t001

Liver samples were obtained by Menghini biopsy needle [11520–19, 19swg (1.0mm) x 70mm; Dixons Surgical Instruments Ltd, Wickford, Essex, UK] for all 37 liver grafts at two stages of OLT: T₁ –after organ retrieval and transportation to the implanting centre while the liver was prepared on the bench while immersed in an ice bath at a temperature of 0–4°C T1and T₂ (post reperfusion biopsy) once the graft had been reconnected and the patient haemodynamics stabilised, and usually towards the end of the recipient procedure, after warm ischaemic period and reperfusion injury. In addition, for six DBD grafts an additional biopsy was taken at T₀, during the donor surgical phase while the liver was still in the donor’s body cavity. There were 80 liver allograft biopsies and these were snap frozen in liquid nitrogen and stored at -80°C until sample preparation for the direct infusion FT-ICR mass spectrometry based metabolomics.

Direct infusion FT-ICR mass spectrometry based metabolomics

Samples were prepared for the metabolomics analysis as described previously [19]. Briefly, biopsies were extracted using a methanol:chloroform:water method, separating the extracts into polar and non-polar fractions [20]. In total, 80 samples liver biopsy samples were extracted and from these one quality control (QC) sample was prepared by pooling a fraction of each of the 80 extracts (which was then aliquoted into 11 identical fractions). The polar metabolite fraction of each sample was analysed by USA; LTQ FT Ultra) from m/z 70 to 590, in positive ion mode, using the SIM-stitching approach [21]. Each sample was analysed in triplicate. To minimise false positive metabolites in the data matrix (due to noise), only peaks present in at least 2 of the 3 replicate measurements of each sample were retained, and then only peaks present in at least 75% of all the samples were retained for further analysis [22]. This data processing also served to exclude any peaks in the mass spectra that arose from the drugs that were known to be administrated to the donors and recipients. The final data matrix consisted of 1260 reproducibly detected peaks (rows) and 91 variables (80 biopsies and 11 quality control samples; columns). The matrix contained 9.29% of missing data which was imputed using a weighted k-nearest neighbours algorithm (k = 5) [23]. Data were then normalized using the probabilistic quotient method [24] and subjected to a generalised log transformation (prior multivariate analysis) to stabilise the technical variance across the peaks and hence to avoid the highest abundance peaks from dominating the multivariate analysis [25]. Putative metabolite names were assigned to the peaks based on their mass-to-charge ratio and taking into account commonly detected ions forms, including [M-e]⁺, [M+H]⁺, [M+Na]^+, [M+³⁹K]⁺, [M+⁴¹K]+, [M+2Na-H]⁺, [M+2³⁹K-H]⁺ and [M+NH₄]⁺.

Statistical analyses

Potential clinical differences between DBD and DCD grafts were evaluated by testing each OLT variable, including cold ischaemia time, warm ischaemia time, hours in ITU, number of days-in hospital following OLT, peak aspartate transaminase (AST), and incidence of primary non-function. Non-parametric 2-sample Wilcoxon rank-sum test was applied for continuous numerical variables (e.g., CIT) and Fisher’s exact tests was used for binary variables (e.g., the occurrence of primary non-function). The obtained p values were adjusted for multiple hypothesis testing using the Benjamini and Hochberg method to control the false discovery rate [26].

Principal component analysis (PCA) was used to represent the multivariate FT-ICR mass spectral metabolomics data in 2-dimensional space in terms of principal components PC1 and PC2. Univariate statistical analysis, on a per peak basis, was used to discover if any metabolites differed significantly a) between donor (T₀), cold (T₁) and post-reperfusion (T₂) phases across all patients, and b) between DCD and DBD grafts in the cold phase (T₁) and, separately, post-reperfusion (T₂) phase. Here, the Anderson-Darling test was used to evaluate normality assumptions, and since ca. 40% of the peaks did not follow a normal distribution, non-parametric statistical methods were used, specifically a two-sample Wilcoxon rank-sum test (DCD and DBD comparison) and its extension to more groups, Kruskal-Wallis one-way analysis of variance (T₀, T₁ and T₂ comparison), both with a Benjamini and Hochberg false discovery correction. In addition, Gain Scores Analysis (Kruskal-Wallis on gain scores) was used to discover those metabolites that changed in a significantly different manner from cold-phase to post-reperfusion between the DCD and DBD grafts. FDR of 5% was used as cut-off values to identify statistically significant features. All statistical analyses were carried out using R version 3.0.2, a free programming language and software for statistical computing and graphics. Putatively identified metabolites were assigned to KEGG metabolic pathways, as defined in the KEGG database [27,28]

Results and Discussion

Clinical outcomes

The liver transplantation procedures were carried out in a similar manner for DBD and DCD grafts and we did not notice any significant differences neither at the procedure level nor with the short-term outcomes. The mean CIT was 484.52 ± 143.59 minutes (DBD) and 461.00 ± 116.97 minutes (DCD), whereas the mean implantation time, when grafts were exposed to further warm ischaemia until the circulation was restored, was 41.85 ± 6.94 minutes (DBD) and 41.80 ± 9.77minutes (DCD). The recipients spent on average 95.04 ± 120.80 (DBD) and 123.22 ± 174.48 (DCD) hours in the intensive care unit. The majority of OLTs were successful; three patients in the entire study group had perioperative mortality (n = 2 in the DBD group). The causes of death were related to PNF in two patients (one patient each in the DCD and DBD groups) and related to hepatic artery thrombosis in the third patient (DBD) (Table 1).

Changes in hepatic metabolism during transplantation

FT-ICR mass spectra of the extracted biopsies contained 1260 reproducibly detected peaks of which 448 (35.56%) were putatively annotated based upon accurate mass measurements and the Kyoto Encyclopaedia of Genes and Genomes database (S1 Table). Principal component analysis verified the high technical reproducibility of the mass spectra, evidenced by the clustering of the measurements of the QC sample on the PCA scores plot (Fig 1). Furthermore, the PCA scores showed a clear separation between the biopsies from the cold phase (T₁) and post reperfusion (T₂). The clustering of the donor biopsies (T₀; obtained from six DBD grafts while the organs were still perfused with warm circulation) close to the post-reperfusion biopsies (T₂), with both groups having very distinctive metabolic profiles compared to the biopsies originating from the cold phase sampling is a striking result (Fig 1). This signifies the distinctively different metabolism of hepatocytes in the perfused state compared to those in cold storage. This metabolic separation was confirmed by univariate testing that detected 688 (54.60%) significantly different peaks between T1 and T₂, 293 peaks (23.25%) between T₀ and T1, and only 124 peaks (9.84%) between T₀ and T₂ (Kruskal-Wallis test, p < 0.05) (Table 2). In our previous proof-of-principle OLT study, we identified a plethora of metabolic responses in the post-reperfused grafts compared to their cold-phase state and concluded that these changes reflected the rapid resumption of the biochemical functions of hepatocytes following reperfusion, including increased urea production, bile acid synthesis and clearance of the preservation solution. Here, in addition to verifying these expected metabolic responses, we observed additional key metabolic changes including, amongst others, putatively annotated essential (threonine and valine) and non-essential (tyrosine, serine and proline) amino acids, taurine (a major constituent of bile), and kynurenine (a central compound in the tryptophan metabolism pathway.

[Figure omitted. See PDF.]

Fig 1. Principal component analysis scores plot showing the similarities and differences between the metabolic profiles of the grafts in donor (T₀), cold (T₁) and post-reperfusion phases (T₂).

The variance explained by PC1 = 39.82% and PC2 = 29.58. The close grouping of the donor and post-reperfusion biopsies along with their clear separation from the cold-phase biopsies (along PC1) is further supportive of the rapid resumption of the biochemical functions in the reperfused grafts and shows the direction of metabolic changes through the patient journey and the OLT procedure—from ‘healthy’ donor grafts through cold-phase to almost fully functional grafts post reperfusion.

https://doi.org/10.1371/journal.pone.0165884.g001

[Figure omitted. See PDF.]

Table 2. Top putatively annotated metabolic fold-changes (FC) in the liver grafts between donor phase (T₀), cold phase (T₁) and post-reperfusion phase (T₂), considering the DBD and DCD biopsies as one group.

The average absolute ppm error was 0.3631, range: 0.0029–0.9828.

https://doi.org/10.1371/journal.pone.0165884.t002

Hepatic metabolism in DBD compared to DCD grafts

We identified a small subset of peaks that distinguished the DCD and DBD grafts at the metabolic level. In particular, we detected 50 peaks including 11 putatively annotated compounds that differed between DCD and DBD in the cold phase (T₁), 64 peaks (10 putatively annotated) that differed between DCD and DBD following reperfusion (T₂), and 72 peaks (10 putatively annotated) that changed from T₁ to T₂ in a significantly different manner between DCD and DBD grafts (Table 3). The PCA scores plots, based only on these sub-selections of peaks, provide a visualisation of the clear separation of the DCD vs. DBD grafts along PC1 for the cold-phase biopsies and almost as clear separation for the post-reperfusion biopsies (Fig 2)

[Figure omitted. See PDF.]

Fig 2. Principal component analysis scores plots highlighting the metabolic separation of the DBD and DCD grafts in the cold phase (T₁) and separately the post-reperfusion phase (T₂), based on analyses of just the 50 and 64 peaks identified as being significantly different (between DBD and DCD) for the T₁ and T₂ groups, respectively.

Variance explained for T₁, PC1 = 36.75% and PC2 = 19.39% and for T₂, PC1 = 25.22% and PC2 = 17.41%.

https://doi.org/10.1371/journal.pone.0165884.g002

[Figure omitted. See PDF.]

Table 3. Top putatively annotated metabolic fold-changes between the DCD and DBD grafts: combined results for the comparison, (i) in the cold phase (T₁), (ii) following reperfusion (T₂), and (iii) in response from going from T₁ to T₂.

The average absolute ppm error was 0.2758, range: 0.0069–0.6835.

https://doi.org/10.1371/journal.pone.0165884.t003

The key metabolic differences between DBD and DCD grafts in the cold phase (T₁) included increased levels (in DCD) of the putatively annotated metabolites tryptophan, kynurenine, glucose and leucine and decreased levels of adenylosuccinate, GMP, ADP, malate, O-acetyl carnitine and pantothenate (Table 3). The observed putative metabolites are involved mainly in tryptophan metabolism, purine metabolism, oxidative phosphorylation and a set of carbohydrate metabolic pathways, including the TCA cycle, pyruvate metabolism, glycolysis/gluconeogenesis and the pentose phosphate pathway (Table 4). Among these findings, tryptophan and its metabolism have received earlier attention in the liver transplantation field. Tryptophan is an essential amino acid that, amongst other roles, serves as a precursor of the neurotransmitter serotonin and vitamin B₃. Histidine-tryptophan-ketoglutarate (HTK) solution, which contains tryptophan to prevent membrane injury, was proposed as an alternative liver preservation solution to the current gold standard, University of Wisconsin solution. Interestingly the organ preservation solutions used in this study cohort (UW solution and Hyper-osmolar citrate/Marshall’s solution) do not have added tryptophan. The systematic review to compare the efficacy and safety of these two solutions did not show overall significant differences, yet in some cases HTK was believed to perform better, especially in terms of biliary tract flush and prevention of biliary complications [29]. Furthermore, tryptophan can be catabolised either via the kynurenine or serotonin pathways, and hence kynurenine was studied previously to investigate tryptophan metabolism in potential cirrhotic liver transplant recipients. The pre-transplant serum levels of kynurenine as well as the kynurenine/tryptophan ratios were positively correlated with the disease severity, while serum levels of tryptophan and serotonin showed no correlation[30]. The significantly higher levels of tryptophan and kynurenine in DCD grafts (vs. DBD grafts) during the cold phase in our study appears to support the previous studies that identified tryptophan metabolism via kynurenine pathways as a key metabolic change in liver transplantation. Although it was not a key objective to analyse biomarkers related to primary non-function in the present study given the small sample size, the two failed allografts due to primary non-function both had abundantly higher levels of tryptophan and kynurenine (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Selected top putatively identified metabolites changed in a significantly different way from the cold phase (T₁) to post-reperfusion (T₂) between DCD and DBD grafts.

Gain Scores Analysis (Kruskal-Wallis on gain scores) was used to discover those metabolites that changed in a significantly different manner from cold-phase to post-reperfusion between the DCD and DBD grafts.

https://doi.org/10.1371/journal.pone.0165884.g003

[Figure omitted. See PDF.]

Table 4. Metabolic pathways discovered to differ significantly between DCD and DBD grafts, including the associated putatively annotated metabolites in those pathways.

https://doi.org/10.1371/journal.pone.0165884.t004

The putatively annotated metabolites that significantly differed between DCD and DBD grafts following reperfusion (T₂) included increased levels (in DCD) of glutathione, threonine, leucine, glutamate, creatine, glutamate, proline and decreased levels of pantothenate (Table 2). In addition, four of these metabolites including glutathione, threonine, glutamate and glutamate were changed in a significantly different manner in DCD and DBD grafts while they were removed from cold storage (T₁) and following reperfusion (T₂). The remaining six metabolites identified as different in the Gain Scores Analyses included tryptophan and glucose (previously observed as significantly different in the cold-phase) as well as ornithine, serine, S-adenosyl methionine (SAM) and glucose (Table 2 and Fig 4). All of these putatively annotated metabolites are primarily involved in amino acid metabolism and translation (aminoacyl-tRNA biosynthesis) and to a lesser extent in metabolism of cofactors and vitamins and bile secretion (Table 3). Interestingly, in DBD grafts, there is a large increase in SAM between the cold phase and post-reperfusion, which does not occur in DCD grafts. While the biochemical implications of this lack of recovery are unknown, we hypothesise that this could affect methylation reactions including DNA methylation, given the importance of this metabolite as a methyl donor [31]. Further targeted investigations of the SAM (and related S-adenosylhomocysteine) metabolic pathway are therefore warranted.

[Figure omitted. See PDF.]

Fig 4. The differences between tryptophan and kynurenine in failed allografts due to Primary non-function/PNF (n = 2) vs. non-PNF (n = 36) in the cold phase and post reperfusion.

The data show the relative abundances of the metabolites with 95% confidence intervals (statistics not applied due to limited sample size).

https://doi.org/10.1371/journal.pone.0165884.g004

The increased level of glutathione in DCD grafts, which occurs in the cold phase (1.13 times higher in DCD compared to DBD grafts) but is considerably more pronounced (1.88 times higher) following reperfusion, is also an intriguing finding. Glutathione, owing to its thiol group, is known as one of the most effective antioxidants preventing cellular damage caused by reactive oxygen species, as occurs in ischaemia/reperfusion injury. One of its precursors, acetylcysteine, has been studied as a protective molecule in the perioperative treatment of patients undergoing liver transplantation [32]. Our findings show that not only were glutathione levels increased but so were other intermediates of glutathione metabolism such as glutamate and ornithine, indicative of disruption to the glutathione pathway. While glutathione and glutamate shared similar responses at the post-reperfusion time point, i.e. a rapid increase of levels in DCD compared to DBD, ornithine levels were higher in the cold-phase and similar following reperfusion (Fig 4). This could be due to ornithine being used up for the biosynthesis of glutathione via the intermediate by-product of glutamate. These changes may explain the increased oxidative stress incurred by the DCD grafts, owing to the increased ischaemic injury through warm ischaemic damage.

Although metabolomics studies have caused initial enthusiasm in the transplant community as a potential investigative tool to differentiate poorly functioning grafts from those that perform well after transplantation certain limitations have halted the wider applicability. Only one recent publication proposes a metabolic biosignature, based on combination of metabolites altered in liver allografts that perform unsatisfactorily following liver transplantation[4] and this by far the most explicit application of metabolomics in the field of solid organ transplantation. These metabolites responsible for early allograft dysfunction and graft failure include phospholipid metabolism, bile production, ammonia and urea cycles as well as and glutathione metabolism. Interesting some of the metabolites have also been identified in our study, the key metabolite changes related to tryptophan and kynurenine metabolism have not been identified that could be explained by differences in allograft allograft types and analytical techniques involved.

Conclusions

Previous studies on larger cohorts (n > 300) of DBD and DCD liver transplantations provide evidence for the inferiority of the DCD donations, including increased incidence of primary non-function, biliary complications as well as lower graft- and patient survival [12,33]. Our study has for the first time identified, in un-targeted manner, key metabolic differences between DCD and DBD liver grafts. Since we did not observe any significant differences between the DCD and DBD procedures and outcomes in our study, we believe that these metabolic differences are reflective of the inherent molecular dissimilarities between DCD and DBD grafts. Some of the identified metabolic alterations correlate with our current understanding of the physiological changes surrounding DCD organ donation, including an impact on glucose metabolism by donor warm ischaemia in DCD grafts. However, our observed changes to the tryptophan/kynurenine axis in the DCD grafts are novel findings. Both of these metabolites were observed at ca. 2-fold higher concentration in the DCD grafts (compared to DBD grafts) in the cold phase, suggesting the possibility that these metabolites are responsible for, or at least could be indicators of, the reported higher incidences of increased graft failures in DCD grafts in the literature. In fact, we observed increased levels of tryptophan/kynurenine in allografts with PNF in our study; however, due to few cases we were not able to verify statistically verify this association. Given the role of metabolomics as a hypothesis generating tool, and not to determine whether this metabolic pathway is indeed the cause of graft failure, we conclude that the subsequent clinical investigations of DCD versus DBD transplantations should employ a targeted analytical approach to robustly quantify the metabolites in the tryptophan/kynurenine pathway in the pursuit of more reliable biomarkers of graft function.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Minimal dataset—FTICR mass spectromtry analysis.xlsx.

https://doi.org/10.1371/journal.pone.0165884.s001

(XLSX)

Author Contributions

1. Conceptualization: MTPRP MAS NM DFM MRV.

2. Data curation: OH MTPRP.

3. Formal analysis: OH JK MRV.

4. Funding acquisition: MTPRP MAS DFM.

5. Investigation: OH JK MRV.

6. Methodology: MTPRP OH JK DAR MRV.

7. Project administration: MTPRP.

8. Resources: MTPRP MAS DAR NM DFM MRV.

9. Software: OH JK MRV.

10. Supervision: DFM MRV.

11. Validation: OH JK MRV.

12. Writing – original draft: OH MTPRP.

13. Writing – review & editing: MAS JK DAR NM DFM RL MRV.

Citation: Hrydziuszko O, Perera MTPR, Laing R, Kirwan J, Silva MA, Richards DA, et al. (2016) Mass Spectrometry Based Metabolomics Comparison of Liver Grafts from Donors after Circulatory Death (DCD) and Donors after Brain Death (DBD) Used in Human Orthotopic Liver Transplantation. PLoS ONE 11(11): e0165884. https://doi.org/10.1371/journal.pone.0165884

References

1. Wishart DS (2005) Metabolomics: the principles and potential applications to transplantation. American journal of transplantation 5: 2814–2820. pmid:16302993

2. Bohra R, Klepacki J, Klawitter J, Klawitter J, Thurman JM, Christians U (2013) Proteomics and metabolomics in renal transplantation-quo vadis? Transplant International 26: 225–241. pmid:23350848

3. Naesens M, Sarwal MM (2010) Molecular diagnostics in transplantation. Nature Reviews Nephrology 6: 614–628. pmid:20736923

4. Cortes M, Pareja E, Garcia-Canaveras JC, Donato MT, Montero S, Mir J et al(2014) Metabolomics discloses donor liver biomarkers associated with early allograft dysfunction. Journal of hepatology 61: 564–574. pmid:24798621

5. Serkova NJ, Zhang Y, Coatney JL, Hunter L, Wachs ME, Niemann CU et al.(2007) Early detection of graft failure using the blood metabolic profile of a liver recipient. Transplantation 83: 517. pmid:17318087

6. Serkova NJ, Niemann CU (2006) Pattern recognition and biomarker validation using quantitative 1H-NMR-based metabolomics 6:717–731 pmid:17009906

7. Amathieu R, Nahon P, Triba M, Bouchemal N, Trinchet JC, Beaugrand M et al. (2011) Metabolomic approach by 1H NMR spectroscopy of serum for the assessment of chronic liver failure in patients with cirrhosis. Journal of proteome research 10: 3239–3245. pmid:21568267

8. Papaspyridonos K, Garcia-Perez I, Angulo S, Domann PJ, Vilca-Melendez H, Heaton N et al.(2008) Fingerprinting of human bile during liver transplantation by capillary electrophoresis. Journal of separation science 31: 3058–3064. pmid:18693311

9. Xu J, Casas-Ferreira AM, Ma Y, Sen A, Kim M, Proitsi P et al.(2015) Lipidomics comparing DCD and DBD liver allografts uncovers lysophospholipids elevated in recipients undergoing early allograft dysfunction. Sci Rep 5: 17737. srep17737 [pii]; pmid:26635289

10. Monbaliu D, Pirenne J, Talbot D (2012) Liver transplantation using donation after cardiac death donors. Journal of hepatology 56: 474–485. pmid:21782762

11. Foley DP, Fernandez LA, Leverson G, Chin LT, Krieger N, Cooper JT et al.(2005) Donation after cardiac death: the University of Wisconsin experience with liver transplantation. Annals of surgery 242: 724. pmid:16244547

12. De Vera ME, Lopez-Solis R, Dvorchik I, Campos S, Morris W, Demetris AJ et al.(2009) Liver Transplantation Using Donation After Cardiac Death Donors: Long-Term Follow-Up from a Single Center. American journal of transplantation 9: 773–781. pmid:19344466

13. Grewal HP, Willingham DL, Nguyen J, Hewitt WR, Taner BC, Cornell D et al.(2009) Liver transplantation using controlled donation after cardiac death donors: An analysis of a large single-center experience. Liver Transplantation 15: 1028–1035. pmid:19718636

14. Reich DJ, Mulligan DC, Abt PL, Pruett TL, Abecassis MMI, D'Alessandro A et al.(2009) ASTS recommended practice guidelines for controlled donation after cardiac death organ procurement and transplantation. American Journal of Transplantation 9: 2004–2011. pmid:19624569

15. Mathur AK, Heimbach J, Steffick DE, Sonnenday CJ, Goodrich NP, Merion RM (2010) Donation after cardiac death liver transplantation: predictors of outcome. American Journal of Transplantation 10: 2512–2519. pmid:20977642

16. Totsuka E, Fung JJ, Hakamada K, Ohashi M, Takahashi K, Nakai M et al.(2004) Synergistic effect of cold and warm ischemia time on postoperative graft function and outcome in human liver transplantation. Transplantation proceedings 36: 1955–1958. pmid:15518710

17. Kupiec-Weglinski JW, Busuttil RW (2005) Ischemia and reperfusion injury in liver transplantation. Transplantation proceedings 37: 1653–1656. pmid:15919422

18. Callaghan CJ, Charman SC, Muiesan P, Powell JJ, Gimson AE, van der Meulen JH (2013) Outcomes of transplantation of livers from donation after circulatory death donors in the UK: a cohort study. BMJ open 3: e003287. pmid:24002984

19. Hrydziuszko O, Silva MA, Perera MT, Richards DA, Murphy N, Mirza D et al.(2010) Application of metabolomics to investigate the process of human orthotopic liver transplantation: a proof-of-principle study. OMICS A Journal of Integrative Biology 14: 143–150. pmid:20210660

20. Wu H, Southam AD, Hines A, Viant MR (2008) High-throughput tissue extraction protocol for NMR-and MS-based metabolomics. Analytical biochemistry 372: 204–212. pmid:17963684

21. Southam AD, Payne TG, Cooper HJ, Arvanitis TN, Viant MR (2007) Dynamic range and mass accuracy of wide-scan direct infusion nanoelectrospray fourier transform ion cyclotron resonance mass spectrometry-based metabolomics increased by the spectral stitching method. Anal Chem 79: 4595–4602. pmid:17511421

22. Payne TG, Southam AD, Arvanitis TN, Viant MR (2009) A signal filtering method for improved quantification and noise discrimination in fourier transform ion cyclotron resonance mass spectrometry-based metabolomics data. J Am Soc Mass Spectrom 20: 1087–1095. pmid:19269189

23. Hrydziuszko O, Viant MR (2012) Missing values in mass spectrometry based metabolomics: an undervalued step in the data processing pipeline. Metabolomics 8: 161–174.

24. Dieterle F, Ross A, Schlotterbeck G, Senn H (2006) Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in 1H NMR metabonomics. Analytical Chemistry 78: 4281–4290. pmid:16808434

25. Parsons HM, Ludwig C, Gunther UL, Viant MR (2007) Improved classification accuracy in 1-and 2-dimensional NMR metabolomics data using the variance stabilising generalised logarithm transformation. BMC bioinformatics 8: 234. pmid:17605789

26. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological) 289–300.

27. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M et al.(2008) KEGG for linking genomes to life and the environment. Nucleic acids research 36: D480–D484. pmid:18077471

28. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research 28: 27–30. pmid:10592173

29. Feng L, Zhao N, Yao X, Sun X, Du L, Diao X et al.(2007)Histidine-tryptophan-ketoglutarate solution vs. University of Wisconsin solution for liver transplantation: A systematic review. Liver Transplantation 13: 1125–1136. pmid:17665493

30. Rossouw JE, Labadarios D, Davis M, Williams R (1977) The degradation of tryptophan in severe liver disease. International journal for vitamin and nutrition research 48: 281–289.

31. Mato JM, Lu SC (2007) Role of S-adenosyl-L-methionine in liver health and injury. Hepatology 45: 1306–1312. pmid:17464973

32. Clavien PA, Petrowsky H, DeOliveira ML, Graf R (2007) Strategies for safer liver surgery and partial liver transplantation. New England Journal of Medicine 356: 1545–1559. pmid:17429086

33. Mateo R, Cho Y, Singh G, Stapfer M, Donovan J, Kahn J et al.(2006) Risk factors for graft survival after liver transplantation from donation after cardiac death donors: an analysis of OPTN/UNOS data. American journal of transplantation 6: 791–796. pmid:16539637