Next Article in Journal
Omics Approaches in Pancreatic Adenocarcinoma
Previous Article in Journal
PD-L1 Blockade by Atezolizumab Downregulates Signaling Pathways Associated with Tumor Growth, Metastasis, and Hypoxia in Human Triple Negative Breast Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 11
Issue 8
10.3390/cancers11081051
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Preliminary Investigation towards the Risk Stratification of Allogeneic Stem Cell Recipients with Respect to the Potential for Development of GVHD via Their Pre-Transplant Plasma Lipid and Metabolic Signature

by
Daniel Contaifer, Jr.
1,
Catherine H. Roberts
2,
Naren Gajenthra Kumar
3,
Ramesh Natarajan
2,
Bernard J. Fisher
2,
Kevin Leslie
4,
Jason Reed
4,
Amir A. Toor
2,5,*,†,‡ and
Dayanjan S. Wijesinghe
1,6,7,*,†,§
1
Department of Pharmacotherapy and Outcomes Sciences, School of Pharmacy, Virginia Commonwealth University, Richmond, VA 23298, USA
2
Department of Internal Medicine, School of Medicine, Virginia Commonwealth University, Richmond, VA 23298, USA
3
Department of Microbiology, School of Medicine, Virginia Commonwealth University, Richmond, VA 23298, USA
4
Department of Physics, College of Humanities and Sciences, Virginia Commonwealth University, Richmond, VA 23284, USA
5
Massey Cancer Center, Virginia Commonwealth University, Richmond, VA 23298, USA
6
Institute for Structural Biology Drug Discovery and Development (ISB3D), VCU School of Pharmacy, Richmond, VA 23298, USA
7
Da Vinci Center, Virginia Commonwealth University, Richmond, VA 23284, USA
*
Authors to whom correspondence should be addressed.
†
Co-senior authors, equal contribution.
‡
Current address: Department of Internal Medicine, Division of Hematology, Oncology & Palliative Care, School of Medicine, Virginia Commonwealth University, P.O. Box 980663, 1250 E Marshall St, Richmond, VA 23298, USA.
§
Current address: Department of Pharmacotherapy and Outcomes Sciences, School of Pharmacy, Virginia Commonwealth University, PO Box 980533, Robert Blackwell Smith Bldg. Rm. 652A, 410 North 12th Street, Richmond, VA 23298-5048, USA.
Cancers 2019, 11(8), 1051; https://doi.org/10.3390/cancers11081051
Submission received: 16 June 2019 / Revised: 12 July 2019 / Accepted: 18 July 2019 / Published: 25 July 2019
Browse Figures
Review Reports Versions Notes

Abstract

:
The clinical outcome of allogeneic hematopoietic stem cell transplantation (SCT) may be influenced by the metabolic status of the recipient following conditioning, which in turn may enable risk stratification with respect to the development of transplant-associated complications such as graft vs. host disease (GVHD). To better understand the impact of the metabolic profile of transplant recipients on post-transplant alloreactivity, we investigated the metabolic signature of 14 patients undergoing myeloablative conditioning followed by either human leukocyte antigen (HLA)-matched related or unrelated donor SCT, or autologous SCT. Blood samples were taken following conditioning and prior to transplant on day 0 and the plasma was comprehensively characterized with respect to its lipidome and metabolome via liquid chromatography/mass spectrometry (LCMS) and gas chromatography/mass spectrometry (GCMS). A pro-inflammatory metabolic profile was observed in patients who eventually developed GVHD. Five potential pre-transplant biomarkers, 2-aminobutyric acid, 1-monopalmitin, diacylglycerols (DG 38:5, DG 38:6), and fatty acid FA 20:1 demonstrated high sensitivity and specificity towards predicting post-transplant GVHD. The resulting predictive model demonstrated an estimated predictive accuracy of risk stratification of 100%, with area under the curve of the ROC of 0.995. The likelihood ratio of 1-monopalmitin (infinity), DG 38:5 (6.0), and DG 38:6 (6.0) also demonstrated that a patient with a positive test result for these biomarkers following conditioning and prior to transplant will be at risk of developing GVHD. Collectively, the data suggest the possibility that pre-transplant metabolic signature may be used for risk stratification of SCT recipients with respect to development of alloreactivity.

1. Introduction

Transplantation of hematopoietic progenitors from an HLA-matched donor is a curative procedure for many patients with hematologic malignancies and disorder of hematopoiesis. Graft vs. host disease (GVHD) is a frequently observed complication of stem cell transplantation (SCT), which contributes to transplant-related mortality and adversely impacts clinical outcomes following transplantation. GVHD after allogeneic hematopoietic stem cell transplantation (HSCT) is a reaction of donor immune cells to recipient tissues. An inflammatory cascade triggered by the preparative regimen causes activated donor T cells to target and destroy epithelial cells. About 35–50% of HSCT recipients will develop acute GVHD. It is mediated by donor-derived T cells responding to minor histocompatibility antigens encountered in the recipient. The T cells encounter these alloantigens, undergo activation, and perform functions such as cytokine secretion (IL-2, IL-4, IL-10, IL-12, IL-17, and interferon gamma by helper T cell subsets) and target lysis (granzyme and perforin secretion by cytotoxic T cells). These functions are accompanied by significant metabolic adaptations in the T cells, including increased glycolysis and oxygen consumption as well as cytokine production [1]. As an example, higher levels of GLUT 1 (glucose transporter 1) expression have been observed in activated T cells, suggesting increased metabolic and biosynthetic rates [2]. Supporting this, correlation has been shown between intracellular ATP concentration in T cells and severity of clinical GVHD in humans, and between increasing glycolysis and GVHD in murine models [3,4]. These observations suggest that donor T cell activation and consequent metabolic and biosynthetic changes may correlate with clinical events as immune reconstitution occurs following SCT, and alloreactivity is triggered.
Just as metabolic changes in the T cell are crucial to the onset of immune reactions, the metabolic milieu in which the T cells find themselves influences their function. In this respect, the lipid molecules constitute a family of important functional mediators. The effects of some lipid molecules on the T cells have been studied recently. These effects include lysophosphatidylserine (lysoPS)-mediated suppression of IL-2 production and suppression of T cell proliferation [5]. This effect is mediated via LPS3/G protein coupled receptor 174, which triggers IL-2 production in CD4+ T cell. Another enzyme, acid sphingomyelinase (ASMase), generates ceramide and modulates signaling cascades involving CD3 and CD28. It is involved in Th1 and Th17 responses through its effect on signal transducer and activator of transcription 3 (STAT 3) and the mammalian target of rapamycin (mTOR) [6]. An acid sphingomyelinase deficient mouse model experiences attenuation of GVHD [7]. Along the same lines, leukotriene C4 has been shown to be important for airway inflammation when administered to murine models along with IL33 [8]. Consistent with such results are the observations that T cells in acetyl-CoA carboxylase deficient mice are resistant to induction of GVHD [9]. T effector cells have been shown to increase their reliance on fatty acid metabolism during GVHD as well [10,11], for example, Prostaglandin E2 (PGE2) has been implicated in modulating T cell effects of mesenchymal stromal cells on T cell populations [12] and PGE2 priming of T cells reduces GVHD when localized to the site of alloreactivity [13]. Further, bone marrow stromal cells also exert an ameliorating influence on GVHD through indoleamine 2,3-dioxygenase (IDO) and PGE2 expression [14]. These lipid-mediated effects have been targeted in the treatment of GVHD, for instance, the effectiveness of leukotriene inhibitor monteleukast has long been recognized in managing GVHD of the lung [15,16]. Additionally, prostaglandins have been studied in GVHD prevention strategies [17], particularly PGE2 [13]. A leukotriene inhibitor, eicosapentanoic acid (EPA), has also been studied in the treatment of GVHD as well as prophylaxis [18,19]. Prostaglandin mediators FT1050 (16,16-dimethyl PGE2) and FT4145 (dexamethasone) are also being studied in clinical trials of GVHD prophylaxis using ex-vivo modification of the allograft [20,21]. These observations make it crucial to gain an understanding regarding the lipid and metabolic changes that come about following dose-intensive myeloablative conditioning, and how the ensuing metabolome and the lipidome might impact T cell function following SCT. Modern methods of lipidomics and metabolomics allow us to study such changes in detail [22,23,24,25,26]. In this paper, we describe the lipidomic and metabolomic profiles of patients undergoing myeloablative conditioning and stem cell transplantation and try to understand the role of these metabolites in mediating alloreactivity, and for potentially predicting GVHD.

2. Materials and Methods

2.1. Patients

