The use of outbreak information in the interpretation of clustering of reported cases of Escherichia coli O157 in space and time in Alberta, Canada, 2000-2002

Temporal, spatial and space-time distribution of infections caused by five major enteric pathogens, Ontario, Canada, 2010-2017

AIMS

In Canada, enteric diseases pose substantial health and economic burdens. The distribution of these diseases is uneven across both geography and time and understanding these patterns is therefore important for the prevention of future outbreaks. We evaluated temporal, spatial and space-time clustering of laboratory-confirmed cases of Campylobacter spp. (n = 28,728), non-typhoidal Salmonella spp. (n = 22,640), Shiga toxin-producing Escherichia coli (STEC; n = 1340), Yersinia spp. (n = 1674) and Listeria monocytogenes (n = 471) infections, reported between 2010 and 2017 inclusive in Ontario, the most populous province in Canada (population ~ 13,500,000 in 2016).

METHODS AND RESULTS

For each enteric pathogen, we calculated the mean incidence rates (IRs) for Ontario's 35 public health unit (PHU) areas and visualized them using choropleth maps. We identified temporal, spatial and space-time high infection rate clusters using retrospective Poisson scan statistics. Campylobacter and Salmonella infections had the highest IRs, while Listeria infections had the lowest. Campylobacter, Salmonella, STEC and Listeria mostly clustered temporally in the spring/summer and sometimes extended into fall, while Yersinia showed a less clear seasonal pattern. The IR visualizations and spatial and space-time scan statistics showed geographic heterogeneity of infection rates with high infection rate clusters detected mainly in PHUs across the southwestern and central-western regions of Ontario for Campylobacter, Salmonella and STEC infections, and mainly in PHUs located in the central-eastern regions for Yersinia and Listeria. A high proportion of cases in some of the significant Salmonella, STEC and Listeria infection clusters were linked to disease outbreaks.

CONCLUSIONS

Results from this study will inform heightened public health surveillance, and prevention and control programmes, in populations and regions of high infection rates. Further research is needed to determine the pathogen-specific socioeconomic, environmental and agricultural risk factors that may be related to the temporal and spatial disease patterns we observed in our study.

Incidence, Demographic, and Seasonal Risk Factors of Infections Caused by Five Major Enteric Pathogens, Ontario, Canada, 2010-2017

In Canada, enteric infections cause significant health and economic burden. We evaluated the individual characteristics of laboratory-confirmed cases of Campylobacter spp. (n = 28,728), non-typhoidal Salmonella spp. (n = 22,640), Yersinia spp. (n = 1674), Verotoxin-producing Escherichia coli (VTEC; n = 1340), and Listeria monocytogenes (n = 471), reported between 2010 and 2017 inclusive, in Ontario, Canada (population ∼13,500,000). We calculated overall and pathogen-specific annual and mean incidence rates (IRs) for Ontario. We used multivariable Poisson and negative binomial regression models to estimate incidence rate ratios (IRRs) for years, seasons, age groups, and sexes, and we included two-way age and sex interaction terms in the models. Campylobacter and Salmonella infections had the highest IRs whereas Listeria infections had the lowest IRs. None of the infections showed long-term trends over the 8-year study period; however, rates of all five infections were elevated in the summer. More Salmonella, VTEC, and Listeria infections were linked to disease outbreaks than were Campylobacter and Yersinia infections. Overall, mean IRs of Campylobacter, Salmonella, Yersinia, and VTEC infections were highest in children 0-4 years old, whereas Listeria IRs peaked in adults 60 years and older. Higher mean IRs of Campylobacter were observed in males. No other differences by sex were statistically significant. The same mean rate was observed in both sexes for Listeria. Adjusting for all other factors, significant age- and sex-specific differences in IRs were observed in Campylobacter, Salmonella, and VTEC infection rates. No significant interactions of age and sex were found for Yersinia and Listeria infections. Future research should focus on the pathogen-specific socioeconomic, environmental, or agricultural risk factors that might be responsible for these infections.

Detection of spatial and spatio-temporal Salmonella Heidelberg and Salmonella Typhimurium human case clusters focused around licensed abattoirs in Ontario in 2015, and their potential relation to known outbreaks

Salmonellosis is one of several zoonotic diseases for which individuals with occupational animal contact, including abattoir workers, are at an increased risk. If meat is contaminated during slaughter, this can increase the risk of enteric illness for consumers. In this study, we investigated whether reported cases of Salmonella Heidelberg and Typhimurium were clustered around abattoirs in Ontario in 2015 and whether there was any evidence (laboratory/exposure) to suggest an abattoir at the centre of a cluster might be the source of exposure. Data for each reported case of S. Heidelberg and S. Typhimurium in Ontario in 2015 were collected. Multi-focused and non-focused spatial and space-time cluster detection tests were performed for each serotype, with and without cases linked to known outbreaks, using Poisson and space-time permutation models. Focused tests included the location of abattoirs operational in all or part of 2015. Laboratory data and exposure information were used to explore the relatedness of cases within identified clusters. Focused spatial tests identified clusters of S. Heidelberg and S. Typhimurium around abattoirs. Focused space-time permutation tests identified 2 significant space-time clusters of S. Heidelberg; one cluster (n = 11 cases) included 8 of 9 cases associated with a known outbreak and the other cluster (n = 18 cases) was not part of a previously identified outbreak. Review of laboratory and risk factor information suggested that cases within each cluster shared a common exposure. Cases were not asked about goat or sheep meat consumption. The focused cluster test, particularly with the space-time permutation model, could assist in identifying outbreaks associated with a particular physical location, such as an abattoir. Improvements to the current case investigation process, such as consistent collection and reporting of high-risk occupation information and more detailed food consumption history, could assist in outbreak identification when coupled with this statistic.

