- Research article
- Open Open Peer Review
Agricultural, socioeconomic and environmental variables as risks for human verotoxigenic Escherichia coli(VTEC) infection in Finland
BMC Infectious Diseasesvolume 11, Article number: 275 (2011)
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.
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.
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.
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.
Verotoxigenic E. coli (VTEC) is a highly infectious causative agent for severe gastrointestinal infection with haemorrhagic diarrhoea, which can lead to haemolytic uraemic syndrome (HUS) or thrombocytopaenic purpura (TTP) . Most of the diagnosed VTEC infections in developed countries, including Finland, are of the serogroup O157, but an increasing number of cases are found with non-O157 serogroups O26, O111, O91, O103 and O145 [2, 3]. The majority of VTEC infections occur among young children, predominantly during the summer months [2–4]. VTEC infection is rare in Finland, with only 10-20 cases reported annually to the National Infectious Disease Register (NIDR) .
Several risk factors for VTEC infections have been identified from outbreak data, modelling studies, and descriptive and analytical epidemiological studies. Cattle are considered the main reservoir for VTEC [6–8]. Living, visiting or working in close proximity to a farm with cattle or having a household contact working in such a farm have been found to be important risk factors [9–11]. Cattle density or VTEC-positive cattle density have been established as important explanatory variables for VTEC infections [4, 12–16] or cases of HUS . Also farm density per area of land  and animal manure application to soil are important . Undercooked beef meat or meat products, water, vegetables and unpasteurized milk, ice cream, juice and cider can be sources of VTEC infections , and consumption of burgers and cold meat cuts are established risk factors [9, 18]. Pathogenic strains of VTEC have been detected in minced beef meat, beef meat cubes and cuts, raw cattle, sheep and goat milk and sprouts . Lake water or contaminated municipal water systems have caused outbreaks of VTEC O157 [17, 19, 20] and consumption of untreated surface water is a risk factor for non-O157 .
Poisson distribution based regression models are standard models for count data, including for infectious disease surveillance data . However, if the data has an excess of zeros, the assumptions of the Poisson distributed models are no longer valid. Models based on negative binomial distribution may slightly correct for this problem, but when there is a severe inflation of zeros, more specific models are needed. Zero-inflated and hurdle models have been developed to overcome this problem [21, 22]. With the use of a complementary log-log link function, the hurdle model reduces to a standard Poisson regression if the data is Poisson distributed . These types of models have been applied recently for modelling sleeping sickness , cholera prevalence  and bacterial counts .
The aim of this study was to identify explanatory variables for VTEC infections reported to the NIDR in Finland between 1997 and 2006. VTEC infection is a rare disease in Finland and we used Poisson distribution based hurdle models with variable selection based on in-depth questionnaires and known risk factors. We chose the hurdle model based on its ability to handle an over-dispersion of zeros and applicability to scarce data. The incidence rate of VTEC infection at municipality level was used as an outcome variable and explanatory variables included various agricultural, socioeconomic and environmental factors. These explanatory variables also included several factors that have not been studied with VTEC infections previously. The null hypothesis that the populations in different municipalities in Finland are at the same risk of contracting VTEC infection was tested using a Poisson distribution based hurdle model.
The presence, isolation and biochemical identification of VTEC bacteria from stool cultures was done as previously described . The isolates were O grouped as previously described  with the specific O antisera obtained from THL (National Institute for Health and Welfare), Statens Serum Institute (SSI, Denmark) and Oxoid (Oxoid, Basingstoke, UK). The isolates of the O157 strains were phage typed as previously described . Chromosomal DNA of the isolates was characterized by the internationally standardized pulsed-field gel electrophoresis (PFGE) method .
Definitions for the microbiologically confirmed VTEC cases
A case of VTEC infection with no history of foreign travel during the two weeks prior to onset of symptoms.
The VTEC case with "earliest onset of symptoms"- date or isolation date within a household or outbreak.
A VTEC case living in the same household with the index case.
Two or more epidemiologically and microbiologically linked VTEC cases caused by strains of identical subtype and not of the same household.
A VTEC case with no epidemiological or microbiological link to another VTEC case.
Index or sporadic case.
Surveillance of VTEC in Finland
For all cases reported to the NIDR since 1997, date of birth, sex, place of residence, the earliest reported date, sampling date, travel information, and death notification were reported. An in-depth interview as a part of the surveillance has been conducted since 1998 , with the questionnaire reformulated in 2003. The interviews were conducted either by THL or the local health district. The information and question types that differed between the questionnaires were either used in conjunction or information was combined from the two forms. However, in some cases data were insufficient for a full analysis of the in-depth questionnaires. The study material is part of routine enhanced surveillance carried out in Finland. The data presented does not include any individual data so that any person could be identified. The interviewed case patients have given informed consent for the surveillance. The data was accessed from the NIDR based on the justification on the infectious disease law (583/1986) and regulation (786/1986).
The information that was available from both forms or could be obtained by combining information included: Place of residence (postal code), attending school or day care and having VTEC positives in the same household. Clinical information included presence of diarrhoea, visible blood in the stool, HUS or TTP as a complication, admittance to a hospital, and death. Certain risk factors were asked in one or both forms, including living or visiting farms with domestic animals, consumption of raw or undercooked meat, unpasteurized milk or untreated water, and swimming in natural water.
As the infectious disease surveillance data is count data and there was an excess number of zeros in the data (mean number of cases 0.31 per municipality with variance of 1.2), we developed a Bayesian Poisson complementary-log log (clog log) hurdle model [22, 26] with a spatial correlation factor for taking into account the neighbouring municipalities [27, 28]. If a Poisson or negative binomial distribution based count model is applied to data with an excess of zeros without addressing these mixtures, the model can be strongly affected. The hurdle model was chosen due to its ability to handle the over-dispersion of zeros. In the zero-inflated models the count of zeros is mixed into two components, while in the hurdle models the zero and positive counts are modelled by separate processes. Zero-inflated models are designed to handle excessive zero counts, while hurdle models can also be adjusted for too few zeros. As the data were quite scarce and there was no need to distinguish between structural and sampling zeros, the hurdle model was chosen for the present study [21, 22]. The hurdle model is composed of two parts, a binary model (logistic regression) generating positive counts versus zero counts, and a zero truncated Poisson model generating positive counts only. In a hurdle model, the data is considered as being generated by a process that generates positive counts only after clearing a zero barrier, or hurdle. Until the hurdle is cleared, the process generates a binary response and the covariates associate with the probability for occurrence of the cases in the logistic regression part. After the hurdle is cleared, the process generates only a positive counts response and the covariates estimate the magnitude for the incidence of VTEC cases. We used a Bayesian full-likelihood approach to modelling the data, taking into account also the missingness in the covariate data (using imputation and missingness indicators).
Domestic primary cases of VTEC were included in the analysis as outcome variables for the 10-year follow-up period. These were assigned with home municipality according to information on the in-depth questionnaires or NIDR information using a municipality division of 416 municipalities in Finland in 2007. The explanatory variables for modelling included various agricultural, socioeconomic and environmental factors. The results from the descriptive epidemiology of the present study and literature were used to identify potential variables to be included into the model. We performed a univariable analysis in the hurdle model with all potential explanatory variables and in the multivariable analysis we tested those variables that had p-values < 0.20 in the univariable analysis. Of the correlated variables, only the most significant ones and one per correlation group were included in the final model within each group, as marked in Additional file 1. Correlation coefficients within each group (1-5) were significant and > 0.80. To identify explanatory variables in the multivariable model, Gibbs' variable selection was used  with all explanatory variables as listed in Additional file 1. Those variables with a posterior inclusion probability of > 0.5 in the Gibbs' variable selection were included in the final model. The significant variables were tested with cross validation in the univariable analysis to identify influential values. For details of the model, see Additional file 2 and Additional file 3.
The statistical packages used included R (version 2.11.1), SPSS (version 18), Arcgis (version 9.3) (for data pre-processing), WinBUGS (version 1.4.3) and Stata (version 9.2). The hurdle models are described in the Additional file 2 and Additional file 3.
Description of the cases
Between 1997 and 2006, 247 cases of microbiologically confirmed VTEC infections were reported to the NIDR in Finland. There was an average annual incidence rate of 4.8 (range 1.9-12.0) per million inhabitants and 25 (10-62) cases on average were identified annually. Of all cases, 46 cases (19%) had a history of foreign travel; these cases were excluded from further analyses. Outbreak cases accounted for 29 (12%) and household cases for 77 (31%) of all the cases, and of these, 36 domestic index cases were enrolled. Altogether 131 domestically acquired primary cases were enrolled. Of the 131 cases, 95 (73%) cases were sporadic.
Of all 247 cases, 249 strains were isolated from primary stool cultures of patients or their close contacts with VTEC infection. Two patients had a dual infection caused by two VTEC strains. Of the 131 domestically acquired primary cases enrolled in this study, 132 VTEC strains were isolated: 74 (56%) were serogroup O157 and 58 (44%) were of the non-O157 serogroup; one case had dual infections of both serogroups. Of the O157 strains, 60 (81%) typed to sorbitol negative and 14 (19%) to the sorbitol positive O157 group. Within the O157 group, the most predominant phage types (PT) were PT2, PT88 or its variant, PT8 and PT4. The genomic DNA of the strains gave 33 and 50 different PFGE patterns for the VTEC O157 and non-O157 strains, respectively. The non-O157 strains typed into 20 serogroups, the most common O groups being O103, O145 and O26.
Of the 131 primary cases, 51 (39%) were aged younger than 5 years and 23 (18%) were older than 40 years. Half (66, 50%) were male, with the gender distribution also being equal in the age groups (data not shown). Interviews were conducted for 83 (63%) cases. The postcode for the home address was known for 25 (30%) of the cases. Secondary household VTEC positives were reported in 24 cases (29%).
The symptoms were typical of VTEC infection: 74 (94%) had diarrhoea and 69 (90%) visible blood in the stools. HUS developed in 21 (28%) cases, and of these, 5 (7%) had TTP; in total, 61 (74%) were admitted to a hospital, and none of the cases died. Information on the reported risk factors was scarce: of the 37 cases aged 0-6 years, 8 (22%) attended day care. Living or visiting a farm with domestic animals was reported in 17 (22%) cases. Unpasteurized milk products were consumed in 16 (20%) cases, raw or undercooked meat was consumed in 17 (24%) cases. Consumption of untreated water was reported in 6 (17%) and swimming in natural water in 3 (8%) cases.
Modelling of the data
The modelling unit was the municipality and between 1 and 11 cases were identified in 70/416 (17%) municipalities (Figure 1). Of all the cattle, socioeconomic and environmental variables, 17 variables were included in the multivariable models (Additional file 1).
In the multivariable model, bulls per human population and level of higher education per population were significant in the zero part of the hurdle models, indicating the importance of these variables for the occurrence of the infection (Table 1). The number of bulls per human population, the number of farms under cultivation per area, the proportion of fresh water to total surface area, the level of higher education per population and the proportion of low income households with children were positively significant in the zero-truncated Poisson of the hurdle models, indicating significance relating to the incidence of VTEC cases. The zero part of the hurdle model should be interpreted so that 1% change in bulls per population in a certain municipality or 1% change in level of higher education per population increases the approximated odds for the occurrence of VTEC by exp(0.01*βbulls/1.386)-1 = 4% and exp(0.01*βeduc/1.386)-1 = 4%, respectively. Likewise, the 1% change in bulls, farms, fresh water, educational level and in low income households in the Poisson part increases the standardised incidence ratio (SIR) by exp(0.01*βbulls) = 5%, exp(0.01*βfarms) = 2%, exp(0.01*βfresh_water) = 2%, exp(0.01*βeduc) = 6%, exp(0.01*βlow_ income) = 7%, respectively .
Interestingly, using a cross validation procedure in the univariable or multivariable hurdle models, several more influential values were found. When two influential observations (in both zero and Poisson parts of the model) were removed from the data, the fresh water variable became negatively significant. When examining these municipalities in more detail, it became plausible that there may have been prolonged, undetected outbreaks. One of these municipalities was a medium size city with nine cases, with six cases during one calendar year and of those, four had identical PFGE profiles. The other was a small municipality with a population of just over 10 000 that had four cases during a two-year period with different PFGE patterns. Furthermore, these municipalities were among those with the highest proportion of fresh water, suggestive of water borne outbreaks. Also the number of cattle was high in these municipalities.
VTEC infection has a low average annual incidence rate of 4.8 per million inhabitants in Finland. The cases described in the present study were predominantly children with typical symptoms of VTEC infection: Most cases had haemorrhagic diarrhoea and some cases of HUS were identified. Half of the cases were infected with O157 serogroup strains. Detailed microbiological analysis including genotyping was essential to distinguish between outbreak and non-outbreak strains. Approximately a third of cases had a household contact who was VTEC infection positive, indicating direct transmission of the infection. Many of the cases reported visits to farms or consuming undercooked meat or unpasteurized milk; however, it was difficult to assess the impact of the risk factors obtained from the in-depth questionnaires as the information was scarce. Bias due to under-reporting of the cases is considered to be low for VTEC; however, diagnostic differences may occur between health care districts.
Hurdle modelling has not yet been widely used on infectious disease surveillance data [23, 25, 31]. One reason for using the hurdle model in the present study was the fact that the variance was much higher than the mean in the observed number of cases per municipality. The data was also characterised by an excess of zeros. As the zero and the truncated Poisson parts of the model gave different sets of significant variables, this further suggested the need for this type of model. If the coefficients of the binary part and the truncated Poisson part are close to each other, the hurdle model may be replaced by a standard Poisson regression model . From the technical point of view, the prior distributions used were as weakly informative as possible, which guaranteed the convergence of the simulations. The results from the simulations were robust, especially to changes in the variances of the spatial factors. We also recommend using a full-likelihood approach in the multivariable analysis to explain the mechanism for the missing values.
We identified the proportion of beef cattle to human population and the proportion of population with higher education as being associated with increased occurrence and incidence of VTEC infections. The proportion of the population with higher education is likely to represent an indirect indicator of consumption habits, possibly eating undercooked beef steak. In addition, the number of farms under cultivation per area, fresh water per area and the proportion of low income households with children were associated with increased incidence of the infection. These are likely to represent routes for the spread of the infection; the low income households may be an indirect indicator of consumption habits. The model without the spatial variable essentially gave the same results. Our results concur with a US study on ecological factors for gastrointestinal infections, where the proportion of the population living on a farm (positive association), low education (negative association) and living in the south (negative association) were the most important socioeconomic factors for acquiring a VTEC infection .
The finding of beef cattle as a strong explanatory variable was not unexpected as cattle have been found to be an important risk factor for VTEC infections in epidemiological and modelling studies [9–11, 13, 14]. The VTEC prevalence is low in cattle in the Nordic countries, including Finland [33, 34]. The overall cattle density  or dairy cattle or calves density  have also been found be significant, but these differences may reflect correlations between the cattle variables, variation in agricultural practices, and the ways in which data are collected. The significance of beef cattle over the other cattle variables may be due to a susceptibility of young animals to VTEC, raising the bull cattle in finishing units with multiple sources of animals, and the faecal contamination of the bulls and the environment during the finishing period . The beef cattle density was an important factor both for the occurrence of the disease--confirming the role of cattle as a reservoir of the infection--and the incidence of the infection, indicating that contact with cattle may be an important route for spreading the infection. A relationship between human VTEC infections and VTEC positive cattle density has been found . Also indirect cattle variables like the proportion of land where manure is applied has been found to be significant ; in the present study the proportion of cultivated land was also a significant risk factor for the incidence of the VTEC infections. Spreading animal manure on fields is common practice in Finland, with most of the manure coming from cattle .
The result for the proportion of fresh water surface per area is intriguing. Previous studies suggest drinking and swimming water as a source of VTEC infections and outbreaks [6, 18–20, 37]. However, two influential values had a strong effect for this variable in the present study. It would be interesting to examine the risk factors separately for water borne cases; this information was not widely available. There is recent evidence that contaminated water is an important route of infection in cattle  Campylobacter infections have been found to have a positive association with water pipe-length per population and a negative association with proportion of population receiving household water from a public supply . It is noteworthy that the proportion of salted water was not a significant variable in the multivariable analysis. Finland has thousands of small lakes, but it also has a long sea coast along the western and southern borders. In a Swedish study, it was found that most of the VTEC cases lived along the coastline, near lakes or along rivers; however, the type of water was not specified . VTEC prevalence studies have indicated that cases often have a private well or have experienced a water failure recently [12, 40].
The finding that the proportion of the population with a higher education and the proportion of low income households with children being associated with increased incidence of the infection is somewhat contradictory. However, we feel that these two variables may explain different phenomena. There is conflicting evidence about educational level as a risk factor for VTEC infections [32, 41]. However, it is known that people with higher education eat raw or undercooked ground beef more frequently , a known risk factor for VTEC infection . The proportion of households where the main supporter of the household works in agriculture did not remain in the final model; however, there was a correlation between this and other agricultural variables. There may also be immunity: It has been found that in outbreaks, long-term exposure can lead to immunity . In other studies, having a household member with an agricultural occupation was found to be a risk factor . The most likely risk factor found for sporadic VTEC infection is having contact with a farm .
The study was limited by the low number of VTEC cases; therefore spatiotemporal analysis was not used. This was compensated for by using specific models designed for scarce data and performing a purely spatial analysis. The seasonal trends could not be evaluated properly due to the low number of cases. Another limitation of the study was the unavailability of data regarding consumption of beef or other food items per municipality. It has been found in Scotland that of all VTEC infections, 16% of cases had a significant food-related exposure, such as eating high-risk foods or having a contact working with meat . Considering the VTEC infections in our study and the strong evidence of domestically acquired infections, transmission via imported food stuff would be plausible, based on the similar pheno- and genotypic characteristics of the non-sorbitol fermenting, sorbitol fermenting and non-O157 strains isolated in Finland and other European countries [7, 43, 44].
In conclusion, cattle are the main reservoir for VTEC infections as well as being also an important factor for the incidence of the disease. In addition, higher education per population was found to be associated both with the occurrence and the incidence of the disease, implying that food exposures may play a role. The role of fresh water and other environmentally related variables warrants further studies. Socioeconomic factors like low income and educational level expedite the incidence of VTEC infection. The recommendations for control measures for the prevention of VTEC infections given by the WHO include good hygiene in animal and slaughter process and food retail, consumer education on proper handling of foods, and up-to-date international and national legislation and regulation .
We enrolled 131 cases of VTEC with 74 VTEC O157 and 58 non-O157 strains isolated (one person had dual infections). With Poisson hurdle models we were able to distinguish between risk factors for the occurrence of the disease (zero part of the model) and incidence of the disease (zero truncated Poisson part of the model) for data characterised by an excess of zeros in the case counts. The number of bulls per human population and the proportion of population with a higher education were identified as being associated with increased occurrence and incidence of human VTEC O157 and non-O157 infections. 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. We recommend using a Bayesian framework with imputation for the missing values of the explanatory variables in conjunction with the missingness equations.
Henderson H: Direct and indirect zoonotic transmission of Shiga toxin-producing Escherichia coli. J Am Vet Med Assoc. 2008, 232 (6): 848-859. 10.2460/javma.232.6.848.
Eklund M, Nuorti JP, Ruutu P, Siitonen A: Shigatoxigenic Escherichia coli (STEC) infections in Finland during 1998-2002: a population-based surveillance study. Epidemiol Infect. 2005, 133 (5): 845-852. 10.1017/S0950268805004450.
Brooks JT, Sowers EG, Wells JG, Greene KD, Griffin PM, Hoekstra RM, Strockbine NA: Non-O157 Shiga toxin-producing Escherichia coli infections in the United States, 1983-2002. J Infect Dis. 2005, 192 (8): 1422-1429. 10.1086/466536.
Pearl DL, Louie M, Chui L, Dore K, Grimsrud KM, Leedell D, Martin SW, Michel P, Svenson LW, McEwen SA: 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. Epidemiol Infect. 2006, 134 (4): 699-711. 10.1017/S0950268805005741.
National Infectious Disease Register (NIDR). [http://www3.ktl.fi/stat]
Hancock DD, Besser TE, Rice DH, Ebel ED, Herriott DE, Carpenter LV: Multiple sources of Escherichia coli O157 in feedlots and dairy farms in the northwestern USA. Prev Vet Med. 1998, 35 (1): 11-19. 10.1016/S0167-5877(98)00050-6.
Heinikainen S, Pohjanvirta T, Eklund M, Siitonen A, Pelkonen S: Tracing shigatoxigenic Escherichia coli O103, O145, and O174 infections from farm residents to cattle. J Clin Microbiol. 2007, 45 (11): 3817-3820. 10.1128/JCM.00198-07.
Lahti E, Eklund M, Ruutu P, Siitonen A, Rantala L, Nuorti P, Honkanen-Buzalski T: Use of phenotyping and genotyping to verify transmission of Escherichia coli O157:H7 from dairy farms. Eur J Clin Microbiol Infect Dis. 2002, 21 (3): 189-195. 10.1007/s10096-001-0682-0.
Parry SM, Salmon RL, Willshaw GA, Cheasty T: Risk factors for and prevention of sporadic infections with vero cytotoxin (shiga toxin) producing Escherichia coli O157. Lancet. 1998, 351 (9108): 1019-1022. 10.1016/S0140-6736(97)08376-1.
Locking M, Allison L, Rae L, Pollock K, Hanson M: VTEC infections and livestock-related exposures in Scotland, 2004. Euro Surveill. 2006, 11 (2): E060223-060224.
O'Brien SJ, Adak GK, Gilham C: Contact with farming environment as a major risk factor for Shiga toxin (Vero cytotoxin)-producing Escherichia coli O157 infection in humans. Emerg Infect Dis. 2001, 7 (6): 1049-1051. 10.3201/eid0706.010626.
Michel P, Wilson JB, Martin SW, Clarke RC, McEwen SA, Gyles CL: Temporal and geographical distributions of reported cases of Escherichia coli O157:H7 infection in Ontario. Epidemiol Infect. 1999, 122 (2): 193-200. 10.1017/S0950268899002083.
Valcour JE, Michel P, McEwen SA, Wilson JB: Associations between indicators of livestock farming intensity and incidence of human Shiga toxin-producing Escherichia coli infection. Emerg Infect Dis. 2002, 8 (3): 252-257. 10.3201/eid0803.010159.
Kistemann T, Zimmer S, Vagsholm I, Andersson Y: GIS-supported investigation of human EHEC and cattle VTEC O157 infections in Sweden: geographical distribution, spatial variation and possible risk factors. Epidemiol Infect. 2004, 132 (3): 495-505. 10.1017/S0950268803001729.
Innocent GT, Mellor DJ, McEwen SA, Reilly WJ, Smallwood J, Locking ME, Shaw DJ, Michel P, Taylor DJ, Steele WB, et al: Spatial and temporal epidemiology of sporadic human cases of Escherichia coli O157 in Scotland, 1996-1999. Epidemiol Infect. 2005, 133 (6): 1033-1041. 10.1017/S0950268805003687.
Haus-Cheymol R, Espie E, Che D, Vaillant V, H DEV, Desenclos JC: Association between indicators of cattle density and incidence of paediatric haemolytic-uraemic syndrome (HUS) in children under 15 years of age in France between 1996 and 2001: an ecological study. Epidemiol Infect. 2006, 134 (4): 712-718. 10.1017/S095026880500542X.
Pebody RG, Furtado C, Rojas A, McCarthy N, Nylen G, Ruutu P, Leino T, Chalmers R, de Jong B, Donnelly M, Fisher I, Gilham C, Graverson L, Cheasty T, Willshaw G, Navarro M, Salmon R, Leinikki P, Wall P, Bartlett C: An international outbreak of Vero cytotoxin-producing Escherichia coli O157 infection amongst tourists; a challenge for the European infectious disease surveillance network. Epidemiol Infect. 1999, 123 (2): 217-223. 10.1017/S0950268899002940.
Hadler JL, Clogher P, Hurd S, Phan Q, Mandour M, Bemis K, Marcus R: Ten-Year Trends and Risk Factors for Non-O157 Shiga Toxin-Producing Escherichia coli Found Through Shiga Toxin Testing, Connecticut, 2000-2009. Clin Infect Dis. 2011, 53 (3): 269-276. 10.1093/cid/cir377.
Paunio M, Pebody R, Keskimaki M, Kokki M, Ruutu P, Oinonen S, Vuotari V, Siitonen A, Lahti E, Leinikki P: Swimming-associated outbreak of Escherichia coli O157:H7. Epidemiol Infect. 1999, 122 (1): 1-5. 10.1017/S0950268898001927.
Olsen SJ, Miller G, Breuer T, Kennedy M, Higgins C, Walford J, McKee G, Fox K, Bibb W, Mead P: A waterborne outbreak of Escherichia coli O157:H7 infections and hemolytic uremic syndrome: implications for rural water systems. Emerg Infect Dis. 2002, 8 (4): 370-375. 10.3201/eid0804.000218.
Hilbe JM: Negative binomial regression - understanding and modeling overdispersed count data. 2011, Cambridge University Press, 2
Ridout M, Demétrio CGB, Hinde J: Models for count data with many zeros. International Biometric Conference. 1998, Cape Town
Berrang-Ford L, Berke O, Sweeney S, Abdelrahman L: Sleeping sickness in southeastern Uganda: a spatio-temporal analysis of disease risk, 1970-2003. Vector Borne Zoonotic Dis. 2010, 10 (10): 977-988. 10.1089/vbz.2008.0196.
Carrel M, Voss P, Streatfield PK, Yunus M, Emch M: Protection from annual flooding is correlated with increased cholera prevalence in Bangladesh: a zero-inflated regression analysis. Environ Health. 2010, 9: 13-10.1186/1476-069X-9-13.
Gonzales-Barron U, Kerr M, Sheridan JJ, Butler F: Count data distributions and their zero-modified equivalents as a framework for modelling microbial data with a relatively high occurrence of zero counts. Int J Food Microbiol. 2010, 136 (3): 268-277. 10.1016/j.ijfoodmicro.2009.10.016.
Spiegelharter D, Thomas A, Best N: WinBUGS Version 1.4. User Manual. MRC Biostatistics Unit. 1999
Besag J, York J, Mollié A: Bayesian image restoration, with two applications in spatial statistics. Ann Inst Statist Math. 1991, 43 (1): 1-59. 10.1007/BF00116466.
Lunn DJ, Thomas A, Best N, Spiegelharter D: WinBUGS -- a Bayesian modelling framework: concepts, structure, and extensibility. Statistics and computing. 2000, 10: 325-337. 10.1023/A:1008929526011.
Congdon P: Bayesian models for categorical data. 2005, John Wiley & Sons Ltd
Ntzoufras I: Bayesian modelling using WinBUGS. 2009, Wiley
Erickson MC, Doyle MP: Food as a vehicle for transmission of Shiga toxin-producing Escherichia coli. J Food Prot. 2007, 70 (10): 2426-2449.
Chang M, Groseclose SL, Zaidi AA, Braden CR: An ecological analysis of sociodemographic factors associated with the incidence of salmonellosis, shigellosis, and E. coli O157:H7 infections in US counties. Epidemiol Infect. 2009, 137 (6): 810-820. 10.1017/S0950268808001477.
Lahti E, Keskimaki M, Rantala L, Hyvonen P, Siitonen A, Honkanen-Buzalski T: Occurrence of Escherichia coli O157 in Finnish cattle. Vet Microbiol. 2001, 79 (3): 239-251. 10.1016/S0378-1135(00)00355-2.
Boqvist S, Aspan A, Eriksson E: Prevalence of verotoxigenic Escherichia coli O157:H7 in fecal and ear samples from slaughtered cattle in Sweden. J Food Prot. 2009, 72 (8): 1709-1712.
Lahti E, Ruoho O, Rantala L, Hanninen ML, Honkanen-Buzalski T: Longitudinal study of Escherichia coli O157 in a cattle finishing unit. Appl Environ Microbiol. 2003, 69 (1): 554-561. 10.1128/AEM.69.1.554-561.2003.
Finland's environmental administration. [http://www.ymparistokeskus.fi/default.asp?node=4032&lan=en]
Garvey P, McKeown P, Carroll A, McNamara E: Epidemiology of Verotoxigenic E. coli in Irland, 2009. EpiInsight. 2010, 11 (9): 1-9.
Gautam R, Bani-Yaghoub M, Neill WH, Dopfer D, Kaspar C, Ivanek R: Modeling the effect of seasonal variation in ambient temperature on the transmission dynamics of a pathogen with a free-living stage: example of Escherichia coli O157:H7 in a dairy herd. Prev Vet Med. 2011, 102 (1): 10-21. 10.1016/j.prevetmed.2011.06.008.
Nygard K, Andersson Y, Rottingen JA, Svensson A, Lindback J, Kistemann T, Giesecke J: Association between environmental risk factors and campylobacter infections in Sweden. Epidemiol Infect. 2004, 132 (2): 317-325. 10.1017/S0950268803001900.
Coia JE, Sharp JCM, Curnow J, Reilly WJ: Ten years of experience of Escherichia coli O157 in Scotland (1984-1993). Recent advances in Verocytotoxin-producing Escherichia Coli O157 in Scotland (1984-1993). Edited by: Karmali MA, Goglio AG. 1994, Elsevier Science B.V, 41-44.
Simonsen J, Frisch M, Ethelberg S: Socioeconomic risk factors for bacterial gastrointestinal infections. Epidemiology. 2008, 19 (2): 282-290. 10.1097/EDE.0b013e3181633c19.
Patil SR, Cates S, Morales R: Consumer food safety knowledge, practices, and demographic differences: findings from a meta-analysis. J Food Prot. 2005, 68 (9): 1884-1894.
Cheasty T, Allerberger F, Beutin L, Caprioli A, Heuvelink A, Karch H, Lofdahl S, Pierard D, Scheutz F, Siitonen A, et al: VTEC O157 phage types isolated from 15 European countries: 1997-1999. 4th International Symposium and Workshop on "Shiga toxin (verocytotoxin) - producing Escherichia coli infections", VTEC 2000: October 29 - November 2, 2000 2000; Kyoto, Japan. 2000, 126-
Feng PC, Monday SR, Lacher DW, Allison L, Siitonen A, Keys C, Eklund M, Nagano H, Karch H, Keen J, Whittam TS: Genetic diversity among clonal lineages within Escherichia coli O157:H7 stepwise evolutionary model. Emerg Infect Dis. 2007, 13 (11): 1701-1706.
Reilly A: Prevention and control of enterohaemorrhagic Escherichia coli (EHEC) infections: memorandum from a WHO meeting. WHO Consultation on Prevention and Control of Enterohaemorrhagic Escherichia coli (EHEC) Infections. Bull World Health Organ. 1998, 76 (3): 245-255.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/11/275/prepub
We thank Raili Ronkainen for assistance with the data handling and Tarja Heiskanen, Aino Kyyhkynen and Anna Liimatainen for their laboratory assistance.
The authors declare that they have no competing interests.
KJ performed the analysis of the descriptive epidemiology data and the statistical analysis and drafted the manuscript. JO designed the statistical analysis and helped to perform it. ME performed the microbiological analysis and participated in the analysis of the descriptive epidemiology. AS participated in the design and coordination of the microbiological analysis. MK participated in the design and coordination of the epidemiological analysis. All authors read and approved the final manuscript.