Skip to main content
  • Research article
  • Open access
  • Published:

Risk prediction in medically treated chronic thromboembolic pulmonary hypertension



At present, there is no generally accepted comprehensive prognostic risk prediction model for medically treated chronic thromboembolic pulmonary hypertension (CTEPH) patients.


Consecutive medically treated CTEPH patients were enrolled in a national multicenter prospective registry study from August 2009 to July 2018. A multivariable Cox proportional hazards model was utilized to derive the prognostic model, and a simplified risk score was created thereafter. Model performance was evaluated in terms of discrimination and calibration, and compared to the Swedish/COMPERA risk stratification method. Internal and external validation were conducted to validate the model performance.


A total of 432 patients were enrolled. During a median follow-up time of 38.73 months (IQR: 20.79, 66.10), 94 patients (21.8%) died. The 1-, 3-, and 5-year survival estimates were 95.5%, 83.7%, and 70.9%, respectively. The final model included the following variables: the Swedish/COMPERA risk stratum (low-, intermediate- or high-risk stratum), pulmonary vascular resistance (PVR, ≤ or > 1600 dyn·s/cm5), total bilirubin (TBIL, ≤ or > 38 µmol/L) and chronic kidney disease (CKD, no or yes). Compared with the Swedish/COMPERA risk stratification method alone, both the derived model [C-index: 0.715; net reclassification improvement (NRI): 0.300; integrated discriminatory index (IDI): 0.095] and the risk score (C-index: 0.713; NRI: 0.300; IDI: 0.093) showed improved discriminatory power. The performance was validated in a validation cohort of 84 patients (C-index = 0.707 for the model and 0.721 for the risk score).


A novel risk stratification strategy can serve as a useful tool for determining prognosis and guide management for medically treated CTEPH patients.

Trial registration: (Identifier: NCT01417338).

Peer Review reports


Chronic thromboembolic pulmonary hypertension (CTEPH) is characterized by organized thromboembolic obstruction of the pulmonary arteries, which leads to progressively elevated pulmonary vascular resistance, pulmonary hypertension (PH), right heart failure and ultimately death [1]. Pulmonary endarterectomy (PEA) is recommended as the first-line treatment for CTEPH [1, 2]. As reported by the International Registry, PEA can increase the 3-year survival of incident CTEPH patients to 89%, in contrast to the 70% for non-operated patients [3]. However, despite the widely acknowledged benefits of PEA, approximately 40% of CTEPH patients are considered inoperable due to surgical inaccessibility of the thrombi, pulmonary arterial pressure disproportionate with the morphological lesions or the presence of severe comorbidities [4, 5]. For those patients, balloon pulmonary angioplasty (BPA) and riociguat, the only PH-targeted drug for CTEPH approved to date, should be considered alternative treatment options [1, 2]. Meanwhile, for non-operated patients, several studies have ventured to explore predictors of prognosis that can help clinicians identify patients at high risk, determine appropriate treatment strategies, and evaluate the efficacy of the possible treatments [3, 6,7,8,9,10,11,12,13]. As reported, numerous variables, such as World Health Organization (WHO) functional class (FC) [3, 6, 8], 6-min walk distance (6MWD) [6, 8, 13, 14], right atrial pressure (RAP) [6, 13], cardiac index [13], pulmonary vascular resistance (PVR) [12, 14] and brain natriuretic peptide (BNP)/N-terminal pro-BNP (NT-proBNP) [8], have been reported to have prognostic value for CTEPH patients.

However, at present, there is no generally accepted comprehensive prognostic risk prediction model for non-operated CTEPH patients. The 2015 European Society of Cardiology (ESC)/European Respiratory Society (ERS) guidelines proposed a risk stratification strategy for pulmonary arterial hypertension (PAH) with a range of clinical characteristics, biochemical markers and cardiac function and hemodynamics evaluations [2]. This strategy has been further abbreviated and validated in several PAH cohorts [11, 15,16,17] and two medically treated CTEPH cohorts [10, 11], with all studies demonstrating that patients in the low-risk stratum tend to have better outcomes. Meanwhile, the REVEAL risk score (RRS), another widely used risk assessment tool for PAH, has also had its utility validated for patients with inoperable and persistent/recurrent CTEPH [9]. Nevertheless, the most substantial limitation of the strategies mentioned above is that they are first derived from PAH cohorts without consideration of the prognostic factors specific to CTEPH patients [10, 11]. Furthermore, the inclusion of a broad panel of data, including relevant comorbidities and clinically available biomarkers, would also be helpful in improving the accuracy of risk prediction.

