Andrea Sucasas-Alonso1
Sonia Pértega-Díaz2,3*
Vanesa Balboa-Barreiro2,3,4
Fermín García-Muñoz Rodrigo5
Alejandro Avila-Alvarez1
- 1Neonatology Department, Complexo Hospitalario Universitario de A Coruña, A Coruña, Spain
- 2Rheumatology and Health Research Group, Department of Health Sciences, Universidade da Coruña, Ferrol, Spain
- 3Nursing and Health Care Research Group, Instituto de Investigación Biomédica de A Coruña (INIBIC), A Coruña, Spain
- 4Research Support Unit, Complexo Hospitalario Universitario A Coruña, A Coruña, Spain
- 5Division of Neonatology, Complejo Hospitalario Universitario Insular Materno-Infantil, Las Palmas de Gran Canaria, Las Palmas, Spain
Objective: To develop predictive clinical models of bronchopulmonary dysplasia (BPD) through competing risk analysis.
Methods: Retrospective observational cohort study, including preterm newborns ≤32 weeks gestational age, conducted between January 1, 2013 and September 30, 2022 in a third-level Neonatal Intensive Care Unit in Spain. A prediction study was carried out using competing risk models, where the event of interest was BPD and the competing event was death. A multivariate competing risk model was developed separately for each postnatal day (days 1, 3, 7 and 14). Nomograms to predict BPD risk were developed from the coefficients of the final models and internally validated.
Results: A total of 306 patients were included in the study, of which 73 (23.9%) developed BPD and 29 (9.5%) died. On day 1, the model with the greatest predictive capacity was that including birth weight, days since rupture of membranes, and surfactant requirement (area under the receiver operating characteristic (ROC) curve (AUC), 0.896; 95% CI, 0.792–0.999). On day 3, the final predictive model was based on the variables birth weight, surfactant requirement, and Fraction of Inspired Oxygen (FiO2) (AUC, 0.891; 95% CI, 0.792–0.989).
Conclusions: Competing risk analysis allowed accurate prediction of BPD, avoiding the potential bias resulting from the exclusion of deceased newborns or the use of combined outcomes. The resulting models are based on clinical variables measured at bedside during the first 3 days of life, can be easily implemented in clinical practice, and can enable earlier identification of patients at high risk of BPD.
1 Introduction
Bronchopulmonary dysplasia (BPD) is a common chronic lung disease in premature infants that leads to significant morbidity and mortality (1–3). Although multiple studies have sought to develop effective preventive interventions for BPD, including postnatal corticosteroids and avoidance of invasive mechanical ventilation (4), the incidence of BDP does not appear to have declined in recent years (5, 6).
Inflammation and oxidative stress that ultimately lead to BPD appear to occur early in the neonatal period, preceding clinical symptoms of BPD and diagnosis based on current definitions (7). Moreover, the description of different BPD phenotypes highlights the need to individualize care and, potentially, to develop more effective preventive strategies (8). Therefore, early and accurate identification of babies who will later develop BPD is among the highest priorities in modern neonatal care. Good-quality predictive models for BPD aim to establish an individual risk of BPD, helping clinicians to identify high-risk patients and individualize care, and hopefully enabling more effective preventive strategies. While a plethora of predictive models have been developed in recent decades (9), none are fully incorporated into routine clinical practice (5, 10) and all suffer from a high risk of bias (11).
Heterogeneity in clinical practices across units and countries, the unavailability of certain variables used in models, the complexity of some models, their development based on population data from decades ago, and a lack of external validation are just some of the factors that may have contributed to the apparent failure in developing an optimal prediction tool. Furthermore, certain methodological issues need to be critically evaluated. Most models were developed using logistic regression analysis of datasets that either exclude infants who died before BPD diagnosis, or use a combined outcome of mortality and BPD (12–16).
This approach attributes equal weight to outcomes that are clearly different for both clinicians and parents (death and BPD) and, more problematically, analyzes predictors of BPD in a lower-risk infant population, eliminating from the analysis those children who die before 36 weeks of postmenstrual age (PMA). Although BPD is usually defined at a fixed point in time (typically 36 weeks PMA), it could be considered a time-to-event outcome, diagnosed at some point after birth, meaning that it could be correctly studied using survival analysis methods. Competing risks analysis, a special type of survival analysis, is particularly useful when patients are at risk of more than one mutually exclusive event, such as BPD and death (17).
In the present study we sought to develop an up-to-date model for early prediction of BPD, based on routinely collected variables, and to explore the utility of competing risk survival analysis methods in this context. We hypothesized that this methodology may constitute a superior means of studying the risk of BPD, overcoming the methodological limitations of other approaches.
2 Materials and methods
This single-center retrospective analysis of data recorded in a prospective registry was conducted between January 1, 2013 and September 30, 2022, at a third-level Neonatal Intensive Care Unit of a hospital of the Spanish Public Health System. Preterm infants ≤32 weeks gestational age (GA) were included in the study. Newborns with major congenital anomalies were excluded.
Variables related to prenatal and obstetric characteristics, demographics, resuscitation in the delivery room, respiratory support, and neonatal evolution during first hospital admission were collected. The main prematurity-related outcomes were recorded.
Gestational age was estimated based on the last menstrual period date and/or obstetric and ultrasound parameters recorded in the maternal medical record. BPD was defined as the need for supplemental oxygen therapy for 28 days of life and was classified based on oxygen requirements and respiratory support at 36 weeks of PMA (18). Patent ductus arteriosus was diagnosed by cardiac ultrasound and managed according to local protocols based on echocardiographic and clinical data (only ductus >1.5 mm in patients on respiratory support were considered) (19). Maternal chorioamnionitis was defined by clinical diagnosis in the medical record or by histologic diagnosis by placental pathology (20–22). Nosocomial infection was defined as positive microbial growth on one or more bloodstream cultures or any sterile body fluid obtained after 72 h of life with accompanying clinical signs of sepsis (23). NEC grade ≥2 according to the Bell classification was considered (24, 25). Intraventricular hemorrhage (IVH) was defined and graded according to Volpe (26, 27). All newborns were screened for retinopathy of prematurity and its grade was classified according to international guidelines (28). Intrauterine growth restriction was described as birth weight ≤1.5 z-score according to Fenton growth charts (29).
There were standard guidelines for respiratory management in place in the NICU during the study period, which were in keeping with European and national recommendations (30, 31). In brief, stabilization in the delivery room as initiated with mask and T-resuscitator applying CPAP or NIPPV depending on the presence of spontaneous breathing and heart rate. NIPPV was used in the NICU as early rescue therapy before considering orotracheal intubation or as initial support at clinical discretion. Surfactant was administered using the INSURE (intubation-surfactant-extubation) procedure in the first hours of life when FiO2 requirements were greater than 30% on NIV or if the infant required orotracheal intubation for any other reason. For the INSURE technique, premedication with caffeine, fentanyl and atropine was administered. Caffeine citrate was given prophylactically to all infants <28 weeks GA and in those <32 weeks who developed apneas (20 mg/Kg bolus followed by 5 mg/Kg per day, which can be increased to 10 mg/Kg if clinically needed). Hydrocortisone administration was considered in infants on IMV and FiO2 >0.3 after the first week of life. According with the European guidelines, the oxygen saturation target was established between 90% and 94%.
No relevant changes were made in clinical practice during the study period.
The study was approved by the local Research Ethics Committee (code 2017/360).
2.1 Statistical analysis
Maternal and perinatal characteristics were described and compared between neonates without BPD, neonates with BPD, and neonates who died. A BPD prediction study was carried out using competing risk analysis, where the event of interest was BPD and the competing event was death (17). PMA was considered as the time variable, and each newborn was followed until death or hospital discharge.
The cumulative incidence function (CIF) of both BPD and death was estimated, using the method proposed by Kalbfleish and Prentice (32). Fine-Gray univariate and multivariate competing risk models were used to identify predictive factors for BPD and death separately (17). Variables considered clinically important and those for which a statistically significant difference (p < 0.05) was observed in the univariate analysis were included in the multivariate models, using Akaike information criterion (AIC) as the selection criterion in a stepwise forward selection strategy. A multivariate competing risk model was developed separately for each timepoint (days 1, 3, 7 and 14) using data from infants who survived to the prediction day. For each day, the final prediction model was the simplest model that included variables significantly associated with BPD among those with the greatest predictive ability.
Nomograms for predicting BPD risk were developed from the coefficients of the final models and were internally validated. The primary event was BPD (defined as the need for supplemental oxygen therapy for 28 days of life), and the estimated cumulative probability of BPD at 36 weeks was considered the prediction of interest, since all the babies who develop BPD do so before this time horizon. The time-dependent area under the receiver operating characteristic (ROC) curve (AUC) for t = 36 weeks was used to quantify the discrimination performance of the competing risk nomograms (33), since the concordance index is inadequate when the aim is to predict the risk of an event at a fixed timepoint. A 10-fold cross-validation was used to obtain an optimism-adjusted AUC.
To further analyze the models' discriminative ability, newborns were divided into four risk groups based on their estimated probability of BPD according to both nomograms. Gray's test was used to compare the cumulative incidence of BPD between groups (34). Bootstrap cross-validation calibration plots (B = 100 resamples) were also used to measure the degree of consistency between predicted and observed BPD probabilities. Finally, the clinical utility and overall advantages of the nomograms were assessed using decision curve analysis (DCA) (35).
All analyses were performed using the statistical software SPSS v.28 and R v4.2. The “cmprsk”, “riskRegression”, and “pec” packages were used for competing risk analysis and internal validation of the models. Nomograms were built with the “QHScrnomo” package. DCA was implemented using the “dcurves” package.
All reported significance levels are two-sided. Statistical significance was defined as a p-value < 0.05.
3 Results
During the study period, 322 preterm newborns ≤32 weeks GA born at the study hospital were evaluated. Sixteen patients with major congenital anomalies were excluded. Ultimately, 306 patients were included in the study, of whom 73 (23.9%) developed BPD and 29 died (9.5%) (Figure 1, Supplementary Figure S1).
Figure 1. Flowchart of study population selection and diagnosis of bronchopulmonary dysplasia (BPD).
Mean GA and birth weight were 29.3 ± 2.3 weeks and 1,135.8 ± 309.7 g, respectively. Table 1 shows the most relevant characteristics of the study population. Results of the univariate analysis of factors associated with either BPD or death are shown in Table 2 and Supplementary Table S1.
Table 1. Maternal and neonatal characteristics of the study cohort.
Table 2. Univariate competing-risk regression to identify risk factors for bronchopulmonary dysplasia (BPD).
Table 3 shows the final models for BPD prediction. On day 1 the greatest predictive capacity was obtained using the variables birth weight, days since rupture of membranes, and surfactant requirement. On day 3, the final predictive model was established using the variables birth weight, surfactant requirement, and FiO2. For days 7 and 14, only birth weight and surfactant requirement were statistically significant in the multivariate analysis. In particular, FiO2 on day of life (DOL) 7 and 14 did not achieve statistical significance.
Table 3. Multivariate competing risk analysis of bronchopulmonary dysplasia (BPD) for different postnatal days of life.
Figure 2 shows the nomograms developed based on the multivariate competing risk models for postnatal days 1 and 3. The AUC scores for the predictions at day of birth and DOL 3 were 0.896 (95% CI, 0.792–0.999) and 0.891 (95% CI, 0.792–0.989), respectively. Optimism adjusted AUC values for the same models were 0.881 (95% CI, 0.767–0.995) and 0.863 (95% CI, 0.760–0.967), showing good predictive accuracy (Supplementary Figures S2A,B), especially in babies 28–32 weeks GA, being the discriminative ability of the models lower for children <28 weeks GA (Supplementary Table S2).
Figure 2. Nomograms predicting bronchopulmonary dysplasia (BPD) at postnatal days 1 (A) and 3 (B). Specific values for individual patients are located on each variable axis. A line is drawn upward to the “Points” axis to determine the number of points corresponding to each variable. The sum of points for all variables is calculated, and a line drawn downwards to the probability axis to determine the risk of BPD.
For each of the nomograms, the estimated probability of BPD was obtained for all cases in the database, and newborns were categorized into four risk groups according to the quartiles of the estimated BPD probabilities. Differences between CIF of BPD curves by risk subgroups were statistically significant (Supplementary Figure S3). The calibration curves of both nomograms are also shown in Supplementary Figures 2C,D. Subgroup analysis (Supplementary Figure S3) showed that the observed percentages of BPD cases were consistent with the estimated probabilities.
The decision curves of the competing risk nomograms are also displayed in Supplementary Figures 2E,F.
4 Discussion
The competing risk analysis presented here shows that BPD can be effectively predicted using simple models that include clinical variables collected at the bedside during the first 3 days of life. To our knowledge, this is the first study to use this methodology to develop predictive models for BPD.
The importance of the availability of good prediction tools in neonatology has been previously highlighted, both for clinical practice and research (36). However, despite the efforts made by researchers and clinicians, BPD is still a major concern for health care providers and families, partly due to a lack of widely used prediction tools. Critical analysis of previous research on predictive models identified the exclusion of infants that die before BPD diagnosis and the use of combined outcomes as areas of potential improvement (10).
The explanation of this problem entails the introduction of two statistical concepts: time-to-event analysis and competing risks. In many studies, the results are assessed longitudinally, and each patient in the cohort is evaluated for a period of time until the event occurs. The study objectives can be to estimate the probability of occurrence of the event or its association with other variables of interest. An event that hinders or modifies the possibility of observing the event of interest is known as a competing event (37). In the case of BPD, patients who die prior to BPD diagnosis are no longer at risk and cannot be treated either as individuals with BPD or as healthy patients, and therefore death constitutes a competing risk for BPD.
Excluding patients who die, as was previously done (14), may not be correct. In fact, the most serious cases of lung disease of prematurity are probably found among those who died in the first days of life, and therefore may have been candidates for early treatment. For this reason, some studies include deaths from respiratory causes among patients with severe BPD (38). However, this also may be methodologically incorrect, as some patients are assigned a theoretical and somewhat arbitrary diagnosis. By contrast, other studies include infants who die in the non-BPD group (39), and consider them as members of the “healthy” controls.
Importantly, a recent systematic review and meta-analysis of BPD prediction models (10) discussed the importance of taking BPD and mortality into account as competing risks, although the authors propose using combined outcomes to address this limitation. In our opinion, the use of combined outcomes results in the attribution of common risk factors to two very different events. The literature includes studies of risk-factors for mortality in premature infants, the results of which differ from those focusing on BPD risk-factors (40, 41). This indicates that, although both variables share several common risk factors, there are many others that are unique to one or the other, since there are multiple causes of death that do not share the same pathophysiological pathways of BPD (40–42).
Competing risk analysis is a special type of survival analysis that aims to correctly estimate the marginal probability of an event in the presence of competing events. To the best of our knowledge, this methodology has not been previously used to predict BPD. A similar analysis focusing on another relevant complication of preterm birth, retinopathy of prematurity, was carried out by Miller et al. (43). This disease has similarities with BPD, since its development is determined by variables in the first days of life, before diagnosis is made weeks later (diagnosis cannot be established in patients who die earlier) (44). Moreover, in a recent clinical trial of hydrocortisone in BPD, the combined variable death or BPD at 36 weeks was used as the main study outcome, but an analysis using competitive risk models was also planned by the investigators, showing the relevance of this methodology in neonatal clinical trials (45).
The results of our models showed that predictive capacity was best established on DOL 1 using simple and readily accessible variables: birth weight, surfactant requirement, and time since rupture of membranes. On DOL 3, the variables ultimately included in the model were birth weight, surfactant requirement, and FiO2. The high discriminatory capacity of the proposed models, in terms of AUC values, could be surprising given the reduced number of variables included. It should be noted that the use of a competing risks methodology could provide a more unbiased estimate of the probability of BPD and the HR associated with each of the covariates considered, which may lead to a better predictive ability. On the other hand, all variables finally included in the models have been previously described in the literature as risk factors for the disease and are biologically plausible (9). In this respect, we would like to highlight the simplicity of our models, which consist of a limited number of objective clinical variables, recorded at bedside during the first days of life. Our data show that this simplicity does not imply inferior predictive capacity relative to more complex models that include more variables or biochemical or even ultrasound markers (46–49).
The aforementioned systematic review by Peng et al. provided some recommendations for the development of future models (9). The authors highlighted the importance of models using variables that are recorded sufficiently early, thereby enabling initiation of preventive treatments and/or recruitment of patients in clinical trials.
The present findings are clinically relevant, as the sooner high-risk patients are identified, the sooner preventive treatment can be instituted. Hallmark studies indicate that the benefits of treatment outweigh the risks in high-risk patients (50) when the estimated probability of experiencing BPD is high [>40% (95% CI 33%–46%)] (51). Our findings indicate that if a patient has received surfactant in the first hours of life, the individual risk of BPD can be effectively predicted based on birth weight and time since rupture of membranes or FiO2 on DOL 3 (Figure 2), without needing to wait for subsequent respiratory progression, when many treatments would be futile and heterogeneity in clinical practice may be greater.
Limitations of our study include the smaller sample size relative to other studies: this implies a limited number of events, making it difficult to explore the utility of more complex models, and may explain why we did not obtain better results on days 7 and 14. Moreover, our study was carried out using data from a retrospective cohort from a single center, without external validation. In particular, the sample studied is relatively mature with a mean gestational age of 29.3 weeks, so results could be not generalizable to more immature babies, as suggested by the AUC values obtained in the subset of babies <28 weeks GA, albeit with a very small sample size. This may also be related to the results of the calibration of the models, which show better performance in children at low risk of BPD and poorer calibration in babies at medium to high BPD risk, probably coinciding with those with a lower GA. External validation is currently underway using a national population-based cohort, and we hope to have data in the coming months. Finally, some limitations of the proposed methodology cannot be overlooked. A time-to-event analysis with PMA (which includes GA) as the time variable prevents us from analysing the effect of GA as an additional covariate in the proposed models. However, it can be adjusted for birth weight, which is a closely related variable. On the other hand, the proposed models include variables collected on specific days of life that may correspond to different PMA for different children, so it would be interesting to explore how to incorporate these time-dependent variables in a more efficient way using an analogous methodology.
We believe that the methodology proposed and the results produced using simple models based on early clinical data could be of great interest to the scientific community in general, and the field of clinical neonatology in particular. In conclusion, our results show that BPD can be effectively predicted using models based on variables that can be measured during the first hours of life, including birth weight, surfactant requirement, time since membrane rupture, and FiO2 in the first 3 days of life. Further research is needed to externally validate our findings and continue the quest to develop high-quality tools to predict BPD.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by A Coruña-Ferrol Territorial Research Ethics Committee (code 2017/360). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
AS: Conceptualization, Investigation, Resources, Writing – original draft, Writing – review & editing. SP: Supervision, Writing – review & editing, Data curation, Formal Analysis, Methodology, Software, Visualization, Writing – original draft. VB: Writing – review & editing, Data curation, Formal Analysis, Methodology, Software. FG: Supervision, Writing – review & editing. AA: Supervision, Writing – original draft, Writing – review & editing, Conceptualization, Funding acquisition, Resources.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article.
This study was funded by Instituto de Salud Carlos III (ISCIII) through the project PI21/01601 and co-funded by the European Union.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2024.1335891/full#supplementary-material
References
1. Lui K, Lee SK, Kusuda S, Adams M, Vento M, Reichman B, et al. Trends in outcomes for neonates born very preterm and very low birth weight in 11 high-income countries. J Pediatr. (2019) 215:32–40.e14. doi: 10.1016/j.jpeds.2019.08.020
2. Cheong JLY, Doyle LW. An update on pulmonary and neuro- developmental outcomes of bronchopulmonary dysplasia. Semin Fetal Neonatal Med. (2018) 42:478–84. doi: 10.1053/j.semperi.2018.09.013
3. Jensen EA, Edwards EM, Greenberg LT, Soll RF, Ehret DEY, Horbar JD. Severity of bronchopulmonary dysplasia among very preterm infants in the United States. Pediatrics. (2021) 148(1):1–8. doi: 10.1542/peds.2020-030007
4. Doyle LW. Postnatal corticosteroids to prevent or treat bronchopulmonary dysplasia. Neonatology. (2021) 118(2):244–51. doi: 10.1159/000515950
5. Avila-Alvarez A, Zozaya C, Pértega-Diaz S, Sanchez-Luna M, Iriondo-Sanz M, Elorza MD, et al. Temporal trends in respiratory care and bronchopulmonary dysplasia in very preterm infants over a 10-year period in Spain. Arch Dis Child Fetal Neonatal Ed. (2022) 107(2):143–9. doi: 10.1136/archdischild-2021-322402
6. Geetha O, Rajadurai VS, Anand AJ, Dela Puerta R, Huey Quek B, Khoo PC, et al. New BPD-prevalence and risk factors for bronchopulmonary dysplasia/mortality in extremely low gestational age infants ≤28 weeks. J Perinatol. (2021) 41(8):1943–50. doi: 10.1038/s41372-021-01095-6
7. Leroy S, Caumette E, Waddington C, Hébert A, Brant R, Lavoie PM. A time-based analysis of inflammation in infants at risk of bronchopulmonary dysplasia. J Pediatr. (2018) 192:60–5.e1. doi: 10.1016/j.jpeds.2017.09.011
8. Pierro M, Van Mechelen K, van Westering-Kroon E, Villamor-Martínez E, Villamor E. Endotypes of prematurity and phenotypes of bronchopulmonary dysplasia: toward personalized neonatology. J Pers Med. (2022) 12(5):687. doi: 10.3390/jpm12050687
9. Peng HB, Zhan YL, Chen Y, Jin ZC, Liu F, Wang B, et al. Prediction models for bronchopulmonary dysplasia in preterm infants: a systematic review. Front Pediatr. (2022) 10(May):1–12. doi: 10.3389/fped.2022.856159
10. Kwok TC, Batey N, Luu KL, Prayle A, Sharkey D. Bronchopulmonary dysplasia prediction models: a systematic review and meta-analysis with validation. Pediatr Res. (2023) 94(1):43–54. doi: 10.1038/s41390-022-02451-8
11. Romijn M, Dhiman P, Finken MJJ, van Kaam AH, Katz TA, Rotteveel J, et al. Prediction models for bronchopulmonary dysplasia in preterm infants: a systematic review and meta-analysis. J Pediatr. (2023) 258:113370. doi: 10.1016/j.jpeds.2023.01.024
12. Laughon MM, Langer JC, Bose CL, Smith PB, Ambalavanan N, Kennedy KA, et al. Prediction of bronchopulmonary dysplasia by postnatal age in extremely premature infants. Am J Respir Crit Care Med. (2011) 183(12):1715–22. doi: 10.1164/rccm.201101-0055OC
13. Baud O, Laughon M, Lehert P. Survival without bronchopulmonary dysplasia of extremely preterm infants: a predictive model at birth. Neonatology. (2021) 118(4):385–93. doi: 10.1159/000515898
14. Greenberg RG, McDonald SA, Laughon MM, Tanaka D, Jensen E, Van Meurs K, et al. Online clinical tool to estimate risk of bronchopulmonary dysplasia in extremely preterm infants. Arch Dis Child Fetal Neonatal Ed. (2022) 107(6):638–43. doi: 10.1136/archdischild-2021-323573
15. Valenzuela-Stutman D, Marshall G, Tapia JL, Mariani G, Bancalari A, Gonzalez Á. Bronchopulmonary dysplasia: risk prediction models for very-low- birth-weight infants. J Perinatol. (2019) 39(9):1275–81. doi: 10.1038/s41372-019-0430-x
16. Baker EK, Davis PG. Bronchopulmonary dysplasia outcome estimator in current neonatal practice. Acta Paediatr Int J Paediatr. (2021) 110(1):166–7. doi: 10.1111/apa.15427
17. Fine JP, Gray R. A proportional hazards regression model for the sub-distribution with covariates adjusted censoring weight for competing risks data. J Am Stat Assoc. (1999) 94:496–509. doi: 10.1080/01621459.1999.10474144
18. Jobe AH, Bancalari E. NICHD/NHLBI/ORD workshop summary. Am J Respir Crit Care Med. (2001) 163:1723–9. doi: 10.1164/ajrccm.163.7.2011060
19. Avila-Alvarez A, Serantes Lourido M, Barriga Bujan R, Blanco Rodriguez C, Portela-Torron F, Bautista-Hernandez V. Cierre quirúrgico del ductus arterioso persistente del prematuro: ¿influye la técnica quirúrgica en los resultados? (Surgical closure of patent ductus arteriosus in premature neonates: does the surgical technique affect the outcome?). An Pediatr (Barc). (2017) 86(5):277–83. (Spanish). doi: 10.1016/j.anpedi.2015.12.011
20. Gibbs RS, Blanco JD, St Clair PJ, Castaneda YS. Quantitative bacteriology of amniotic fluid from women with clinical intraamniotic infection at term. J Infect Dis. (1982) 145(1):1–8. doi: 10.1093/infdis/145.1.1
21. Kim CJ, Romero R, Chaemsaithong P, Chaiyasit N, Yoon BH, Kim YM. Acute chorioamnionitis and funisitis: definition, pathologic features, and clinical significance. Am J Obstet Gynecol. (2015) 213(4 Suppl):S29–52. doi: 10.1016/j.ajog.2015.08.040
22. Kachikis A, Eckert LO, Walker C, Bardají A, Varricchio F, Lipkind HS, et al. Chorioamnionitis: case definition & guidelines for data collection, analysis, and presentation of immunization safety data. Vaccine. (2019) 37(52):7610–22. doi: 10.1016/j.vaccine.2019.05.030
23. Wynn JL. Defining neonatal sepsis. Curr Opin Pediatr. (2016) 28(2):135–40. doi: 10.1097/MOP.0000000000000315
24. Bell MJ, Ternberg JL, Feigin RD, Keating JP, Marshall R, Barton L, et al. Neonatal necrotizing enterocolitis. Therapeutic decisions based upon clinical staging. Ann Surg. (1978) 187:1–7. doi: 10.1097/00000658-197801000-00001
25. Neu J. Necrotizing enterocolitis: the search for a unifying pathogenic theory leading to prevention. Pediatr Clin North Am. (1996) 43(2):409–32. doi: 10.1016/s0031-3955(05)70413-2
26. Intracranial hemorrhage: germinal matrix intraventricular hemorrhage. In: Volpe JJ, editor. Neurology of the Newborn. 5th Edn. Philadelphia, PA: Saunders Elsevier (2008). p. 517–88.
27. Volpe JJ. Brain injury in the premature infant: overview of clinical aspects, neuropathology, and pathogenesis. Semin Pediatr Neurol. (1998) 5(3):135–51. doi: 10.1016/s1071-9091(98)80030-2
28. Section on Ophthalmology American Academy of Pediatrics; American Academy of Ophthalmology; American Association for Pediatric Ophthalmology and Strabismus. Screening examination of premature infants for retinopathy of prematurity. Pediatrics. (2006) 117(2):572–6. doi: 10.1542/peds.2005-2749 [Erratum in: Pediatrics. 2006 118(3):1324].16452383
29. Fenton TR, Kim JH. A systematic review and meta-analysis to revise the Fenton growth chart for preterm infants. BMC Pediatr. (2013) 13:59. doi: 10.1186/1471-2431-13-59
30. Sweet DG, Carnielli VP, Greisen G, Hallman M, Klebermass-Schrehof K, Ozek E, et al. European consensus guidelines on the management of respiratory distress syndrome: 2022 update. Neonatology. (2023) 120(1):2–23. doi: 10.1159/000528914
31. Zeballos Sarrato G, Avila-Alvarez A, Escrig Fernández R, Izquierdo Renau M, Ruiz Campillo CW, Gómez Robles C, et al. Guía española de estabilización y reanimación neonatal 2021. Análisis, adaptación y consenso sobre las recomendaciones internacionales. An Pediatría. (2022) 96(2):145.e1–9. doi: 10.1016/j.anpedi.2021.06.003
32. Kalbfleisch JD, Prentice R. The Statistical Analysis of Failure Time Data. New York: John Willey and Sons (1980).
33. Blanche P, Kattan MW, Gerds TA. The c-index is not proper for the evaluation of t-year predicted risks. Biostatistics. (2019) 20(2):347–57. doi: 10.1093/biostatistics/kxy006
34. Gray RJ. A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann Stat. (1988) 16:1141–54. doi: 10.1214/aos/1176350951
35. Vickers AJ. Decision analysis for the evaluation of diagnostic tests, prediction models and molecular markers. Am Stat. (2008) 62(4):314–20. doi: 10.1198/000313008X370302
36. Crilly CJ, Haneuse S, Litt JS. Predicting the outcomes of preterm neonates beyond the neonatal intensive care unit: what are we missing? Pediatr Res. (2021) 89(3):426–45. doi: 10.1038/s41390-020-0968-5
37. Pintilie M. An introduction to competing risks analysis. Rev Esp Cardiol. (2011) 64(7):599–605. doi: 10.1016/j.recesp.2011.03.017
38. Fairchild KD, Nagraj VP, Sullivan BA, Moorman JR LD. Oxygen desaturations in the early neonatal period predict development of bronchopulmonary dysplasia. Pediatr Res. (2019) 85(7):987–93. doi: 10.1038/s41390-018-0223-5
39. Ryan SW, Nycyk J, Shaw BN. Prediction of chronic neonatal lung disease on day 4 of life. Eur J Pediatr. (1996) 155:668–71. doi: 10.1007/BF01957150
40. del Río R, Thió M, Bosio M, Figueras J, Iriondo M. Predicción de mortalidad en recién nacidos prematuros. Revisión sistemática actualizada. An Pediatría. (2020) 93(1):24–33. doi: 10.1016/j.anpedi.2019.11.003
41. Iriondo M, Thio M, Del Rió R, Baucells BJ, Bosio M, Figueras-Aloy J. Prediction of mortality in very low birth weight neonates in Spain. PLoS One. (2020) 15(7 July 2020):1–13. doi: 10.1371/journal.pone.0235794
42. Ge WJ, Mirea L, Yang J, Bassil KL, Lee SK, Shah PS. Prediction of neonatal outcomes in extremely preterm neonates. Pediatrics. (2013) 132(4):e876–85. doi: 10.1542/peds.2013-0702
43. Miller MM, Revenis ME, Lai MM, Meleth AD, Jeffress ES, Carrera A, et al. Risk and clinical course of retinopathy of prematurity in 78 infants of gestational age 22–25 weeks. JAAPOS. (2014) 18(3):266–70. doi: 10.1016/j.jaapos.2014.01.016
44. Castro Conde JR, Echániz Urcelay I, Botet Mussons F, Pallás Alonso CR, Narbona E, Sánchez Luna M. Retinopathy of prematurity. Prevention, screening and treatment guidelines. An Pediatr. (2009) 71(6):514–23. doi: 10.1016/j.anpedi.2009.07.006
45. Onland W, Merkus MP, Nuytemans DH, Jansen-van der Weide MC, Holman R, van Kaam AH, et al. Systemic hydrocortisone to prevent bronchopulmonary dysplasia in preterm infants (the SToP-BPD study): statistical analysis plan. Trials. (2018) 19(1):1–10. doi: 10.1186/s13063-018-2505-y
46. Pezza L, Alonso-Ojembarrena A, Elsayed Y, Yousef N, Vedovelli L, Raimondi F, et al. Meta-analysis of lung ultrasound scores for early prediction of bronchopulmonary dysplasia. Ann Am Thorac Soc. (2022) 19(4):659–67. doi: 10.1513/AnnalsATS.202107-822OC
47. Trittmann JK, Bartenschlag A, Zmuda EJ, Frick J, Stewart WCL, Nelin LD. Using clinical and genetic data to predict pulmonary hypertension in bronchopulmonary dysplasia. Acta Paediatr Int J Paediatr. (2018) 107(12):2158–64. doi: 10.1111/apa.14600
48. Veneroni C, Wallström L, Sindelar R, Dellaca RL. Oscillatory respiratory mechanics on the first day of life improves prediction of respiratory outcomes in extremely preterm newborns. Pediatr Res. (2019) 85(3):312–7. doi: 10.1038/s41390-018-0133-6
49. Verder H, Heiring C, Ramanathan R, Scoutaris N, Verder P, Jessen TE, et al. Bronchopulmonary dysplasia predicted at birth by artificial intelligence. Acta Paediatr Int J Paediatr. (2021) 110(2):503–9. doi: 10.1111/apa.15438
50. Baud O, Lehert P. The beneficial effect of prophylactic hydrocortisone treatment in extremely preterm infants improves upon adjustment of the baseline characteristics. Pediatr Res. (2024) 95(1):251–6. doi: 10.1038/s41390-023-02785-x
51. Jensen EA, Wiener LE, Rysavy MA, Dysart KC, Gantz MG, Eichenwald EC, et al. Assessment of corticosteroid therapy and death or disability according to pretreatment risk of death or bronchopulmonary dysplasia in extremely preterm infants. JAMA Netw Open. (2023) 6(5):e2312277. doi: 10.1001/jamanetworkopen.2023.12277
Keywords: bronchopulmonary dysplasia, prediction, death, mechanical ventilation, preterm
Citation: Sucasas-Alonso A, Pértega-Díaz S, Balboa-Barreiro V, García-Muñoz Rodrigo F and Avila-Alvarez A (2024) Prediction of bronchopulmonary dysplasia in very preterm infants: competitive risk model nomogram. Front. Pediatr. 12:1335891. doi: 10.3389/fped.2024.1335891
Received: 9 November 2023; Accepted: 8 February 2024;
Published: 20 February 2024.
Edited by:
Minesh Khashu, University Hospitals Dorset NHS Foundation Trust, United KingdomReviewed by:
Anne Greenough, King’s College London, United KingdomSimón Lubian Lopez, Hospital Universitario Puerta del Mar, Spain
© 2024 Sucasas-Alonso, Pértega-Díaz, Balboa-Barreiro, García-Muñoz Rodrigo and Avila-Alvarez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sonia Pértega-Díaz s.pertega@udc.es