Published in Volume
118, Issue 10
(October 1, 2008)J Clin Invest.
Copyright © 2008, American Society for Clinical Investigation
Metabolite profiling of blood from individuals undergoing planned myocardial infarction reveals early markers of myocardial injury
1Cardiology Division and
2Center for Immunology and Inflammatory Diseases, Massachusetts General Hospital, Boston, Massachusetts, USA.
3Donald W. Reynolds Cardiovascular Clinical Research Center on Atherosclerosis at Harvard Medical School, Boston, Massachusetts, USA.
4Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA.
5Center for Human Genetic Research, Massachusetts General Hospital, Boston, Massachusetts, USA.
6Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, Massachusetts, USA.
7Cardiovascular Institute, Beth Israel Deaconess Medical Center, Boston, Massachusetts, USA.
8Research Group for Clinical Bioinformatics, University for Health Sciences, Medical Informatics and Technology (UMIT), Tirol, Austria.
9Cardiovascular Division, Brigham and Women’s Hospital, Boston, Massachusetts, USA.
Address correspondence to: Robert E. Gerszten, Cardiology Division and Center for Immunology and Inflammatory Diseases, Massachusetts General Hospital, Room 8307, 149 13th Street, Charlestown, Massachusetts 02129, USA. Phone: (617) 724-8322; Fax: (617) 726-1544; E-mail:
First published September 2, 2008
Submitted: January 22,
2008; Accepted: July 9,
Emerging metabolomic tools have created the opportunity to establish metabolic signatures of myocardial injury. We applied a mass spectrometry–based metabolite profiling platform to 36 patients undergoing alcohol septal ablation treatment for hypertrophic obstructive cardiomyopathy, a human model of planned myocardial infarction (PMI). Serial blood samples were obtained before and at various intervals after PMI, with patients undergoing elective diagnostic coronary angiography and patients with spontaneous myocardial infarction (SMI) serving as negative and positive controls, respectively. We identified changes in circulating levels of metabolites participating in pyrimidine metabolism, the tricarboxylic acid cycle and its upstream contributors, and the pentose phosphate pathway. Alterations in levels of multiple metabolites were detected as early as 10 minutes after PMI in an initial derivation group and were validated in a second, independent group of PMI patients. A PMI-derived metabolic signature consisting of aconitic acid, hypoxanthine, trimethylamine N-oxide, and threonine differentiated patients with SMI from those undergoing diagnostic coronary angiography with high accuracy, and coronary sinus sampling distinguished cardiac-derived from peripheral metabolic changes. Our results identify a role for metabolic profiling in the early detection of myocardial injury and suggest that similar approaches may be used for detection or prediction of other disease states.
Novel metabolomics technologies have enhanced the feasibility of acquiring high-throughput snapshots of a whole organism’s metabolic status. Metabolomics studies to date have focused principally on model organisms (1–4). However, the profiling of low–molecular weight biochemicals, including lipids, sugars, and amino acids, that serve as substrates and products in metabolic pathways is particularly relevant to human conditions such as myocardial ischemia (5). In addition to serving as disease biomarkers, circulating metabolites may themselves participate in previously unanticipated roles as regulatory signals with hormone-like functions. In this regard, a recent publication suggests that the tricarboxylic acid (TCA) cycle intermediate succinate, released from ischemic tissues, may activate a novel G protein–coupled receptor in the kidney to modulate blood pressure (6).
However, metabolomics techniques still suffer from signal-to-noise issues, and applications to humans have been limited by the profound degree of interindividual variability. Studies to identify novel disease-related pathways have been further restricted by the inherent unpredictability of the onset of pathological states. Human metabolomics studies have also been complicated by potentially confounding clinical variables such as diet or drug effects, particularly if NMR or mass spectrometry–based (MS-based) profiling techniques are used without unambiguously identifying metabolites associated with peaks of interest. Indeed, published findings of 1H-NMR spectral variation correlating with coronary artery disease severity (7) were later found to be confounded by cholesterol-lowering (i.e., HMG-CoA reductase inhibition) therapy (8). The analysis of samples from large patient cohorts, stratified by known risk factors or exposures, may minimize the impact of clinical confounding variables (8). However, the throughput of many technologies, particularly those that are MS based, currently precludes the analysis of large numbers of samples.
To help circumvent these problems, we applied a targeted MS-based metabolomics platform, which provides high analyte specificity (9, 10), to serial blood samples obtained from carefully phenotyped patients undergoing alcohol septal ablation for hypertrophic obstructive cardiomyopathy (HOCM). This human model of planned myocardial infarction (PMI; refs. 11, 12) faithfully reproduces spontaneous myocardial infarction (SMI; refs. 13–15). The serial sampling before and after the controlled myocardial insult allows each patient to serve as his or her own biological control. To assess whether metabolic alterations were localized to the heart after myocardial injury, we simultaneously sampled the venous effluent from the heart as well as peripheral blood. The specificity of metabolic alterations was also verified by studying patients undergoing catheterization without the induction of myocardial infarction (MI). Finally, to assess generalizability, we tested our metabolic biomarkers in patients with SMI.
PMI recapitulates SMI. Clinical characteristics of the HOCM patients and other study populations are detailed in Table 1. The septal ablation recapitulated important clinical features of SMI, including substernal chest pain and electrocardiographic changes as well as development of echocardiographic evidence of septal wall motion abnormalities, as previously described by our group and others (13–15). The standard biochemical metrics of myocardial injury, creatine kinase (CK) muscle and brain subunits (CK-MB) and troponin T, were within normal limits prior to septal ablation and increased to 200 ± 98 ng/ml and 4.5 ± 2.6 ng/ml, respectively. CK-MB peaked at 8.9 ± 4.5 hours and cardiac troponin T at 14.9 ± 8.0 hours after PMI, time courses consistent with SMI (16).
Baseline clinical characteristics of study cohorts
Derivation and validation of early metabolic changes in peripheral plasma. We first studied peripheral blood samples across the range of time points available (10 minutes to 24 hours) to characterize metabolic alterations associated with PMI in a derivation cohort of 20 patients. Table 2 describes the most significantly changed metabolites 10 minutes after the onset of myocardial injury. We observed a total of 7 metabolites that were significantly changed (nominal P < 0.005), with an estimated false discovery rate of approximately 11% (see Methods). Based on this false discovery rate, we expected that 6 of the 7 metabolites were truly differentially changed.
Metabolite changes in the peripheral blood 10 minutes after myocardial injury
Because biomarker discovery studies are vulnerable to unintentional overfitting of data (17, 18), we then assessed the performance of metabolites that were significantly changed in the derivation cohort in an independent validation cohort of patients also undergoing the septal ablation procedure (n = 16; Table 2). Significant changes (P < 0.05) were noted in the validation cohort for 6 of the 7 metabolites observed in the derivation cohort. These metabolites all showed concordance with the derivation cohort in the direction of change, and the magnitude of change in the 2 cohorts was highly correlated (r2 = 0.87, P = 0.01). Metabolic changes included products of purine and pyrimidine catabolism (hypoxanthine and malonic acid, respectively) as well as several amino acids. Of note, the alterations in these metabolites were seen when no significant rises in the clinically available biomarkers, CK-MB and troponin T, were detectable in the plasma (P = NS for both).
By 60 minutes after PMI, more metabolic changes were noted (Table 3). All of the metabolic changes documented at 10 minutes were also observed at 60 minutes, underscoring the consistency of our findings. At 60 minutes, 13 metabolites were significantly changed (nominal P < 0.005), with an estimated false discovery rate of 6.5%. Significant changes in 10 of the 13 metabolites observed in the derivation cohort were noted in the validation cohort as well (P < 0.05), with strong trends for the remaining 3 metabolites. Again, the magnitude of change in the 2 cohorts was highly correlated at this time point (r2 = 0.93, P < 0.0001). By 60 minutes, additional changes in purine metabolites (xanthine and inosine), an inflammatory mediator (methylhistamine), and other amines and amino acids were detected.
Metabolite changes in the peripheral blood 60 minutes after myocardial injury
We explored the specificity of the findings observed in the PMI cohort by examining blood samples from patients undergoing routine cardiac catheterization, without the induction of MI that occurs in the unique ablation injury model. Whereas the metabolic changes in the derivation and validation cohorts were highly correlated, as noted above, there was no correlation between the derivation cohort and the catheterization control group overall (P = 0.47 at 10 minutes; P = 0.76 at 60 minutes). However, cardiac catheterization alone was associated with changes in 3 metabolites, tryptophan, tyrosine, and phenylalanine, at either the 10-minute or the 60-minute time point (P < 0.01). These 3 metabolites were therefore not included in Tables 2 and 3 and were excluded from further analysis. Thus, metabolites with changes that were not specific to myocardial injury and that may instead reflect procedural events such as arteriotomy or catheter manipulation were eliminated using the appropriate patient controls.
Kinetic analysis of metabolic changes in peripheral plasma. Figure 1 shows representative metabolites across a spectrum of time points after PMI. Some metabolic changes were relatively transient (≤120 minutes), as seen with malonic acid and alanine, whereas other early metabolic changes persisted (≥240 minutes), as seen with hypoxanthine, aconitic acid, trimethylamine N-oxide (TMNO), and threonine. Later-appearing metabolites included anthranilic acid and creatine, the latter of which was detected in a time course in our study consistent with that of CK-MB release from cardiomyocytes undergoing necrosis (19).
Kinetics of metabolic changes in peripheral human plasma after planned myocardial injury. Representative metabolites demonstrating transient (≤120 minutes after injury; see Table 2), sustained (≥240 minutes after injury), and late (1,440 minutes after injury, i.e., 1 day) changes in plasma levels compared with baseline are shown. Data represent the entire cohort of 36 patients: boxes denote interquartile range, lines denote median, “+” denotes mean, whiskers denote fifth and ninety-fifth percentiles, and outliers are shown by points.
Early metabolic changes in coronary sinus plasma. We further examined the specificity of our findings by exploring the anatomic origin of the early metabolic changes. In a subgroup of 13 patients, we compared metabolite levels in samples obtained simultaneously from the peripheral blood and from a catheter placed in the coronary sinus, the venous outflow of the heart. This simultaneous sampling allowed us to identify transmyocardial enrichment of metabolites. As a reference for comparison, B-type natriuretic peptide, a protein released from the heart in response to left ventricular wall stress, was enriched approximately 1.3-fold in the coronary sinus samples 60 minutes after injury (data not shown). We assessed whether the metabolites that increased at 10 or 60 minutes in the peripheral blood (1-methylhistamine, aminoisobutyric acid, hypoxanthine, inosine, isoleucine/leucine, malonic acid, and xanthine; Tables 2 and 3) were also enriched in the coronary sinus at 10 minutes (Table 4). We evaluated relative changes in metabolite levels in the coronary sinus compared with the periphery for the rest of the metabolites in our platform as well. Additional biomarkers that showed greater than 1.3-fold transmyocardial enrichment 10 minutes after PMI are shown in Table 5. Peripheral and coronary data for the 60-minute time point are shown in Tables 6 and 7. All of the metabolites that increased in the peripheral samples also increased in the coronary sinus samples, underscoring the consistency of our findings. Of the 7 metabolites that increased in the peripheral blood, 5 showed at least a 1.3-fold enrichment in the coronary sinus at 10 minutes after myocardial injury. As expected, a transmyocardial enrichment pattern was evident for metabolites related to myocardial anaerobic metabolism, including lactic acid and succinic acid, as well as ATP degradation products, such as hypoxanthine and AMP. Furthermore, changes in the levels of certain metabolites were first apparent in the coronary sinus, and were only later detected in the peripheral samples. For example, inosine and 1-methylhistamine levels were significantly elevated by 10 minutes in the coronary sinus, but were not significantly elevated until 60 minutes in the periphery (Tables 3 and 4). A similar pattern was observed for malic acid and glycerol-3-phosphate (Figure 2). The pattern of early enrichment in the coronary sinus with later detection in the periphery supports the cardiac origin of these observed metabolic changes. Coronary sinus samples unmasked some metabolic changes that were not statistically significantly changed in peripheral plasma, such as with taurine (Figure 2) as well as succinic acid, orotic acid, and ribose-5-phosphate (data not shown), perhaps because of overall limited sample numbers, rapid catabolism, excretion, or wide distribution — and consequent dilution — once circulated. Finally, other metabolites, such as glutamic acid (Figure 2), were concomitantly elevated in both the peripheral and coronary sinus samples at 10 minutes, although to a greater degree in the coronary sinus.
Kinetic analyses of representative metabolites that are enriched in the coronary sinus after myocardial injury. Data (mean ± SEM) are restricted to the 13 patients who had simultaneous peripheral blood and coronary sinus sampling performed. Black and white bars show changes in coronary sinus and peripheral levels, respectively. G3P, glycerol-3-phosphate. *P < 0.05 vs. baseline in peripheral blood. #P < 0.05 between coronary sinus and peripheral samples.
Relative changes in plasma metabolite levels in coronary sinus versus periphery 10 minutes after PMI
Additional metabolites enriched in the coronary sinus 10 minutes after PMI
Relative changes in plasma metabolite levels in coronary sinus versus periphery 60 minutes after PMI
Additional metabolites enriched in the coronary sinus 60 minutes after PMI
Validation of metabolic markers in SMI. We examined whether our findings were applicable to a cohort of patients with SMI presenting for acute coronary angiography. Because the exact time of onset of SMIs relative to sample collection was heterogeneous (162 ± 102 minutes), we focused on the group of metabolites that changed significantly in the peripheral blood in a sustained pattern after PMI (Figure 1). These 4 metabolites in the platform — hypoxanthine, aconitic acid, TMNO, and threonine — were significantly changed at each of the 60-, 120-, and 240-minute time points in the entire PMI cohort (P < 0.05 for each). The median changes across these 3 time points compared with baseline are shown in Figure 3A. We next examined the difference in levels of each of these metabolites in the patients presenting with SMI compared with control patients presenting to the cardiac catheterization suite with nonacute cardiovascular disease (Figure 3A). There was concordance of both the direction and magnitude of differences in the SMI cohort with changes in the PMI cohort. We generated a simple composite score by summing the equally weighted intensities for each of these metabolites and confirmed that absolute MS intensity units in addition to relative changes successfully distinguished SMI patients from controls (P = 0.004; Figure 3B). This revealed excellent discriminatory ability, with receiver-operating-characteristic area under the curve (AUC) of 0.88 (Figure 3C). The composite score further confirmed that metabolic biomarkers derived in the PMI model were similarly altered in the SMI samples.
Metabolites whose levels significantly changed from baseline at 60, 120, and 240 minutes in PMI patients assessed in an independent SMI cohort. (A) Black bars denote average of median changes across the 3 time points compared with baseline values in PMI patients. White bars denote relative levels of each metabolite in patients presenting with SMI compared with control patients presenting to the cardiac catheterization suite with nonacute cardiovascular disease. (B) Composite scores for metabolites in A, derived by taking the weighted sum of metabolites that increased in PMI minus the sum of metabolites that decreased in PMI. Boxes denote interquartile range, lines denote median, and whiskers denote tenth and ninetieth percentiles. (C) Receiver-operating-characteristic curve for composite metabolite scores for each patient in the SMI cohort (cases) and the second control cohort (controls). AUC (mean ± SEM) is shown.
Validation of metabolic changes using isotope standards. Suppression or enhancement of the MS signal can occur from constituents present in a complex biological matrix such as plasma. We therefore performed additional studies using isotope-labeled metabolite standards for key metabolites in Tables 2–7 to provide unambiguous identification and quantification (see Methods). Supplemental Table 2 (supplemental material available online with this article; doi:
10.1172/JCI35111DS1) shows absolute and relative concentration changes for the metabolites with available isotope standards, relevant to Tables 2 and 3. The magnitude of changes in the original and follow-up studies were highly correlated (r2 = 0.9, P = 0.01). Similarly, the results in Supplemental Table 3 compare peripheral and coronary sinus samples using the isotopically labeled standards, relevant to Tables 6 and 7. These separate analyses also underscored our initial findings; the magnitude of changes in the original and follow-up studies were again highly correlated (r2 = 0.87, P = 0.01). Data for malic acid and glutamic acid, 2 metabolites that change early in the coronary sinus, are provided in Supplemental Table 4. Isotope standards were available for 3 of the 4 metabolites in the multimarker panel: TMNO, threonine, and hypoxanthine. Repeat analyses incorporating the isotope standards for these 3 markers in the PMI subjects versus cardiac catheterization controls (Supplemental Figure 3) were highly consistent with our initial findings (Figure 3A).
Functional consequences of metabolic changes. We used an in vitro cardiomyocyte model to investigate whether metabolites altered in PMI modulate the response to hypoxic challenge. Baseline apoptosis of cultured rat cardiomyocytes was 12.0% ± 1.9%, which increased to 24.2% ± 2.2% (P = 0.0001) after 3 hours of hypoxic challenge (see Methods) as assessed by FITC-conjugated annexin V staining, consistent with prior studies (20). Cell growth media was supplemented with candidate coronary sinus–enriched metabolites at the time of onset of hypoxia at concentrations approximately 2- to 10-fold those previously reported in the plasma (21–25). Figure 4A shows the percent change in apoptosis relative to the hypoxia control for the lowest concentration of each metabolite screened. Consistent with prior studies, pretreatment of cardiomyocytes with inosine diminished apoptosis (26), while hypoxanthine potentiated apoptosis following hypoxia (27). Xanthine, glycerol-3-phosphate, and malonic acid also diminished the hypoxia-triggered apoptotic response. Figure 4B shows a dose-dependent inhibition of apoptosis at physiologically relevant concentrations of inosine more than 100-fold lower than previously reported (26).
Cultured rat cardiomyocytes were subjected to a 3-hour hypoxic challenge as described in Methods. Cells were incubated with the indicated exogenous metabolites across a range of concentrations. Upon reoxygenation, apoptosis was assessed by FITC-labeled annexin V staining. (A) Percent change in apoptosis relative to the hypoxia control for the lowest concentration of each metabolite screened: xanthine, 14 μM; inosine, 10 μM; glycerol-3-phosphate, 1.7 mM; malonic acid, 6.24 μM; orotic acid, 19 mM; methylhistamine, 1.38 μM; succinic acid, 200 mM; taurine, 2 mM; hypoxanthine, 40 μM. (B) Percent change in apoptosis relative to the hypoxia control at physiologically relevant concentrations of inosine. Representative data (mean ± SEM) from 14 total experiments are shown. *P < 0.05, **P < 0.01 vs. hypoxia control.
Although emerging metabolomics profiling technologies hold enormous promise for illuminating biology, applications to human disease are still lacking. We therefore chose a clinical model of MI that allows for precise kinetic analysis in patients who served as their own controls. This experimental design markedly enhanced our power to identify statistically meaningful changes, while the targeted MS platform, in which metabolites are analyzed with high sensitivity and structural specificity, allowed us to unambiguously identify metabolites and thus avoid potentially confounding clinical variables (7). We hypothesized that small biochemicals may leak from injured myocardial cells before cellular damage would permit egress of macromolecules, potentially allowing for earlier detection of disease (28). Presently used indicators of myocardial injury (i.e., the myocardial isoform of CK-MB and cardiac troponins) reflect leakage of cardiac enzymes or structural proteins from irreversibly damaged cardiomyocytes, and are not reliably detected until at least 4 hours after myocardial injury (16). In contrast, the metabolic changes we identified in the PMI derivation cohort were readily apparent as early as 10 minutes after injury and were serially validated in distinct clinical cohorts. Coronary sinus catheterization provided the opportunity to sample directly from the organ of interest and served as another layer of validation that underscored the biological plausibility of the observed association between metabolic changes and MI. Finally, metabolite changes observed in humans modulated the response to hypoxic injury in vitro, which suggests that these metabolites not only serve as markers, but may also modulate the evolving response to ischemia in vivo.
The metabolic, structural, and functional consequences of ischemia and/or reperfusion have been examined in a wide variety of experimental animal models, including regional ischemia after coronary vessel occlusion in dogs (29–31) and swine (32–35) as well as in humans with coronary disease (36–40) or undergoing heart surgery (41–43). However, prior studies have generally assayed relatively limited subsets of metabolites in focused approaches. Some prior studies have also relied on invasive instrumentation, such as the implantation of dialysate catheters into the myocardium (32). In support of the applicability of a targeted liquid chromatography–MS (LC-MS) platform to human cardiovascular disease, we confirmed changes in metabolites that had been previously reported as appearing in either plasma or effluent from isolated, perfused hearts (32). These changes included the sequential purine degradation products (ADP, AMP, inosine, hypoxanthine, and xanthine), intermediates of anaerobic glycolysis (lactic acid) and the TCA cycle (succinic acid), and amino acids known to be concentrated within cardiomyocytes (taurine and glutamic acid) (32).
Analysis of patient plasma using a comprehensive LC-MS platform also identified additional metabolites in these pathways that had not to our knowledge been previously reported, such as the TCA cycle intermediates malic acid and aconitic acid. We have also identified what we believe to be novel changes in human plasma levels of pyrimidine metabolites, including orotic acid, malonic acid, and aminoisobutyric acid. Alterations in TMNO, in contrast, suggest injury-mediated modulation of dietary compounds in the plasma. The platform also identified the stable circulating derivative of the proinflammatory metabolite of histamine, methylhistamine (44).
While the application of metabolomics directly to humans has important advantages, functional studies must be pursued in complementary model systems. Our in vitro studies identified roles for several of the metabolites discussed above, including roles for glycerol-3-phosphate and malonic acid in diminishing hypoxia-triggered apoptosis in cardiomyocytes. Future studies will be necessary to compare metabolite profiles in intracellular versus extracellular compartments, to define tissue-specific responses to ischemia and injury, and ultimately to further define the functional roles for the observed metabolic alterations in animal models.
From a clinical perspective, metabolic responses arising uniquely from an abnormal hypertrophic myocardial substrate or the direct effects of ethanol must be considered; these could potentially restrict broader application to general populations. Ultimately, the validation of our findings in SMI patients mitigates the concern that ethanol-induced thrombosis and downstream injury generates a different metabolic profile from infarction as a result of coronary plaque rupture. However, further testing in larger cohorts will be necessary to confirm these preliminary findings, to define the specificity of the metabolic changes in the presence of comorbid conditions, and to extend their relevance to subgroups based on sex and race. Furthermore, while the findings in the PMI cohort suggest that our early markers are particularly useful in diagnosing very recent onset of injury, cohorts of “normal” MI patients presenting soon after symptom onset will be needed to further assess its clinical performance. Complementary studies of other body fluids, particularly urine, may also provide insight into metabolites that are rapidly cleared from the circulation.
It is unlikely that the metabolic changes we observed herein are attributable to medication exposures. Aspirin, beta blockers, and heparin were used almost uniformly in each of the cohorts studied (Table 1), removing them as potential confounders. When subjects were stratified by calcium channel blocker or statin use, no significant differences in metabolic responses were observed for any of the metabolites in Tables 2–5 (P > 0.05 for each). Large-scale studies will ultimately be necessary to identify any subtle potential metabolic alterations secondary to medications.
The recognized advantages of coupling LC with targeted MS/MS analysis include greatly increased sensitivity and highly specific identification of analytes. Targeted MS methods provide absolute measurement of analyte concentrations when appropriate standards are added, as we demonstrated in the present study. A targeted approach, however, is blind to changes in metabolites whose retention times and MS characteristics have not been incorporated into the analysis method. As efforts to define the human metabolome grow (25), we anticipate increasingly comprehensive targeted platforms for biomarker and pathway discovery.
We previously reported initial observations of the application of a more limited LC-MS/MS–based platform (10), in which many of the analytes had not been unambiguously identified, to patients undergoing cardiac exercise testing. Both the present work and the prior study detected circulating purine breakdown products, including hypoxanthine, which may not always be of cardiac muscle–specific origin. However, even with deletion of hypoxanthine from the composite score, the remaining 3 metabolites (aconitic acid, TMNO, and threonine) derived from the PMI subjects have excellent discrimination of SMI cases from controls, with AUC of 0.83. We hypothesize that in certain clinical contexts, such as myocardial ischemia at rest, elevations in hypoxanthine would be informative. In contrast, in the setting of vigorous exercise or diffuse muscle injury, the specificity of hypoxanthine for myocardial injury would be diminished. The present results also suggested the modulation of a distinct subset of TCA cycle intermediaries, amino acids, and amines not observed in the exercise cohort, which may reflect differences in the human models or distinct metabolic signatures of myocardial ischemia and MI. Ultimately, the integration of multiple data sets will be necessary for large-scale validation in patients to identify plasma biomarkers of myocardial ischemia and early injury, for which no circulating biomarkers currently exist (45).
In summary, we demonstrated the application of a metabolomics platform to a carefully phenotyped patient cohort for the discovery of blood markers with the potential to detect the presence of very early myocardial injury. We identified abnormalities in circulating metabolites as soon as 10 minutes after myocardial injury, a time frame in which, to our knowledge, no currently used biomarkers are elevated. Beyond the diagnostic utility we demonstrate here, other metabolic signatures may predict disease, establish a reference for return to normality, or aid in the design of new therapeutics for metabolic modulation.
Patients with HOCM undergoing septal ablation. A total of 36 patients undergoing PMI using alcohol septal ablation for the treatment of symptomatic HOCM were included in this study. Inclusion criteria for this cohort were (a) primary HOCM; (b) septal thickness of 16 mm or greater; (c) resting outflow tract gradient of greater than 30 mmHg, or an inducible outflow tract gradient of at least 50 mmHg; (d) symptoms refractory to optimal medical therapy; and (e) appropriate coronary anatomy. The most proximal accessible septal branch was instrumented using standard angioplasty guiding catheters and guidewires and 1.5 or 2.0 mm × 9 mm Maverick balloon catheters. Radiographic and echocardiographic contrast injections confirmed proper selection of the septal branch and balloon catheter position. Ethanol was infused through the balloon catheter at 1 ml/min. Additional injections in the same or other septal branches were administered as needed in order to cease blood flow to the isolated myocardium and reduce the gradient to less than 20 mmHg (46).
All protocols for obtaining blood from patients, including the coronary sinus sampling protocol, were approved by the Massachusetts General Hospital Institutional Review Board, and all subjects gave written informed consent. Blood was drawn at baseline, just prior to the onset of the ablation, and at 10 minutes and 1, 2, 4, and 24 hours after the onset of injury (Supplemental Figure 1). Of the 36 patients, 13 consented to the placement of a coronary sinus catheter during the ablation, allowing for the simultaneous sampling from coronary sinus and femoral catheters at the baseline, 10-minute, and 60-minute time points. The coronary sinus catheter was subsequently removed prior to the patient leaving the catheterization suite.
A control cohort of 16 patients undergoing elective, diagnostic cardiac catheterization for cardiovascular disease, but not acute myocardial ischemia, was also enrolled. Blood was drawn prior to the onset of cardiac catheterization and at 10 and 60 minutes after the procedure was begun.
Patients with SMI. We enrolled a cohort of 12 patients undergoing emergent cardiac catheterization for acute ST-segment elevation SMI within 8 hours of symptom onset. For this cohort, blood samples were obtained in the coronary catheterization suite. For the validation studies in SMI subjects, a second control cohort of 9 patients undergoing elective, diagnostic cardiac catheterization for cardiovascular disease, but not acute myocardial ischemia, was also enrolled.
Blood sample processing. Samples were obtained from femoral venous catheters during the procedure, or, where indicated, from a catheter placed in the coronary sinus. Samples were collected in K2EDTA-treated tubes (BD). All blood samples were centrifuged at 2,000 g for 10 minutes to pellet cellular elements. The supernatant plasma was then aliquoted and immediately frozen at –80°C to minimize freeze-thaw degradation. Additional blood samples were sent to the clinical chemistry laboratory for evaluation of the standard cardiac markers CK, CK-MB, and troponin T (Roche Diagnostics).
Metabolite selection for platform. We incorporated metabolites into our platform based on the following criteria. (a) With regard to biological diversity, we chose metabolites to provide sentinels in as many of the known metabolic pathways as possible (i.e., TCA cycle intermediates, purine biosynthesis, etc.). (b) With regard to biological relevance, we specifically added some metabolites with specific potential relevance to human metabolic and/or cardiovascular diseases. (c) With regard to biophysical characteristics, we eliminated metabolites that were not readily detectable by LC-MS methods or would require the addition of highly specific methods that could not be readily generalized to other classes of compounds. We also consulted available databases with reported concentrations to assess the feasibility of identifying specific metabolites in plasma.
HPLC and MS analysis. For the analysis of sugars, ribonucleotides, and organic acids, 200 μl plasma was subjected to ethanol precipitation (80% ethanol, 19.9% H2O, and 0.1% formic acid) at 4°C for 2 hours and centrifuged at 15, 000 g for 15 minutes, and the supernatant was evaporated in a nitrogen chamber at 30°C (Caliper Life Sciences). Samples were reconstituted in 60 μl H2O, and aliquots were separated sequentially and automatically by injection onto 3 HPLC columns with orthogonal separation characteristics as previously described (10). Sugars and ribonucleotides were separated on a Luna amino column (Phenomenex) under normal phase using acetonitrile, water, 0.25% ammonium hydroxide, and 10 mM ammonium acetate (pH 11) in a run time of 3.5 minutes. Organic acids were separated using a Synergi Polar-RP column (Phenomenex) under reverse phase using acetonitrile, water, and 5 mM ammonium acetate (pH 5.6–6.0) in a run time of 3.5 minutes. For the analysis of amino acids, plasma was diluted 10-fold with H2O and then subjected to reverse phase chromatography on a Luna phenyl-hexyl column (Phenomenex) using acetonitrile, water, and 0.1% acetic acid (pH 3.5–4.0) in a run time of 1.5 minutes.
The 3 HPLC columns were connected in parallel via an automated switching valve on a robotic sample loader (Leap Technologies) to a triple quadrupole mass spectrometer (AB4000Q; Applied Biosystem/Sciex), operated in electrospray ionization mode using a turbo ion spray LC-MS interface. Either positive or negative ions were selected for targeted MS/MS analysis using selective reaction monitoring conditions. A total of 210 known metabolites were monitored for each sample. Precursor and product ions of metabolites are shown in Supplemental Table 5. Metabolite quantification was performed by integrating peak areas for parent/daughter ion pairs using Multiquan software (version 1.0; Applied Biosystem/Sciex), and subsequently, all metabolite peaks were manually reviewed for peak quality in a blinded manner prior to statistical analysis.
LC-MS analyses with isotope standards. Representative metabolites from each of the tables were purchased from Cambridge Isotope Laboratories (Supplemental Table 1). LC-MS analyses ensured purity of greater than 98% (data not shown), and dose-response curves were established first in solvent, and subsequently in plasma. Endogenous levels of all metabolites were more than 2 orders of magnitude above the limit of quantitation and fell well within the linear range of the dose-response relationship (Supplemental Figure 2). The indicated amount of each metabolite (Supplemental Table 1) was added to plasma prior to sample preparation, and absolute quantitation was then determined by the ratio of the endogenous analyte peak to the isotope standard peak.
In vitro cardiomyocyte hypoxia model. Protocols for cardiomyocyte procurement and use in this study were approved by the IACUC of Beth Israel Deaconess Medical Center. Cardiomyocytes were isolated from 1-day-old rats and cultured for 72 hours, as previously described (20), in 100 mm dishes. A 3-hour hypoxic challenge was performed in a chamber (Billups-Rothenberg) filled with less than 0.5% O2, 95% nitrogen, and 5% CO2 at 37°C for 3 hours. Cells were incubated with exogenous metabolites across a range of concentrations. Upon reoxygenation, apoptosis was assessed by staining with FITC-labeled annexin V (Invitrogen). Counterstaining with Hoechst 33258 (4 mg/ml, 1:10 dilution in PBS of 40 mg/ml stock solution in DMSO; Invitrogen) allowed for normalization of annexin-labeled cells to total cell number.
Statistics. For clinical characteristics, values for continuous variables are presented as mean ± SD, and comparisons between groups were performed using 2-tailed, 2-sample Student’s t tests. Association between categorical variables was assessed using Fisher’s Exact test.
We performed preliminary studies using sample preparation and MS replicates of pooled human samples to assess the coefficient of variation (calculated as SD/mean of data set; expressed as a percentage) for the metabolites in the platform. This analysis showed that the aggregate coefficient of variation was approximately 20%. From the 36 patients in whom peripheral samples were collected in the PMI study, 20 patients were randomly selected for analysis as a derivation set. Levels of metabolites were tested for statistically significant percent change from baseline using the Wilcoxon signed-rank test. We used a significance threshold of P < 0.005 in the derivation cohort, because this threshold would be expected to yield approximately 1 false positive discovery of 210 metabolites analyzed, assuming independent hypotheses. False discovery rate for a given metabolite is the fraction of metabolites with the observed P value or lower that correspond to the null distribution. Conservative (i.e., upper bound) estimates of false discovery rate were obtained via the Benjamini-Hochberg procedure (47). The procedure was applied to the set of 210 measured metabolites for the comparisons of 10- and 60-minute time points with baseline.
Metabolites that changed significantly at the 10- or 60-minute time points in the PMI derivation cohort, and did not change significantly (P > 0.01) between the same time points in the control cohort of patients undergoing diagnostic catheterization without MI, were selected as candidate early biomarkers for testing in the PMI validation cohort that consisted of 16 patients. Criteria for validation was P < 0.05 by Wilcoxon signed-rank test with the direction of change concordant with that observed in the derivation cohort. The relationship between change in metabolites in the derivation and validation cohorts was assessed with a Spearman correlation coefficient. Data in all tables indicate median percent change from baseline and interquartile ranges.
Further analysis was carried out in the subgroup of 13 PMI patients with matched coronary sinus and peripheral samples. Metabolites were considered to be enriched in the coronary sinus if they changed significantly in the coronary sinus at 10 or 60 minutes compared with baseline (P < 0.05, Wilcoxon signed-rank test) and changed to a greater extent in the coronary sinus than in the periphery (median change in coronary sinus at least 1.3-fold that in the periphery; P < 0.05, Wilcoxon signed-rank test).
To evaluate whether metabolic changes observed in the PMI patients were generalizable to SMI, we selected all of those metabolites that displayed significant changes from baseline at 1, 2, and 4 hours in the derivation and validation PMI cohorts (P < 0.05 at all 3 time points). A Wilcoxon rank-sum test was used to examine levels of these individual metabolites in the patients presenting with SMI compared with control patients presenting to the cardiac catheterization suite with nonacute cardiovascular disease. These metabolites were also compiled into a composite MS intensity unit score for SMI and control patients. To ensure equal weighting of each metabolite in this composite score, the intensity values of each metabolite were rescaled to have a common median intensity of 1.0 × 106 arbitrary units. The composite score was defined as the sum of metabolites that increased in PMI minus the sum of metabolites that decreased in PMI.
The authors gratefully acknowledge support from the NIH (grants R01 HL072872 to R.E. Gerszten and M.S. Sabatine; U01HL083141 to R.E. Gerszten, M.S. Sabatine, and F.P. Roth; and R01DK081572 to R.E. Gerszten), the Donald W. Reynolds Foundation (to R.E. Gerszten and M.S. Sabatine), Fondation Leducq (to A. Rosenzweig and R.E. Gerszten), Heart Failure Society of America (to G.D. Lewis), an American Heart Association Fellow-to-Faculty Award (to G.D. Lewis), the Harvard/MIT Clinical Investigator Training Program (to G.D. Lewis), and the Broad Institute Scientific Planning and Allocation of Resources Committee to establish the metabolomics platform. A. Asnani was supported by a predoctoral award from the Sarnoff Cardiovascular Research Foundation. F.P. Roth was also supported in part by National Human Genome Research Institute, NIH, grants HG003224, HG0017115, NS054052, and HG004233 and by the Keck Foundation.
Nicholson, J.K., Wilson, I.D. 2003. Opinion: understanding ‘global’ systems biology: metabonomics and the continuum of metabolism. Nat. Rev. Drug Discov. 2:668-676.
Raamsdonk, L.M., et al. 2001. A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations. Nat. Biotechnol. 19:45-50.
Allen, J., et al. 2003. High-throughput classification of yeast mutants for functional genomics using metabolic footprinting. Nat. Biotechnol. 21:692-696.
An, J., et al. 2004. Hepatic expression of malonyl-CoA decarboxylase reverses muscle, liver and whole-animal insulin resistance. Nat. Med. 10:268-274.
Beecher, C.W.W. 2003. Chap. 17 inThe human metabolome.
Kluwer Academic Publishers. Boston, Massachusetts, USA. 311–319.
He, W., et al. 2004. Citric acid cycle intermediates as ligands for orphan G-protein-coupled receptors. Nature. 429:188-193.
Brindle, J.T., et al. 2002. Rapid and noninvasive diagnosis of the presence and severity of coronary heart disease using 1H-NMR-based metabonomics. Nat. Med. 8:1439-1444.
Kirschenlohr, H.L., et al. 2006. Proton NMR analysis of plasma is a weak predictor of coronary artery disease. Nat. Med. 12:705-710.
Lee, M.S., Kerns, E.H. 1999. LC/MS applications in drug development. Mass. Spectrom. Rev. 18:187-279.
Sabatine, M.S., et al. 2005. Metabolomic identification of novel biomarkers of myocardial ischemia. Circulation. 112:3868-3875.
Sigwart, U. 1995. Non-surgical myocardial reduction for hypertrophic obstructive cardiomyopathy. Lancet. 346:211-214.
Knight, C., et al. 1997. Nonsurgical septal reduction for hypertrophic obstructive cardiomyopathy: outcome in the first series of patients. Circulation. 95:2075-2081.
Lakkis, N.M., et al. 1998. Echocardiography-guided ethanol septal reduction for hypertrophic obstructive cardiomyopathy. Circulation. 98:1750-1755.
Lakkis, N.M., Nagueh, S.F., Dunn, J.K., Killip, D., Spencer (3rd.), W.H. 2000. Nonsurgical septal reduction therapy for hypertrophic obstructive cardiomyopathy: one-year follow-up. J. Am. Coll. Cardiol. 36:852-855.
Yoerger, D.M., et al. 2006. Time course of pressure gradient response after first alcohol septal ablation for obstructive hypertrophic cardiomyopathy. Am. J. Cardiol. 97:1511-1514.
Zimmerman, J., et al. 1999. Diagnostic marker cooperative study for the diagnosis of myocardial infarction. Circulation. 99:1671-1677.
Ransohoff, D.F. 2005. Bias as a threat to the validity of cancer molecular-marker research. Nat. Rev. Cancer. 5:142-149.
Ransohoff, D.F. 2004. Rules of evidence for cancer molecular-marker discovery and validation. Nat. Rev. Cancer. 4:309-314.
Roberts, R., Gowda, K.S., Ludbrook, P.A., Sobel, B.E. 1975. Specificity of elevated serum MB creatine phosphokinase activity in the diagnosis of acute myocardial infarction. Am. J. Cardiol. 36:433-437.
Matsui, T., et al. 1999. Adenoviral gene transfer of activated PI 3-kinase and Akt inhibits apoptosis of hypoxic cardiomyocytes in vitro. Circulation. 100:2373-2379.
Albert, J.D., et al. 1986. Extremity amino acid metabolism during starvation and intravenous refeeding in humans. Am. J. Physiol. 251:E604-E610.
Scriver, C.R., Gregory, D.M., Sovetts, D., Tissenbaum, G. 1985. Normal plasma free amino acid values in adults: the influence of some common physiological variables. Metabolism. 34:868-873.
Traut, T.W. 1994. Physiological concentrations of purines and pyrimidines. Mol. Cell. Biochem. 140:1-22.
Wevers, R.A., Engelke, U., Heerschap, A. 1994. High-resolution 1H-NMR spectroscopy of blood plasma for metabolic studies. Clin. Chem. 40:1245-1250.
Wishart, D.S., et al. 2007. HMDB: The Human Metabolome Database. Nucleic Acids Res. 35:D521-D526.
Goldhaber, S.Z., et al. 1982. Inosine: a protective agent in an organ culture model of myocardial ischemia. Circ. Res. 51:181-188.
Xia, Y., Zweier, J.L. 1995. Substrate control of free radical generation from xanthine oxidase in the postischemic heart. J. Biol. Chem. 270:18797-18803.
Song, D., O’Regan, M.H., Phillis, J.W. 1998. Mechanisms of amino acid release from the isolated anoxic/reperfused rat heart. Eur. J. Pharmacol. 351:313-322.
Dorheim, T.A., Wang, T., Mentzer (Jr.), R.M., Van Wylen, D.G. 1990. Interstitial purine metabolites during regional myocardial ischemia. J. Surg. Res. 48:491-497.
Delyani, J.A., Van Wylen, D.G. 1994. Endocardial and epicardial interstitial purines and lactate during graded ischemia. Am. J. Physiol. 266:H1019-H1026.
Mei, D.A., Gross, G.J., Nithipatikom, K. 1996. Simultaneous determination of adenosine, inosine, hypoxanthine, xanthine, and uric acid in microdialysis samples using microbore column high-performance liquid chromatography with a diode array detector. Anal. Biochem. 238:34-39.
Backstrom, T., et al. 2003. Cardiac outflow of amino acids and purines during myocardial ischemia and reperfusion. J. Appl. Physiol. 94:1122-1128.
Zemgulis, V., et al. 2001. Energy-related metabolites during and after induced myocardial infarction with special emphasis on the reperfusion injury after extracorporeal circulation. Acta. Physiol. Scand. 171:129-143.
Hassel, B., Ilebekk, A., Tonnessen, T. 1998. Cardiac accumulation of citrate during brief myocardial ischaemia and reperfusion in the pig in vivo. Acta. Physiol. Scand. 164:53-59.
Osterlund, B., et al. 2002. Myocardial ischemia induces coronary t-PA release in the pig. Acta. Anaesthesiol. Scand. 46:271-278.
Mudge (Jr.), G.H., Mills (Jr.), R.M., Taegtmeyer, H., Gorlin, R., Lesch, M. 1976. Alterations of myocardial amino acid metabolism in chronic ischemic heart disease. J. Clin. Invest. 58:1185-1192.
Srinivasan, K.N., et al. 1998. Comparison of glycoprotein components, tryptophan, lipid peroxidation and antioxidants in borderline and severe hypertension and myocardial infarction. Clin. Chim. Acta. 275:197-203.
Wirleitner, B., et al. 2003. Immune activation and degradation of tryptophan in coronary heart disease. Eur. J. Clin. Invest. 33:550-554.
Turgan, N., et al. 1999. Urinary hypoxanthine and xanthine levels in acute coronary syndromes. Int. J. Clin. Lab. Res. 29:162-165.
Pisarenko, O.I., et al. 1989. Features of myocardial metabolism of some amino acids and ammonia in patients with coronary artery disease. Eur. Heart J. 10:209-217.
Smolenski, R.T., et al. 1993. Formation and breakdown of uridine in ischemic hearts of rats and humans. J. Mol. Cell. Cardiol. 25:67-74.
Kennergren, C., et al. 1999. Extracellular amino acids as markers of myocardial ischemia during cardioplegic heart arrest. Cardiology. 91:31-40.
Svedjeholm, R., et al. 1991. Myocardial uptake of amino acids and other substrates in relation to myocardial oxygen consumption four hours after cardiac operations. J. Thorac. Cardiovasc. Surg. 101:688-694.
Keyzer, J.J., Breukelman, H., Wolthers, B.G., Richardson, F.J., de Monchy, J.G. 1985. Measurement of N tau-methylhistamine concentrations in plasma and urine as a parameter for histamine release during anaphylactoid reactions. Agents Actions. 16:76-79.
Morrow, D.A., de Lemos, J.A., Sabatine, M.S., Antman, E.M. 2003. The search for a biomarker of cardiac ischemia. Clin. Chem. 49:537-539.
Baggish, A.L., et al. 2006. Pathological effects of alcohol septal ablation for hypertrophic obstructive cardiomyopathy. Heart. 92:1773-1778.
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. 57:289-300.