Author + information
- Received August 12, 2015
- Revision received October 6, 2015
- Accepted October 13, 2015
- Published online January 19, 2016.
- Gregory D. Lewis, MD∗,†,‡∗ (, )
- Debby Ngo, MD†,
- Anna R. Hemnes, MD§,
- Laurie Farrell, RN∗,
- Carly Domos, MS∗,
- Paul P. Pappagianopoulos, MEd†,
- Bishnu P. Dhakal, MD∗,
- Amanda Souza, MS‡,
- Xu Shi, PhD∗,
- Meredith E. Pugh, MD, MSCI§,
- Arkadi Beloiartsev, MD∗,
- Sumita Sinha, PhD∗,
- Clary B. Clish, PhD‡ and
- Robert E. Gerszten, MD∗,‡,‖∗∗ ( )()
- ∗Cardiology Division, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts
- †Pulmonary and Critical Care Unit, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts
- ‡Broad Institute of MIT and Harvard, Cambridge, Massachusetts
- §Vanderbilt University Pulmonary Unit, Nashville, Tennessee
- ‖Cardiovascular Research Center, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts
- ↵∗Reprint requests and correspondence:
Dr. Gregory D. Lewis, Cardiology Division, Massachusetts General Hospital, Bigelow 800, Fruit Street, Boston, Massachusetts 02114.
- ↵∗∗Dr. Robert E. Gerszten, Cardiovascular Research Center, Massachusetts General Hospital, CPZ 185, 185 Cambridge Street, Boston, Massachusetts 02114.
Background Pulmonary hypertension and associated right ventricular (RV) dysfunction are important determinants of morbidity and mortality, which are optimally characterized by invasive hemodynamic measurements.
Objectives This study sought to determine whether metabolite profiling could identify plasma signatures of right ventricular-pulmonary vascular (RV-PV) dysfunction.
Methods We measured plasma concentrations of 105 metabolites using targeted mass spectrometry in 71 individuals (discovery cohort) who underwent comprehensive physiological assessment with right-sided heart catheterization and radionuclide ventriculography at rest and during exercise. Our findings were validated in a second cohort undergoing invasive hemodynamic evaluations (n = 71), as well as in an independent cohort with or without known pulmonary arterial (PA) hypertension (n = 30).
Results In the discovery cohort, 21 metabolites were associated with 2 or more hemodynamic indicators of RV-PV function (i.e., resting right atrial pressure, mean PA pressure, pulmonary vascular resistance [PVR], and PVR and PA pressure-flow response [ΔPQ] during exercise). We identified novel associations of RV-PV dysfunction with circulating indoleamine 2,3-dioxygenase (IDO)–dependent tryptophan metabolites (TMs), tricarboxylic acid intermediates, and purine metabolites and confirmed previously described associations with arginine–nitric oxide metabolic pathway constituents. IDO-TM levels were inversely related to RV ejection fraction and were particularly well correlated with exercise PVR and ΔPQ. Multisite sampling demonstrated transpulmonary release of IDO-TMs. IDO-TMs also identified RV-PV dysfunction in a validation cohort with known risk factors for pulmonary hypertension and in patients with established PA hypertension.
Conclusions Metabolic profiling identified reproducible signatures of RV-PV dysfunction, highlighting both new biomarkers and pathways for further functional characterization.
There is an unmet clinical need to develop new circulating biomarkers of pulmonary hypertension (PH). Novel biomarkers of right ventricular (RV)–pulmonary vascular (PV) dysfunction might aid in earlier PH detection and inform selection of patients for PH-specific therapies, while reducing reliance on subjective measures (e.g., functional class) and invasive hemodynamic assessment. The profiling of metabolites, including lipids, sugars, nucleotides, amino acids, and amines, is particularly relevant to the understanding of RV-PV dysfunction. Studies in both experimental models and patients have suggested a maladaptive imbalance between pulmonary vasodilators and vasoconstrictors in RV-PV dysfunction. Several small-molecule mediators of PV tone have been identified, including the nitric oxide (NO) second messenger cyclic guanosine monophosphate, prostaglandins, and other biogenic amines (1–3). PV remodeling and right ventricular (RV) hypertrophy have been associated with a potential switch from fatty acid oxidation to glycolysis (4,5), which provides further motivation to apply a more global metabolite profiling approach to PH. Ultimately, metabolite profiling might highlight measures to redress metabolic derangements in RV-PV dysfunction (6,7).
In this study, we hypothesized that simultaneous invasive hemodynamic measurements during right-sided heart catheterization and blood sampling in dyspneic subjects would permit derivation of metabolic signatures of RV-PV dysfunction. We recruited subjects with dyspnea and a wide range of right-sided heart catheterization hemodynamic profiles, as well as predefined groups of people with and without known World Health Organization Group 1 pulmonary arterial hypertension (PAH), to identify metabolite profiles of RV-PV dysfunction. We performed blood sampling in multiple anatomic sites to localize the source of metabolic signatures. To optimize the precision of our phenotyping of RV-PV dysfunction, we also measured the degree of impairment in exercise capacity (peak oxygen uptake [Vo2]) and performed minute-by-minute hemodynamic measurements across varying RV loading conditions during exercise to probe the capacity of the RV-PV unit to accommodate increased blood flow (8). Furthermore, we used a well-described animal model to evaluate the expression of key enzymes and metabolites in an emerging pathway associated with PH in our human studies.
We enrolled 2 consecutive cohorts of patients who underwent cardiopulmonary exercise testing (CPET) with invasive hemodynamic monitoring at Massachusetts General Hospital to evaluate dyspnea on exertion. Patients provided informed consent to participate in this study. Exclusion criteria consisted of the following: 1) left ventricular (LV) systolic dysfunction (left ventricular ejection fraction [LVEF] <0.45) or severe valvular heart disease; 2) chronic kidney disease (worse than stage 3); 3) parenchymal lung disease with a pulmonary mechanical limitation to exercise as defined by minute ventilation (VE)/(forced expiratory volume in 1 s [FEV1] × 35) >0.7 at the ventilatory anaerobic threshold (9); or 4) incomplete pulmonary arterial catheter pressure measurements.
A third cohort was prospectively recruited from patients referred to the Vanderbilt Pulmonary Hypertension Clinic who were diagnosed with right-sided heart catheterization–proven PAH attributable to idiopathic, heritable, drug- or toxin-associated, or connective tissue disease–associated PAH and who provided consent for having their blood drawn. The control subjects were recruited prospectively from an e-mail distribution to all Vanderbilt employees seeking healthy individuals of similar age, sex, and body mass index (BMI) to enroll as control subjects in their PAH clinic blood sampling study. The human studies were approved by the Partners Healthcare System and Vanderbilt University Institutional Review Boards.
Subjects in the first 2 cohorts were mandated to be fasting for at least 8 h before CPET and blood draws for metabolite measurements. Invasive CPET with first-pass radionuclide ventriculography was performed as described previously (7,10) (Online Appendix, Online Figure 1). In the third PAH versus control validation cohort, peripheral venous plasma samples were obtained in the fasting state at 3 (median interquartile range: 0 to 6) days from the time of right-sided heart catheterization documenting PAH.
Profiles of plasma metabolites were obtained by use of a multiplexed liquid chromatography mass spectrometry (LC-MS) system as described previously (11) (Online Appendix). Additionally, commercially available enzyme-linked immunosorbent assays were used to measure N-terminal pro–B-type natriuretic peptide (NT-proBNP), cystatin C, and interleukin (IL)-6 levels from separate aliquots of the same samples used for LC-MS analysis of metabolites.
Animal model of PH
Male C57Bl/6 mice (8 to 10 weeks old) were obtained from the Jackson Laboratory (Bar Harbor, Maine). Mice were housed in specially designed cage racks and exposed to chronic hypoxia (Fio2 0.10) or room air (control) for 3 weeks. All animal studies were approved by the Institutional Animal Care and Use Committee of the Massachusetts General Hospital and conformed to the position of the American Heart Association on animal use. After mice were euthanized, blood was collected via LV puncture, and a thoracotomy was performed for heart and lung harvesting. RV hypertrophy was calculated as the ratio of the RV wall weight to LV wall plus septum weight (RV/[LV + S] − Fulton ratio).
Blood samples from the mice were centrifuged at 8,000 revolutions/min for 5 min at 4°C. Plasma was immediately stored at −80°C. Lung tissue was snap frozen in liquid nitrogen immediately after harvest and stored at −80°C. For homogenization of the lung, 25 mg of tissue sample was mixed with 250 μl of 50:50 methanol:water solution. The tissue was homogenized for 4 min at 25 Hz in the Tissue Lyser II (Qiagen, Hilden, Germany). Metabolite profiling was performed on plasma and lung homogenate as described previously for the human plasma samples. Indoleamine 2,3-dioxygenase-dependent (IDO) ribonucleic acid (RNA) content was measured by quantitative real-time polymerase chain reaction (Online Appendix).
All metabolite values and NT-proBNP levels were log-transformed because of non-normal distribution. Age- and sex-adjusted Pearson correlation coefficients were calculated to determine correlations between metabolites known to cluster within well-defined groups (e.g., arginine-NO metabolites and tryptophan metabolites). Linear regression analyses were performed to examine the relation of each metabolite with hemodynamic measurements that reflect RV-PV function: right atrial pressure (RAP), pulmonary arterial pressure (PAP), and pulmonary vascular resistance (PVR) at rest, and PVR and change in mean PAP relative to change in cardiac output (ΔPQ) during exercise. Relationships between metabolites and pulmonary arterial wedge pressure (PAWP) are included for the sake of comparison of metabolite relationships to left-sided heart filling pressures.
Multimarker scores were created for metabolites that cluster into known metabolic pathways if multiple constituents demonstrated significant associations with hemodynamic indicators of RV-PV function. Scores were derived by calculating the sum of standardized biomarker values (z-scores) weighted according to the value of their corresponding β-coefficients, as described previously (11). Multimarker scores fit to the derivation cohort that were significantly associated with RV-PV dysfunction were then evaluated by linear regression in the second validation cohort. Receiver operating coefficient analyses were performed to identify the ability of metabolites to predict RV-PV dysfunction in the validation cohort. The C-statistic was calculated by assessing the proportion of cases (RV-PV dysfunction) for which the biomarker value exceeded that in controls (no RV-PV dysfunction).
To further validate novel circulating metabolites in an independent cohort, we compared levels of IDO tryptophan metabolites (IDO-TMs) in patients with known PAH and control subjects. To have 80% power to detect a 20% difference in 4 IDO-TMs with p < 0.0125 (Bonferroni correction p = 0.05 for 4 metabolites) between groups, we performed metabolite profiling in 11 patients with PAH and 19 control subjects. Between-group comparisons were performed for the PAH versus control cohort with Student t test or Wilcoxon rank sum test, as appropriate. Additional information regarding adjustment for multiple hypothesis testing and regression models is available in the Online Appendix. The STATA version 12.0 software package (StataCorp LP, College Station, Texas) was used for statistical analysis.
Study participants undergoing CPET had an average age of 62 years in both cohorts, with a slight female predominance (Table 1). Average resting right-sided heart catheterization measurements were in the high-normal range. No patients demonstrated LV systolic dysfunction, as defined by either an LVEF <0.45 at rest or a fall in LVEF to <0.45 during exercise. Exercise unmasked impaired RV-PV reserve function in both the derivation and validation groups, as evidenced by multipoint changes in ΔPQ and peak exercise PVR (Table 1). In the setting of these hemodynamic profiles, average exercise capacity was reduced in both the derivation and validation cohorts (peak Vo2 = 62 ± 13% predicted and 66 ± 14% predicted (12), respectively).
Analytes associated with RV-PV dysfunction
We measured 5 hemodynamic indicators of RV-PV dysfunction (resting RAP, PAP, PVR, exercise PVR, and ΔPQ). We integrated the hemodynamic measurements with mass spectrometry–based analyses of metabolites (Figure 1), identifying 21 metabolites that were significantly associated with 2 or more measurements of RV-PV dysfunction in a regression analysis (p < 0.0095 for each analyte).
Many metabolites clustered within previously defined pathways. Levels of arginine-NO metabolites (Arg-Ms: arginine, ornithine, citrulline, asymmetric dimethylarginine, and symmetric dimethylarginine) were related to indexes of RV-PV dysfunction (Figure 1, Table 2). The ratio of arginine to ornithine + citrulline, which reflects global arginine bioavailability and has emerged as a potential PH biomarker in select PH populations with sickle cell disease (13) and heart failure (1), was inversely related to PAP, PVR, and ΔPQ (all p < 0.005) (Table 2). A cluster of purine degradation products (i.e., purine-Ms: allantoin, xanthosine, inosine, xanthine, uric acid) was also closely related to RV-PV dysfunction. Several of these purine-Ms are associated with oxidative stress, although only urate has been associated with PH (14).
The amino acid tryptophan can be metabolized via 2 enzymatic pathways: tryptophan hydroxylase (TH), which yields serotonin (5-hydroxytryptophan) and the main metabolite of serotonin (5-hydroxyindoleacetic acid), or via IDO, which generates kynurenine, anthranilate, kynurenate, and quinolinate. We identified striking associations between RV-PV dysfunction and IDO-TMs (Figure 1, Table 2) that have not previously been associated with PH. By contrast, tryptophan metabolites downstream from TH were not related to hemodynamic indexes of RV-PV dysfunction (Figure 1). Additionally, plasma levels of tricarboxylic acid (TCA) cycle metabolites (TCA-Ms) were closely associated with hemodynamic indexes of RV-PV dysfunction. Although TCA intermediates generally function intracellularly, we have previously documented that these metabolites are detectable in plasma (15).
We next assessed interpathway and intrapathway correlations. Relationships across annotated pathways are displayed in Online Table 1. Moderate correlations within groups of related metabolites were observed for IDO-TMs, TCA-Ms, purine-Ms, and Arg-Ms, whereas TH-TMs were not strongly interrelated. Although many of the metabolites highlighted by the LC-MS based platform fell into previously described pathways, other metabolites had less obvious metabolic connections (Figure 1). We identified associations between indexes of RV-PV dysfunction and the catecholamine vanillylmandelic acid, as well as aminoisobutyric acid and cyclic adenosine monophosphate, 2 additional metabolites potentially related to sympathetic nervous system activation (16,17). Plasma levels of the homocysteine metabolite cystathionine, as well as glyceraldehyde and glucoronate, were also related to indexes of RV-PV dysfunction.
To assess how clusters of metabolites from a given metabolic pathway relate to hemodynamic measurements, composite scores were generated by summing z-scores of individual metabolites weighted by β-coefficients within clusters (18) (Table 2). Composite scores from each group, adjusted for age and sex, were closely related to PAP and PVR at rest (Table 2). By contrast, none of the metabolite scores were significantly associated with systemic MAP (data not shown). Although purine-Ms and Arg-Ms demonstrated similar β-coefficients for PAP and PAWP, IDO-TMs and TCA-Ms demonstrated β-coefficients <0.3 with PAWP and ≥0.4 for PAP, which suggests relative specificity for pre-capillary PAP as opposed to PAP driven by left-sided filling pressure. Moreover, IDO scores remained associated with PAP, PVR, and ΔPQ (p < 0.05) after adjustment for each of the other pathway scores. Individuals in the top quartile of the IDO-TM score had a 2.9-fold greater ΔPQ (Figure 2) and a 2.2-fold greater PVR at rest and with exercise (p for trend < 0.001 for all). IDO-TMs remained significantly related to indexes of RV-PV dysfunction after further adjustment for NT-proBNP, a marker of myocardial wall stress and PH (19); cystatin C, indicative of renal function; and IL-6, an inflammatory cytokine previously associated with PH (Table 3). After stepwise adjustment for the combination of age, sex, NT-proBNP, cystatin C, and IL-6, IDO-TM scores remained significantly associated with RAP, mPAP, PVR, exercise PVR, and ΔPQ (p ≤ 0.01 for all) (Table 3). Of the 4 IDO-TMs, only kynurenate was significantly associated with serum creatinine (β = 0.41; p = 0.003). Furthermore, IDO-TM scores were not significantly associated with BMI, diabetes, hypertension, smoking status, or exposure to beta-blockers, angiotensin-converting enzyme inhibitors/angiotensin receptor blockers, calcium channel blockers, diuretic agents, statins, aspirin, or warfarin (p > 0.1 for all) (Online Table 2).
The novel finding in humans that IDO-dependent metabolites are closely related to indexes of RV-PV dysfunction led us to further examine their relationship with other indicators of RV-PV dysfunction beyond hemodynamic variables. IDO-TM composite scores were inversely associated with right ventricular ejection fraction (RVEF) (β = 0.38, p = 0.007), and the association with RVEF remained significant after adjustment for age, sex, and NT-proBNP levels (p = 0.015). Furthermore, because exercise intolerance is a cardinal manifestation of RV-PV dysfunction, and reduced peak Vo2 is associated with poor prognosis in various forms of PH (20–22), we examined whether IDO-TMs also predicted impaired exercise capacity. We found that higher levels of IDO-TMs were associated with reduced percent-predicted peak Vo2, even after adjusting for NT-proBNP (β = 0.36; p = 0.009). IDO-TMs thus provide markers of RV-PV dysfunction that signal both hemodynamic impairment and reduced functional capacity.
Localization of metabolic changes in human studies
To examine whether there was net uptake or release of IDO-TMs across the pulmonary circulation, we examined relative concentrations of the 4 IDO-TMs in radial arterial and proximal pulmonary arterial samples. In subjects with normal resting PVR, there were no significant differences in pulmonary and radial arterial levels of IDO-TMs (Figure 3); however, subjects with elevated PVR (>2 Wood units) demonstrated significantly higher radial arterial levels of kynurenine, anthranilate, and quinolinate than pulmonary arterial levels (p < 0.05 for all), which indicates release of these metabolites into the pulmonary circulation. The observed gradient in IDO-TMs across the pulmonary circulation in individuals with elevated PVR suggests a pulmonary origin of the observed increases in circulating IDO-TMs in patients with PV dysfunction.
Localization of metabolic changes in an animal model
To add information regarding tissue specificity, IDO and IDO-TMs were examined in lung tissue from an experimental animal model of hypoxia-induced PH (Figure 4). RV hypertrophy was confirmed by an elevated Fulton ratio of 0.35 ± 0.02 in the hypoxia group versus 0.22 ± 0.02 in controls (p < 0.001). Lung messenger RNA expression of IDO was 2.5-fold higher in mice with PH than in controls (p < 0.001) (Figure 4A); the IDO substrate-product ratio was elevated in lung tissue (p < 0.001) (Figure 4B); plasma levels of IDO-TMs were elevated in experimental PH (mean z-score: 0.42 vs. −0.50, p = 0.004) (Figure 4C), and their degree of elevation was closely correlated to Fulton ratios (r = 0.48 to 0.72; all p < 0.05).
Additional validation studies
We assessed the performance of our metabolite biomarkers in a second set of CPET subjects who served as a validation cohort. Multimarker combinations of metabolites associated with RV-PV dysfunction, identical to those identified in the derivation cohort (Table 2), demonstrated significant associations with indexes of RV-PV dysfunction in the validation cohort, with the majority of β-coefficients >0.25 and p < 0.05 (Table 4). We subsequently stratified the validation group according to whether or not abnormal pulmonary vascular function was present (+PH; n = 35), as evidenced by mPAP ≥21 mm Hg with PVR >2.0 Wood Units and ΔPQ >3.0 mm Hg/l/min (8–11,23,24), or absent (−PH; n = 36). Receiver operating characteristic curves generated for each IDO-TM in the validation cohort demonstrated C-statistics that ranged from 0.66 to 0.74 (Figure 5A). IDO-TMs retained their significance in logistic regression models that included variables known to be associated with PH (age, sex, BMI, RVEF, and PAWP), with C-statistic increments ranging from 0.02 to 0.09 (Table 5).
To further assess the performance of IDO-TMs as potential markers of RV-PV dysfunction, we also examined IDO-TM levels in an independent cohort of patients with known World Health Organization Group 1 PAH and control subjects of similar age, sex, and BMI (Table 6). The 3 IDO-TMs with transpulmonary release in the high-PVR group (anthranilate, quinolinate, and kynurenine) demonstrated significantly higher levels in PAH, with C-statistics of 0.79 to 0.86 (Figure 5B); kynurenate did not differ between groups.
Metabolite profiling permits comprehensive, relatively unbiased investigation of metabolic pathways that potentially underlie complex phenotypes such as RV-PV dysfunction. Using a mass spectrometry–based metabolite profiling platform, we identified a panel of metabolites for which the plasma concentrations were closely associated with RV-PV dysfunction. Our platform highlighted small molecules previously associated with RV-PV dysfunction, including arginine-NO pathway intermediates, catecholamines, and indicators of oxidative stress. Furthermore, we identified circulating metabolites (TCA cycle intermediates and IDO-TMs) not previously associated with RV-PV dysfunction in humans (Central Illustration), which represent novel potential biomarkers of PH. IDO-TMs are particularly attractive potential markers of RV-PV dysfunction because of their transpulmonary release in patients with high PVR, their relationship with RVEF and steep PQ, and their elevation in human and experimental PAH. Taken together, our data indicated that metabolite profiling may provide biomarkers of PV disease and a window into the metabolic pathways altered in RV-PV dysfunction.
A growing number of studies have evaluated potential circulating biomarkers of PV disease (25,26); however, studies to date have been largely limited to 1 or a few analytes that are compared to single resting noninvasive estimates or measurements of right-sided heart hemodynamics (1,14,27). A strength of our investigation was the application of a well-established targeted LC-MS metabolite profiling platform to measure more than 100 metabolites in carefully phenotyped derivation and validation cohorts. Our analysis benefited from blood sampling performed in parallel with “gold standard” invasive hemodynamic measurements at rest. We also characterized multipoint pressure-flow (PQ) relationships during exercise to further refine patient phenotypes related to RV-PV function (Online Figure 1). Our group and others have recently reported that ΔPQ responses to exercise predict functional capacity and outcomes (8,23,28,29). By linking metabolic signatures to PQ responses indicative of PV reserve capacity, we may gain insight into early markers of RV-PV dysfunction. Ultimately, biomarkers of RV-PV dysfunction will provide the greatest utility if they can differentiate early forms of PH before overt disease ensues.
We found that the ratio of arginine to ornithine + citrulline (Arg/[Orn+Cit]), which reflects global arginine bioavailability, was strongly inversely related to PAP, PVR during rest and exercise, and PQ (Figure 1, Table 2). These findings extend previous work demonstrating that Arg/(Orn+Cit) ratio is a marker of PH severity in select PH-populations with sickle cell disease (13) and LV systolic dysfunction (1). Our study documented a novel association of these metabolites with pulmonary hemodynamic measurements in an unselected population with exertional dyspnea. The particularly strong relationship between Arg/(Orn+Cit) levels and ΔPQ suggested that this ratio may signal impaired NO-mediated vasodilation.
In contrast to arginine pathway metabolites, the association of IDO-TMs with PH in humans has not been established. The consistent, strong correlations between the 4 IDO-TMs and indexes of RV-PV dysfunction, coupled with their transpulmonary release and association with RVEF, suggest that these markers have specificity for the pulmonary circulation. IDO has a known immunomodulatory role upon induction by proinflammatory stimuli, including interferon-γ, and is expressed in fibroblasts, macrophages, and vascular endothelium (30). The IDO-product kynurenine was identified as a potential endothelium-derived vasodilator in a murine sepsis model (31). In vitro studies by the same group demonstrated that kynurenine triggered dilation of previously constricted porcine arteries. Recently, IDO expression was localized to pulmonary endothelium in mouse lung, and mice lacking the IDO gene demonstrated exaggerated hypoxic pulmonary vasoconstriction (32). Enhanced endothelial IDO ameliorated PH in preclinical models, modulating PH-derived pulmonary artery smooth muscle cells toward a proapoptotic and less proliferative/synthetic state. However, to the best of our knowledge, there have been no studies linking this pathway to human disease.
IDO-TMs may be released as vasodilators in response to PH just as natriuretic peptides are secreted in the setting of myocardial wall stretch. As a heme-containing enzyme, IDO may be inhibited by NO (33), and it is plausible that production of IDO-TMs represents a compensatory pathway activated in states of NO deficiency, such as RV-PV dysfunction. Our finding that IDO-TMs were inversely correlated with plasma markers of Arg-NO bioavailability (Online Table 1) supported this hypothesis, although further investigation of the relationship between NO-deficient states and heightened IDO activity is needed. In contrast to our findings related to IDO-TMs and RV-PV dysfunction, we observed no significant relationship between levels of TH-TMs and RV-PV dysfunction. Serotonin (5-hydroxytryptamine) and the serotonin metabolite 5-hydroxyindoleacetic acid are tryptophan metabolites formed in pulmonary artery endothelial cells by TH (Figure 4A). These metabolites may contribute to pulmonary vasoconstriction and proliferation of pulmonary artery smooth muscle cells, but their storage in platelets and low abundance in plasma make them challenging and often unreliable to measure as PH biomarkers (2,34).
Plasma levels of TCA cycle intermediates, like IDO-TMs, were closely associated with hemodynamic indexes of RV-PV function in our study. Although the TCA cycle is operative within intracellular compartments, we have previously demonstrated that TCA cycle intermediates are detectable in plasma (35,36). Until recently, however, a potential role of TCA cycle intermediates in signaling biologic activity in a hormone-like role was unrecognized. A seminal study by He et al. (18) demonstrated that the span 1 TCA cycle intermediate succinate was sufficient to augment blood pressure in mice by activating a novel G protein–coupled receptor in the kidney. This group hypothesized that release of succinate in the setting of tissue ischemia may serve as a compensatory, protective mechanism to maintain systemic blood pressure. Recently, Zhao et al. (36) utilized pathological specimens from explanted or resected lungs and found elevated TCA cycle intermediates in lungs with severe PAH. Our study extends these findings by showing directionally concordant elevations in TCA cycle intermediates measured in plasma in patients with less advanced PV disease. Ours is the first study to relate circulating levels of succinate and other TCA cycle intermediates to PH, which raises the possibility that succinate or potentially other TCA-related metabolites could modulate vascular tone in humans.
There is mounting evidence of abnormal mitochondrial function, both within and outside the pulmonary circulation in PH. Beyond a well-recognized glycolytic shift that occurs in PH (37), down-regulation of the mitochondrial deacetylase sirtuin 3 has also been implicated in mitochondrial dysfunction in PH (38). The extent to which circulating TCA cycle or glycolytic intermediates accurately reflect mitochondrial dysfunction in PH will be an important area of future investigation.
We also found that metabolites that may relate to sympathetic nervous system activation tracked with RV-PV dysfunction (i.e., vanillylmandelic acid, cyclic adenosine monophosphate, and aminoisobutyric acid). These findings are consistent with previous studies relating heightened catecholamine levels to pulmonary vasoconstriction (17).
We used a liquid chromatography tandem mass spectrometry platform that provides key “sentinels” across known metabolic pathways but does not fully cover the entire plasma metabolome. Our study design focused on subjects with dyspnea on exertion in the main derivation and validation cohorts, as opposed to a narrowly defined patient population with PH, which resulted in heterogeneity in pathogenesis and degree of PH in our cohort. However, because PH purports a poor prognosis regardless of pathogenesis, biomarkers will be most useful in their application to patients with suspected, as opposed to confirmed, pulmonary vascular disease. Thus, our focus on dyspneic patients is advantageous in this regard. Identification of patients with clear evidence of pulmonary vascular disease in our second validation cohort demonstrated that these biomarkers performed particularly well in separating patients with clear evidence of PAH from those without RV-PV dysfunction (Figure 5). Evaluation of metabolite profiles during exercise warrants future study as a potential way to further improve the specificity of metabolite signatures for impaired RV-PV reserve capacity.
Accounting for potential confounders is of critical importance in a biomarker discovery study. Renal dysfunction in particular is associated with increases in several circulating metabolites, including tryptophan metabolites that are known to be increased in the setting of uremia (39). We excluded patients with advanced kidney disease, and average creatinine levels were normal in both our derivation and validation cohorts (1.1 mg/dl in each). Furthermore, we adjusted for cystatin C concentrations and found that IDO-TM, but not other TMs such as indoxysulfate, remained strongly predictive of pulmonary hemodynamic measurements. Simultaneous peripheral venous and central venous sampling was not available to permit comparison of metabolite levels from the 2 sampling sites. The small sample size of the validation cohort limited our power to identify additional metabolites that may have differed between PAH and control subjects. Finally, in our animal model studies, we acknowledge that the cellular source of metabolites within the lung was not localized, and we did not include invasive hemodynamic measurements to document the degree of PH elicited by hypoxia. However, Fulton ratios were similar to those in other studies in which a >2-fold increment in PAP was obtained (33).
The present study demonstrated the potential value of integrating metabolomic and clinical data in human cohorts who underwent detailed phenotyping while identifying small-molecule reporters of RV-PV function. This study begins to address an unmet clinical need to develop noninvasive markers of RV-PV function and reserve capacity. Novel markers of RV-PV functional reserve may aid in earlier detection of RV-PV dysfunction and help to quantify the burden of RV-PV dysfunction in patients with cardiopulmonary diseases.
COMPETENCY IN MEDICAL KNOWLEDGE: In patients with exertional dyspnea, plasma levels of tryptophan breakdown products, tricarboxylic acid cycle intermediaries and purine metabolites correlate with the severity of pulmonary hypertension and right ventricular dysfunction.
TRANSLATIONAL OUTLOOK: Additional research is needed to characterize the pathophysiological pathways by which these and other metabolites contribute to pulmonary hypertension and dyspnea of hemodynamic origin. Circulating biochemical signatures may ultimately be useful for earlier patient diagnosis or to identify those who might benefit from an intervention.
For an expanded Methods section as well as additional tables and figures, please see the online version of this article.
This work was supported by National Institutes of Health (NIH) grants 1R01 HL098280 and R01 DK081572 (Bethesda, Maryland) and the Leducq Foundation. The project was also supported by Clinical and Translational Science Award No. UL1TR000445 (Dr. Hemnes and Dr. Pugh) from the National Center for Advancing Translational Sciences (Bethesda, Maryland). Its contents are solely the responsibility of the authors and do not necessarily represent official views of the National Center for Advancing Translational Sciences or the NIH. Dr. Lewis was supported by NIH grants K23HL091106, NIH R01 HL131029, and NIH R01 HL119154 (Bethesda, Maryland), American Heart Association (AHA) grant 15GPSGC24800006 (Dallas, Texas), and the Hassenfeld Clinical Scholar Award (Boston, Massachusetts). Dr. Pugh has served as an advisory board member for Gilead. Dr. Gerszten was supported by an AHA Established Investigator Award (Dallas, Texas). All other authors have reported that they have no relationships relevant to the contents of this paper to disclose.
- Abbreviations and Acronyms
- body mass index
- cardiopulmonary exercise test
- change in mean pulmonary arterial pressure relative to change in cardiac output
- indoleamine 2,3-dioxygenase
- indoleamine 2,3-dioxygenase–dependent tryptophan metabolite
- liquid chromatography mass spectrometry
- left ventricular
- left ventricular ejection fraction
- nitric oxide
- N-terminal pro–B-type natriuretic peptide
- pulmonary arterial hypertension
- pulmonary arterial pressure
- pulmonary arterial wedge pressure
- pulmonary hypertension
- pulmonary vascular
- pulmonary vascular resistance
- right atrial pressure
- ribonucleic acid
- right ventricular
- right ventricular ejection fraction
- right ventricular-pulmonary vascular
- tricarboxylic acid
- tryptophan hydroxylase
- tryptophan metabolite
- Received August 12, 2015.
- Revision received October 6, 2015.
- Accepted October 13, 2015.
- American College of Cardiology Foundation
- Shao Z.,
- Wang Z.,
- Shrestha K.,
- et al.
- Xu W.,
- Koeck T.,
- Lara A.R.,
- et al.
- Takeyama D.,
- Kagaya Y.,
- Yamane Y.,
- et al.
- Stasch J.P.,
- Pacher P.,
- Evgenov O.V.
- Lewis G.D.,
- Shah R.,
- Shahzad K.,
- et al.
- Lewis G.D.,
- Bossone E.,
- Naeije R.,
- et al.
- Lewis G.D.,
- Lachmann J.,
- Camuso J.,
- et al.
- Lewis G.D.,
- Farrell L.,
- Wood M.J.,
- et al.
- Guazzi M.D.,
- De Cesare N.,
- Fiorentini C.,
- et al.
- Guazzi M.D.,
- Alimento M.,
- Fiorentini C.,
- Pepi M.,
- Polese A.
- Galiè N.,
- Hoeper M.M.,
- Humbert M.,
- et al.
- Wensel R.,
- Opitz C.F.,
- Anker S.D.,
- et al.
- Inuzuka R.,
- Diller G.P.,
- Borgia F.,
- et al.
- Saggar R.,
- Lewis G.D.,
- Systrom D.M.,
- Champion H.C.,
- Naeije R.,
- Saggar R.
- Kovacs G.,
- Berghold A.,
- Scheidl S.,
- Olschewski H.
- Lewis G.D.,
- Murphy R.M.,
- Shah R.V.,
- et al.
- Thomas S.R.,
- Mohr D.,
- Stocker R.