Email updates

Keep up to date with the latest news and content from BMC Infectious Diseases and BioMed Central.

Open Access Highly Accessed Research article

Environmental and demographic risk factors for campylobacteriosis: do various geographical scales tell the same story?

Julie Arsenault15*, Olaf Berke3, Pascal Michel25, André Ravel25 and Pierre Gosselin45

Author affiliations

1 Faculté de médecine vétérinaire, Université de Montréal, Saint-Hyacinthe, Québec, Canada

2 Laboratoire de lutte contre les zoonoses d’origine alimentaire, Agence de la santé publique du Canada, Saint-Hyacinthe, Québec, Canada

3 Department of Population Medicine, University of Guelph, Guelph, Ontario, Canada

4 Institut national de santé publique du Québec (INSPQ), Beauport, Québec, Canada, Centre hospitalier universitaire de Québec (CHUQ), Sainte-Foy, Québec, Canada

5 Groupe de recherche en épidémiologie des zoonoses et santé publique, Faculté de médecine vétérinaire, Saint-Hyacinthe, Québec, Canada

For all author emails, please log on.

Citation and License

BMC Infectious Diseases 2012, 12:318  doi:10.1186/1471-2334-12-318


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2334/12/318


Received:18 November 2011
Accepted:15 November 2012
Published:22 November 2012

© 2012 Arsenault et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Campylobacter is a common cause of bacterial gastro-enteritis characterized by multiple environmental sources and transmission pathways. Ecological studies can be used to reveal important regional characteristics linked to campylobacteriosis risk, but their results can be influenced by the choice of geographical units of analysis. This study was undertaken to compare the associations between the incidence of campylobacteriosis in Quebec, Canada and various environmental characteristics using seven different sets of geographical units.

Methods

For each set of geographical unit, a conditional autoregressive model was used to model the incidence of reported cases of campylobacteriosis according to environmental (poultry density, ruminant density, slaughterhouse presence, temperature, and precipitation) and demographic (population density, level of education) characteristics. Models were compared in terms of number of significant predictors, differences in direction and magnitude of predictors, and fit of the models.

Results

In general, the number of significant predictors was reduced as the aggregation level increased. More aggregated scales tend to show larger but less precise estimates for all variables, with the exception of slaughterhouse presence. Regional characteristics associated with an increased regional risk of campylobacteriosis, for at least some geographical units, were high ruminant density, high poultry density, high population density, and presence of a large poultry slaughterhouse, whereas a reduction in risk was associated with a lower percentage of people with diplomas, a lower level of precipitation, and warmer temperature. Two clusters of elevated residual risk were observed, with different location and size depending on the geographical unit used.

Conclusions

Overall, our results suggest that the use of municipality or census consolidated subdivision were the most optimal scales for studying environmental determinants of campylobacteriosis at a regional level. This study highlights the need for careful selection and analysis of geographical units when using ecological study designs.

Keywords:
Campylobacteriosis; Quebec; Spatial; Poultry; Ruminants; Climate; Modifiable areal unit problem

Background

In Canada, as in many industrialized countries, infection with Campylobacter spp. is a leading cause of bacterial gastro-enteritis, with an annual average of 39 cases reported per 100,000 people over the last decade [1]. Many case–control studies have been conducted to identify risk factors for campylobacteriosis and have consistently revealed that some factors, including foreign travel, consumption of raw milk, eating in a restaurant, barbecuing, and coming into contact with raw poultry meat are associated with higher risk of disease [2]. Consumption of contaminated food is believed to explain approximately half of the reported cases [3-5]. In the last decade, environmental hypotheses have been put forward to explain campylobacteriosis risks that are not directly attributed to food [5,6]. Wild birds, poultry, sheep and cattle seem to be of particular importance in the natural life cycle of the bacteria, because they are easily infected and can excrete the bacteria in high numbers into the environment [7-14]. Aquatic environments are often contaminated by Campylobacter, which could make recreational and drinking water significant pathways of transmission between animals and humans [15-17]. Domestic flies could act as mechanical vectors for the transmission of the bacteria to humans as well [5]. Finally, meteorological factors likely affect the survival of the bacteria, and have been reported to influence the risk of disease [6,18-20].

Because of the regional-level intrinsic influence of environmental characteristics, the use of an ecological study design (i.e. the unit of analysis is a region) is a valid choice when studying these factors in relation to campylobacteriosis occurrence in populations. In the past, few ecological studies had been conducted to identify regional factors associated with the incidence of human campylobacteriosis. In Sweden, a positive association was found between the incidence of campylobacteriosis and both ruminant density and drinking water quality measures [17]. In the United States, counties with high poultry densities were reported to have higher incidence rates of campylobacteriosis [21,22]. In the province of Manitoba, Canada, the incidence of campylobacteriosis was reported to be highest among populations living in areas with high densities of farm animals, including cows, pigs, and chickens [23]. Although ecological studies such as these may reveal important characteristics linked to campylobacteriosis incidence, their conclusions may be influenced by the geographical scales used in the analysis, an issue known as the “modifiable areal unit problem” [24,25]. Results may also be biased if the boundaries of the geographical units of analysis do not follow the ones at which the process under study operates, an issue known as the zoning effect [24]. Although theoretical and empirical works have been published to better understand these design issues [25-29], we are not aware of any study directly addressing these issues for infectious diseases with environmental reservoirs. Exploring the impact of geographical scale on epidemiological inferences is also a valuable approach to better identify the various scales involved in the process [30].