Patients were enrolled prospectively in a Virginia Commonwealth University (VCU) Institutional Review Board (IRB)-approved observational pilot study; patients provided written informed consent prior to enrollment (Ethical code: 45 CFR 46.108(b) and 45 CFR 46.109(e) and 45 CFR 46.110 by VCU IRB Panel A, permission date: 27 March 2019). Patients underwent myeloablative conditioning followed by either HLA-matched related (MRD), or HLA-matched unrelated donor (URD); patients undergoing autologous stem cell transplantation (auto) were also included as controls (Table 1). HLA matching was at the allelic level; allogeneic HCT recipients received ATG, calcineurin inhibitors, and either mycophenolate mofetil or methotrexate for GVHD prophylaxis. Blood samples were drawn after completion of myeloablative conditioning therapy on day 0 prior to SCT. Blood was processed for plasma collection and frozen at −80 °C until mass spectroscopy-based metabolomic and lipidomic analysis. Given the small patient cohort, acute and chronic GVHD data were pooled, and Glucksberg and NIH Consensus criteria were used to diagnose and stage GVHD.

2.2. Lipid and Metabolite Extraction for LC-MS/MS Analyses

Blood plasma lipids extraction was carried out using a biphasic solvent system of cold methanol, methyl tertiary butyl ether (MTBE), and water with some modifications (Matyash et al., 2008). In detail, 225 µL of cold methanol containing a mixture of odd chain and deuterated lipid internal standards [lysoPE(17:1), lysoPC(17:0), PC(12:0/13:0), PE(17:0/17:0), PG(17:0/17:0), sphingosine (d17:1), d7 cholesterol, SM(17:0), C17 ceramide, d3 palmitic acid, MG(17:0/0:0/0:0), DG(18:1/2:0/0:0), DG(12:0/12:0/0:0), and d5 TG(17:0/17:1/17:0)] were added to a 20 µL blood plasma aliquot in a 1.5 mL polypropylene tube and then vortexed. Next, 750 µL of cold MTBE was added, followed by vortexing and shaking with an orbital mixer. Phase separation was induced by adding 188 µL of MS-grade water. Upon vortexing (20 s), the sample was centrifuged at 12,300 rpm for 2 min. The upper organic phase was collected in two 300 µL aliquots and evaporated with a vapor trap. Dried extracts were resuspended using 110 µL of a methanol/toluene (9:1, v/v) mixture containing CUDA (50 ng/mL; internal standard for quality control of injection) with support of vortexing (10 s), and centrifuged at 800 rpm for 5 min, followed by transferring 100 µL of the supernatant into auto-sampler vial with an insert. The lower polar layer was collected and an aliquot of 125 µL was evaporated to dryness in a SpeedVac, and resuspended in acetonitrile for polar metabolite analysis via HILIC LC-MS/MS method.

2.3. Metabolomics: GC-MS Metabolite Extraction

A total of 30 μL of plasma sample was added to a 1.0 mL of pre-chilled (−20 °C) extraction solution composed of acetonitrile, isopropanol, and water (3:3:2, v/v/v). Samples were vortexed and shaken for 5 min at 4 °C using the orbital mixing chilling/heating plate. Next, the mixture was centrifuged for 2 min at 14,000 rcf. Then, 450 μL of the supernatant was dried with cold trap concentrator. The dried aliquot was then reconstituted with 450 μL acetonitrile:water (50:50) solution and centrifuged for 2 min at 14,000 rcf. The supernatant was transferred to a polypropylene tube, and subjected to drying in a cold trap. The process of derivatization began with the addition of 10 μL of 40 mg/mL Methoxyamine hydrochloride solution to each dried sample and standard. Samples were shaken at maximum speed at 30 °C for 1.5 h. Then, 91 μL of MSTFA + FAME mixture was added to each sample and standard, and capped immediately. After shaking at maximum speed at 37 °C, the content was transferred to glass vials with micro-inserts inserted and capped immediately.

2.4. Lipids: LC-MC/MC Conditions

Untargeted lipid analysis was undertaken with Sciex TripleTOF 6600 coupled to Agilent 1290 LC. Lipids were separated on an Acquity UPLC CSH C18 column (100 × 2.1 mm; 1.7 µm) (Waters, Milford, MA, USA). The column was maintained at 65 °C and the flow-rate of 0.6 mL/min. The mobile phases consisted of (A) 60:40 (v/v) acetonitrile:water with 10 mM ammonium acetate and (B) 90:10 (v/v) isopropanol:acetonitrile with 10 mM ammonium acetate. The separation was conducted following a stepwise gradient: 0–2 min 15–30% (B), 2–2.5 min 48% (B), 2.5–11 min 82% (B), 11–11.5 min 99% (B), 11.5–12 min 99% (B), 12–12.1 15% (B), and 12–14 min 15% (B). Negative and positive electrospray ionization (ESI) modes were applied with nitrogen serving as the desolvation gas and the collision gas. The parameters for detection of lipids were as follows: Curtain Gas: 35; CAD: High; Ion Spray Voltage: 4500 V; Source Temperature: 350 °C; Gas 1: 60; Gas 2: 60; Declustering Potential: +/−80 V, and collision energies +/−10.

2.5. Metabolites: GC-MS Conditions

A Leco Pegasus IV time of flight mass spectrometer coupled with Agilent 6890 GC equipped with a Gerstel automatic liner exchange system (ALEX) that included a multipurpose sample (MPS2) dual rail, and a Gerstel CIS cold injection system (Gerstel, Muehlheim, Germany) was used to complement HILIC metabolite analysis. The transfer line was maintained at 280 °C. Chromatography separation was achieved on a 30 m long, 0.25 mm i.d. Rtx-5Sil MS column (0.25 μm 95% dimethyl 5% diphenyl polysiloxane film) with the addition of a 10 m integrated guard column was used (Restek, Bellefonte, PA, USA) with helium (99.999%; Airgas, Radnor, PA, USA) at a constant flow of 1 mL/min. The oven temperature was held constant at 50 °C for 1 min and then ramped at 20 °C/min to 330 °C at which it was held constant for 5 min. The GC temperature program was set as follows: 50 °C to 275 °C final temperature at a rate of 12 °C/s and hold for 3 min. The injection volume was 1 μL in splitless mode at 250 °C. Electron impact ionization at 70 V was employed with an ion source temperature of 250 °C. The scan mass ranged from 85 to 500 Da with an acquisition rate of 17 spectra/second.

2.6. Statistical Analysis

Prior to statistical analysis, metabolomic and lipidomic data were subjected to preprocessing as follows. Data were first normalized by a variant of a ‘vector normalization’ by calculating the sum of all peak heights for all identified metabolites for each sample and thereafter normalizing each compound by the total average of the sum. Multivariate statistical analysis tends to focus on metabolites with high intensities. To avoid this tendency, log scaling was applied to reduce the effect of large peaks and scale the data into a more normally distributed pattern. Pareto scaling, which uses the square root of the standard deviation as the scaling to change the emphasis from metabolites with high concentrations to those with moderate or small abundances, was also used for analyzing parameters with large variation.
The patient cohort was divided into patients that developed either acute or chronic GVHD and patients that did not develop GVHD. Because of the small sample size and since the analysis was performed before transplantation, autologous stem cell transplant patients that did not develop GVHD were grouped with allograft recipients that also did not develop GVHD. Partial least squares discriminant analysis (PLS-DA) was used to create a bilinear model to fit the data (Figure 1). Multivariate statistical methods such as PLS-DA have been introduced to reduce the complexity of metabolic spectra and help identify meaningful patterns in high-resolution mass spectrometric data. In this method, the PLS-DA scores can be filtered through calculation of the variable importance in projection (VIP), and used to estimate the contribution of lipids and metabolites for class separation. The top 20 most important variables in the PLS-DA model were selected for further investigation (Figure 1).
Because of the relatively small sample size, cross validation was employed to further evaluate the classification model performance. An algorithm based on support vector machine (SVM) was used to identify the top 20 variable importance in projection (VIP) to further select the best hyperplane that represents the largest separation between the two groups (Figure 1). This method was coupled with Monte-Carlo cross validation (MCCV) through balanced subsampling to support the validation in the small sample cohort (Figure 1). In each MCCV, two thirds of the samples are used to evaluate the feature importance and the remaining one-third is used as test population. The SVM and MCCV allows a multivariate area under the curve of the receiver operating characteristic (AUCROC) analyses to estimate the success of the classification model (Figure 1), creating several AUCROC models to test performances with different numbers of predictors. This procedure was repeated multiple times to calculate the performance estimates and build a confidence interval for each model. Based on the performance in the multivariate AUCROC, the best number of predictors is a reference for selection of potential biomarkers [27].
Univariate AUCROC analyses was used to find potential biomarkers with sufficient power to separate the groups (Figure 1). The result shows the t-test p-value and the AUCROC value with the confidence interval computed using 500 bootstrap replications. The criteria to choose the stronger potential biomarkers for GVHD was to select compounds with the highest AUCROC performance and lowest p-value (p < 0.05). The calculated optimal cutoff was used to estimate the associated sensitivity and specificity values. Positive and negative likelihood ratios are calculated from the sensitivity and specificity output. The data were analyzed using MetaboAnalyst 3.5 maintained by Xia Lab at McGill University. The statistical approach outlined here (Figure 1) was designed to provide meaningful and validated results optimized for the small sample size of the cohort, aiming to find potential biomarkers of GVHD to support future studies. To obtain an adequately powered sample size, patients who did not experience GVHD following an allograft were combined with patients who underwent an autologous SCT in this pilot project.

