ORIGINAL RESEARCH article

Front. Pediatr., 27 June 2025

Sec. Pediatric Pulmonology

Volume 13 - 2025 | https://doi.org/10.3389/fped.2025.1566747

Clinical characteristics of bronchopulmonary dysplasia and the risk of sepsis onset prediction via machine learning models

  • 1. The Second School of Clinical Medicine, Southern Medical University, Guangzhou, China

  • 2. Institute of Pediatrics, Faculty of Pediatrics, The Seventh Medical Center of PLA General Hospital, Beijing, China

  • 3. National Engineering Laboratory for Birth Defect Prevention and Control of Key Technology, Beijing, China

  • 4. Beijing Key Laboratory of Pediatric Organ Failure, Beijing, China

  • 5. Department of Pediatric Cardiology, Faculty of Pediatrics, the Seventh Medical Center of PLA General Hospital, Beijing, China

Article metrics

View details

1

Citations

1,6k

Views

746

Downloads

Abstract

Bronchopulmonary dysplasia (BPD), also known as chronic lung disease, is the most common cause of respiratory morbidity in preterm infants. Sepsis plays a significant role in the pathogenesis of BPD, and the systemic inflammatory response caused by sepsis is associated with lung development, leading to simplified alveoli and abnormal vascular development, which are the histological hallmarks of BPD. In this study, we conducted a retrospective analysis of the clinical characteristics of 306 preterm infants with BPD treated at our hospital from December 2019 to December 2022. We subsequently utilized ten machine learning (ML) algorithms and used clinical features to acquire models to predict BPD with sepsis. The performance of the model was evaluated according to the mean area under the receiver operating characteristic curve (AUC), sensitivity, specificity, and accuracy. The mean area under the curve (AUC) of the best predictive model was 0.93. A nomogram for sepsis onset was developed in the primary cohort with four factors: invasive respiratory support, CRIB II score, NEC, and chorioamnionitis. By including clinical features, ML algorithms can predict BPD with sepsis, and the random forest (RF) model (sorted by the mean AUC) performs the best. Our prediction model and nomogram can help clinicians make early diagnoses and formulate better treatment plans for preterm infants with BPD.

Introduction

Bronchopulmonary dysplasia (BPD) is usually caused by mechanical ventilation and long-term use of oxygen (1). BPD, resulting from lung injury that disrupts alveolar and pulmonary vascular development, is one of the most common causes of morbidity and mortality in preterm infants. Owing to the increased survival of extremely low-gestational-age newborns, BPD remains the most common complication associated with prematurity, and its prevalence is increasing (13).

BPD is a multifactorial pathology influenced by a variety of prenatal and postnatal factors that affect mothers and infants (4). BPD pathology worsens and further promotes reactive oxygen species (ROS) production, and subsequent inflammation leads to sepsis (5). If children with BPD develop an infection, this may further exacerbate the inflammatory response, thereby increasing the risk of sepsis (6).

Sepsis is a clinical syndrome involving organ dysfunction caused by a disordered host response to infection (7). Preterm infant sepsis refers to an infection involving the bloodstream in infants aged <28 days of age (8). Sepsis is divided into early-onset sepsis (EOS) and late-onset sepsis (LOS) on the basis of age at presentation after birth. EOS refers to sepsis in neonates occurring within 72 h (h) of birth (some experts use 7 days), and LOS is defined as sepsis occurring at or after 72 h of life (810). EOS is generally caused by an in utero infection or by vertical bacterial transmission from the mother during vaginal delivery, whereas LOS usually occurs not only by vertical bacterial transmission but also by horizontal bacterial transmission from healthcare providers and the environment (11). The clinical manifestations of neonatal sepsis can range from nonspecific or vague symptoms to hemodynamic collapse. The early symptoms may include irritability, lethargy, or poor feeding. Others may quickly progress to respiratory distress, fever, hypothermia, or hypotension, accompanied by poor perfusion and shock (12).

Preterm infants with BPD are more susceptible to developing sepsis because of chronic lung and airway damage. Predicting the risk of BPD with sepsis is crucial for early intervention and improving patient prognosis. Few studies have reported the establishment of predictive models in adults (1315). However, in neonatal medicine, statistical or machine learning (ML) models for predicting patients who may develop BPD with sepsis are relatively rare. Therefore, such models need to be established to help doctors identify risks early, thereby allowing them to take preventive measures or undergo early treatment.

In this study, we conducted a retrospective analysis of the clinical characteristics of preterm infants with BPD at our hospital from December 2019 to December 2022. We subsequently used ML algorithms to identify high-risk factors for the co-occurrence of BPD and sepsis.

Methods

Study design and study participants

This study was approved by the Ethics Committee of the Seventh Medical Center of the PLA General Hospital (S2024-046-01) and was conducted in accordance with the Declaration of Helsinki. The studies were conducted in accordance with local legislation and institutional requirements. Written informed consent for participation was not required from the participants or their legal guardians/next of kin because of the retrospective design of the study.

The study population included extremely preterm infants with BPD who were born or admitted to the Seventh Medical Center of PLA General Hospital (Beijing, China). Trained neonatologists at our center identified the patients. The following criteria were applied to construct the initial dataset:

The inclusion criteria were as follows: very preterm infants (gestational age ≤32 weeks, >28 weeks) and extremely preterm infants (gestational age ≤28 weeks) who met the diagnostic criteria for BPD.

The exclusion criteria were as follows: very preterm infants and extremely preterm infants who died were withdrawn from treatment or discharged within 28 days after birth. Patients with genetic metabolic diseases. Infants without BPD.