Area-Level Clustering of Shiga Toxin-Producing Escherichia coli Infections and Their Socioeconomic and Demographic Factors in Ontario, Canada: An Ecological Study

Shiga toxin-producing Escherichia coli (STEC) infections are an important health burden for human populations in Ontario and worldwide. We assessed 452 STEC cases that were reported to Ontario's reportable disease surveillance system between 2015 and 2017. A retrospective scan statistic using a Poisson model was used to detect high-rate STEC clusters at the forward sortation area (FSA; the first three digits of a postal code) level. A significant spatial cluster in the southwest region of Ontario was identified. A case-case logistic regression analysis was applied to compare FSA-level socioeconomic and demographic characteristics among STEC cases included inside the spatial cluster with cases outside of the cluster. Cases included in the spatial cluster had higher odds of living in FSAs with a low median family income, low proportion of lone-parent families, and low proportion of the visible minority population. In addition, STEC cases inside the cluster had higher odds of coming from rural FSAs. Our study demonstrated that STEC cases were spatially clustered in Ontario and their clustering was associated with FSA-level socioeconomic and demographic determinants of cases.

Detection of spatial, temporal and space-time Salmonella Heidelberg and Salmonella Typhimurium clusters in Ontario in 2015, and comparisons to known outbreaks

PURPOSE

Salmonellosis is one of several reportable diseases in Ontario (ON). Two or more cases of the same serotype that are linked to a common exposure or related to one another in time and/or space are considered a potential outbreak. While laboratory data can help to determine the molecular relatedness of cases, results may take up to several weeks. This study aimed to assess the utility of the retrospective spatial scan statistic in detecting clusters of Salmonella Heidelberg and Salmonella Typhimurium cases using data from ON in 2015. Identified clusters were validated by laboratory data (where available) to determine whether identified clusters were likely outbreaks.

METHODS

Data representing the location of each reported S. Heidelberg or S. Typhimurium case in 2015, responsible serotype and symptom onset date were exported to SaTScan for retrospective spatial, temporal, and space-time analyses using the spatial scan statistic with Bernoulli models and a space-time permutation model. Analyses were performed with and without those cases linked to known outbreaks. Laboratory subtyping data (i.e. pulsed field gel electrophoresis (PFGE) and/or phage type) and food and environmental exposure information (e.g. travel, animal contact, poultry and other food item consumption) were used to explore the relatedness of cases within identified clusters.

RESULTS

Spatial, temporal and space-time analyses identified a known outbreak of S. Heidelberg in 2015 (n = 9 cases) and a previously unidentified cluster of S. Heidelberg cases. Most cases (94%) within a cluster detected via a space-time permutation model of S. Heidelberg cases shared an identical PFGE pattern and appeared to represent a true outbreak.

CONCLUSIONS

The spatial scan statistic, and particularly the space-time permutation model, could assist in outbreak identification before laboratory data are available, allowing for faster cluster identification and implementation of control measures.

A case-case study comparing the individual risk factors and symptomatology of Salmonella Heidelberg and Salmonella Typhimurium in Ontario in 2015, following implementation of the Ontario Investigation Tools

Salmonella Heidelberg and Salmonella Typhimurium are among the most common serotypes responsible for human salmonellosis in Ontario. Introduction of the Ontario Investigation Tools (OIT) in 2014 allowed for standardized case investigation and reporting. This study compared the risk factors and symptomatology for sporadic S. Heidelberg and S. Typhimurium cases reported in Ontario in 2015, following implementation of the OIT. Multilevel logistic regression models were applied to assess associations between serotype and individual-level demographic characteristics, exposures and symptoms for sporadic confirmed cases of S. Heidelberg and S. Typhimurium in Ontario in 2015. There were 476 sporadic cases of S. Typhimurium (n = 278) and S. Heidelberg (n = 198) reported in Ontario in 2015. There were significant associations between the odds of the isolate from a case being one of these serotypes, and travel, consumption of sprouts (any type), contact with reptiles and development of malaise, fever or bloody diarrhoea. The S. Typhimurium and S. Heidelberg cases differed in both symptom presentation and risk factors for illness. Case-case comparisons of Salmonella serotypes have some advantages over case-control studies in that these are less susceptible to selection and recall bias while allowing for rapid comparison of cases to identify potential high-risk exposures that are unique to one of the serotypes when compared to the other. Comparing cases of two different Salmonella serotypes can help to highlight risk factors that may be uniquely associated with one serotype, or more strongly associated with one serotype compared to another. This information may be useful for understanding relative source attribution between common serotypes of Salmonella.

Evaluating the risks associated with Shiga-toxin-producing Escherichia coli (STEC) in private well waters in Canada

Shiga-toxin-producing Escherichia coli (STEC) represent a major concern for waterborne disease outbreaks associated with consumption of contaminated groundwater. Over 4 million people rely on private groundwater systems as their primary drinking water source in Canada; many of these systems do not meet current standards for water quality. This manuscript provides a scoping overview of studies examining STEC prevalence and occurrence in groundwater, and it includes a synopsis of the environmental variables affecting survival, transport, persistence, and overall occurrence of these important pathogenic microbes in private groundwater wells used for drinking purposes.

Total coliform and Escherichia coli contamination in rural well water: analysis for passive surveillance

With increasing stress on our water resources and recent waterborne disease outbreaks, understanding the epidemiology of waterborne pathogens is crucial to build surveillance systems. The purpose of this study was to explore techniques for describing microbial water quality in rural drinking water wells, based on spatiotemporal analysis, time series analysis and relative risk mapping. Tests results for Escherichia coli and coliforms from private and small public well water samples, collected between 2004 and 2012 in Alberta, Canada, were used for the analysis. Overall, 14.6 and 1.5% of the wells were total coliform and E. coli-positive, respectively. Private well samples were more often total coliform or E. coli-positive compared with untreated public well samples. Using relative risk mapping we were able to identify areas of higher risk for bacterial contamination of groundwater in the province not previously identified. Incorporation of time series analysis demonstrated peak contamination occurring for E. coli in July and a later peak for total coliforms in September, suggesting a temporal dissociation between these indicators in terms of groundwater quality, and highlighting the potential need to increase monitoring during certain periods of the year.