3. Results

Fourteen patients were analyzed in this study (Table 1). Of these, 10 underwent a myeloablative allograft and 4 underwent an autologous SCT. The entire study cohort was composed of three HLA-matched related donors, seven HLA-matched unrelated donors, and four autologous stem cell transplant recipients, with mean age of 50 (±10) years old; 57% patients were women, and 11 were Caucasians and 3 were African Americans.
Following preprocessing and filtering to remove low-quality data, in the individual patient datasets, the final aggregate, analyzable dataset consisted of 225 plasma lipids and 139 non lipid small-molecule metabolites derived from the patients.

3.1. Pre-Transplant Plasma Lipid and Metabolite Profiles Reveals Class Separation between Those Patients Who Ultimately Developed GVHD and Those Who Did Not

To estimate the potential of the lipids and metabolites to predispose SCT recipients to the development of alloreactivity in the form of either acute or chronic GVHD, post-conditioning and immediate pre-transplant plasma lipid and metabolite data were analyzed via PLS-DA. The degree of separation of patients with future GVHD against patients with no GVHD was visualized by the scores plot of the two principal components. The distance of class separation suggests that metabolic variation may correlate with the development of GVHD in patients undergoing SCT (Figure 2).
The top 20 variables that contributed to the separation observed in the PLS-DA model were represented by 2 metabolites and 18 lipids (Table 2). These characterized the metabolic variation present in the patients, pre-transplant, and correlates strongly with GVHD development post-transplant. Compared to patients with no GVHD, GVHD patients had decreased levels of 2-aminobutyric acid, hexose, monounsaturated fatty acids (14:1, 16:1, 18:1, 19:1, and 20:1), and poly-unsaturated fatty acid (20:3), as well as plasmenyl-ethanolamine (PE(p-34:1) or PE(o-34:2)). Further, GVHD patients presented elevation of monoacylglycerols (1-monoolein and 1-monopalmitin), diacylglycerols (38:5 and 38:6), along with elevated lysophosphocholine (14:0, 20:0), phosphocholines (28:0, 14:0/16:1, 16:0/18:3), and phosphoethanolamines (16:0/18:1, 18:0/22:5).

3.2. The More Important Variables for Class Separation Suggest Metabolic Pathway Tendencies Predispoising to Alloreactivity

So far, our study has identified 20 metabolites, the presence of which correlates with the eventual development of GVHD, suggesting an inherent metabolic disturbance that predisposes a patient towards alloreactivity as early as the day of SCT. Further examination of these metabolites indicates that they modulate three related metabolic pathways. Activated phospholipid metabolism appears to be one of the main alterations associated with GVHD pre-transplant. 1-monopalmitin is a monoacylglycerol that had a high VIP score (2.54), and its elevation in the GVHD group indicates phospholipid degradation in cell membranes to produce diacylglycerol, the precursor for MAGs. Alternatively, the non-GVHD state appears to be associated with elevation of monounsaturated and polyunsaturated fatty acids (Figure 3A). Hexose, more commonly called glucose, also had a high VIP score (2.67), indicating that elevated glucose uptake is on demand for energy production in the tricarboxylic acid cycle (TCA), as well as increased aerobic glycolysis required for hematopoietic cell proliferation. These processes increase the production of the reduced form of nicotinamide adenine dinucleotide (NADH) used in the electron transport chain, but its upregulation induces excessive reducing power that triggers processes such as fatty acid unsaturation and ROS production (Figure 3B). The levels of oxidative stress originated by excessive ROS production is controlled by the glutathione metabolism, where NADH is also used to reduce glutathione for its antioxidant action over ROS. The excessive demand for antioxidative process can deplete glutathione and its precursor cysteine, increasing the demand of 2-aminobutyric acid that can either modulate the glutathione synthesis or be used in ophtalmate synthesis, a tripeptide analog of glutathione, with similar compensative antioxidative actions (Figure 3C).

3.3. The More Important Variables for Class Separation Can Be Used to Build Models for GVHD Association

To evaluate the potential to build models to possibly predict future GVHD, the top 20 highest VIP scores were used in an exploratory analysis. Hence, the analysis showed that plasma metabolites and lipids obtained post-conditioning on day 0 prior to SCT, may be used to build predictive models for GVHD (Figure 4). Models with 2, 3, or 5 variables demonstrated the same level of performance, fitting the selection criteria, with AUCROC ranging from 0.915 to 0.935. With the criteria of finding the lowest number of predictors that can also physiologically explain the metabolic profile of the classes, a model with five predictors was chosen for further exploration as potential biomarkers for risk stratification of patients with the potential for development of GVHD following SCT.

3.4. Univariate ROC Curve Analysis Finds Potential Biomarkers of GVHD With Plasmatic Data Pre-Transplant

The optimal cutoff found for each predictor from the ROC curve was used to estimate the sensitivity and specificity, and to calculate the positive and negative likelihood ratios.
The five best biomarkers and their respective estimates are depicted in Table 3, and the comparison plot for each potential biomarker is presented in Figure 5 where presence of outliers is of notice. These data support the previous finding that a model with five metabolic biomarkers may provide a robust model for GVHD prediction.
All potential biomarkers show high sensitivity and specificity, except for FA 20:1, which, while it has maximum sensitivity, has low specificity. The positive likelihood ratio of 1-monopalmitin and DG (38:5) and DG (38:6) shows that a patient with a positive test result for these biomarkers will have very high odds of developing GVHD. 1-monopalmitin, 2-aminobutyric acid, and DG 38:5 showed CI ranging inside an acceptable AUROC values (0.6–1.0). Yet, the bootstrap confidence interval for DG 36:6 and FA 20:1 are indicative that these two compounds must be taken with caution despite their high sensitivity estimate.
The evaluation of SVM model performance analyzed by AUROC shows the five specific biomarkers model appears to be accurate in predicting the future development of GVHD (AUC = 0.995) with patient’s day 0 plasma drawn pre-transplant (Figure 6A). The predictive power of class probability of the SVM method was also tested and the confusion matrix represented by the probability plot in Figure 6B shows the method is robust without any misclassification.

4. Discussion