The patients were divided into BPD with sepsis and BPD without sepsis groups on the basis of Doppler echocardiography results after at least 36 weeks of corrected gestation.

Clinical feature collection

Clinical data were collected from the electronic medical records of the Seventh Medical Center of the PLA General Hospital, including maternal pregnancy factors [amniotic fluid disorders, hypertension during pregnancy, gestational diabetes mellitus (GDM), preeclampsia, placental abnormality, delivery method, placenta pathology], newborn clinical data [gestational age (GA), birth weight (BWt), sex, multifetal gestations, 1-min Apgar score (5APGAR), 10 min Apgar score (10APGAR), severity of BPD, pulmonary hypertension (PH) (early-PH and BPD-related PH), patent ductus arteriosus (PDA), severe PDA (defined as needing surgery ligation), intraventricular hemorrhage (IVH), necrotizing enterocolitis (NEC), retinopathy of prematurity (ROP), clinical risk index for babies score II (CRIB II), hospital length of stay (day), invasive respiratory support (day), etc.], related laboratory tests (platelet (PLT), C-reactive protein (CRP), white blood cell count (WBC), neutrophil count (N), hemoglobin (HGB)), signs and symptoms of infection episodes such as difficulty breathing, fever (>37.5°C), hypothermia (<36.5°C), abdominal distension, feeding intolerance, etc., and microbiological characteristics. The diagnosis of chorioamnionitis was confirmed via placental histopathology. Echocardiograms at 4–7 days of age were used to assess “early-PH” (16, 17). Other diagnostic criteria were based on related criteria (18, 19).

Statistical methods

To filter for missing data, the missing data module in Python 3.9.12 was used. In Supplementary Figure S1, each column represents a clinical variable, and the white line represents missing data. The denser the lines in each column are, the greater the number of missing values for that variable. Detailed information regarding missing values is provided in Supplementary Table S1. We removed the antenatal corticosteroid (ANC) variable, which was missing in >25% of the observations, to facilitate and ensure study accuracy.

Continuous data are presented as the mean ± standard deviation (SD) or median (interquartile range, Q1, Q3), and intergroup comparisons of normally distributed continuous data were made via two-sample t tests, with the test value being the t value. Intergroup comparisons of nonnormally distributed continuous data via nonparametric tests were performed via the Mann‒Whitney U test, with the test value being the Z value. Categorical data are represented by the number of cases and percentage (%), and intergroup comparisons were made via the rank sum test, with the test value being χ2. To identify the predictors of sepsis, a variance inflation factor (VIF) was first used to test all predictors for multicollinearity, followed by the inclusion of all predictors in the model via multivariate logistic regression analysis. The odds ratios (ORs) for independent risk factors for sepsis were estimated via a stepwise selection method with a 95% confidence interval (CI). A two-tailed P value less than 0.05 indicated a statistically significant difference. Statistical analyses were performed via SPSS version 27.0.

Binary classification was performed using Scikit-learn (version 0.24.1) in Python (version 3.9.12). Ten ML algorithms were used to differentiate between septic BPD patients and nonseptic BPD patients. The following 10 models are used: logistic regression (LR), random forest (RF), support vector machine (SVM), decision tree (DTREE), AdaBoost (ADB), Gaussian naive Bayes (NB), linear discriminant analysis (LDA), k-nearest neighbors (KNN), gradient boosting classifier (GB), and multilayer perceptron (MLP). Considering the limited number of data samples, each classification model for the ML algorithms was built on all the data via the default parameters in the scikit-learn library. The synthetic minority oversampling technique (SMOTE) expands the number of samples in a minority class to ensure equal representation among groups in ML. The bootstrap method was used 1,000 times for internal validation. The performance of each algorithm was assessed on the basis of average sensitivity, specificity, the mean area under the receiver operating characteristic (ROC) curve, and the mean F1 of the resampled samples 1,000 times for pediatric patients with combined septic and nonseptic BPD. Receiver operating characteristic (ROC) curves were plotted via the matplotlib library (version 3.3.4) in Python as part of the internal validation process. The precise contribution (magnitude and direction) of the feature output by each classifier was determined via Shapley additive explanations (SHAPs). The SHAP values were calculated via the RF algorithm for each classifier. SHAP summary plots were visualized in Python via the Sharp library (version 0.39.0).

Nomogram charts (rms package) were drawn using the selected risk factors. The concordance statistic (C statistic) and calibration curve (rms package) were used to distinguish and calibrate the nomograms. Decision curve analysis (DCA) and a clinical impact curve (CIC) were used to evaluate the clinical utility of the model (2022). The DCA curve and CIC (ggDCA package) were used to evaluate the effectiveness and clinical applicability of the risk prediction nomogram. Statistical analyses were performed via version 4.5.0 of the R statistical software.

Results

Demographic and clinical features of sepsis BPD patients and nonsepsis BPD patients

A total of 306 patients with BPD were enrolled in this study, including 177 men (57.8%) and 129 women (42.2%). The demographic and clinical features of the patients are summarized in Table 1.

Table 1

