Abstract
Quantifying metabolic derangements in pulmonary hypertension (PH) by plasma metabolomics could identify biomarkers useful for diagnosis and treatment. The objective of this paper is to test the hypotheses that circulating metabolites are differentially expressed in PH patients compared with controls and among different hemodynamic subtypes of PH associated with left heart disease. We studied patients enrolled in the CATHGEN biorepository with PH (right heart catheterization mPAP ≥ 25 mmHg; n = 280). Of these, 133 met criteria for postcapillary PH, 82 for combined precapillary and postcapillary PH (CpcPH), and 65 for precapillary PH. Targeted profiling of 63 metabolites (acylcarnitines, amino acids, and ketones) was performed using tandem flow injection mass spectrometry. Multivariable linear regression was used to determine differences in metabolite factors derived from a principal components analysis between PH cases, PH subtypes, and non-PH controls. In adjusted models, the metabolite factor loaded with long-chain acylcarnitines was higher in all PH cases versus non-PH controls (P = 0.00008), but did not discriminate between CpcPH and postcapillary PH (P = 0.56). In analyses of subtypes, CpcPH patients had lower levels of factors loaded with urea cycle amino acids and short chain acylcarnitines as compared to controls (P = 0.002 and P = 0.01, respectively) and as compared to postcapillary PH (P = 0.04 and P = 0.02, respectively). Compared to controls, PH was strongly associated with greater concentrations of long-chain acylcarnitines. Postcapillary PH and CpcPH were weakly associated with distinct metabolomic profiles. These findings suggest the presence of unique metabolic abnormalities in subtypes of PH and may reflect underlying pathophysiology.
Additionally, PH in left heart disease is further divided into two broad hemodynamic subtypes that are defined by different levels of precapillary (pulmonary arterial) and postcapillary (pulmonary venous) PH (Table 1). As a whole, PH is defined by a mean pulmonary artery pressure (mPAP) ≥ 25 mmHg. 2 An abnormal pulmonary capillary wedge pressure differentiates PH in left heart disease (PCWP > 15 mmHg) from precapillary PH (PCWP ≤ 15 mmHg). In “passive” or postcapillary PH, elevated pressures result purely from postcapillary (pulmonary venous) congestion from elevated left side pressures, with no precapillary component of PH (pulmonary vascular resistance [PVR] of less than 3 Woods units [WU]). However, patients occasionally develop both precapillary and postcapillary PH, where the PCWP > 15 mmHg and PVR ≥ 3 WU. This combined precapillary and postcapillary PH—defined subsequently as “CpcPH”—can develop in the setting of prolonged or severe postcapillary pressure elevations. 3
While postcapillary PH may occur in any pathology raising left atrial pressure, the prevalence and severity of CpcPH is highly variable. Factors leading to pulmonary vascular remodeling depend on individual genetic, biologic, and environmental influences. Pulmonary vascular remodeling in CpcPH exhibits features seen in PAH, including intimal fibrosis and medial hypertrophy. 4 The management of patients with CpcPH is unclear: 2 while CpcPH has features of PAH, it is unknown whether using therapies targeting pathogenic pathways in PAH is beneficial. Identifying reliable biomarkers that could distinguish between postcapillary PH and CpcPH could lead to earlier diagnosis and give insight to novel therapeutic pathways.
Abnormalities in energy metabolism are common features in experimental models of both right and left ventricular failure. 5 In both failing left and right hearts, as mitochondrial metabolism is actively suppressed,6,7 a shift occurs from the highly efficient fatty acid oxidation to glucose metabolism and glycolysis. The resultant increase in lactate production and acidosis further impairs cardiac function. Advances in metabolomic profiling have enabled more comprehensive snapshots of dysregulated metabolism; in fact, metabolomic profiling has already been used to identify circulating markers of this impaired energy metabolism in coronary artery disease and heart failure.8–10 Similar profiling may increase understanding of PH and especially PH associated with left heart disease and provide long-sought biomarkers that provide earlier diagnosis of PH and discriminate CpcPH and postcapillary PH. Thus, we hypothesized that circulating metabolites reflecting mitochondrial energetics—including amino acids, acylcarnitines, and ketones—would discriminate PH from non-PH controls, subtypes of PH associated with left heart disease from non-PH controls, and between postcapillary PH and CpcPH in a discovery cohort.
Methods
Study sample
The CATHGEN biorepository consists of biological samples collected on 9300 sequential consenting individuals undergoing cardiac catheterization at Duke University Medical Center (Durham, NC, USA) between 2001 and 2010. 11 The majority of the participants enrolled in the biorepository were undergoing a catheterization for the evaluation of ischemic cardiovascular disease. At the time of arterial access for cardiac catheterization, blood was obtained from the femoral artery, immediately processed to separate plasma, and frozen at −80℃ for later use. All patients had fasted at least 6 h prior to sample collection. The Duke Databank for Cardiovascular Disease provided data on demographics, medical history, angiographic findings, and longitudinal follow-up. For this study, for PH cases, we identified study participants who had undergone right heart catheterization which showed a mPAP ≥ 25 mmHg, and who also had targeted metabolomic profiling data available (n = 280). Of these, postcapillary PH was defined as PCWP > 15 mmHg, PVR < 3 WU; CpcPH as PCWP > 15 mmHg, PVR ≥ 3 WU; and precapillary PH as PCWP ≤ 15 mmHg (Table 1). A cohort of 280 participants matched to reflect a similar distribution of clinical variables—including age deciles, race, sex, left ventricular ejection fraction, presence of coronary artery disease—and without PH based on echocardiogram estimated RVSP < 35 mmHg served as controls. The protocol for CATHGEN was approved by the Duke Institutional Review Board and all participants provided written informed consent for future analysis of stored samples.
Metabolite profiling
As previously described, we used targeted, quantitative tandem flow injection mass spectrometry (MS/MS) to determine levels of 60 metabolites (45 acylcarnitines and 15 amino acids) and standard assays for measurement of ketones, β-hydroxybutyrate, and non-esterified fatty acids.12,13 For MS/MS analyses, plasma samples were spiked with a cocktail of heavy-isotope internal standards (Cambridge Isotope Laboratories, MA, USA; CDN Isotopes, Canada) and deproteinated with methanol. The aliquoted supernatant was dried and esterified with hot, acidic methanol (acylcarnities) or n-butanol (amino acids). Tandem MS/MS using a Waters TQD (Milford, MA, USA) was employed to quantitatively assess acylcarnitine and amino acid ester concentrations. Total ketones (KET), 3-hydroxybutyrate (HBUT), and total non-esterified fatty acids (NEFA) were measured on a Beckman DxC600 autoanalyzer (Brea, CA, USA) using reagents from Wako (Richmond, VA, USA). All assays were tested in random batch order by the Metabolomics Core Laboratory at the Duke Molecular Physiology Institute, Duke University, Durham, NC, USA. Laboratory personnel were blinded to the clinical status of patients.
Statistical analysis
Given the large number of metabolites, principal components analysis (PCA) with varimax rotation was used in the larger CATHGEN cohort with metabolomics data (n = 3899) to reduce the large number of co-linear metabolites into 14 uncorrelated components9,13 (Table 4). As participants in this study were taken from the larger PCA, scoring coefficients from the original analysis were used to calculate weighted factor scores for each individual in this PH substudy. Metabolites with an absolute factor load ≥0.4 are listed as composing a given factor. To test the significance of differences in metabolite factors between overall PH case/control status and between PH subtypes, we constructed generalized multivariable linear regression models for each factor, adjusted for relevant clinical variables (age, presence/absence of coronary artery disease, sex, race, diabetes, and left ventricular ejection fraction) and metabolite assay batch number. Factor 5 (ketone-related analytes) was also adjusted for preoperative heparin status, as heparin is known to affect β-hydroxybutyrate, ketone, fatty acid, and triglyceride levels. 14 Separate regression models were used to test for mean metabolite factor levels between: (1) all PH participants and non-PH controls (primary analysis); (2) left heart disease-associated PH subgroups (postcapillary PH and CpcPH participants) each with non-PH controls (secondary analysis); and (3) postcapillary PH participants and CpcPH participants (secondary analysis). As these analyses were exploratory in nature and given co-linearity of the metabolites, we considered two-sided P values with significance defined as P ≤ 0.05. Although results are presented in the context of Bonferroni correction for multiple comparisons for the primary analysis, secondary analyses are presented without this adjustment. Statistical analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA).
Results
Baseline characteristics
Baseline patient characteristics.
Values are mean or percentages.