The identification of biomarkers that might enable the risk stratification of SCT patients with respect to the potential for development of GVHD has significant clinical utility in the execution of SCT protocols [28,29,30]. We have previously demonstrated that combined lipidomics and metabolomics approached are extremely useful in identifying putative metabolic biomarkers of disease [31]. Here, we demonstrate the discovery of 18 lipids and 2 metabolites separating the GVHD from the non-GVHD cohorts, also implicating specific metabolic pathways involved in the pathogenesis of alloreactivity.
2-aminobutyric acid is a byproduct in the cysteine biosynthesis pathway and relates to glutathione (GSH) metabolism. The presence of an alpha-amino acid as a predictor of GVHD after SCT is substantiated by the reported role of amino acid metabolic changes pre-transplantation [32]. Moreover, Reikvam et al. [33] showed that the altered metabolism of branched chain amino acids, as well as of isobutyryl-carnitine and propyonyl-carnitinine in lipids, before pre-conditioning are associated with the development of acute GVHD. Our results suggest that even as early as after conditioning therapy and before transplantation, the effects of amino acid and lipid metabolism may be associated with GVHD onset.
The metabolic implication of 2-aminobutyric acid in modulation of GSH homeostasis by production of ophthalmic acid is also well known [34]. Ophtalmate is an analog of GSH in which the cysteine group is replaced by L-2-aminobutyrate. It has been proposed that oxidative stress leads to intracellular depletion of GSH, depletion of cysteine, and consequent activation of ophthalmate synthesis [35]. 2-aminobutiric acid increases intracellular GSH levels by regulation of AMP-activated protein kinase and increase of the reduced form of nicotinamide adenine dinucleotide phosphate (NADPH). Moreover, it has been demonstrated that chemotherapy agents are capable of inactivating glutathione reductase, the enzyme that catalyzes the reduction of glutathione disulfide to the sulfhydryl form glutathione [36,37].
Particularly in the immune system, activated T cells undergoing clonal expansion have increased energy demands that increase production of ROS by the mitochondrial electron transport chain [38]. GSH is not necessary for cell activation, but activated T cells regulate their oxidative stress by using GSH, a key component for metabolic reprograming for cell differentiation [39]. Also, TCR ligation and binding with costimulatory molecules induces metabolic remodeling of the naive T cell to anabolic growth and biomass accumulation, and increases aerobic glycolysis [40].
A novel observation from the metabolic profile of GVHD-prone patients in our study is that monoacylglycerol (MAG), diacylglycerol (DAG), fatty acids, phospholipids, and plasmalogens metabolism are significantly altered, identifying lipids as potential mediators of GVHD. Lipid modulation is expected due to metabolic demands of compromised hematopoietic tissue, immunologic response, and underlined inflammatory profile of the patients related to the conditioning regimen. The effects of MUFA and PUFA from the omega-3 family in decreasing inflammation have been extensively studied in efforts to introduce intake of dietary lipids toward treatment of several diseases [41]. These lipids have the potential to decrease production of cytokines in response to LPS and increase the concentration of the anti-inflammatory cytokine IL-10 [42]. The ratio of saturated to monounsaturated fatty acids in membrane phospholipids is critical to normal cellular function. Alterations in this ratio have been correlated with cancer, and the oxidative stress characteristic of the pathology [43,44].
In a study analyzing plasma phospholipids changes in patients with acute leukemia, it was demonstrated that all phospholipids’ concentrations found in patients at the time of diagnosis were significantly lower than in the reference group [45]. Endogenous lipids are important not only in regulation of inflammation, but also in expressing antitumor functions in several types of cancer [46,47,48]. Free fatty acids may be esterified in cell membrane phospholipids undergoing hydrolysis by phospholipases to generate bioactive lipid mediators, including lysophosphocholine (LPC), diacylglycerol (DAG), monoacylglycerol (MAG), and unsaturated fatty acids (MUFA and PUFA). LPC is known to exert immune-regulatory activity by increases in the numbers of T cells, monocytes and neutrophils, and also inducing protein kinase C activation in bone marrow derived mast cells [45].
DAGs are not only a precursor for free fatty acids but also an important signaling molecule in cells. Protein kinase C (PKC) is the major cellular target of DAG. The protein kinase D (PKD) is a substrate of PKC responsible for several cell responses as proliferation, differentiation, apoptosis, and immune response through TCR signaling [49]. Under oxidative stress, ROS induces activation of PKD to protect the cell from oxidative-stress-induced cell death [50].
1-monopalmitin is a MAG formed via release of a fatty acid from DAG by diacylglycerol lipase. The hydrolysis of MAG to FFA and glycerol is conducted by the MAG lipase (MAGL) in different tissues, although ABHD6, a MAG hydrolase, has also been implicated in the pathogenesis of metabolic syndrome [51], inflammation [52], and in cancer [53]. Deletion or inhibition of ABHD6 activity has been shown to be beneficial in certain cancers [54]. The importance of 1-monoacylglycerols with a saturated fatty acid group is demonstrated by its accumulation upon ABHD6 suppression, and direct binding to the ligand binding domain of the peroxisome proliferator-activated receptors PPARα and PPARγ, and activating these transcription factors [55]. MAGs metabolism is also related to the effects of endocannabinoids in the immune system. The endocannabinoid receptor CB2 is identified as a peripheral receptor preferably present in B cells, T cells, macrophages, monocytes, natural killers, and polymorphonuclear. The enzyme MAGL that hydrolyzes MAGs also hydrolyzes 2-arachidonoylglycerol (2-AG), an endogenous endocannabinoid acting through CB2 receptors in the immune system with immune suppression effects [56,57].

5. Limitations

The small size of the cohort used in this study is a limitation that prevents definitive confirmation that the biomarkers identified in this study may be used as predictors to risk stratifying patients with respect to the potential for the development of GVHD. The small sample size also compelled us to combine patients with acute and chronic GVHD, to look at cumulative GVHD incidence and to pool autologous SCT controls with those allograft recipients who did not experience GVHD. Despite these limitations, this study allows an understanding of the metabolic milieu in the study population at the time of transplantation and can be used to direct future studies similar to other studies described in literature [58]. Furthermore, significance following the use of the statistical strategy to overcome the small sample size using Monte Carlo cross validation with bootstrap resampling to create confidence intervals further indicate the study results are of significance. Such approaches have been used previously in the literature as appropriate tools to deal with the small sample size problem [59].

6. Conclusions

Our study demonstrates that the pre-transplant lipidome and the metabolome of SCT recipients has significant potential towards their risk stratification with respect to the development of GVHD, indicating potential use as biomarkers for this purpose. The identified potential biomarkers indicate a pro-inflammatory metabolic profile in patients that will eventually develop GVHD. The role of GSH and its association with 2-aminobutyric acid, coupled with signs of altered glucose metabolism, support the hypothesis that in patients that will develop GVHD, GSH levels are excessively depleted due to elevated oxidative stress related to glucose metabolism, in response to chemotherapy treatment. In this scenario, both groups have elevated 2-aminobutyric acid as a compensatory mechanism, although in patients susceptible to develop GVHD, it is not produced in a large enough amount to compensate for the increased antioxidant demand. Furthermore, the decreased levels of plasma hexose indicate excessive glucose uptake for glycolysis, oxidative phosphorylation, and consequent ROS production. The protective effect of MUFA and PUFA as anti-inflammatory agents is decreased in patients that will develop GVHD, and elevation of phospholipids and DAG and MAG indicates an increased traffic of inflammatory lipids. This pro-inflammatory profile in patients with risk of GVHD is associated with immune suppression, and characterize an unfavorable environment coupled with the overwhelming physiologic impact of undergoing the donor graft, predisposing the host to GVHD.

Author Contributions

Conceptualization, A.A.T. and D.S.W.; data curation, D.C.J., N.G.K., and D.S.W.; formal analysis, D.C.J. and D.S.W.; funding acquisition, D.S.W.; investigation, A.A.T. and D.S.W.; methodology, R.N. and B.J.F.; project administration, C.H.R., A.A.T., and D.S.W.; resources, A.A.T. and D.S.W.; supervision, D.S.W.; validation, A.A.T.; visualization, D.C.J.; writing—original draft, D.C.J.; writing—review and editing, C.H.R., N.G.K., R.N., B.J.F., K.L., J.R., A.A.T., and D.S.W.

Funding

Research reported in this publication was supported by research grants from National Institutes of Health under grant numbers HD087198 (to DSW). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. This work also received support via a Young Investigator Award from SCIEX for clinical lipidomic research (DSW).

Acknowledgments

