Article Text

Download PDFPDF

Association between perinatal interventional activity and 2-year outcome of Swiss extremely preterm born infants: a population-based cohort study
  1. Mark Adams1,2,
  2. Thomas M Berger3,
  3. Cristina Borradori-Tolsa4,
  4. Myriam Bickle-Graz5,
  5. Sebastian Grunt6,
  6. Roland Gerull7,
  7. Dirk Bassler1,
  8. Giancarlo Natalucci1
  9. For the Swiss Neonatal Network & Follow-Up Group
    1. 1 Department of Neonatology, Universitätsspital Zürich, Zürich, Switzerland
    2. 2 Epidemiology, Biostatistics and Prevention Institute, University Zürich, Zürich, Schweiz, Switzerland
    3. 3 Neonatal and Paediatric Intensive Care Unit, Kinderspital Luzern, Luzern, Switzerland
    4. 4 Department of Paediatrics, Geneva University Hospital, Geneva, Switzerland
    5. 5 Department Woman-Mother-Child, University Hospital Lausanne, Lausanne, Switzerland
    6. 6 Division of Neuropaediatrics, Development and Rehabilitation, Children’s University Hospital, Inselspital, University of Bern, Bern, Switzerland
    7. 7 Department of Neonatology, University of Basel Children’s Hospital (UKBB), Basel, Switzerland
    1. Correspondence to Dr Mark Adams; mark.adams{at}


    Objectives To investigate if centre-specific levels of perinatal interventional activity were associated with neonatal and neurodevelopmental outcome at 2 years of age in two separately analysed cohorts of infants: cohort A born at 22–25 and cohort B born at 26–27 gestational weeks, respectively.

    Design Geographically defined, retrospective cohort study.

    Setting All nine level III perinatal centres (neonatal intensive care units and affiliated obstetrical services) in Switzerland.

    Patients All live-born infants in Switzerland in 2006–2013 below 28 gestational weeks, excluding infants with major congenital malformation.

    Outcome measures Outcomes at 2 years corrected for prematurity were mortality, survival with any major neonatal morbidity and with severe-to-moderate neurodevelopmental impairment (NDI).

    Results Cohort A associated birth in a centre with high perinatal activity with low mortality adjusted OR (aOR 0.22; 95% CI 0.16 to 0.32), while no association was observed with survival with major morbidity (aOR 0.74; 95% CI 0.46 to 1.19) and with NDI (aOR 0.97; 95% CI 0.46 to 2.02). Median age at death (8 vs 4 days) and length of stay (100 vs 73 days) were higher in high than in low activity centres. The results for cohort B mirrored those for cohort A.

    Conclusions Centres with high perinatal activity in Switzerland have a significantly lower risk for mortality while having comparable outcomes among survivors. This confirms the results of other studies but in a geographically defined area applying a more restrictive approach to initiation of perinatal intensive care than previous studies. The study adds that infants up to 28 weeks benefited from a higher perinatal activity and why further research is required to better estimate the added burden on children who ultimately do not survive.

    • neonatology
    • fetal medicine
    • quality in health care

    This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See:

    Statistics from

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Strengths and limitations of this study

    • Geographically defined cohort study comparing 2-year outcome of extremely preterm infants between hospitals with high and low perinatal interventional activity based on a quantitative score of 3 obstetric and 4 neonatal indicators.

    • Crude and risk-adjusted ORs for outcome between high and low activity centres were calculated using regression models and measuring predictive validity of risk adjustment.

    • Missing outcome data of surviving infants lost to follow-up were addressed by multiple imputation and the validity of data imputation was tested by sensitivity analysis using non-imputed data.

    • The inter-relatedness of some of the perinatal activity score components (eg, surfactant and respiratory support) form a limitation. However, our aim was to quantify activity and not measure quality of care. As such, the components quantify at best the level of proactive care by focusing on the whole repertoire of first perinatal supportive interventions.

    • Further minor limitations are represented by the composite nature of the outcome measures ‘any major morbidity’ and ‘neurodevelopmental impairment’, which could have led to a loss of information.


    Over recent decades, progress in perinatal care has led to vastly improved survival rates for infants born at extremely low gestational age (GA). Consequently, the GA limit for initiating life-sustaining intensive care at birth was lowered to 22 weeks, below which treatment is generally not recommended.1 This, however, raised ethical concerns regarding the added burden to infants who ultimately do not survive and the increased risk of neurosensory impairment among those that do survive.2 There is little consensus about the policy of treating these infants. Instead, hospital practices regarding the initiation of intensive care have dramatically influenced rates of survival and survival without impairment.3 This has led to large outcome variation among hospitals. In 2015, Rysavy et al quantified the effect of perinatal interventional activity on outcome between a group of 24 centres of the Eunice Kennedy Shriver National Institute of Child Health and Human  Development (NICHD) network in the USA. Among infants born at 22–23 weeks gestation, centre rates of active treatment accounted for 78% of the variation in survival between centres. For those born at 24 weeks gestation, it accounted for 22% of the variation.4 There is a similar variation in decision-making approaches in Europe, which appears to reflect local philosophy and practice rather than economic or demographic markers, individual variations in health expenditure or rates of preterm birth.5

    Swiss obstetricians and neonatologists are known to be generally restrictive with regard to initiation of intensive care for infants at the limit of viability. A recent study revealed that the decision-making regarding initiation of treatment was based almost exclusively on GA and that most infants born <24 weeks gestation received a priori comfort care.6 Nevertheless, centre-to-centre outcome variability has repeatedly been reported and is assumed to be in part associated with the local ethical decision-making.7–9 What remained unclear was whether the lower mortality in some centres leads to higher risk of neurodevelopmental impairment (NDI) in survivors.

    This study aimed to quantify the outcome differences between centres with a higher propensity towards providing active perinatal care at birth from those that provide more a priori comfort care. We assessed infants in cohort A (22+0 to 25+6 weeks GA) and compared their outcome variability with that in cohort B (26+0 to 27+6 weeks GA) where we expected no difference in either perinatal care provision or outcome based on previous studies.4 10–13


    Study population and procedure

    Geographically defined, retrospective cohort study including all live-born infants between 22+0 and 27+6 weeks of gestation in Switzerland from 2006 to 2013. We excluded infants born with a major congenital malformation, defined as being the primary cause of death or requiring surgery or a chromosomal anomaly. We also excluded four infants born and treated exclusively in three step-down units. We extracted electronically recorded data (challenged for plausibility and completeness) from the prospective national database of the Swiss Neonatal Network & Follow-up Group (SNN). Data from delivery room deaths were audited by a researcher visiting the hospitals’ maternity wards.

    Participating centres were obliged to inform parents about the scientific use of anonymised data.

    Neonatal data

    Perinatal and neonatal variables were defined as follows: GA as the best estimate available based on prenatal ultrasound examination during the first trimester of pregnancy. Birthweight z-scores, intraventricular haemorrhage grade 3 or higher (IVH), cystic periventricular leukomalacia (cPVL), bronchopulmonary dysplasia (BPD), severe retinopathy of prematurity stages 3–4 (ROP), necrotising enterocolitis (NEC), late-onset neonatal sepsis (LOS) and socioeconomic status (SES) were defined as previously published.14

    Neurodevelopmental outcome at 2 years of age

    Experienced neurodevelopmental or developmental paediatricians performed the standardised Swiss follow-up assessment15 at 18–24 months corrected age in one of the Swiss follow-up centres. Between 2006 and 2012, examinations were based on Bayley Scales of Infant Development, 2 edition (BSID-II)16 and afterwards on Bayley Scales of Infant and Toddler Development, 3 edition (Bayley-III).17 Up until 2006, some of the children were tested using the Griffiths Mental Development Scales-Revised (GMDS)18 during the transition period to BSID-II. According to recent literature,19–22 Bayley-III scores below 85 (-1SD) were considered as equivalent of indices below 70 (-2SD) in the BSID-II and the GMDS. The GMDS development quotient was considered equivalent to BSID-II Mental Developmental Index.23 Infants with significant disability precluding completion of the development test were assigned a development score of 1 below −3SD. Cerebral palsy was defined according to Rosenbaum et al 24 and was graded according to the Gross Motor Function Classification System (GMFCS) for children aged ≤2 years.25 Vision and hearing were assessed by either direct examination or caregiver report.

    Perinatal interventional activity score

    In order to differentiate between Swiss centres with a higher propensity to initiate intensive care for infants born below 26 weeks gestation, we calculated a perinatal activity score per centre based on the model presented by Serenius et al.10 Obstetric and neonatal activity scores reflecting the intensity of care in each centre were calculated on the basis of the rates of three key obstetric indicators (delivery at level III hospitals, that is, inborn; complete course of antenatal steroids; caesarean section), and four key neonatal indicators (surfactant within 2 hours after birth; any of the following activities in the delivery room: bag/mask ventilation, continuous positive airway pressure, endotracheal intubation, epinephrine or epinephrine supply, cardiac compression; no recorded decision for primary non-intervention in the delivery room (available for all infants); infants admitted for intensive care (out of infants alive at 30 min after birth)). The mean obstetric and neonatal activity scores comprised the perinatal activity score (henceforth, activity score). Each step included normalisation by assigning the centre with the highest rate for each indicator a score of 100; the remaining centres were assigned proportional scores. We calculated an activity score for infants born between 22–25 weeks gestation (cohort A) and a separate one for infants born 26–27 weeks gestation (cohort B). The first score was used to classify the nine Swiss perinatal centres into a ‘low’ and a ‘high’ activity group using the arbitrary threshold of an activity score of 80 in order to generate two groups of similar size with clear difference in perinatal activity. To allow comparison between cohort A and cohort B, cohort B was split into the same centres with low and high activity as for cohort A. To confirm the dichotomous association of perinatal activity with outcome, we also tested for a linear association. For this, we generated a variable that increases by one for every 5-point increase in perinatal activity using the two perinatal activity scores for cohort A and B.


    Outcome parameters were defined as follows: mortality at 2 years corrected age, death or moderate-to-severe NDI, death or major morbidity, survival with moderate-to-severe NDI and survival with major morbidity. Moderate-to-severe NDI was defined as either mental or motor development index below 70 (-2SD) in the BSID-II; cognitive or motor composite score below 85 (see below) in the Bayley-III; a global score of the GMDS below 70 (-2SD); cerebral palsy with GMFCS above 1; the absence of useful hearing even with aids (ie, >90 dB hearing level); blindness or only perception of light. Major morbidities encompassed IVH, cPVL, BPD, ROP, NEC or LOS.

    Data completeness

    Patient population coverage was assessed by comparison with the birth registry of the Swiss Federal Statistical Office and yielded 86.4% of all live births between 22 and 27 weeks gestation. 1.2% of the datasets were missing information on ‘major morbidity’ and were therefore eliminated in the non-imputed analyses including this outcome. Other data were missing in less than 0.1% of cases. As 2-year follow-up data were missing in 11% of all surviving infants between 22 and 25 weeks gestation and 17% between of 26 and 27 weeks gestation, we performed a fivefold imputation with chained equations for the cohorts A (all infants), A (survivors), B (all infants) and B (survivors), respectively, allowing the calculation of pooled adjusted ORs (aORs) for the above-listed outcomes.26 Data for imputation were determined as missing at random. As a sensitivity analysis, we performed the same tests with crude, that is, non-imputed data. As outborn infants were by concept unevenly distributed between high and low activity centres, they cannot be assumed to be missing at random for the imputation. We, therefore, performed another sensitivity analysis excluding all outborn infants for both the activity score as well as for the regression analysis.

    Statistical analysis

    We performed multivariable-adjusted logistic regression to compare outcome between centres with lower and higher perinatal activity. Adjustment was made for GA and GA2 (as the combination better models the non-linear dependency of most neonatal outcomes on GA and therefore results in risk adjustment with higher predictive validity), birthweight z-score, male sex, multiple births and SES. To estimate the validity of the adjustment, we calculated how well the above listed variables predicted outcome. If this predicted validity is low, the adjustment made has limited explanatory power.27 We calculated the area under receiver operating characteristics curve (AUC). AUC values between 0.7 and 0.8 were considered to represent moderate, and >0.8 to represent high predictive validity, respectively. All statistical analyses were performed using R V.3.4.28


    Study population and neonatal data

    A total of 2063 infants were born alive in Switzerland between 2006 and 2013 below 28 weeks gestation (figure 1). A total of 1839 of those infants were registered by SNN between 2007 and 2013. They correspond to 1900 life births registered at the Swiss Federal Statistical Office in the same time period which results in a population coverage of 97% for the study period (reference data for 2006 were not available as GA was not included in the national register prior to 2007). Of the 2063 infants, we excluded four infants because they were born in and treated at three step-down units and 85 infants because of major congenital malformations. Among the 1974 included infants, 712 (36%) died before discharge and 6 (0.3%) after discharge but before reaching 2 years of age. In cohorts A and B, 319 (89% of survivors) and 743 (83%) were assessed at 2 years corrected age, respectively. One hundred ninety-four (15%) children were lost to follow-up (109 refused follow-up, 57 could not be reached/moved away, 28 unknown loss to follow-up). In both cohorts, difference between baseline characteristics of children with and without follow-up was negligible (table 1), particularly concerning the factor ‘activity score’ which determines the comparison groups for the subsequent analysis. Outborn children were less frequently assessed in both cohorts.

    Figure 1

    Study population overview. GA, gestational age; NICU, neonatal intensive care unit.

    Table 1

    Baseline characteristics cohorts A/B: children with and without 2-year follow-up (FU)

    Perinatal interventional activity score

    Table 2 displays the proportion of infants receiving either of the treatments determining the obstetrical activity (items1–3) and the neonatal activity (items4–7) for cohort A in each centre. Activity scores were determined out of the mean of the obstetric and neonatal activity scores and ranged from 57 to 98. Four centres with an activity score <80 were classified as ‘low activity centres’ encompassing 461 infants in cohort A and 554 in cohort B. The remaining five centres were classified as ‘high activity centres’ and combined 466 infants in cohort A and 493 in cohort B. A separate centre-specific activity score was determined for cohort B using the same routine (online supplementary table 1). Each 5-point increase in perinatal activity score reduced the risk for mortality within 2 years of life for cohort A (aOR 0.74, 95% CI 0.68 to 0.80) and cohort B (aOR 0.82, 95% CI 0.70 to 0.97).

    Supplementary file 1

    Table 2

    Centre activity grouping cohort A (standardised ratios)

    After risk adjustment, the aOR in cohort A for mortality (0.22, 95% CI 0.15 to 0.31) and death or NDI (0.30, 95% CI 0.21 to 0.42) were significantly lower (figure 2) for high compared with low activity centres. After restricting the cohort to survivors, aOR for any major neonatal morbidities was 0.75 (95% CI 0.47 to 1.19) and for moderate-to-severe NDI 0.95 (95% CI 0.46 to 1.98). The crude baseline characteristics of low versus high activity centres in cohort A (table 3) were comparable for size, GA range, birthweight z-score, multiple births and SES. However, high activity centres had more male infants and less outborn infants. The crude outcome proportions reflect the aOR results. Repeating the same evaluation using non-imputed data resulted in almost the same point estimates and CIs (dotted lines in figure 2). As outborn infants were unevenly distributed between children with and without follow-up, we repeated the analysis excluding all 75 outborn infants from the beginning for both the activity score and the regression analysis. The resulting aOR was almost equivalent to the ones reported above (online supplementary figure 1). Excluding the delivery room deaths, median neonatal intensive care unit (NICU) age at death in high activity centres was higher (8 vs 4 days) as was length of stay (100 vs 73 days). As can be seen in figure 3, the distribution of age at death was different between high and low activity centres. In high activity centres, the highest density was reached above the median, whereas it was below the median in low activity centres. The violin plot thereby allows taking the full distribution of the data into account rather than relying on median and IQR alone.

    Figure 2

    ORs of outcome when born in a unit with high perinatal interventional activity. ORs with 95% CI of outcome when born in high activity units, crude or adjusted for GA, GA2, birthweight z-score, male sex, multiple births and socioeconomic status. Dotted lines represent results received when using non-imputed data. AUC, area under receiver operator characteristic curve representing predictive validity of risk adjustment. GA, gestational age; NDI, neurodevelopmental impairment.

    Figure 3

    Age at death in NICU (days).  Log-transformed age at death in NICU (days) for low and high activity centres in cohort A and cohort B. Crude y-axis scale on right-hand side. Violin plot: box-plot with median, IQR and density (width=frequency). GA, gestational age; NICU, neonatal intensive care unit.

    Table 3

    Low versus high activity centres cohorts A and B

    Centres classified as having ‘high activity’ for cohort A, still had a higher perinatal activity score (91 vs 84) in cohort B with matching outcome variability. aOR for mortality (0.49, 95% CI 0.34 to 0.71) and death or NDI (0.72, 95% CI 0.53 to 0.97) were significantly lower in high activity centres (figure 2). In contrast to cohort A, the aOR for any major morbidities was also significantly lower (0.63, 95% CI 0.48 to 0.84) whereas moderate-to-severe NDI outcome in survivors was higher without reaching significance (1.17, 95% CI 0.77 to 1.78). Again, crude baseline characteristics (table 3) were comparable with more male and fewer outborn infants in high activity centres. Crude outcomes mirror the aOR result. In cohort B, median age at death in NICU in high activity centres was again higher (11 vs 6 days), whereas length of stay was only somewhat longer (85 vs 78 days).

    Neurodevelopmental outcome at 2 years of age

    Follow-up assessment was performed at a mean (SD) corrected age of 22.6 (3.1) months (table 4). Of 1062 examined children, 729 (69%), 265 (25%) and 46 (4%) were assessed with the Bayley-II, Bayley-III and GMDS, respectively, while in 19 (2%) children only a part of a developmental assessment could be performed because of severe cerebral palsy or poor compliance. A total of 220 (21%) infants suffered from moderate-to-severe NDI and 842 (79%) infants showed favourable outcome.

    Table 4

    Outcome at 2 years of age (corrected)


    Swiss perinatal centres with a higher score of perinatal activity for infants at the borderline of viability (cohort A) had significantly lower aOR for mortality and for the composite outcome death or NDI. After restricting the cohort to survivors at 2 years of age corrected for prematurity, there was no significant difference between centres with high and low intensity regarding risk for major morbidity or NDI. Cohort B mirrored the outcome differences between high and low activity perinatal centres determined for cohort A, that is, significantly lower aOR for mortality and for death or NDI. In contrast to cohort A, the lower mortality in high activity centres in cohort B was based exclusively on NICU mortality. High activity centres in cohort B had a significantly lower aOR for any major neonatal morbidity. Except for the aOR for mortality, death or NDI or any major morbidity in cohort A, the predictive validity measured by AUC was low to negligible. For those aORs, we would expect the true value to lie between the crude and the adjusted level, which confirm the results.29

    The association between perinatal interventional activity and outcome has been reported before. We based our centre-specific activity score on the model of the Swedish regional activity score published by Serenius et al.10 In their study, live-born infants between 22 and 24 weeks GA treated in regions with high versus low activity also had lower aOR for mortality (0.43, 95% CI 0.26 to 0.73) and for death or NDI (aOR 0.48, 95% CI 0.27 to 0.84) at 2.5 years corrected age. As in our study, survivors in high activity regions did not have a higher risk for any NDI (aOR 0.63, 95% CI 0.31 to 1.28). This reflects the difficulty in assessing outcome at 2 years in general and using NDI as a measure of efficacy.30 31 Other reports between intensity of care and survival date back as far as 2004 when Håkansson et al for the first time reported how a proactive perinatal strategy increases the number of live births and improves the infant’s postnatal condition in Sweden.32 Rysavy et al were also able to show a clear association between the outcome of infants receiving active care and those that did not for infants between 22 and 24 weeks GA, but not for those between 25 and 26 weeks GA.4 A study determining the perinatal outcome for extremely preterm infants in relation to place of birth in England reported lower odds for mortality when they were born in a level 3 centre with higher perinatal activity based on staffing and activity data collected by questionnaire.11 In 2013, an NICHD study could show that ‘centre intervention rates’ significantly predicted mortality rates for infants <25 weeks but not for infants ≥25 weeks gestation.12 Centre intervention rates included parameters similar or identical to the ones chosen for the calculation of the activity score in this study. Using yet another version of a similar intervention score, the MOSAIC study of 2009 (Models of OrganiSing Access to Intensive Care for very preterm births) was able to show an association between the score and mortality for infants <26 weeks but not for infants at 26–27 weeks gestation.13 Swiss centres rarely initiated intensive care below 24 weeks gestation during the study period.6 8 This is why we chose the same GA ranges for cohorts A and B as the MOSAIC study but not as the other studies. Yet, even with the more restrictive approach, we confirm the results of previous studies in infants <26 weeks gestation.

    Swiss centres classified as having ‘high activity’ for cohort A still had a higher perinatal activity score in cohort B. In contrast to the studies reported above,4 10–13 the resulting outcome variability in cohort B is congruent to the variability in cohort A. As variability in mortality in cohort B was exclusively observed in NICU but not in delivery room-mortality, it is difficult to determine an association to either higher perinatal activity or quality of care, or a combination of both. Arguments exist for each: centres of the high activity group also have a higher perinatal activity score in cohort B and the continuous perinatal activity score reduces risk of mortality per each 5-point increase in activity score. As variability in mortality in cohort B exists exclusively for NICU and not for delivery room mortality, variability in the decision-making process regarding redirection of care based on futility may exist. On the other hand, centres with a high perinatal activity also had a significantly lower aOR for any major morbidities which would favour the assumption of difference in quality of care. This would be in line with the observations by Smith et al whose findings suggest that the approach taken with infants at the limits of viability is associated with improved outcomes of more-mature infants.33

    As reported in other studies,10 34 median age at death was higher in high activity centres. Serenius et al argue that the longer duration of life allows for a trial of life, whereas Costeloe et al see the reason in different principle causes of death with a tendency towards later complications of prematurity in the case of the group of infants dying later.10 34  Figure 3 displays that the age at death for both cohorts are comparable in their IQR but, together with the longer median length of stay for all infants, indicate an added burden to infants who ultimately do not survive in high activity centres.


    A 5-point increase in perinatal interventional activity significantly decreases the aOR for mortality of infants born 22–25 weeks gestation in Switzerland. High activity centres have a significantly lower aOR for mortality and death or NDI and outcome among survivors at 2 years of age is comparable. The effect observed was repeated for infants between 26 and 27 weeks gestation, indicating that all extremely preterm infants benefited. Age at death and duration of hospitalisation are higher suggesting a higher burden for patients who ultimately do not survive in high activity centres. Although the results favour the high activity approach, further research is required to better estimate the accompanying added burden on the children. Outcome data at 5 years of age or later are urgently required to confirm these findings with a more predictive long-term assessment.


    We gratefully thank all children and their parents who participated in this study. We also thank the following members of the Swiss Neonatal Network & Follow-Up Group.


    1. 1.
    2. 2.
    3. 3.
    4. 4.
    5. 5.
    6. 6.
    7. 7.
    8. 8.
    9. 9.
    10. 10.
    11. 11.
    12. 12.
    13. 13.
    14. 14.
    15. 15.
    16. 16.
    17. 17.
    18. 18.
    19. 19.
    20. 20.
    21. 21.
    22. 22.
    23. 23.
    24. 24.
    25. 25.
    26. 26.
    27. 27.
    28. 28.
    29. 29.
    30. 30.
    31. 31.
    32. 32.
    33. 33.
    34. 34.


    • Contributors MA performed the biostatistical evaluation and together with GN was involved in study design, data collection, analysis and interpretation of the results and wrote, reviewed and revised the manuscript. TMB, DB, CB-T, MB-G, SG and RG were involved in data collection, study design, analysis and interpretation of data and reviewed the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

    • Funding MA receives a salary as network coordinator for the Swiss Neonatal Network. GN received financial support by the Swiss National Science Foundation; grant PZOOP3_161146.

    • Disclaimer The authors received no specific funding for this work. The sponsors were not involved in the following matters: study design; collection, analysis, and interpretation of data; writing of the report and the decision to submit the manuscript for publication.

    • Competing interests MA receives a salary as network coordinator for the Swiss Neonatal Network.

    • Ethics approval Data collection and evaluation for this study has been approved by the institutional ethical review board (KEK-ZH-Nr 2014–0552).

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data sharing statement Additional data are available by emailing

    • Collaborators Swiss Neonatal Network & Follow-up Group: Aarau: Cantonal Hospital Aarau, Children’s Clinic, Department of Neonatology (Ph. Meyer, C. Anderegg), Department of Neuropaediatrics (A. Capone Mori, D. Kaeppeli); Basel: University of Basel Children’s Hospital (UKBB), Department of Neonatology (S. Schulzke), Department of Neuropaediatrics and Developmental Medicine (P. Weber); Bellinzona: San Giovanni Hospital, Department of Paediatrics (G.P. Ramelli, B. Goeggel Simonetti); Berne: University Hospital Berne, Department of Neonatology (M. Nelle), Department of Paediatrics (B. Wagner), Department of Neuropaediatrics (M. Steinlin, S. Grunt); Biel: Children’s Hospital Wildermeth, Department of Pediatrics (M. Gebauer), Development and Paediatric Neurorehabilitation Center (R. Hassink); Chur: Children’s Hospital Chur, Department of Neonatology (W. Bär), Department of Neuropaediatrics (E. Keller, Ch. Killer); Fribourg: Cantonal Hospital Fribourg, Department of Neuropaediatrics (K. Fuhrer); Geneva: Department of child and adolescent, University Hospital (HUG), Neonatology Units (R. E. Pfister), Division of Development and Growth (P. S. Hüppi, C. Borradori-Tolsa); Lausanne: University Hospital (CHUV), Department of Neonatology (J.-F. Tolsa, M. Roth-Kleiner), Department of Child Development (M. Bickle-Graz); Lucerne: Children’s Hospital of Lucerne, Neonatal and Paediatric Intensive Care Unit (T. M. Berger), Department of Neuropaediatrics (T. Schmitt-Mechelke, F. Bauder); Lugano: Regional Hospital Lugano, Department of Paediatrics (V. Pezzoli); Muensterlingen: Cantonal Hospital Muensterlingen, Department of Paediatrics (B. Erkert, A. Mueller); Neuchatel: Cantonal Hospital Neuchatel, Department of Paediatrics (M. Ecoffey); St. Gallen: Cantonal Hospital St. Gallen, Department of Neonatology (A. Malzacher), Children’s Hospital St. Gallen, Neonatal and Paediatric Intensive Care Unit (J. P. Micallef), Department of Child Development (A. Lang-Dullenkopf); Winterthur: Cantonal Hospital Winterthur, Department of Neonatology (L. Hegi), Social Paediatrics Center (M. von Rhein); Zurich: University Hospital Zurich (USZ), Department of Neonatology (D. Bassler, R. Arlettaz), University Children’s Hospital Zurich, Department of Neonatology (V. Bernet) and Child Development Center (B. Latal, G. Natalucci); Zollikerberg: Hospital Zollikerberg, Neonatology Clinic (M. Moenkhoff).

    • Patient consent for publication Not required.