The objective of this study was to estimate the effect of selected environmental and demographic characteristics on the regional incidence of human campylobacteriosis in Quebec using different geographical segregations of the study area, and to assess the impact of the choice of geographical units on epidemiological inference.

Methods

An ecological study was conducted using human cases of campylobacteriosis reported in the province of Quebec, Canada, between 1996 and 2006 inclusively. The area under study was defined as the populated areas of the province, with the exclusion of non-organized territories, incompletely enumerated Indian reserves and settlements, and northern areas (Nunavik). Populated areas were defined as regions covered by census blocks in which at least one person was living according to the 2001 census of Statistics Canada.

Data acquisition on campylobacteriosis cases

Following approval of the research protocol by the research ethics board of the Faculty of Medicine at the University of Montreal and by the research ethics board of the Agency of Health and Social Services of Montreal, available data on all laboratory-confirmed cases of Campylobacter infections reported in the province of Quebec for the study period were obtained from local health authorities. Cases were geocoded to the full 6-digit postal code areas, municipalities, and local community service center (CLSC). Recurrent cases occurring within less than 5 years of the first episode were excluded [31].

Selection of geographical frameworks and scales

For the purpose of our study, we defined a geographical framework as a set of boundaries delineating an administrative or natural organization of the study area. These frameworks may include different subsets at various scales, which relates to the “aggregation effect” of the modifiable areal unit problem [24]. The different areas defined by the boundary of a geographical framework at a defined scale were termed geographical units, and refer to the notion of “zoning” effects [24]. Geographical sets of units used in this study are presented and defined in Table 1. These units were chosen among those previously described as applicable and of potential interest for the ecological study of campylobacteriosis in Quebec [32].

Table 1. Geographical units

Definition and acquisition of environmental and demographic variables

Demographic and environmental variables were selected based on current literature and availability of data. They are defined in Table 2. Data on animal productions for the years 1998, 2001 and 2004 were obtained for each registered farm from the Ministry of Agriculture, Fisheries and Food of Quebec. Data were geocoded at the centroid of the main production site. For the calculation of ruminant density, only farms with enough ruminants to be registered—estimated to be 15 for small ruminants and 10 for bovine herds—were considered. The same rule was applied to the poultry density variable, with a minimum of 57 to 2100 birds required on the farms (depending on the type of flock) for inclusion. The total number of ruminants and poultry were then averaged over the 3 years of data using the weights 4, 3, and 4, respectively, before calculating densities by km2 of populated areas.

Table 2. Definition of risk factors

Slaughterhouse data for the year 2006 were obtained from the websites of the Canadian Food Inspection Agency and the Ministry of Agriculture, Fisheries and Food of Quebec. Slaughterhouses were geocoded using their civic number and street address within GeoPinPoint 2008 software [33]; manual geocoding at the street level was done whenever the automatic procedure was unsuccessful.

Demographic information was obtained for the census years 1996, 2001, and 2006 from Statistics Canada at the Dissemination Area level (i.e. smallest hierarchical division of the territory at which socio-economic data are released, comprising an average of 400 people). For each census year, data from each Dissemination Area within a geographical unit were averaged, weighted by their population size. For each geographical unit, data were then averaged over the 3 census years using weights of 2.5, 5, and 2.5, respectively.

Climate variables were obtained from the National Land and Water Information Service of Agriculture and Agri-Food Canada for 1996 to 2003, inclusively. These values represent a 10x10 km cell size, generated by the interpolation (ANUSPLIN v4 algorithm) of meteorological data collected from weather stations across Canada. The rasters were first transformed at a cell size of 0.5 km using the Spatial Analyst extension of ArcInfo 9.2 [34], and then average values of raster cells within each geographical unit were computed. This was necessary as some smaller units did not cover any pixel center of the original raster, creating missing values when extracting data.

Statistical modeling

All covariates were categorized prior to modeling in order to get a constant set of predictors for all models and/or because the variable to outcome relationship was assumed to be non-linear in form. Animal density variables were put into 3 categories: absence, medium levels, and high levels of production. The cut-off for high levels was set at the 85th percentile of the respective empirical distribution at the level of municipalities. The slaughterhouse variable was organized into 3 categories: absence of any kind of slaughterhouse, presence of large poultry slaughterhouse(s), and presence of small poultry slaughterhouse(s) or slaughterhouse(s) for other species. A poultry slaughterhouse belonging to one of the four largest poultry transformation companies in Quebec, which have a combined market share of over 90%, was defined as “large” [35]. For demographic and climatic variables, 3 categories were defined using the 15th and 85th percentile of their empirical distributions at the level of municipalities.