The authors would like to thank the West Coast Metabolomics Center for helping with the acquisition of some of the data as a fee for service.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jones, N.; Cronin, J.G.; Dolton, G.; Panetti, S.; Schauenburg, A.J.; Galloway, S.A.E.; Sewell, A.K.; Cole, D.K.; Thornton, C.; Francis, N.J. Metabolic Adaptation of Human CD4+ and CD8+ T-Cells to T-Cell Receptor-Mediated Stimulation. Front. Immunol. 2017. [Google Scholar] [CrossRef] [PubMed]
  2. Cretenet, G.; Clerc, I.; Matias, M.; Loisel, S.; Craveiro, M.; Oburoglu, L.; Kinet, S.; Mongellaz, C.; Dardalhon, V.; Taylor, N. Cell surface Glut1 levels distinguish human CD4 and CD8 T lymphocyte subsets with distinct effector functions. Sci. Rep. 2016, 6, 24129. [Google Scholar] [CrossRef] [PubMed]
  3. Martínez, C.; Millán, O.; Rovira, M.P.; Fernández-Avilés, F.; López, A.; Suárez-Lledó, M.; Carreras, E.; Urbano-Ispízua, Á.; Brunet, M. Harmacodynamics of T cell function for monitoring pharmacologic immunosuppression after allogeneic hematopoietic stem cell transplantation. Int. J. Hematol. 2017, 105, 497–505. [Google Scholar] [CrossRef] [PubMed]
  4. Nguyen, H.D.; Chatterjee, S.; Haarberg, K.M.; Wu, Y.; Bastian, D.; Heinrichs, J.; Fu, J.; Daenthanasanmak, A.; Schutt, S.; Shrestha, S.; et al. Metabolic reprogramming of alloantigen-activated T cells after hematopoietic cell transplantation. J. Clin. Investig. 2016, 126, 1337–1352. [Google Scholar] [CrossRef] [PubMed]
  5. Shinjo, Y.; Makide, K.; Satoh, K.; Fukami, F.; Inoue, A.; Kano, K.; Otani, Y.; Ohwada, T.; Aoki, J. Lysophosphatidylserine suppresses IL-2 production in CD4 T cells through LPS3/GPR174. Biochem. Biophys. Res. Commun. 2017, 494, 332–338. [Google Scholar] [CrossRef] [PubMed]
  6. Bai, A.; Guo, Y. Acid sphingomyelinase mediates human CD4+ T-cell signaling: Potential roles in T-cell responses and diseases. Cell Death Dis. 2017, 8, e2963. [Google Scholar] [CrossRef]
  7. Rotolo, J.A.; Stancevic, B.; Lu, S.X.; Zhang, J.; Suh, D.; King, C.G.; Kappel, L.W.; Murphy, G.F.; Liu, C.; Fuks, Z.; et al. Cytolytic T cells induce ceramide-rich platforms in target cell membranes to initiate graft-versus-host disease. Blood 2009, 114, 3693–3706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Lund, S.J.; Portillo, A.; Cavagnero, K.; Baum, R.E.; Naji, L.H.; Badrani, J.H.; Mehta, A.; Croft, M.; Broide, D.H.; Doherty, T.A. Leukotriene C4 Potentiates IL-33–Induced Group 2 Innate Lymphoid Cell Activation and Lung Inflammation. J. Immunol. 2017, 199, 1096–1104. [Google Scholar] [CrossRef] [PubMed]
  9. Raha, S.; Raud, B.; Oberdörfer, L.; Castro, C.N.; Schreder, A.; Freitag, J.; Longerich, T.; Lochner, M.; Sparwasser, T.; Berod, L.; et al. Disruption of de novo fatty acid synthesis via acetyl-CoA carboxylase 1 inhibition prevents acute graft-versus-host disease. Eur. J. Immunol. 2016, 46, 2233–2238. [Google Scholar] [CrossRef]
  10. Chiaranunt, P.; Ferrara, J.L.; Byersdorfer, C.A. Rethinking the paradigm: How comparative studies on fatty acid oxidation inform our understanding of T cell metabolism. Mol. Immunol. 2015, 68, 564–574. [Google Scholar] [CrossRef]
  11. Byersdorfer, C.A.; Tkachev, V.; Opipari, A.W.; Goodell, S.; Swanson, J.; Sandquist, S.; Glick, G.D.; Ferrara, J.L.M. Effector T cells require fatty acid metabolism during murine graft-versus-host disease. Blood 2013, 122, 3230–3237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Najar, M.; Raicevic, G.; Boufker, H.I.; Kazan, H.F.; De Bruyn, C.; Meuleman, N.; Bron, D.; Toungouz, M.; Lagneaux, L. Mesenchymal stromal cells use PGE2 to modulate activation and proliferation of lymphocyte subsets: Combined comparison of adipose tissue, Wharton’s Jelly and bone marrow sources. Cell. Immunol. 2010, 264, 171–179. [Google Scholar] [CrossRef] [PubMed]
  13. Highfill, S.L.; Kelly, R.M.; O’Shaughnessy, M.J.; Zhou, Q.; Xia, L.; Panoskaltsis-Mortari, A.; Taylor, P.A.; Tolar, J.; Blazar, B.R. Multipotent adult progenitor cells can suppress graft-versus-host disease via prostaglandin E2 synthesis and only if localized to sites of allopriming. Blood 2009, 114, 693–701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, M.; Sun, X.; Kuang, X.; Liao, Y.; Li, H.; Luo, D. Mesenchymal stem cells suppress CD8+T cell-mediated activation by suppressing natural killer group 2, member D protein receptor expression and secretion of prostaglandin E2, indoleamine 2, 3-dioxygenase and transforming growth factor-β. Clin. Exp. Immunol. 2014, 178, 516–524. [Google Scholar] [CrossRef] [PubMed]
  15. Or, R.; Gesundheit, B.; Resnick, I.; Bitan, M.; Avraham, A.; Avgil, M.; Sacks, Z.; Shapira, M.Y. Sparing Effect by Montelukast Treatment for Chronic Graft Versus Host Disease: A Pilot Study. Transplantation 2007, 83, 577–581. [Google Scholar] [CrossRef] [PubMed]
  16. Bolanos-Meade, J.; Vogelsang, G.B. Montelukast in Chronic Graft-Versus-Host Disease: A Breath of Fresh Air? Transplantation 2007, 83, 527–528. [Google Scholar] [CrossRef] [PubMed]
  17. Reading, J.L.; Vaes, B.; Hull, C.; Sabbah, S.; Hayday, T.; Wang, N.S.; DiPiero, A.; A Lehman, N.; Taggart, J.M.; Carty, F.; et al. Suppression of IL-7-dependent Effector T-cell Expansion by Multipotent Adult Progenitor Cells and PGE2. Mol. Ther. 2015, 23, 1783–1793. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Takatsuka, H.; Takemoto, Y.; Yamada, S.; Wakae, T.; Mori, A.; Okada, M.; Iwata, N.; Okamoto, T.; Kanamaru, A.; Kakishita, E. Oral eicosapentaenoic acid for acute colonic graft-versus-host disease after bone marrow transplantation. Drugs Under Exp. Clin. Res. 2002, 28, 121–125. [Google Scholar]
  19. Takatsuka, H.; Takemoto, Y.; Iwata, N.; Suehiro, A.; Hamano, T.; Okamoto, T.; Kanamaru, A.; Kakishita, E. Oral eicosapentaenoic acid for complications of bone marrow transplantation. Bone Marrow Transplant. 2001, 28, 769–774. [Google Scholar] [CrossRef] [Green Version]
  20. Cutler, C.S.; Shoemaker, D.; Ballen, K.K.; Robbins, D.; Desponts, C.; Kao, G.S.; Chen, Y.B.A.; Dey, B.R.; McAfee, S.L.; Alyea, E.P.; et al. FT1050 (16-dimethyl Prostaglandin E2)-Enhanced Umbilical Cord Blood Accelerates Hematopoietic Engraftment After Reduced Intensity Conditioning and Double Umbilical Cord Blood Transplantation. Blood 2011, 118, 653. [Google Scholar]
  21. Guerrettaz, L.; Sahaf, N.; Mitchell, L.; Lynn, C.; Raynel, S.; Le, T.; Foster, H.; Truong, C.; Rezner, B.; Robbins, D.; et al. Ex Vivo Modulation of Donor Cells Results in Enhanced Survival and Reduced Gvhd Mortality. Blood 2015, 126, 1884. [Google Scholar]
  22. Kumar, N.G.; Contaifer, D.; Baker, P.R.; Ekroos, K.; Jefferson, K.K.; Wijesinghe, D.S. Untargeted lipidomic analysis to broadly characterize the effects of pathogenic and non-pathogenic staphylococci on mammalian lipids. PLoS ONE 2018, 13, e0206606. [Google Scholar]
  23. Contaifer, D.; Carl, D.E.; Warncke, U.O.; Martin, E.J.; Mohammed, B.M.; Van Tassell, B.; Brophy, D.F.; Chalfant, C.E.; Wijesinghe, D.S. Unsupervised analysis of combined lipid and coagulation data reveal coagulopathy subtypes among dialysis patients. J. Lipid Res. 2016, 58, 586–599. [Google Scholar] [CrossRef] [PubMed]
  24. Jayaraman, S.P.; Anand, R.J.; De Antonio, J.H.; Mangino, M.; Aboutanos, M.B.; Kasirajan, V.; Ivatury, R.R.; Valadka, A.B.; Glushakova, O.; Hayes, R.L.; et al. Metabolomics and Precision Medicine in Trauma: The State of the Field. Shock Augusta Ga 2017, 50, 5–13. [Google Scholar] [CrossRef] [PubMed]
  25. Hoeferlin, L.A.; Huynh, Q.K.; Mietla, J.A.; Sell, S.A.; Tucker, J.; Chalfant, C.E.; Wijesinghe, D.S. The Lipid Portion of Activated Platelet-Rich Plasma Significantly Contributes to Its Wound Healing Properties. Adv. Wound Care 2015, 4, 100–109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wijesinghe, D.S.; Chalfant, C.E. Systems-Level Lipid Analysis Methodologies for Qualitative and Quantitative Investigation of Lipid Signaling Events During Wound Healing. Adv. Wound Care 2013, 2, 538–548. [Google Scholar] [CrossRef]
  27. Xia, J.; Broadhurst, D.I.; Wilson, M.; Wishart, D.S. Translational biomarker discovery in clinical metabolomics: An introductory tutorial. Metabolomics 2013, 9, 280–299. [Google Scholar] [CrossRef]
  28. Ali, A.M.; DiPersio, J.F.; Schroeder, M.A. The role of biomarkers in the diagnosis and risk stratification of Acute Graft vs. Host Disease (aGvHD): A systematic review. Biol. Blood Marrow Transplant. 2016, 22, 1552–1564. [Google Scholar] [CrossRef]
  29. Stikvoort, A.; Chen, Y.; Rådestad, E.; Törlén, J.; Lakshmikanth, T.; Björklund, A.; Mikes, J.; Achour, A.; Gertow, J.; Sundberg, B.; et al. Combining Flow and Mass Cytometry in the Search for Biomarkers in Chronic Graft-versus-Host Disease. Front. Immunol. 2017. [Google Scholar] [CrossRef]
  30. Paczesny, S.; Levine, J.E.; Braun, T.M.; Ferrara, J.L.M. Plasma Biomarkers in Graft-versus-Host Disease: A New Era? Biol. Blood Marrow Transplant. 2009, 15, 33–38. [Google Scholar] [CrossRef] [Green Version]
  31. Contaifer, D., Jr.; Buckley, L.F.; Wohlford, G.; Kumar, N.G.; Morriss, J.M.; Ranasinghe, A.D.; Carbone, S.; Canada, J.M.; Trankle, C.; Abbate, A.; et al. Metabolic modulation predicts heart failure tests performance. PLoS ONE 2019, 14, e0218153. [Google Scholar]
  32. Alborghetti, M.R.; Correa, M.E.P.; Whangbo, J.; Shi, X.; Aricetti, J.A.; Da Silva, A.A.; Miranda, E.C.M.; Sforca, M.L.; Caldana, C.; Gerszten, R.E.; et al. Clinical Metabolomics Identifies Blood Serum Branched Chain Amino Acids as Potential Predictive Biomarkers for Chronic Graft vs. Host Disease. Front. Oncol. 2019. [Google Scholar] [CrossRef] [PubMed]
  33. Reikvam, H.; Hatfield, K.; Bruserud, Ø. The pretransplant systemic metabolic profile reflects a risk of acute graft versus host disease after allogeneic stem cell transplantation. Metabolomics 2016, 12, 12. [Google Scholar] [CrossRef] [PubMed]
  34. Soga, T.; Baran, R.; Suematsu, M.; Ueno, Y.; Ikeda, S.; Sakurakawa, T.; Kakazu, Y.; Ishikawa, T.; Robert, M.; Nishioka, T.; et al. Differential Metabolomics Reveals Ophthalmic Acid as an Oxidative Stress Biomarker Indicating Hepatic Glutathione Consumption. J. Biol. Chem. 2006, 281, 16768–16776. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Irino, Y.; Toh, R.; Nagao, M.; Mori, T.; Honjo, T.; Shinohara, M.; Tsuda, S.; Nakajima, H.; Satomi-Kobayashi, S.; Shinke, T.; et al. 2-Aminobutyric acid modulates glutathione homeostasis in the myocardium. Sci. Rep. 2016, 6, 36749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Frischer, H.; Ahmad, T. Severe generalized glutathione reductase deficiency after antitumor chemotherapy with BCNU [1,3-bis(chloroethyl)-1-nitrosourea]. J. Lab. Clin. Med. 1977, 89, 1080–1091. [Google Scholar] [PubMed]
  37. Babson, J.R.; Reed, D.J. Inactivation of glutathione reductase by 2-chloroethyl nitrosourea-derived isocyanates. Biochem. Biophys. Res. Commun. 1978, 83, 754–762. [Google Scholar] [CrossRef]
  38. Buck, M.D.; Pearce, E.L.; O’Sullivan, D.; O’Sullivan, D. T cell metabolism drives immunity. J. Exp. Med. 2015, 212, 1345–1360. [Google Scholar] [CrossRef] [Green Version]
  39. Mak, T.W.; Grusdat, M.; Duncan, G.S.; Dostert, C.; Nonnenmacher, Y.; Cox, M.; Binsfeld, C.; Hao, Z.; Brüstle, A.; Itsumi, M.; et al. Glutathione Primes T Cell Metabolism for Inflammation. Immunity 2017, 46, 675–689. [Google Scholar] [CrossRef] [Green Version]
  40. Leslie, K.A.; Rasheed, M.; Sabo, R.T.; Roberts, C.C.; Toor, A.A.; Reed, J. Reconstituting donor T cells increase their biomass following hematopoietic stem cell transplantation. Analyst 2018, 143, 2479–2485. [Google Scholar] [CrossRef]
  41. Calder, P.C. Omega-3 polyunsaturated fatty acids and inflammatory processes: Nutrition or pharmacology? Br. J. Clin. Pharmacol. 2013, 75, 645–662. [Google Scholar] [CrossRef] [PubMed]
  42. Calder, P.C. Omega-3 fatty acids and inflammatory processes: From molecules to man. Biochem. Soc. Trans. 2017, 45, 1105–1115. [Google Scholar] [CrossRef] [PubMed]
  43. Ntambi, J.M. Regulation of stearoyl-CoA desaturase by polyunsaturated fatty acids and cholesterol. J. Lipid Res. 1999, 40, 1549–1558. [Google Scholar] [PubMed]
  44. Di Meo, S.; Reed, T.T.; Venditti, P.; Victor, V.M. Role of ROS and RNS Sources in Physiological and Pathological Conditions. Oxidative Med. Cell. Longev. 2016, 2016, 1–44. [Google Scholar] [CrossRef] [PubMed]
  45. Sevastou, I.; Kaffe, E.; Mouratis, M.A.; Aidinis, V. Lysoglycerophospholipids in chronic inflammatory disorders: The PLA2/LPC and ATX/LPA axes. Biochim. Biophys. Acta (BBA)-Mol. Cell Biol. Lipids 2013, 1831, 42–60. [Google Scholar] [CrossRef] [PubMed]
  46. Chiurchiù, V.; Leuti, A.; Maccarrone, M. Bioactive Lipids and Chronic Inflammation: Managing the Fire Within. Front. Immunol. 2018. [Google Scholar] [CrossRef]
  47. Murray, M.; Hraiki, A.; Bebawy, M.; Pazderka, C.; Rawling, T. Anti-tumor activities of lipids and lipid analogues and their development as potential anticancer drugs. Pharmacol. Ther. 2015, 150, 109–128. [Google Scholar] [CrossRef]
  48. Beloribi-Djefaflia, S.; Vasseur, S.; Guillaumond, F. Lipid metabolic reprogramming in cancer cells. Oncogenesis 2016, 5, e189. [Google Scholar] [CrossRef]
  49. Wang, Q.J. PKD at the crossroads of DAG and PKC signaling. Trends Pharmacol. Sci. 2006, 27, 317–323. [Google Scholar] [CrossRef]
  50. Storz, P.; Toker, A. Protein kinase D mediates a stress-induced NF-κB activation and survival pathway. EMBO J. 2003, 22, 109–120. [Google Scholar] [CrossRef]
  51. The Serine Hydrolase ABHD6 Is a Critical Regulator of the Metabolic Syndrome: Cell Reports. Available online: http://www.cell.com/cell-reports/fulltext/S2211-1247(13)00507-X?_returnURL=http%3A%2F%2Flinkinghub.elsevier.com%2Fretrieve%2Fpii%2FS221112471300507X%3Fshowall%3Dtrue (accessed on 6 February 2018).
  52. Alhouayek, M.; Masquelier, J.; Cani, P.D.; Lambert, D.M.; Muccioli, G.G. Implication of the anti-inflammatory bioactive lipid prostaglandin D2-glycerol ester in the control of macrophage activation and inflammation by ABHD6. Proc. Natl. Acad. Sci. 2013, 110, 17558–17563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Li, F.; Fei, X.; Xu, J.; Ji, C. An unannotated α/β hydrolase superfamily member, ABHD6 differentially expressed among cancer cell lines. Mol. Biol. Rep. 2009, 36, 691–696. [Google Scholar] [CrossRef] [PubMed]
  54. Grüner, B.M.; Schulze, C.J.; Yang, D.; Ogasawara, D.; Dix, M.M.; Rogers, Z.N.; Chuang, C.H.; McFarland, C.D.; Chiou, S.H.; Brown, J.M.; et al. An in vivo multiplexed small-molecule screening platform. Nat. Methods 2016, 13, 883–889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Poursharifi, P.; Madiraju, S.R.M.; Prentki, M. Monoacylglycerol signalling and ABHD6 in health and disease. Diabetes Obes. Metab. 2017, 19, 76–89. [Google Scholar] [CrossRef] [Green Version]
  56. Chiurchiù, V.; Battistini, L.; Maccarrone, M. Endocannabinoid signalling in innate and adaptive immunity. Immunology 2015, 144, 352–364. [Google Scholar] [CrossRef] [PubMed]
  57. Cabral, G.A.; Griffin-Thomas, L. Emerging Role of the CB2 Cannabinoid Receptor in Immune Regulation and Therapeutic Prospects. Expert Rev. Mol. Med. 2009, 11, e3. [Google Scholar] [CrossRef]
  58. Bacchetti, P.; Deeks, S.G.; McCune, J.M. Breaking Free of Sample Size Dogma to Perform Innovative Translational Research. Sci. Transl. Med. 2011, 3, 87ps24. [Google Scholar] [CrossRef] [PubMed]
  59. Xi, B.; Gu, H.; Baniasadi, H.; Raftery, D. Statistical Analysis and Modeling of Mass Spectrometry-Based Metabolomics Data. Breast Cancer 2014, 1198, 333–353. [Google Scholar] [Green Version]
