Ibuprofen for hemodynamically significant patent ductus arteriosus in extremely preterm infants: a target trial emulation
Highlight box
Key findings
• In this multicenter prospective cohort emulating a target trial, ibuprofen treatment for hemodynamically significant patent ductus arteriosus (hsPDA) showed heterogeneous effects, with an association with lower mortality observed among clinically intolerant, high-risk extremely preterm infants.
What is known and what is new?
• Randomized trials have not demonstrated improved outcomes with pharmacological treatment for patent ductus arteriosus (PDA) in preterm infants. Persistent controversy in PDA management has contributed to a lack of clinical equipoise and challenges in patient recruitment, potentially leading to the underrepresentation of high-risk infants in randomized trials.
• By stratifying infants according to clinical tolerance to PDA-related hemodynamic changes, this study demonstrated heterogeneity in the association between ibuprofen treatment and mortality across clinical tolerance subgroups.
What is the implication, and what should change now?
• Risk stratification based on clinical tolerance may help inform decision-making for PDA management.
• Future randomized trials may consider specifically targeting this clinically intolerant phenotype to further evaluate these findings.
Introduction
Patent ductus arteriosus (PDA) is one of the most common cardiovascular conditions in extremely preterm infants. A hemodynamically significant PDA (hsPDA) can lead to pulmonary overcirculation and systemic hypoperfusion and has been associated with an increased risk of prematurity-related morbidities and mortality (1-3). The optimal management of PDA in preterm infants—particularly whether early pharmacological treatment should be initiated—has challenged neonatologists for more than five decades and remains one of the most controversial issues in neonatology (4,5).
Numerous randomized controlled trials (RCTs) have evaluated pharmacological treatment for PDA but have failed to demonstrate that ductal closure improves clinical outcomes in preterm infants (6-13). Although recent RCTs have adopted increasingly rigorous designs, marked and strongly held differences in clinical opinion regarding PDA treatment persist, resulting in a lack of sufficient clinical equipoise and substantial challenges in patient recruitment. Consequently, a considerable proportion of eligible infants—particularly those at higher risk—have not been enrolled in randomization (6,8,14). As a result, infants enrolled in prior RCTs may have disproportionately represented those for whom PDA had a relatively limited impact on clinical outcomes, such that ductal closure was unlikely to confer a measurable benefit (14).
Clinically, extremely preterm infants exhibit interindividual variability in their tolerance to PDA-related hemodynamic changes. While most infants are able to tolerate the hemodynamic effects of PDA, a subset develop cardiopulmonary instability due to inadequate compensation and require higher levels of respiratory or hemodynamic support. These clinically intolerant, high-risk infants may have been underrepresented in prior RCTs, and whether they benefit from pharmacological treatment remains unclear.
Accordingly, we used prospective cohort data from the Chinese Multicenter Collaboration Platform for PDA in Extremely Preterm Infants (CMCP-PDA) to emulate a target trial evaluating the association between early ibuprofen treatment for hsPDA and clinical outcomes in extremely preterm infants. We further conducted subgroup analyses based on infants’ clinical tolerance to PDA-related hemodynamic changes, with the aim of identifying subpopulations that may derive benefit from pharmacological treatment and of informing targeted therapeutic strategies and the design of future randomized trials. We present this article in accordance with the STROBE reporting checklist (available at https://tp.amegroups.com/article/view/10.21037/tp-2026-1-0017/rc).
Methods
Data source
This study employed a target trial emulation design based on a multicenter prospective cohort. Data were derived from the CMCP-PDA cohort, which has been registered with the Chinese Clinical Trial Registry (ChiCTR2400091371). The CMCP-PDA cohort was established to investigate the epidemiology and management of PDA in extremely preterm infants in China, to address ongoing clinical controversies regarding the identification of infants who may benefit from PDA treatment, and to promote collaborative quality improvement and multicenter clinical research. The CMCP-PDA cohort was initiated on November 1, 2024, and prospectively enrolled all preterm infants with a gestational age ≤29+6 weeks who were admitted to participating centers within 3 days after birth. Clinical data were prospectively collected using the REDCap electronic data capture system. The database comprehensively captured neonatal demographic characteristics, maternal perinatal factors, laboratory findings, therapeutic interventions, and major neonatal morbidities. Echocardiographic assessments were systematically performed on postnatal days 3–5 to evaluate PDA-related hemodynamic characteristics, with concurrent recording of the level of respiratory and hemodynamic support at the time of assessment. A total of 13 tertiary neonatal intensive care units participated in CMCP-PDA, including 3 pediatric specialty hospitals and 10 maternal and child health hospitals, all of which are large tertiary centers with established capacity for the care of extremely preterm infants in China. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study protocol was approved by the Institutional Ethics Committee of the Children’s Hospital of Nanjing Medical University (No. 202406018-1), and informed consent was obtained from all individual participants’ legal guardians.
Quality control
CMCP-PDA employed a standardized data entry manual, and all data abstractors at participating centers received standardized training. The electronic data capture system incorporated automated input validation and error-checking functions to ensure standardized and consistent data collection across sites. The CMCP-PDA coordinating center (Children’s Hospital of Nanjing Medical University) performed regular audits of data accuracy and completeness across participating centers and provided timely feedback to maintain data quality. Additionally, standardized echocardiography protocols were developed based on the 2024 American Society of Echocardiography guidelines for targeted neonatal echocardiography in the neonatal intensive care unit (15) and Practical Neonatal Echocardiography, edited by Siassi (16). These protocols included detailed echocardiographic operating procedures and standardized forms for PDA assessment and measurement (Table S1). Echocardiographic examinations were performed as targeted neonatal echocardiography, primarily by physicians specialized in neonatal cardiac ultrasonography and, in some centers, by trained neonatologists. In all cases, image acquisition and interpretation were performed by the same physician at each center. All participating sonographers underwent dedicated training to ensure consistency and adherence to standardized echocardiographic techniques across centers.
Patients
All preterm infants with a gestational age ≤29+6 weeks admitted to CMCP-PDA-participating centers within 3 days after birth between November 1, 2024, and September 30, 2025 were assessed for eligibility. Infants were eligible for inclusion if they underwent echocardiographic assessment on postnatal days 3–5 and were diagnosed with hsPDA. Infants were excluded if they were transferred to non-CMCP-PDA hospitals, received pharmacological treatment for PDA prior to echocardiographic assessment, died on the day of echocardiography, or had contraindications to pharmacological PDA therapy. Contraindications included a platelet count <50×109/L on admission, intraventricular hemorrhage (IVH) grade ≥3, necrotizing enterocolitis (NEC) stage ≥ II, or spontaneous intestinal perforation occurring before or on the day of echocardiography. hsPDA was defined according to Clyman’s echocardiographic criteria as a ductal diameter ≥1.5 mm with left-to-right shunting in combination with one or more of the following: ductal flow velocity ≤2.5 m/s; a left atrial-to-aortic root diameter ratio ≥1.5; or absent or reversed diastolic flow in the descending aorta (17). The study sample size was determined by accrual within the CMCP-PDA cohort.
Target trial emulation
Time zero
Time zero was defined as the day of echocardiographic assessment on postnatal days 3–5. To emulate the eligibility and exclusion process of a RCT, infants were required to meet all inclusion criteria at time zero. Infants who had received pharmacological treatment for PDA prior to time zero, died at time zero, or had contraindications to pharmacological treatment were excluded.
Grace period
Given that pharmacological treatment for PDA is not routinely initiated immediately after echocardiographic assessment in clinical practice, a 2-day grace period was specified.
Treatment strategies and assignment
Across all participating CMCP-PDA centers, initial pharmacological treatment for PDA consisted of oral ibuprofen, administered at a dose of 10 mg/kg on the first day, followed by 5 mg/kg on each of the subsequent two days. Other pharmacologic agents for PDA, including acetaminophen/paracetamol or indomethacin, were not used as initial treatment during the study period. Infants were classified into the treatment group or the non-treatment group according to whether ibuprofen was initiated within the grace period. Treatment assignment was defined based on the grace period and was not altered by any subsequent pharmacologic treatments. Inverse probability of treatment weighting (IPTW) was applied to balance baseline characteristics between groups and emulate random treatment assignment.
Follow-up and outcomes
All infants were followed until death or discharge from the neonatal intensive care unit. Outcomes included mortality, IVH grade ≥3, NEC stage ≥ II, bronchopulmonary dysplasia (BPD) grade ≥2, death or BPD grade ≥2, and a composite outcome of death, IVH grade ≥3, NEC stage ≥II, or BPD grade ≥2.
Target trial emulation design and participant assignment
A schematic illustration of the target trial design and participant assignment pathways is presented in Figure 1.
Subgroup analysis
Clinical tolerance to hsPDA was prespecified as a potential effect modifier, based on the physiological concept that the clinical consequences of ductal shunting depend not only on echocardiographic characteristics, but also on the infant’s ability to tolerate the resulting cardiopulmonary load. Infants requiring escalation of respiratory or hemodynamic support are more likely to experience clinically significant PDA-related compromise. Therefore, subgroup analyses were performed according to tolerance status (intolerant vs tolerant) to evaluate whether treatment effects differed by clinical tolerance. The intolerant subgroup was defined as infants unable to tolerate the hemodynamic effects of hsPDA at the time of echocardiography, operationalized by a requirement for higher levels of respiratory or hemodynamic support. Specifically, infants were classified as intolerant if they required conventional mechanical ventilation with a fraction of inspired oxygen (FiO2) >30%, high-frequency ventilation regardless of FiO2, or one or more vasoactive-inotropic agents to maintain blood pressure. The tolerant subgroup comprised infants who tolerated the hemodynamic effects of hsPDA and required only noninvasive respiratory support, or conventional mechanical ventilation with FiO2 ≤30%, without the need for vasoactive-inotropic support. Vasoactive-inotropic agents included dopamine, dobutamine, epinephrine, norepinephrine, and vasopressin.
Definitions
IVH was defined and graded according to Papile criteria (18). NEC was defined and staged according to Bell criteria (19). BPD was defined and graded according to Jensen criteria (20). BPD grade ≥2 was defined as a requirement for supplemental oxygen via nasal cannula (>2 L/min), noninvasive positive pressure ventilation, or invasive mechanical ventilation at 36 weeks’ postmenstrual age or at discharge to home, whichever occurred first. Small for gestational age was defined as a birth weight below the 10th percentile for gestational age and sex. Clinical risk index for babies II (CRIB-II) score was calculated using sex, gestational age, birth weight, temperature, and base excess on admission (21). Antenatal corticosteroids were defined as receipt of a partial or complete course of antenatal corticosteroids prior to delivery. Chorioamnionitis was diagnosed based on clinical findings and/or histopathological examination.
Statistical analysis
Categorical variables were presented as frequencies and percentages, and continuous variables were summarized as medians with interquartile ranges (IQRs) because they were not normally distributed. In the unweighted observational cohort, between-group comparisons for categorical variables were performed using the chi-squared test, and comparisons for continuous variables were conducted using the Mann-Whitney U test. In the weighted target trial emulation cohort, descriptive statistics were calculated using survey-weighted analytical methods to account for the weighting structure. Between-group comparisons for categorical variables were performed using the Rao-Scott adjusted Chi-squared test, with P values reported based on the Rao-Scott F statistic. Comparisons of continuous variables were conducted using weighted linear regression models to obtain weight-adjusted P values. In addition, balance in baseline characteristics between groups was assessed using standardized mean differences (SMDs), with an SMD <0.10 considered indicative of adequate covariate balance.
Because treatment allocation in this study was not randomized, estimates of the effect of pharmacological treatment may be subject to confounding. Based on prior literature and clinical expertise, potential confounders were prespecified before study initiation (22,23). These included multiple births, cesarean delivery, premature rupture of membranes >24 hours, chorioamnionitis, antenatal corticosteroids, sex, gestational age, birth weight, small for gestational age, 5-minute Apgar score, CRIB-II score, the number of surfactant administrations within 3 days of life, and the level of respiratory and hemodynamic support on the day of echocardiographic assessment. Echocardiographic variables used to define cohort eligibility were not included in the propensity score model, as they were conditioned on at the design stage.
A propensity score-based approach was applied, and a weighted target trial emulation cohort was constructed using IPTW. Propensity scores were estimated using a logistic regression model, with receipt of pharmacological treatment as the dependent variable and the prespecified confounders as independent variables. IPTW weights were calculated based on the estimated propensity scores, with weights defined as 1/propensity score for treated infants and 1/(1 − propensity score) for untreated infants, to balance baseline characteristics between groups and emulate random treatment assignment. All subsequent effect estimates were derived from the IPTW-weighted cohort. Weighted log-binomial regression was used to estimate risk ratios (RRs) and corresponding 95% confidence intervals (CIs). When log-binomial models failed to converge, weighted Poisson regression was used as an alternative. All regression analyses incorporated robust standard errors clustered by study center to account for inter-facility variation in clinical practices. In analyses of the overall cohort, no additional covariate adjustment was performed beyond weighting to estimate marginal treatment effects. In subgroup analyses, regression models were forced to adjust for four key covariates (sex, gestational age, birth weight, and CRIB-II score), with further adjustment for any additional baseline covariates that remained imbalanced after weighting.
To assess the robustness of the findings, three sensitivity analyses were performed. First, analyses were repeated in the unweighted observational cohort using generalized linear mixed models, specifying a binomial distribution with a logit link function, with study center included as a random intercept to account for clustering. Second, in the non-treatment group, a small number of infants died within 2 days after echocardiographic assessment and therefore did not reach the prespecified grace period (patient 1 in Figure 1). Classifying these infants as untreated could potentially introduce immortal time bias. Accordingly, these infants were excluded, and the analyses were repeated. Third, the length of the grace period was shortened from 2 days to 1 day after echocardiographic assessment to evaluate whether the estimated effect of pharmacological treatment remained consistent across different grace periods.
There was one missing value for chorioamnionitis in this cohort. Because this missing observation was considered likely to represent the absence of the event, the value was imputed as “no”. All statistical analyses were performed using SAS version 9.4. All tests were two-sided, and a P value <0.05 was considered statistically significant.
Results
Study population
During the study period, a total of 641 preterm infants with a gestational age ≤29+6 weeks who were admitted within 3 days after birth were enrolled across CMCP-PDA participating centers. Among these infants, 492 underwent echocardiographic assessment on postnatal days 3–5, of whom 271 (55.1%) were diagnosed with hsPDA. After excluding infants who did not meet the eligibility criteria of the target trial emulation, 239 infants were ultimately included in the analysis (Figure 2).
In the final analytic cohort, the median gestational age was 27.9 weeks (IQR, 26.9–29.0 weeks), and the median birth weight was 1,040 g (IQR, 890–1,235 g). Echocardiographic assessment was performed at a median postnatal age of 4 days (IQR, 3–5 days): 114 infants (47.7%) were assessed on postnatal day 3, 65 (27.2%) on day 4, and 60 (25.1%) on day 5. Within 2 days after echocardiographic assessment, 65 infants received ibuprofen treatment, whereas 174 infants did not. Among 65 treated infants, ibuprofen was initiated at a median postnatal age of 5 days (IQR, 4–5 days): 4 infants (6.2%) initiated treatment on postnatal day 3, 21 (32.3%) on day 4, 24 (36.9%) on day 5, 12 (18.5%) on day 6, and 4 (6.2%) on day 7. After construction of the weighted target trial emulation cohort using IPTW, the effective weighted sample size was 239 infants in both the treatment and non-treatment groups.
Baseline characteristics
The distribution of IPTW weights for the treatment and non-treatment groups is shown in Figure S1. The weights for both groups were distributed within the region of common support, with no evidence of extreme weights.
Maternal and neonatal baseline characteristics of the treatment and non-treatment groups in the unweighted observational cohort and the IPTW-weighted target trial emulation cohort are presented in Table 1. Baseline characteristics between the treatment and non-treatment groups in the target trial emulation cohort were well balanced, with all SMDs <0.10.
Table 1
| Characteristic | Observational cohort (unweighted) | Emulated trial cohort (IPTW-weighted) | |||||||
|---|---|---|---|---|---|---|---|---|---|
| Non-treatment (n=174) | Treatment (n=65) | P value | SMD | Non-treatment (n=239) | Treatment (n=239) | P value | SMD | ||
| Maternal characteristics | |||||||||
| Multiple births | 61 (35.1) | 25 (38.5) | 0.62 | 0.07 | 87 (36.4) | 93 (38.9) | 0.75 | 0.05 | |
| Cesarean delivery | 103 (59.2) | 35 (53.8) | 0.45 | 0.11 | 137 (57.3) | 134 (56.1) | 0.86 | 0.03 | |
| Premature rupture of membranes >24 h | 30 (17.2) | 10 (15.4) | 0.73 | 0.05 | 40 (16.7) | 39 (16.3) | 0.90 | 0.02 | |
| Chorioamnionitis | 42 (24.1) | 21 (32.3) | 0.20 | 0.18 | 62 (25.9) | 59 (24.7) | 0.85 | 0.03 | |
| Antenatal corticosteroids | 136 (78.2) | 55 (84.6) | 0.26 | 0.17 | 190 (79.5) | 191 (79.9) | 0.98 | 0.00 | |
| Neonatal characteristics | |||||||||
| Male sex | 93 (53.4) | 31 (47.7) | 0.42 | 0.12 | 124 (51.9) | 122 (51.0) | 0.90 | 0.02 | |
| Gestational age, weeks | 28.0 [27.0–29.0] | 27.7 [26.7–29.0] | 0.33 | 0.12 | 27.9 [26.7–29.0] | 27.7 [26.8–29.0] | 0.88 | 0.02 | |
| <28 | 84 (48.3) | 36 (55.4) | 119 (49.8) | 126 (52.7) | |||||
| 28–29+6 | 90 (51.7) | 29 (44.6) | 120 (50.2) | 113 (47.3) | |||||
| Birth weight, g | 1,060 [880–1,280] | 1,020 [900–1,180] | 0.42 | 0.15 | 1,027 [870–1,239] | 1,033 [897–1,195] | 0.91 | 0.01 | |
| Small for gestational age | 8 (4.6) | 4 (6.2) | 0.87 | 0.07 | 12 (5.0) | 13 (5.4) | 0.91 | 0.02 | |
| Apgar 5 min score | 8.5 [8–9] | 8 [8–9] | 0.96 | 0.01 | 8 [7–9] | 8 [7–9] | 0.66 | 0.07 | |
| CRIB-II score† | 8 [6–10] | 8 [7–10] | 0.25 | 0.18 | 8 [6–10] | 8 [6–9] | 0.90 | 0.02 | |
| ≥2 surfactant doses within 3 days of life | 36 (20.7) | 12 (18.5) | 0.70 | 0.06 | 47 (19.7) | 41 (17.2) | 0.66 | 0.06 | |
| Respiratory support at echocardiography | 0.16 | 0.32 | >0.99 | 0.00 | |||||
| Noninvasive ventilation | 83 (47.7) | 21 (32.3) | 104 (43.5) | 104 (43.5) | |||||
| Conventional ventilation, FiO2 ≤30% | 54 (31.0) | 26 (40.0) | 80 (33.5) | 80 (33.5) | |||||
| Conventional ventilation, FiO2 >30% | 16 (9.2) | 6 (9.2) | 22 (9.2) | 21 (8.8) | |||||
| High-frequency ventilation, any FiO2 | 21 (12.1) | 12 (18.5) | 33 (13.8) | 34 (14.2) | |||||
| Inotropes at echocardiography | 21 (12.1) | 7 (10.8) | 0.78 | 0.04 | 28 (11.7) | 26 (10.9) | 0.87 | 0.02 | |
Data are presented as n (%) or median [IQR]. †, higher values represent greater illness severity. CRIB-II, clinical risk index for babies II; FiO2, fraction of inspired oxygen; IPTW, inverse probability of treatment weighting; IQR, interquartile range; SMD, standardized mean difference.
Outcomes in the target trial emulation cohort
In the IPTW-weighted target trial emulation cohort, ibuprofen treatment was associated with a lower, but not statistically significant, risk of mortality (RR, 0.51; 95% CI: 0.25–1.02; P=0.055). No statistically significant differences were observed between the treatment and non-treatment groups in the risks of BPD, IVH, NEC, or their related composite outcomes (Table 2).
Table 2
| Outcomes | Observational cohort (unweighted) | Emulated trial cohort (IPTW-weighted) | |||||||
|---|---|---|---|---|---|---|---|---|---|
| Non-treatment (n=174) | Treatment (n=65) | χ2 | P value | Non-treatment (n=239) | Treatment (n=239) | RR (95% CI)† | P value | ||
| Death | 24 (13.8) | 5 (7.7) | 1.652 | 0.19 | 33 (13.8) | 17 (7.1) | 0.51 (0.25–1.02) | 0.055 | |
| IVH ≥ grade 3 | 16 (9.2) | 5 (7.7) | 0.133 | 0.71 | 23 (9.6) | 14 (5.9) | 0.59 (0.23–1.52) | 0.27 | |
| NEC ≥ stage II | 11 (6.3) | 4 (6.2) | 0.000 | >0.99 | 15 (6.3) | 13 (5.4) | 0.88 (0.38–2.01) | 0.75 | |
| BPD ≥ grade 2‡ | 41/151 (27.2) | 25/60 (41.7) | 4.208 | 0.040 | 58/207 (28.0) | 81/222 (36.5) | 1.31 (0.82–2.11) | 0.26 | |
| Death/BPD | 64 (36.8) | 30 (46.2) | 1.742 | 0.18 | 89 (37.2) | 98 (41.0) | 1.10 (0.69–1.74) | 0.69 | |
| Death/IVH/NEC/BPD | 73 (42.0) | 34 (52.3) | 2.052 | 0.15 | 102 (42.7) | 110 (46.0) | 1.08 (0.77–1.51) | 0.66 | |
Data are presented as n (%) or n/N (%). †, weighted log-binomial regression with cluster-robust standard errors by study center, estimating marginal risk ratios without additional covariate adjustment after IPTW; ‡, among infants alive at 36 weeks of postmenstrual age, or discharge to home, if earlier. BPD, bronchopulmonary dysplasia; CI, confidence interval; IPTW, inverse probability of treatment weighting; IVH, intraventricular hemorrhage; NEC, necrotizing enterocolitis; RR, risk ratio.
Outcomes stratified by clinical tolerance subgroups
According to the prespecified clinical tolerance criteria, the study population was stratified into a tolerant subgroup (n=174, 72.8%) and an intolerant subgroup (n=65, 27.2%). Baseline characteristics in the tolerant and intolerant subgroups before and after weighting are presented in Tables S2,S3.
In the target trial emulation analyses, after IPTW with additional adjustment for key clinical covariates and any residual imbalances, ibuprofen treatment was not significantly associated with mortality in the tolerant subgroup (RR, 0.96; 95% CI: 0.23–4.00; P=0.95). In contrast, in the intolerant subgroup, ibuprofen treatment was associated with a significantly reduced risk of mortality (RR, 0.18; 95% CI: 0.06–0.51; P=0.001) (Table 3).
Table 3
| Outcomes | Observational cohort (unweighted) | Emulated trial cohort (IPTW-weighted) | |||||||
|---|---|---|---|---|---|---|---|---|---|
| Non-treatment | Treatment | χ2 | P value | Non-treatment | Treatment | RR (95% CI) | P value | ||
| Tolerant subgroup | n=129 | n=45 | n=173 | n=178 | † | ||||
| Death | 10 (7.8) | 4 (8.9) | 0.000 | >0.99 | 13 (7.5) | 13 (7.3) | 0.96 (0.23–4.00) | 0.95 | |
| IVH ≥ grade 3 | 5 (3.9) | 3 (6.7) | 0.127 | 0.72 | 7 (4.0) | 9 (5.1) | 1.17 (0.38–3.65) | 0.78 | |
| NEC ≥ stage II | 8 (6.2) | 3 (6.7) | 0.000 | >0.99 | 10 (5.8) | 9 (5.1) | 0.92 (0.24–3.50) | 0.90 | |
| BPD ≥ grade 2‡ | 36/120 (30.0) | 17/41 (41.5) | 1.819 | 0.17 | 50/161 (31.1) | 60/165 (36.4) | 1.17 (0.71–1.91) | 0.54 | |
| Death/BPD | 45 (34.9) | 21 (46.7) | 1.967 | 0.16 | 61 (35.3) | 73 (41.0) | 1.18 (0.81–1.73) | 0.38 | |
| Death/IVH/NEC/BPD | 49 (38.0) | 23 (51.1) | 2.370 | 0.12 | 66 (38.2) | 79 (44.4) | 1.18 (0.86–1.64) | 0.31 | |
| Intolerant subgroup | n=45 | n=20 | n=66 | n=61 | § | ||||
| Death | 14 (31.1) | 1 (5.0) | 3.949 | 0.047 | 20 (30.3) | 4 (6.6) | 0.18 (0.06–0.51) | 0.001 | |
| IVH ≥ grade 3 | 11 (24.4) | 2 (10.0) | 1.016 | 0.31 | 16 (24.2) | 5 (8.2) | 0.30 (0.11–0.83) | 0.02 | |
| NEC ≥ stage II | 3 (6.7) | 1 (5.0) | – | >0.99 | 5 (7.6) | 4 (6.6) | 0.97 (0.16–5.69) | 0.97 | |
| BPD ≥ grade 2‡ | 5/31 (16.1) | 8/19 (42.1) | 2.892 | 0.08 | 8/46 (17.4) | 21/57 (36.8) | 2.21 (0.71–6.89) | 0.17 | |
| Death/BPD | 19 (42.2) | 9 (45.0) | 0.044 | 0.83 | 28 (42.4) | 25 (41.0) | 1.12 (0.55–2.28) | 0.74 | |
| Death/IVH/NEC/BPD | 24 (53.3) | 11 (55.0) | 0.016 | 0.90 | 36 (54.5) | 31 (50.8) | 1.04 (0.76–1.42) | 0.81 | |
Data are presented as n (%) or n/N (%). †, weighted Poisson regression (employed after weighted log-binomial regression failed to converge), with cluster-robust standard errors by study center, and additionally adjusted for sex, gestational age, birth weight, CRIB-II score, premature rupture of membranes >24 h, chorioamnionitis, and antenatal corticosteroids; ‡, among infants alive at 36 weeks of postmenstrual age, or discharge to home, if earlier; §, weighted Poisson regression (employed after weighted log-binomial regression failed to converge), with cluster-robust standard errors by study center, and additionally adjusted for sex, gestational age, birth weight, CRIB-II score, cesarean delivery, premature rupture of membranes >24 h, chorioamnionitis, antenatal corticosteroids, and surfactant doses within 3 days of life. BPD, bronchopulmonary dysplasia; CI, confidence interval; IPTW, inverse probability of treatment weighting; IVH, intraventricular hemorrhage; NEC, necrotizing enterocolitis; RR, risk ratio.
Sensitivity analyses
Across all sensitivity analyses—including analyses using generalized linear mixed models in the unweighted cohort, exclusion of two infants in the non-treatment group who died within 2 days after echocardiographic assessment, and shortening of the grace period from 2 days to 1 day—the direction of the results was consistent with the primary analysis. Specifically, the association between ibuprofen treatment and a reduced risk of mortality remained primarily evident in the intolerant subgroup. Detailed results of the sensitivity analyses are presented in Tables S4-S8.
Discussion
Using prospective cohort data from the CMCP-PDA, we emulated a target trial to evaluate the effects of ibuprofen treatment for hsPDA on clinical outcomes in extremely preterm infants. Our findings suggest that the treatment effect of ibuprofen among infants with hsPDA identified between postnatal days 3 and 5 was strongly modified by clinical tolerance to PDA-related hemodynamic changes. Specifically, ibuprofen treatment was not associated with improved clinical outcomes among infants who were clinically tolerant of hsPDA, whereas it was associated with a reduced risk of mortality among clinically intolerant infants.
In our cohort, the majority of infants (174/239, 72.8%) demonstrated good clinical tolerance to hsPDA-related hemodynamic changes and required only low levels of respiratory or circulatory support. Consistent with the findings of most previous RCTs (6-13) and systematic reviews (24-26), early pharmacological treatment aimed at promoting ductal constriction was not associated with improved clinical outcomes in this population. Although these infants exhibited structurally significant ductal shunting, their apparent compensatory capacity for PDA-related hemodynamic changes prevented functional impairment or clinical instability, thereby limiting the potential benefit of pharmacological treatment. Moreover, while pharmacological treatment facilitates ductal constriction, it may also be accompanied by potential adverse effects, including reduced renal and intestinal perfusion and delayed enteral feeding (1). Ultimately, this pattern mirrors that observed in most studies, with increased rates of ductal closure but without corresponding improvements in clinical outcomes.
In the intolerant subgroup, ibuprofen treatment was associated with a significantly reduced risk of mortality. Nevertheless, this observed association should be interpreted with caution. Mortality is a complex outcome influenced by multiple clinical and contextual factors, and residual confounding cannot be fully excluded in an observational setting, particularly within small subgroups. Although we applied a target trial emulation framework with IPTW to mitigate confounding, unmeasured factors related to illness severity, clinician decision-making, or center-level practices may still have contributed to the observed mortality differences. Therefore, this finding should be regarded as hypothesis-generating rather than definitive evidence of a causal survival benefit.
At first glance, this finding appears inconsistent with the conclusions of most previous studies, which have generally reported overall null effects or even suggested potential harm associated with pharmacological treatment. However, this apparent discrepancy is more likely attributable to differences in effect estimation frameworks and study population selection, rather than to fundamental contradictions in the underlying evidence.
First, few previous studies were designed to explicitly differentiate infants’ clinical tolerance to PDA shunting. Clinical intolerance was rarely incorporated as a core inclusion criterion or prespecified as an effect modifier, potentially diluting benefits among those most likely to respond (6-9). Consequently, study cohorts often comprised substantial proportions of clinically stable infants with compensated PDA-related hemodynamic changes. If only a small fraction of infants experienced cardiopulmonary instability truly attributable to PDA shunting, any treatment benefit in that subgroup could be masked in marginal (overall) effect estimates, leading to an apparent null effect despite meaningful heterogeneity. In contrast, within the hsPDA population, we further operationalized intolerance using the intensity of respiratory and circulatory support. This approach more closely aligns with bedside assessment of whether PDA is causing functional compromise, and it may therefore better identify a clinically relevant target population most likely to benefit from treatment.
Second, beyond considerations related to effect modification, substantial clinical controversy and strong practice preferences have long existed in the management of PDA. A lack of clinical equipoise among clinicians has been a major barrier to recruitment in RCTs, particularly for high-risk infants (6,14). Multiple randomized trials have shown that a large proportion of otherwise eligible infants were not enrolled in randomization due to a lack of clinical equipoise (6,8,9,13). For example, in the recent NICHD (6) and Baby-OSCAR (8) trials, the proportions of eligible infants who were not enrolled were as high as 77.3% (1,662/2,149) and 48.6% (618/1,271), respectively. Notably, in the NICHD trial, two-thirds of eligible but non-enrolled infants received pharmacological treatment outside the trial (6), suggesting that infants perceived to be at higher risk may have been underrepresented in randomized studies. Furthermore, a post hoc analysis of the PDA-TOLERATE trial revealed that eligible infants who were not enrolled due to clinicians’ lack of equipoise—but who received treatment outside the trial—had lower gestational age, required greater respiratory support, and were more severely ill overall. Despite this higher baseline risk, however, these infants exhibited significantly lower mortality than those who were actually enrolled (27). This finding further highlights the potential impact of a lack of clinical equipoise on the composition and risk profile of trial populations. Taken together, these observations suggest that, to some extent, participants enrolled in prior RCTs may have been preferentially drawn from clinically more stable and PDA-tolerant infants, rather than from those experiencing true cardiopulmonary instability attributable to ductal shunting. In contrast, our study leveraged a real-world, multicenter prospective cohort encompassing the overall population of extremely preterm infants encountered in routine clinical practice. By further conducting targeted analyses among infants who were clinically intolerant of hsPDA, we were able to reveal clinically meaningful heterogeneity in treatment effects.
From a physiological perspective, infants in the intolerant subgroup required high levels of respiratory or hemodynamic support, indicating a state of pulmonary overcirculation, systemic hypoperfusion, or markedly increased myocardial workload. This constellation likely reflects an acute hemodynamic burden imposed by PDA shunting (28). In this context, ibuprofen-induced ductal constriction may reduce left-to-right shunting and pulmonary blood flow, while also improving systemic perfusion and blood pressure stability. Among infants whose hemodynamic compensation is already limited, such changes may be more likely to translate into an early survival benefit.
At present, there is no universally accepted echocardiographic definition or optimal timing for the assessment of PDA, and substantial heterogeneity exists across previous studies (1). The echocardiographic criteria used in this study incorporated ductal characteristics, markers of pulmonary overcirculation, and indicators of systemic hypoperfusion, allowing for a comprehensive assessment of PDA-related hemodynamics. Importantly, these parameters are readily measurable in routine clinical practice and demonstrate good inter-center consistency. This study focused on PDA-related hemodynamics and treatment decisions after the ultra-acute transitional period and did not evaluate pharmacological interventions initiated within the first 72 hours of life. We defined time zero as postnatal days 3–5 for echocardiographic assessment and evaluation of clinical tolerance, based on two key considerations. First, this timing reduced the inclusion of infants undergoing spontaneous ductal constriction or closure. As observed in our cohort, approximately one-third of infants (170/492) had a small or already closed ductus at the time of echocardiographic assessment. Second, by postnatal days 3–5, pulmonary vascular resistance has typically declined in most preterm infants, facilitating the development of left-to-right shunting (29). Consequently, cardiopulmonary instability observed during this period is more likely to be attributable to PDA-related hemodynamic changes. In addition to timing-related considerations, gestational age-related physiological variability may also influence PDA natural history and treatment outcomes. Infants born at the lowest gestational ages, particularly those <25 weeks’ gestation, exhibit markedly lower rates of spontaneous ductal closure and prolonged transitional circulation. Although such infants were included in our cohort, this study was not powered to evaluate gestational age-specific treatment effects, and future studies should further explore PDA management strategies in these extremely immature populations.
A major strength of this study is the use of a multicenter, prospective, real-world cohort combined with a target trial emulation approach to construct a randomized trial-like analytic framework. By explicitly defining time zero, a grace period, and treatment strategies, this design helps mitigate time-related and selection biases that commonly affect traditional observational studies. More importantly, we further incorporated stratification by clinical tolerance, enabling focused analyses of infants who may have been underrepresented in prior studies, particularly those who are clinically intolerant of PDA-related hemodynamic changes. This allowed us to identify a potential subgroup that may derive benefit from pharmacological treatment. Several limitations should also be acknowledged. First, ibuprofen was administered orally in this study. Although no increase in NEC was observed in our cohort, other gastrointestinal adverse outcomes were not systematically assessed. As oral ibuprofen, compared with intravenous ibuprofen, may raise concerns regarding gastrointestinal complications, the absence of an intravenous ibuprofen arm limits the generalizability of our findings to other routes of administration. Second, although a target trial emulation framework was applied, residual confounding and the influence of unmeasured confounders cannot be fully excluded. Third, clinical intolerance was defined based on the intensity of respiratory and hemodynamic support, which reflects real-world clinical decision-making but may still be influenced by subjective judgment and inter-center variability. Fourth, the relatively small sample size of the clinically intolerant subgroup warrants cautious interpretation of mortality and other major morbidities, and these findings should be confirmed in future studies. Finally, mortality was analyzed as an all-cause outcome without adjudication of individual causes.
Conclusions
Our findings suggest that the effects of ibuprofen treatment for hsPDA are heterogeneous, with an association with lower mortality observed among high-risk extremely preterm infants who are clinically intolerant to PDA-related shunting, but not among clinically tolerant infants. Risk stratification based on clinical tolerance may help inform decision-making for PDA management, and future randomized trials may consider specifically targeting this clinically intolerant phenotype to further evaluate these findings.
Acknowledgments
We thank all the data abstractors from the Chinese Multicenter Collaboration Platform for PDA in Extremely Preterm Infants (CMCP-PDA). Group information of the CMCP-PDA: Children’s Hospital of Nanjing Medical University; Changzhou Maternity and Child Health Care Hospital; Nanjing Women and Children’s Healthcare Hospital; General Hospital of Ningxia Medical University; The Affiliated Hospital of Qingdao University; The First Hospital Affiliated with Shandong First Medical University; Beijing Children’s Hospital, Capital Medical University; Affiliated Women’s Hospital of Jiangnan University; Northwest Women’s and Children’s Hospital; The First Affiliated Hospital of Xinjiang Medical University; The Affiliated Hospital of Xuzhou Medical University; Children’s Hospital, Zhejiang University School of Medicine; The First Affiliated Hospital with Nanjing Medical University.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://tp.amegroups.com/article/view/10.21037/tp-2026-1-0017/rc
Data Sharing Statement: Available at https://tp.amegroups.com/article/view/10.21037/tp-2026-1-0017/dss
Peer Review File: Available at https://tp.amegroups.com/article/view/10.21037/tp-2026-1-0017/prf
Funding: This work was supported by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tp.amegroups.com/article/view/10.21037/tp-2026-1-0017/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study protocol was approved by the Institutional Ethics Committee of the Children’s Hospital of Nanjing Medical University (No. 202406018-1) and endorsed by all participating hospitals. Informed consent was obtained from all individual participants’ legal guardians.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Ambalavanan N, Aucott SW, Salavitabar A, et al. Patent Ductus Arteriosus in Preterm Infants. Pediatrics 2025;155:e2025071425. [Crossref] [PubMed]
- Hamrick SEG, Sallmon H, Rose AT, et al. Patent Ductus Arteriosus of the Preterm Infant. Pediatrics 2020;146:e20201209. [Crossref] [PubMed]
- Dani C, Pratesi S. Patent ductus arteriosus and oxidative stress in preterm infants: a narrative review. Transl Pediatr 2020;9:835-9. [Crossref] [PubMed]
- Benitz WE, Backes CH. At a crossroads for early medical treatment of persistent patent ductus arteriosus in preterm infants. J Perinatol 2024;44:1534-7. [Crossref] [PubMed]
- Shelton EL, Reese J. A historical perspective of investigations into the mechanisms and management of patent ductus arteriosus. Semin Perinatol 2023;47:151723. [Crossref] [PubMed]
- Laughon MM, Thomas SM, Watterberg KL, et al. Expectant Management vs Medication for Patent Ductus Arteriosus in Preterm Infants: The PDA Randomized Clinical Trial. JAMA 2026;335:588-99. [Crossref] [PubMed]
- Potsiurko S, Dobryanskyy D, Sekretar L, et al. Randomized Noninferiority Trial of Expectant Management versus Early Treatment of Patent Ductus Arteriosus in Preterm Infants. Am J Perinatol 2024;41:730-8. [Crossref] [PubMed]
- Gupta S, Subhedar NV, Bell JL, et al. Trial of Selective Early Treatment of Patent Ductus Arteriosus with Ibuprofen. N Engl J Med 2024;390:314-25. [Crossref] [PubMed]
- Hundscheid T, Onland W, Kooi EMW, et al. Expectant Management or Early Ibuprofen for Patent Ductus Arteriosus. N Engl J Med 2023;388:980-90. [Crossref] [PubMed]
- Rozé JC, Cambonie G, Le Thuaut A, et al. Effect of Early Targeted Treatment of Ductus Arteriosus with Ibuprofen on Survival Without Cerebral Palsy at 2 Years in Infants with Extreme Prematurity: A Randomized Clinical Trial. J Pediatr 2021;233:33-42.e2. [Crossref] [PubMed]
- de Waal K, Phad N, Stubbs M, et al. A Randomized Placebo-Controlled Pilot Trial of Early Targeted Nonsteroidal Anti-Inflammatory Drugs in Preterm Infants with a Patent Ductus Arteriosus. J Pediatr 2021;228:82-86.e2. [Crossref] [PubMed]
- Sung SI, Lee MH, Ahn SY, et al. Effect of Nonintervention vs Oral Ibuprofen in Patent Ductus Arteriosus in Preterm Infants: A Randomized Clinical Trial. JAMA Pediatr 2020;174:755-63. [Crossref] [PubMed]
- Clyman RI, Liebowitz M, Kaempf J, et al. PDA-TOLERATE Trial: An Exploratory Randomized Controlled Trial of Treatment of Moderate-to-Large Patent Ductus Arteriosus at 1 Week of Age. J Pediatr 2019;205:41-48.e6. [Crossref] [PubMed]
- Clyman RI. Early Treatment of the Presymptomatic Patent Ductus Arteriosus. JAMA Pediatr 2025;179:818-20. [Crossref] [PubMed]
- McNamara PJ, Jain A, El-Khuffash A, et al. Guidelines and Recommendations for Targeted Neonatal Echocardiography and Cardiac Point-of-Care Ultrasound in the Neonatal Intensive Care Unit: An Update from the American Society of Echocardiography. J Am Soc Echocardiogr 2024;37:171-215. [Crossref] [PubMed]
- Siassi B, Noori S, Acherman R, et al. Practical Neonatal Echocardiography. New York: McGraw-Hill; 2019:163-78.
- Nawaytou H, Hills NK, Clyman RI. Patent ductus arteriosus and the risk of bronchopulmonary dysplasia-associated pulmonary hypertension. Pediatr Res 2023;94:547-54. [Crossref] [PubMed]
- Papile LA, Burstein J, Burstein R, et al. Incidence and evolution of subependymal and intraventricular hemorrhage: a study of infants with birth weights less than 1,500 gm. J Pediatr 1978;92:529-34. [Crossref] [PubMed]
- Bell MJ, Ternberg JL, Feigin RD, et al. Neonatal necrotizing enterocolitis. Therapeutic decisions based upon clinical staging. Ann Surg 1978;187:1-7. [Crossref] [PubMed]
- Jensen EA, Dysart K, Gantz MG, et al. The Diagnosis of Bronchopulmonary Dysplasia in Very Preterm Infants. An Evidence-based Approach. Am J Respir Crit Care Med 2019;200:751-9. [Crossref] [PubMed]
- Parry G, Tucker J, Tarnow-Mordi W, et al. CRIB II: an update of the clinical risk index for babies score. Lancet 2003;361:1789-91. [Crossref] [PubMed]
- Qian A, Jiang S, Gu X, et al. Treatment of patent ductus arteriosus and short-term outcomes among extremely preterm infants: a multicentre cohort study. EClinicalMedicine 2024;67:102356. [Crossref] [PubMed]
- Relangi D, Somashekar S, Jain D, et al. Changes in Patent Ductus Arteriosus Treatment Strategy and Respiratory Outcomes in Premature Infants. J Pediatr 2021;235:58-62. [Crossref] [PubMed]
- Buvaneswarran S, Wong YL, Liang S, et al. Active Treatment vs Expectant Management of Patent Ductus Arteriosus in Preterm Infants: A Meta-Analysis. JAMA Pediatr 2025;179:877-85. [Crossref] [PubMed]
- Mitra S, Scrivens A, von Kursell AM, et al. Early treatment versus expectant management of hemodynamically significant patent ductus arteriosus for preterm infants. Cochrane Database Syst Rev 2020;12:CD013278. [Crossref] [PubMed]
- Mitra S, Florez ID, Tamayo ME, et al. Association of Placebo, Indomethacin, Ibuprofen, and Acetaminophen With Closure of Hemodynamically Significant Patent Ductus Arteriosus in Preterm Infants: A Systematic Review and Meta-analysis. JAMA 2018;319:1221-38. [Crossref] [PubMed]
- Liebowitz M, Katheria A, Sauberan J, et al. Lack of Equipoise in the PDA-TOLERATE Trial: A Comparison of Eligible Infants Enrolled in the Trial and Those Treated Outside the Trial. J Pediatr 2019;213:222-226.e2. [Crossref] [PubMed]
- Backes CH, Hill KD, Shelton EL, et al. Patent Ductus Arteriosus: A Contemporary Perspective for the Pediatric and Adult Cardiac Care Provider. J Am Heart Assoc 2022;11:e025784. [Crossref] [PubMed]
- Hillman NH. Increasing pulmonary blood flow at birth: the nerve of the baby. J Physiol 2017;595:1437-8. [Crossref] [PubMed]