For statistical modeling, the dependent variable was defined as the average number of reported cases of campylobacteriosis per 100,000 people per year. The incidence was directly standardized for age group (0–4 yrs, 5–15 yrs, 16–44 yrs, >45 yrs), using the study population for the year 2001 as the standard population. Prior to modeling, the incidence rate was smoothed using the empirical Bayesian estimation in order to improve normality and variance homogeneity [36]. This procedure also has the advantage of avoiding the use of transformed scales, which can be difficult to interpret [37]. Ordinary regression models were built for each geographical unit without consideration of the spatial structure of the data. All covariates were included as fixed effects, with no interactions. No variable selection was done in order to avoid any influence of the selection method or covariate(s) exclusion in the comparison of models. Models were built using SAS software [38]. Studentized residuals were computed from the final models to assess the spatial dependence. An empirical semi-variogram was estimated based on Euclidean distances, with a 95% confidence band estimated in R (package geoR [39]). The choice of Euclidean distances to evaluate the spatial dependence in residuals was based on preliminary analyses (see Additional file 1).

Additional file 1. Choice of distance matrix.

Format: DOC Size: 190KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

In the presence of spatial structure in studentized residuals from the ordinary regression models, a conditional autoregressive (CAR) model was fit in R using the “spdep” package [40] for each set of geographical units. A binary neighbor matrix was used to account for spatial autocorrelation. The cut-off within which surrounding areas were considered neighbors was determined as the practical range of the semi-variogram using studentized residuals from ordinary regression models. The practical range was estimated by maximum likelihood in R using the package “geoR” [39], based on an exponential model. All variables were included as fixed effects. Residuals from the final CAR models were visually assessed for normality using QQ-plots.

The geographical sets of units were first compared based on the estimated associations between outcome and exposure variables. More specifically, the number of exposure variables with statistically significant associations was compared between sets of units. Among significant variables, differences in direction and magnitude in risk estimates were also described. For the comparisons in terms of magnitude, the maximal percentage of variation between statistically significant coefficient estimates for each explanatory variable was calculated as follows:

<a onClick="popup('http://www.biomedcentral.com/1471-2334/12/318/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2334/12/318/mathml/M1">View MathML</a>

The geographical sets of units were then compared in term of model fitting. Three criteria were used: 1) Pearson’s coefficient of correlation (r2) between observed and predicted outcome (a good fit is indicated by strong correlation); 2) spatial dependence in model residuals as measured by Moran’s I using the same neighborhood definition as used in the conditional autoregressive models; and 3) presence of areas of unexplained higher risk. This was evaluated with a spatial scan test scanning for primary and secondary clusters of high residual values based on the normal model, performed in SaTScan version 8.1.1. [41]. Significant clusters (p<0.05) were mapped.

Results

A total of 28,521 cases of campylobacteriosis reported in the study area between 1996 and 2006 were included in the analysis. The average annual population within the study area was 7,367,517 people, giving an overall estimate of annual incidence of campylobacteriosis of 35.2 cases per 100,000 people. The size of the study area was 218 668 km2.

The distribution of exposure variables included in regression models is presented in Table 3. For all ordinary regression models, there was evidence of spatial dependence among studentized residuals. Thus, CAR models were built for all sets of geographical units. Results from CAR models are presented in Table 4. Neighbors were defined as units with population centers within the practical range estimated from the semi-variogram, ranging from 15 km to 34 km depending on the model (see Table 5). Residuals from all CAR models were normally distributed according to visual inspection.

Table 3. Distribution of risk factors

Table 4. Regression coefficients of CAR models

Table 5. Cluster and clustering in residuals

The number of exposure variables significantly associated with campylobacteriosis incidence in CAR models tended to decrease as the level of aggregation increased (Table 4). For all variables, the direction of the association for statistically significant estimates was consistent across all sets of geographical units. A large variability was detected in point estimates of statistically significant variables obtained from different geographical units, ranging from 16% to 453% (Table 4). In general, the estimates of regression coefficients increased with a high level of aggregation, with the exception of slaughterhouses for which a reverse trend was noted. The “poultry density” was particularly variable, and this high degree of variation was driven by the coefficients estimated at the watershed level, which were three times higher compared to others.

The Pearson r2 between the observed and predicted values of the CAR models ranged from 21% to 52%, with an increasing trend as the data got more aggregated (Table 4). No clustering was detected in residuals according to Moran’s I test (Table 5). The scan test detected the same small hotspot area in residuals for the three sets of geographical units at a smaller scale, whereas a very large hotspot was identified for two of the geographical units at a larger scale (Table 5 and Figure 1). This large area was located in a zone with a high predicted incidence of campylobacteriosis (see Figure 2). No secondary cluster was found to be significant (all p>0.05).

thumbnailFigure 1. A-G. Distribution of residuals. Residuals from a CAR model predicting the annual incidence of reported cases of campylobacteriosis per 100,000 people in Quebec (1996–2006) according to various geographical units. Significant clusters (p<0.05) in residuals according to the scan test are illustrated. Classification was done using Jenk’s natural breaks. Dark grey areas represent the unpopulated areas, non-organized territories or incompletely enumerated Indian reserves and settlements within Quebec, whereas light grey shows the frontier area of Quebec.

