Abstract
The prognostic importance of follow-up haemodynamics and the validity of multidimensional risk assessment are not well established for systemic sclerosis (SSc)-associated pulmonary arterial hypertension (PAH).
We assessed incident SSc-PAH patients to determine the association between clinical and haemodynamic variables at baseline and first follow-up right heart catheterisation (RHC) with transplant-free survival. RHC variables included cardiac index, stroke volume index (SVI), pulmonary arterial compliance and pulmonary vascular resistance. Risk assessment was performed according to the number of low-risk criteria: functional class I or II, 6-min walking distance (6MWD) >440 m, right atrial pressure <8 mmHg and cardiac index ≥2.5 L·min−1·m−2.
Transplant-free survival from diagnosis (n=513) was 87%, 55% and 35% at 1, 3 and 5 years, respectively. At baseline, 6MWD was the only independent predictor. A follow-up RHC was available for 353 patients (median interval 4.6 months, interquartile range 3.9–6.4 months). The 6MWD, functional class, cardiac index, SVI, pulmonary arterial compliance and pulmonary vascular resistance were independently associated with transplant-free survival at follow-up, with SVI performing better than other haemodynamic variables. 1-year outcomes were better with increasing number of low-risk criteria at baseline (area under the curve (AUC) 0.63, 95% CI 0.56–0.69) and at first follow-up (AUC 0.71, 95% CI 0.64–0.78).
Follow-up haemodynamics and multidimensional risk assessment had greater prognostic significance than at baseline in SSc-PAH.
Abstract
For newly diagnosed patients with systemic sclerosis-associated PAH, some haemodynamic variables, particularly the stroke volume index, and a multidimensional risk assessment were more useful during early follow-up than at baseline http://ow.ly/LUuj30lrvfW
Introduction
Pulmonary arterial hypertension (PAH) develops in ∼10% of individuals with systemic sclerosis (SSc) [1] and, when present, is one of the most frequent causes of death in these patients [2]. The diagnosis of PAH requires haemodynamic measurements during right heart catheterisation (RHC), with a mean pulmonary artery pressure (mPAP) ≥25 mmHg and pulmonary artery wedge pressure ≤15 mmHg [3].
While essential for the diagnosis of PAH, haemodynamic variables such as the right atrial pressure (RAP), cardiac index and stroke volume provide important information about right ventricular (RV) function, and predict survival in SSc-PAH [4–6]. Several observational studies, predominantly of idiopathic PAH (IPAH) patients, have noted the importance of clinical variables, such as the 6-min walking distance (6MWD), New York Heart Association functional class (NYHA-FC) and haemodynamic variables achieved during early follow-up as they reflect response to initial PAH therapy [7–10]. The long-term prognostic significance of early follow-up clinical and haemodynamic variables after initial treatment has not been specifically studied in the SSc-PAH population.
Periodic risk assessment using a multidimensional approach is recommended by recent European Society of Cardiology (ESC) and European Respiratory Society (ERS) guidelines, with thresholds for a “low-risk” patient defined by several clinical, exercise, imaging-derived and haemodynamic parameters [3]. A study from the French PAH registry demonstrated that the number of low-risk criteria present was a simple and valid method of assessing risk in IPAH and heritable and drug-induced PAH [11]. However, the validity of this approach to risk assessment has not been studied specifically in SSc-PAH.
The main objective of this study was to evaluate the prognostic value of haemodynamic variables at baseline and after initial treatment in newly diagnosed SSc-PAH patients. A secondary objective was to determine whether the number of low-risk criteria at baseline and first follow-up RHC was associated with transplant-free survival in SSc-PAH.
Methods
This study complied with the Declaration of Helsinki. Although French law does not require ethics committee approval or informed consent for retrospective data collection, the data collected in the French Pulmonary Arterial Hypertension Network registry were anonymised and complied with the requirements of the Commission Nationale Informatique et Liberté (CNIL). CNIL, the organisation dedicated to privacy, information technology and civil rights in France, approved the methods used to collect and analyse data on May 24, 2003 (approval number 842063).
Study population
Incident, treatment-naïve patients in the French Pulmonary Arterial Hypertension Network registry with SSc-PAH were prospectively enrolled between January 2006 and March 2017. Patients were eligible for inclusion if they had a diagnosis of SSc, were ≥18 years of age, and had newly diagnosed group 1 precapillary PAH on RHC, defined as resting mPAP ≥25 mmHg, pulmonary artery wedge pressure ≤15 mmHg and pulmonary vascular resistance (PVR) >3 Wood units [3]. The diagnosis of SSc was made according to current criteria by the treating physician according to American College of Rheumatology/European League Against Rheumatism criteria [12]. Patients with a recorded history of severe interstitial lung disease on computed tomography in the registry, those coded as having group 2 (secondary to left heart disease), group 3 (due to lung disease or chronic hypoxaemia), group 4 (chronic thromboembolic pulmonary hypertension) or with missing baseline haemodynamic variables were excluded [3]. Classification of patients as having group 1 PAH was at the discretion of clinicians at each site.
Measurements
Variables were assessed at baseline and at the time of the first follow-up RHC within 12 months of diagnosis. We recorded standard haemodynamic variables such as RAP, mPAP, cardiac output, cardiac index, PVR and mixed venous oxygen saturation (SvO2), and calculated variables: stroke volume index (SVI) from the cardiac index divided by heart rate, and pulmonary arterial compliance (PCa), calculated by stroke volume divided by pulse pressure (the difference between systolic and diastolic pulmonary arterial pressure). The initial treatment strategy was defined according to number of PAH medications prescribed within 4 months of the initial RHC. Low-risk criteria were defined as NYHA-FC I or II, 6MWD >440 m, RAP <8 mmHg and cardiac index ≥2.5 L·min−1·m−2 [3]. The primary outcome was death or lung transplantation. Transplant-free survival time was defined from the date of initial RHC until the occurrence of a primary outcome event for the baseline analysis and from the date of first follow-up RHC until the occurrence of a primary outcome event for the follow-up analysis. Patients without an outcome event were censored at the time of last clinical contact.
Statistical analysis
Continuous variables are expressed as mean±sd for normally distributed variables or median (interquartile range (IQR)) for non-normally distributed variables. Normality was assessed using the Kolmogorov–Smirnov test. Categorical variables are presented as n (%). Changes from baseline to follow-up were assessed using a paired t-test, Wilcoxon signed-rank test or Chi-squared test where appropriate. All comparisons were two-sided, with a p-value <0.05 considered significant. Cox proportional-hazards regression and Kaplan–Meier analysis were used to assess the association between variables at baseline, at first follow-up RHC and the number of low-risk criteria with transplant-free survival. The proportional-hazards assumption was tested using log-minus-log plots and collinearity of haemodynamic variables was assessed using linear regression or Spearman's rank correlation. Variables with a p-value ≤0.1 in the univariable analysis were eligible for entry into the multivariable models only if they were not highly correlated (absolute value of Pearson's r or Spearman's ρ <0.6) with other continuous variables and if <25% of individuals had missing values for that variable. Case-wise deletion was used for multivariable modelling with no imputation for missing values. Statistical analyses were performed using STATA version 13.1 for Mac (College Station, TX, USA).
Role of the funding source
The funding bodies supporting this work had no role in the study design, data collection, analysis, interpretation of the data or writing of the manuscript.
Results
Study population and treatment strategies
Of 573 potentially eligible patients in the registry between 2006 and 2017, 513 were included in the baseline analysis (figure 1). Patient characteristics are shown in table 1. The majority were female (78.2%) and had limited systemic sclerosis (69.2%). Most patients were in NYHA-FC III (59.6%) and 13.4% were in NYHA-FC IV. The most frequent initial treatment strategy was monotherapy (56.1%) followed by initial combination therapy with two PAH therapies (31.4%). Six (1.2%) patients received initial triple combination therapy. Death occurred in 256 (49.9%) patients during follow-up and 12 (2.3%) received lung transplantation. Total follow-up from diagnosis was 1331.88 person-years with a median (IQR) follow-up time of 2.09 (0.96–3.69) years. Median transplant-free survival from the time of diagnosis was 3.42 years (41 months) in the total population, and 1-, 3- and 5-year rates were 87%, 55% and 35%, respectively (online supplementary figure S1).
Study inclusion diagram. SSc-PAH: systemic sclerosis-associated pulmonary arterial hypertension; ILD: interstitial lung disease; mPAP: mean pulmonary arterial pressure; PVR: pulmonary vascular resistance; RHC: right heart catheterisation.
Baseline characteristics
Changes in clinical and haemodynamic variables with initial treatment
There were 353 patients with a follow-up RHC performed within a year of diagnosis (median interval 4.6 months, IQR 3.9–6.4 months). Changes in functional and haemodynamic variables from baseline to first follow-up are shown in table 2. There were significant improvements in the proportion of patients in NYHA-FC I or II (24% versus 43.5%, p<0.001), median 6MWD (300 m versus 318 m, p<0.001) and most haemodynamic variables.
Changes in exercise, functional capacity and haemodynamic variables at follow-up
Baseline variables and transplant-free survival
The univariable associations between baseline characteristics and transplant-free survival are shown in table 3. In multivariable models that included age, BMI, 6MWD, NYHA-FC, RAP, mPAP and either cardiac index, SVI, PVR or PCa, the 6MWD was the only variable independently associated with transplant-free survival (adjusted hazard ratio (HR) 0.97 per 10 m, 95% CI 0.95–0.99 per 10 m; p<0.001). Although several baseline haemodynamic variables were significant in the univariable analysis, none were independently associated, regardless of which variable (cardiac index, SVI, PCa or PVR) was entered, in the multivariable model.
Univariable analyses for baseline and follow-up factors and transplant-free survival
Risk assessment at baseline
At baseline, the proportions of patients with 0, 1, 2, 3 or 4 low-risk criteria were 20%, 35%, 30%, 11% and 3%, respectively (figure 2). There was a significantly lower risk of death or transplantation with increasing number of low-risk criteria (online supplementary table S1). Figure 3a shows transplant-free survival according to the number of low-risk criteria present at the time of PAH diagnosis. The ability of low-risk criteria at baseline to discriminate 1-year survival was modest with an area under the curve (AUC) of 0.63 (95% CI 0.56–0.69) (online supplementary table S2).
Number of low-risk criteria present at baseline and first follow-up right heart catheterisation. NYHA-FC: New York Heart Association functional class; 6MWD: 6-min walking distance; RAP: right atrial pressure; CI: cardiac index.
Transplant-free survival according to the number of low-risk criteria present at a) baseline and b) first follow-up right heart catheterisation.
Follow-up variables and transplant-free survival
After the first follow-up RHC, median (IQR) additional follow-up was 1.71 (0.70–3.15) years and total observation time was 802.35 person-years. At the time of first follow-up RHC, the 6MWD, NYHA-FC and most haemodynamic variables were associated with transplant-free survival in univariable analysis (table 3). As cardiac index, SVI, PCa and PVR were highly correlated at follow-up (Spearman ρ ranged from 0.62 to 0.87), separate multivariable models incorporating each haemodynamic variable with the other variables with p≤0.1 in table 4. In these multivariable models, the cardiac index, SVI, PVR and PCa at follow-up were independently associated with the primary outcome after adjusting for age, sex, 6MWD, NYHA-FC, RAP and mPAP. Because of a high degree of collinearity, mPAP was not included in the models with PVR (r=0.73, p<0.001) or PCa (r=−0.67, p<0.001). The SVI at first follow-up RHC had the lowest −log likelihood and Akaike's information criterion value. Figure 4 shows transplant-free survival according to quartiles for cardiac index, SVI, PVR and PCa at follow-up RHC (log-rank test p<0.001 for all). Among patients who achieved or maintained cardiac index ≥2.5 L·min−1·m−2 at the follow-up RHC, 28% had a SVI <38 mL·m−2 [10] and these patients had significantly worse survival than those with an SVI ≥38 mL·m−2 (online supplementary figure S2).
Multivariable analyses for clinical and haemodynamic variables at first follow-up right heart catheterisation
Transplant-free survival according to follow-up haemodynamic variable quartiles for a) cardiac index, b) stroke volume index,c) pulmonary arterial compliance and d) pulmonary vascular resistance.
Risk assessment at follow-up
At first follow-up RHC (n=353), the proportions of patients with 0, 1, 2, 3 or 4 low-risk criteria were 15%, 26%, 30%, 22% and 6%, respectively (figure 2). Transplant-free survival was better with an increasing number of low-risk criteria present at first follow-up (figure 3b, online supplementary table S3) All low-risk criteria at follow-up were independently associated with outcomes (online supplementary table S4). Discriminating ability for 1-year outcome was reasonable (AUC 0.71, 95% CI 0.64–0.78), such that having at least three low-risk criteria had 31.4% sensitivity and 92.7% specificity for surviving without lung transplantation in the subsequent year (online supplementary table S5). In addition, we performed an exploratory analysis of the 189 (54%) patients with brain natriuretic peptide (BNP) or N-terminal pro-BNP (NT-proBNP) at follow-up to determine whether achieving the low-risk thresholds for these biomarkers (NT-proBNP <300 ng·L−1 and BNP <50 ng·L−1) was associated with transplant-free survival. Although having NT-proBNP/BNP low-risk criteria at follow-up was associated with better outcomes in univariable analysis (HR 0.31, 95% CI 0.12–0.83; p=0.02), it was not independently associated when added to the other risk assessment criteria (online supplementary table S6). These results were similar whether or not the invasive haemodynamic criteria were included in the model (data not shown).
Discussion
The main findings from this large, multicentre cohort study of incident SSc-PAH patients were: 1) baseline haemodynamic variables were not independent predictors of long-term transplant-free survival; 2) baseline risk assessment using the number of low-risk criteria, as suggested in the ESC/ERS guidelines, allowed modest discrimination of future risk; 3) follow-up haemodynamic variables within the first year of diagnosis were significantly and independently associated with transplant-free survival, with the SVI providing better explanation of the data than cardiac index, PVR or PCa; and 4) risk assessment using the number of low-risk criteria present at first follow-up RHC discriminated transplant-free survival better than at baseline. To our knowledge, this is the largest study to assess haemodynamics as prognostic factors in SSc-PAH and the first to compare the prognostic importance of haemodynamic variables achieved after initial treatment in this specific subgroup.
We confirmed our primary hypothesis that the SVI provides more useful clinical information than the cardiac index in patients with PAH, which was based upon our recent findings in patients with IPAH and drug-induced and heritable PAH [10]. In our SSc-PAH cohort, SVI at follow-up predicted transplant-free survival when adjusted for demographic, functional and RV loading variables (i.e. RAP and mPAP; table 4). Unadjusted transplant-free survival for patients in the highest two quartiles of SVI (>37 mL·m−2) at first follow-up were similar (figure 4), which is strikingly similar to the optimal cut-point of 38 mL·m−2 in patients with IPAH, heritable and drug-induced PAH [10], and which corresponds to the lower limit of normal for RV function [13]. Furthermore, we found that an important proportion of SSc-PAH patients who achieved a cardiac index goal of ≥2.5 L·min−1·m−2 still had a low SVI <38 mL·m−2 (28%) and a worse prognosis (online supplementary figure S 2), confirming that cardiac index can be falsely reassuring when the SVI is not also taken into account, and that SVI is likely a better haemodynamic indicator of RV function. While PCa reflects the pulsatile component of RV load and is calculated from the stroke volume [14], and PCa has been associated with prognosis in SSc-PAH patients [5, 9], our results suggest that the prognostic information provided by PCa is largely related to the stroke volume in SSc-PAH. Interestingly, in the study by Campo et al. [5], SVI and PCa at baseline were both associated with survival, but the model performance for each variable was quite similar, suggesting again that stroke volume is the main determinant of prognosis. Given the consistency of the association between SVI and prognosis across several PAH aetiologies, serial measures of SVI might be better in early phase or proof-of-concept clinical trials, to follow after treatment initiation or when assessing prognosis during follow-up.
RV function is a major determinant of survival in PAH and contractile function is worse in SSc-PAH than in IPAH patients for a similar degree of resistive or pulsatile afterload [15, 16]. This discrepancy is likely related to intrinsic systolic dysfunction in SSc-PAH patients due to myocardial fibrosis and impaired sarcomere contractility as compared to IPAH patients [17]. Thus, markers of RV function are essential in assessing risk and therapeutic effect in SSc-PAH patients. While RHC remains mandatory for the diagnosis of PAH, non-invasive measures of RV function would be preferred during follow-up. The SVI and other measures of RV function can be obtained non-invasively with cardiac magnetic resonance imaging (MRI); however, it is not well established whether serial MRI measurements outperforms invasive prognostic variables in SSc-PAH. One study, which included connective tissue disease (CTD)-associated PAH patients, suggested that changes in RV ejection fraction distinguished survivors from non-survivors, whereas changes in cardiac index or PVR did not [18]. There is an increasing body of literature on the prognostic utility of MRI-derived indices of RV function in PAH, including the RV ejection fraction [19, 20], RV end-systolic volume [21], RV end-diastolic volume [22], SVI [22] and RV–pulmonary arterial coupling [21, 23]. A recent MRI study by Swift et al. [21], which included 147 patients with CTD-PAH, found that RV end-systolic volume was prognostic in the overall cohort, but the ratio of RV elastance to arterial elastance (Ea), a load independent measure of RV–pulmonary arterial coupling, was independently associated with outcomes in the CTD-PAH group. In that study, the Ea was calculated from RHC and MRI measurements as (mPAP – PAWP (pulmonary arterial wedge pressure)) divided by SVI. Another study by Vanderpool et al. [23] found that entirely non-invasive MRI estimates of RV–pulmonary arterial coupling (calculated as stroke volume divided by end-systolic volume) were independently associated with survival [23]. Therefore, when repeated invasive haemodynamics are not possible, it may be useful to serially assess SVI, RV–pulmonary arterial coupling and other measures of RV function non-invasively for prognostication with MRI where it is readily available. Three-dimensional echocardiography-derived estimates of RV–pulmonary arterial coupling correlate well with MRI, which may be a future non-invasive prognostic tool [24]. In addition, the tricuspid annular plane systolic excursion derived from echocardiography reflects RV function [25] and has been associated with survival in SSc-PAH patients, albeit in a single small study [26].
PVR was independently predictive of transplant-free survival, although the model fit less well than those evaluating PCa and SVI. Others have also found that baseline PVR independently predicts survival in SSc-PAH [5, 27] although this is inconsistent [28]. Although RAP was not independently predictive in our multivariable models, it was associated in univariable analysis and this is a consistent finding in other studies of SSc-PAH [5, 29], which contrasts with the frequent observation that RAP is an independent prognostic factor in IPAH [10, 30–32]. Interestingly, having a low-risk RAP of <8 mmHg at follow-up was still independent of the other low-risk criteria (online supplementary table S4), and was an important variable in a previous meta-analysis [4], suggesting that a normal RAP is still an appropriate treatment goal for these patients.
Additionally, we confirmed the association between 6MWD and prognosis in SSc-PAH, even though the utility and validity of the 6MWD in SSc-PAH trials has been questioned due to other factors that affect exercise capacity in these patients, such as musculoskeletal limitations [33, 34] and because 6MWD changes correlate poorly with haemodynamic changes with treatment [35]. Nevertheless, in our study, the 6MWD was independently related to prognosis when measured at baseline or during follow-up. This suggests that the 6MWD still provides useful prognostic information, consistent with the meta-analysis by Lefèvre et al. [4]. Another recent study by Gadre et al. [36] derived a clinical score including the 6MWD, Borg dyspnoea score and oxygen saturation during the walk test, which was associated with haemodynamic severity and survival in SSc-PAH. Interestingly, while baseline 6MWD independently predicted long-term outcomes, haemodynamics and NYHA-FC at baseline did not. It may be that initial treatment decisions are more strongly influenced by the severity of haemodynamics or NYHA-FC rather than 6MWD in this population, with more severe haemodynamics or NYHA-FC III–IV invoking more aggressive combination treatment or use of parenteral therapies, which could mask or obviate the prognostic influence of baseline variables. The presence of NYHA-FC IV at follow-up was strongly related to mortality in our study, as would be expected for patients who had persistent severe functional impairment despite treatment.
When we assessed the number of low-risk criteria from the ESC/ERS guidelines [3] in our SSc-PAH population, we found that this simple tool is useful to assess risk as has been recently reported for IPAH, drug-induced and heritable PAH [11]. Similar to several previous studies, baseline assessment of risk did not discriminate short-term survival as well as risk assessment during follow-up, which emphasises the importance of response to treatment over baseline indices of disease severity [7, 11, 37, 38]. It is a sobering reminder of the dismal prognosis for SSc-PAH patients as only 28% of patients achieved three or four low-risk criteria after initial treatment and this “low-risk” group still had a 22% mortality rate at 3 years. Our results are similar to a recent study by Hoeper et al. [38] from the COMPERA registry, which included 347 patients with CTD-PAH at baseline and 213 at follow-up. They also demonstrated that few of the CTD-PAH patients achieved a low-risk profile after initial treatment, and these ‘low-risk’ SSc-PAH patients had worse 5-year (55.5%) survival compared to idiopathic, heritable and drug-induced PAH patients (72.4%) [38]. In contrast to idiopathic, heritable and drug-induced PAH in the French registry [11] and idiopathic PAH patients from the COMPERA registry [39], we did not find that the NT-proBNP/BNP low-risk criteria added prognostic value in SSc-PAH patients. This may argue for the importance of invasive haemodynamic variables in SSc-PAH patients; however, a lack of association could be due to insufficient power as a result of missing data for NT-proBNP or BNP measurements at follow-up or the lack of data on renal function during follow-up, which could confound the association between these biomarkers and outcomes.
We must acknowledge some limitations to this study. It is important to recognise that the initial treatment strategy (monotherapy versus combination therapy) may explain the lack of independent association between baseline variables and outcomes, given that disease severity likely influenced treatment decisions, as has been proposed for other PAH subgroups [10]. The extent of missing data for certain variables that were significant in univariable analysis and which have been previously shown to be related to mortality in SSc-PAH, such as impaired renal function [5], SvO2 [40] and diffusion capacity of the lung for carbon monoxide [6] precluded their inclusion in multivariable analyses. As discussed earlier, the NT-proBNP and BNP analysis is limited by missing data for 46% of the follow-up cohort and lack of adjustment for renal function. Thus it remains unclear whether a non-invasive risk assessment using NYHA, 6MWD and NT-proBNP/BNP in SSc-PAH is accurate. Given a lack of echocardiographic data, we were unable to determine whether some patients had signs or features of left heart disease. Even though patients with PAWP >15 mmHg were excluded, the effect of occult left heart disease on haemodynamic variables cannot be ruled out. Since we limited our analysis to patients without any record of interstitial lung disease (ILD) in the registry, in order to assess a more homogenous population, it is not certain whether these results apply to those with mild ILD. However, another recent study of SSc patients with precapillary pulmonary hypertension from France and the United States suggested that outcomes in patients with a limited extent of ILD are similar to patients without ILD and PAH [41]. Lastly, although we were particularly interested in the value of haemodynamics during follow-up, the indications for RHC within the first year were not always clear and some patients did not undergo a follow-up RHC at all (n=99; figure 1). These patients who did not have any follow-up haemodynamic testing had similar baseline haemodynamics measurements and NYHA-FC as the rest of the cohort, but were significantly older (mean age 70.24 years versus 65.0 years, p<0.001) and had lower baseline 6MWD (193±145 m versus 276±137 m, p<0.001), thus, their exclusion could potentially have introduced a selection bias in our follow-up analysis.
Conclusions
In summary, we found that the haemodynamic variables during follow-up were independently associated with transplant-free survival in a large cohort of patients with SSc-PAH, whereas 6MWD was the only independent predictor at baseline. While traditional haemodynamic measures within the first year of treatment, including cardiac index and mPAP, were related to outcome, the SVI, PCa and PVR may represent more useful haemodynamic indicators of RV function and treatment targets. The number of low-risk criteria present at baseline and first follow-up RHC also predicted long-term outcomes, supporting the utility of this risk assessment tool in the SSc-PAH population. However, despite the fact that nearly a third of patients received combination therapy, only a minority of patients (28%) achieved three or four low-risk criteria and long-term survival remains decidedly poor in this group of patients.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary figures and tables ERJ-00678-2018_Supplement
Acknowledgements
We wish to acknowledge Laurence Rottat (Hôpital Bicêtre, Le Kremlin-Bicêtre, France) for her help in obtaining the data for this study and her hard work in managing data in the French PAH Registry.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Other contributors from the French PH network and registry (alphabetical order): Kaïs Ahmad (Lyon), Fabrice Bauer (Rouen), Emmanuel Bergot (Caen), Laurent Bertoletti (Saint Etienne), Arnaud Bourdin (Montpellier), Mathieu Canuet (Strasbourg), Céline Chabanne (Rennes), Bruno Degano (Besançon), Nicolas Favrolt (Dijon), Gilbert Habib (Marseille, La Timone), Delphine Horeau-Langlard (Nantes), Sylvie Leroy (Nice), Pascal Magro (Tours), Pierre Mauran (Reims), Jean-François Mornex (Lyon), Sylvain Palat (Limoges), Christophe Pison (Grenoble), Patrice Poubeau (La Réunion), Martine Reynaud-Gaubert (Marseille, Nord), Pascal Roblot (Poitiers), Olivier Sanchez (Paris-HEGP), Cécile Tromeur (Brest).
Author contributions: J. Weatherald and A. Boucly contributed equally to literature search, study design, data collection, analysis, figures, interpretation, writing. D. Launay: data collection, interpretation, writing. V. Cottin: data collection, interpretation, writing. G. Prévôt: data collection, interpretation, writing. D. Bourlier: data collection, interpretation, writing. C. Dauphin: data collection, interpretation, writing. A. Chaouat: data collection, interpretation, writing. L. Savale: data collection, interpretation, writing. X. Jaïs: data collection, interpretation, writing. M. Jevnikar: data collection, interpretation, writing. J. Traclet: data collection, interpretation, writing. P. De Groote: data collection, interpretation, writing. G. Simonneau: data collection, interpretation, writing. E. Hachulla: data collection, interpretation, writing. L. Mouthon: data collection, interpretation, writing. D. Montani: data collection, interpretation, writing. M. Humbert and O. Sitbon contributed equally to study design, data collection, analysis, interpretation, writing.
Conflict of interest: J. Weatherald reports grants from the European Respiratory Society and the Canadian Thoracic Society, during the conduct of the study; personal fees and non-financial support from Actelion Pharmaceuticals and Bayer, personal fees from Novartis, grants from Canadian Vascular Network, outside the submitted work.
Conflict of interest: A. Boucly reports non-financial support from GSK, Bayer and MSD, and personal fees and non-financial support from Actelion, outside the submitted work.
Conflict of interest: D. Launay reports grants from Pfizer, during the conduct of the study; and personal fees from Actelion, grants and personal fees from GSK, outside the submitted work.
Conflict of interest: V. Cottin reports personal fees for consultancy, lecturing and travel to medical meetings from Actelion and Roche, personal fees for development of educational presentations, consultancy, lecturing and travel to medical meetings from Boehringer Ingelheim, personal fees for consultancy from Bayer and GSK, personal fees for adjudication committee work from Gilead, personal fees consultancy and travel to medical meetings from MSD, personal fees for consultancy and lecturing from Novartis and Sanofi, institutional grants from Boehringer Ingelheim and Roche, personal fees for chairing a DSMB from Promedior amd Celgene, and personal fees for consultancy and chairing a DSMB from Galapagos, outside the submitted work.
Conflict of interest: G. Prévôt reports personal fees for consultancy, lecturing, development of educational presentations and travel to medical meetings from Actelion, Bayer, Boehringer Ingelheim, GSK and Roche, outside the submitted work.
Conflict of interest: D. Bourlier reports personal fees from Novartis, outside the submitted work.
Conflict of interest: C. Dauphin has nothing to disclose.
Conflict of interest: A. Chaouat reports grants from Direction Générale de l'Offre de Soins (DGOS), during the conduct of the study, personal fees and non-financial support from Actelion Pharmaceuticals, Bayer and MSD, and non-financial support from GSK, outside the submitted work.
Conflict of interest: L. Savale reports grants, personal fees and non-financial support from Actelion Pharmaceuticals, Bayer, GlaxoSmithKline and MSD and grants and personal fees from Pfizer, outside the submitted work.
Conflict of interest: X. Jaïs reports personal fees and non-financial support from Actelion Pharmaceuticals, grants and personal fees from Bayer, and grants, personal fees and non-financial support from MSD and GSK, outside the submitted work.
Conflict of interest: M. Jevnikar has nothing to disclose.
Conflict of interest: J. Traclet reports personal fees from Actelion Pharmaceuticals, Boehringer Ingelheim and Roche, outside the submitted work.
Conflict of interest: P. De Groote has nothing to disclose.
Conflict of interest: G. Simonneau reports grants, personal fees and non-financial support from Actelion Pharmaceuticals, Bayer, GlaxoSmithKline and MSD and grants and personal fees from Pfizer, outside the submitted work.
Conflict of interest: E. Hachulla has nothing to disclose.
Conflict of interest: L. Mouthon has nothing to disclose.
Conflict of interest: D. Montani reports grants and personal fees from Actelion and Bayer, and personal fees from BMS, GSK, MSD and Pfizer, outside the submitted work.
Conflict of interest: M. Humbert has relationships with drug companies including Actelion, Bayer, GSK, Novartis and Pfizer. In addition to being investigator in trials involving these companies, relationships include consultancy service and membership of scientific advisory boards.
Conflict of interest: O. Sitbon reports grants, personal fees and non-financial support from Actelion Pharmaceuticals, Bayer, GlaxoSmithKline and MSD, and grants and personal fees from Pfizer, outside the submitted work.
Support statement: This work was supported by a public grant overseen by the French National Research Agency (ANR) as part of the second Investissements d'Avenir programme (reference ANR-15-RHUS-0002). Jason Weatherald was the recipient of a joint European Respiratory Society/Canadian Thoracic Society Long-Term Research Fellowship (LTRF 2015–4780). Funding information for this article has been deposited with the Crossref Funder Registry.
- Received April 10, 2018.
- Accepted August 12, 2018.
- Copyright ©ERS 2018