Cancers 11 01051 g001 550
Figure 1. Statistical approach used for the identification of potential lipid and metabolite-based biomarkers for the prediction of the alloreactivity following stem cell transplantation (SCT). A supervised statistical approach in the form of partial least squares discriminant analysis (PLS-DA) was applied to the consolidated metabolomic and lipidomic data as a first step to find potential biomarkers and to detect the presence of class separation, if any, between graft vs. host disease (GVHD) and non-GVHD patients. The top 20 most important variables, which separates the groups, were selected by the variable importance in projection (VIP). A support vector machine (SVM) used these VIP variables to further find the binary classification of patients in the two groups. The result was cross-validated with Monte Carlo cross-validation (MCCV) and was used in a multivariate area under the curve of the receiver operating characteristic (AUCROC) analysis to select the best model based on low dimensionality and high accuracy. Model estimates were validated with bootstrap CI. The best model was used as reference to select potential biomarkers of GVHD through univariate analysis of the top 20 VIP that passed the criteria of high AUCROC estimate and low t-test p-value, and validated with bootstrap CI.
Figure 1. Statistical approach used for the identification of potential lipid and metabolite-based biomarkers for the prediction of the alloreactivity following stem cell transplantation (SCT). A supervised statistical approach in the form of partial least squares discriminant analysis (PLS-DA) was applied to the consolidated metabolomic and lipidomic data as a first step to find potential biomarkers and to detect the presence of class separation, if any, between graft vs. host disease (GVHD) and non-GVHD patients. The top 20 most important variables, which separates the groups, were selected by the variable importance in projection (VIP). A support vector machine (SVM) used these VIP variables to further find the binary classification of patients in the two groups. The result was cross-validated with Monte Carlo cross-validation (MCCV) and was used in a multivariate area under the curve of the receiver operating characteristic (AUCROC) analysis to select the best model based on low dimensionality and high accuracy. Model estimates were validated with bootstrap CI. The best model was used as reference to select potential biomarkers of GVHD through univariate analysis of the top 20 VIP that passed the criteria of high AUCROC estimate and low t-test p-value, and validated with bootstrap CI.
Cancers 11 01051 g001
Cancers 11 01051 g002 550
Figure 2. Supervised method PLS-DA identified the direction of maximum covariance between the data et and the class membership, extracting features in the form of latent variables. The model separation of GVHD (YES) and no-GVHD (NO) pre-transplant indicates metabolic differences for prediction of GVHD phenomena with no overlap in the 95% confidence intervals (shaded areas). In the NO group, the allogeneic donor type (allo SCT) is indicated by the red circle, and the autologous donor type (auto SCT) is indicated by pink circles. In the YES group, acute GVHD is indicated by the dark green circle, and the chronic GVHD is indicated by the light green circle.
Figure 2. Supervised method PLS-DA identified the direction of maximum covariance between the data et and the class membership, extracting features in the form of latent variables. The model separation of GVHD (YES) and no-GVHD (NO) pre-transplant indicates metabolic differences for prediction of GVHD phenomena with no overlap in the 95% confidence intervals (shaded areas). In the NO group, the allogeneic donor type (allo SCT) is indicated by the red circle, and the autologous donor type (auto SCT) is indicated by pink circles. In the YES group, acute GVHD is indicated by the dark green circle, and the chronic GVHD is indicated by the light green circle.
Cancers 11 01051 g002
Cancers 11 01051 g003 550
Figure 3. The more important lipids and metabolites associated with GVHD may indicate the main compensatory metabolic pathways modulated in hematopoietic cells pre-transplant. (A) Phospholipids metabolism are affected suggesting Phospholipase C activity to produce diacylglycerols (DAGs), and diacylglycerol lipase and monoacylglycerol lipase activity to produce monoacylglycerols (MAGs) and free fatty acids (FFA), respectively. Predominance of monounsaturated (MUFA) and polyunsaturated (PUFA) free fatty acids are result of increased fatty acyl-CoA desaturases activity, increasing the anticancer activity of ω-3 fatty acids. This metabolic pathway is linked to elevated levels of the reduced form of nicotinamide adenine dinucleotide (NADH), produced in the tri-carboxylic acid cycle (TCA), and important on desaturases activity. (B) Glucose uptake is increased due to increased energy demands in hematopoietic cells, and increased aerobic glycolysis with increased NADH, providing NADH for the electron transport chain and resulting in reactive oxygen species (ROS). (C) Oxidative stress compensation is achieved by activation of cysteine-glutathione pathway, and antioxidant action. Exacerbation of these mechanism leads to depletion of cysteine and glutathione, causing increased 2-aminobutyric acid-ophtalmate compensatory pathway activation originated from amino acids metabolism.
Figure 3. The more important lipids and metabolites associated with GVHD may indicate the main compensatory metabolic pathways modulated in hematopoietic cells pre-transplant. (A) Phospholipids metabolism are affected suggesting Phospholipase C activity to produce diacylglycerols (DAGs), and diacylglycerol lipase and monoacylglycerol lipase activity to produce monoacylglycerols (MAGs) and free fatty acids (FFA), respectively. Predominance of monounsaturated (MUFA) and polyunsaturated (PUFA) free fatty acids are result of increased fatty acyl-CoA desaturases activity, increasing the anticancer activity of ω-3 fatty acids. This metabolic pathway is linked to elevated levels of the reduced form of nicotinamide adenine dinucleotide (NADH), produced in the tri-carboxylic acid cycle (TCA), and important on desaturases activity. (B) Glucose uptake is increased due to increased energy demands in hematopoietic cells, and increased aerobic glycolysis with increased NADH, providing NADH for the electron transport chain and resulting in reactive oxygen species (ROS). (C) Oxidative stress compensation is achieved by activation of cysteine-glutathione pathway, and antioxidant action. Exacerbation of these mechanism leads to depletion of cysteine and glutathione, causing increased 2-aminobutyric acid-ophtalmate compensatory pathway activation originated from amino acids metabolism.
Cancers 11 01051 g003
Cancers 11 01051 g004 550
Figure 4. The top 20 variables from PLS-DA method were used to build appropriate classification models with their respective confidence interval. Predictive models pre-transplant demonstrated robust classification even with only two variables.
Figure 4. The top 20 variables from PLS-DA method were used to build appropriate classification models with their respective confidence interval. Predictive models pre-transplant demonstrated robust classification even with only two variables.
Cancers 11 01051 g004
Cancers 11 01051 g005 550
Figure 5. Comparison of five biomarkers of GVHD (YES) against no GVHD (NO). The comparative analysis revealed elevated 1-monopalmitin, DG (38:5), and DG (38:6) and decreased 2-aminobutyric acid and FA (20:1) in patients that will develop GVHD.
Figure 5. Comparison of five biomarkers of GVHD (YES) against no GVHD (NO). The comparative analysis revealed elevated 1-monopalmitin, DG (38:5), and DG (38:6) and decreased 2-aminobutyric acid and FA (20:1) in patients that will develop GVHD.
Cancers 11 01051 g005
Cancers 11 01051 g006 550
Figure 6. Final model with five specific biomarkers predicting GVHD with cross validation shows strong performance accuracy to predict GVHD (A). The cross-validation for the SVM method for model building and performance evaluation is performed multiple times to test a model with five specific biomarkers and it shows that the model can predict GVHD patients without misclassification (B).
Figure 6. Final model with five specific biomarkers predicting GVHD with cross validation shows strong performance accuracy to predict GVHD (A). The cross-validation for the SVM method for model building and performance evaluation is performed multiple times to test a model with five specific biomarkers and it shows that the model can predict GVHD patients without misclassification (B).
Cancers 11 01051 g006
Table
Table 1. Demographic data for the cohort study.
Table 1. Demographic data for the cohort study.
Donor TypeAgeGenderRaceBMIPrep Regimen;
GVHD Prophylaxis		DiseaseGVHD (Grade/Severity)GVHD Organs/System
MRD60FemaleC21.41Bu/Cy; Cyls/MTXMDSAcute (1)Skin, liver, GI
URD50FemaleC24.42Bu/Cy; Tac/MMFAML *Acute (4)Skin, GI
URD31FemaleC34.99Bu/Cy; Tac/MTXAMLChronic (mod)Skin, GI
URD40FemaleC42.08TBI/Cy; Tac/MTXALLAcute (2)GI
URD57MaleC22.09Bu/Flu; Tac/MMFMDSAcute (2)Skin
URD50MaleC26.13Flu/Mel; Tac/MTXALLAcute (1)Skin
MRD59MaleC21.80Bu/Flu; Cyls/MMFMDSChronic (mod)Oral, GI
URD50FemaleC28.85Bu/Cy; Tac/MTXCML &no-
MRD63FemaleC27.81Bu/Flu; Cyls/MTXET #no-
URD66MaleC29.45Flu/Mel; Tac/MMFMCLno-
auto48FemaleC17.60BEAM; n/aPTCLn/a-
auto35FemaleAA25.04Mel 200; n/aMMn/a-
auto44MaleAA45.22BEAM; n/aNHLn/a-
auto40MaleAA51.22BEAM; n/aALLn/a-
BMI = Body Mass Index; MRD = match related donor; URD = unrelated donor; auto = autologous; C = Caucasian; AA = African American; HTN = hypertension; MI = myocardial infarction; n/a = not applicable. Bu = Busulfan; Cy = Cytoxan; Cyls = cyclosporine; MTX = methotrexate; Tac = tacrolimus; MMF = mycopehnolate mofetil; TBI = total body irradiation; Flu = fludarabine; Mel = melphalan; BEAM = carmustine, etoposide, cytarabine and melphalan; MDS = myelodysplasia; AML = acute myeloid leukemia; HL = Hodgkin lymphoma; CML = chronic myeloid leukemia; ET = essential thrombocythemia; MCL = mantle cell lymphoma; PTCL = Peripheral T-cell lymphoma NOS; MM = multiple myeloma; NHL = non-Hodgkin lymphoma; ALL = acute lymphoblastic leukemia; GI= gastrointestinal system * Secondary to chemotherapy; & blast crisis; # w/progression to myelofibrosis.
Table
Table 2. The metabolites and lipids constituting the top 20 VIP’s that predict class separation on the day of transplant between those patients who eventually went on to develop GVHD as opposed to those who did not develop GVHD following SCT. Patients that went onto develop GVHD had decreased 2-aminobutyric acid, hexose, unsaturated fatty acids, and plasmenyl-ethanolamine PE (p-34:1) or PE (o-34:2), along with elevated monoacylglycerols and diacylglycerols, lysophosphocholines, phosphocholines. and phosphoethanolamines. GVHD’ patients (YES) and non-GVHD (NO) are represented by dark grey (high) or light grey (low), respectively.
Table 2. The metabolites and lipids constituting the top 20 VIP’s that predict class separation on the day of transplant between those patients who eventually went on to develop GVHD as opposed to those who did not develop GVHD following SCT. Patients that went onto develop GVHD had decreased 2-aminobutyric acid, hexose, unsaturated fatty acids, and plasmenyl-ethanolamine PE (p-34:1) or PE (o-34:2), along with elevated monoacylglycerols and diacylglycerols, lysophosphocholines, phosphocholines. and phosphoethanolamines. GVHD’ patients (YES) and non-GVHD (NO) are represented by dark grey (high) or light grey (low), respectively.
ClassPredictorsVIPGVHD
NOYES
Alpha-amino acid2-aminobutyric acid1.80
MonosaccharideHexose2.67
Monoacylglycerol1-monoolein2.11
1-monopalmitin2.54
DiacylglycerolDG 38:51.97
DG 38:61.91
Fatty acidFA 14:11.90
FA 16:12.08
FA 18:11.92
FA 19:12.09
FA 20:12.01
FA 20:32.07
LysophosphatidylcholineLPC 14:02.79
LPC 20:02.46
PhosphatidylcholinePC 28:02.29
PC 14:0/16:11.95
PC 16:0/18:31.95
PhosphatidylethanolaminePE 16:0/18:12.02
PE 18:0/22:52.23
Plasmenyl-ethanolaminePE (p-34:1) or PE (o-34:2)1.75
Table
Table 3. Potential biomarkers performance for a model with top five predictors based on accuracy performance (area under the curve (AUC) estimate).
Table 3. Potential biomarkers performance for a model with top five predictors based on accuracy performance (area under the curve (AUC) estimate).
PredictorAUCp-ValueSensitivitySpecificityLR+LR−
1-monopalmitin0.96
(0.82–1.00)		0.00050.86
(0.71–1.00)		1.00
(1.00–1.00)		infinity0.14
DG 38:50.96
(0.80–1.00)		0.0030.86
(0.57–1.00)		0.86
(0.64–1.00)		6.00.17
DG 38:60.86
(0.57–1.00)		0.0350.86
(0.43–1.00)		0.86
(0.64–1.00)		6.00.17
2-aminobutyric acid0.86
(0.61–1.00)		0.0290.86
(0.57–1.00)		0.71
(0.29–1.00)		3.00.20
FA 20:10.82
(0.49–0.97)		0.0391.00
(1.00–1.00)		0.57
(0.29–0.86)		2.30.00