thumbnailFigure 2. A-B. Predicted incidence of campylobacteriosis. Predicted annual incidence of reported cases of campylobacteriosis per 100,000 people in Quebec (1996–2006) according to a conditional autoregressive (CAR) model at the level of municipality (A) or census division (B). Classification was done using Jenk’s natural breaks. Dark grey areas represent the unpopulated areas, non-organized territories or incompletely enumerated Indian reserves and settlements within Quebec, whereas light grey shows the frontier area of Quebec.

Discussion

The various geographical sets of units gave various insights into the spatial distribution of campylobacteriosis incidence. These differences can be attributed to a change in scale (e.g. for administrative units) or to a combination of changes in scale and zoning. Different criteria were used to compare the results obtained from different sets of geographical units, selected for their epidemiological relevance. First, the difference in estimated associations between campylobacteriosis incidence and exposure variables was explored, since different results can lead to a different understanding of the environmental factors influencing the spatial distribution of the disease. Next, residual clustering or hotspots of unexplained risk were explored. Local zones of unexplained risk were detected, which potentially identify different scales at which disease processes are occurring and specific areas and scales warranting further investigation.

Associations between campylobacteriosis and exposure variables

The number of significant variables in the CAR models was reduced as the level of aggregation increased. This was likely driven by a reduction of statistical power due to lower sample size rather than reduced strength of association as the point estimate in regression coefficients were in general larger in more aggregated data. It should be noted that, with the exception of covariates referring to animal density, significant effects on regression coefficients were only observed for those covariates showing less variability at lower level geographical units [32]. Choice of unit based on minimal intra-zonal variance was suggested as a potential solution to the modifiable areal unit problem, as recently reviewed [25]. The agriculture-based custom framework was created in order to increase the size of units to improve stability in rates, while preserving high intra-unit homogeneity in animal production. However, it has the drawback of being heterogeneous in unit size and many units had a non-compact shape, which might have hampered the ability to find statistically significant associations.

The direction of the associations was consistent and in agreement with current biological knowledge. Poultry and ruminants are frequently colonized with Campylobacter and can shed the bacteria in high numbers [10,42-44]. Contamination in humans can then occur through direct contact with these animals or following indirect contact with a contaminated environment. We thus expected a positive association with a likely dose–response relationship, as we observed. For climate variables, the directionality of the association was also compatible with Campylobacter biology, which has reduced survival in a warmer environment and is sensitive to desiccation [45-49]. The inclusion of education as predictor variable was justified by a previous study reporting its influence on some steps needed for a case to be reported [50], with potential bias on spatial patterns. We observed a reduction in the incidence of campylobacteriosis with a lower level of education, which has also been reported by others and might represent a proxy for socio-economic status [51]. Traveling and eating in a restaurant were reported as risk factors for campylobacteriosis in other countries, which are also associated with socio-economic variables including the level of education [51,52]. To our knowledge, the importance of traveling as a risk factor for campylobacteriosis has not been evaluated in Canada, but is supported by a study conducted in a Canadian community, where about 22% of reported cases of campylobacteriosis were associated with international traveling [53]. We could not exclude travel-related cases from our study because this information was missing for most cases. The biological pathway linking slaughterhouse presence to campylobacteriosis incidence is twofold: professional exposure of poultry slaughterhouse workers was previously reported to increase the risk of campylobacteriosis [54]; and slaughterhouse effluents were reported to harbor large quantities of Campylobacter [55], which could then contaminate the surrounding drinking or recreational water. However, we do not have a clear explanation for positive association between campylobacteriosis incidence and population density observed for some geographical units, apart from some effects due to differences in case reporting.

The slaughterhouse variable was associated with campylobacteriosis incidence, with a significant decrease in magnitude with aggregation and non-overlapping confidence intervals between geographical units. This met our expectation for this particular variable, since its effect would likely be limited to the nearby areas where workers live or environmental contamination occurs, and a reduction of the coefficient in the aggregation process would occur due to a dilution effect. We further explored these associations by multiplying the point estimate of large poultry slaughterhouse coefficients with the total populations at risk living in areas with large poultry slaughterhouses. In summary, the case load due to the presence of large poultry slaughterhouses in the study area ranged from 9 to 172 cases, depending on the framework. This illustrates the impact of choice of geographical units on potential public health conclusions.

For all other variables, the estimated coefficients of regression tend to increase in magnitude as the units become more aggregated, but this increase was not consistent across all geographical sets of units and thus not predictable. This increase is in accordance with previous work based on simple linear regression models, although another study conducted in a multivariate regression setting reported that the results were essentially unpredictable [29]. For the animal density variables, the estimated regression coefficients with the highest values also had large standard errors, and were not statistically significant. Among the statistically significant ones, the variation in point estimates was too small in our perspective to affect epidemiological conclusions on the biological relevance of these variables. The only exception was perhaps for poultry density at the watershed level. At this level of aggregation, we suspected the presence of colinearity between poultry density and ruminant density based on our knowledge of the study area. Also, according to a simulation study, the aggregation process can introduce colinearity between exposure variables where none existed before [26]. In order to explore this possibility, alternative models were built by excluding either the poultry density or the ruminant density from the model. For the watershed units, the exclusion of the poultry density variable led to larger and significant coefficients for the ruminant density variable (from 16.0 [p=0.22] to 23.1 [p=0.08] for ruminant density >20 per km2 and from 8.9 [p=0.12] to 15.9 [p<0.01] for ruminant density ≤250 per km2). This is suggestive of colinearity between these two variables, and so the effect attributed to poultry or ruminant density for frameworks at a larger scale is probably a cumulative effect of both types of production. Indeed, any epidemiological study aiming to study the distinct effect of poultry and ruminant density should avoid any use of geographical units at a scale larger than census consolidated subdivisions. From another perspective, the large regression coefficients for animal densities observed at the watershed scale could be indicative of real biological processes occurring at that specific scale, which is not impossible considering the importance of aquatic environments in the ecology of Campylobacter [56].

For the demographic and climatic variables, the point estimates of regression coefficients also increased with aggregation, although confidence intervals increased and point estimates were not statistically different between geographical units (i.e. point estimates for one geographical unit were included in the confidence intervals of estimates from other geographical units). It should be noted that high values of estimated regression coefficients for climate variables as seen for the “smallest” and “CLSC” geographical units compared to “municipality” are probably misleading. In fact, 51 and 39 of the divisions constituting these units were located on the islands of Montreal and Laval, respectively, urban areas characterized by warm temperatures and low amounts of precipitation. At the scale where it was measured, almost no variation in weather variables was seen within this area; indeed, these areas were probably over-represented in the estimated relationship between campylobacteriosis and weather variables.

Evaluation of the fit of models

The fit of models was explored in different ways. As previously reported, an increase of the correlation coefficient between observed and predicted data was seen as the level of aggregation increased. This is attributed to the smoothing effect of aggregation [25,29] and highlights the fact that too much emphasis should not be placed on the interpretation of this measure. The absence of clustering in residuals was suggestive of a good adjustment for spatial dependence in all models. And finally, a small hotspot of unexplained residual risk was detected by the scan test for the three sets of geographical units at the smallest scales, i.e. the “smallest,” “municipality,” and “census consolidated subdivision” geographical units. Considering the small size of this hotspot, it was likely smoothed out in the aggregation process for geographical units at a large scale. A community episode of campylobacteriosis from a common source is a potential explanation for this high risk zone, or it could have been caused by a locally-acting factor not measured in our study. Molecular typing data on Campylobacter strains were not available from our dataset, but would have been very useful in distinguishing between these two processes, as we expect a higher similarity of strains isolated within a cluster in the presence of an outbreak from a common source [57]. For two of the geographical units at a large scale (i.e. CLSC and census division), a large area of significant unexplained risk was detected. This area was characterized by high poultry production and relatively high predicted risk of campylobacteriosis. We do not have a clear explanation for this finding. It is possible that in areas of high animal density, the effect of environmental variables on campylobacteriosis risk is amplified in a non-linear way due to the dynamic of bacterial transmission between reservoirs. In more aggregated units, the areas included in the average animal density categories were less homogeneous, which might lead to poorer adjustment of the model. Geographically-weighted regression might be useful in exploring this possibility, but this was beyond the scope of our study.

Limitations

The choice of framework and scales included in the study were limited to the ones at which cases could be located. Additionally, only cases reported through the surveillance system were included, limiting the generalization of results to all cases occurring in the population. In fact, case reporting involves multiple steps, which could be influenced by health care system, physician, or patient characteristics. In the province of Quebec, health care services are accessible to all, with universal health insurance plan allowing the entire population to receive free hospital and medical services, which should limit any related reporting bias. Regarding physicians, their propensity to ask for a stool sample have been reported to vary according to their speciality and severity of illness of their patients [58]. Also, the decision of a diarrheic patient to consult a physician has been reported to depend upon severity of illness, urbanicity, and various socio-economic indicators of patients [50,59]. Our models only offer partial adjustment related to patient characteristics, with the inclusion of an educational status variable. Thus, we cannot exclude that the observed association with environmental risk factors were biased by spatially-varying factors influencing steps conducting to the reporting of a case. However, as we used the same dataset for all geographical sets of units, the underreporting bias was kept constant and this should not invalidate our conclusions regarding the influence of geographical units on risk factor estimation. Also, our results could be dependent on the method used for statistical modeling. Using the Poisson regression model with spatially correlated random effects was explored for our study, but no convergence was obtained for most models. Our study design did not allow us to point out which mechanisms were exactly responsible for the differences observed. Theoretical studies on modifiable areal units assume that unit size does not vary much, but in our case we had small urban units with large rural units in most frameworks, and this could complicate interpretation.

Conclusions

Our study allowed us to evaluate to what extent the modifiable areal unit problem influences the epidemiological conclusions on the spatial distribution of campylobacteriosis. In general, the various geographical units gave consistent results for the direction of the associations between incidence and various environmental and demographic characteristics, which were in agreement with current knowledge on campylobacteriosis epidemiology. However, the strength and statistical significance of the associations observed at one scale should be used with caution when inference is made at another scale. In general, municipality or census consolidated subdivision scales are recommended for the study of the spatial distribution of campylobacteriosis in our study area, based on theoretical grounds [32] and low degree of colinearity in agriculture variables at this scale. These were also associated with significant and biologically meaningful associations between campylobacteriosis incidence and the exposure variables under study.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JA designed the methodology, performed the statistical analyses and drafted the manuscript. OB supervised the analyses. PM, AR and PG assisted in designing the methodology and drafting the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was made possible through doctoral research awards from the Canadian Institutes of Health Research and the Faculté de Médecine vétérinaire de l’Université de Montréal attributed to J. Arsenault. We are most grateful to the provincial health units of Quebec, to the Ministère de l’Agriculture, des Pêcheries et de l’Alimentation du Québec, and to the National Land and Water Information Service led by Agriculture and Agri-Food Canada for sharing their data.

References

  1. Public Health Agency of Canada:

    Notifiable Diseases On-Line.

    http://dsol-smed.phac-aspc.gc.ca/dsol-smed/ndis/index_e.html webcite

    OpenURL

  2. Engberg J: Contributions to the epidemiology of Campylobacter infections. A review of clinical and microbiological studies.

    Dan Med Bull 2006, 53(4):361-389. PubMed Abstract OpenURL

  3. Evers EG, Van Der Fels-Klerx HJ, Nauta MJ, Schijven JF, Havelaar AH: Campylobacter source attribution by exposure assessment.

    Int J Risk Assess Manag 2008, 8(1/2):174-190. Publisher Full Text OpenURL

  4. Vellinga A, Van Loock F: The dioxin crisis as experiment to determine poultry-related Campylobacter enteritis.

    Emerg Infect Dis 2002, 8(1):19-22. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Ekdahl K, Normann B, Andersson Y: Could flies explain the elusive epidemiology of campylobacteriosis?

    BMC Infect Dis 2005, 5(1):11. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  6. Skelly C, Weinstein P: Pathogen survival trajectories: an eco-environmental approach to the modeling of human campylobacteriosis ecology.

    Environ Health Perspect 2003, 111(1):19-28. PubMed Abstract | PubMed Central Full Text OpenURL

  7. Abulreesh HH, Paget TA, Goulder R: Campylobacter in waterfowl and aquatic environments: incidence and methods of detection.

    Environ Sci Technol 2006, 40(23):7122-7131. PubMed Abstract | Publisher Full Text OpenURL

  8. Kassa H, Harrington B, Bisesi MS: Risk of occupational exposure to Cryptosporidium, Giardia, and Campylobacter associated wiih the feces of giant Canada geese.

    Appl Occup Environ Hyg 2001, 16(9):905-909. OpenURL

  9. Busato A, Hofer D, Lentze T, Gaillard C, Burnens A: Prevalence and infection risks of zoonotic enteropathogenic bacteria in Swiss cow-calf farms.

    Vet Microbiol 1999, 69(4):251-263. PubMed Abstract | Publisher Full Text OpenURL

  10. Inglis GD, Kalischuk LD, Busz HW: A survey of Campylobacter species shed in faeces of beef cattle using polymerase chain reaction.

    Can J Microbiol 2003, 49(11):655-661. PubMed Abstract | Publisher Full Text OpenURL

  11. Inglis GD, Kalischuk LD, Busz HW, Kastelic JP: Colonization of cattle intestines by Campylobacter jejuni and Campylobacter lanienae.

    Appl Environ Microbiol 2005, 71(9):5145-5153. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Stanley KN, Wallace JS, Currie JE, Diggle PJ, Jones K: Seasonal variation of thermophilic campylobacters in lambs at slaughter.

    J Appl Microbiol 1998, 84(6):1111-1116. PubMed Abstract | Publisher Full Text OpenURL

  13. Jones K, Howard S, Wallace JS: Intermittent shedding of thermophilic campylobacters by sheep at pasture.

    J Appl Microbiol 1999, 86(3):531-536. PubMed Abstract | Publisher Full Text OpenURL

  14. Wallace JS, Stanley KN, Jones K: The seasonal incidence of thermophilic campylobacters in sheep. In ampylobacters, Helicobacters, and Related Organisms. Edited by Newell DG, Ketley JM, Feldman RA. New York: Plenum publishing corporation; 1997:359-362. OpenURL

  15. Adak GK, Cowden JM, Nicholas S, Evans HS: The Public Health Laboratory Service national case–control study of primary indigenous sporadic cases of Campylobacter infection.

    Epidemiol Infect 1995, 115(1):15-22. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Koenraad PMFJ, Rombouts FM, Notermans SHW: Epidemiological aspects of thermophilic Campylobacter in water-related environments: A review.

    Water Environ Res 1997, 69(1):52-63. Publisher Full Text OpenURL

  17. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Fleury M, Charron DF, Holt JD, Allen OB, Maarouf AR: A time series analysis of the relationship of ambient temperature and common bacterial enteric infections in two Canadian provinces.

    Int J Biometeorol 2006, 50(6):385-391. PubMed Abstract | Publisher Full Text OpenURL

  19. Louis VR, Gillespie IA, O’Brien SJ, Russek-Cohen E, Pearson AD, Colwell RR: Temperature-driven Campylobacter seasonality in England and Wales.

    Appl Environ Microbiol 2005, 71(1):85-92. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Tam CC, Rodrigues LC, O’Brien SJ, Hajat S: Temperature dependence of reported Campylobacter infection in England, 1989–1999.

    Epidemiol Infect 2006, 134(1):119-125. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Potter RC, Kaneene JB, Gardiner J: A comparison of Campylobacter jejuni enteritis incidence rates in high- and low-poultry-density counties: Michigan 1992–1999.

    Vector Borne Zoonotic Dis 2002, 2(3):137-143. PubMed Abstract | Publisher Full Text OpenURL

  22. Sandberg M, Nygard K, Meldal H, Valle PS, Kruse H, Skjerve E: Incidence trend and risk factors for Campylobacter infections in humans in Norway.

    BMC Publ Health 2006, 6:179. BioMed Central Full Text OpenURL

  23. Green CG, Krause D, Wylie J: Spatial analysis of Campylobacter infection in the Canadian province of Manitoba.

    Int J Health Geogr 2006, 5(1):2. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Waller LA, Gotway CA: Applied spatial statistics for public health data. Hoboken, New Jersey: John Wiley & Sons Inc; 2004. OpenURL

  25. Wong D: The modifiable areal unit problem (MAUP). In The SAGE handbook of spatial analysis. Edited by Fotheringham AS, Rogerson PA. London: SAGE publications; 2009:105-123. OpenURL

  26. Reynolds HD: The modifiable area unit problem: empirical analysis by statistical simulation. Toronto: University of Toronto; 1998. OpenURL

  27. Gotway CA, Young LJ: Combining incompatible spatial data [Review].

    J Am Stat Assoc 2002, 97(458):632-648. Publisher Full Text OpenURL

  28. Flowerdew R, Manley DJ, Sabel CE: Neighbourhood effects on health: does it matter where you draw the boundaries?

    Soc Sci Med 2008, 66(6):1241-1255. PubMed Abstract | Publisher Full Text OpenURL

  29. Fotheringham AS, Wong DWS: The modifiable areal unit problem in multivariate statistical analysis.

    Environ Plan A 1991, 23:1025-1044. Publisher Full Text OpenURL

  30. Manley D, Flowerdew R, Steel D: Scales, levels and processes: Studying spatial patterns of British census variables.

    Computers, Environments and Urban Systems 2006, 30(2):143-160. Publisher Full Text OpenURL

  31. Arsenault J, Ravel A, Michel P, Berke O, Gosselin P: Do patients with recurrent reported episodes of campylobacteriosis differ from those with a single disease event?

    BMC Publ Health 2011, 11(1):32. OpenURL

  32. Arsenault J: Épidémiologie spatiale de la campylobactériose au Québec [Spatial epidemiology of campylobacteriosis in Quebec]. Saint-Hyacinthe, Qc, Canada: Université de Montréal; 2010. [PhD thesis] OpenURL

  33. DMTI Spatial Inc: GeoPinpoint Suite version 2008.3, release version 6.4. Canada: Markham (Ontario); 2008. OpenURL

  34. ESRI: ArcGIS Desktop 9.3.1 (ArcInfo). Redlands, CA, USA; 2009. OpenURL

  35. Ministère de l’Agriculture des Pêcheries et de l’Alimentation du Québec: Portrait sommaire de l’industrie de la volaille au Québec [Overview of poultry industry in Quebec]. Canada: Gouvernement du Québec: Direction des politiques commerciales et intergouvernementales, Québec; 2006. OpenURL

  36. Berke O: Choropleth mapping of regional count data of Echinococcus multilocularis among red foxes in Lower Saxony, Germany.

    Prev Vet Med 2001, 52(2):119-131. PubMed Abstract | Publisher Full Text OpenURL

  37. Berke O: Exploratory disease mapping: kriging the spatial risk function from regional count data.

    Int J Health Geogr 2004, 3(1):18. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  38. SAS Institute Inc: SAS version 9.2. NC, USA; [Cary] OpenURL

  39. Ribeiro PJJ, Diggle PJ: geoR: a package for geostatistical analysis.

    R-NEWS 2001, 1(2):15-18. OpenURL

  40. Bivand R:

    spdep: Spatial dependence: weighting schemes, statistics and models. R package version 0.5-4. 2010. OpenURL

  41. Kulldorff M, Huang L, Konty K: A scan statistic for continuous data based on the normal probability model.

    Int J Health Geogr 2009, 8:58. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  42. Park SF: The physiology of Campylobacter species and its relevance to their role as foodborne pathogens.

    Int J Food Microbiol 2002, 74(3):177-188. PubMed Abstract | Publisher Full Text OpenURL

  43. Waldenstrom J, Broman T, Carlsson I, Hasselquist D, Achterberg RP, Wagenaar JA, Olsen B: Prevalence of Campylobacter jejuni, Campylobacter lari, and Campylobacter coli in different ecological guilds and taxa of migrating birds.

    Appl Environ Microbiol 2002, 68(12):5911-5917. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Stanley KN, Wallace JS, Jones K: The seasonality of thermophilic campylobacters in beef and dairy cattle. In Campylobacters, Helicobacters, and Related Organisms. Edited by Newell DG, Ketley JM, Feldman RA. New York: Plenum publishing corporation; 1997:163-167. OpenURL

  45. Easton J: Fate and transport of campylobacters in soil arising from farming practices. In Campylobacters, Helicobacters, and Related Organisms. Edited by Newell DG, Ketley JM, Feldman RA. New York: Plenum publishing corporation; 1997:461-465. OpenURL

  46. Korhonen LK, Martikainen PJ: Comparison of the survival of Campylobacter jejuni and Campylobacter coli in culturable form in surface water.

    Can J Microbiol 1991, 37(7):530-533. PubMed Abstract | Publisher Full Text OpenURL

  47. Buswell CM, Herlihy YM, Lawrence LM, McGuiggan JT, Marsh PD, Keevil CW, Leach SA: Extended survival and persistence of Campylobacter spp. in water and aquatic biofilms and their detection by immunofluorescent-antibody and -rRNA staining.

    Appl Environ Microbiol 1998, 64(2):733-741. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Thomas C, Hill DJ, Mabey M: Evaluation of the effect of temperature and nutrients on the survival of Campylobacter spp. in water microcosms.

    J Appl Microbiol 1999, 86(6):1024-1032. PubMed Abstract | Publisher Full Text OpenURL

  49. Fernandez H, Vergara M, Tapia F: Dessication resistance in thermotolerant Campylobacter species.

    Infection 1985, 13(4):197. PubMed Abstract | Publisher Full Text OpenURL

  50. Tam CC, Rodrigues LC, O’Brien SJ: The study of infectious intestinal disease in England: what risk factors for presentation to general practice tell us about potential for selection bias in case–control studies of reported cases of diarrhoea.

    Int J Epidemiol 2003, 32(1):99-105. PubMed Abstract | Publisher Full Text OpenURL

  51. Simonsen J, Frisch M, Ethelberg S: Socioeconomic risk factors for bacterial gastrointestinal infections.

    Epidemiology 2008, 19(2):282-290. PubMed Abstract | Publisher Full Text OpenURL

  52. Unicomb LE, Dalton CB, Gilbert GL, Becker NG, Patel MS: Age-specific risk factors for sporadic Campylobacter infection in regional Australia.

    Foodborne Pathog Dis 2008, 5(1):79-85. PubMed Abstract | Publisher Full Text OpenURL

  53. Ravel A, Nesbitt A, Marshall B, Sittler N, Pollari F: Description and burden of travel-related cases caused by enteropathogens reported in a Canadian community.

    J Travel Med 2011, 18(1):8-19. PubMed Abstract | Publisher Full Text OpenURL

  54. Wilson IG: Airborne Campylobacter infection in a poultry worker: case report and review of the literature.

    Commun Dis Public Health 2004, 7(4):349-353. PubMed Abstract OpenURL

  55. Höller C: Long-term study of occurrence, distribution and reduction of Campylobacter spp. in the sewage system and wastewater treatment plant of a big town.

    Water Sci Technol 1988, 20:529-531. OpenURL

  56. Thomas C, Gibson H, Hill DJ, Mabey M: Campylobacterepidemiology: an aquatic perspective.

    J Appl Microbiol 1999, 85:168S-177S. OpenURL

  57. Clark CG, Taboada E, Grant CC, Blakeston C, Pollari F, Marshall B, Rahn K, Mackinnon J, Daignault D, Pillai D, et al.: Comparison of molecular typing methods useful for detecting clusters of Campylobacter jejuni and C. coli isolates through routine surveillance.

    J Clin Microbiol 2012, 50(3):798-809. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Hennessy TW, Marcus R, Deneen V, Reddy S, Vugia D, Townes J, Bardsley M, Swerdlow D, Angulo FJ: Survey of physician diagnostic practices for patients with acute diarrhea: clinical and public health implications.

    Clin Infect Dis 2004, 38(Suppl 3):S203-S211. PubMed Abstract | Publisher Full Text OpenURL

  59. Herikstad H, Yang S, Van Gilder TJ, Vugia D, Hadler J, Blake P, Deneen V, Shiferaw B, Angulo FJ: A population-based estimate of the burden of diarrhoeal illness in the United States: FoodNet, 1996–7.

    Epidemiol Infect 2002, 129(1):9-17. PubMed Abstract | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2334/12/318/prepub