key: cord-0813083-93486y6e authors: Magnusson, Amanda; Ahle, Margareta; Swolin‐Eide, Diana; Elfvin, Anders; Andersson, Roland E. title: Population‐based study showed that necrotising enterocolitis occurred in space–time clusters with a decreasing secular trend in Sweden date: 2017-04-24 journal: Acta Paediatr DOI: 10.1111/apa.13851 sha: 04e4c6f40520ee04e8ce8a8dc19bc874e4b411be doc_id: 813083 cord_uid: 93486y6e AIM: This study investigated space–time clustering of neonatal necrotising enterocolitis over three decades. METHODS: Space–time clustering analyses objects that are grouped by a specific place and time. The Knox test and Kulldorff's scan statistic were used to analyse space–time clusters in 808 children diagnosed with necrotising enterocolitis in a national cohort of 2 389 681 children born between 1987 and 2009 in Sweden. The municipality the mother lived in and the delivery hospital defined closeness in space and the time between when the cases were born – seven, 14 and 21 days – defined closeness in time. RESULTS: The Knox test showed no indication of space–time clustering at the residential level, but clear indications at the hospital level in all the time windows: seven days (p = 0.026), 14 days (p = 0.010) and 21 days (p = 0.004). Significant clustering at the hospital level was found during 1987–1997, but not during 1998–2009. Kulldorff's scan statistic found seven significant clusters at the hospital level. CONCLUSION: Space–time clustering was found at the hospital but not residential level, suggesting a contagious environmental effect after delivery, but not in the prenatal period. The decrease in clustering over time may reflect improved routines to minimise the risk of contagion between patients receiving neonatal care. Necrotising enterocolitis (NEC) is the most common gastrointestinal emergency among neonates, and it mainly affects preterm infants, with mortality rates ranging from 10% to 50%. The highest mortality rate is found among infants requiring surgery (1) (2) (3) (4) (5) . The overall incidence of NEC varies between studies, from 0.3 to 1.0 per 1000 live births (1, 6, 7) . However, in extremely preterm and very low birth weight infants, the incidence is approximately 7% (5, 8) . The pathogenesis of NEC is multifactorial, and there are some factors that remain unknown (4, 5) . Most cases of NEC occur sporadically. Nevertheless, reports of clusters or outbreaks suggest that an infectious element could be a causal factor in NEC (3, (9) (10) (11) (12) (13) . This hypothesis is supported by the fact that improvements in infection-control procedures have stopped outbreaks of NEC (11, 14) . Seasonal variations in the incidence of NEC have been described, which also indicate that an infectious agent may contribute to the clustering of the disease (6, 12, 15) . Several microbial organisms have been proposed as possible causes of NEC, for example Klebsiella pneumonia, Staphylococcus Aureus, Escherichia coli, Clostridium difficile, norovirus and rotavirus, but no specific causative organism was identified in some outbreaks (9, 11, 14, 16, 17) . It has also been suggested that overcrowding in neonatal intensive care units (NICUs) has contributed to clusters of NEC (14) . Nevertheless, the majority of reports describing outbreaks of NEC are retrospective and based on observed suspected outbreaks that could just be random. Abbreviations NEC, Necrotising enterocolitis; NICU, Neonatal intensive care unit. This study investigated space-time clustering of necrotising enterocolitis from 1987 to 2009 using national Swedish data on nearly 2.4 million births. Clustering was found at the hospital level during 1987-1997, but not during 1998-2009, and not at the residential level. The decrease in clustering over time could be related to enhanced routines to minimise the spread of any potential necrotising enterocolitis inducing contagion between patients in the neonatal intensive care unit. Furthermore, most of the described outbreaks have been on a hospital level, while clustering based on the mother's residential municipality has not been addressed (9, (11) (12) (13) . In reports on NEC outbreaks, the cluster concept tends to be used subjectively without a standard definition (18) . Our group previously presented a national, populationbased study on NEC epidemiology and trends in Sweden, which described an increase in the incidence of NEC between 1987 and 2009 (6). The same cohort was used in the present study to investigate space-time clusters of NEC on two levels for closeness in space: the mother's residential municipality and the delivery hospital. Furthermore, the present study was designed to examine whether there had been any change in the occurrence of space-time clusters over time, by studying two subperiods: 1987-1997 and 1998-2009 . A cohort of newborn infants with a diagnosis of NEC was identified from the following registers held by the Swedish National Board of Health and Welfare: the National Patient Register, the Swedish Medical Birth Register and the National Cause of Death Register. All children born between 1987 and 2009 in Sweden with a discharge diagnosis of NEC according to the 9th or 10th revision of the International Classification of Diseases -ICD-9 code 777F or ICD-10 code P77were identified. The NEC diagnosis was introduced to ICD-9 in 1987 and is based on the modified Bell NEC staging criteria (19, 20) . As it was not possible to identify the exact date for the NEC diagnosis, the date of birth of the study subjects was used for time comparisons in the cluster analysis. Further details about the identification process were previously described (6) . An anonymised extract covering the background population of all children born in Sweden during the same time period as the NEC cases was also obtained from the Birth Register. This extract contained perinatal information and demographic data, including the municipality the mother lived in and the delivery hospital. Sweden has a highly centralised care policy for very preterm and extremely preterm infants, based on intention to transfer mothers with a high risk of preterm delivery to a regional level three hospital before they give birth. As a result, most of the infants diagnosed with NEC are admitted to the NICU at the hospital in which they were born. Two methods were used to analyse for space-time interactions between NEC cases: the Knox space-time cluster analysis and Kulldorff's space-time permutation scan statistic (21, 22) . The Knox test is based on an analysis of the proximity in space and time of all possible n(n À 1)/2 distinct pairs of cases (23) . Each individual pair is classified into one of four cells in a 2 9 2 table, with distance (close/not close) and time (close/not close) on the two axes, according to whether the two parts are close or not close to each other in terms of geographical distance and time. A pair of cases is regarded as being in close proximity if their dates of birth are close and if their geographical locations at the time of birth are close. Closeness in the date of birth was divided into time windows of seven, 14 and 21 days apart. Two geographical levels were used to define closeness in space: the mother's residential municipality and the delivery hospital. The number of pairs of cases observed in close proximity was compared with the expected number of pairs, which was obtained from the cross-products of the column and row totals. If the observed number of pairs of cases exceeded the expected number of pairs, there was evidence of space-time clustering. The magnitude of the excess, or deficit, was estimated by calculating the strength of clustering using the equation S = [(O-E)/E] 9 100, where S was the strength, O was the number of pairs of cases observed and E was the expected number of pairs. To study any changes over time in NEC clustering, the population was divided into two cohorts according to the subjects' year of birth: 1987-1997 and 1998-2009. The Knox test was used to compare the two time periods, and the binomial test was used to compare the change in incidence of NEC in the two time periods. In addition, Kulldorff's scan statistic, based on a spacetime permutation model, was used to identify the presence of space-time and purely temporal clusters of cases (21). Kulldorff's scan statistic is based on the number of observed cases among all births that have taken place within a circle of varying radius in space in one dimension and in a time window with a varying duration in the other dimension. The statistic is centred at all geographical locations to look for possible clusters. Thus, the circular window is flexible in location, size and time. For the analyses of clustering on the residential level, we used the geographical coordinates of the centre of the mothers' residential municipality. For the analyses of clustering at each delivery hospital, we used a purely temporal scan statistic, with a time window of varying duration. The number of observed cases in a cluster was compared to what would have been expected if the spatial and temporal locations of all cases were independent of each other, so that there was no space-time interaction. As described by Kulldorff et al., the scan statistic makes minimal assumptions about the time, geographical location or size of the cluster and can be adjusted for both purely spatial and purely temporal variations (21) . The Poisson distribution was used for testing the statistical significance of the difference between the observed and expected number of pairs in the Knox test. Kulldorff's scan statistic was assessed by Monte Carlo hypothesis testing in 999 simulations, which meant that the smallest p value we could get was 0.001 (24) . Statistical significance was set at p < 0.05. The study used Stata Statistical Software, version 13 (StataCorp LP, College Station, TX, USA) and SaTScan, version 9.4.2 (Kulldorff M. and Information Management Services Inc., MA, USA) for the statistical analyses (25) . The study was approved by the Regional Ethical Review Board of Link€ oping (Dnr 2010/405-32). The study was based on a total of 2 389 681 births from 1987 to 2009, and the patient characteristics are described in Table 1 . Information about the delivery hospital and the mothers' residential municipality was missing for 5,621 and 19,130 children, respectively. We identified 808 cases of NEC, including 27 pairs of twins. Each twin pair with NEC was counted as one instance of NEC for the cluster analyses. Information about the mother's residential municipality and delivery hospital was missing for 12 and seven of the 808 cases, respectively. After we excluded the 27 second twins and the births with missing information on municipality or delivery hospital, there were 769 cases for the analyses based on municipality and 774 cases for the analyses based on delivery hospital. Due to the centralised care of preterm infants in Sweden, 422 of the 774 NEC cases (54%) occurred at a hospital that did not match the residential municipality of the mother. To be specific, 58% of all the NEC cases among extremely preterm births, with a gestational age under 28 weeks, and 22% of all the NEC cases among term births, with a gestational age over 36 weeks, occurred at a hospital that was not the closest to the mother's municipality. The cohort in the first time period, 1987-1997, consisted of 1 113 946 infants and 289 cases of NEC, resulting in an NEC incidence of 0.26 per 1000 live births. During the second time period, 1998-2009, the cohort consisted of 1 275 735 infants and 519 cases of NEC, giving an NEC incidence of 0.41 per 1000 live births. There was a significant increase in the incidence of NEC in the second time period compared to the first time period (p < 0.001) ( Table 1 ). The Knox test did not indicate any space-time clustering at a residential level in any of the studied time windows of seven, 14 or 21 days. There was a significant space-time clustering at a hospital level, with the strongest clustering at a time window of seven days (S = 36.3, p = 0.026) ( Table 2 ). The Knox test is sensitive to time-related shifts in the background population, which can give biased results. We therefore performed separate analyses for each of the two time periods. The first time period showed significant space-time clustering of NEC in the time windows of seven and 14 days, with the strongest clustering at seven days (S = 122.2, p = 0.003) ( Table 2) . During the second time period, there was no significant clustering at a hospital level in any of the studied time windows. Kulldorff's scan statistic At a residential level, Kulldorff's scan statistic only identified one single space-time cluster of four cases during 17 days in January 1990. The four cases came from four different municipalities within a radius of 34 kilometres. At a hospital level, the purely temporal cluster analysis identified seven instances of temporal clusters at seven different hospitals (Table 3 ). In four of the seven clusters identified by Kulldorff 0 s scan statistic, the cluster only consisted of two patients. However, several of these clusters occurred in hospitals with a low number of deliveries and few expected cases of NEC in the given time interval. Of the seven statistically significant clusters, five occurred during November to April and only two clusters occurred during May to October. The present study showed that NEC occurred in clusters at a hospital level, as found with both the Knox test and Kulldorff's scan statistic. When we compared two different time periods -1987-1997 and 1998-2009, using the Knox test, significant clustering was only found in the early time period and the strongest significance was found using seven days as the time window. Our results showed no signs of space-time clustering related to the mother's residential municipality with the Knox test and only one single cluster with Kulldorff 0 s scan statistic. Several possible explanations for clustering on a hospital level have previously been described. One explanation is that NEC is associated with a nosocomial infection spread from one child to another in a NICU. Hill et al. described an outbreak of NEC associated with Klebsiella pneumoniae in all cases at one NICU (26) . Han et al. and Alfa et al. described outbreaks of NEC associated with the Clostridium species (27, 28) . A second possible mechanism for clustering on a hospital level could be transmission from the healthcare workers to the infants, as suggested by Harbarth et al., who described an outbreak of Enterobacter cloacae during a period of overcrowding and understaffing in the NICU (29) . As the present study was a retrospective register study, no investigations could be carried out into whether the bacteria in the infants and among the staff contributed to the clusters. Contamination of human milk fortifier or formula is a third possibility for clustering on a hospital level. Van Acker et al. described an outbreak of NEC where the same bacteria were isolated from both the neonates with NEC and the powdered milk formula (10) . In Sweden, most infants receive human breast milk in NICUs, either from their mother or from a milk bank, but this milk is frequently enriched with human milk fortifier. A fourth possible explanation for clustering on a hospital level may be an accumulation of preterm births at referral hospitals due to referrals of at-risk pregnancies. This could theoretically lead to an overestimation of the number of clusters. The results from the Knox test showed significant clustering during 1987-1997, but not during 1998-2009, which did not support an overestimation of clusters due to centralised care, as the centralisation of neonatal intensive care in Sweden has increased over the last few decades. The finding of a decrease in clustering over time could be related to improvements in the neonatal intensive care of preterm infants. In this study, it was not possible to analyse whether the decrease in clusters was related to improved control of infection in the NICU, less overcrowding, better routines in the NICU or other reasons for reduced transmission of NEC between patients. Even though the Knox test showed no significant clustering of NEC during the last decade, Kulldorff's scan statistic indicated that clusters of NEC do still occur. Clustering on a residential level would, as described above, indicate that NEC is associated with causative agents, such as infections in the community. Stuart et al. described a strong association with the norovirus in an outbreak of NEC (11). Chany et al. showed a significant association between coronavirus infections and NEC (30) . These findings could indicate that the virus had its origin in the community and was then transmitted to the infants. The findings in the present study, in which the Knox test found no clustering on a residential level and Kulldorff's scan statistic found only one cluster, are strong indications against the theory that there is a connection between NEC and infections spread in the community. However, when studying the clusters on a hospital level with the Kulldorff 0 s scan statistic, it was noticed that the majority of the clusters at a hospital level occurred during November to April, which is also the season when most infections in the community occur. Our group has previously described this seasonal variation, with a peak in incidence of all cases of NEC in November and a decrease in May (6). The present study showed indications of space-time clustering of NEC on a hospital level in Sweden, but not at the level of the mother's residential municipality, suggesting a contagious environmental effect after delivery. The decrease in clustering on a hospital level over the last few decades may indicate that improved routines in modern neonatal care are effective in minimising the transfer of agents involved in the development of NEC between patients in the NICU. However, continued awareness of signs of clusters is still warranted to further minimise the risk of environmental factors for NEC being transferred from one patient to another. Necrotising enterocolitis hospitalisations among neonates in the United States Low birthweight, gestational age, need for surgical intervention and gramnegative bacteraemia predict intestinal failure following necrotising enterocolitis A cluster of necrotizing enterocolitis in term infants undergoing open heart surgery Necrotising enterocolitis Necrotizing enterocolitis Epidemiology and trends of necrotizing enterocolitis in Sweden Epidemiology of neonatal necrotising enterocolitis: a population-based study Variations in incidence of necrotizing enterocolitis in Canadian neonatal intensive care units Cluster of necrotizing enterocolitis in a neonatal intensive care unit: New Mexico Outbreak of necrotizing enterocolitis associated with Enterobacter sakazakii in powdered milk formula An outbreak of necrotizing enterocolitis associated with norovirus genotype GII.3 Epidemic occurrence of neonatal necrotizing enterocolitis Cluster of late preterm and term neonates with necrotizing enterocolitis symptomatology: descriptive and case-control study A decrease in the number of cases of necrotizing enterocolitis associated with the enhancement of infection prevention and control measures during a Staphylococcus aureus outbreak in a neonatal intensive care unit Seasonal variation in the incidence of necrotizing enterocolitis Nosocomial necrotising enterocolitis outbreaks: epidemiology and control measures Necrotising enterocolitis: is there a relationship to specific pathogens? Epidemiology of necrotizing enterocolitis temporal clustering in two neonatology practices Neonatal necrotizing enterocolitis. Therapeutic decisions based upon clinical staging Necrotizing enterocolitis in neonates fed human milk A space-time permutation scan statistic for disease outbreak detection The Knox method and other tests for space-time interaction The detection of space-time interactions Modified randomization tests for nonparametric hypotheses SaTScan is a trademark of Martin Kulldorff. The SaTScanTM software was developed under the joint auspices of Martin Kulldorff, the National Cancer Institute, and Farzad Nosocomial colonization with Klebsiella, type 26, in a neonatal intensive-care unit associated with an outbreak of sepsis, meningitis, and necrotizing enterocolitis An outbreak of necrotizing enterocolitis associated with a novel clostridium species in a neonatal intensive care unit An outbreak of Clostridium difficile necrotizing enterocolitis: a case for oral vancomycin therapy? Outbreak of Enterobacter cloacae related to understaffing, overcrowding, and poor hygiene practices Association of coronavirus infection with neonatal necrotizing enterocolitis We would like to thank Nils-Gunnar Pehrsson and Henrik Eriksson at Statistiska Konsultbyr an for their statistical expertise. This study was financed by grants from the ALF agreement between the Swedish government and county councils to Sahlgrenska University Hospital. The authors have no conflict of interests to declare.