Share and Cite

MDPI and ACS Style

Contaifer, D., Jr.; Roberts, C.H.; Kumar, N.G.; Natarajan, R.; Fisher, B.J.; Leslie, K.; Reed, J.; Toor, A.A.; Wijesinghe, D.S. A Preliminary Investigation towards the Risk Stratification of Allogeneic Stem Cell Recipients with Respect to the Potential for Development of GVHD via Their Pre-Transplant Plasma Lipid and Metabolic Signature. Cancers 2019, 11, 1051. https://doi.org/10.3390/cancers11081051

AMA Style

Contaifer D Jr., Roberts CH, Kumar NG, Natarajan R, Fisher BJ, Leslie K, Reed J, Toor AA, Wijesinghe DS. A Preliminary Investigation towards the Risk Stratification of Allogeneic Stem Cell Recipients with Respect to the Potential for Development of GVHD via Their Pre-Transplant Plasma Lipid and Metabolic Signature. Cancers. 2019; 11(8):1051. https://doi.org/10.3390/cancers11081051

Chicago/Turabian Style

Contaifer, Daniel, Jr., Catherine H. Roberts, Naren Gajenthra Kumar, Ramesh Natarajan, Bernard J. Fisher, Kevin Leslie, Jason Reed, Amir A. Toor, and Dayanjan S. Wijesinghe. 2019. "A Preliminary Investigation towards the Risk Stratification of Allogeneic Stem Cell Recipients with Respect to the Potential for Development of GVHD via Their Pre-Transplant Plasma Lipid and Metabolic Signature" Cancers 11, no. 8: 1051. https://doi.org/10.3390/cancers11081051

APA Style

Contaifer, D., Jr., Roberts, C. H., Kumar, N. G., Natarajan, R., Fisher, B. J., Leslie, K., Reed, J., Toor, A. A., & Wijesinghe, D. S. (2019). A Preliminary Investigation towards the Risk Stratification of Allogeneic Stem Cell Recipients with Respect to the Potential for Development of GVHD via Their Pre-Transplant Plasma Lipid and Metabolic Signature. Cancers, 11(8), 1051. https://doi.org/10.3390/cancers11081051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop