Abstract
To assess the relationship of cytokines with functional and clinical outcomes in pulmonary arterial hypertension (PAH). Endothelial dysfunction and vascular inflammation are characteristic of PAH. We investigated whether markers of angiogenesis and inflammation associated with functional, hemodynamic parameters, and clinical outcomes in PAH. PAH patients (n = 206) were pooled from two clinical trials: TRUST-1 and FREEDOM-C2. Baseline and post-treatment cytokine levels were correlated to baseline clinical and hemodynamic parameters, were assessed in clinical subgroups, and were associated with clinical outcomes. In 206 patients (mean age = 48 years; 74% women) with WHO group-1 PAH, most cytokine levels were higher in those with 6-min walking distance (6MWD) < median (335 m) vs. those above median, including Ang-1 (11.9 ± 10.1 vs. 5.9 ± 6.0 ng/mL), Ang-2 (14.3 ± 11.8 vs. 12.2 ± 11.2 ng/mL), and MMP-9 (221 ± 262.3 vs. 119 ± 171 ng/mL). Baseline 6MWD inversely correlated with Ang-1 (r = −0.27, P < 0.0001), Ang-2 (r = −0.20, P = 0.004), and MMP-9 (r = −0.27, P < 0.0001). MMP-9 levels differed significantly by NYHA functional class (P = 0.001) suggesting an association between MMP-9 and subjective PAH severity. Mean Ang-2 levels were higher in those with baseline right atrial pressure (RAP) > 15 mmHg compared to those with RAP < 15 mmHg (23,841 vs. 11,020 pg/mL). Baseline RAP was associated with change in MMP-9 levels (r = −0.53, P = 0.03). Finally, baseline Ang-1, VEGF and MMP-9 levels were associated with risk of death and hospitalization at 16-week follow-up. Inflammatory cytokines and vascular angiogenesis markers are associated with baseline functional, hemodynamic parameters in PAH, and predict death and hospitalization. Larger prospective studies are needed to confirm the utility of cytokines in PAH.
Pulmonary arterial hypertension (PAH) is a progressive, debilitating, and often fatal disease characterized by a variety of physiologic changes in the distal pulmonary vasculature, ultimately resulting in a rise in pulmonary vascular resistance (PVR) and right ventricular afterload.1,2 Smooth muscle cell proliferation, endothelial dysfunction, and vascular inflammation are some of the underlying mechanisms thought to be responsible for the pathogenesis of PAH. 3 In recent years, greater attention has been focused on vascular inflammation in patients with PAH, ranging from idiopathic PAH (IPAH) to connective tissue disease-associated PAH.4,5
Three major molecular pathways have already been implicated in the pathogenesis of PAH specifically involving endothelin-1, nitric oxide, and prostacyclin. 6 All currently available PAH-specific therapies target one of these pathways. 7 Among the available therapeutic agents for the treatment of PAH, prostacyclin analogs, such as Epoprostenol and Treprostinil, are potent agents that have proven efficacy with regard to symptomatic improvement and, in the case of intravenous Epoprostenol, a demonstrated survival benefit.8–11 As novel treatments for this incurable disease are being developed, surrogate markers to evaluate treatment response are required; these may in turn shed light on the pathophysiology of the disease itself. 12 Emerging evidence suggests the potential for cytokine signatures as important circulating biomarkers in PAH pathology as well a means to measure treatment response.13–16
In this study, we sought to investigate the association between circulating markers of vascular angiogenesis and inflammatory cytokine levels with PAH severity at baseline and in response to treatment with a prostacyclin analog. To achieve this, we utilized the treatment arms of two previously published randomized trials of Treprostinil in PAH, TRUST-1 17 and FREEDOM-C2. 18 We hypothesized that serum markers of vascular angiogenesis, such as angiopoietin-2, and of vascular inflammation, such as IL-8, along with matrix metalloproteinases, such as matrix metalloproteinase-9 (MMP-9), would be associated with baseline clinical, hemodynamic, and functional parameters in PAH. Moreover, we sought to investigate the effect of oral and intravenous treprostinil therapy on the levels of these biomarkers. Finally, we examined whether baseline levels of these markers might be associated with clinical outcomes such as death and hospitalization after PAH therapy with oral and intravenous treprostinil.
Methods
Data collection
Data were pooled from the treatment arms of two clinical trials assessing the effect of treprostinil on PAH, TRUST-1 17 and FREEDOM-C2. 18 Specific details of the inclusion and exclusion criteria as well as the study results have been previously reported. Baseline data from 206 patients (28 patients from TRUST-1 and 178 from FREEDOM-C2) were pooled while follow-up data were available for a total of 193 patients (16 from TRUST-1 and 177 from FREEDOM-C2).
In brief, the TRUST-1 trial comprised 28 WHO Group I PAH patients who were evaluated at baseline and after 12 weeks of treatment with IV treprostinil with 6-min walking distance (6MWD), New York Heart Association (NYHA) functional class (FC) assessment, and N-terminal pro-B type natriuretic peptide (NT-proBNP). These patients all had invasive hemodynamic assessment at baseline before starting IV treprostinil. Similarly, data from 178 patients enrolled in the treatment arm of the FREEDOM-C2 trial were available for analysis. These patients underwent clinical assessment for 6MWD, NYHA FC, and NT-proBNP at baseline and at 16-week follow-up visits after treatment with oral treprostinil.
Cytokine levels were measured for all patients at baseline and when possible at 12- or 16-week follow-up (TRUST-1 and FREEDOM C-2, respectively). Biomarkers were assessed in EDTA plasma using a multiplexed immunoassay (DiscoveryMAP® v. 3.0 assay, Myriad RBM, Inc., Austin, TX, USA). Previous publications of these data have included detailed information regarding biomarker measurements and the assays used.17,18 Pre-specified cytokines that were measured included angiopoietin-1 (Ang-1), angiopoietin-2 (Ang-2), basic fibroblast growth factor (BFGF), interleukin-1 beta (IL-1b), IL-6, IL-8, IL-13, MMP-9, and vascular endothelial growth factor (VEGF). In addition, in a subset of 95 patients, levels of NT-proBNP levels were also measured.
Statistical analysis
Descriptive analyses were performed to summarize the baseline characteristics of the combined group. Normality of the data was assessed by considering skewness and kurtosis. Continuous variables were reported as mean ± SD or median and range for parametric and non-parametric descriptive summaries, respectively. Categorical variables were represented as frequencies. Stratified analyses based on baseline 6MWD, baseline NYHA FC, death, and hospitalization were performed, in which differences in continuous variables were assessed using the Student’s t-test for normally distributed data whereas non-normal data were compared between stratified groups utilizing the Mann–Whitney U test. Categorical variable comparisons were performed utilizing Pearson’s chi-square test. The Holm-Bonferroni step-down procedure to test the significance of multiple comparisons of analysis of biomarkers was performed. Furthermore, Spearman correlation analyses were carried out to examine correlations between baseline cytokine levels and baseline 6MWD as well as between change in cytokine levels and baseline 6MWD. ANOVA models were built to evaluate associations between cytokines and NYHA FC. Finally, Kaplan–Maier analysis and log-rank testing was used to determine the association between baseline cytokine levels and subsequent composite outcome including death and/or hospitalization. P values ≤ 0.05 were considered statistically significant. All analyses were carried out using SAS statistical software, version 9.4.
Results
Characteristics of the overall study group
Baseline characteristics of the combined patient population (n = 206).
Values are represented as mean ± SD or median (range) for continuous variables and as n (%) for categorical variables.
HIV, human immunodeficiency virus; ERA, endothelin receptor antagonist; PDE-5i, phosphodiesterase-5 inhibitor.
Stratified analyses of baseline characteristics
Baseline cytokine profile comparison between groups stratified by 6-min walk distance (6MWD).
Between-group comparisons were performed by Mann–Whitney U-test.
Baseline cytokine profile comparison between patients who died during the study and patients who were alive at the end of the study.
Between-group comparisons were performed by Mann–Whitney U-test.
Baseline cytokine profile comparison between patients who were hospitalized during study and those who were not.
Between-group comparisons were performed by Mann–Whitney U-test.
Baseline cytokine profile comparison between patients with baseline NYHA FC 2 and those with NYHA FC 3/4.
Between-group comparisons were performed by Mann–Whitney U-test.
Spearman correlation reveals significant inverse associations between pro-inflammatory cytokine levels and baseline 6MWD
Spearman correlation analyses between cytokines and 6-min walk distance (6MWD).
Correlation between baseline hemodynamic parameters and cytokine levels
Among patients enrolled in TRUST-1, both Ang-2 and IL-8 levels were associated with baseline RA pressure > 15 mmHg vs. those with RA pressure < 15 mmHg (23,841 pg/mL vs. 11,020 pg/mL, P = 0.065) and (229 pg/mL vs. 14 pg/mL, P = 0.03), respectively. Conversely, baseline cardiac index (CI) (r = 0.69, P = 0.005) and PVR index (PVRI) (R = −0.68, P = 0.005) were associated with change from baseline Ang-2 levels with IV treprostinil therapy. Furthermore, baseline RA pressure was associated with change from baseline MMP-9 levels (R = −0.53, P = 0.035). Additionally, the change from baseline MMP-9 levels was associated with change from baseline CI and PVRI after IV treprostinil therapy (R = −0.57, P = 0.054 and R = 0.61, P = 0.036, respectively).
Correlation between baseline cytokine levels and NT-ProBNP
In a subset of 95 patients, NT-proBNP levels were measured at baseline and correlated to baseline cytokine levels. From this sub-dataset, there was a significant positive correlation between NT-proBNP and Ang-2 at baseline (R = 0.67, P < 0.0001). There was also significant positive correlation noted between baseline levels of IL8 and NT-proBNP (R = 0.30, P = 0.0027) as well as baseline levels of MMP-9 and NT-proBNP (R = 0.26, P = 0.012).
ANOVA model demonstrates significant dose response in MMP-9 levels in association with worsening NYHA class
ANOVA model derived significance levels for comparisons of cytokine levels among NYHA FCs.
Differentiation between cytokines in patients with collagen vascular disease versus IPAH
Baseline cytokine levels in patients with IPAH vs. those with collagen vascular disease-related PAH.
IPAH, idiopathic/familial pulmonary arterial hypertension; CVD PAH, collagen vascular disease-related pulmonary arterial hypertension.
Baseline cytokine levels are associated with composite outcome of death and hospitalization
Finally, we assessed the association between baseline cytokine levels and the composite outcome of death and/or hospitalization in the two studies using Kaplan–Maier analysis (Fig. 1). As such, we found that higher than median levels of Ang-1 (log rank P < 0.001), VEGF (log rank P = 0.001), and MMP-9 (log rank P < 0.001) were associated with higher rates of hospitalizations and death. However, Ang-2 failed to show association with these outcomes.
Kaplan–Maier curves showing association between baseline Ang-1 levels (a), VEGF levels (b), mmp9 levels (c), and composite of death and hospitalization at 16 weeks.
Discussion
In this study, pooling data from the treatment arms of two randomized clinical trials of treprostinil therapy in PAH, we presented several important findings: (1) in cross-sectional analysis, we showed that increased levels of these cytokine signatures were associated with lower 6MWD, hospitalizations, and death; (2) further, we demonstrated that baseline 6MWD correlated inversely with multiple cytokines potentially involved in the pathophysiology of PAH; (3) additionally, we demonstrated an association between markers of inflammation including Ang-2 and MMP-9 and hemodynamic characteristics in patients with PAH, and revealed a direct relationship between baseline MMP-9 levels and NYHA FC in PAH patients; and (4) we demonstrated that baseline levels of pro-inflammatory cytokines were associated with subsequent risk of death and hospitalizations over a short-term follow-up of 16 weeks. Collectively, these findings demonstrate the potential utility of these cytokines as biological markers in PAH patients which may correlate with functional capacity and PAH-related hospitalizations and deaths.
PAH is an often-fatal disease characterized by increased PVR associated with endothelial dysfunction and vascular inflammation. 3 The three known pathways implicated in the pathophysiology of PAH are endothelin-1, nitric oxide, and prostacyclin. With the progress in our understanding of the complex interplay at the molecular between these pathways and the structural impact of this interplay, resulting in increases in PA pressure leading to PAH, research investigating the role of pro-inflammatory cytokines is gaining greater attention.13,15,16 In view of these recent advances, our findings that elevated levels of markers of vascular angiogenesis and pro-inflammatory cytokines are associated with both surrogate outcome measures such as 6MWD and clinical outcomes such hospitalizations are of importance. Together, these results suggest that these cytokines may have clinical and subclinical implications as biomarkers of baseline assessment, gauging treatment response, and may suggest future outcomes in patients with PAH.
As the molecular pathophysiology of PAH becomes further elucidated, more surrogate outcomes to evaluate disease severity and treatment response are required. 10 We demonstrated that Ang-2 levels were inversely associated with baseline 6MWD. Building upon the previous evidence suggesting Ang-2 as a biomarker for PAH,17,19 we showed further potential for Ang-2 in a larger sample set. Additionally, matrix metalloproteinases are zinc-dependent, multi-functional, soluble, or cell-membrane anchored proteases. 20 Out of this family of proteolytic enzymes, MMP-9 may play a crucial role in the initiation, development, and progression of PAH. Specifically, MMP-9 plays an important role in endothelial cell migration, smooth muscle cell migration and proliferation, adventitial fibroblast trans-differentiation, and inflammatory cell recruitment and infiltration. 20 Through its multitude of effects, MMP-9 facilitates the process of pulmonary vascular remodeling thereby leading to PAH. Our study revealed that MMP-9 was the only biomarker associated with baseline NYHA FC, further supporting its potential role as a composite biomarker of PAH.
Invasive hemodynamic assessment via right heart catheterization is the gold standard for diagnosis of pulmonary hypertension (PH). As such, using right-sided invasive hemodynamic measures, we demonstrated an association of Ang-2, IL-8, and MMP-9 with CI as well as PVRI. Furthermore, utilizing a composite outcome of death and hospitalization, we established predictive values for Ang-1, VEGF, and MMP-9. Taken together, these results suggest that measurement of these inflammatory cytokines at baseline may be useful to incorporate into future risk stratification algorithms to assess the severity of PAH and risk for future hospitalization and death. However, our observations with regards to associations between cytokines and hospitalization/death should be cautiously interpreted since they are over a short time period and obviously are non-randomized. However, they may provide a “proof of concept” and should be evaluated further in larger studies with longer duration of follow-ups.
Recently, a study analyzing the role of circulating Ang-1 as a biomarker of disease severity in PAH 21 suggested that Ang-1 was not reliable. However, this study had patients with PH secondary to multiple underlying pathophysiological conditions, rendering the results less generalizable. While the results from our study are somewhat at odds with these recently published data, we had a larger sample size, all of whom had WHO Group I PAH. Nonetheless, since both these studies were cross-sectional and observational in nature, larger population-based studies are required to elucidate the role of Ang-1 in detail.
Strengths and limitations
One of the strengths of this study is that it included a large dataset of rigorously characterized WHO Group I patients who were assessed with biomarker assessment at baseline and at a follow-up after treprostinil therapy. Regardless of the size of this study, the nature of these analyses remains cross-sectional and observational, leaving the possibility for residual confounding, which is a limitation of all observational studies. Furthermore, despite including over 200 WHO Group I PAH patients, our study sample size is still limited and future population-based studies are necessary to examine the role of these inflammatory biomarkers in PAH pathophysiology. Though we did have access to baseline hemodynamic data from TRUST-1, neither trial included right heart imaging data with echocardiography or cardiac MRI and this we could not correlate cytokine levels with right ventricular function. Finally, given the cross-sectional nature of most of our study analyses, we cannot evaluate causality and hence future research should focus on studying the mechanistic roles of these cytokines in PAH. However, despite these limitations, our study is the largest study to date to assess the role of inflammatory cytokines and markers of angiogenesis in PAH and to inspect their relationships with the characteristics of PAH disease severity in a systematic manner.
In summary, utilizing the pooled biomarker data from two PAH clinical trials, we conclude that the cytokine levels and markers of vascular angiogenesis are associated with baseline 6MWD in patients with PAH and with death and hospitalizations. Moreover, we established an association between markers of inflammation and hemodynamic characteristics of the patients with PAH. Finally, we established that baseline levels of Ang-1 and MMP-9 helped assess future risk for hospitalization and deaths. As such, inflammatory cytokines may be useful biomarkers in PAH, but larger prospective studies are required to confirm the utility of these findings in future risk prediction algorithms.
Footnotes
Conflict of interest
AR receives research support and is a principal investigator with United Therapeutics. AR is also a principal investigator for Actelion and Bellerophon, a speaker for Bayer, and consultant for Sr. Jude. RB is a PI and has grant support from Bellerophon, Actelion, Eiger, St Jude, and Bayer. YR and KS are United Therapeutics employees.
Funding
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
