Skip to main content

Multi-centre, randomised non-inferiority trial of early treatment versus expectant management of patent ductus arteriosus in preterm infants (the BeNeDuctus trial): statistical analysis plan

A Study protocol to this article was published on 04 August 2018



Controversy exists about the optimal management of a patent ductus arteriosus (PDA) in preterm infants. A persistent PDA is associated with neonatal mortality and morbidity, but causality remains unproven. Although both pharmacological and/or surgical treatment are effective in PDA closure, this has not resulted in an improved neonatal outcome. In most preterm infants, a PDA will eventually close spontaneously, hence PDA treatment potentially increases the risk of iatrogenic adverse effects. Therefore, expectant management is gaining interest, even in the absence of convincing evidence to support this strategy.


The BeNeDuctus trial is a multicentre, randomised, non-inferiority trial assessing early pharmacological treatment (24–72 h postnatal age) with ibuprofen versus expectant management of PDA in preterm infants in Europe. Preterm infants with a gestational age of less than 28 weeks and an echocardiographic-confirmed PDA with a transductal diameter of > 1.5 mm are randomly allocated to early pharmacological treatment with ibuprofen or expectant management after parental informed consent.

The primary outcome measure is the composite outcome of mortality, and/or necrotizing enterocolitis Bell stage ≥ IIa, and/or bronchopulmonary dysplasia, all established at a postmenstrual age of 36 weeks. Secondary short-term outcomes are comorbidity and adverse events assessed during hospitalization and long-term neurodevelopmental outcome assessed at a corrected age of 2 years. This statistical analysis plan focusses on the short-term outcome and is written and submitted without knowledge of the data.

Trial registration NTR5479. Registered on October 19, 2015, with the Dutch Trial Registry, sponsored by the United States National Library of Medicine NCT02884219 (registered May 2016) and the European Clinical Trials Database EudraCT 2017-001376-28.

Peer Review reports


A patent ductus arteriosus (PDA) is common in preterm infants [1] and is associated with neonatal mortality and morbidity, such as bronchopulmonary dysplasia (BPD), necrotizing enterocolitis (NEC) and intraventricular haemorrhage (IVH). Controversy exists about its optimal management, as early PDA treatment induces PDA closure but does not improve overall outcome [2,3,4].

Ductus arteriosus closure is delayed in prematurity, since it is not yet programmed for prompt postnatal closure [5]. This is supported by the observation of a high amount of spontaneous closure [6] even after ‘failed’ pharmacological treatment [7]. Since early pharmacological treatment has not been proven to improve outcome, it potentially increases the risk of iatrogenic adverse effect in patients in whom the PDA would have closed spontaneously. These observations have led to an increased interest in expectant PDA management [8]. Evidence to support expectant PDA management is scarce and conflicting [9], due to a high amount of open label treatment in placebo-controlled randomised controlled trials (RCTs) [4, 10] and a heterogeneous definition of (haemodynamic significant) PDA [11].

The BeNeDuctus trial assesses early pharmacological treatment within 24–72 h postnatal age (PNA) with ibuprofen versus expectant management of PDA in preterm infants in Europe. The study protocol was published previously [10]. This paper describes the statistical analysis plan (SAP) for short-term outcomes in detail, which is written and submitted without knowledge of the data.


The primary aim of the BeNeDuctus trial is to investigate whether an expectant management is non-inferior to early treatment (24–72 h PNA) for PDA (diameter > 1.5mm) in preterm infants, born at a gestational age less than 28 weeks, with regard to the composite outcome of mortality and/or NEC (Bell stage ≥ IIa) and/or BPD at a postmenstrual age (PMA) of 36 weeks.


Design and setting

The BeNeDuctus trial is a multicentre, non-inferiority trial in preterm infants with an echocardiographic-confirmed PDA > 1.5 mm before 72 h of PNA. Patients are recruited from 17 neonatal intensive care units (NICUs) in the Netherlands, Belgium and Denmark.

Study protocol development and conduct

The BeNeDuctus trial is registered with the Dutch Trial Register NTR5479 (registered on 19 October 2015), the registry sponsored by the United States National Library of Medicine NCT02884219 (registered May 2016) and the European Clinical Trials Database EudraCT 2017-001376-28. This study has been approved by the medical ethics committee of the Radboud University (CMO Arnhem-Nijmegen; Number 2016-2552/NL57885.091.16). The local ethics committee of each participating hospital approved the local feasibility of the study protocol.

