Learn more: PMC Disclaimer | PMC Copyright Notice
Metabolomic profile of patients with left ventricular assist devices: a pilot study
Abstract
Background
Metabolomic profiling has important diagnostic and prognostic value in heart failure (HF). We investigated whether left ventricular assist device (LVAD) support has an impact on the metabolomic profile of chronic HF patients and if specific metabolic patterns are associated with the development of adverse events.
Methods
We applied untargeted metabolomics to detect and analyze molecules such as amino acids, sugars, fatty acids and other metabolites in plasma samples collected from thirty-three patients implanted with a continuous-flow LVAD. Data were analyzed at baseline, i.e., before implantation of the LVAD, and at long-term follow-up.
Results
Our results reveal significant changes in the metabolomic profile after LVAD implant compared to baseline. In detail, we observed a pre-implant reduction in amino acid metabolism (aminoacyl-tRNA biosynthesis) and increased galactose metabolism, which reversed over the course of support [median follow-up 187 days (63–334 days)]. These changes were associated with improved patient functional capacity driven by LVAD therapy, according to NYHA functional classification of HF (NYHA class I-II: pre-implant =0% of the patients; post-implant =97% of the patients; P<0.001). Moreover, patients who developed adverse thromboembolic events (n=4, 13%) showed a pre-operative metabolomic fingerprint mainly associated with alterations of fatty acid biosynthesis and mitochondrial beta-oxidation of short-chain saturated fatty acids.
Conclusions
Our data provide preliminary evidence that LVAD therapy is associated with changes in the metabolomic profile of HF and suggest the potential use of metabolomics as a new tool to stratify LVAD patients in regard to the risk of adverse events.
Introduction
Myocardial energy substrate metabolism plays an important role in the regulation of cardiac function because continuous energy generation is required for contraction (1). Metabolic remodeling and bioenergetic insufficiency are common features of heart failure (HF), mainly due to impaired mitochondrial function and energy substrate adaptation in the failing heart (1,2). Previous studies reported impaired mitochondrial biogenesis, proteomic remodeling of mitochondrial oxidative pathways, and altered fatty acids and glucose metabolism in HF patients with idiopathic dilated cardiomyopathy or ischemic heart disease (3-6). Moreover, mitochondrial function has been identified as a potential therapeutic target in HF (7).
To date, data regarding the metabolic profile of chronic HF patients who receive mechanical circulatory support with left ventricular assist devices (LVADs) are scarce. Shahinian et al. reported proteome remodeling associated with post-LVAD cardiac remodeling in patients with ischemic cardiomyopathy (8). Weitzel et al. measured different cardiac metabolites in failing hearts post-LVAD support, and showed that while metabolites related to glucose metabolism recovered their values post-LVAD, amino acid levels, creatine, and the citric acid cycle intermediate succinate, did not show a similar pattern (9). Ahmad et al. showed that greater circulating levels of long-chain acylcarnitines in patients with chronic systolic HF were independently predictive of functional status and mortality, and that abnormalities were modifiable with LVAD support (10). These studies contributed to advancing our understanding of the effect of LVAD therapy upon progression of HF and left ventricular biology; however, they did not correlate changes in the metabolic patterns of LVAD patients with clinical outcomes.
Metabolomics is the systematic study of small-molecule metabolite profiles that specific cellular processes leave behind, and it allows us a powerful new approach to the system biology applied to medicine. The term metabolome is derived from the word genome and refers to the complete set of metabolites in an organism and its organelles, or the total complement of metabolites in a cell (11). In this way, metabolomics refers to the use of analytical methods to identify and quantify all metabolites in a biological system, as well as the monitoring of changes in the metabolome of a biofluid or tissue sample following perturbation (11). In parallel with genomics, transcriptomics and proteomics, the application of metabolomic technologies holds the potential to increase our understanding of pathophysiological processes and consequently the potential to identify novel biomarkers to develop new therapeutic strategies (11). Indeed, the identification and quantification of these low molecular weight molecules (e.g., lipids, amino acids, and sugars) might allow for the definition of the phenotype of a given disease. Accordingly, from a clinical perspective, metabolomics might enhance the understanding of the mechanisms underlying a disease, with potential important implications in terms of personalized medicine (11). Metabolomics began to develop in the 1990s with the pioneering work by Lindon, Nicholson and Holmes (12), and has received inputs from different sciences, including informatics, statistics, analytical sciences, chemistry and system biology. Over the years, it has been applied to several medical disciplines, such as neurology, cardiology, and blood diseases (11-13). Metabolic profiling amplifies changes both in the proteome and the genome, and represents an accurate approximation to the phenotype of an organism in health and disease (13).
With regard to LVAD therapy, the outcome of a patient is often strictly dependent on the ability to adopt a patient-tailored therapeutic treatment. Given this, we hypothesized that metabolomics profiling of LVAD patients could further contribute to identifying novel patient-specific therapeutic targets. To test our hypothesis, in this study we applied untargeted metabolomics to: (I) characterize the metabolic profile of chronic HF before LVAD implantation, (II) evaluate the metabolic signature of the clinical effects of LVAD therapy, and (III) identify potential diagnostic fingerprint biomarkers in the patients’ metabolic profile that may indicate patient-specific susceptibility to adverse events.
Methods
Study design
This is a retrospective study conducted in thirty-three patients with end-stage HF implanted with a continuous-flow LVAD at San Raffaele Scientific Institute in Milan from November 2016 to October 2018. Clinical data were recorded at baseline, i.e., before LVAD implantation, and during follow-up visits. Similarly, to evaluate the metabolic profile, patients had blood samples collected prior to the LVAD implantation, forty-eight hours before the intervention, and post-implant, following at least fifty days of LVAD support. Adverse events were diagnosed according to Interagency Registry for Mechanically Assisted Circulatory Support (INTERMACS) criteria.
Approved by the San Raffaele Scientific Institute Institutional Review Board, the study was performed in accordance with the ethical guidelines of the Declaration of Helsinki; all patients provided written informed consent.
Plasma metabolites measurement
We used the Gas Chromatography coupled to Mass Spectrometry analytical technique (GC-MS) to identify and quantify metabolites in platelet-free plasma samples (PFP). Data mining was performed by Multivariate Chemometrics methods. Models from discriminant analysis were calculated and differences in metabolites between pre- vs. post-intervention were evaluated.
The analytical method has been described elsewhere (14). Briefly, 400 µL of PFP were treated with 1,200 µL of cold methanol in 2 mL Eppendorf tubes, vortex mixed, and centrifuged for ten minutes (min) at 14,000 rpm (16,900 G); 400 µL of the upper phase were transferred in glass vials (1.5 mL) and evaporated to dryness overnight; 50 µL of a 0.24 M (20 mg/mL) solution of methoxylamine hydrochloride in pyridine was added to each vial, samples were vortex mixed, and left to react for seventeen hours at room temperature in the dark. Then, 50 µL of N-Methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) were added and left to react for one hour at room temperature. Samples were then diluted with hexane (100 µL) with tetracosane (0.01 mg/mL) as the internal standard, just before GC-MS analysis. Analyses were performed on an Agilent 5977B GC/MS interfaced to the GC 7890B (Agilent Technologies, Palo Alto, CA, USA), equipped with a DB-5ms column (Agilent J&W Scientific, Folsom, CA, USA). Injector temperature was 230 °C, detector temperature 280 °C, helium carrier gas flow rate was set at one mL/min. GC oven temperature program was the following: 90 °C initial temperature, one minute hold time, increasing 10 °C/min to a final temperature of 270 °C, 7 minutes hold time. Samples (1 µL) were injected in split mode (1:4). After a solvent delay of 3 minutes, mass spectra were acquired in full scan mode using 2.28 scans/s with a mass range of 50–700 Amu. Each acquired chromatogram was analyzed by means of the free software AMDIS (Automated Mass spectral Deconvolution and Identification System; http://chemdata.nist.gov/mass-spc/amdis); the software automatically identifies each chromatographic peak by comparison of the relative mass spectra and the retention times with those stored in an in-house library comprising 255 metabolites (15). Other metabolites were identified using NIST14 (National Institute of Standards and Technology’s mass spectral database) and the GOLM Metabolome Database (GMD; http://gmd.mpimp-golm.mpg.de/). Through this approach, 108 compounds were detected and quantified: ninety-seven were accurately identified and eleven compounds were not identified and were defined as unknown.
Statistical analysis
The AMDIS data matrix including 108 metabolites was processed with the SIMCA12plus software. Missing values were replaced with half of the minimum positive values in the original data, and after normalization by sum, data were log transformed and categorized using Pareto scaling for the purposes of analysis, including univariate analysis, partial least square discriminant analysis (PLS-DA), orthogonal PLS-DA (OPLS-DA) (16), and its associated Variable Importance in Projection (VIP) score. PLS-DA models were tested with the leave-one-out cross validation (LOOCV) method for the evaluation of statistical parameters (correlation coefficient R2, cross validation coefficient Q2), which allowed us to determine the optimal number of components for the model description.
Results
Patients’ demographics and pre-operative clinical characteristics are reported in Table 1. Median age at implant was sixty-eight years (IQR: 65–72 years); 31 (94%) patients were male; etiology of HF was ischemic in seventeen patients (51%); seven patients (21%) were on INTERMACS profile one-two, according to INTERMACS classification with TCS-modifier (17); all patients were on NYHA functional class III-IV; intention to treat was destination therapy in thirty patients (91%); ten patients (30%) had diabetes mellitus; thirty patients (91%) were stabilized with short-term mechanical circulatory support prior to LVAD implantation (intra-aortic balloon pump, Impella® or extra-corporeal membrane oxygenation). Fifteen patients (45%) were implanted with the HeartWare HVAD (Medtronic Inc., USA), while eighteen (55%) had the HeartMate3 (Abbott Laboratories, USA). Median time of follow-up at the time of the post-implant analysis was 187 days (IQR: 63–334 days).
Table 1
Variable | Value |
---|---|
Age at implant (year) | 68 [65–72] |
Male sex | 31 (94%) |
Ischemic HF etiology | 17 (51%) |
BMI (kg/m2) | 26.02±5.85 |
INTERMACS profile (*) | |
1 | 4 (12%) |
2 | 3 (9%) |
3 | 18 (55%) |
4 | 8 (24%) |
NYHA functional class | |
I | 0 |
II | 0 |
III | 17 (52%) |
IV | 16 (48%) |
Intention to treat | |
DT | 30 (91%) |
BTT/BTC | 3 (9%) |
Diabetes | 10 (30%) |
Creatinine (mg/dL) | 1.48 (1.16–1.80) |
Hemoglobin (g/dL) | 11.5 (10.3–12.2) |
Temporary MCS | 30 (91%) |
HVAD | 15 (45%) |
HM3 | 18 (55%) |
*, INTERMACS with TCS-modifier. HF, heart failure; BMI, body mass index; INTERMACS, Interagency Registry for Mechanically Assisted Circulatory Support; NYHA, New York Heart Association; DT, destination therapy; BTT, bridge to transplant; BTC, bridge to candidacy; MCS, mechanical circulatory support (Intra Aortic Balloon Pump, Impella, Extra-Corporeal Membrane Oxygenation); HVAD, HeartWare Ventricular Assist Device; HM3, HeartMate3 Ventricular Assist Device.
Figure 1 shows the results of the PLS-DA model. PLS-DA allowed identifying the differences in the metabolomic profile between the pre- vs. post-implant samples; Table 2 summarizes metabolomic differences and changes between the pre- and post-implant samples. According to our results, different important variables (VIP) separating the pre- vs. post-implant samples were identified, belonging to different chemical classes, including: (I) amino acids (alanine, proline and their carboxylation adducts, 4-hydroxyproline, cysteine, threonine, ornithine, glutamine), (II) sugars and reduced sugars (mannose, fructose, glucose, xylitol, sorbitol), (III) carboxylic acids (2-hydroxybutyric, 2-ketobutyric, aminomalonic, gluconic γ-lactone, quinic), and (IV) a purinic derivative. In particular, consistent with the severity of HF, our results show pre-implant reduced amino acid metabolism and increased galactose metabolism, which reversed post-implant (Table 2).
Results of the PLS-DA model discriminating the pre-implant (green dots) vs. post-implant (red dots) metabolomic profile of LVAD patients. Model parameters: R2X=0.198; R2Y=0.819; Q2=0.565. LVAD, left ventricular assist device; PLS-DA, partial least square discriminant analysis.
Table 2
Metabolites | Pre | Post | Chemical family |
---|---|---|---|
Sorbitol | é | ê | Polyol |
Mannose | é | ê | Sugar |
Xylitol | é | ê | Polyol |
2-hydroxybutyric acid | é | ê | Carboxylic acid |
Fructose | é | ê | Sugar |
Benzoic acid | é | ê | Carboxylic acid |
Gluconate acid γ-lactone | é | ê | Carboxylic acid |
Glucose | é | ê | Sugar |
2-ketobutyric acid | é | ê | Carboxylic acid |
Aminomalonic acid | ê | é | Carboxylic acid |
Threonine | ê | é | Amino acid |
Alanine CO2 | ê | é | Amino acid derived |
Proline | ê | é | Amino acid |
Proline CO2 | ê | é | Amino acid derived |
Glutamine | ê | é | Amino acid |
4-hydroxyproline | ê | é | Amino acid |
Cysteine | ê | é | Amino acid |
Alanine | ê | é | Amino acid |
Ornithine | ê | é | Amino acid derived |
Quinic acid | ê | é | Carboxylic acid |
Purine | ê | é | Purine |
Pre, baseline values prior to LVAD implantation; Post, post-implant values, i.e., following LVAD implantation; é: increased; ê: decreased. LVAD, left ventricular assist device; PLS-DA, partial least square discriminant analysis.
Data measured post-implant provide a metabolic signature of the clinical effects of LVAD therapy. These changes are associated with improvement of post-implant patients’ functional status; in detail, we observed a significant shift of patients from NYHA class III-IV to class I-II following LVAD implant (NYHA class I-II: pre-implant =0% of the patients; post-implant =97% of the patients; P<0.001).
Moreover, the comparison of the pre-implant metabolomic profile of patients who developed a thrombotic event (stroke; n=4, 12% of the patients) vs. a matching control population revealed a metabolomic fingerprint in this group mainly associated with alteration of fatty acid biosynthesis and mitochondrial beta-oxidation of short-chain saturated fatty acids (Figure 2), and related to altered profiles of octanoic acid, kynurenine, xylitol, inositol, glycine, pseudouridine, 3,4,5-trihydroxypentanoic acid, aminomalonic acid and gluconic acid γ-lactone.
Results of the OPLS-DA model discriminating the baseline metabolomic fingerprint of LVAD patients who suffered a thromboembolic complication vs. a matching control population. Model parameters: R2X=0.524; R2Y=0.981; Q2=0.854. LVAD, left ventricular assist device; OPLS-DA, orthogonal partial least square discriminant analysis.
Discussion
We investigated the diagnostic value of metabolomics in the setting of durable mechanical circulatory support in a real-world population of patients implanted with a continuous-flow LVAD. Our results show that, beyond producing the starting point of the metabolic profile of chronic HF patients requiring LVAD implantation, metabolomics has allowed characterization of significant changes in this profile over the course of support. In this regard, this study provides, for the first time, preliminary evidence that LVAD therapy triggers significant changes in the metabolomics of HF concomitant to the clinical improvement of functional status (i.e., NYHA class). Moreover, with this study we describe specific metabolic changes driven by an intervention rather than different baseline disorders or disease severity. While it is reasonable to assume that these changes are mainly associated with hemodynamic recovery offered by durable LVAD support, we also speculate that prolonged exposure of patients’ blood to the artificial titanium surface of the LVAD might further impact metabolomics with respect to the baseline pre-operative profile. Similarly, we speculate that prolonged exposure of blood cells to supraphysiological shear stress might also contribute to metabolomic changes over time. In turn, this might expose susceptible patients to higher risk of developing adverse events; this is consistent with previous studies by our group, where we suggested that shear stress triggers and promotes thromboembolic complications in patients with a baseline prothrombotic profile (18-23).
Furthermore, our results suggest that metabolomics have the potential to effectively identify pre-operative diagnostic biomarkers associated with patient-specific susceptibility to develop adverse events, consistent with the recognition of a metabolomic fingerprint in patients who developed a thromboembolic complication (stroke). This peculiar feature of metabolomics might have important implications as far as (I) the identification of new criteria for patients’ risk stratification, and (II) the definition of optimized, patient-tailored post-operative therapeutic strategies are concerned. Of note, our analysis simply required the collection of peripheral blood samples, which is a fast and relatively safe procedure, and supports effective translation of metabolomics in routine clinical practice.
Our results are in accordance with and further expand upon previous studies, which showed that metabolomic analysis has discriminative value in the prognosis and diagnosis of HF (24,25). It is worth noting that the sampling of blood was performed during routine follow-up, before any adverse event occurred, thus reinforcing the prognostic value of our assumption; indeed, our analysis revealed early signals of equilibrium disruption triggering adverse events. Importantly, the baseline “thrombotic fingerprint” eliminates any possible inference of post-implant therapeutic approach, including (I) pharmacological therapy, (II) differences in hemodynamic characteristics of the patients (MAP, aortic insufficiency, right ventricular dysfunction, etc.) and (III) different hemodynamic characteristics of the LVAD support itself (rpm, flow). Our results warrant further studies aimed at exploring the potential of metabolomics to identify further fingerprints prognostic of different scenarios of LVAD adverse events.
Study limitations
This is a single-center retrospective analysis limited by relatively low numbers of patients; in particular, the prognostic capability of metabolomics to stratify patients at risk of developing thromboembolic complications—as well as other scenarios of adverse events—needs to be validated in larger cohort studies. On the other hand, the single-center nature of this study excludes possible biases associated with inhomogeneous post-operative therapeutic approaches. Moreover, we did not analyze differences in the baseline metabolic profile of patients with different etiologies of HF (ischemic vs. non-ischemic). Furthermore, evaluation of the impact of temporary mechanical circulatory support on metabolomic changes warrants further studies. Indeed, we cannot exclude that temporary mechanical circulatory support may have affected the pre-operative metabolic profile. We also acknowledge that post-implant data were collected over a variable duration of follow-up. Finally, the occurrence of specific metabolic changes right before the clinical manifestation of the adverse event was not investigated, and requires future analysis, as it might significantly improve the prognostic power of the method.
Conclusions
This study investigated changes in the metabolomic profile of chronic HF patients implanted with a continuous-flow LVAD. Post-implant, the metabolic profile was altered as compared to baseline, making the metabolic signature a marker of the clinical effects of LVAD therapy. Furthermore, metabolomic profiling allowed for the identification of diagnostic fingerprint biomarkers indicative of patient-specific susceptibility to adverse events, namely thromboembolic complications. Our findings suggest that the analysis of the metabolic profile can complement currently available clinical tools in managing patients implanted with LVADs.
Acknowledgments
Funding: This work was supported by Fondazione Cariplo (Milano, Italy) grant No. 2015-1044 (awarded to FC).
Notes
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
Footnotes
Conflicts of Interest: The authors have no conflicts of interest to declare.