Distributions of hemodynamics in PH subgroups for (a) mean pulmonary artery pressure (mPAP, mmHg), (b) pulmonary capillary wedge pressure (PCWP, mmHg), and (c) pulmonary vascular resistance (PVR, Woods units). Due to the cutoffs used to define precapillary, postcapillary, and CpcPH, there are skewed distributions of the PCWP and PVR.
Baseline patient hemodynamics.
Values are mean (standard deviation).
Principal components analysis.
Included in the table are 12 factors identified from PCA. Column 2 lists the overall description of the factor, column 3 the individual metabolites comprising that factor. ALA, alanine; ARG, arginine; ASX, aspartic acid/asparagine; CIT, citrulline; GLX, glutamine/glutamate; GLY, glycine; HBUT, beta-hydroxybutyrate; HIS, histidine; KET, ketones; LEU/ILE, leucine/isoleucine; MET, methionine; NEFA, non-esterified fatty acids; ORN, ornithine; PHE, phenylalanine; PRO, proline; SER, serine; TYR, tyrosine; VAL, valine.
Metabolomic profiling of PH cases versus no-PH controls
Mean values for metabolite factor levels, stratified by PH subtype.
Values are mean (standard deviation) and are unitless.
*P value < 0.05 in comparison of listed categories after adjustment for age, race, sex, diabetes, ejection fraction, number of diseased coronary arteries, and batch number.
†Correction for multiple comparisons, with 14 components, is P < 0.0036.
Metabolomic profiling of PH-LHD subtypes
In comparing PH associated with left heart disease subtypes to no-PH controls, the most significant results revealed that both postcapillary PH and CpcPH groups individually continued to show higher mean levels of factor 4 than no-PH controls (long-chain acylcarnitines, P = 0.00001 and P = 0.002, respectively). Individuals with postcapillary PH also had lower mean levels of factor 5 (ketone-related metabolites, P = 0.04) and factor 14 (C22 acylcarnitine, P = 0.01). Individuals with CpcPH also demonstrated higher factor 2 levels (long-chain dicarboxyacylcarnitines, P = 0.01) and lower factor 8 (urea cycle amino acids, P = 0.005), factor 9 (C3-C5 acylcarnitines, P = 0.01), factor 12 (glycine and valine, P = 0.04), and factor 14 (C22 acylcarnitine, P = 0.04). Factor 4, however, was not significantly different between CpcPH and postcapillary PH patients. Instead, when compared to postcapillary PH patients, CpcPH patients exhibited higher factor 5 levels (ketone-related metabolites, P = 0.04), but lower factor 8 (urea cycle amino acids, P = 0.04) and factor 9 (C3-C5 acylcarnitines, P = 0.02) levels.
Discussion
In this study, we hypothesized that metabolite profiles reporting on mitochondrial energetics would differentiate PH from controls without PH, and different hemodynamic subtypes of PH associated with left heart disease from one another. Our most significant finding was for higher levels of circulating long-chain acylcarnitines (factor 4) in patients with PH as compared with controls without PH even after adjustment for many clinical variables, a finding that met rigorous adjustment for multiple comparisons overall and was significant in comparisons of each PH subtype with controls. These results suggest a shared mitochondrial dysfunction common to all types of PH. Our second most significant finding was for urea cycle metabolites (factor 8), with lower levels in PH patients than no-PH controls, and differentiating PH subtypes in left heart disease with CpcPH having the lowest levels overall.
While the metabolomics changes that we found were limited to a discovery cohort, they are broadly consistent with other findings that have been described in the literature. The non-failing heart primarily utilizes long-chain fatty acids as an energy source, and an early step in their oxidation in the mitochondria is the formation of long-chain acylcarnitines.
15
In our study, the increased levels of long-chain acylcarnitines (factor 4, Fig. 2) seen across the PH subtypes as compared with no-PH potentially reflects the inefficiencies in mitochondrial fuel metabolism previously described in both left and right ventricular failure. Recently, these abnormalities in fatty acid metabolism have been identified in the plasma and RV in patients with PAH.
16
Similar abnormalities in long-chain acylcarnitine levels have been noted in left heart failure including heart failure with preserved ejection fraction and improve with ventricular assist device therapy.8,10 These findings are consistent with abnormalities in metabolic reprogramming in animal models of heart failure.
17
This suggests that abnormalities in long-chain acylcarnitines reflect a common signature associated with myocardial mitochondrial dysfunction.
Normalized concentrations of long-chain acylcarnitines in PH hemodynamic subgroups. Shown are values normalized to control. *P < 0.05 by two-way ANOVA with Tukey test for multiple comparisons.
The other nominally significant results included results for factor 2, comprised of long-chain dicarboxylacylcarnitines, unique in being the only factor which was only different in CpcPH compared to controls. Dicarboxylated acylcarnitines are derivatives of mitochondrial lipid or amino acid oxidation and a product of omega-oxidation by microsomal cytochrome p450 enzymes in peroxisomes. We have previously described a strong association of similar dicarboxyl acylcarnitines with cardiovascular events in a cardiovascular disease cohort. 21 The precise substrates and pathways that generate the dicarboxylated acylcarnitines and factors that foster their production require further investigation.
Urea cycle amino acids (factor 8) were lower in CpcPH compared to postcapillary PH patients as well as controls. This factor is strongly weighted by glycine and may be reporting on subsequent degradation to ammonia and entry into the urea cycle. Abnormalities in these metabolites has been noted in human pulmonary endothelium with BMPR2 mutations. 18 Recent data have linked increased glycine levels and decreased urea cycle metabolites with a chromosomal locus and identified a strong association with decreased risk of atherosclerosis. 19 Glycine may have cardioprotective properties against coronary arterial endothelial cell inflammation, further supporting the hypothesis that genetic susceptibility can lead to a differential response to cardiopulmonary vascular inflammation and stress. 20
A number of other metabolic abnormalities have been reported in pulmonary vascular disease and right heart failure. In a recent study using an alternative targeted metabolomics platform, differences in indoleamine 2,3-dioxygenase-dependent tryptophan metabolites, tricarboxylic acid intermediates, purine metabolites, and arginine-nitric oxide metabolic pathway constituents were associated with right ventricular–pulmonary vascular dysfunction. 22 In future studies, combining multiple biomarkers reflecting the many pathways dysregulated in pulmonary vascular disease—including the nitric oxide pathway, mitochondrial bioenergetics, and vasoactive mediators—could provide additional insight into left heart disease associated PH.
Our results simultaneously highlight metabolic pathways dysregulated in PH while identifying circulating, easily measurable biomarkers that distinguish PH from no-PH and PH subtypes associated with left heart disease. There are some limitations to our study, however. The cross-sectional design of our analysis allows only for an examination of associations, not causation. Given the exploratory nature of this small sample size and multiple comparisons made, we may be subject to both false-positive and false-negative findings. The use of PCA mitigates the number of multiple comparisons by reducing a large number of metabolites into unrelated groupings in an unbiased manner, and our primary finding for long-chain acylcarnitines survives adjustment for a conservative correction for multiple comparisons. While we adjusted for potential confounders such as age, race, diabetes, left ventricular ejection fraction, and coronary artery disease, there may be unmeasured confounders. Also, as a result of the inclusion criteria to the CATHGEN biorepository, patients with precapillary PH had very mild disease by hemodynamic criteria. Thus, our findings do not likely extend to those with severe PAH. Finally, while some results were nominally significant for factors with lower eigenvalues (e.g. factor 12 which was nominally significant in comparisons of PH to no-PH and CpcPH to no-PH) likely are false positives representing statistical noise of higher number factors that explain a much smaller proportion of variance in the dataset.
Conclusions
We observed unique metabolite profiles in PH compared to controls and in different hemodynamic subtypes of PH seen in left heart disease. These metabolites likely reflect changes in mitochondrial energetics associated with pulmonary vascular disease and heart failure, but could also report on global metabolic disturbances in this disease state.
Footnotes
Conflict of interest
The author(s) declare that there is no conflict of interest.
Funding
This project was supported in part by grant number HL095987 from the National Institutes of Health (Bethesda, MD, USA). SR is supported by HL114643 and a Burroughs Wellcome Career Award for Medical Scientists (Durham, NC, USA).