During the course of the study the ethics committee approved two amendments. On 25 September 2017, the study protocol was adapted for use by participating sites outside The Netherlands, the data monitoring plan was updated and cyclo-oxygenase inhibitor (i.e. indomethacin and/or ibuprofen) was changed to ibuprofen due to unavailability of indomethacin. On 26 February 2018, reimbursement of investigational medicinal product in Belgium was added and the Data Safety Monitoring Board (DSMB) paragraph was edited according to an updated DSMB charter.

The trial is conducted according to the principles of the Declaration of Helsinki [12], good clinical practice guidelines and Dutch, Belgian and Danish legislation regarding medical research involving human subjects [13]. Informed consent from both parents needs to be obtained before patients can be randomised in the trial.

An independent DSMB monitors the progress of the study and performs interim-analysis on safety. The trial will be reported according to the Consolidated Standards of Reporting Trials (CONSORT) guidelines [14].

In- and exclusion criteria are described in detail in the previously published study protocol [10]. In summary, preterm infants born at a gestational age of less than 28 weeks with an echocardiographic-confirmed PDA with a diameter > 1.5 mm and transductal (predominantly) left-to-right shunting at a PNA of less than 72 h are eligible.

Included patients are randomised to either early pharmacological treatment with ibuprofen, dosage according to local protocol or expectant management. The use of ‘open-label’ pharmacological treatment in the expectant management group is most definitely discouraged and only allowed when predefined open-label treatment criteria are met [10]. In a strict sense open-label treatment is a misnomer, since in this non-placebo controlled nor blinded clinical trial, no label could be opened in the expectant management study group. Due to the widespread use of the term open label in PDA research, we will continue to use this term for pharmacological treatment in the expectant treatment group.

Randomisation and data collection

Eligible patients are randomly allocated to early pharmacological treatment or expectant management per centre and stratified according to gestational age in two strata (stratum A less than 260/7 weeks’ gestation or stratum B 260/7–276/7 weeks’ gestation). Stratification per centre was done for pragmatic reasons and is not likely to influence the primary outcome. Therefore, the analyses will not be adjusted for centre. In an additional analysis, the primary outcome will be adjusted for gestational age, as this is a prognostic covariate, using Poisson regression. Block size varies from four to eight. Multiple birth infants are preferably randomised independently, unless there is an explicit request from the parents to expose all siblings to the same treatment group as the first patient is allocated to. The randomisation process is centrally controlled and web-based.

Data management is implemented according to good clinical practice guidelines. Patient data, both until hospital discharge and long-term follow-up data, are entered via an electronic case record form in a central web-based Castor® Database to facilitate on-site data entry [15]. The most important data is shown in Table 1 for baseline characteristics and Table 2 for outcome measures.

Table 1 Baseline maternal and patient characteristics of randomised patients
Table 2 Outcome parameters per analysis

Primary outcome

The primary outcome is the dichotomous composite outcome of mortality, and/or NEC (Bell stage ≥ IIa) [16], and/or BPD at a PMA of 36 weeks. BPD is defined as the need for supplemental oxygen or positive pressure ventilatory support at a PMA of 36 weeks and diagnosed following international standard criteria by Bancalari [17], including an oxygen reduction test according to Walsh [18] (Table 2).

Short-term secondary outcomes

Secondary outcomes are adverse events till discharge, which are shown in Table 2.

Short-term secondary outcome definitions

Hyperlactataemia is defined as a serum concentration > 2.5 mmol/L, either capillary or arterial [19]. Renal failure is defined as creatinine > 120 μmol/L or urine output < 0.5 ml/kg/h. Hypertension is defined as mean arterial blood pressure > 95th percentile, according to Zubrow or Dionne [20, 21].

Respiratory support is divided in invasive, encompassing conventional and high-frequency ventilation, and non-invasive, encompassing non-invasive positive pressure ventilation, nasal continuous positive airway pressure and high-flow nasal cannula therapy. High-flow nasal cannula therapy is defined as a flow > 1 l per minute. Low flow, defined as flow < 1 l per minute, is considered as respiratory support for supplemental oxygen. Pulmonary air leakage is defined as pneumothorax, pulmonary interstitial emphysema or pneumomediastinum.

IVH is graded according to Volpe [22]. Periventricular echogenicity is graded according to Hashimoto [23].

Time to full enteral feeding is defined as the first moment full enteral feeding is reached, even in cases in which enteral feeding needs to be stopped afterwards, for example due to NEC.