A descriptive analysis of the spatio-temporal distribution of enteric diseases in New Brunswick, Canada

Background

Enteric diseases affect thousands of Canadians annually and several large outbreaks have occurred due to infection with enteric pathogens. The objectives of this study were to describe the spatial and temporal distributions of reportable Campylobacter, Escherichia coli, Giardia, Salmonella and Shigella from 1994 to 2002 in New Brunswick, Canada. By examining the spatial and temporal distributions of disease incidence, hypotheses as to potential disease risk factors were formulated.

Methods

Time series plots of monthly disease incidence were examined for seasonal and secular trends. Seasonality of disease incidence was evaluated using the temporal scan statistic and seasonal–trend loess (STL) decomposition methods. Secular trends were evaluated using negative binomial regression modeling. The spatial distribution of disease incidence was examined using maps of empirical Bayes smoothed estimates of disease incidence. Spatial clustering was examined by multiple methods, which included Moran’s I and the spatial scan statistic.

Results

The peak incidence of Giardia infections occurred in the spring months. Salmonella incidence exhibited two peaks, one small peak in the spring and a main peak in the summer. Campylobacter and Escherichia coli O157 disease incidence peaked in the summer months. Moran’s I indicated that there was significant positive spatial autocorrelation for the incidence of Campylobacter, Giardia and Salmonella. The spatial scan statistic identified clusters of high disease incidence in the northern areas of the province for Campylobacter, Giardia and Salmonella infections. The incidence of Escherichia coli infections clustered in the south-east and north-east areas of the province, based on the spatial scan statistic results. Shigella infections had the lowest incidence rate and no discernable spatial or temporal patterns were observed.

Conclusions

By using several different spatial and temporal methods a robust picture of the spatial and temporal distributions of enteric disease in New Brunswick was produced. Disease incidence for several reportable enteric pathogens displayed significant geographic clustering indicating that a spatially distributed risk factor may be contributing to disease incidence. Temporal analysis indicated peaks in disease incidence, including previously un-reported peaks.

Area-level global and local clustering of human Salmonella Enteritidis infection rates in the city of Toronto, Canada, 2007-2009

BACKGROUND

Salmonella enterica serotype Enteritidis (S. Enteritidis) remains a major foodborne pathogen in North America yet studies examining the spatial epidemiology of salmonellosis in urban environments are lacking. Our ecological study combined a number of spatial statistical methods with a geographic information system to assess area-level heterogeneity of S. Enteritidis infection rates in the city of Toronto.

METHODS

Data on S. Enteritidis infections between January 1, 2007 and December 31, 2009 were obtained from Ontario's surveillance system, and were grouped and analyzed at the forward sortation area (FSA)-level (an area signified by the first three characters of the postal code). Incidence rates were directly standardized using the FSA-level age- and sex-based standard population. A spatial empirical Bayes method was used to smooth the standardized incidence rates (SIRs). Global clustering of FSAs with high or low non-smoothed SIRs was evaluated using the Getis-Ord G method. Local clustering of FSAs with high, low, or dissimilar non-smoothed SIRs was assessed using the Getis-Ord Gi* and the Local Moran's I methods.

RESULTS

Spatial heterogeneity of S. Enteritidis infection rates was detected across the city of Toronto. The non-smoothed FSA-level SIRs ranged from 0 to 16.9 infections per 100,000 person-years (mean = 6.6), whereas the smoothed SIRs ranged from 2.9 to 11.1 (mean = 6.3). The global Getis-Ord G method showed significant (p ≤ 0.05) maximum spatial clustering of FSAs with high SIRs at 3.3 km. The local Getis-Ord Gi* method identified eight FSAs with significantly high SIRs and one FSA with a significantly low SIR. The Local Moran's I method detected five FSAs with significantly high-high SIRs, one FSA with a significantly low-low SIR, and four significant outlier FSAs (one high-low, and three low-high).

CONCLUSIONS

Salmonella Enteritidis infection rates clustered globally at a small distance band, suggesting clustering of high SIRs in small distinct areas. This finding was supported by the local cluster analyses, where distinct FSAs with high SIRs, mainly in downtown Toronto, were detected. These areas should be evaluated by future studies to identify risk factors of disease in order to implement targeted prevention and control programs. We demonstrated the usefulness of combining several spatial statistical techniques with a geographic information system to detect geographical areas of interest for further study, and to evaluate spatial processes that influenced S. Enteritidis infection rates. Our study methodology could be applied to other foodborne disease surveillance data.

Typhoid fever acquired in the United States, 1999-2010: epidemiology, microbiology, and use of a space-time scan statistic for outbreak detection

Although rare, typhoid fever cases acquired in the United States continue to be reported. Detection and investigation of outbreaks in these domestically acquired cases offer opportunities to identify chronic carriers. We searched surveillance and laboratory databases for domestically acquired typhoid fever cases, used a space-time scan statistic to identify clusters, and classified clusters as outbreaks or non-outbreaks. From 1999 to 2010, domestically acquired cases accounted for 18% of 3373 reported typhoid fever cases; their isolates were less often multidrug-resistant (2% vs. 15%) compared to isolates from travel-associated cases. We identified 28 outbreaks and two possible outbreaks within 45 space-time clusters of ⩾2 domestically acquired cases, including three outbreaks involving ⩾2 molecular subtypes. The approach detected seven of the ten outbreaks published in the literature or reported to CDC. Although this approach did not definitively identify any previously unrecognized outbreaks, it showed the potential to detect outbreaks of typhoid fever that may escape detection by routine analysis of surveillance data. Sixteen outbreaks had been linked to a carrier. Every case of typhoid fever acquired in a non-endemic country warrants thorough investigation. Space-time scan statistics, together with shoe-leather epidemiology and molecular subtyping, may improve outbreak detection.