Accordingly, the objectives of the current study were to identify prognostic predictors from a broad range of data, including clinical assessments, comorbid conditions, routinely available biomarkers, evaluations of cardiac/pulmonary function and hemodynamic parameters, in a national prospective multicenter CTEPH registry dataset, and to further establish a risk assessment tool specific to medically treated CTEPH patients.


Study design

In this national multicenter prospective registry study, patients with CTEPH were consecutively recruited from 18 participating medical centers throughout China. We performed retrospective data analysis using the abovementioned prospective CTEPH registry data in the current study. The study protocol was approved by the Institutional Review Board (IRB) of Fuwai Hospital (Approval No. 2009-208), complies with the Declaration of Helsinki, and is registered on (Identifier: NCT01417338). Written informed consent was obtained from all enrolled patients.

Study participants

Patients were enrolled in the registry according to the following criteria: 1) right heart catheterization (RHC) performed within one month before enrollment between August 2009 and July 2018; 2) PH confirmed by RHC with a pulmonary arterial pressure (mPAP) ≥ 25 mmHg and pulmonary arterial wedge pressure (PAWP) ≤ 15 mmHg at rest; 3) CTEPH diagnosed based on mismatch on ventilation/perfusion (V/Q) scintigraphy with at least one large perfusion defect in one segment or two subsegments or evidence of pulmonary vascular lesions on computed tomography (CT)/magnetic resonance imaging/pulmonary angiography; 4) administration of at least 3 months of effective anticoagulation; and 5) age ranging from 14 to 85 years. Patients complicated with systemic vasculitis, severe pulmonary disease such as interstitial fibrosis or other comorbid conditions that could have caused nonthromboembolic pulmonary hypertension or those with a life expectancy of less than half a year were excluded. In the current study, we selected medically treated patients who had not undergone PEA or BPA before enrollment or during follow-up. The validation cohort consisted of medically treated CTEPH patients retrospectively enrolled from four centers across China from October 2006 to July 2009 according to the same criteria above.

Measurements and data collection

Electrocardiography (ECG), chest X-ray, transthoracic echocardiography, pulmonary function tests, V/Q scintigraphy lung scan, high-resolution CT, pulmonary angiography (if necessary), RHC and laboratory tests were performed to evaluate cardiac and pulmonary function, aid in the diagnosis and guide the treatments of CTEPH. Operability for PEA was assessed by an experienced multidisciplinary team (MDT) in the operation centers. Surgically inoperable CTEPH was defined as CTEPH in which the thrombus was in a peripheral location. For enrolled, medically treated CTEPH patients, the following data were collected: (1) demographics, medical history, clinical symptoms and vital signs; (2) examination results; and (3) treatments.

Endpoint and follow-up

The primary endpoint of this study was all-cause mortality. Overall survival was measured from the date of RHC to the date of death from any cause. Follow-up was performed by telephone calls, outpatient visits or inpatient admissions every 6 months ± 2 weeks. At each follow-up, vital status was confirmed, as well as surgical events, interventions, and instances of cardiac hospitalization. Patients were followed until death or until the cutoff date of the current study (March 2019).

Statistical analysis