Sepsis is defined as positive blood culture and need for antibiotics. Meningitis is defined as sepsis with antibiotic treatment regimen (dosage and duration) for meningitis, irrespective of spinal fluid culture. Thrombocytopenia is defined as a platelet count < 100 x 109/L. Hyperglycaemia is defined as elevated glucose level necessitating insulin treatment according to local protocol.

Retinopathy of prematurity (ROP) is graded according to the international classification [24].

Long-term secondary outcomes

Neurodevelopmental outcome is assessed in all Dutch and Belgian children in the National Neonatal Follow Up Program at a corrected age of 24 months by (a) paediatric and neurologic examination; (b) cognitive assessment with Bayley Scales of Infant and Toddler Development, Third Dutch Edition (BSID-III-NL); (c) behavioural assessment with Child Behaviour Check List (CBCL), Teacher Report Form (TRF) questionnaire. Despite the Movement Assessment Battery for Children, Second Dutch Edition (Movement ABC 2-NL) was mentioned in the study protocol, this will not be included in the neurodevelopmental outcome at 24 months corrected age as it is assessed at five years corrected age. Motor scales of the BSID-III-NL will be used to assess motor function. Equivalent assessments may be used for non-Dutch or Belgian children.

Long-term secondary outcome definitions

Participants who are alive and have a BSID III cognitive and motor composite score at 2 years corrected age ≥ 85 (−1 standard deviation) will be defined as having survived without neurodevelopmental impairment.

General health is considered (a) good, i.e. no readmissions in first 2 years after discharge home; (b) moderate, i.e. 1–3 readmissions; or (c) poor, i.e. >3 readmissions.

Neurologic examination is subdivided in (a) normal, i.e. no symptoms of pathology; (b) mild abnormal, i.e. mild muscle tone abnormalities and coordination problems that do not hamper functioning and/or development; and (c) abnormal, i.e. muscle tone abnormalities and/or coordination problems that are hampering development, cerebral palsy and/or infantile encephalopathy.

Vision is scored as (a) normal, i.e. no problems with sight; (b) mildly abnormal, i.e. treated by ophthalmologist/orthoptist for abbreviation (goggles) or strabismus/amblyopia; (c) abnormal, i.e. limited vision, but ability to see anything; and (d) severely abnormal, i.e. blind.

Hearing is scored as (a) normal, i.e. no hearing problems; (b) mildly abnormal, i.e. light hearing loss for which control or treatment; (c) abnormal, i.e. neurosensory hearing loss (partially) corrected with hearing aid; and (d) severely abnormal, i.e. neurosensory hearing loss, deafness.

Adverse events (AEs) and serious adverse events (SAE)

As this trial is undertaken in a fragile patient group, context-specific (S)AEs are not individually reported to the medical ethics committee, as described in the study protocol [10, 25, 26]. Non-context specific (S)AEs are reported to the medical ethics committee and will be listed.

Statistical methods specified in the study protocol

Sample size calculation

As described in the study protocol, the a priori risk of mortality, NEC and BPD is estimated to be 20%, 10% and 15%, respectively, for preterm infants born at a gestational age of less than 28 weeks, based on data from The Netherlands Perinatal Registry (PRN foundation) in the period 2008–2012. Non-inferiority is defined as a significant difference in the primary outcome parameter between the two groups of less than 10%. In other words, the one-side 95% confidence interval of the observed difference between an expectative approach and ibuprofen treatment should not exceed the non-inferiority margin of 10%.

With an estimated a priori risk for the composite of mortality and/or NEC and/or BPD at 36 weeks PMA of 35%, a one-sided type I error of 5% and a power of 80%, the sample size to exclude a non-inferiority margin of 10% for the difference of proportion of participants reaching the primary outcome parameter is 564 patients, being 282 patients in each group. This sample size is calculated using PASS 2008, version 08.0.8 NCSS.

Original proposed analyses

The originally proposed analyses are described in the published study protocol [10]. For the primary outcome, a 95% one-sided confidence interval for the risk difference will be calculated, and when based on this interval, a difference of 10% or more can be excluded, non-inferiority will be concluded.

Treatment effects for the dichotomous clinical outcomes will be reported using risk differences with 95% confidence interval. Normally distributed data will be presented as mean ± standard deviations, uneven distributed data as medians with interquartile ranges. Categorical data will be analysed using the chi-square for two and multiway tables. Continuous data will be analysed using the Student’s t test. Both intention-to-treat and per-protocol analyses will be employed. Statistical significance is defined as a p value < 0.05. Here, we present our final and further detailed SAP.

