- Research article
- Open Access
Prognostic value of perioperative NT-proBNP after corrective surgery for pediatric congenital heart defects
BMC Pediatrics volume 19, Article number: 497 (2019)
It is critically important to assess the prognostic value of NT-proBNP in the form of repeated measures among children undergoing surgery for congenital heart defects (CHD). The aim of the present study is to assess the value of repeated perioperative NT-proBNP in evaluating the time dependent and temporal trajectory in prognostics diagnosis during the perioperative period in a large series of children with CHD.
Repeated measures of NT-proBNP from 329 consecutive children with CHD were obtained before and 1, 12, and 36 h after surgery, respectively. For fully utilizing longitudinal characteristics, we employed parallel cross-sectional logistic regression, a two stage mixed effect model and trajectories over time analysis to mine the predictive value of perioperative NT-proBNP on the binary outcome of prolonged intensive care unit (ICU) stay.
The two stage mixed effects model confirmed that both the mean NT-proBNP level (aOR = 1.46, P = 0.001) and the time trends had prognostic value on the prediction of prolonged ICU stay. In the fully adjusted logistic regression analyses based on gaussian distributions, “rapidly rising NT-proBNP” put the subjects at 5.4-times higher risk of prolonged ICU stay compared with “slowly rising” group (aOR = 5.40, P = 0.003).
Comprehensive assessment of the time dependent and temporal trajectory in perioperative NT-proBNP, indicated by repeated measurements, can provide more accurate identification of children with higher risk of prolonged ICU stay after CHD surgery.
Congenital heart disease (CHD) accounts for nearly 1/3 of all major congenital anomalies (9 per 1000 live births) . Although the improvements in surgical intervention have significantly benefited children with CHD, the incidence of adverse postoperative outcomes remains about 27% . Thus, the identification of prognostic biomarkers for adverse outcome could confer multiple benefits. The guidelines published by the European Society of Cardiology(2016) recommend the clinical biomarkers Brain natriuretic peptide (BNP) and N-terminal pro-B-type natriuretic peptide (NT-proBNP) for diagnosis and prognosis of heart failure . At present, some evidence have pointed out the prognostic value of BNP after surgical intervention for CHD in pediatrics [4,5,6,7]. Compared to BNP, NT-proBNP transportation is more stable, and its half-life in serum is longer than that of BNP . Therefore, NT-proBNP becomes more important and better marker of cardiac strain over time.
Some researches pointed out the prognostic value of NT-proBNP after surgical intervention for CHD. Walsh  et al. found the preoperative NT-proBNP level was a significant predictor of duration of intensive care unit (ICU) stay and peak postoperative level was a predictor of the intensity of overall medical management. They attained NT-proBNP preoperatively and at 2, 12, 24, 48 and 72 h after surgery and just used related analysis to explore the relation between NT-proBNP level and outcome. Goei  et al. demonstrated that the difference in NT-proBNP between pre and postoperative was the strongest independent predictor of cardiac outcome. They applied multivariate Cox regression analyses to evaluate the relation between the change in NT-proBNP levels and the study end point. And Michael  et al. also assessed a reduction in NT-proBNP had a lower subsequent rate of cardiovascular death or HF hospitalization. Qu  et al. have confirmed that time-varying NT-proBNP level, particularly 1-h postoperative levels, had prognostic value on the prediction of outcome after surgery.
Compared with a single measurement, perioperative serial NT-proBNP measurements were more accurate in predicting postoperative adverse events, since it can take into account the hemodynamic stress caused by anesthesia and surgery [10, 13]. However, there is no general consensus on methods for examining data with a longitudinal repeated variable and non-time-varying outcome. Chronological order is often ignored in routine characteristics. In order to overcome the shortcomings from routine characteristics and in-depth dig the values of repeated measurement data, the use of time dependent assessment based on the approach of two stage mixed effects model for clinical trial analyses has recently gained broad supports [14, 15]. This model is more flexible than the previous approaches, and does not require sphericity of the covariance data structure . And this model characterizes within-subject patterns of longitudinal measurements, and the association between features of the longitudinal measurements process and the duration of outcome events. Therefore, our study would use this model to assess NT-proBNP prognostic value of the time dependent effect.
In addition, temporal trajectory of NT-proBNP levels in patients has been confirmed to be the powerful diagnostic and prognostic indicator available in treatment of kidney or heart disease [17,18,19]. Trajectories over time allow to assess how the measured data has changed through two or more assessment points. In this study, we would combined the temporal trajectory analysis to dig more prognostic value from the repeated perioperative NT-proBNP.
Therefore, the aim of the present study is to assess the value of repeated perioperative NT-proBNP in evaluating the time dependent and temporal trajectory in prognostics diagnosis.
Subjects and data collection
This retrospective study was conducted at Guangzhou Women and Children’s Medical Center. We derived the data from the clinical system databases, 364 eligible consecutive patients undergoing CHD surgery at this hospital between June and December 2014. Patients younger than 18 years with CHD who underwent cardiac surgery were eligible for the study. Children who had preoperative arrhythmia (potentially malignant ventricular arrhythmias and high-degree atrioventricular block), rheumatic heart disease, infective endocarditis, myocarditis, pericardial disease, renal dysfunction, or neoplasms, or who could not be separated from cardiopulmonary bypass (CPB) were excluded. Thus, results from 329 cases were presented in this study, which with complete clinical information were included.
Plasma NT-proBNP levels were determined in each patient at preoperative, 1 h, 12 h and 36 h postoperatively, respectively. The NT-proBNP levels were measured with a commercially available fluorescence immunoassay (competitive Enzyme Immuno Assay; ReLIA II, Shenzhen, China) and Multi-Detection Microlpate Reader (VICTOR X5; PerkinElmer, Waltham, Mass). Clinical and biochemical data were collected retrospectively form the medical records, including patient demographics (age, weight, gender), the Risk Adjustment in Congenital Heart Surgery, version 1(RACHS-1) score, Cardiopulmonary bypass (CPB) duration, aorta cross clamp (ACC) time and follow-up measurements after surgery: the duration of Intensive Care Unit (ICU). The primary outcome measure was the prolonged ICU stay, as the short-term adverse outcome to estimate prognostic .. And the patients were divided two groups: ICU stay <= 3 days and ICU stay > 3 days according to the mean duration of ICU stay, which was 3.0 days with a median (interquartile range [IQR]) of (2–5) days (As shown in Additional file 1: Figure S1).
Continuous variables were showed as mean and standard deviation and compared by using the 1-way analysis of variance for 2-group comparisons. Categorical variables were expressed as percentages and compared by using the Fisher exact test or the X2 test. NT-proBNP were compared between the two groups and was log-transformed to obtain normality. A parallel cross-sectional logistic regression model was built to evaluate the relation between NT-proBNP and prolonged ICU stay. Multivariate regression analyses were adjusted for cardiac risk factors and factors recognized to influence NT-proBNP levels: age, gender, RACHS-1, Weight, CPB time and ACC time.
Statistical analysis was performed using R software, version 3.3.2. The package of nlme was used for two stage mixed effects analysis. Mclust package was the optimal model for parameterized Gaussian mixture clustering by hierarchical clustering. The results of the logistic regression model analyses are presented as odds ratios (ORs) with 95% CIs. P < 0.05 considered statistically significant.
In addition, we applied two methods, time dependent assessment and trajectories over time, to explore whether the change of NT-proBNP could predict the postoperative recovery of children patients undergoing CHD surgery. For time dependent assessment, we proposed a two stage mixed effects model as followed. In stage 1, the longitudinally time-varying NT-proBNP level was first modeled using a linear mixed effects (LME) model. In the second stage, best linear unbiased predictor (BLUP) estimated of the random coefficients from LME model were used as predictors in a logistic regression model.
For trajectories over time, the statistical analysis also included two steps. First, subject-specific NT-proBNP levels were clustered based on Gaussian distributions with a predetermined number of clusters (K). Correspondingly, each subject was grouped into one of the K clustered trajectories. Second, to examine the association between trajectory groups over 4 peri-operative assessment points and the binary outcome of prolonged ICU stay, we regressed Yi on this trajectory index along with the covariates as:
where Ci is the clustering trajectory index, which may be more than one-dimensional when K > 2.
Second, we added clusters of NT-proBNP into the logistic regression model adjusted other risk factors to assess the effect of trajectories over time of NT-proBNP with the CHD pediatric prognostics after surgery.
Study group characteristics
The patients’ clinical characteristics were listed in Table 1. Patients were grouped according to length of stay in the ICU: group1, 3 days or less (n = 181; 55.02%) and group 2, more than 3 days (n = 148; 44.98%). With the exception of gender, all the other clinical baseline characteristics showed significant differences between the two groups (all P < 0.001) (Table 2).
The higher level of perioperative NT-proBNP was a potential prognostic factor for longer stay in ICU. In both groups, the level of postoperative NT-proBNP showed a rising trend regardless of preoperative condition and reached the peak at 12 h after operation (Fig. 1).
Prognostic value of time independent assessment of perioperative NT-proBNP on prolonged ICU stay
Parallel cross-sectional logistic regression models for prognostic value of time independent assessment are presented in Table 2. After adjustment for potential confounding factors, the values of NT-proBNP at each of the four times were non-statistically significant associated with prolonged ICU stay (Table 3).
Prognostic value of time dependent assessment of perioperative NT-proBNP on prolonged ICU stay
We used two stage mixed effects model to further investigate the time dependent prognostic value of NT-proBNP. As exhibited in Table 4, patient-specific intercepts and visit-specific slopes were first extracted from Stage 1 and then continuously modeled in Stage 2. The results suggested that there was an association between the patient-specific intercepts and prolonged ICU stay(β = 0.38, SE = 0.11, aOR = 1.46, P = 0.001). Controlled factors mentioned in Table 4 and the time effect, this model showed the higher level of NT-proBNP predicted the higher risk of prolonged ICU stay. Furthermore, compared with before surgery, two effect estimates for both 12 h and 36 h postoperative time points were statistically significant (P = 0.01 and 0.004, respectively). This demonstrated that both the mean NT-proBNP exposure level and the time trends, as indicated by the patient-specific intercept and visit-specific slopes, respectively, were associated with prolonged ICU stay.
The prognostic value of perioperative NT-proBNP levels trajectories on prolonged ICU stay
According to Bayesian information criterion (BIC), the number of clusters chose was two, based on Post−/Pre-(log NT-proBNP) ratio. The average ratios of the NT-proBNP levels within each cluster by time point were displayed in Fig. 2. The trajectory of NT-proBNP levels of cluster 1 elevated slowly over time, whereas rose rapidly in cluster 2. In detail, 87% of subjects fell into cluster 1 and 13% of them fell into cluster 2.
We then modeled each clustering index in relation to odds of suffering a prolonged ICU stay. As indicated in Table 5, compared with the odds for children in cluster 1, the odds of suffering a prolonged ICU stay for children in cluster 2 was 5.4 times (β = 1.67, SE = 0.56, aOR = 5.40, P = 0.003). It means that patients with rising rapidly NT-proBNP levels postoperative might have poorer postoperative recovery.
The principal findings of this study are as follows: (1) the higher perioperative NT-proBNP level prolonged ICU stay for pediatric patients with CDH after surgery; (2) both the mean NT-proBNP exposure level and the time trends were associated with prolonged ICU stay; (3) patients with rising rapidly perioperative NT-proBNP level suffer higher risk of prolonged ICU stay.
Firstly, we carried out assessment of prognostic value within other influence factors (CPB time, ACC time, gender, age, weight and RACHS-1 score). Some studies have pointed out that CPB duration is an independent predictor of adverse outcomes after cardiac surgery . CPB also influenced BNP/NT-proBNP levels proportionally, which was confirmed in our first sub-model of two stage mixed effects model. The inconsistencies among previous findings that prolonged ACC time significantly correlated with major post-operative adverse events , were due to strong correlation between CPB time and ACC time . Here, we found the correlation coefficient was 0.84 (P < 0.001) in our results (Additional file 1: Figure S2). Regardless of the variable of CPB time, length of ACC time restored predictive ability correctly in multiple models (data not shown). RACHS-1 was created to compare in-hospital mortality for children undergoing CHD surgery and demonstrated available for prediction of ICU stay length . To a certain extent, our results supported that, RACHS-1 risk category was positively correlated with length of ICU stay in two stage mixed effect model. Other variables, such as age and gender almost showed not significant association with adverse outcomes.
Increasing evidence pointed out that the change in NT-proBNP have been shown to play key roles as prognostic biomarkers in patient undergoing cardiac surgery for CHD correction [10, 25]. Similar results were reported that the higher NT-proBNP could predict the longer duration of ICU stay  . Our results found that children stay in ICU longer than 3 days showed higher level NT-proBNP at each time point, which confirmed NT-proBNP has predictive value in pediatric CHD patient undergoing surgery. Otherwise, our result showed the peak level of NT-proBNP at 12 h after surgery, which Jiangbo had discuss before . Then we further developed a parallel cross-sectional logistic regression model with covariate, like age or gender, to analysis the predictive value of each time point. However, the result of this model showed that NT-proBNP levels at each time point had no significant effect on the incidence of prolonged ICU stay. According to previous study [11, 12], we assumed that the time effect might play a role. So, we chose two stage mixed effects model to investigate time dependent prognostic value of perioperative NT-proBNP. Results from this model indicated that the mean level of NT-proBNP, especially at the 12 h and 36 h postoperative time points, significantly correlated to duration of ICU stay. Surprisingly, compared to the level of before surgery, NT-proBNP levels at 12 h after surgery, the peak level after surgery, was negatively associated with outcome event. After considering the time effect of the NT-proBNP with prolong.
ed. ICU stay, the peak NT-proBNP level become negative effect with the risk of longer prolonged ICU. It suggested that it is not accurate enough to judge the prognosis of CHD children solely based on the indicators of a single time point, like the peak time point.
So we further used temporal trajectory analysis to dig more prognostic value of temporal trajectories of perioperative NT-proBNP levels in patients. Although NT-proBNP has been confirmed to be the powerful diagnostic and prognostic indicator available in treatment of kidney or heart disease [17,18,19], to our knowledge, this is the first study to assess the prognostic value of temporal trajectory of longitudinal NT-proBNP levels by clustering model in children undergoing cardiac surgery for CHD. When patients undergo surgical procedures, the NT-proBNP levels at specific time points was closely related to individual basal secretion, so using conventional characteristic directly to predict the incidence of heart defects was highly susceptible. In our study, the time trajectory of NT-proBNP was proven to be a better prognostic indicator, since it can reflect an acute change in volume state or left ventricular function that influenced subsequent risk of adverse event. Clinicians might could judge the prognosis of patients by considering the speed of changes of NT-proBNP to adjust their medical decision.
In conclusion, compared with single measurement, repeated measures of perioperative NT-proBNP levels were regarded as a more valuable predictor of ICU stay time length for CHD children undergoing open heart operation. Both the mean NT-proBNP exposure levels and the time trends were associated with prolonged ICU stay, and patients with rising rapidly perioperative NT-proBNP level suffer higher risk of prolonged ICU stay. Further studies are needed to confirm how the time trajectory of NT-proBNP can be used in clinical practice to predict postoperative risk and if repeated measurements of perioperative NT-proBNP should be performed during postoperative management to identify CHD children at high risk of adverse events.
Availability of data and materials
Aorta cross clamp
Best linear unbiased predictor
Brain natriuretic peptide
Congenital heart defects
Cardiopulmonary bypass (CPB) duration
Intensive care unit
Linear mixed effects
N-terminal pro-B-type natriuretic peptide
the Risk Adjustment in Congenital Heart Surgery, version 1
van der Linde D, Konings EE, Slager MA, Witsenburg M, Helbing WA, Takkenberg JJ, Roos-Hesselink JW. Birth prevalence of congenital heart disease worldwide: a systematic review and meta-analysis. J Am Coll Cardiol. 2011;58(21):2241–7.
Parker DM, Everett AD, Stabler ME, Vricella L, Jacobs ML, Jacobs JP, Parikh CR, Pasquali SK, Brown JR. Novel biomarkers improves prediction of 365-day readmission after pediatric congenital heart surgery. Ann Thorac Surg. 2019.
Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JGF, Coats AJS, Falk V, Gonzalez-Juanatey JR, Harjola VP, Jankowska EA, et al. 2016 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure: the task force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC) developed with the special contribution of the heart failure association (HFA) of the ESC. Eur Heart J. 2016;37(27):2129–200.
Hsu JH, Keller RL, Chikovani O, Cheng H, Hollander SA, Karl TR, Azakie A, Adatia I, Oishi P, Fineman JR. B-type natriuretic peptide levels predict outcome after neonatal cardiac surgery. J Thorac Cardiov Sur. 2007;134(4):939–45.
Niedner MF, Foley JL, Riffenburgh RH, Bichell DP, Peterson BM, Rodarte A. B-type natriuretic peptide: perioperative patterns in congenital heart disease. Congenit Heart Dis. 2010;5(3):243–55.
Nahum E, Pollak U, Dagan O, Amir G, Frenkel G, Birk E. Predictive value of B-type natriuretic peptide level on the postoperative course of infants with congenital heart disease. Isr Med Assoc J. 2013;15(5):216–20.
Cantinotti M, Giordano R, Scalese M, Molinaro S, della Pina F, Storti S, Arcieri L, Murzi B, Marotta M, Pak V, et al. Prognostic role of BNP in children undergoing surgery for congenital heart disease: analysis of prediction models incorporating standard risk factors. Clin Chem Lab Med. 2015;53(11):1839–46.
Mir TS, Haun C, Lilje C, Laer S, Weil J. Utility of N-terminal brain natriuretic peptide plasma concentrations in comparison to lactate and troponin in children with congenital heart disease following open-heart surgery. Pediatr Cardiol. 2006;27(2):209–16.
Walsh R, Boyer C, LaCorte J, Parnell V, Sison C, Chowdhury D, Ojamaa K. N-terminal B-type natriuretic peptide levels in pediatric patients with congestive heart failure undergoing cardiac surgery. J Thorac Cardiovasc Surg. 2008;135(1):98–105.
Goei D, van Kuijk JP, Flu WJ, Hoeks SE, Chonchol M, Verhagen HJM, Bax JJ, Poldermans D. Usefulness of repeated N-terminal pro-B-type natriuretic peptide measurements as incremental predictor for long-term cardiovascular outcome after vascular surgery. Am J Cardiol. 2011;107(4):609–14.
Zile MR, Claggett BL, Prescott MF, McMurray JJV, Packer M, Rouleau JL, Swedberg K, Desai AS, Gong JJ, Shi VC, et al. Prognostic implications of changes in N-terminal pro-B-type natriuretic peptide in patients with heart failure. J Am Coll Cardiol. 2016;68(22):2425–36.
Qu JB, Liang HY, Zhou N, Li LJ, Wang YF, Li JB, Cui YQ. Perioperative NT-proBNP level: potential prognostic markers in children undergoing congenital heart disease surgery. J Thorac Cardiov Sur. 2017;154(2):631–40.
Sandhu S, Har BJ, Aggarwal SG, Kavanagh KM, Ramadan D, Exner DV, Investigators R. Predictive value of repeated versus single N-terminal pro B-type natriuretic peptide measurements early after-myocardial infarction. Pace. 2009;32:S86–9.
Donohue MC, Aisen PS. Mixed model of repeated measures versus slope models in Alzheimer's disease clinical trials. J Nutr Health Aging. 2012;16(4):360–4.
Andersen SW, Millen BA. On the practical application of mixed effects models for repeated measures to clinical trial data. Pharm Stat. 2013;12(1):7–16.
Maurissen JP, Vidmar TJ. Repeated-measure analyses: which one? A survey of statistical models and recommendations for reporting. Neurotoxicol Teratol. 2017;59:78–84.
Roberts MA, Hare DL, Sikaris K, Ierino FL. Temporal trajectory of B-type natriuretic peptide in patients with CKD stages 3 and 4, dialysis, and kidney transplant. Clin J Am Soc Nephrol. 2014;9(6):1024–32.
Beltrami M, Ruocco G, Ibrahim A, Lucani B, Franci B, Nuti R, Palazzuoli A. Different trajectories and significance of B-type natriuretic peptide, congestion and acute kidney injury in patients with heart failure. Intern Emerg Med. 2017;12(5):593–603.
Greene SJ, Maggioni AP, Fonarow GC, Solomon SD, Bohm M, Kandra A, Prescott MF, Reimund B, Hua TA, Lesogor A, et al. Clinical profile and prognostic significance of natriuretic peptide trajectory following hospitalization for worsening chronic heart failure: findings from the ASTRONAUT trial. Eur J Heart Fail. 2015;17(1):98–108.
Qu J, Liang H, Zhou N, Li L, Wang Y, Li J, Cui Y. Perioperative NT-proBNP level: potential prognostic markers in children undergoing congenital heart disease surgery. J Thorac Cardiovasc Surg. 2017;154(2):631–40.
Salis S, Mazzanti VV, Merli G, Salvi L, Tedesco CC, Veglia F, Sisillo E. Cardiopulmonary bypass duration is an independent predictor of morbidity and mortality after cardiac surgery. J Cardiothor Vasc An. 2008;22(6):814–22.
Al-Sarraf N, Thalib L, Hughes A, Houlihan M, Tolan M, Young V, McGovern E. Cross-clamp time is an independent predictor of mortality and morbidity in low- and high-risk cardiac patients. Int J Surg. 2011;9(1):104–9.
Nissinen J, Biancari F, Wistbacka JO, Peltola T, Loponen P, Tarkiainen P, Virkkila M, Tarkka M. Safe time limits of aortic cross-clamping and cardiopulmonary bypass in adult cardiac surgery. Perfusion-Uk. 2009;24(5):297–305.
Larsen SH, Pedersen J, Jacobsen J, Johnsen SP, Hansen OK, Hjortdal V. The RACHS-1 risk categories reflect mortality and length of stay in a Danish population of children operated for congenital heart disease. Eur J Cardio-Thorac. 2005;28(6):877–81.
Cuypers JAAE, Eindhoven JA, Slager MA, Opic P, Utens EMWJ, Helbing WA, Witsenburg M, van den Bosch AE, Ouhlous M, van Domburg RT, et al. The natural and unnatural history of the Mustard procedure: long-term outcome up to 40 years. Eur Heart J. 2014;35(25):1666.
The authors thank the doctors and nurses at the Department of Heart Center, Guangzhou Women and Children’s Medical Center, in particular Xinxin Chen, MD, and his team of surgeons.
We are grateful to all of the subjects and their caretakers.
This work was supported by the Guangzhou Institute of Pediatrics/Guangzhou Women and Children’s Medical Center [NO: KCP-2016-002] and the National Natural Science Foundation of China (Grant no.81401755). The first funding supported us to carry out this study including design and data collection. And the second funding mainly provided us sources to learn the statistics methods and scientific research. Both of the funding supported us to accomplish the study and the manuscript.
Ethics approval and consent to participate
This work was approved by the commit of Guangzhou Women and Children’s Medical Center. The use of human subjects in this research has complied with all the relevant national regulations and institutional policies. All participants in the study gave their written informed consent. The research was approved by the author’s equivalent committee.
Consent for publication
All of authors have agreed to publication on BMC Pediatrics.
No conflict of interest exits in the submission of this manuscript. All authors have read and approved this version of the article, and due care has been taken to ensure the integrity of the work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The vioplot of duration of ICU stay among 329 patients. Black box shows the interquartile range, solid white point shows the median of NT-proBNP levels, and black lines extend out from the box are up and down whiskers lines. The external shape with pink color is the kernel density estimation, which presents the population distribution among ICU stay time. Figure S2. The dash line showed CPB duration strongly correlated with AACC time. Cardiopulmonary bypass (CPB) duration strongly correlated with Aortic cross-clamping (ACC) time. Dash line is the regression between CPB duration and ACC time. Correlation coefficient is 0.84 and p value is less than 0.001.
About this article
Cite this article
Lin, F., Zheng, L., Cui, Y. et al. Prognostic value of perioperative NT-proBNP after corrective surgery for pediatric congenital heart defects. BMC Pediatr 19, 497 (2019). https://doi.org/10.1186/s12887-019-1830-y
- Congenital heart defects
- Repeated measures