Abstract
Background Although there is an ever growing number of adult patients with congenital heart disease (ACHD), many are still afflicted by premature death. Previous reports suggested that natriuretic peptides may identify ACHD patients with adverse outcome. We investigated the prognostic power of brain natriuretic peptide (BNP) across the spectrum of ACHD in a large contemporary cohort.
Methods We retrospectively studied 3,392 consecutive and well-characterised ACHD patients under long-term follow-up at a tertiary ACHD centre between 2006-2019. The primary study endpoint was all-cause mortality.
Results A total of 11,974 BNP measurements were analysed. The median BNP at baseline was 47 [24-107] ng/L. During a median follow-up of 8.6 years (29,115 patient-years), 615 (18.1%) patients died. On univariate and multivariate analysis both baseline BNP and temporal changes in BNP levels were predictive of mortality (p<0.001 for both) independent of congenital heart disease diagnosis, complexity, anatomic/haemodynamic features, and/or systolic systemic ventricular function. Patients within the highest quartile of baseline BNP (>107 ng/L) and those within the highest quartile of temporal BNP change (>35 ng/L) had 5.8 and 3.6-fold increased risk of death, respectively.
Conclusion Baseline BNP and temporal BNP changes are both significantly associated with all-cause mortality in ACHD independent of congenital heart disease diagnosis, complexity, anatomic/haemodynamic features, and/or systolic systemic ventricular function. BNP levels represent an easy to obtain and inexpensive marker conveying prognostic information and should be used for the routine surveillance of patients with ACHD.
Introduction
Continuous improvements in the diagnosis and treatment of congenital heart disease (CHD) have brought about an ever-growing population of adults with CHD (ACHD), increasingly including older patients, all with unique healthcare needs. However, ACHD patients are still afflicted by premature death, with heart failure being the leading cause of morbidity and mortality.1 This has important implications, including increasing rates of hospitalisation, higher healthcare utilisation, and higher overall health resource use and costs.2, 3 Risk stratification of all patients at the time of diagnosis and during periodic life-long follow-up is a fundamental part of ACHD practice, aimed at identifying patients at increased risk of adverse outcomes, who require escalation of treatment and closer surveillance.
Natriuretic peptides, either brain natriuretic peptide (BNP) or N-terminal pro-B-type natriuretic peptide (NTproBNP), are established biomarkers in acquired heart failure, used both for diagnosis and assessment of disease severity course. Natriuretic peptide levels and their change over time are also important predictors of heart failure hospitalization and mortality.4, 5
Similar to acquired heart disease, natriuretic peptides are considered clinically relevant in ACHD.6, 7 Previous studies have shown that an elevated single natriuretic peptide value6 and its increase over time is associated with the risk of adverse events in ACHD.7-10 However, variability in their levels across the spectrum of different CHD with different underlying cardiac abnormalities and haemodynamics, exercise limitation, and degree of neurohormonal activation9, 11 limit their implementation as a screening tool for such a heterogenous population. The aim of our study was, therefore, to examine the prognostic role of BNP and its serial measurements in a large, contemporary, and well-characterized cohort of ACHD patients.
Methods
We retrospectively examined data on 3,392 consecutive, clinically stable ACHD patients under periodic long-term follow-up at our large tertiary ACHD centre (Royal Brompton Hospital, London) who had BNP measurements performed between 2006 and 2019. We studied patients’ data including demographics, primary diagnosis, previous interventions, clinical characteristics, laboratory parameters, echocardiographic findings, and cardiopulmonary exercise testing from dedicated electronic health records.12 Exclusion criteria were defined as age < 18 years, acute decompensated heart failure, current systemic infection, or malignancy at the time of BNP sampling, severe renal impairment (serum creatinine level >200 umol/l) and/or liver dysfunction (liver aspartate transaminase more than twice the upper limit of normal (>82 IU/L). Laboratory, echocardiographic, and cardiopulmonary exercise data were linked to clinical visits within 6 months of the respective tests. Systemic ventricular function was assessed by echocardiography and graded in a semi-quantitative manner into normal, mildly, moderately, or severely impaired. Patients were manually categorised according to I) CHD diagnosis, II) disease complexity based on the established Bethesda classification, III) distinct morphologic and hemodynamic features (univentricular heart, systemic right ventricle, Eisenmenger syndrome, Tetralogy of Fallot spectrum, and remaining others with systemic left ventricle who were defined as the reference group). 9, 13 14 Cyanosis was defined as oxygen saturations <90% at rest. New York Heart Association (NYHA) functional class and medications were manually extracted from the clinical records. Heart failure medications were defined by combination usage of more than one of the following medications: β-blocker, angiotensin converting enzyme inhibitor (ACEI) or angiotensin receptor blocker (ARB) or angiotensin receptor neprilysin inhibitor (ARNI), mineralocorticoid antagonist (MRA), and diuretics. The temporal change of BNP levels (ΔBNP) was defined as the difference of BNP level between one measurement level minus the previous BNP level. Data on overall mortality were retrieved from the Office for National Statistics, which registers all United Kingdom deaths, and all-cause mortality served as the study primary endpoint. As this was a retrospective analysis based on data collected for routine clinical care and administrative purposes (UK National Research Ethics Service guidance), individual informed consent was not required. The study was approved by the Institutional Research and Governance Committee.
Statistical Analysis
Continuous variables are presented as mean ± standard deviation (SD) when normally distributed or as median with interquartile range (IQR) if the distribution departed from normality, whereas categorical variables are presented as numbers (%). Comparison of continuous variables between groups was performed using the non-parametric Wilcoxon rank sum test or Kruscal–Wallis test, while categorical variables were compared using the Fisher’s exact test. For the prognostic model, overall survival data was included in univariate and multivariate time dependent Cox proportional hazard survival models after testing the proportional hazards assumption (by assessing the relationship between scaled Schoenfeld residuals and survival time). To assess the effect of change in BNP after the first BNP measurement, stratified by outcome, a random effects mixed linear model was employed using the lme4 and lmerTest packages. Analyses were performed using R-package version 4.1.0. A two-sided P-value of <0.05 was considered indicative of statistical significance.
Results
I. Study population
Baseline demographics, clinical characteristics, and laboratory data are summarised in Table 1; during the study period, 3,392 ACHD patients (male n=1714, 50.5%) had at least one BNP measurement. The mean age at first BNP measurement (baseline) was 40.0 ± 17.0 years, whereas median BMI was 24.5 kg/m2 (IQR 21.9-28.3). Moreover, 2,509 (74.0%), 665 (19.6%), and 218 (6.4%) patients were in the NYHA functional class D, D, and D/D, respectively. There were 1,450 (42.7%) patients with simple, 1,225 (36.1%) moderate, and 717 (21.1%) with complex CHD. During a median follow-up time of 8.6 years (IQR 5.9-11.0, corresponding to a total of 29,115 patient-years of follow-up), 615 (18.1%) patients died.
II. BNP and ΔBNP
Overall, a total of 11,974 BNP measurements were available in the study population with a median BNP of 47 ng/L (IQR 24-107 ng/L). Amongst the 3,392 overall patients with BNP assessment, 2,020 (59.6%) patients had multiple BNP measurements during the study period. The median ΔBNP was 3.5 ng/L (IQR -20 to 35 ng/L) and the median time interval between BNP measurements was 189 days (IQR 89-460 days). The distribution of BNP at baseline and the distribution of ΔBNP is shown in Supplementary Figure 1. There was a significant correlation between BNP levels and NYHA functional class, as well as systemic ventricular function impairment (both p<0.001). Higher BNP levels were seen with increasing disease severity, including worse NYHA functional class and more impaired systemic ventricular function, as shown in Figure 1A and 1B. In contrast, ΔBNP levels did not significantly correlate with change in functional class (p=0.74), but did correlate with worsening systemic ventricular function (p=0.002) as shown in Figure 1C and 1D.
III. Prognostic value of BNP and ΔBNP
BNP at baseline. On univariate Cox analysis age (HR 1.05, 95% CI 1.05-1.06), higher NYHA functional class (NYHA D-D; HR 3.58, 95% CI 3.05-4.21), cyanosis (HR 2.76, 95% CI 2.29-3.33), creatinine concentration (HR 1.02, 95% CI 1.02-1.03), systemic ventricular dysfunction (more than mildly reduced systolic function; HR 2.95, 95% CI 2.19-3.96), treatment with a β-blocker (HR 1.63, 95% CI 1.45-1.83), or heart failure medication (HR 3.59, 95% CI 3.06-4.20), and BNP (HR per /100-ng/L increase 1.16, 95% CI 1.15-1.18) were significantly associated with a higher risk of all-cause mortality (p<0.001 for all). Disease complexity (HR 1.67, 95% CI 1.41-1.97, p=0.001) and Eisenmenger syndrome (HR 2.02, 95% CI 1.70-2.61, p<0.001) were also associated with a higher mortality risk (Table 2).
On multivariate Cox analysis, BNP independently predicted all-cause mortality after adjusting for confounding factors including age, distinct morphological and haemodynamic features, disease complexity, heart rhythm, NYHA functional class D-D, cyanosis, creatinine, systemic ventricular dysfunction, peak VO2, and prescription of β-blocker and heart failure medication (HR per 100-ng/L increase 1.13, 95% CI 1.08-1.18, p<0.001) (Table 2).
The estimated average baseline BNP at study inclusion (intercept) was significantly higher in patients who reached the primary endpoint during the study compared to the remainder (p<0.001), 274 (95%CI 255-293) ng/L vs 90 (95%CI 86-94) ng/L, respectively. Moreover, the average increase of BNP (ΔBNP slope) was significantly higher in patients who reached the primary endpoint compared to the remaining patients, 47 (95%CI 41-53) ng/L versus 3 (95%CI 2-4) ng/L (p<0.001) as presented in Figure 2 based on the results of the mixed effects linear models. When stratifying patients by quartiles of BNP, the risk of all-cause mortality increased with higher BNP level as shown in Figure 3A. Patients within the highest quartile of BNP (>107 ng/L), had a hazard ratio of 5.77 (95% CI 4.91-6.79, p<0.001) of dying compared to the remaining patients (Figure 4).
Serially assessed BNP. On univariate time dependent Cox proportional hazard analysis age (HR 1.06, 95% CI 1.05-1.06), higher NYHA functional class (NYHA D-D; HR 3.47, 95% CI 2.80-4.31), disease complexity (HR 1.53, 95% CI 1.24-1.88), cyanosis (HR 2.68, 95% CI 2.15-3.35), creatinine concentration (HR 1.03, 95% CI 1.02-1.03), systemic ventricular dysfunction (more than mildly reduced systolic function; HR 2.99, 95% CI 2.18-4.10), prescription of heart failure medications (HR 3.37, 95% CI 2.74-4.14) and ΔBNP (HR per /Δ100-ng/L increase 1.22, 95% CI 1.19-1.26) were all significantly associated with a higher risk of all-cause mortality (p<0.001 for all). Eisenmenger syndrome (HR 1.84, 95% CI 1.43-2.37) and prescription of β-blocker (HR 1.27 95% CI 1.03-1.57) were also associated with a higher mortality risk (Table 3).
On multivariate time dependent Cox proportional hazard analysis ΔBNP independently predicted all-cause mortality (HR per Δ100ng/L increase 1.19, 95% CI 1.12-1.26, p<0.001) after adjusting for confounding factors including age, distinct morphological and haemodynamic features, disease complexity, heart rhythm, NYHA functional class D-D, cyanosis, creatinine, systemic ventricular dysfunction, peak VO2, and prescription of β-blocker and HF medication. Similar to BNP, patients within the highest quartile of ΔBNP (>35 ng/L) had the highest risk of mortality (Figure 3B) and a hazard ratio of 3.59 (95% CI 2.93-4.40, p<0.001) of death compared to the remaining (Figure 5).
Discussion
We report herewith the largest study to date regarding baseline and serial BNP measurements in a contemporary cohort of 3,392 ACHD patients with 29,115 patient-years follow-up. Both isolated BNP measurements and temporal BNP change correlated closely with all-cause mortality independently of the underlying CHD diagnosis, disease complexity, specific anatomic features, and systolic ventricular function. Patients with a BNP value in the highest quartile of measurements had an approximately 6-fold increased mortality risk compared to the remainder. A similar effect was seen for change in BNP, where patients within the highest quartile of ΔBNP had 3.6-fold increased risk of death.
BNP and underlying CHD
Our data supports earlier evidence linking natriuretic peptides levels with outcome in ACHD. 7-11 However, potential variability in BNP levels across the spectrum of different CHD diagnoses with different underlying cardiac abnormalities limited the implementation of these earlier reports into clinical practice.13, 15 In our, much larger ACHD contemporary cohort, we demonstrate unequivocally that both BNP and ΔBNP relate to all-cause mortality, independently of underlying congenital heart disease diagnoses, severity and morphological features in the very heterogeneous spectrum of ACHD. Furthermore, our data shows that serial BNP assessment is a valuable adjunct to absolute BNP measurements 7 and in combination, identifies the cohort at the highest risk of mortality.
BNP and NYHA functional class
Patients with ACHD are afflicted by a life-long disease, with multiple adaptations occurring during their trajectory16. Nonetheless, the NYHA classification remains a valuable clinical tool with correlation to exercise capacity, severity of underlying cardiac disease and mid-to long-term mortality.16 Our data showed that baseline BNP values correlate well with NYHA functional class, thus identifying more symptomatic and higher-risk patients. Interestingly, however, there was no correlation with ΔBNP and change in NYHA class, albeit they well identify independently patients at risk of mortality. We submit, therefore, ΔBNP might be more sensitive than NYHA functional class in identifying patients at risk of clinical deterioration; ΔBNP might manifest itself before NYHA functional class change is appreciated, begging for additional medical attention and potential further therapeutic intervention/s.
BNP and systemic ventricular function
ACHD patients have multiple factors influencing their prognosis and, for some, the progression of heart failure. In the follow-up assessment of these patients, we routinely employ standardized transthoracic echocardiography to assess ventricular function, an established prognostic marker.17 Our study demonstrated that baseline BNP levels and ΔBNP were significantly associated with systemic ventricular dysfunction. This is even more relevant in the complex and heterogenous underlying anatomy of CHD patients, as echocardiographic assessment of ventricular function is often not straightforward, requires tertiary expertise, and even then, subtle changes in function can be overlooked. Moreover, serial echocardiograms are time consuming, often subject to interobserver variability, and relatively costly. Cardiac magnetic resonance, the ‘gold standard’ for assessing systolic ventricular function in such a heterogeneous ACHD population is costly and not universally available, thus unlikely to serve as a routine follow up tool. Furthermore there is a real challenge in assessing diastolic function of the systemic, pulmonary, or single ventricle in such a heterogeneous group in a uniform way by any imaging modality, albeit this impacts on prognosis and outcome in ACHD cohorts.18 Hence the value of BNP and ΔBNP - as surrogate markers of global heart function - which our data suggest, are predictive of mortality and correlates with gross assessment of systemic ventricular function. We advocate, therefore, that periodic BNP measurements, being ‘objective’, inexpensive, and easy to obtain, should be employed in the serial assessment of patients at risk of heart failure and adverse outcome from within the ACHD cohort.
Higher BNP and ΔBNP levels and increased mortality risk
Our study showed that patients with the highest BNP and the highest temporal change of BNP levels were at the highest risk of death. The change of natriuretic peptide levels correlates with prognosis in acquired cardiovascular disease,19, 20 and the annual change of natriuretic peptide levels may also be associated with a poorer prognosis in ACHD.10 Our data substantiates this; ACHD patients prone to heart failure (such as patients with a systemic right ventricle or single ventricle) and patients with higher baseline BNP levels warrant serial BNP measurements and closer surveillance. If temporal BNP change is observed, earlier intervention (in the form of up titration of heart failure therapy, bringing forward haemodynamic intervention/s, and/or referral for advance heart failure assessment including transplantation) should be considered. The current study illustrates how an objective and simple marker may be well applied in the setting of a life-long chronic disease such as ACHD.
BNP and its applicability to current ACHD and heart failure guidelines recommendations
It has been suggested that CHD may represent the original heart failure syndrome characterized by a triad comprising cardiac abnormality, exercise limitation, and neurohormonal activation,11 affecting younger patients compared to acquired heart disease cardiovascular cohorts. Different classes of biomarkers have been reported to be associated with adverse events in CHD, including markers of myocardial injury and/or inflammation. 11, 21, 22 Current guidelines for the management of ACHD state that serial testing of natriuretic peptides14 plays a role in identifying patients at risk for adverse events, and our current data unequivocally supports such an approach. Moreover, similar to acquired heart disease, we submit that ACHD patients on a heart failure trajectory should be monitored by serial natriuretic peptide levels according to the current American College of Cardiology (ACC)/American Heart Association (AHA) and European Society of Cardiology (ESC) guidelines for the management of heart failure. 23 24 It could also be of value to obtain a baseline natriuretic peptide level on all patients with ACHD, as even some patients with simple CHD lesions may encounter late complications,25, 26 so that comparisons can be made.
Strength and limitations of the current study
To the best of our knowledge, this is the largest cohort study assessing BNP and its prognostic implication in ACHD. Our study is also the first study where temporal BNP change is shown to be prognostic beyond baseline BNP levels, independently of other established risk factors. Establishing the merits for serial BNP assessment in this population, we submit is a potential game-changer in ACHD clinical practice, namely to assess BNP as a reproducible, prognostic marker with ‘low cost’, easy to measure, and widespread availability.
Although our study provides solid representation of the ACHD spectrum, this remains a single-centre retrospective cohort study. We included patients who had BNP measurements only, as per study design and this might induce a selection bias towards patients with heart failure. Future, prospective studies of neurohormonal biomarkers may validate our data and show the impact of intervention on their levels and, thus outcome. Multicentre, prospective registries may shed additional light on the BNP as a surrogate marker in terms of diagnostics and prognosis of heart failure in ACHD and, in turn, whether BNP reduction improves prognosis in these patients.
Conclusions
Baseline BNP measurements and temporal BNP change in our large, contemporary ACHD cohort were both significantly associated with all-cause mortality independently of underlying diagnosis, disease complexity, specific anatomic and haemodynamic features, including systolic ventricular function. The mortality of patients with BNP values in the highest quartile of measurements (>107 ng/L) was approximately 6-fold higher to the remaining patients. A similar effect was seen for change in BNP, where patients within the highest quartile of ΔBNP had 3.6-fold increased risk of death. BNP measurements are widely available, reproducible, simple and inexpensive to obtain, and should be incorporated in the routine, life-long surveillance of these patients.
Data Availability
The data that support the findings of this study are available from the first author YY, upon reasonable request.
Conflicts of interest
The authors declare no conflicts of interest.
Source of funding
None
Figure Legends
Figure 1. Boxplot of BNP and ΔBNP according to NYHA functional class and systemic ventricular function. A. BNP and NYHA functional class. B. BNP and systemic ventricular dysfunction. C. ΔBNP and NYHA functional class. D. ΔBNP and systemic ventricular dysfunction.
Figure 2. Scatterplot with regression line for the BNP or ΔBNP during follow up period. A. Time point zero is denoted as the date of study inclusion and the first BNP measurement. B. Time point zero is denoted as the date when patients had primary endpoint or when they were censored without primary endpoint.
Figure 3. Kaplan-Meier curves depicting survival free from the all-cause mortality according to the quartile of baseline BNP level and ΔBNP level. A. Kaplan-Meier curves of survival free from all-cause mortality by the quartile of baseline BNP level. Shaded regions indicate 95% CI. B. Kaplan-Meier curves of survival free from all-cause mortality by the quartile of ΔBNP level. Shaded regions indicate 95% CI.
Figure 4. Kaplan-Meier curves depicting survival free from the all-cause mortality according to baseline BNP level. Kaplan-Meier curves of survival free from all-cause mortality by the higher BNP (>107 ng/L) and the lower BNP (≤107 ng/L). Shaded regions indicate 95% CI.
Figure 5. Kaplan-Meier curves depicting survival free from the all-cause mortality according to ΔBNP level. Kaplan-Meier curves of survival free from all-cause mortality by the higher ΔBNP (>35 ng/L) and the lower ΔBNP (≤35 ng/L). Shaded regions indicate 95% CI.
Supplementary Figure 1.
A. Distribution of BNP at baseline. B. Distribution of ΔBNP at serial assessment.
Acknowledgements
None
References
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.