Continuous variables are presented as the mean ± standard deviation or median [interquartile range (IQR)]. Differences were compared by Student’s t test or the Mann–Whitney U test for two groups and 1-way analysis of variance or the Kruskal–Wallis test for multiple groups, as appropriate. Categorical variables are shown as frequencies and percentages and were compared with the chi-square test. Multiple imputation was used to replace missing values for corresponding variables. Cox proportional hazards analyses were performed to compute hazard ratios (HRs) with 95% confidence intervals (CIs). The proportional hazards assumption was examined by the Schoenfeld residuals method. Univariable Cox analyses were first conducted to screen candidate variables, which were based on the literature and clinical expertise and included demographics, clinical assessments, comorbidities, clinically assessed biomarkers, and variables obtained from pulmonary function tests, echocardiography and RHC. In the univariable analyses, 6MWD, WHO-FC, NT-proBNP, cardiac index, RAP and mixed venous oxygen saturation (SvO2), which were included as variables in the Swedish/COMPERA risk stratification method, showed significant prognostic value for mortality [15, 16]. Because the Swedish/COMPERA risk stratification method has previously been validated in CTEPH patients [10], we integrated the six parameters into a composite variable according to the specified rules. Hence, we graded the risks for the six variables from 1 to 3 (1 = low, 2 = intermediate, 3 = high) using the recommended thresholds in the guidelines. The rounded means of these grades were then used to define the risk stratum, which was further included in the multivariable model as a priority. Because categorical specifications are favored in daily practice, other continuous risk factors were dichotomized according to the optimal thresholds determined by maximally selected rank statistics before being entered in the multivariable Cox model [18]. Thereafter, stepwise variable selection with entry and exit criteria (P < 0.05) was used to obtain the final model. For routine prognostic assessments, a simplified risk score was then derived by assigning integer numbers to each variable according to the adjusted HRs, and the overall risk score was defined as the sum of the risk points.

Model performance was evaluated for discrimination and calibration. Harrell’s C-index was used to assess the discriminatory power, while the calibration of the 5-year risk prediction was visually evaluated by plotting and comparing the predicted and observed risk. For model comparison, net reclassification improvement (NRI) and integrated discriminatory index (IDI) were used. For NRI, we classified the 5-year mortality into low (< 25%), intermediate (25–50%) and high (> 50%) based on previous literature [11]. Survival was estimated by means of Kaplan–Meier analysis, and difference were compared by the log-rank test. Sensitivity analyses were performed in the following subgroups: (1) newly diagnosed patients, (2) surgically inoperable patients, and (3) a subset that excluded patients with chronic liver disease. Internal validation was evaluated with bootstrapping. External validation was performed to justify whether the derived model and the risk score were also predictive of death in a validation cohort. Differences were considered statistically significant when the two-sided P value was < 0.05. All analyses were performed with the R statistical package (version 4.0.0, R Foundation for Statistical Computing, Vienna, Austria).


Study cohort

A total of 432 medically treated CTEPH patients were enrolled in the current study (Fig. 1). The mean age of the cohort was 53.54 ± 12.25 years, and 52.3% of the patients were males. A total of 59.3% of these patients were newly diagnosed, and 54.2% were surgically inoperable. A total of 251 (58.1%) patients had a history of pulmonary embolism, and chronic kidney disease, defined as an estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73 m2, was the most common comorbidity (21.3%). In addition to conventional medical therapies, 55.3% of the patients received at least one PH-targeted drug, while only 9.3% of them received combination therapy. However, as riociguat, the only approved targeted drug for CTEPH, was not available in China before June 2018, it was not in used by any of the patients in the current study cohort. Baseline patient characteristics are presented in Table 1.

Fig. 1
figure 1

Flowchart of patient selection. CTEPH: chronic thromboembolic pulmonary hypertension; PEA: pulmonary endarterectomy; BPA: balloon pulmonary angiography

Table 1 Baseline characteristics of the overall analyzed cohort, survivors and non-survivors

During a median follow-up time of 38.73 months (IQR: 20.79, 66.10), 94 patients (21.8%) died. The leading cause of death was right heart failure, and other frequent causes included infection, hemoptysis, and respiratory failure; sudden death was also commonly reported in the cohort. The survival estimates at 1, 3 and 5 years were 95.5% (95% CI: 93.5–97.5%), 83.7% (79.9–87.8%), and 70.9% (65.5–76.7%), respectively.

Prognostic variables

In univariable Cox analysis, WHO-FC, 6MWD, SvO2, RAP, cardiac index, mPAP, PVR, NT-proBNP, total bilirubin (TBIL), creatinine, uric acid, blood urea nitrogen, coronary heart disease and CKD were significant predictors for survival (Additional file 1: Table S1). In addition, the Swedish/COMPERA risk stratum also showed significant predictive value for survival: compared with the low-risk group, the intermediate- (HR 4.008; 95% CI 1.746–9.199) and high-risk (HR 9.740; 95% CI 3.537–26.820) groups demonstrated increased risks for mortality (Additional file 1: Table S1 and Figure S1).

Risk model and simplified risk score

According to multivariable Cox proportional hazards analysis, the final model included the following variables: the Swedish/COMPERA risk stratum (low-, intermediate- or high-risk stratum), PVR (≤ or > 1600 dyn s/cm5), TBIL (≤ or > 38 µmol/L) and CKD (no or yes). The adjusted HRs and CIs are presented in Table 2. The model showed a C-index of 0.715 (95% CI 0.660–0.770) with good calibration (Fig. 2a).

Table 2 Multivariate Cox proportional hazards model for predicting mortality in medically treated CTEPH patients
Fig. 2
figure 2

Calibration of the derived model (a) and the risk score (b). Five-year rates of mortality as predicted versus the observed rates. The diagonal line indicates perfect calibration

To simplify its use in clinical practice, a risk score was then derived based on the adjusted coefficients of the above model (Table 2). The overall risk score, which ranged from 0 to 13, demonstrated good discriminatory power similar to that of the model, with a C-index of 0.713 (95% CI 0.658–0.768). Figure 2b shows that the risk score was also well calibrated.

Compared with the model consisting only of the Swedish/COMPERA risk stratum, both the full model [NRI 0.300 (0.100–0.542); IDI 0.095 (0.041–0.191)] and the risk score [NRI 0.300 (0.091–0.480); IDI 0.093 (0.044–0.169)] demonstrated improved discriminatory performance (Table 3). Based on the distribution of the risk score, we divided the patients into 3 risk groups, with low-, intermediate- and high-risk groups determined by risk scores of 0–3, 4–5 or ≥ 6 points, respectively.

Table 3 C-index, NRI and IDI of the derived model and risk score

The characteristics of the three subgroups are shown in Additional file 1: Table S2. The variables included in the derived risk model differed significantly between the three risk score groups; furthermore, the survival differences between the three risk categories were statistically significant, with P < 0.0001 for all comparisons, P = 0.006 for the 0–3 versus 4–5 points group comparison, P = 0.007 for the 4–5 versus ≥ 6 points group comparison and P < 0.001 for the 0–3 versus ≥ 6 points group comparison (Fig. 3). The estimated 5-year survival rates for the three groups were 91.9% (84.9–99.5%), 79.0% (71.9–868%) and 49.3% (40.0–60.7%), respectively (Table 4).

Fig. 3
figure 3

Comparison of Kaplan–Meier survival curves according to the derived risk score (0–3, 4–5 or ≥ 6 points)

Table 4 Estimated 1-, 3- and 5-year survival of the three risk groups

Sensitivity analysis and external validation

Sensitivity analyses were performed in three subgroups of patients: patients who were newly diagnosed, those who were surgically inoperable and those without chronic liver disease (excluding 5 patients with chronic hepatitis). The baseline characteristics are shown in Additional file 1: Tables S3-S5. Both the derived model [C-index in newly diagnosed patients: 0.709 (95% CI 0.638–0.780), in surgically inoperable patients: 0.737 (0.649–0.825), in patients without chronic liver disease: 0.718 (0.663–0.773)] and the risk score [C-index: 0.709 (0.640–0.778), 0.743 (0.655–0.831), 0.716 (0.661–0.771), respectively] showed consistent significant predictive value for survival in subgroup analyses. The three risk groups defined by the new risk score also showed significant discriminatory power (Additional file 1: Tables S6-S8 and Figures S2-S4).

The validation cohort consisted of 84 medically treated CTEPH patients with a mean age of 52.01 ± 12.94 years. Twenty-nine patients (34.5%) died during a median follow-up time of 103.56 months (IQR 48.70, 120.00). The survival estimates at 1, 3 and 5 years were 96.4%, 84.2%, and 75.7, respectively. The baseline data of this cohort are presented in Additional file 1: Table S9. Both the derived model [C-index: 0.707 (95% CI: 0.623–0.791)] and the risk score [C-index: 0.721 (0.633–0.809)] demonstrated consistent performance in the validation cohort (Additional file 1: Figure S5). Due to the limited sample size, no subgroup analyses were further performed in the validation cohort.