Spatial-temporal clustering of companion animal enteric syndrome: detection and investigation through the use of electronic medical records from participating private practices

There is interest in the potential of companion animal surveillance to provide data to improve pet health and to provide early warning of environmental hazards to people. We implemented a companion animal surveillance system in Calgary, Alberta and the surrounding communities. Informatics technologies automatically extracted electronic medical records from participating veterinary practices and identified cases of enteric syndrome in the warehoused records. The data were analysed using time-series analyses and a retrospective space–time permutation scan statistic. We identified a seasonal pattern of reports of occurrences of enteric syndromes in companion animals and four statistically significant clusters of enteric syndrome cases. The cases within each cluster were examined and information about the animals involved (species, age, sex), their vaccination history, possible exposure or risk behaviour history, information about disease severity, and the aetiological diagnosis was collected. We then assessed whether the cases within the cluster were unusual and if they represented an animal or public health threat. There was often insufficient information recorded in the medical record to characterize the clusters by aetiology or exposures. Space–time analysis of companion animal enteric syndrome cases found evidence of clustering. Collection of more epidemiologically relevant data would enhance the utility of practice-based companion animal surveillance.

Influenza surveillance and incidence in a rural area in China during the 2009/2010 influenza pandemic

BACKGROUND

Most influenza surveillance is based on data from urban sentinel hospitals; little is known about influenza activity in rural communities. We conducted influenza surveillance in a rural region of China with the aim of detecting influenza activity in the 2009/2010 influenza season.

METHODS

The study was conducted from October 2009 to March 2010. Real-time polymerase chain reaction was used to confirm influenza cases. Over-the-counter (OTC) drug sales were daily collected in drugstores and hospitals/clinics. Space-time scan statistics were used to identify clusters of ILI in community. The incidence rate of ILI/influenza was estimated on the basis of the number of ILI/influenza cases detected by the hospitals/clinics.

RESULTS

A total of 434 ILI cases (3.88% of all consultations) were reported; 64.71% of these cases were influenza A (H1N1) pdm09. The estimated incidence rate of ILI and influenza were 5.19/100 and 0.40/100, respectively. The numbers of ILI cases and OTC drug purchases in the previous 7 days were strongly correlated (Spearman rank correlation coefficient [r] = 0.620, P = 0.001). Four ILI outbreaks were detected by space-time permutation analysis.

CONCLUSIONS

This rural community surveillance detected influenza A (H1N1) pdm09 activity and outbreaks in the 2009/2010 influenza season and enabled estimation of the incidence rate of influenza. It also provides a scientific data for public health measures.

Spatial epidemiology of Escherichia coli O157:H7 in dairy cattle in relation to night roosts Of Sturnus vulgaris (European Starling) in Ohio, USA (2007-2009)

The goal of our study was to use spatial scan statics to determine whether the night roosts of European starlings (Sturnus vulgaris) act as point sources for the dissemination of Escherichia coli O157:H7 among dairy farms. From 2007 to 2009, we collected bovine faecal samples (n = 9000) and starling gastrointestinal contents (n = 430) from 150 dairy farms in northeastern Ohio, USA. Isolates of E. coli O157:H7 recovered from these samples were subtyped using multilocus variable-number tandem repeat analysis (MLVA). Generated MLVA types were used to construct a dendrogram based on a categorical multistate coefficient and unweighted pair-group method with arithmetic mean (UPGMA). Using a focused spatial scan statistic, we identified statistically significant spatial clusters among dairy farms surrounding starling night roosts, with an increased prevalence of E. coli O157:H7-positive bovine faecal pats, increased diversity of distinguishable MLVA types and a greater number of isolates with MLVA types from bovine-starling clades versus bovine-only clades. Thus, our findings are compatible with the hypothesis that starlings have a role in the dissemination of E. coli O157:H7 among dairy farms, and further research into starling management is warranted.

The use of the temporal scan statistic to detect methicillin-resistant Staphylococcus aureus clusters in a community hospital

Background

In healthcare facilities, conventional surveillance techniques using rule-based guidelines may result in under- or over-reporting of methicillin-resistant Staphylococcus aureus (MRSA) outbreaks, as these guidelines are generally unvalidated. The objectives of this study were to investigate the utility of the temporal scan statistic for detecting MRSA clusters, validate clusters using molecular techniques and hospital records, and determine significant differences in the rate of MRSA cases using regression models.

Methods

Patients admitted to a community hospital between August 2006 and February 2011, and identified with MRSA > 48 hours following hospital admission, were included in this study. Between March 2010 and February 2011, MRSA specimens were obtained for spa typing. MRSA clusters were investigated using a retrospective temporal scan statistic. Tests were conducted on a monthly scale and significant clusters were compared to MRSA outbreaks identified by hospital personnel. Associations between the rate of MRSA cases and the variables year, month, and season were investigated using a negative binomial regression model.

Results

During the study period, 735 MRSA cases were identified and 167 MRSA isolates were spa typed. Nine different spa types were identified with spa type 2/t002 (88.6%) the most prevalent. The temporal scan statistic identified significant MRSA clusters at the hospital (n = 2), service (n = 16), and ward (n = 10) levels (P ≤ 0.05). Seven clusters were concordant with nine MRSA outbreaks identified by hospital staff. For the remaining clusters, seven events may have been equivalent to true outbreaks and six clusters demonstrated possible transmission events. The regression analysis indicated years 2009–2011, compared to 2006, and months March and April, compared to January, were associated with an increase in the rate of MRSA cases (P ≤ 0.05).

Conclusions

The application of the temporal scan statistic identified several MRSA clusters that were not detected by hospital personnel. The identification of specific years and months with increased MRSA rates may be attributable to several hospital level factors including the presence of other pathogens. Within hospitals, the incorporation of the temporal scan statistic to standard surveillance techniques is a valuable tool for healthcare workers to evaluate surveillance strategies and aid in the identification of MRSA clusters.

Detection of Clostridium difficile infection clusters, using the temporal scan statistic, in a community hospital in southern Ontario, Canada, 2006-2011

BACKGROUND

In hospitals, Clostridium difficile infection (CDI) surveillance relies on unvalidated guidelines or threshold criteria to identify outbreaks. This can result in false-positive and -negative cluster alarms. The application of statistical methods to identify and understand CDI clusters may be a useful alternative or complement to standard surveillance techniques. The objectives of this study were to investigate the utility of the temporal scan statistic for detecting CDI clusters and determine if there are significant differences in the rate of CDI cases by month, season, and year in a community hospital.

METHODS

Bacteriology reports of patients identified with a CDI from August 2006 to February 2011 were collected. For patients detected with CDI from March 2010 to February 2011, stool specimens were obtained. Clostridium difficile isolates were characterized by ribotyping and investigated for the presence of toxin genes by PCR. CDI clusters were investigated using a retrospective temporal scan test statistic. Statistically significant clusters were compared to known CDI outbreaks within the hospital. A negative binomial regression model was used to identify associations between year, season, month and the rate of CDI cases.

RESULTS

Overall, 86 CDI cases were identified. Eighteen specimens were analyzed and nine ribotypes were classified with ribotype 027 (n = 6) the most prevalent. The temporal scan statistic identified significant CDI clusters at the hospital (n = 5), service (n = 6), and ward (n = 4) levels (P ≤ 0.05). Three clusters were concordant with the one C. difficile outbreak identified by hospital personnel. Two clusters were identified as potential outbreaks. The negative binomial model indicated years 2007-2010 (P ≤ 0.05) had decreased CDI rates compared to 2006 and spring had an increased CDI rate compared to the fall (P = 0.023).

CONCLUSIONS

Application of the temporal scan statistic identified several clusters, including potential outbreaks not detected by hospital personnel. The identification of time periods with decreased or increased CDI rates may have been a result of specific hospital events. Understanding the clustering of CDIs can aid in the interpretation of surveillance data and lead to the development of better early detection systems.

Weather and livestock risk factors for Escherichia coli O157 human infection in Alberta, Canada

This study investigated the extent to which proximity to cattle and weather events in Alberta predispose human populations to E. coli O157 disease. Cases of human E. coli O157 infection in Alberta between 2004 and 2011 were obtained from the province's Communicable Disease Reporting System and Discharge Abstract Database. Regression models based on spatial area incorporated human infection data with livestock and weather covariates. A variety of regression models were applied (i.e. least squares, spatial lag/error, Poisson, negative binomial) to test the most appropriate approach. Ratios for the total number of calves, bulls and beef cows to human population were highlighted as significant cattle density variables in all final best-fitting models. Weather variables were not significant in final regression models averaged over the full study period. Our results provide evidence of a significant association between measures of cattle density and human E. coli O157 disease in Alberta.

Comparison of covariate adjustment methods using space-time scan statistics for food animal syndromic surveillance

BACKGROUND

Abattoir condemnation data show promise as a rich source of data for syndromic surveillance of both animal and zoonotic diseases. However, inherent characteristics of abattoir condemnation data can bias results from space-time cluster detection methods for disease surveillance, and may need to be accounted for using various adjustment methods. The objective of this study was to compare the space-time scan statistics with different abilities to control for covariates and to assess their suitability for food animal syndromic surveillance. Four space-time scan statistic models were used including: animal class adjusted Poisson, space-time permutation, multi-level model adjusted Poisson, and a weighted normal scan statistic using model residuals. The scan statistics were applied to monthly bovine pneumonic lung and "parasitic liver" condemnation data from Ontario provincial abattoirs from 2001-2007.

RESULTS

The number and space-time characteristics of identified clusters often varied between space-time scan tests for both "parasitic liver" and pneumonic lung condemnation data. While there were some similarities between isolated clusters in space, time and/or space-time, overall the results from space-time scan statistics differed substantially depending on the covariate adjustment approach used.

CONCLUSIONS

Variability in results among methods suggests that caution should be used in selecting space-time scan methods for abattoir surveillance. Furthermore, validation of different approaches with simulated or real outbreaks is required before conclusive decisions can be made concerning the best approach for conducting surveillance with these data.

A spatial and temporal analysis of notifiable gastrointestinal illness in the Northwest Territories, Canada, 1991-2008

Background

This is the first study to describe the geographical and temporal distribution of notifiable gastrointestinal illness (NGI) in the Northwest Territories (NWT), Canada. Understanding the distribution of NGI in space and time is important for identifying communities at high risk. Using data derived from the Northwest Territories Communicable Disease Registry (NWT CDR), a number of spatial and temporal techniques were used to explore and analyze NGI incidence from the years 1991 to 2008. Relative risk mapping was used to investigate the variation of disease risk. Scan test statistics were applied to conduct cluster identification in space, time and space-time. Seasonal decomposition of the time series was used to assess seasonal variation and trends in the data.

Results

There was geographic variability in the rates of NGI with higher notifications in the south compared to the north. Incidence of NGI exhibited seasonality with peaks in the fall months for most years. Two possible outbreaks were detected in the fall of 1995 and 2001, of which one coincided with a previously recognized outbreak. Overall, incidence of NGI fluctuated from 1991 to 2001 followed by a tendency for rates to decrease from 2002 to 2008.

Conclusions

The distribution of NGI notifications varied widely according to geographic region, season and year. While the analyses highlighted a possible bias in the surveillance data, this information is beneficial for generating hypotheses about risk factors for infection.

Improving resolution of public health surveillance for human Salmonella enterica serovar Typhimurium infection: 3 years of prospective multiple-locus variable-number tandem-repeat analysis (MLVA)

Background

Prospective typing of Salmonella enterica serovar Typhimurium (STM) by multiple-locus variable-number tandem-repeat analysis (MLVA) can assist in identifying clusters of STM cases that might otherwise have gone unrecognised, as well as sources of sporadic and outbreak cases. This paper describes the dynamics of human STM infection in a prospective study of STM MLVA typing for public health surveillance.

Methods

During a three-year period between August 2007 and September 2010 all confirmed STM isolates were fingerprinted using MLVA as part of the New South Wales (NSW) state public health surveillance program.

Results

A total of 4,920 STM isolates were typed and a subset of 4,377 human isolates was included in the analysis. The STM spectrum was dominated by a small number of phage types, including DT170 (44.6% of all isolates), DT135 (13.9%), DT9 (10.8%), DT44 (4.5%) and DT126 (4.5%). There was a difference in the discriminatory power of MLVA types within endemic phage types: Simpson's index of diversity ranged from 0.109 and 0.113 for DTs 9 and 135 to 0.172 and 0.269 for DTs 170 and 44, respectively. 66 distinct STM clusters were observed ranging in size from 5 to 180 cases and in duration from 4 weeks to 25 weeks. 43 clusters had novel MLVA types and 23 represented recurrences of previously recorded MLVA types. The diversity of the STM population remained relatively constant over time. The gradual increase in the number of STM cases during the study was not related to significant changes in the number of clusters or their size. 667 different MLVA types or patterns were observed.

Conclusions

Prospective MLVA typing of STM allows the detection of community outbreaks and demonstrates the sustained level of STM diversity that accompanies the increasing incidence of human STM infections. The monitoring of novel and persistent MLVA types offers a new benchmark for STM surveillance.

A part of this study was presented at the MEEGID × (Molecular Epidemiology and Evolutionary Genetics of Infectious Diseases) Conference, 3-5 November 2010, Amsterdam, The Netherlands

Gastrointestinal disease outbreak detection using multiple data streams from electronic medical records

BACKGROUND

Passive reporting and laboratory testing delays may limit gastrointestinal (GI) disease outbreak detection. Healthcare systems routinely collect clinical data in electronic medical records (EMRs) that could be used for surveillance. This study's primary objective was to identify data streams from EMRs that may perform well for GI outbreak detection.

METHODS

Zip code-specific daily episode counts in 2009 were generated for 22 syndromic and laboratory-based data streams from Kaiser Permanente Northern California EMRs, covering 3.3 million members. Data streams included outpatient and inpatient diagnosis codes, antidiarrheal medication dispensings, stool culture orders, and positive microbiology tests for six GI pathogens. Prospective daily surveillance was mimicked using the space-time permutation scan statistic in single and multi-stream analyses, and space-time clusters were identified. Serotype relatedness was assessed for isolates in two Salmonella clusters.

RESULTS

Potential outbreaks included a cluster of 18 stool cultures ordered over 5 days in one zip code and a Salmonella cluster in three zip codes over 9 days, in which at least five of six cases had the same rare serotype. In all, 28 potential outbreaks were identified using single stream analyses, with signals in outpatient diagnosis codes most common. Multi-stream analyses identified additional potential outbreaks and in one example, improved the timeliness of detection.

CONCLUSIONS

GI disease-related data streams can be used to identify potential outbreaks when generated from EMRs with extensive regional coverage. This process can supplement traditional GI outbreak reports to health departments, which frequently consist of outbreaks in well-defined settings (e.g., day care centers and restaurants) with no laboratory-confirmed pathogen. Data streams most promising for surveillance included microbiology test results, stool culture orders, and outpatient diagnoses. In particular, clusters of microbiology tests positive for specific pathogens could be identified in EMRs and used to prioritize further testing at state health departments, potentially improving outbreak detection.

Agricultural, socioeconomic and environmental variables as risks for human verotoxigenic Escherichia coli (VTEC) infection in Finland

Background

Verotoxigenic E. coli (VTEC) is the cause of severe gastrointestinal infection especially among infants. Between 10 and 20 cases are reported annually to the National Infectious Disease Register (NIDR) in Finland. The aim of this study was to identify explanatory variables for VTEC infections reported to the NIDR in Finland between 1997 and 2006. We applied a hurdle model, applicable for a dataset with an excess of zeros.

Methods

We enrolled 131 domestically acquired primary cases of VTEC between 1997 and 2006 from routine surveillance data. The isolated strains were characterized by virulence type, serogroup, phage type and pulsed-field gel electrophoresis. By applying a two-part Bayesian hurdle model to infectious disease surveillance data, we were able to create a model in which the covariates were associated with the probability for occurrence of the cases in the logistic regression part and the magnitude of covariate changes in the Poisson regression part if cases do occur. The model also included spatial correlations between neighbouring municipalities.

Results

The average annual incidence rate was 4.8 cases per million inhabitants based on the cases as reported to the NIDR. Of the 131 cases, 74 VTEC O157 and 58 non-O157 strains were isolated (one person had dual infections). The number of bulls per human population and the proportion of the population with a higher education were associated with an increased occurrence and incidence of human VTEC infections in 70 (17%) of 416 of Finnish municipalities. In addition, the proportion of fresh water per area, the proportion of cultivated land per area and the proportion of low income households with children were associated with increased incidence of VTEC infections.

Conclusions

With hurdle models we were able to distinguish between risk factors for the occurrence of the disease and the incidence of the disease for data characterised by an excess of zeros. The density of bulls and the proportion of the population with higher education were significant both for occurrence and incidence, while the proportion of fresh water, cultivated land, and the proportion of low income households with children were significant for the incidence of the disease.

Seasonal variation of HUS occurrence and VTEC infection in children with acute diarrhoea from Argentina

In order to study the seasonality of haemolytic uraemic syndrome (HUS) and verotoxigenic Escherichia coli (VTEC) infection in children, 437 patients under 6 years of age with acute diarrhoea were studied, 8% of whom progressed to HUS. VTEC was found in 10% of all of the stool samples analysed and seasonal occurrence of HUS (p < 0.01) was confirmed. VTEC infection was more prevalent in warm months, although the differences were not statistically significant. Moreover, a significant difference in the detection of O157:H7 serotype and in the vt profile between cold and warm months (autumn and winter; spring and summer, respectively) was established. The O157:H7 serotype was isolated more frequently during warm months. Moreover, a predominance of vt (2) was noted, which was partially replaced by the combination of vt (1) with vt (2) in the cold season. The results of this study indicate the seasonal variation of the disease and the presence of serotype O157:H7 and the vt types. They also reinforce the need to develop prevention programmes considering the seasonal pattern of the disease, which would generate an impact on public health. Control strategies of the pathogen in cattle in the most risky season of the year would also be of benefit.

Utilizing spatiotemporal analysis of influenza-like illness and rapid tests to focus swine-origin influenza virus intervention

In the spring of 2009, a novel strain of H1N1 swine-origin influenza A virus (S-OIV) emerged in Mexico and the United States, and soon after was declared a pandemic by the World Health Organization. This work examined the ability of real-time reports of influenza-like illness (ILI) symptoms and rapid influenza diagnostic tests (RIDTs) to approximate the spatiotemporal distribution of PCR-confirmed S-OIV cases for the purposes of focusing local intervention efforts. Cluster and age adjusted relative risk patterns of ILI, RIDT, and S-OIV were assessed at a fine spatial scale at different time and space extents within Cameron County, Texas on the US-Mexico border. Space-time patterns of ILI and RIDT were found to effectively characterize the areas with highest geographical risk of S-OIV within the first two weeks of the outbreak. Based on these results, ILI and/or RIDT may prove to be acceptable indicators of the location of S-OIV hotspots. Given that S-OIV data is often difficult to obtain real-time during an outbreak; these findings may be of use to public health officials targeting prevention and response efforts during future flu outbreaks.

Spatial and temporal clustering of calcium oxalate and magnesium ammonium phosphate uroliths in dogs living in Ontario, Canada between 1998 and 2006

Using the spatial scan statistic with a Bernoulli model, in a comparison of the two most common canine uroliths, calcium oxalate (CaOx) and magnesium ammonium phosphate (struvite) we determined whether there was evidence of spatial and/or temporal clustering of each urolith type based on canine submissions from Ontario to the Canadian Veterinary Urolith Centre (CVUC) between 1998 and 2006. During this period, there were 10,478 canine submissions, excluding cases that were identified as recurrent. We were able to georeference approximately 93% of these incident cases. After adjusting for spatial and temporal distributions of dogs based on the demographic risk factors of age, sex, and breed-type, statistically significant spatial and temporal clusters were present for both CaOx and struvite urolith types. A purely temporal struvite cluster occurred between February 10, 1998 and December 20, 2001, whereas, a purely temporal CaOx cluster occurred between September 2, 2005 and December 21, 2006. Hypotheses to explain the spatial clustering of uroliths include variation in the spatial distribution of water hardness, diet-type, access to veterinary care, and the use of surgical versus medical therapies to treat these uroliths. Based on the cluster locations, water hardness was unlikely to explain the spatial difference between the two cluster types whereas variables related to human population density were more consistent with our findings; the CaOx cluster occurred in the highest population density area of Ontario, and the struvite cluster occurred in the lowest population density area of southern Ontario. The temporal struvite cluster at the beginning and CaOx cluster at the end of the study period reflect a similar trend away from struvite towards CaOx urolithiasis among both canines and humans in the developed countries of North America and Europe.

A new application of spatiotemporal analysis for detecting demographic variations in AIDS mortality: an example from Florida

The purpose of the present study was to characterize, geographically and temporally, the patterns of acquired immune deficiency syndrome (AIDS) death disparity in 67 Florida jurisdictions, and to determine if the detected trends varied according to age, race, and sex. The space-time scan statistic proposed by Kulldorff et al was used to examine the excess AIDS deaths that occurred between 1987 and 2004. Results were geographically referenced in maps using EpiInfo and EpiMap made available by the Centers for Disease Control. Miami-Dade and the nearby counties including Broward, Martin, and Palm Beach are the most likely clusters (observed/expected: 1505.16) with temporal dimension (also called cluster's age) persisting from 1996 to the present. Union county had the longest cluster for the cluster period 1987-1998, but not for 1999-2004. African-Americans contributed to more clusters compared with whites. Time trends indicated that AIDS mortality peaked in 1995 and then sharply dropped until 1998, when the decrease stopped. By accounting for the temporal dimension of disease clustering, the present study revealed the persistence of geographic clusters, which is not often provided by other geographic detection methods. These findings may be informative for medical resource allocation and better focus public health intervention strategies for AIDS care.

Differences between notifiable and administrative health information in the spatial-temporal surveillance of enteric infections

PURPOSE

The purpose of this study is to compare the spatial and temporal information generated from two distinct health data sources available for the surveillance of intestinal infections associated with Escherichia coli O157:H7.

METHODS

Our study area is the province of Alberta, Canada. Data are from two sources: a fee-for-service administrative health data system and a notifiable disease data reporting system. The study period is between 1999 and 2005. We compare the systems by observing correlations in the infections over time, the variability in the overall distribution of cases (as measured by a geographic dissimilarity index), and the relative locations of spatial-temporal clusters of infection.

RESULTS

Our results indicate considerable variability in information generated from these two systems. The geographic distribution of cases varies considerably, with annual indices of dissimilarity suggesting considerable variation in the geographic distribution of cases throughout the study period (D=0.445). The temporal patterns identified by these two sources of information are negatively correlated (-0.40, p<0.001). Notifiable disease clusters occur in the summer in southern regions of the province, whereas cases identified from administrative health data system cluster in the winter season, and further to the north.

CONCLUSIONS

Notifiable disease data may suffer from selection bias; administrative health data may be insufficiently precise without laboratory confirmation. Our results illustrate differences in the spatial and temporal information generated from these two systems of case identification. Future surveillance of gastrointestinal illness of infectious origin may benefit from case ascertainment algorithms based on both sources of data.

A multi-level approach for investigating socio-economic and agricultural risk factors associated with rates of reported cases of Escherichia coli O157 in humans in Alberta, Canada

Using negative binomial and multi-level Poisson models, the authors determined the statistical significance of agricultural and socio-economic risk factors for rates of reported disease associated with Escherichia coli O157 in census subdivisions (CSDs) in Alberta, Canada, 2000-2002. Variables relating to population stability, aboriginal composition of the CSDs, and the economic relationship between CSDs and urban centres were significant risk factors. The percentage of individuals living in low-income households was not a statistically significant risk factor for rates of disease. The statistical significance of cattle density, recorded at a higher geographical level, depended on the method used to correct for overdispersion, the number of levels included in the multi-level models, and the choice of using all reported cases or only sporadic cases. Our results highlight the importance of local socio-economic risk factors in determining rates of disease associated with E. coli O157, but their relationship with individual risk factors requires further evaluation.

Tracking the spatial diffusion of influenza and norovirus using telehealth data: a spatiotemporal analysis of syndromic data

BACKGROUND

Telehealth systems have a large potential for informing public health authorities in an early stage of outbreaks of communicable disease. Influenza and norovirus are common viruses that cause significant respiratory and gastrointestinal disease worldwide. Data about these viruses are not routinely mapped for surveillance purposes in the UK, so the spatial diffusion of national outbreaks and epidemics is not known as such incidents occur. We aim to describe the geographical origin and diffusion of rises in fever and vomiting calls to a national telehealth system, and consider the usefulness of these findings for influenza and norovirus surveillance.

METHODS

Data about fever calls (5- to 14-year-old age group) and vomiting calls (> or = 5-year-old age group) in school-age children, proxies for influenza and norovirus, respectively, were extracted from the NHS Direct national telehealth database for the period June 2005 to May 2006. The SaTScan space-time permutation model was used to retrospectively detect statistically significant clusters of calls on a week-by-week basis. These syndromic results were validated against existing laboratory and clinical surveillance data.

RESULTS

We identified two distinct periods of elevated fever calls. The first originated in the North-West of England during November 2005 and spread in a south-east direction, the second began in Central England during January 2006 and moved southwards. The timing, geographical location, and age structure of these rises in fever calls were similar to a national influenza B outbreak that occurred during winter 2005-2006. We also identified significantly elevated levels of vomiting calls in South-East England during winter 2005-2006.

CONCLUSION

Spatiotemporal analyses of telehealth data, specifically fever calls, provided a timely and unique description of the evolution of a national influenza outbreak. In a similar way the tool may be useful for tracking norovirus, although the lack of consistent comparison data makes this more difficult to assess. In interpreting these results, care must be taken to consider other infectious and non-infectious causes of fever and vomiting. The scan statistic should be considered for spatial analyses of telehealth data elsewhere and will be used to initiate prospective geographical surveillance of influenza in England.

Epidemiological characteristics of reported sporadic and outbreak cases of E. coli O157 in people from Alberta, Canada (2000-2002): methodological challenges of comparing clustered to unclustered data

Using multivariable models, we compared whether there were significant differences between reported outbreak and sporadic cases in terms of their sex, age, and mode and site of disease transmission. We also determined the potential role of administrative, temporal, and spatial factors within these models. We compared a variety of approaches to account for clustering of cases in outbreaks including weighted logistic regression, random effects models, general estimating equations, robust variance estimates, and the random selection of one case from each outbreak. Age and mode of transmission were the only epidemiologically and statistically significant covariates in our final models using the above approaches. Weighing observations in a logistic regression model by the inverse of their outbreak size appeared to be a relatively robust and valid means for modelling these data. Some analytical techniques, designed to account for clustering, had difficulty converging or producing realistic measures of association.

The use of randomization tests to assess the degree of similarity in PFGE patterns of E. coli O157 isolates from known outbreaks and statistical space-time clusters

Using isolates from reported cases of Escherichia coli O157 from Alberta, Canada in 2002, we applied randomization tests to determine if cases associated with an outbreak or statistical space-time cluster had more similar pulsed-field gel electrophoresis patterns, based on Dice coefficients, than expected by chance alone. Within each outbreak and space-time cluster, we assessed the mean, median, 25th percentile, 75th percentile, standard deviation, coefficient of variation, and interquartile range of the Dice coefficients of each pairwise comparison among the isolates. To assess the statistical significance of measures of location (e.g. mean) and variation (e.g. standard deviation) we created randomization distributions using all isolates or only isolates from sporadic cases. We determined that randomization tests are an appropriate tool for evaluating the similarity among isolates from cases that have been linked epidemiologically or statistically. We found little difference between using all cases or only sporadic cases when creating our randomization distributions.