Interim analyses and safety reporting

To protect patients and to assist and advise the principal investigator in protecting the safety, validity and credibility of the trial, a DSMB is installed. The DSMB consists of a neonatologist, paediatric cardiologist and biostatistician. The members have no competing interests and are not involved in the trial. Their composition, tasks, responsibilities and working procedures are described in a charter (available on request).

Safety analysis is scheduled when 15%, 30%, 50% and 75% of the data are gathered. Before every interim analysis the DSMB receives a report from the trial coordinator and data manager, that includes data on (1) context-specific (S)AEs, subdivided in (a) hemodynamic, (b) pulmonary, (c) central nervous system, (d) gastrointestinal, (e) infectious, (f) renal, (g) ophthalmologic, (h) metabolic and (i) mortality; (2) non-context-specific (S)AEs; and (3) open-label treatment in expectant study group.

The DSMB charter states that there are two possible reasons for stopping the trial early, namely concerns for safety and futility. In principle, the trial will not be stopped early before the minimum number of evaluable patients required (n = 564) are included for beneficial effect of ibuprofen treatment on the primary outcome. The interim analyses are not associated with alpha spending.

Statistical analysis plan

Overall principles

As described above the first phase of data analysis and reporting will include all outcome data up to discharge to home. Analyses will start once all data to discharge of the last included patient have been obtained, the database has been checked and locked and this SAP has been submitted for publication. The analyses will be performed by TH, WdB and RD.

Intention-to-treat analysis will be performed. Since open-label treatment might have an effect on both PDA closure and might be associated with adverse events, analyses will also be done in a per-protocol population as described in detail below. For all relevant parameters, outcome with 95% confidence intervals will be presented. A p value < 0.05 will be considered statistically significant. For the primary analysis, no adjustment for multiple testing will be applied. The secondary analyses are to be considered as exploratory. For statistical programming and analysis, we will use SPSS version 25.0 (IBM Corp. 2017) and the R environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria).

The short-term outcome analyses will be performed and published before the assessment of the long-term outcomes at 2 years of corrected age of all patients have been completed, since it is considered unethical by the neonatal community to withhold the results of the primary and short-term secondary outcomes for another 2 years. In the second phase, the neurodevelopmental outcomes after 2 years will be analysed, reported and published afterwards.


As stated in the study protocol [10], the investigator or attending physician can decide to withdraw a subject from the study for urgent medical reasons. If they wish, parents or caregivers can withdraw their consent for the study at any time for any reason. Patients in the expectant management arm that meet the open-label criteria and patients in the early pharmacological treatment arm that meet the surgical ligation criteria will remain in follow up and are therefore not withdrawn from the study, as are patients that develop contraindications for continuation of ibuprofen (Fig. 2). Other reasons for withdrawal will be presented as supplementary material with reason and count per treatment group.

Handling of missing data

To minimize the amount of missing data, multiple attempts will be undertaken to retrieve data. Both level III and II (referral) hospitals will be contacted, since most patients will be transferred to referral hospitals after the initial hospitalization. As the primary outcome will be assessed at a PMA of 36 weeks, so generally before discharge to home, we anticipate minimal missing values. Missing data till discharge will therefore not be imputed.

For the long-term outcome at 2 years of corrected age, a higher amount of missing data is expected due to loss to follow up. In case of loss to follow-up the reason for this will be recorded. Their available trial results and characteristics will be analysed and sensitivity analyses will be performed as described later.

Definition of analysis sets

Intention-to-treat population

The intention-to-treat population consists of all randomised infants, regardless of protocol deviations or use of open-label pharmacological treatment or surgical ligation (Table 3). This population includes randomised patients who died before the first ibuprofen dosage was administered.

Table 3 Definition of population analysis sets

Per-protocol population

In the per-protocol population, only patients included and treated in accordance with the study protocol will be included. Patients included in the expectant management group receiving pharmacological treatment after fulfilling the open-label criteria will be included in the per-protocol expectant management group. Patients in whom other protocol deviations occurred will be excluded from the per-protocol analysis (Table 3).

Statistical analyses

Patient flow

The flow of trial participants is displayed in the CONSORT flow diagram (Fig. 1). As the echocardiographic screening for the presence of a PDA is common practice in several of the study sites, while in other centres this is considered a study procedure, patient flow differs. Therefore, informed consent is needed before an echocardiography can be performed as study procedure, while in case echocardiographic screening is common practice informed consent can be obtained after echocardiographic assessment for eligibility, as is shown in Fig.1. Reasons why patients are not eligible, and why eligible patients are not included will be summarized. Patient characteristics will be collected from all eligible infants that are not included in this study in order to assess any potential recruitment bias.