The present study identified predictors for survival from a broad range of data collected for the study group of medically treated CTEPH patients and thereafter established a new risk assessment tool specific to CTEPH patients. To the best of our knowledge, this is the first risk prediction model specifically derived from the data from CTEPH patients.

Previous studies have reported that numerous variables are associated with the outcomes of medically treated CTEPH patients [3, 6,7,8,9,10,11,12,13]. However, regarding risk prediction tools, published and recommended risk stratification strategies have only been derived for PAH, including the RRS [19, 20], the recommended strategy from the ESC/ERS guidelines [2] and its three abbreviated versions [15,16,17]. Although these strategies have been validated in CTEPH patients [9,10,11], it is unknown whether other variables relevant to CTEPH patients would add incremental value. The current study derived a new risk prediction strategy based on real-life registry data from medically treated CTEPH patients and demonstrated that the newly derived prediction model, which combined the Swedish/COMPERA risk stratification method with PVR, serum TBIL and CKD, performed well in predicting survival in those patients. It should be noted that despite the validation of the Swedish/COMPERA risk stratification method, the six variables included in the stratum have also been reported to be individually associated with the outcomes of non-operated CTEPH patients [6, 8, 13, 14]. In addition, PVR is also widely acknowledged as a prognostic factor for CTEPH patients [12, 14].

In contrast to the above-listed variables, TBIL and CKD have been infrequently reported as risk factors in CTEPH patients. In a small sample of 77 inoperable CTEPH patients, the serum concentration of TBIL was found to be an independent prognostic predictor for mortality, with patients whose TBIL ≥ 23.7 µmol/L having markedly worse survival [21]. Similarly, hyperbilirubinemia (serum TBIL > 1.2 mg/dL) was also reported as a predictor of mortality in PAH patients [22]. Furthermore, among biomarkers related to hepatic function, elevated total bilirubin could be the strongest predictor for the adverse outcome of cardiovascular death, superior to transaminases in sensitivity to hemodynamic abnormalities [23]. As chronic liver disease can also affect variables concerning liver function, such as total bilirubin, we further performed sensitivity analyses based on a subgroup consisting of patients without chronic liver disease, where 5 patients with chronic hepatitis were excluded from the cohort. Notably, after excluding these patients, the derived risk model and the risk score showed consistent significant prognostic power. Similar to hepatic dysfunction, right heart failure may also be the potential link between CTEPH and CKD. Renal insufficiency is included in the original RRS [19, 20]; the latest version, RRS 2.0, has updated this category as eGFR < 60 mL/min/1.73 m2 or renal insufficiency, as renal function is an important risk predictor for PAH patients [24, 25]. Our study further supports the prognostic use of renal insufficiency in CTEPH patients. Regarding utility, the risk prediction tools derived in our study also emphasize the importance of controlling these comorbid conditions, which could have significant effects on the outcomes.

The 1-, 3- and 5-year survival estimates for medically treated CTEPH patients in the current study cohort (95.5%, 83.7%, 70.9%) were higher than those of the International Registry (1- and 3-year survival estimates: 88% and 70%, respectively)3 or those reported by Delcroix et al. (1-, 3-, and 5-year survival estimates: 92.0%, 74.7%, and 59.8%, respectively) [10]. Other differences with this latter study, such as the younger cohort age (53 vs. 69 years), the much lower percentages of comorbidities (any comorbidities: 70% versus 91%), and of patients with intermediate risk (18% versus 68%) and the smaller number of patients who received PH-targeted therapy–especially the unavailability of stimulators of soluble guanylate cyclase (sGCs) in our cohort (versus 37% use in Delcroix et al.)–should be elucidated [10]. The potential reason for these differences may be the disparity in the enrolled patients, as the current registry also enrolled previously diagnosed patients, which could lead to potential survival bias. Therefore, the survival estimates are more similar to those of Spanish Registry of Pulmonary Arterial Hypertension (REHAP) registry (1-, 3-, and 5-year survival estimates: 92.6%, 80.7%, and 64.9%, respectively), which included both newly and previously diagnosed patients [26, 27]. As the cohort in the current study included both newly and previously diagnosed patients, similar to daily practice, our results can be more broadly generalized to any clinical scenario. Furthermore, it should also be noted that we performed sensitivity analysis in the subgroup of newly diagnosed patients, and the model illustrated consistently good performance in risk prediction.

Despite the fact that PEA is recommended as the first-line treatment for CTEPH, only 14% of the patients in our study underwent this procedure, while approximately 46% who were surgically operable did not. The low rate of PEA could be largely attributed to the unbalanced development of medical centers between different areas, as there are only three surgical centers in China, all located in Beijing. Furthermore, the high financial cost of PEA could be another barrier for patients to undergo the procedure. Regarding its effect on the study, as shown in the sensitivity analysis, the model performed consistently well for surgically inoperable patients. Additionally, it should be noted that although the surgery rate was much lower than that in Western countries, it reflects real-world experience regarding the treatment of CTEPH in developing countries to some extent.

We regarded the six variables included in the Swedish/COMPERA method as a composite variable without redefining the existing categories. Due to the relatively small number of events, the significant predictive value of the six parameters as continuous variables in the univariable analyses, and the fact that the risk strata have been previously validated and have shown consistent discriminative power in our cohort [10], we did not include each variable separately with new cutoffs, which may have helped avoid overfitting. However, it is possible that the stratification strategy with the cutoff values first derived for PAH might not be suitable for patients with CTEPH. Therefore, further studies are still needed to investigate appropriate cutoff values specific to CTEPH.

Our study has several limitations. First, we only evaluated variables at baseline to establish and validate the model and did not utilize the data from follow-up visits, as these data were incomplete; only 34 patients (7.9%) had data on follow-up RHC. Among these patients, 2 met the primary endpoint during the follow-up period. Therefore, both the small sample size and the low event rate in the follow-up data prevented us from performing further analyses. However, as reevaluations at follow-ups are necessary for risk assessment and treatment guidance, further studies are required to include information obtained at follow-up to achieve better assessments. Second, as mentioned above, we did not include the six variables in the Swedish/COMPERA risk strata separately and did not find new cutoff values, which could lead to potential inconsistency with the actual categories for the CTEPH patients. Finally, as we only performed external validation in a retrospective cohort with a small sample size, further external validation in independent larger cohorts is required.


Our novel risk stratification strategy can serve as a useful tool for determining prognosis and guiding management for medically treated CTEPH patients.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. Kim NH, Delcroix M, Jais X, Madani MM, Matsubara H, Mayer E, et al. Chronic thromboembolic pulmonary hypertension. Eur Respir J. 2019;53(1):D92–9.

    Google Scholar 

  2. Galie N, Humbert M, Vachiery JL, Gibbs S, Lang I, Torbicki A, et al. 2015 ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT). Eur Heart J. 2016;37(1):67–119.

    Article  Google Scholar 

  3. Delcroix M, Lang I, Pepke-Zaba J, Jansa P, D’Armini AM, Snijder R, et al. Long-term outcome of patients with chronic thromboembolic pulmonary hypertension: results from an international prospective registry. Circulation. 2016;133(9):859–71.

    Article  CAS  Google Scholar 

  4. Pepke-Zaba J, Delcroix M, Lang I, Mayer E, Jansa P, Ambroz D, et al. Chronic thromboembolic pulmonary hypertension (CTEPH): results from an international prospective registry. Circulation. 2011;124(18):1973–81.

    Article  Google Scholar 

  5. Mayer E, Jenkins D, Lindner J, D’Armini A, Kloek J, Meyns B, et al. Surgical management and outcome of patients with chronic thromboembolic pulmonary hypertension: results from an international prospective registry. J Thorac Cardiovasc Surg. 2011;141(3):702–10.

    Article  Google Scholar 

  6. Taniguchi Y, Jaïs X, Jevnikar M, Boucly A, Weatherald J, Brenot P, et al. Predictors of survival in patients with not-operated chronic thromboembolic pulmonary hypertension. J Heart Lung Transplant. 2019;38(8):833–42.

    Article  Google Scholar 

  7. Lewczuk J, Piszko P, Jagas J, Porada A, Wójciak S, Sobkowicz B, et al. Prognostic factors in medically treated patients with chronic pulmonary embolism. Chest. 2001;119(3):818–23.

    Article  CAS  Google Scholar 

  8. Simonneau G, D’Armini AM, Ghofrani HA, Grimminger F, Jansa P, Kim NH, et al. Predictors of long-term outcomes in patients treated with riociguat for chronic thromboembolic pulmonary hypertension: data from the CHEST-2 open-label, randomised, long-term extension trial. Lancet Respir Med. 2016;4(5):372–80.

    Article  CAS  Google Scholar 

  9. Benza RL, Farber HW, Frost A, Grunig E, Hoeper MM, Busse D, et al. REVEAL risk score in patients with chronic thromboembolic pulmonary hypertension receiving riociguat. J Heart Lung Transplant. 2018;37(7):836–43.

    Article  Google Scholar 

  10. Delcroix M, Staehler G, Gall H, Grünig E, Held M, Halank M, et al. Risk assessment in medically treated chronic thromboembolic pulmonary hypertension patients. Eur Respir J. 2018;52(5):1800248.

    Article  Google Scholar 

  11. Humbert M, Farber HW, Ghofrani HA, Benza RL, Busse D, Meier C, et al. Risk assessment in pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension. Eur Respir J. 2019;53(6):1802004.

    Article  Google Scholar 

  12. Scholzel BE, Post MC, Thijs Plokker HW, Snijder RJ. Clinical worsening during long-term follow-up in inoperable chronic thromboembolic pulmonary hypertension. Lung. 2012;190(2):161–7.

    Article  Google Scholar 

  13. Condliffe R, Kiely DG, Gibbs JS, Corris PA, Peacock AJ, Jenkins DP, et al. Prognostic and aetiological factors in chronic thromboembolic pulmonary hypertension. Eur Respir J. 2009;33(2):332–8.

    Article  CAS  Google Scholar 

  14. Saouti N, de Man F, Westerhof N, Boonstra A, Twisk J, Postmus PE, et al. Predictors of mortality in inoperable chronic thromboembolic pulmonary hypertension. Respir Med. 2009;103(7):1013–9.

    Article  Google Scholar 

  15. Kylhammar D, Kjellstrom B, Hjalmarsson C, Jansson K, Nisell M, Soderberg S, et al. A comprehensive risk stratification at early follow-up determines prognosis in pulmonary arterial hypertension. Eur Heart J. 2018;39(47):4175–81.

    Article  Google Scholar 

  16. Hoeper MM, Kramer T, Pan Z, Eichstaedt CA, Spiesshoefer J, Benjamin N, et al. Mortality in pulmonary arterial hypertension: prediction by the 2015 European pulmonary hypertension guidelines risk stratification model. Eur Respir J. 2017;50(2):1700740.

    Article  Google Scholar 

  17. Boucly A, Weatherald J, Savale L, Jais X, Cottin V, Prevot G, et al. Risk assessment, prognosis and guideline implementation in pulmonary arterial hypertension. Eur Respir J. 2017;50(2):1700889.

    Article  Google Scholar 

  18. Kassambara A KM, Biecek P, Fabian S. Package “survminer.”.

  19. Benza RL, Miller DP, Gomberg-Maitland M, Frantz RP, Foreman AJ, Coffey CS, et al. Predicting survival in pulmonary arterial hypertension: insights from the Registry to Evaluate Early and Long-Term Pulmonary Arterial Hypertension Disease Management (REVEAL). Circulation. 2010;122(2):164–72.

    Article  Google Scholar 

  20. Benza RL, Gomberg-Maitland M, Miller DP, Frost A, Frantz RP, Foreman AJ, et al. The REVEAL Registry risk score calculator in patients newly diagnosed with pulmonary arterial hypertension. Chest. 2012;141(2):354–62.

    Article  Google Scholar 

  21. Gong JN, Zhai ZG, Yang YH, Liu Y, Gu S, Kuang TG, et al. Serum bilirubin and 6-min walk distance as prognostic predictors for inoperable chronic thromboembolic pulmonary hypertension: a prospective cohort study. Chin Med J. 2015;128(23):3125–31.

    Article  CAS  Google Scholar 

  22. Takeda Y, Takeda Y, Tomimoto S, Tani T, Narita H, Kimura G. Bilirubin as a prognostic marker in patients with pulmonary arterial hypertension. BMC Pulm Med. 2010;10:22.

    Article  Google Scholar 

  23. Allen LA, Felker GM, Pocock S, McMurray JJ, Pfeffer MA, Swedberg K, et al. Liver function abnormalities and outcome in patients with chronic heart failure: data from the Candesartan in Heart Failure: Assessment of Reduction in Mortality and Morbidity (CHARM) program. Eur J Heart Fail. 2009;11(2):170–7.

    Article  Google Scholar 

  24. Chakinala MM, Coyne DW, Benza RL, Frost AE, McGoon MD, Hartline BK, et al. Impact of declining renal function on outcomes in pulmonary arterial hypertension: A REVEAL registry analysis. J Heart Lung Transplant. 2018;37(6):696–705.

    Article  Google Scholar 

  25. Benza RL, Gomberg-Maitland M, Elliott CG, Farber HW, Foreman AJ, Frost AE, et al. Predicting survival in patients with pulmonary arterial hypertension: the REVEAL risk score calculator 2.0 and comparison with ESC/ERS-based risk assessment strategies. Chest. 2019;156(2):323–37.

    Article  Google Scholar 

  26. Escribano-Subias P, Blanco I, López-Meseguer M, Lopez-Guarch CJ, Roman A, Morales P, et al. Survival in pulmonary hypertension in Spain: insights from the Spanish registry. Eur Respir J. 2012;40(3):596.

    Article  Google Scholar 

  27. Escribano-Subías P, Del Pozo R, Román-Broto A, Domingo Morera JA, Lara-Padrón A, Elías Hernández T, et al. Management and outcomes in chronic thromboembolic pulmonary hypertension: From expert centers to a nationwide perspective. Int J Cardiol. 2016;203:938–44.

    Article  Google Scholar 

Download references




This study was supported by National Key Technology R&D Program of China (Project Number: 2011BAI11B15) and National Key Research and Development Program of China (Project Number: 2016YFC1304400). The funding body had no role in the design of the study, and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



Contributing to the conception and design: RLQ, CMX and JGH. Patient enrollment and data collection: RLQ, YHY, ZWY, HYT, SQL, JYS, YQJ, GCZ, CJZ, GYW, YHL, ZZC, ZXY, ZYS, ZQZ, WC, YCC, SL, XXC, YLQ. Data analysis, interpretation: RLQ, CMX, GLS. Drafting the article: RLQ. Revising the article: RLQ and JGH. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jianguo He.

Ethics declarations

Ethics approval and consent to participate

This study complied with the Declaration of Helsinki. The study protocol was approved by the Institutional Review Board (IRB) of Fuwai Hospital (Approval No. 2009-208). Written informed consent was obtained from each enrolled patient.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Figure S1. Comparison of Kaplan-Meier survival curves according to the Swedish/COMPERA risk stratum (Low-, intermediate- and high-risk groups). Figure S2. Comparison of Kaplan-Meier survival curves according to the derived risk score (0-3, 4-5 or ≥ 6 points) in newly diagnosed CTEPH patients. Figure S3. Comparison of Kaplan-Meier survival curves according to the derived risk score (0-3, 4-5 or ≥ 6 points) in surgically inoperable CTEPH patients. Figure S4. Comparison of Kaplan-Meier survival curves according to the derived risk score (0-3, 4-5 or ≥ 6 points) in patients without chronic liver disease. Figure S5. Calibration of the derived model (A) and the risk score (B) in the validation cohort. Table S1. Univariate Cox proportional hazards analyses of candidate variables for all-cause mortality in the overall analyzed cohort. Table S2. Baseline characteristics of the 3 risk groups classified by the new derived risk score. Table S3. Baseline characteristics of the newly diagnosed CTEPH patients. Table S4. Baseline characteristics of the surgically inoperable CTEPH patients. Table S5. Baseline characteristics of patients without chronic liver disease. Table S6. Estimated 1-, 3- and 5-year survival of the three risk groups in newly diagnosed CTEPH patients. Table S7. Estimated 1-, 3- and 5-year survival of the three risk groups in surgically inoperable patients. Table S8. Estimated 1-, 3- and 5-year survival of the three risk groups in patients without chronic liver disease. Table S9. Baseline characteristics of the validation cohort.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Quan, R., Yang, Y., Yang, Z. et al. Risk prediction in medically treated chronic thromboembolic pulmonary hypertension. BMC Pulm Med 21, 128 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: