Introduction: Milrinone is administered after patent ductus arteriosus (PDA) ligation to prevent and treat postoperative hemodynamic instability (i.e., postligation cardiac syndrome). We aimed to explore the effectiveness of milrinone on in-hospital outcomes in infants who underwent PDA ligation using a nationwide inpatient database in Japan. Methods: Using the Japanese Diagnosis Procedure Combination database, we identified patients who received milrinone after PDA ligation (n = 428) in neonatal intensive care units between July 2010 and March 2021 and those who did not (n = 3,392). We conducted a 1:4 propensity score-matched analysis with adjustment for background characteristics (e.g., gestational age, birth weight, comorbidities, preoperative treatments, and hospital background) to compare morbidities (bronchopulmonary dysplasia, intraventricular hemorrhage, necrotizing enterocolitis, and retinopathy of prematurity), mortality, total hospitalization costs, and other outcomes. For sensitivity analysis, we performed an overlap propensity score-weighted analysis. Results: In-hospital morbidity, bronchopulmonary dysplasia, intraventricular hemorrhage, and necrotizing enterocolitis occurred in 58%, 48%, 9.5%, and 7.1% of patients, respectively; the in-hospital mortality was 5.4%. After 1:4 propensity score matching, no significant difference was observed regarding mortality (7.1 vs. 5.7%), in-hospital morbidity (55 vs. 50%), bronchopulmonary dysplasia (44 vs. 41%), intraventricular hemorrhage (7.8 vs. 9.1%), necrotizing enterocolitis (8.5 vs. 8.9%), retinopathy of prematurity (21 vs. 22%), or total hospitalization costs (median: approximately 86,000 vs. 82,000 US dollars) between milrinone users (n = 425) and nonusers (n = 1,698). Sensitivity analyses yielded consistent results. Conclusions: Milrinone use after PDA ligation was not associated with improved in-hospital outcomes, such as mortality and morbidity.

Patent ductus arteriosus (PDA) – a major cardiovascular complication among preterm infants – is associated with mortality and severe comorbidities, including bronchopulmonary dysplasia (BPD), necrotizing enterocolitis (NEC), and retinopathy of prematurity (ROP) [1‒4]. Pharmacotherapy or surgery for PDA is required in 20–70% of very preterm infants [4‒7]; if pharmacotherapy with ibuprofen, indomethacin, or paracetamol is ineffective, PDA ligation is performed [1, 6, 8]. From 2006 to 2015, 2.9–8.4% of preterm infants in the USA required PDA ligation [7].

Postligation cardiac syndrome (PLCS) is a phenomenon of hemodynamic instability that manifests as postoperative hypoxemia and hypotension and requires prolonged respiratory support [9]. PLCS occurs in 7.2–44% of patients who undergo PDA ligation [10‒14] and can result in mortality, morbidities (such as BPD and NEC), and the need for home oxygen therapy [10, 13, 14]. Milrinone, a selective phosphodiesterase III inhibitor, is administered following PDA ligation for the prevention and treatment of PLCS [15‒18]. A retrospective study involving 52 patients who underwent PDA ligation reported that milrinone reduced the development of PLCS, catecholamine use, respiratory support, and steroid use [17]. Although some observational studies have reported that milrinone therapy after PDA ligation can reduce oxygenation or ventilation failure, cardiac function, and pharmacokinetics, there is little evidence regarding postoperative mortality, BPD, NEC, and ROP [18, 19]. Therefore, using a nationwide inpatient database in Japan, we aimed to explore the effectiveness of milrinone on in-hospital outcomes in infants who underwent PDA ligation.

Database

This nationwide, retrospective cohort study was conducted using the Diagnostic Procedure Combination database. This database includes the administrative claims data and discharge abstracts of approximately 8,000,000 inpatients per year from >1,200 hospitals throughout Japan [20]. All 82 university hospitals in Japan are required to participate in the database, whereas the participation of community hospitals was voluntary. The database covers approximately 90% of hospitals with neonatal intensive care units (NICUs) in Japan [21]. The requirement for informed consent was waived due to the anonymity of the patient database, and the Institutional Review Board of the University of Tokyo approved this study (approval number: 3501-[5] [May 19th, 2021]).

The database comprised the following patient data: sex, gestational age, birth weight, main diagnoses, comorbidities on admission, and complications after admission recorded in Japanese text data and using the International Classification of Disease, Tenth Revision (ICD-10) codes indexed by the original Japanese procedure codes, patient transportation from other hospitals or in-hospital births, length of stay, discharge status, and hospitalization costs. Hospitalization costs were based on reference prices in the Japanese national fee schedule for item-by-item prices for surgical, pharmaceutical, laboratory, and other inpatient services. All abstract discharge data for each patient were recorded at discharge by attending physicians. Favorable sensitivities and specificities for diagnosis and procedure records in the database were demonstrated in several previous validation studies [22‒24]; additionally, several perinatal and neonatal studies were conducted using this database [25‒27].

Protocol

We retrospectively identified patients admitted to the NICU who underwent PDA ligation between July 2010 and March 2021. The original Japanese procedure codes for PDA ligation and NICU admission were used to identify the patients. Patients who simultaneously underwent heart surgery other than PDA ligation were excluded. We divided patients into two groups to investigate the effectiveness of early-phase initiation of milrinone after PDA ligation. The milrinone group comprised patients who were administered milrinone on the first or second day after PDA ligation, and the non-milrinone group comprised those who were not administered milrinone after PDA ligation.

The primary outcomes were in-hospital morbidity (BPD, intraventricular hemorrhage [IVH], NEC, and ROP) and mortality. BPD was defined as oxygen use, noninvasive positive pressure ventilation, or ventilation at 36 weeks of corrected gestational age in preterm infants born under 32 weeks of gestational age. NEC was defined as undergoing abdominal surgery coded with the original Japanese procedure codes or diagnoses recorded with ICD-10 codes. ROP was defined as laser photocoagulation or antivascular endothelial growth factor treatment. Definitions of the ICD-10 codes for IVH and NEC are listed in Supplementary Table 1 (for all online suppl. material, see https://doi.org/10.1159/000533598). Secondary outcomes were postoperative steroid and blood transfusion use, discharge destination, post-discharge device use (including home oxygen therapy, tracheostomy, and home mechanical ventilation), length of hospital stay, and total hospitalization cost. We defined the currency exchange rate as 130 Japanese yen per 1 US dollar.

We examined the following patient characteristics: sex, gestational age, birth weight, in-hospital birth, general comorbidities, and the presence of neonatal asphyxia. Gestational age was categorized into four groups: <28 weeks (extremely preterm), 28–31 weeks (moderately preterm), 32–36 weeks (late preterm), and ≥37 weeks (term). Birth weight was categorized into four groups: <1,000 g (extremely low), 1,000–1,499 g (very low), 1,500–2,499 g (low), and ≥2,500 g (normal). General comorbidities were defined according to the pediatric complex chronic conditions (CCCs) classification system, version 2 [28], and categorized into four groups according to the number of CCCs (0, 1, 2, and ≥3). Additionally, we investigated the following treatment characteristics: preoperative drug use (cyclooxygenase inhibitors [indomethacin and ibuprofen], catecholamines [dopamine and dobutamine], adrenalin, milrinone, steroids, antibiotics, sedative agents [fentanyl, midazolam, phenobarbital, and dexmedetomidine], and muscle relaxants [rocuronium and vecuronium]), preoperative transfusions (red blood cells, platelets, fresh frozen plasma, albumin, and gamma-globulin), preoperative procedures (central venous catheterization, inhaled nitric oxide, intratracheal surfactant administration, intubation, and jaundice phototherapy), hospital type, hospital volume, and fiscal year. Regarding the hospital type, we investigated whether patients were admitted to a general perinatal medical center or a teaching hospital. The NICU in a general perinatal medical center in Japan is similar to level 3 NICUs in the USA [29]. The hospital volume was defined as the number of PDA ligations performed annually at each hospital and was categorized into tertiles (low, medium, and high). The fiscal years were categorized into two halves: 2010–2015 and 2016–2021. These background characteristics were adjusted in the subsequent analyses.

Statistical Analysis

We conducted propensity score matching at a 1:4 ratio to compare outcomes between the two groups [30, 31]. We calculated the propensity scores of the patients and their treatment characteristics using a logistic regression model. Each patient in the milrinone group was matched with 4 patients with the closest estimated propensity scores within a caliper (≤0.2 of the pooled standard deviation) in the non-milrinone group. Nearest-neighbor matching with replacement was used to obtain matched samples by matching the subjects. Additionally, we calculated the standardized differences to examine the balance in the baseline characteristics of patients between the two groups for all patients, as well as for the 1:4 propensity score-matched paired cohort. A standardized difference of 10% implied a negligible difference between the two groups [31]. In the 1:4 propensity score-matched paired cohort, we performed a χ2 test to compare proportions for categorical outcomes and the Kruskal-Wallis test to compare medians for continuous outcomes. A subgroup analysis of patients with a gestational age <28 weeks was also conducted because preterm infants with <28 weeks of gestational age often require PDA ligation and intensive cardiotropic support after PDA ligation [1, 2, 12].

For sensitivity analysis, we conducted an overlap propensity score-weighted analysis for the main and subgroup analyses to examine the robustness of our analysis. Overlap propensity score-weighted analysis is an extension of propensity score analysis that balances covariates between two groups and imitates a randomized trial without excluding most observations from the study sample [32, 33]. This analysis can compare highly comparable populations without being overly influenced by outliers that are almost always or never treated [33].

The threshold for significance was set at a level of 0.05, and all reported p values were two-sided. All statistical analyses were conducted using Stata/MP 17.0 (StataCorp, College Station, TX, USA).

We identified 4,117 patients who underwent PDA ligation in the NICU between July 2010 and March 2021. We excluded 297 patients who underwent other cardiac surgeries during PDA ligation. Of the 3,820 eligible patients, the milrinone and non-milrinone groups comprised 428 and 3,392 patients, respectively (Fig. 1).

Fig. 1.

Patient flow. NICU, neonatal intensive care unit; PDA, patent ductus arteriosus.

Fig. 1.

Patient flow. NICU, neonatal intensive care unit; PDA, patent ductus arteriosus.

Close modal

Table 1 describes the background characteristics of all patients, as well as the 1:4 propensity score-matched patients. Before matching, the milrinone group was more likely to have a higher gestational age, birth weight, and number of CCCs than the non-milrinone group. After 1:4 propensity score matching, the milrinone and non-milrinone groups comprised 1,698 patients; the background characteristics were well-balanced between the two groups.

Table 1.

Background characteristics of all patients and 1:4 propensity score-matched patients

All patients1:4 propensity score-matched patients
milrinonenon-milrinoneASD*milrinonenon-milrinoneASD*
n = 428n = 3,392(%)n = 425n = 1,698(%)
Patient characteristics 
Male 198 (46) 1,571 (46) 0.1 196 (46) 799 (47) 1.9 
Gestational age category, weeks 
 <28 224 (52) 2,209 (65) 26.2 224 (53) 933 (55) 4.5 
 28–31 46 (11) 404 (12) 3.7 46 (11) 193 (11) 1.7 
 32–36 61 (14) 263 (7.8) 20.9 59 (14) 198 (12) 6.7 
 ≥37 73 (17) 309 (9.1) 23.7 72 (17) 291 (17) 0.5 
 Missing data 24 (5.6) 207 (6.1) 2.1 24 (5.6) 83 (4.9) 3.4 
Birth weight category, g 
 <1,000 239 (56) 2,341 (69) 27.5 239 (56) 998 (59) 5.1 
 1,000–1,499 55 (13) 454 (13) 1.6 55 (13) 212 (12) 1.4 
 1,500–2,499 78 (18) 359 (11) 21.9 76 (18) 278 (16) 4.0 
 ≥2,500 52 (12) 214 (6.3) 20.3 51 (12) 195 (11) 1.6 
 Missing data (0.93) 24 (0.71) 2.5 (0.94) 15 (0.88) 0.6 
In-hospital birth 284 (66) 2,494 (74) 15.7 283 (67) 1,116 (66) 1.8 
Number of CCC 
 0 129 (30) 1,051 (31) 1.8 129 (30) 488 (29) 3.5 
 1 232 (54) 1,929 (57) 5.4 230 (54) 976 (57) 6.8 
 2 50 (12) 363 (11) 3.1 49 (12) 183 (11) 2.4 
 ≥3 17 (4.0) 49 (1.4) 15.6 17 (4.0) 51 (3.0) 5.4 
Neonatal asphyxia 235 (55) 1,994 (59) 7.8 234 (55) 942 (55) 0.8 
Treatment characteristics 
Preoperative drug use           
 Cyclooxygenase inhibitors 244 (57) 2,407 (71) 29.4 244 (57) 955 (56) 2.4 
 Catecholamines 243 (57) 2,085 (61) 9.6 241 (57) 912 (54) 6.0 
 Adrenalin 20 (4.7) 80 (2.4) 12.6 19 (4.5) 63 (3.7) 3.8 
 Milrinone 52 (12) 74 (2.2) 39.4 49 (12) 210 (12) 2.6 
 Steroids 202 (47) 1,638 (48) 2.2 201 (47) 841 (50) 4.5 
 Antibiotics 290 (68) 2,470 (73) 11.1 287 (68) 1,136 (67) 1.3 
 Sedative agents 190 (44) 1,336 (39) 10.2 189 (45) 771 (45) 1.9 
 Muscle relaxant 20 (4.7) 104 (3.1) 8.3 19 (4.5) 78 (4.6) 0.6 
Preoperative transfusion 
 Red blood cell 247 (58) 1,893 (56) 3.8 245 (58) 986 (58) 0.9 
 Platelets 49 (25) 377 (24) 1.7 49 (12) 178 (11) 3.3 
 Fresh frozen plasma 110 (26) 910 (27) 2.6 109 (26) 415 (24) 2.8 
 Albumin 116 (27) 935 (28) 1.0 115 (27) 478 (28) 2.4 
 Gamma-globulin 140 (33) 1,313 (39) 12.5 140 (33) 559 (33) 0.0 
Preoperative procedures 
 Central venous catheterization 323 (75) 2,556 (75) 0.3 321 (76) 1,253 (74) 4.0 
 Inhaled nitric oxide 26 (6.1) 252 (7.4) 5.4 26 (6.1) 116 (6.8) 2.9 
 Intratracheal surfactant 214 (50) 2,038 (60) 20.4 214 (50) 899 (53) 5.2 
 Intubation 357 (83) 2,970 (88) 11.8 355 (84) 1,398 (82) 3.2 
 Phototherapy for jaundice 241 (56) 2,201 (65) 17.6 240 (56) 904 (53) 6.5 
Type of intervention 
 Laparoscopy (0.23) 47 (1.4) 12.9 (0.24) (0.24) 0.0 
 Catheter closure 12 (2.8) 76 (2.2) 3.6 12 (2.8) 46 (2.7) 0.7 
Teaching hospital 395 (92) 3,167 (93) 4.2 393 (92) 1,594 (94) 5.6 
General perinatal medical center 421 (98) 3,344 (99) 1.8 418 (98) 1,673 (99) 1.4 
Hospital volume 
 <4 168 (39) 1,068 (31) 16.3 166 (39) 619 (36) 5.4 
 4–6 105 (25) 1,182 (35) 22.7 105 (25) 412 (24) 1.0 
 >6 155 (36) 1,142 (34) 5.3 154 (36) 667 (39) 6.3 
Fiscal year of surgery 
 2010–2015 207 (48) 1,937 (57) 17.6 206 (48) 846 (50) 2.7 
 2016–2021 221 (52) 1,455 (43) 17.6 219 (52) 852 (50) 2.7 
All patients1:4 propensity score-matched patients
milrinonenon-milrinoneASD*milrinonenon-milrinoneASD*
n = 428n = 3,392(%)n = 425n = 1,698(%)
Patient characteristics 
Male 198 (46) 1,571 (46) 0.1 196 (46) 799 (47) 1.9 
Gestational age category, weeks 
 <28 224 (52) 2,209 (65) 26.2 224 (53) 933 (55) 4.5 
 28–31 46 (11) 404 (12) 3.7 46 (11) 193 (11) 1.7 
 32–36 61 (14) 263 (7.8) 20.9 59 (14) 198 (12) 6.7 
 ≥37 73 (17) 309 (9.1) 23.7 72 (17) 291 (17) 0.5 
 Missing data 24 (5.6) 207 (6.1) 2.1 24 (5.6) 83 (4.9) 3.4 
Birth weight category, g 
 <1,000 239 (56) 2,341 (69) 27.5 239 (56) 998 (59) 5.1 
 1,000–1,499 55 (13) 454 (13) 1.6 55 (13) 212 (12) 1.4 
 1,500–2,499 78 (18) 359 (11) 21.9 76 (18) 278 (16) 4.0 
 ≥2,500 52 (12) 214 (6.3) 20.3 51 (12) 195 (11) 1.6 
 Missing data (0.93) 24 (0.71) 2.5 (0.94) 15 (0.88) 0.6 
In-hospital birth 284 (66) 2,494 (74) 15.7 283 (67) 1,116 (66) 1.8 
Number of CCC 
 0 129 (30) 1,051 (31) 1.8 129 (30) 488 (29) 3.5 
 1 232 (54) 1,929 (57) 5.4 230 (54) 976 (57) 6.8 
 2 50 (12) 363 (11) 3.1 49 (12) 183 (11) 2.4 
 ≥3 17 (4.0) 49 (1.4) 15.6 17 (4.0) 51 (3.0) 5.4 
Neonatal asphyxia 235 (55) 1,994 (59) 7.8 234 (55) 942 (55) 0.8 
Treatment characteristics 
Preoperative drug use           
 Cyclooxygenase inhibitors 244 (57) 2,407 (71) 29.4 244 (57) 955 (56) 2.4 
 Catecholamines 243 (57) 2,085 (61) 9.6 241 (57) 912 (54) 6.0 
 Adrenalin 20 (4.7) 80 (2.4) 12.6 19 (4.5) 63 (3.7) 3.8 
 Milrinone 52 (12) 74 (2.2) 39.4 49 (12) 210 (12) 2.6 
 Steroids 202 (47) 1,638 (48) 2.2 201 (47) 841 (50) 4.5 
 Antibiotics 290 (68) 2,470 (73) 11.1 287 (68) 1,136 (67) 1.3 
 Sedative agents 190 (44) 1,336 (39) 10.2 189 (45) 771 (45) 1.9 
 Muscle relaxant 20 (4.7) 104 (3.1) 8.3 19 (4.5) 78 (4.6) 0.6 
Preoperative transfusion 
 Red blood cell 247 (58) 1,893 (56) 3.8 245 (58) 986 (58) 0.9 
 Platelets 49 (25) 377 (24) 1.7 49 (12) 178 (11) 3.3 
 Fresh frozen plasma 110 (26) 910 (27) 2.6 109 (26) 415 (24) 2.8 
 Albumin 116 (27) 935 (28) 1.0 115 (27) 478 (28) 2.4 
 Gamma-globulin 140 (33) 1,313 (39) 12.5 140 (33) 559 (33) 0.0 
Preoperative procedures 
 Central venous catheterization 323 (75) 2,556 (75) 0.3 321 (76) 1,253 (74) 4.0 
 Inhaled nitric oxide 26 (6.1) 252 (7.4) 5.4 26 (6.1) 116 (6.8) 2.9 
 Intratracheal surfactant 214 (50) 2,038 (60) 20.4 214 (50) 899 (53) 5.2 
 Intubation 357 (83) 2,970 (88) 11.8 355 (84) 1,398 (82) 3.2 
 Phototherapy for jaundice 241 (56) 2,201 (65) 17.6 240 (56) 904 (53) 6.5 
Type of intervention 
 Laparoscopy (0.23) 47 (1.4) 12.9 (0.24) (0.24) 0.0 
 Catheter closure 12 (2.8) 76 (2.2) 3.6 12 (2.8) 46 (2.7) 0.7 
Teaching hospital 395 (92) 3,167 (93) 4.2 393 (92) 1,594 (94) 5.6 
General perinatal medical center 421 (98) 3,344 (99) 1.8 418 (98) 1,673 (99) 1.4 
Hospital volume 
 <4 168 (39) 1,068 (31) 16.3 166 (39) 619 (36) 5.4 
 4–6 105 (25) 1,182 (35) 22.7 105 (25) 412 (24) 1.0 
 >6 155 (36) 1,142 (34) 5.3 154 (36) 667 (39) 6.3 
Fiscal year of surgery 
 2010–2015 207 (48) 1,937 (57) 17.6 206 (48) 846 (50) 2.7 
 2016–2021 221 (52) 1,455 (43) 17.6 219 (52) 852 (50) 2.7 

ASD, absolute standardized difference; CCC, complex chronic condition (pediatric complex chronic conditions classification system, version 2).

*An ASD of ≤10% denotes a negligible difference between the two groups.

The crude outcomes of all patients are shown in Table 2. Among all patients, in-hospital morbidity occurred in 2,215 (58%) patients; 206 (5.4%) patients died. After 1:4 propensity score matching, no significant differences were observed in mortality (7.1 vs. 5.7%), in-hospital morbidity (55 vs. 50%), length of hospital stay (median: 113 vs. 102 days), hospitalization costs (median: approximately 86,000 vs. 82,000 US dollars), or other outcomes between the two groups. The results of the overlap propensity score-weighted analysis, wherein the weighted background characteristics were balanced (online suppl. Table 2), were comparable to those of the main analyses (Table 3).

Table 2.

Comparisons of outcomes between the milrinone group and non-milrinone group among all patients and 1:4 propensity score-matched patients

All patients1:4 propensity score-matched patients
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 428n = 3,392n = 425n = 1,698
In-hospital morbidity 233 (54) 1,982 (58) 0.12 233 (55) 855 (50) 0.099 
 Bronchopulmonary dysplasia 187 (44) 1,648 (49) 0.029 187 (44) 690 (41) 0.31 
 Intraventricular hemorrhage 33 (7.7) 330 (9.7) 0.18 33 (7.8) 154 (9.1) 0.40 
 Necrotizing enterocolitis 36 (8.4) 313 (9.2) 0.58 36 (8.5) 151 (8.9) 0.78 
 Retinopathy of prematurity 91 (21) 864 (25) 0.058 91 (21) 367 (22) 0.93 
In-hospital mortality 30 (7.0) 176 (5.2) 0.12 30 (7.1) 97 (5.7) 0.30 
Postoperative steroid use 218 (51) 1,646 (49) 0.35 217 (51) 846 (50) 0.65 
Postoperative blood transfusion 207 (48) 1,527 (45) 0.19 205 (48) 769 (45) 0.28 
Discharged to home 296 (69) 2,451 (72) 0.18 293 (69) 1,205 (71) 0.41 
Post-discharge device use 65 (15) 518 (15) 0.96 65 (15) 288 (17) 0.41 
All patients1:4 propensity score-matched patients
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 428n = 3,392n = 425n = 1,698
In-hospital morbidity 233 (54) 1,982 (58) 0.12 233 (55) 855 (50) 0.099 
 Bronchopulmonary dysplasia 187 (44) 1,648 (49) 0.029 187 (44) 690 (41) 0.31 
 Intraventricular hemorrhage 33 (7.7) 330 (9.7) 0.18 33 (7.8) 154 (9.1) 0.40 
 Necrotizing enterocolitis 36 (8.4) 313 (9.2) 0.58 36 (8.5) 151 (8.9) 0.78 
 Retinopathy of prematurity 91 (21) 864 (25) 0.058 91 (21) 367 (22) 0.93 
In-hospital mortality 30 (7.0) 176 (5.2) 0.12 30 (7.1) 97 (5.7) 0.30 
Postoperative steroid use 218 (51) 1,646 (49) 0.35 217 (51) 846 (50) 0.65 
Postoperative blood transfusion 207 (48) 1,527 (45) 0.19 205 (48) 769 (45) 0.28 
Discharged to home 296 (69) 2,451 (72) 0.18 293 (69) 1,205 (71) 0.41 
Post-discharge device use 65 (15) 518 (15) 0.96 65 (15) 288 (17) 0.41 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 113 (41–154) 114 (54–149.5) 0.70 113 (41–154) 102 (32–146) 0.071 
Total hospitalization cost, US dollars 85,096 (37,636–119,227) 94,173 (45,568–116,241) 0.16 85,943 (38,039–119,267) 81,968 (35,294–114,142) 0.15 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 113 (41–154) 114 (54–149.5) 0.70 113 (41–154) 102 (32–146) 0.071 
Total hospitalization cost, US dollars 85,096 (37,636–119,227) 94,173 (45,568–116,241) 0.16 85,943 (38,039–119,267) 81,968 (35,294–114,142) 0.15 

IQR, interquartile range.

Table 3.

Overlap propensity score-weighted analyses comparing the milrinone group and the non-milrinone group

All patientsPatients under 28 weeks of GA
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 352n = 352n = 191n = 191
In-hospital morbidity 198 (56) 180 (51) 0.056 147 (77) 139 (73) 0.18 
 Bronchopulmonary dysplasia 161 (46) 148 (42) 0.22 127 (66) 119 (62) 0.43 
 Intraventricular hemorrhage 27 (7.7) 29 (8.2) 0.70 21 (11) 22 (12) 0.78 
 Necrotizing enterocolitis 28 (8.0) 30 (8.5) 0.67 18 (9.4) 21 (11) 0.54 
 Retinopathy of prematurity 77 (22) 76 (22) 0.85 66 (35) 66 (35) 1.00 
In-hospital mortality 24 (6.8) 21 (6.0) 0.56 12 (6.3) 12 (6.3) 0.95 
Postoperative steroid use 181 (51) 168 (48) 0.18 128 (67) 116 (61) 0.08 
Postoperative blood transfusion 168 (48) 153 (43) 0.12 110 (58) 101 (53) 0.21 
Discharged to home 243 (69) 244 (69) 0.87 128 (67) 134 (70) 0.40 
Post-discharge device use 53 (15) 52 (15) 0.88 35 (18) 37 (19) 0.73 
All patientsPatients under 28 weeks of GA
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 352n = 352n = 191n = 191
In-hospital morbidity 198 (56) 180 (51) 0.056 147 (77) 139 (73) 0.18 
 Bronchopulmonary dysplasia 161 (46) 148 (42) 0.22 127 (66) 119 (62) 0.43 
 Intraventricular hemorrhage 27 (7.7) 29 (8.2) 0.70 21 (11) 22 (12) 0.78 
 Necrotizing enterocolitis 28 (8.0) 30 (8.5) 0.67 18 (9.4) 21 (11) 0.54 
 Retinopathy of prematurity 77 (22) 76 (22) 0.85 66 (35) 66 (35) 1.00 
In-hospital mortality 24 (6.8) 21 (6.0) 0.56 12 (6.3) 12 (6.3) 0.95 
Postoperative steroid use 181 (51) 168 (48) 0.18 128 (67) 116 (61) 0.08 
Postoperative blood transfusion 168 (48) 153 (43) 0.12 110 (58) 101 (53) 0.21 
Discharged to home 243 (69) 244 (69) 0.87 128 (67) 134 (70) 0.40 
Post-discharge device use 53 (15) 52 (15) 0.88 35 (18) 37 (19) 0.73 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 114 (41–154) 102 (37–145) 0.70 131 (96–168) 126 (91–159) 0.30 
Total hospitalization cost, US dollars 87,687 (37,232–119,260) 82,908 (36,015–113,627) 0.16 110,124 (79,111–124,797) 105,850 (77,565–121,370) 0.30 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 114 (41–154) 102 (37–145) 0.70 131 (96–168) 126 (91–159) 0.30 
Total hospitalization cost, US dollars 87,687 (37,232–119,260) 82,908 (36,015–113,627) 0.16 110,124 (79,111–124,797) 105,850 (77,565–121,370) 0.30 

IQR, interquartile range; GA, gestational age.

In the subgroup analysis of patients born at <28 weeks of gestational age, the milrinone and non-milrinone groups comprised 224 and 2,209 patients, respectively (Fig. 1). After 1:4 propensity score matching (n = 221 vs. n = 882; online suppl. Table 3), the outcomes did not differ significantly, except for postoperative steroid use; the milrinone group received steroids after ligation more often than the non-milrinone group (68 vs. 60%; p = 0.025) (Table 4). In the overlapping propensity score-weighted analysis, no significant differences were observed in any outcomes, including steroid use (Table 3).

Table 4.

Comparisons of outcomes between milrinone group and non-milrinone group in all patients and 1:4 propensity score-matched patients with a gestational age <28 weeks

All patients1:4 propensity score-matched patients
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 224n = 2,209n = 221n = 882
In-hospital morbidity 173 (77) 1,651 (75) 0.41 170 (77) 625 (71) 0.072 
 Bronchopulmonary dysplasia 148 (66) 1,417 (64) 0.80 145 (66) 532 (60) 0.27 
 Intraventricular hemorrhage 26 (12) 284 (13) 0.59 24 (11) 109 (12) 0.54 
 Necrotizing enterocolitis 23 (10) 248 (11) 0.66 21 (10) 87 (10) 0.87 
 Retinopathy of prematurity 78 (35) 795 (36) 0.73 75 (34) 285 (32) 0.65 
In-hospital mortality 15 (6.7) 125 (5.7) 0.53 15 (6.8) 53 (6.0) 0.67 
Postoperative steroid use 153 (68) 1,294 (59) 0.005 150 (68) 526 (60) 0.025 
Postoperative blood transfusion 131 (59) 1,166 (53) 0.10 128 (58) 464 (53) 0.16 
Discharged to home 149 (67) 1,612 (73) 0.039 146 (66) 584 (66) 0.97 
Post-discharge device use 41 (18) 405 (18) 0.99 39 (18) 181 (21) 0.34 
All patients1:4 propensity score-matched patients
milrinonenon-milrinonep valuemilrinonenon-milrinonep value
n = 224n = 2,209n = 221n = 882
In-hospital morbidity 173 (77) 1,651 (75) 0.41 170 (77) 625 (71) 0.072 
 Bronchopulmonary dysplasia 148 (66) 1,417 (64) 0.80 145 (66) 532 (60) 0.27 
 Intraventricular hemorrhage 26 (12) 284 (13) 0.59 24 (11) 109 (12) 0.54 
 Necrotizing enterocolitis 23 (10) 248 (11) 0.66 21 (10) 87 (10) 0.87 
 Retinopathy of prematurity 78 (35) 795 (36) 0.73 75 (34) 285 (32) 0.65 
In-hospital mortality 15 (6.7) 125 (5.7) 0.53 15 (6.8) 53 (6.0) 0.67 
Postoperative steroid use 153 (68) 1,294 (59) 0.005 150 (68) 526 (60) 0.025 
Postoperative blood transfusion 131 (59) 1,166 (53) 0.10 128 (58) 464 (53) 0.16 
Discharged to home 149 (67) 1,612 (73) 0.039 146 (66) 584 (66) 0.97 
Post-discharge device use 41 (18) 405 (18) 0.99 39 (18) 181 (21) 0.34 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 132 (95–167.5) 129 (99–160) 0.46 131 (95–167) 127 (61–159) 0.071 
Total hospitalization cost, US dollars 110,420 (77,900–125,337) 106,822 (84,287–121,611) 0.35 110,124 (77,334–124,373) 105,996 (67,024–120,261) 0.053 
MedianIQRMedianIQRMedianIQRMedianIQR
Length of hospital stay, days 132 (95–167.5) 129 (99–160) 0.46 131 (95–167) 127 (61–159) 0.071 
Total hospitalization cost, US dollars 110,420 (77,900–125,337) 106,822 (84,287–121,611) 0.35 110,124 (77,334–124,373) 105,996 (67,024–120,261) 0.053 

IQR, interquartile range.

In this nationwide database study, we explored the effects of milrinone administration on short-term outcomes following PDA ligation using two propensity score analyses. Postoperative milrinone use was not significantly associated with reduced in-hospital morbidity or mortality. The subgroup analysis of patients born before 28 weeks of gestation demonstrated results similar to those in the main analysis.

PDA ligation can cause hemodynamic instability (i.e., PLCS) associated with severe morbidity [12, 13, 34]. Previous observational studies demonstrated improved cardiac function in patients who received milrinone after PDA ligation [15, 17‒19] as milrinone can reduce biventricular afterload [35]; however, whether milrinone therapy after PDA ligation improves clinical outcomes (i.e., mortality and severe morbidity) remains unclear. To our knowledge, the present study is the largest retrospective cohort study to compare the outcomes of neonates treated with and without milrinone after PDA ligation.

In the present study, milrinone was used in 11% of patients who underwent PDA ligation, while in a previous study of 85 patients who underwent PDA ligation, 19% were administered milrinone [16]. The patients in our study had a lower gestational age and birth weight than those in the previous study, and patients with a high gestational age may have tended to receive milrinone after ligation. Indeed, the milrinone group had a higher gestational age and birth weight than the non-milrinone group.

Our findings revealed that milrinone use after PDA ligation was not associated with decreased morbidity, mortality, discharge destination, post-discharge device use, length of hospital stay, or hospitalization costs. These findings are inconsistent with those of previous retrospective studies that reported the effectiveness of PDA ligation for hemodynamic instability [15, 17‒19]. However, these previous studies were conducted in a single institute, comprised a small number of patients, and did not adjust for clinical background [15, 17‒19]. The present study involved approximately 4,000 patients and was adjusted for patient backgrounds using propensity score matching. Furthermore, although previous studies have investigated postoperative cardiac function, hemodynamic instability, and pharmacokinetics [15, 17‒19], these outcomes may not have direct clinical relevance. The results of the present study suggest that milrinone administration after PDA ligation may be ineffective in improving short-term clinical outcomes. Clinicians should be aware of the uncertain efficacy regarding milrinone initiation after ligation.

In the subgroup analysis of patients with a gestational age <28 weeks, the results were almost consistent with those of the main analysis, excluding postoperative steroid use. That is, although a previous study reported the decreased administration of postoperative steroids in patients treated with milrinone, the milrinone group received steroids after ligation more often than the non-milrinone group [11]. Since milrinone is associated with decreasing blood pressure due to excessive vasodilation [11], steroids would have been used to raise blood pressure, specifically in preterm infants [36]. Milrinone use after ligation in preterm infants under 28 weeks of gestation may induce excessive vasodilation and thus require steroid administration.

This study had several limitations. First, we did not acquire laboratory data or imaging findings (e.g., echocardiography) to evaluate the severity of PDA and its complications. However, we adjusted for the clinical background, such as preoperative drug use and procedures, to estimate disease severity. Second, we did not describe the actual drug dosages administered to the patients because we were unable to obtain data on the partial use of drug vials, even though the effect of milrinone can differ by dosage [37]. Third, we were unable to obtain previous maternal histories, such as antenatal corticosteroid therapy, maternal comorbidities, and mode of delivery. Still, we believe that adjustment for neonatal conditions after birth can be a proxy for maternal factors. Fourth, we did not set outcomes specific to surgical complications as primary and secondary outcomes in this study. However, we consider that comprehensive outcomes are more relevant for patients and clinicians than outcomes specific to surgery. Therefore, the results of this study will be useful for patients and clinicians. Finally, we were unable to adjust for all confounding factors because this was a retrospective study. However, studies using real-world data can be superior to randomized controlled trials in the representativeness of the population and feasibility of conducting research [38]. Regardless of these limitations, this study is informative in that it describes the relevant outcomes of milrinone use after ligation in a large cohort. Further randomized controlled trials are warranted from the results of this study.

We performed propensity score analysis using a nationwide cohort of patients who underwent PDA ligation to compare the short-term outcomes of milrinone use after PDA ligation. Milrinone use after PDA ligation was not significantly associated with decreased in-hospital mortality or morbidities such as BPD, IVH, NEC, and ROP. Milrinone administration after PDA ligation may, therefore, be ineffective at improving short-term outcomes.

This study was approved by the Institutional Review Board of the University of Tokyo (approval number: 3501-[5] [May 19th, 2021]). The need for informed consent was waived by the Institutional Review Board of the University of Tokyo.

Takaaki Konishi received grants from Pfizer Co. Ltd., Kanzawa Medical Research Foundation, and Japan Kampo Medicines Manufacturers Association outside of the submitted work. There are no other conflicts of interest to disclose.

This work was supported by grants from the Ministry of Health, Labor, and Welfare, Japan (23AA2003 and 22AA2003 to H.Y.).

Hiroki Kitaoka and Takaaki Konishi conceptualized and designed the study, designed the data collection instruments, collected data, carried out the initial analyses, and drafted and revised the manuscript. Yoshihiko Shitara, Atsushi Ito, and Kohei Kashima conceptualized and designed the study and critically revised the manuscript. Asahi Fujita designed the data collection instruments, collected the data, carried out the initial analyses, and reviewed the manuscript. Hiroki Matsui coordinated and supervised data collection and reviewed the manuscript for important intellectual content. Motohiro Kato, Naoto Takahashi, and Hideo Yasunaga conceptualized and designed the study and critically reviewed and revised the manuscript for important intellectual content. All authors approved the final manuscript as submitted and agreed to be accountable for all aspects of the work.

1.
Hamrick SEG, Sallmon H, Rose AT, Porras D, Shelton EL, Reese J, et al. Patent ductus arteriosus of the preterm infant. Pediatrics. 2020 Nov;146(5):351–65.
2.
Schena F, Francescato G, Cappelleri A, Picciolli I, Mayer A, Mosca F, et al. Association between hemodynamically significant patent ductus arteriosus and bronchopulmonary dysplasia. J Pediatr. 2015;166(6):1488–92.
3.
Clyman RI, Wickremasinghe A, Merritt TA, Solomon T, McNamara P, Jain A, et al. Hypotension following patent ductus arteriosus ligation: the role of adrenal hormones. J Pediatr. 2014;164(6):1449–55.e1.
4.
Noori S, McCoy M, Friedlich P, Bright B, Gottipati V, Seri I, et al. Failure of ductus arteriosus closure is associated with increased mortality in preterm infants. Pediatrics. 2009 Jan;123(1):e138–144.
5.
Brooks JM, Travadi JN, Patole SK, Doherty DA, Simmer K. Is surgical ligation of patent ductus arteriosus necessary? The Western Australian experience of conservative management. Arch Dis Child Fetal Neonatal Ed. 2005;90(3):235–9.
6.
Clyman RI. Ibuprofen and patent ductus arteriosus. N Engl J Med. 2000;343(10):728–30.
7.
Bixler GM, Powers GC, Clark RH, Walker MW, Tolia VN. Changes in the diagnosis and management of patent ductus arteriosus from 2006 to 2015 in United States neonatal intensive care units. J Pediatr. 2017;189:105–12.
8.
Weisz DE, Mirea L, Rosenberg E, Jang M, Ly L, Church PT, et al. Association of patent ductus arteriosus ligation with death or neurodevelopmental impairment among extremely preterm infants. JAMA Pediatr. 2017;171(5):443–9.
9.
Noori S, Friedlich P, Seri I, Wong P. Changes in myocardial function and hemodynamics after ligation of the ductus arteriosus in preterm infants. J Pediatr. 2007;150(6):597–602.
10.
Harting MT, Blakely ML, Cox CS, Lantin-Hermoso R, Andrassy RJ, Lally KP. Acute hemodynamic decompensation following patent ductus arteriosus ligation in premature infants. J Invest Surg. 2008;21(3):133–8.
11.
Jain A, Sahni M, El-Khuffash A, Khadawardi E, Sehgal A, McNamara PJ. Use of targeted neonatal echocardiography to prevent postoperative cardiorespiratory instability after patent ductus arteriosus ligation. J Pediatr. 2012;160(4):584–9.e1.
12.
Teixeira LS, Shivananda SP, Stephens D, Van Arsdell G, McNamara PJ. Postoperative cardiorespiratory instability following ligation of the preterm ductus arteriosus is related to early need for intervention. J Perinatol. 2008;28(12):803–10.
13.
Ulrich TJB, Hansen TP, Reid KJ, Bingler MA, Olsen SL. Post-ligation cardiac syndrome is associated with increased morbidity in preterm infants. J Perinatol. 2018;38(5):537–42.
14.
Wang H, Jain A, Weisz DE, Moraes TJ. Trends in patent ductus arteriosus ligation in neonates and changes in outcomes: a 10-year multicenter experience. Pediatr Pulmonol. 2021;56(10):3250–7.
15.
Sehgal A, Francis JV, Lewis AI. Use of milrinone in the management of haemodynamic instability following duct ligation. Eur J Pediatr. 2011;170(1):115–9.
16.
Korkmaz L, Ozdemir A, Pamukçu Ö, Güneş T, Ozturk MA. Which inotropic drug, dobutamine or milrinone, is clinically more effective in the treatment of postligation cardiac syndrome in preterm infants?Am J Perinatol. 2022;39(2):204–15.
17.
El-Khuffash AF, Jain A, Weisz D, Mertens L, McNamara PJ. Assessment and treatment of post patent ductus arteriosus ligation syndrome. J Pediatr. 2014;165(1):46–52.e1.
18.
Ting JY, Resende M, More K, Nicholls D, Weisz DE, EL-Khuffash A, et al. Predictors of respiratory instability in neonates undergoing patient ductus arteriosus ligation after the introduction of targeted milrinone treatment. J Thorac Cardiovasc Surg. 2016;152(2):498–504.
19.
Hallik M, Ilmoja ML, Tasa T, Standing JF, Takkis K, Veigure R, et al. Population pharmacokinetics and dosing of milrinone after patent ductus arteriosus ligation in preterm infants. Pediatr Crit Care Med. 2019;20(7):621–9.
20.
Yasunaga H. Real world data in Japan: chapter II the diagnosis procedure combination database. Ann Clin Epidemiol. 2019;1(3):76–9.
21.
Konishi T, Fujiogi M, Shigemi D, Nishioka K, Matsui H, Fushimi K, et al. Surgical and obstetric outcomes of breast cancer surgery during pregnancy: a nationwide database study in Japan. Breast Cancer Res Treat. 2022 Oct;195(3):289–99.
22.
Konishi T, Yoshimoto T, Fujiogi M, Yamana H, Tanabe M, Seto Y, et al. Validity of operative information in Japanese administrative data: a chart review-based analysis of 1,221 cases at a single institution. Surg Today. 2022 Oct;52(10):1484–90.
23.
Ono S, Ishimaru M, Ida Y, Yamana H, Ono Y, Hoshi K, et al. Validity of diagnoses and procedures in Japanese dental claims data. BMC Health Serv Res. 2021;21(1):1116–0.
24.
Yamana H, Moriwaki M, Horiguchi H, Kodan M, Fushimi K, Yasunaga H. Validity of diagnoses, procedures, and laboratory data in Japanese administrative data. J Epidemiol. 2017;27(10):476–82.
25.
Shigemi D, Ishimaru M, Matsui H, Fushimi K, Yasunaga H. Suicide attempts during pregnancy and perinatal outcomes. J Psychiatr Res. 2021;133:101–5.
26.
Shigemi D, Isogai S, Uda K, Aso S, Matsui H, Fushimi K, et al. Association between rehabilitation during hospitalization and perinatal outcomes among pregnant women with threatened preterm birth. J Matern Neonatal Med. 2021;34(7):1028–33.
27.
Michihata N, Matsui H, Fushimi K, Yasunaga H. Hospital volume and mortality due to preterm patent ductus arteriosus. Pediatr Int. 2016;58(11):1171–5.
28.
Feudtner C, Feinstein JA, Zhong W, Hall M, Dai D. Pediatric complex chronic conditions classification system version 2: updated for ICD-10 and complex medical technology dependence and transplantation. BMC Pediatr. 2014;14(1):199–7.
29.
Isayama T. The clinical management and outcomes of extremely preterm infants in Japan: past, present, and future. Transl Pediatr. 2019;8(3):199–211.
30.
Yasunaga H. Introduction to applied statistics: chapter 1 propensity score analysis. Ann Clin Epidemiol. 2020;2(2):33–7.
31.
Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res. 2011;46(3):399–424.
32.
Desai RJ, Franklin JM. Alternative approaches for confounding adjustment in observational studies using weighting based on the propensity score: a primer for practitioners. BMJ. 2019;367:l5657.
33.
Thomas LE, Li F, Pencina MJ. Overlap weighting: a propensity score method that mimics attributes of a randomized clinical trial. JAMA. 2020;323(23):2417–8.
34.
El-Khuffash AF, Jain A, Dragulescu A, Mcnamara PJ, Mertens L. Acute changes in myocardial systolic function in preterm infants undergoing patent ductus arteriosus ligation: a tissue Doppler and myocardial deformation study. J Am Soc Echocardiogr. 2012;25(10):1058–67.
35.
Noori S, Seri I. Neonatal blood pressure support: the use of inotropes, lusitropes, and other vasopressor agents. Clin Perinatol. 2012;39(1):221–38.
36.
Htun ZT, Schulz EV, Desai RK, Marasch JL, McPherson CC, Mastrandrea LD, et al. Postnatal steroid management in preterm infants with evolving bronchopulmonary dysplasia. J Perinatol. 2021;41(8):1783–96.
37.
Hallik M, Tasa T, Starkopf J, Metsvaht T. Dosing of milrinone in preterm neonates to prevent postligation cardiac syndrome: simulation study suggests need for bolus infusion. Neonatology. 2017;111(1):8–11.
38.
Sherman RE, Anderson SA, Dal Pan GJ, Gray GW, Gross T, Hunter NL, et al. Real-world evidence: what is it and what can it tell us?N Engl J Med. 2016;375(23):2293–7.