Fig. 1
figure 1

Consolidated Standards of Reporting Trials (CONSORT) 2010 flow diagram

The compliance to allocated treatment will be reported for all randomised patients in an additional flow chart, with reasons why allocated treatment was not given (Fig. 2).

Fig. 2
figure 2

Patient flow of received treatment per treatment allocation group

Protocol deviations

Protocol deviations are defined as deviations in eligibility criteria or open-label pharmacological treatment in case the predefined criteria, as described in the study protocol, were not fulfilled. All protocol deviations will be listed as supplementary material on publication.

Baseline characteristics

We will present the following maternal characteristics: (a) maternal age; (b) maternal obstetrical disease, divided in (1) pre-eclampsia; (2) haemolysis, elevated liver enzymes and low platelets (HELLP) syndrome; (3) placental abruption; (4) preterm premature rupture of membranes (PPROM); and (5) clinical chorioamnionitis, defined as maternal fever without other cause or prenatal antibiotics use without other source of infection than chorioamnionitis or PPROM; (c) maternal medication, especially non-steroidal inflammatory drugs (i.e. ibuprofen and/or indomethacin) and magnesium sulphate; (d) tocolysis; (e) administration of antenatal corticosteroids, with incomplete administration defined as < 48 h or > 14 days before birth; (f) type of delivery, either vaginally or caesarean section; and (g) multiple birth (Table 1).

The following baseline characteristics of patients will be presented: (a) gestational age at birth; (b) birth weight; (c) outborn; (d) gender; (e) Apgar scores at five minutes; (f) resuscitation after birth, subdivided in (1) non-invasive respiratory support; (2) invasive respiratory support; and (3) circulatory support, defined as the need for chest compressions and/or epinephrine; (g) respiratory distress syndrome; (h) surfactant administration, and; (i) PDA diameter. All variables will be presented as summary statistics according to allocation group of the trial (Table 1).

Continuous normally distributed variables will be summarized with mean and standard deviation. Continuous non-normally distributed variables will be summarized using median and interquartile range. Categorical variables will be summarized using counts and percentages.

Primary outcome

Crude estimates of the relative risk and absolute risk difference of the primary outcome in the expectant management group compared with the early pharmacological treatment group will be calculated. The main analysis of the effect of expectant management on the primary outcome will be performed using 95% one sided confidence interval for the risk difference. Differences in baseline characteristics are not expected due to randomisation. In case of an unexpected disbalance in baseline characteristics correction will be made with appropriate statistical testing (i.e. logistic regression) if possible.

Short-term secondary outcome

Short-term secondary outcomes will be analysed using a generalized linear model. Secondary outcomes will be presented in tables with counts and percentages per treatment group and presented as absolute risk differences. No formal adjustments or multiple comparisons will be made.

Long-term secondary outcome at 2 years corrected age

The key long-term outcome is neurodevelopmental outcome at 2 years corrected age.

Subgroup analyses

We will perform pre-specified exploratory subgroup analyses for the primary outcome by examining treatment and subgroup interaction effects in logistic regression models. We will perform five analyses, each examining one subgroup: (a) gestational age groups (less than or greater or equal to 260/7 weeks), (b) gender (male/female), (c) multiple pregnancy (yes/no), (d) antenatal corticosteroids (complete versus incomplete or no antenatal corticosteroids); and (e) birth weight (less than or greater or equal to 1000 g).

Sensitivity analyses

For BPD classification, there is uncertainty on how to categorize continuous positive airway pressure or high-flow nasal cannula with low or no supplemental oxygen [27]. The international criteria used will categorize those infants as having severe BPD [17]. In very premature infants, the reason for this support can be an impaired control of breathing rather than chronic parenchymal lung damage. Therefore, infants supported by continuous positive airway pressure or high-flow nasal cannula in room air will be classified as having mild BPD in an auxiliary sensitivity analysis in line with the STOP BPD trial [28].

If possible, we will perform sensitivity analyses on the primary outcome to investigate the impact of correlation between infant outcomes within twin or higher-order multiple births using generalized estimating equations with a logit link function.

For long-term neurodevelopmental outcome, sensitivity analyses will be performed (a) using complete cases only and by applying (b) best-case and (c) worst-case scenarios for the unobserved neurodevelopmental impairment outcome data.

Differences between the study protocol and statistical analysis plan

In accordance with the previously published study protocol [10], we further elaborate on the rationale of the analyses that will be performed and planned subgroup as well as sensitivity analyses.

Trial status

Initially, 15 participating centres started recruitment between December 2016 and April 2019, two additional centres started recruitment between August 2019 and May 2020. In another centre no patient was found to be eligible.

During the study we experienced lower inclusion rates than expected, both due to a lower rate of parental informed consent, but also a lower rate of PDA > 1.5mm than anticipated. As the funding by ZonMw stopped on 15 December 2020 and the conflicting DOXA Trial ( NCT04430790) started recruitment, enrollment was ended when in total 273 (48.4% of projected inclusion) patients were randomized (early pharmacological treatment n=137; expectant management n=136). The power calculation was based on data from 2008-2012, which might not represent current data. Even with the actual number of inclusions the BeNeDuctus Trial is the 3rd largest RCT investigating PDA treatment since 2000 with regard to the study population. As the first two largest studies investigated prophylactic treatment [29, 30], this study is in fact the largest RCT investigating early pharmacological treatment.

The last patient is included on 5 December 2020. The last assessment of the primary outcome will be completed in March 2021. The last follow-up at 2 years corrected age is expected around March 2023. The database for the short-term outcome measurements will be locked only after data has been monitored and this SAP has been submitted for publication in a peer-reviewed journal. The long-term outcome measurements will be locked mid-2023 after monitoring of the data.

Availability of data and materials

The results will be presented at scientific meetings and published in peer reviewed medical journals. The data that support the findings of this study are available from the corresponding author upon reasonable request. There will be an embargo on the data for two to five years.



Bronchopulmonary dysplasia


Bayley scales of infant and toddler development, third Dutch edition


Child behaviour checklist


Consolidated standards of reporting trials


Data safety monitoring board


Intraventricular haemorrhage

Movement ABC-2-NL:

Movement assessment battery for children, second Dutch edition


Necrotizing enterocolitis


Neonatal intensive care unit(s)


Patent ductus arteriosus


Postmenstrual age


Postnatal age


Retinopathy of prematurity


Statistical analysis plan


Teacher report form


  1. Lee JA, Kim MJ, Oh S, Choi BM. Current Status of Therapeutic Strategies for Patent Ductus Arteriosus in Very-Low-Birth-Weight Infants in Korea. J Korean Med Sci. 2015;30(Suppl 1):S59–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Evans P, O'Reilly D, Flyer JN, Soll R, Mitra S. Indomethacin for symptomatic patent ductus arteriosus in preterm infants. Cochrane Database Syst Rev. 2021;1:Cd013133.

    PubMed  Google Scholar 

  3. Mitra S, Florez ID, Tamayo ME, Mbuagbaw L, Vanniyasingam T, Veroniki AA, 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(12):1221–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Jansen EJS, Hundscheid T, Onland W, Kooi EMW, Andriessen P, de Boode WP. Factors Associated With Benefit of Treatment of Patent Ductus Arteriosus in Preterm Infants: A Systematic Review and Meta-Analysis. Frontiers in Pediatrics. 2021;9(45):626262.

    Article  Google Scholar 

  5. Hundscheid T, van den Broek M, van der Lee R, de Boode WP. Understanding the pathobiology in patent ductus arteriosus in prematurity-beyond prostaglandins and oxygen. Pediatr Res. 2019;86(1):28–38.

    Article  CAS  PubMed  Google Scholar 

  6. Semberova J, Sirc J, Miletin J, Kucera J, Berka I, Sebkova S, et al. Spontaneous Closure of Patent Ductus Arteriosus in Infants </=1500 g. Pediatrics. 2017;140(2):e20164258.

    Article  Google Scholar 

  7. Weber SC, Weiss K, Buhrer C, Hansmann G, Koehne P, Sallmon H. Natural History of Patent Ductus Arteriosus in Very Low Birth Weight Infants after Discharge. J Pediatr. 2015;167(5):1149–51.

    Article  PubMed  Google Scholar 

  8. Ngo S, Profit J, Gould JB, Lee HC. Trends in Patent Ductus Arteriosus Diagnosis and Management for Very Low Birth Weight Infants. Pediatrics. 2017;139(4):e20162390.

    Article  Google Scholar 

  9. Hundscheid T, Jansen EJS, Onland W, Kooi EMW, Andriessen P, de Boode WP. Conservative Management of Patent Ductus Arteriosus in Preterm Infants—A Systematic Review and Meta-Analyses Assessing Differences in Outcome Measures Between Randomized Controlled Trials and Cohort Studies. Frontiers in Pediatrics. 2021;9(44):626261.

    Article  Google Scholar 

  10. Hundscheid T, Onland W, van Overmeire B, Dijk P, van Kaam A, Dijkman KP, et al. Early treatment versus expectative management of patent ductus arteriosus in preterm infants: a multicentre, randomised, non-inferiority trial in Europe (BeNeDuctus trial). BMC Pediatr. 2018;18(1):262.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Zonnenberg I, de Waal K. The definition of a haemodynamic significant duct in randomized controlled trials: a systematic literature review. Acta Paediatr. 2012;101(3):247–51.

    Article  PubMed  Google Scholar 

  12. World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. JAMA. 2013;310(20):2191–4.

  13. Central Committee on Research Involving Human Subjects (CCMO). Standard research file [Available from: Accessed 10 Aug 2021.

  14. Kessler KM. The CONSORT statement: explanation and elaboration. Consolidated Standards of Reporting Trials. Ann Intern Med. 2002;136(12):926–7; author reply -7.

    Article  Google Scholar 

  15. Castor EDC. Castor Electronic Data Capture 2019 [2019]. Available from:

  16. Neu J. Necrotizing enterocolitis: the search for a unifying pathogenic theory leading to prevention. Pediatr Clin North Am. 1996;43(2):409–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Bancalari E, Claure N. Definitions and diagnostic criteria for bronchopulmonary dysplasia. Semin Perinatol. 2006;30(4):164–70.

    Article  PubMed  Google Scholar 

  18. Walsh MC, Yao Q, Gettner P, Hale E, Collins M, Hensman A, et al. Impact of a physiologic definition on bronchopulmonary dysplasia rates. Pediatrics. 2004;114(5):1305–11.

    Article  PubMed  Google Scholar 

  19. Fauchère JC, Bauschatz AS, Arlettaz R, Zimmermann-Bär U, Bucher HU. Agreement between capillary and arterial lactate in the newborn. Acta Paediatr. 2002;91(1):78–81.

    Article  PubMed  Google Scholar 

  20. Dionne JM, Abitbol CL, Flynn JT. Hypertension in infancy: diagnosis, management and outcome. Pediatr Nephrol. 2012;27(1):17–32.

    Article  PubMed  Google Scholar 

  21. Zubrow AB, Hulman S, Kushner H, Falkner B. Determinants of blood pressure in infants admitted to neonatal intensive care units: a prospective multicenter study. Philadelphia Neonatal Blood Pressure Study Group. J Perinatol. 1995;15(6):470–9.

    CAS  PubMed  Google Scholar 

  22. Volpe JJ. Neurology of the newborn. 5th ed. Philadelphia, Pa.; London: Saunders; 2008.

    Google Scholar 

  23. Hashimoto K, Hasegawa H, Kida Y, Takeuchi Y. Correlation between neuroimaging and neurological outcome in periventricular leukomalacia: diagnostic criteria. Pediatr Int. 2001;43(3):240–5.

    Article  CAS  PubMed  Google Scholar 

  24. The International Classification of Retinopathy of Prematurity revisited. Arch Ophthalmol. 2005;123(7):991–9.

  25. Cook D, Lauzier F, Rocha MG, Sayles MJ, Finfer S. Serious adverse events in academic critical care research. CMAJ. 2008;178(9):1181–4.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Cook D, Moore-Cox A, Xavier D, Lauzier F, Roberts I. Randomized trials in vulnerable populations. Clinical trials. 2008;5(1):61–9.

    Article  PubMed  Google Scholar 

  27. Poindexter BB, Feng R, Schmidt B, Aschner JL, Ballard RA, Hamvas A, et al. Comparisons and Limitations of Current Definitions of Bronchopulmonary Dysplasia for the Prematurity and Respiratory Outcomes Program. Ann Am Thorac Soc. 2015;12(12):1822–30.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Onland W, Merkus MP, Nuytemans DH, Jansen-van der Weide MC, Holman R, van Kaam AH. Systemic Hydrocortisone To Prevent Bronchopulmonary Dysplasia in preterm infants (the SToP-BPD study): statistical analysis plan. Trials. 2018;19(1):–178.

  29. Schmidt B, Davis P, Moddemann D, Ohlsson A, Roberts RS, Saigal S, et al. Long-term effects of indomethacin prophylaxis in extremely-low-birth-weight infants. N Engl J Med. 2001;344(26):1966–72.

    Article  CAS  PubMed  Google Scholar 

  30. Van Overmeire B, Allegaert K, Casaer A, Debauche C, Decaluwe W, Jespers A, et al. Prophylactic ibuprofen in premature infants: a multicentre, randomised, double-blind, placebo-controlled trial. Lancet. 2004;364(9449):1945–9.

    Article  CAS  PubMed  Google Scholar 

Download references


We would like to thank the members of the Data Safety Monitor Board of the BeNeDuctus trial: Dr. Anne van Kempen, neonatologist; Dr. Caroline van Baal, biostatistician; and Hans Breur, pediatric cardiologist.

We would like to thank Karin Deckers and Wendy Jansen (research coordinators) for their invaluable support.

BeNeDuctus trial study group members:

Peter H. Dijk, MD, PhD, Department of Neonatology, University Medical Centre Groningen, Beatrix Children’s Hospital, Groningen, The Netherlands;

Anton H.L.C. van Kaam, MD, PhD, Emma Children's Hospital Amsterdam University Medical Centers, Department of Neonatology, University of Amsterdam, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands;

Tessa de Baat, MD, Emma Children's Hospital Amsterdam University Medical Centers, Department of Neonatology, University of Amsterdam, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands;

Koen P. Dijkman, MD, Department of Neonatology, Maxima Medical Centre Veldhoven, Veldhoven, The Netherlands;

Eduardo Villamor, MD, PhD, Department of Neonatology, Maastricht University Medical Centre, Maastricht, The Netherlands;

André A. Kroon, MD, PhD, Department of Neonatology, Erasmus MC Sophia Children’s Hospital, Rotterdam, The Netherlands;

Remco Visser, MD, PhD, Department of Neonatology, Leiden University Medical Centre, Willem Alexander Children's Hospital, Leiden, The Netherlands;

Susanne M. de Tollenaer, MD, PhD, Department of Neonatology, Isala Women’s and Children’s Hospital Zwolle, Zwolle, The Netherlands;

Filip Cools, MD, PhD, Department of Neonatology, UZ Brussel – Vrije Universiteit Brussel, Brussels, Belgium;

Marisse Meeus, MD, Department of Neonatology, Antwerp University Hospital, Edegem, Belgium;

Anne-Britt Johansson, MD, Department of Neonatology, Hôpital Universitaire des Enfants Reine Fabiola, Bruxelles/Brussels, Belgium;

Frank Derriks, MD, Department of Neonatology, Cliniques Universitaires de Bruxelles, Hôpital Erasme, Brussels, Belgium;

Catheline Hocq, MD, Department of Neonatology, Cliniques Universitaires St Luc, Brussels, Belgium;

Alexandra Zecic, MD, Department of Neonatology, Ghent University Hospital, Gent, Belgium;

Tine Brink Henriksen, MD, PhD, Neonatal Intensive Care Unit, Department of Pediatrics and Adolescent Medicine, Aarhus University Hospital, Aarhus, Denmark;

Kasper Jacobsen Kyng, MD, PhD, Neonatal Intensive Care Unit, Department of Pediatrics and Adolescent Medicine, Aarhus University Hospital, Aarhus, Denmark.


This trial is funded by ZonMw – The Netherlands Organization for Health Research & Development (project number 843002622).

Author information

Authors and Affiliations




TH, WdB and RD were involved in the conception and design of the statistical analysis plan and drafted the manuscript. WO, DN, DV, AM, EK and WdV are members of the BeNeDuctus steering committee. All individuals mentioned in the Acknowledgements are involved in the BeNeDuctus trial study as local investigator at the participating centres. All authors read, edited and approved the final manuscript.

Corresponding author

Correspondence to Tim Hundscheid.

Ethics declarations

Ethics approval and consent to participate

This study has been approved by the MEC of the Radboud University (CMO Arnhem-Nijmegen; Number 2016-2552/NL57885.091.16). Informed consent has been obtained from the parents of all participants in the trial.

Consent for publication

Not applicable.

Competing interests

Nothing to declare

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hundscheid, T., Donders, R., Onland, W. et al. Multi-centre, randomised non-inferiority trial of early treatment versus expectant management of patent ductus arteriosus in preterm infants (the BeNeDuctus trial): statistical analysis plan. Trials 22, 627 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: