Characteristics of the study population
Baseline characteristics of treatment-naïve PAH patients at the time of diagnosis and healthy controls are summarized in Table 1. During a median follow-up of 42 [interquartile range: 32-58] months, twelve PAH patients (Non-survivors) reached the primary endpoint, the causes included end-stage heart failure (n=5), euthanasia because of end-stage cardiovascular and pulmonary disease (n=2), lung transplantation (n=1), progression of systemic sclerosis (n=1), multi-organ failure (n=1), sudden death presumed cerebral (n=1), and malignancy (n=1). Seven of them reached the endpoint within six months after diagnosis and hence had no follow-up measurements taken. Non-survivors were older, had higher heart rate and shorter 6-minute walking distance than survivors (Table 1).
Table 1
Baseline characteristics of all PAH patients and healthy controls
|
Control
|
PAH
|
|
All
|
Survivors
|
Non-survivors
|
N
|
111
|
43
|
31
|
12
|
Aetiology
|
|
|
|
|
- iPAH, n (%)
- CTD-PAH, n (%)
- CHD-PAH, n (%)
|
|
15 (35)
17 (40)
11 (25)
|
14 (45)
8 (26)
9 (29)
|
1 (8)
9 (75)
2 (17)
|
Clinical characteristics
|
|
|
|
|
Age, years old
|
43±13
|
53±17 **
|
49±14
|
62±19 †
|
Sex, women n (%)
|
59 (53)
|
29 (67)
|
20 (65)
|
9 (75)
|
sBP, mmHg
|
123 [115-128]
|
122 [114-132]
|
123 [116-132]
|
115 [106-132]
|
HR, beats/min
|
68 [62-76]
|
78 [67-90] **
|
76 [63-87]
|
87 [76-99] †
|
BMI, kg/m2
|
23.8±2.9
|
27.0±6.1 ***
|
27.6±6.6
|
25.2±4.5
|
eGFR, mL/min/1.73m2
|
---
|
66.8±18.4
|
68.9±16.9
|
61.3±21.8
|
hs-CRP, mg/L
|
---
|
3.1 [1.5-10.5]
|
3.1 [1.2-9.0]
|
3.4 [2.1-22.3]
|
NYHA, I:II:III:IV
|
---
|
1:13:23:6
|
1:11:16:3
|
0:2:7:3
|
6MWD, m
|
---
|
337±153
|
377±136
|
198±133 ††
|
Right heart catheterization
|
|
|
|
mPAP, mmHg
|
---
|
50.5±16.1
|
50.5±15.9
|
50.6±17.3
|
PAWP, mmHg
|
---
|
11.8±5.6
|
11.6±6.0
|
12.2±4.9
|
PVR, WU
|
---
|
7.1 [5.1-11.8]
|
7.9 [5.6-12.0]
|
6.4 [4.2-11.3]
|
CO, L/min
|
---
|
4.7 [3.9-5.5]
|
4.6 [3.9-5.5]
|
4.9 [3.9-6.5]
|
CI, L/min/m2
|
---
|
2.5 [2.2-3.3]
|
2.5 [2.2-3.2]
|
2.5 [2.0-3.6]
|
PAH therapy types at time of censoring
|
|
|
|
|
- No therapy, n (%)
- Mono-therapy, n (%)
- Dual-therapy, n (%)
- Triple-therapy, n (%)
|
|
4 (9.3)
9 (20.9)
13 (30.2)
17 (39.5)
|
2 (6.4)
4 (12.9)
11 (35.5)
14 (45.2)
|
2 (16.7)
5 (41.6)
2 (16.7)
3 (25.0)
|
Data are presented as mean ± SD, median [IQR], or numbers (percentages). **P<0.01, *** P <0.001 PAH versus control. †P<0.05, ††P<0.01, survivors versus non-survivors. Unpaired T Test, Mann-Whitney U Test, and Chi-Square. PAH: pulmonary arterial hypertension, iPAH: idiopathic PAH, CTD-PAH: connective tissue diseases associated PAH, CHD-PAH: congenital heart diseases associated PAH, sBP: systolic blood pressure, HR: heart rate, BMI: body mass index, eGFR: estimated glomerular filtration rate, hs-CRP: high-sensitivity C-reactive protein, NYHA: New York Heart Association classification, 6MWD: 6-minute walking distance, mPAP: mean pulmonary arterial pressure, PAWP: pulmonary arterial wedge pressure, PVR: pulmonary vascular resistance, CO: cardiac output, CI: cardiac index.
KP-metabolite profile in PAH
At the time of diagnosis (baseline), Trp was significantly lower in treatment-naïve PAH patients compared to controls, while Kyn, 3-HK, QA, KA and AA were significantly higher in treatment-naïve PAH patients and no significant difference in 3-HA was found (Figure 1B). Binary logistic regression analyses showed that KP-metabolites significantly distinguished PAH patients from controls at baseline both in a univariate model and when corrected for age, sex, and body mass index (Table 2). Moreover, manual stepwise logistic regression analyses showed that including the whole panel of altered metabolites in the model predicted better than only one metabolite by significantly increasing the Chi-square of the model (Table 2).
MCT-PH rats showed a similar KP-metabolite profile as that observed in PAH patients (Figure 1C), which was accompanied by an increase in NAD+ in the lungs from these rats (1.26 [0.81-1.46] in 6 healthy rats, vs 1.70 [1.43-1.92] µM/mg in 7 MCT-PH rats, P=0.035). Conversely, the KP-metabolite profile was different in SuHx-PH rats (Figure 1D). As the MCT-PH rats show the most severe inflammatory phenotype, these results suggest a link between inflammation and KP-metabolism in PAH.
Table 2
Prediction of PAH with each 1 µM decrease for Trp and 1 nM increase in other KP metabolites by binary logistic regression.
|
Odds Ratio [95% CI]
|
P value
|
Chi-Square
|
P value
|
Univariable
|
|
|
|
Trp
|
1.047 [1.024-1.070]
|
<0.001
|
21.083
|
<0.001
|
Kyn
|
1.003 [1.002-1.004]
|
<0.001
|
39.370
|
<0.001
|
3-HK
|
1.097 [1.051-1.144]
|
<0.001
|
38.390
|
<0.001
|
QA
|
1.053 [1.032-1.075]
|
<0.001
|
56.114
|
<0.001
|
KA
|
1.076 [1.023-1.132]
|
0.005
|
8.842
|
0.003
|
AA
|
1.924 [1.453-2.550]
|
<0.001
|
40.132
|
<0.001
|
Whole panel
|
|
|
92.492
|
<0.001
|
Adjusted for age, sex and body mass index
|
|
|
Trp
|
1.036 [1.011-1.060]
|
<0.001
|
35.551
|
<0.001
|
Kyn
|
1.003 [1.002-1.005]
|
<0.001
|
57.983
|
<0.001
|
3-HK
|
1.089 [1.043-1.137]
|
<0.001
|
53.224
|
<0.001
|
QA
|
1.047 [1.026-1.069]
|
<0.001
|
63.944
|
<0.001
|
KA
|
1.090 [1.029-1.155]
|
0.003
|
35.368
|
<0.001
|
AA
|
1.870 [1.398-2.502]
|
<0.001
|
54.705
|
<0.001
|
Whole panel
|
|
|
93.187
|
<0.001
|
Trp: Tryptophan, Kyn: Kynurenine, 3-HK: 3-Hydroxykynurenine, QA: Quinolinic acid, KA: Kynurenic acid, AA: Anthranilic acid.
IL-6/IL-6Rα contributed to the activation of KP-metabolism in vitro
Stimulation with IL-6/IL-6Rα complex induced KP activation and mimicked the KP-metabolite profile observed in PAH patients in all three cell types, while stimulation with IL-6 alone failed to induce a similar profile (Figure 2A&B). Moreover, the cells responded differently to other cytokines and hypoxia failed to induce a KP profile similar to that seen in PAH patients (Figure 2C-E). Taken together, these results indicate that inflammation, particularly activation of IL-6/IL-6Rα signaling contributed to the KP activation in PAH patients.
Effects of PAH therapy on KP-metabolism
Following six months of PAH therapy, Trp was still significantly lower in PAH patients compared to controls, while only Kyn, 3-HK, QA, and AA were still significantly higher in PAH patients (Figure 1E-J). After one year of PAH therapy, Trp was still significantly lower in PAH patients than in controls, while only 3-HK and QA were still significantly higher in PAH patients (Figure 1E-J). These data suggest KP-metabolite profile partly normalized in PAH patients after PAH therapy.
When baseline KP-metabolite levels were compared between survivors and non-survivors, only Kyn was significantly higher in non-survivors versus survivors (Figure 3D). However, when comparing these levels at the latest measurement available, Kyn, 3-HK, QA, KA, and AA were all significantly higher in non-survivors compared with survivors (Figure 3, right panels). In these survivors, Wilcoxon matched-pairs signed rank test showed that Kyn, 3-HK, QA, KA, and AA were significantly decreased after one year but not six months of PAH therapy, indicating that only long-term PAH therapy decreased KP-metabolite levels (Figure 3, left panels).
These data indicate that KP-metabolites could be potential predictors of response to PAH therapy in survivors, and regular monitoring of KP-metabolites may be important to evaluate clinical status of PAH patients.
Correlation of KP-metabolites and disease-severity
Significant correlations between different KP-metabolites were seen in healthy controls and PAH patients at baseline. In healthy controls, KP-metabolites correlated with each other, with one exception (lack of correlation between Trp and QA (Figure 4A)). In PAH patients, Trp was not correlated with any other metabolite, while other metabolites still correlated with each other at baseline as well as after PAH therapy for six months and one year (Figure 4B-D).
Significant correlations of KP-metabolite levels with baseline characteristics were observed both in healthy controls and PAH patients (Figure 4E&F). Importantly, in PAH patients, higher levels of Kyn, 3-HK, 3-HA, and KA correlated with higher pulmonary vascular resistance, while higher levels of Kyn and KA correlated with reduced cardiac output and cardiac index, higher levels of QA correlated with a reduced cardiac index. Higher levels of Kyn, 3-HK, QA, KA, and AA correlated with worse functional class, higher levels of Kyn correlated with shorter 6-minute walking distance (Figure 4G). In addition, all KP-metabolites correlated with estimated glomerular filtration rate, and Trp, 3-HK, QA, and AA correlated with high-sensitivity C-reactive protein (Figure 4G), further supporting the potential link between inflammation and KP-metabolism.
Survival analyses
PAH patients were stratified into two groups based on the median level of KP-metabolites measured at baseline. High levels of Kyn (>1.328 µM), 3-HK (>22.71 nM) or QA (>75.23 nM) predicted worse early survival (Breslow Test), while high Kyn levels also predicted worse long-term survival (Log-rank Test, Figure 5). However, there was no difference between patients with low and high levels of Trp, 3-HA, KA, or AA (Figure 5).
Since survivors had lower levels of KP-metabolites at the latest measurement timepoint, which may be associated with a better response to PAH therapy, we compared the survival curves in PAH patients based on the latest available measurement. Again, high levels of Kyn, 3-HK, or QA predicted worse early survival, but also worse long-term survival (Figure 5). In addition, patients with high levels of 3-HA (>8.833 nM), KA (>12.92 nM), or AA (>4.944 nM) had worse early survival, and patients with high levels of 3-HA, or KA also had worse long-term survival (Figure 5).
When considering KP-metabolites as continuous variables, for each 1nM increase in Kyn, 3-HK, QA, or KA was associated with an increased hazard ratio of death (Table 3).
These results indicate that elevations in KP-metabolites are potential predictors of survival for PAH patients, with Kyn being the strongest prognostic biomarker.
Table 3
Cox proportional hazard analyses for death by 1nM increase in KP-metabolite levels in PAH patients
|
Baseline measurement
|
Latest measurement
|
Univariable
|
Hazard ratio [95% CI]
|
P value
|
Hazard ratio [95% CI]
|
P value
|
Kyn
|
1.001 [1.000-1.002]
|
0.012
|
1.001 [1.000-1.002]
|
0.001
|
3-HK
|
1.020 [1.009-1.030]
|
<0.001
|
1.021 [1.010-1.032]
|
<0.001
|
QA
|
1.006 [1.002-1.010]
|
0.001
|
1.007 [1.003-1.011]
|
0.001
|
KA
|
1.045 [0.996-1.096]
|
0.074
|
1.054 [1.007-1.104]
|
0.025
|
3-HA
|
1.119 [1.001-1.251]
|
0.048
|
1.114 [0.992-1.250]
|
0.067
|
AA
|
1.015 [0.906-1.138]
|
0.798
|
1.085 [0.987-1.193]
|
0.093
|
Kyn: kynurenine, 3-HK: 3-hydroxy-kynurenine, QA: quinolinic acid, KA: kynurenic acid, 3-HA: 3-hydroxykynurenic acid, AA: anthranilic acid.