VariablesTotal (N = 306)Total (N = 103)
Sepsis (N = 103)Nonsepsis (N = 203)t/χ2PEOS (N = 47)LOS (N = 56)t2P
Clinical variables
GA(weeks),a mean ± SD26.82 ± 1.5327.63 ± 1.90−3.850<0.00126.85 ± 1.6326.79 ± 1.46−0.0100.992
BWt(g)a, mean ± SD990.65 ± 282.811,104.42 ± 334.13−3.561<0.001997.77 ± 246.79984.34 ± 313.54−0.7150.475
Sex0.1210.7280.1130.737
Male61 (59.2%)116 (57.1%)27 (57.4%)34 (60.7%)
Female42 (40.8%)87 (42.9%)20 (42.6%)22 (39.3%)
Csec7.7930.0050.0820.775
Yes41 (39.8%)114 (56.2%)18 (38.3%)23 (41.1%)
No62 (60.2%)87 (42.8%)29 (61.7%)33 (58.9%)
Singleton0.0320.8580.2010.654
Yes70 (68.0%)140 (69.0%)33 (70.2%)37 (66.1%)
No33 (32.0%)63 (21.0%)14 (29.8%)19 (18.0%)
Stay in the hospital (d)a98.53 ± 41.7584.67 ± 34.05−3.591<0.00197.30 ± 46.0999.57 ± 38.12−0.8680.386
Invasive respiratory support (d)a31.64 ± 38.2028.65 ± 19.63−5.473<0.00133.78 ± 45.4629.89 ± 41.37−0.1560.876
Noninvasive respiratory support (d)a28.65 ± 19.6328.26 ± 23.32−0.6490.51630.13 ± 23.0927.42 ± 17.93−0.1090.913
Nasal catheter oxygen inhalation (d)a14.04 ± 15.6416.79 ± 12.92−2.9570.00312.55 ± 12.3915.29 ± 17.94−0.2880.773
CRIB IIa,c9 (7,11)6 (4,8)−6.560<0.0018 (7,10)9 (8,11)−1.2470.212
1APGARa,c7 (6,8)8 (6,9)−1.9520.0517 (6,8)8 (6,8)−0.0530.958
5APGARa,c9 (8,9)9 (8,10)−2.2200.0269 (8,9)9 (8,9)−1.1020.271
10APGARa,c9 (8,10)9 (9,10)−1.8180.0699 (8,10)9 (8,9.5)−0.4610.654
PIH14 (13.7%)43 (21.2%)2.4840.1156 (13%)8 (14.3%)0.0330.856
GDM23 (22.8%)54 (26.6%)0.5230.47010 (22.2%)13 (23.2%)0.0140.906
Placental abruption19 (18.4%)33 (16.3%)0.2150.6439 (19.1%)10 (17.9%)0.0280.866
Abnormal amniotic fluid20 (19.4%)39 (19.3%)0.0010.9828 (17.0%)12 (21.4%)0.3170.573
Abnormal fetal membranesb2 (1.9%)1 (0.5%)0.2641 (2.1%)1 (1.1%)1
Umbilical cord abnormalities16 (15.7%)23 (11.4%)1.1210.2906 (12.8%)10 (18.2%)0.5620.453
Maternal age(years)32.58 ± 5.1132.11 ± 4.230.8230.41131.30 ± 4.7333.68 ± 5.21−2.2880.024
Outcomes
ROP15.2920.0021.5350.674
Non31 (33.7%)115 (58.4%)15 (34.1%)16 (33.3%)
I12 (13%)16 (8.1%)6 (13.6%)6 (12.5%)
II30 (32.6%)41 (20.8%)12 (27.3%)18 (37.5%)
III19 (20.7%)25 (12.7%)11 (25.0%)8 (9.9%)
Early-PH84 (84.0%)58 (29.1%)6.1760.01338 (45.2%)46 (54.8%)0.0120.913
PDA0.8590.6510.028
Non5 (4.9%)15 (7.5%)1 (2.1%)5 (78.6%)
Without surgical treatment88 (86.3%)169 (84.9%)44 (95.7%)44 (8.9%)
Surgical treatment9 (8.8%)15 (7.5%)2 (4.3%)7 (12.5%)
NEC18 (17.47%)194 (9.56%)14.447<0.0015 (10.6%)13 (23.2%)2.8020.094
WMI6.6760.0360.681
Non72 (69.9%)168 (82.8%)35 (48.6%)37 (51.4%)
I23 (22.3%)26 (12.8%)9 (39.1%)14 (60.9%)
II8 (7.8%)9 (4.4%)3 (37.5%)5 (62.5%)
BPD severity31.646<0.0012.4290.297
Mild12 (11.9%)66 (33.2%)7 (14.9%)5 (9.3%)
Moderate40 (39.6%)94 (47.2%)21 (44.7%)19 (35.2%)
Severe49 (48.5%)39 (19.6%)19 (40.4%)30 (55.6%)
BPD-PH13 (13.1%)10 (5.0%)6.0990.0144 (8.7%)9 (17.0%)1.4820.223
IVH72 (70.6%)121 (61.4%)2.4680.11633 (71.7%)39 (70.6%)0.0530.817
Chorioamnionitisb19 (18.4%)3 (1.5%)<0.00113(27.7%)6(10.7%)4.8770.027

Clinical characteristics of infants with BPD with sepsis vs. those without sepsis and EOS vs. LOS in infants with BPD with sepsis.

Continuous data are presented as the mean ± SD values except for median (interquartile range, Q1, Q3); categorical data as n/N (%); the figures in parentheses are percentages. n: Number of participants. Chi-square, Fisher's exact test, and Mann–Whitney U test were used for discrete and continuous variables.

SD, standard deviation; EOS, early onset sepsis; LOS, late-onset sepsis; BWt, birth weight; GA, gestational age; Csec, cesarean section; 1APGAR, Apgar score at 1 min; 5APGAR, Apgar score at 5 min; 10APGAR, Apgar score at 10 min; CRIB-II, clinical risk index II for infants; PIH, pregnancy-induced hypertension; GDM, gestational diabetes mellitus; IVH, cerebral intraventricular hemorrhage; ROP, retinopathy of prematurity; PDA, patent ductus arteriosus; BPD, bronchopulmonary dysplasia; NEC, necrotizing enterocolitis; WMI, brain white matter injury; BPD-PH, BPD-associated pulmonary hypertension.

a

Mann–Whitney U test, Z value.

b

Fisher's exact test.

c

Data are shown as median (interquartile range, Q1, Q3).

Our analysis revealed that the lower the birth weight was, the greater the chance of developing sepsis (Figure 1A). The lower the gestational age was, the greater the incidence of sepsis, especially in preterm infants born at ≤27 weeks (Figure 1B).

Figure 1

The features of the patients with BPD are described in Table 1. A comparison of the two groups revealed that several characteristics, such as GA, BWt, delivery mode, hospital stay, duration of invasive ventilator ventilation, oxygen intrathecal time, CRIB II score and 5APGAR, were significantly different (p < 0.05). These results suggest that the earlier the gestational age and the lower the birth weight, the greater the proportion of infants with sepsis, the longer the hospital stay, and the longer the duration of invasive mechanical ventilation. Sex, singleton pregnancy, noninvasive ventilator ventilation time, and 1APGAR and 10APGAR levels were not significantly different between the two groups.

The analysis of the clinical characteristics of maternal and postnatal complications, including chorioamnionitis, varying degrees of ROP, NEC, early postnatal pulmonary hypertension (early PH), severity grading of BPD, and BPD-PH (p < 0.05), revealed significant differences between sepsis patients and nonsepsis BPD patients (Table 1).

Machine learning algorithms and the development and evaluation of a nomogram for BPD patients with sepsis

In previous clinical practices, BWt, GA, invasive respiratory support, CRIB II, 5APGAR, NEC, early PH, and chorioamnionitis were found to be independent risk factors for sepsis in preterm infants with BPD. Furthermore, collinearity diagnostic analysis demonstrated that the VIFs of these risk factors, except BWt and GA, were less than 4, indicating that there was no strong indication of multicollinearity among the variables. Considering that data were missing, we included a sample size of 284 cases for model construction after removing missing data. These four variables, invasive respiratory support, CRIB II, NEC, and chorioamnionitis, were incorporated into the final predictive model on the basis of stepwise regression results. Multivariate analysis revealed that invasive respiratory support (OR, 1.02; 95% CI, 1.00–1.04; p < 0.05), CRIB II (OR, 1.28; 95% CI, 1.12–1.46; p < 0.05), NEC (OR, 3.10; 95% CI, 1.23–7.82; p < 0.05), and chorioamnionitis (OR, 10.40; 95% CI, 2.85–38.02; p < 0.05) independently increased the risk for the development of sepsis in BPD infants (Table 2).

Table 2

VariablesOR (95% CI)P value
Invasive respiratory support1.019 (1.003–1.036)0.022
NEC3.098 (1.228–7.815)0.017
Chorioamnionitis10.400 (2.845–38.015)<0.001
CRIB II1.280 (1.121–1.460)<0.001

Binary logistic regression analysis of sepsis and nonsepsis patients with BPD.

NEC, necrotizing enterocolitis; CRIB-II, clinical risk index II for babies. P values were calculated via a stepwise multivariate regression model.

We constructed 10 models of ML by comparing their model performance. In the predictive model built with the above four factors, the mean area under the ROC curve of each model reached or approached 0.8. Among them, the RF model emerged as the most effective predictor, with a mean area under the receiver operating characteristic (ROC) curve of 0.93 (Figure 2A). The mean F1 of the 4/10 models reached 0.8, with the RF model having the best predictive performance, with a mean F1 of 0.87 (Table 3).

Figure 2

Table 3

ModelsAUC Mean, 95% CIF1 Mean, 95% CISensitivity Mean, sdSpecificity Mean, sd
LR0.790 (0.784 0.794)0.679 (0.650 0.706)0.593, 0.0260.846, 0.017
RF0.934 (0.913 0.952)0.877 (0.851 0.902)0.875, 0.0230.879, 0.024
SVM0.794 (0.772 0.804)0.692 (0.658 0.728)0.622, 0.0410.827, 0.033
DTREE0.891 (0.858 0.919)0.865 (0.834 0.894)0.851, 0.0250.885, 0.027
ADB0.848 (0.824 0.867)0.760 (0.718 0.798)0.725, 0.0440.818, 0.042
NB0.778 (0.771 0.785)0.619 (0.484 0.643)0.489, 0.0400.911, 0.016
LDA0.791 (0.788 0.794)0.666 (0.641 0.694)0.568, 0.0260.864, 0.021
KNN0.878 (0.854 0.901)0.810 (0.778 0.839)0.821, 0.0290.794, 0.031
GB0.905 (0.885 0.925)0.840 (0.808 0.867)0.815, 0.0300.874, 0.030
MLP0.801 (0.790 0.810)0.699 (0.663 0.730)0.636, 0.0360.818, 0.027

Predictive performance comparison of the ten types of machine learning algorithms.

LR, logistic regression; RF, random forest; SVM, support vector machine; DTREE, decision tree; ADB, AdaBoost; NB, Gaussian naive Bayes; LDA, linear discriminant analysis; KNN, k-nearest neighbors; GB, gradient boosting classifier; MLP, multilayer perceptron. AUC, area under the curve; SD, standard deviation.

The RF model was used to demonstrate the importance of various features in the model with the best predictive performance. The CRIB II score ranked first, followed by the duration of invasive mechanical ventilation and NEC. These findings indicated that a higher CRIB II score and longer duration of invasive mechanical ventilation were associated with greater risk (Figure 2B). In addition, the coexistence of chorioamnionitis and NEC is a major risk factor for sepsis. These findings suggest that, in the absence of chorioamnionitis in pregnant mothers before birth, actively preventing prenatal infections and NEC after birth can prevent the occurrence of sepsis.

We constructed a nomogram based on the logistic regression model, which allowed for a more intuitive prediction of the risk of sepsis (Figure 2C). A nomogram for sepsis onset was developed in the primary cohort with four factors: invasive respiratory support, CRIB II score, NEC, and chorioamnionitis. These factors were screened via logistic regression analysis. The C-statistic of the nomogram was 0.79. The Brier score was 0.164, which is smaller than 0.25.

The calibration curve of the constructed regression model generated via the bootstrap method (with 1,000 repeated samples) indicated an average absolute error of 0.03, suggesting that the predicted risk of sepsis was quite accurate after calibration and that the model did not overfit (Figure 2D).

According to the DCA curve of the prediction model and the verified DCA curve (Figure 2E), the net benefit corresponding to the curve was above 0 over a wide range of decision thresholds (0.2–0.8). It is far from the two extreme curves of “None” and “All”. The CIC curves of the prediction model validation (Figure 2F) show that the results predicted by the model are close to the actual results in a wide range of risk thresholds (0.5–1.0), which indicates practical application value in a wider range of clinical situations.

Clinical characteristics of BPD patients with EOS and LOS

Next, we analyzed the clinical information regarding birth details, postnatal complications, and symptoms from the BPD complicated with EOS group and BPD complicated with LOS group and found significant differences in the LOS group in terms of biochemical indicators (platelet count, peak values of CRP), symptoms during infection (fever, feeding intolerance, and abdominal distension), and postnatal complications (PDA) (p < 0.05). These findings suggest that infants with BPD and LOS are more likely to have abnormalities in the above indicators and clinical symptoms. Other factors, such as GA, BWt, the Apgar score, sex, the CRIB II score, the total SOFA score, the WBC count, HGB, postnatal complications such as NEC and early postnatal PH, and symptoms during infection, such as oxygen desaturation and bradycardia, were not significantly different between the two groups. However, in both groups of patients with sepsis, the percentage of patients with oxygen desaturation reached over 95%, whereas respiratory arrest and bradycardia were more common in patients with LOS, indicating that these patients require close attention from physicians (Table 1).

The analysis of biochemical indicators from the groups with BPD complicated by EOS and LOS revealed that the peak values of CRP were relatively high in the group with BPD complicated by LOS. In contrast, in the group with BPD complicated by EOS, the proportion of infants with abnormal white blood cell counts was relatively high (Table 4).

Table 4

VariablesEOS (N = 47)LOS (N = 56)t/χ2P
Laboratory variables
nSOFAa1 (0,2)2 (1,4)−1.8370.066
WBCb (109/L)16.85 ± 12.2914.55 ± 7.74−0.3310.741
Nb (109/L)7.13 ± 6.1410.81 ± 20.03−0.1560.876
HGB (g/L)131.08 ± 26.74122.13 ± 31.121.3970.166
PLTb (109/L)200.18 ± 158.55130.43 ± 95.60−2.4120.016
Peak episode CRP, mg/L(median ± SD)b23.71 ± 37.4441.00 ± 37.22−2.9020.004
Clinical signs associated episode
Decreased oxygen saturation46 (97.9%)56 (100%)1.2030.273
Fever or hypothermia5 (10.6%)15 (26.8%)4.2580.039
Abdominal distension/increased abdominal girth8 (17.0%)29 (51.8%)13.417<0.001
Feeding intolerance1 (7.7%)12 (21.4%)8.6320.003
Increased bradycardic episodes8 (17%)15 (26.8%)1.4050.236
Apnea12 (25.5%)31 (55.4%)9.3470.002

Laboratory characteristics and clinical signs of EOS vs. LOS in infants with BPD with sepsis.

Continuous data are presented as the mean ± SD values except for median (interquartile range, Q1, Q3); categorical data as n/N (%); figures in parentheses are percentages.

nSOFA, sequential organ failure assessment score; WBC, white blood cell; N, neutrophil count; HGB, hemoglobin; PLT, platelet; CRP, C-reactive protein; SD, standard deviation.

a

The data was shown as median (interquartile range, Q1, Q3).

b

Mann–Whitney U test, Z value.

The etiological characteristic analysis of blood cultures from the groups with BPD complicated by EOS and BPD complicated by LOS revealed that in the group with BPD complicated by LOS, gram-positive bacteria accounted for 37.8%, with Staphylococcus epidermidis and Enterococcus faecalis being the main bacteria, each accounting for 16.2%, followed by Streptococcus hemolyticus and Staphylococcus aureus. Clinically suspected fungal infections accounted for 23.4%, and Candida spp. accounted for 3.6%. A greater proportion of infants with LOS had coinfections with more than one pathogen. In contrast, the clinically suspected rate was higher in patients with EOS (64.4%), indicating that the positive rate of pathogen detection is lower in patients with EOS and that fungal infections are rare (Table 5).

Table 5

Clinical diagnosisEOS(N = 47)LOS (N = 56)
29 (64.4%)7 (12.7%)
Organisms from blood culture
Gram-nagative bacteria(total)4 (8.7%)16 (28.8%)
Klebsiella pnenmoniae1 (2.1%)8 (14.4%)
Serratia marcescens-3 (5.4%)
Enterobacter cloacae2 (4.2%)1 (1.8%)
Enterobacter aerogenes-3 (5.4%)
Escherichia coli1 (2.1%)1 (1.8%)
Gram-positive bacteria(total)10 (21%)25 (45.0%)
Staphylococcus epidermidis5 (10.6%)9 (16.2%)
Staphylococcus aureus-4 (7.2%)
Enterococcus faecalis2 (4.2%)1 (1.8%)
Enterococcus faecalis group D)3 (6.3%)9 (16.2%)
Hemolytic staphylococcus-1 (1.8%)
Staphylococcus human subspecies-1 (1.8%)
Fungal infection(total)2 (4.2%)12 (21.6%)
Pseudohyphae-2 (3.6%)
Candida guilliermondii1 (1.8%)
Probable fungal infection2 (4.2%)9 (16.2%)
Other organisms3 (6.3%)3 (5.4%)
2 organisms1 (2.1%)3 (5.4%)
3 organisms-2(3.6%)

Microbiological characteristics of infants and manifestations of EOS and LOS episodes [results are reported as N(%) unless otherwise specified].

Discussion

BPD is a common complication in premature infants and a chronic lung disease in neonates. A recent study revealed that in the United States, the incidence of chronic lung disease in extremely preterm infants born between 24 and 28 weeks of gestation has been increasing since 2012 (23). In the past decade, the survival rate of extremely preterm infants in China has improved significantly. However, the incidence of BPD has not decreased significantly, reaching 40.7% (24). In addition to chronic lung disease, the incidence of long-term complications in BPD patients, such as those involving the cardiovascular, nervous, digestive, and endocrine systems, has also increased (24).

Neonatal sepsis, especially early-onset sepsis (EOS), often presents with subtle and nonspecific symptoms, making prompt diagnosis challenging. ML models can identify high-risk neonates before symptoms become apparent. This early detection enables timely intervention, which is crucial for reducing morbidity and mortality associated with sepsis (25). ML models can integrate various clinical data points, such as vital signs and laboratory indicators, to identify risk factors more accurately. This approach significantly improves diagnostic sensitivity and overall accuracy. In addition, ML models can help allocate medical resources more efficiently by rapidly screening high-risk neonates, ensuring that critical interventions are directed toward those who need them most, enhancing the overall efficiency of healthcare delivery (26).

Inflammation/infection, including chorioamnionitis in utero and postnatal systemic infectious inflammation, has been shown to be a risk factor for neonatal sepsis (27, 28). Sepsis-induced systemic inflammation is a cause of neonatal mortality (29, 30).

In this study, we conducted a retrospective analysis of the clinical characteristics of preterm infants with BPD from December 2019 to December 2021 and used ML methods to rank the importance of various features and fit a predictive model. All risk factors were ranked in terms of their importance in deciding whether to include or reject them. The most important factor is the CRIB II score for neonates, followed by invasive mechanical ventilation, intrauterine infection (chorioamnionitis), and the subsequent development of acute NEC. This finding is consistent with the results of a previous questionnaire survey study (31). Early chorioamnionitis was suspected to be the main cause of EOS triggered by intrauterine infection, whereas late-onset NEC in later stages of enteral feeding was the main cause of LOS. In terms of the predictive performance of the ten models adopted, the mean area under the ROC curve of all the models can approach or reach 0.8. The RF model showed the best predictive performance, with a value of 0.93. Its clinical and practical value has been determined by DCA and CIC curves. Moreover, the nomogram constructed from these findings can intuitively predict the risk of sepsis, and when combined with the results of etiological analysis, it is beneficial to further guide clinical doctors in the diagnosis, treatment, and rational use of antimicrobial drugs. In the etiological analysis of blood for EOS and LOS in preterm infants with BPD in our hospital, gram-positive bacteria were commonly observed in LOS, and clinically suspected fungal infections, such as Candida albicans and Candida parapsilosis, were frequently observed. In this study, the changes in the platelet count and CRP level were the most significant for LOS, whereas the change in the WBC count was more pronounced for the early detection of sepsis. The platelet count has good sensitivity and specificity for diagnosing neonatal sepsis and can be used as a diagnostic tool for neonatal sepsis (32, 33). In addition, the PLT can serve as a diagnostic indicator of late-onset neonatal pneumonia (34). However, considering that CRP is nonspecific for diagnosing sepsis and has poor sensitivity, a single standard for diagnosing sepsis must be more accurate. In clinical practice, a comprehensive assessment is needed, which should be combined with clinical symptoms, blood cell counts, and other indicators.

The advent of big data and the artificial intelligence era has also driven medical progress. In recent years, the combination of ML with medical data has not only aided scientific research output but has also been continuously applied in clinical work. Statistical models are the best way to analyze and predict outcomes (35). A recent study revealed that ML models constructed from the vital signs of newborns within 24 h before infection can also predict sepsis quite well, with the area under the ROC curve reaching 0.82 (36).

In summary, our study indicated that sepsis is a risk factor for BPD. The ML algorithm suggests that the CRIB II score, duration of invasive mechanical ventilation, incidence of chorioamnionitis, and incidence of neonatal necrotizing enterocolitis are high-risk factors for the co-occurrence of sepsis. By combining the nomogram and characteristics of etiology, the risk of sepsis can be calculated, which may further reduce the exposure to and duration of antibiotic use in preterm infants and has a certain guiding significance for clinical diagnosis and treatment.

Limitations

Our study was a retrospective case analysis study from a single center with a relatively small number of case samples. The etiological characteristics of the samples included only peripheral blood and did not include samples from other sources, such as oropharyngeal or tracheobronchial secretions or sputum. Further studies are needed to expand the sample size and to conduct prospective multicenter cohort studies. Additionally, in this study, the subjects were preterm infants, making it difficult to collect sufficient blood samples for the detection of inflammation-related indicators. In the future, given the availability of technologies that can measure these indicators, we will also incorporate these indicators into our analysis. Although ML models have shown promise in predicting sepsis, their generalizability and clinical interpretability still require further research and validation. Moreover, to increase predictive accuracy, future studies may need to incorporate more patient data and vital sign indicators as well as further refine the algorithms and structure of the models.

Statements

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 This study was approved by the Ethics Committee of the Seventh Medical Center of the PLA General Hospital (S2024-046-01) and was conducted in accordance with the Declaration of Helsinki. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because of the retrospective design of the study.

Author contributions

YaW: Data curation, Formal analysis, Investigation, Writing – original draft. YiW: Investigation, Project administration, Software, Validation, Writing – original draft. LS: Conceptualization, Formal analysis, Methodology, Resources, Writing – original draft. JL: Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft. YX: Formal analysis, Methodology, Project administration, Software, Writing – original draft. LY: Formal analysis, Investigation, Methodology, Project administration, Writing – original draft. SH: Conceptualization, Investigation, Methodology, Project administration, Validation, Writing – original draft, Writing – review & editing. FZ: Conceptualization, Formal analysis, Methodology, Resources, Supervision, Validation, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the National Natural Science Foundation of China (grant number: 82341095). and the Beijing E-Town Cooperation and Development Foundation (grant number: YCXJ-JZ-2023-017).

Acknowledgments

We wish to thank all the parents and infants who supported our research and pushed us to do our best.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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.2025.1566747/full#supplementary-material

References

  • 1.

    ThébaudBGossKNLaughonMWhitsettJAAbmanSHSteinhornRHet alBronchopulmonary dysplasia. Nat Rev Dis Primers. (2019) 5:78. 10.1038/s41572-019-0127-7

  • 2.

    StollBJHansenNIBellEFShankaranSLaptookARWalshMCet alNeonatal outcomes of extremely preterm infants from the NICHD Neonatal Research Network. Pediatrics. (2010) 126:44356. 10.1542/peds.2009-2959

  • 3.

    WadhawanRVohrBRFanaroffAAPerrittRLDuaraSStollBJet alDoes labor influence neonatal and neurodevelopmental outcomes of extremely-low-birth-weight infants who are born by cesarean delivery?Am J Obstet Gynecol. (2003) 189:5016. 10.1067/S0002-9378(03)00360-0

  • 4.

    BalanyJBhandariV. Understanding the impact of infection, inflammation, and their persistence in the pathogenesis of bronchopulmonary dysplasia. Front Med (Lausanne). (2015) 2:90. 10.3389/fmed.2015.00090

  • 5.

    ChoiYRekersLDongYHolzfurtnerLGoetzMJShahzadTet alOxygen toxicity to the immature lung-part I: pathomechanistic understanding and preclinical perspectives. Int J Mol Sci. (2021) 22:11006. 10.3390/ijms222011006

  • 6.

    BaudOMauryLLebailFRamfulDMoussawiFENicaiseCet alEffect of early low-dose hydrocortisone on survival without bronchopulmonary dysplasia in extremely preterm infants (PREMILOC): a double-blind, placebo-controlled, multicentre, randomised trial. Lancet. (2016) 387:182736. 10.1016/S0140-6736(16)00202-6

  • 7.

    van der PollTShankar-HariMWiersingaWJ. The immunology of sepsis. Immunity. (2021) 54:245064. 10.1016/j.immuni.2021.10.012

  • 8.

    SinghMAlsaleemMGrayCP. Neonatal Sepsis. Petersburg: StatPearls Publishing (2024).

  • 9.

    HornikCPFortPClarkRHWattKBenjaminDKJr.SmithPBet alEarly and late onset sepsis in very-low-birth-weight infants from a large group of neonatal intensive care units. Early Hum Dev. (2012) 88(Suppl 2):S6974. 10.1016/S0378-3782(12)70019-1

  • 10.

    StollBJGordonTKoronesSBShankaranSTysonJEBauerCRet alLate-onset sepsis in very low birth weight neonates: a report from the national institute of child health and human development neonatal research network. J Pediatr. (1996) 129:6371. 10.1016/S0022-3476(96)70191-9

  • 11.

    SongWJungSYBaekHChoiCWJungYHYooS. A predictive model based on machine learning for the early detection of late-onset neonatal sepsis: development and observational study. JMIR Med Inform. (2020) 8:e15965. 10.2196/15965

  • 12.

    CelikIHHannaMCanpolatFEMohanP. Diagnosis of neonatal sepsis: the past, present and future. Pediatr Res. (2022) 91:33750. 10.1038/s41390-021-01696-z

  • 13.

    YangZCuiXSongZ. Predicting sepsis onset in ICU using machine learning models: a systematic review and meta-analysis. BMC Infect Dis. (2023) 23:635. 10.1186/s12879-023-08614-0

  • 14.

    FleurenLMKlauschTLTZwagerCLSchoonmadeLJGuoTRoggeveenLFet alMachine learning for the prediction of sepsis: a systematic review and meta-analysis of diagnostic test accuracy. Intensive Care Med. (2020) 46:383400. 10.1007/s00134-019-05872-y

  • 15.

    WangDLiJSunYDingXZhangXLiuSet alA machine learning model for accurate prediction of sepsis in ICU patients. Front Public Health. (2021) 9:754348. 10.3389/fpubh.2021.754348

  • 16.

    MouraniPMSontagMKYounoszaiAMillerJIKinsellaJPBakerCDet alEarly pulmonary vascular disease in preterm infants at risk for bronchopulmonary dysplasia. Am J Respir Crit Care Med. (2015) 191:8795. 10.1164/rccm.201409-1594OC

  • 17.

    KimHHSungSIYangMSHanYSKimHSAhnSYet alEarly pulmonary hypertension is a risk factor for bronchopulmonary dysplasia-associated late pulmonary hypertension in extremely preterm infants. Sci Rep. (2021) 11:11206. 10.1038/s41598-021-90769-4

  • 18.

    XieXKongBDuanT. Obstetrics and Gynecology. 9th Edition. Beijing: People’s Medical Publishing House (2018).

  • 19.

    ShaoXYeHQiuX. Practice of Neonatology. 5th ed. Beijing: People’s Medical Press (2019).

  • 20.

    VickersAJElkinEB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. (2006) 26:56574. 10.1177/0272989X06295361

  • 21.

    RoussonVZumbrunnT. Decision curve analysis revisited: overall net benefit, relationships to ROC curve analysis, and application to case-control studies. BMC Med Inform Decis Mak. (2011) 11:45. 10.1186/1472-6947-11-45

  • 22.

    KerrKFBrownMDZhuKJanesH. Assessing the clinical impact of risk prediction models with decision curves: guidance for correct interpretation and appropriate use. J Clin Oncol. (2016) 34:253440. 10.1200/JCO.2015.65.5654

  • 23.

    HorbarJDGreenbergLTBuzasJSEhretDEYSollRFEdwardsEM. Trends in mortality and morbidities for infants born 24 to 28 weeks in the US: 1997–2021. Pediatrics. (2024) 153:e2023064153. 10.1542/peds.2023-064153

  • 24.

    ZhuZHeYYuanLChenLYuYLiuLet alTrends in bronchopulmonary dysplasia and respiratory support among extremely preterm infants in China over a decade. Pediatr Pulmonol. (2024) 59:399407. 10.1002/ppul.26761

  • 25.

    LeSHoffmanJBartonCFitzgeraldJCAllenAPellegriniEet alPediatric severe sepsis prediction using machine learning. Front Pediatr. (2019) 7:413. 10.3389/fped.2019.00413

  • 26.

    SahuPRaj StanlyEASimon LewisLEPrabhuKRaoMKunhikattaV. Prediction modelling in the early detection of neonatal sepsis. World J Pediatr. (2022) 18:16075. 10.1007/s12519-021-00505-1

  • 27.

    ThomasWSpeerCP. Chorioamnionitis: important risk factor or innocent bystander for neonatal outcome?Neonatology. (2011) 99:17787. 10.1159/000320170

  • 28.

    AfonsoESmetsKDeschepperMVerstraeteEBlotS. The effect of late-onset sepsis on mortality across different gestational ages in a neonatal intensive care unit: a historical study. Intensive Crit Care Nurs. (2023) 77:103421. 10.1016/j.iccn.2023.103421

  • 29.

    SalimiUDummulaKTuckerMHDela CruzCSSampathV. Postnatal sepsis and bronchopulmonary dysplasia in premature infants: mechanistic insights into “New BPD”. Am J Respir Cell Mol Biol. (2022) 66:13745. 10.1165/rcmb.2021-0353PS

  • 30.

    SpasojevićIObradovićBSpasićS. Bench-to-bedside review: neonatal sepsis-redox processes in pathogenesis. Critical Care. (2012) 16:221. 10.1186/cc11183

  • 31.

    MoftianNSamad SoltaniTMirniaKEsfandiariATabibMSRezaei HachesuP. Clinical risk factors for early-onset sepsis in neonates: an international Delphi study. Iran J Med Sci. (2023) 48:5769. 10.30476/IJMS.2022.92284.2352

  • 32.

    WorkuMAynalemMBisetSWolduBAdaneTTigabuA. Role of complete blood cell count parameters in the diagnosis of neonatal sepsis. BMC Pediatr. (2022) 22:411. 10.1186/s12887-022-03471-3

  • 33.

    YochpazSFriedmanNZirkinSBlumovichAMandelDMaromR. C-reactive protein in early-onset neonatal sepsis—a cutoff point for CRP value as a predictor of early-onset neonatal sepsis in term and late preterm infants early after birth?J Matern Fetal Neonatal Med. (2022) 35:45527. 10.1080/14767058.2020.1856068

  • 34.

    MetwaliWAElmashadAMHazzaaSMEAl-BeltagiMHamzaMB. Salivary C-reactive protein and mean platelet volume as possible diagnostic markers for late-onset neonatal pneumonia. World J Clin Pediatr. (2024) 13:88645. 10.5409/wjcp.v13.i1.88645

  • 35.

    MayM. Eight ways machine learning is assisting medicine. Nat Med. (2021) 27:23. 10.1038/s41591-020-01197-2

  • 36.

    HonoréAForsbergDAdolphsonKChatterjeeSJostKHerleniusE. Vital sign-based detection of sepsis in neonates using machine learning. Acta Paediatr. (2023) 112:68696. 10.1111/apa.16660

Summary

Keywords

bronchopulmonary dysplasia, sepsis, machine learning algorithms, nomogram, prediction model

Citation

Wang Y, Wang Y, Song L, Li J, Xie Y, Yan L, Hu S and Feng Z (2025) Clinical characteristics of bronchopulmonary dysplasia and the risk of sepsis onset prediction via machine learning models. Front. Pediatr. 13:1566747. doi: 10.3389/fped.2025.1566747

Received

25 January 2025

Accepted

10 June 2025

Published

27 June 2025

Volume

13 - 2025

Edited by

Christopher G. Wilson, Loma Linda University, United States

Reviewed by

Samuel Franklin Feng, Paris-Sorbonne University Abu Dhabi, United Arab Emirates

Vladimir Pohanka, Retired, Liptovská Teolička, Slovakia

Updates

Copyright

*Correspondence: Siqi Hu Zhichun Feng

†These authors have contributed equally to this work and share first authorship

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics