INTRODUCTION
Campylobacteriosis is the most commonly reported bacterial enteric infection in humans throughout the world [1]. It causes severe abdominal pain, diarrhoea of several days' duration and can lead to Guillain–Barré syndrome, reactive arthritis, irritable bowel syndrome and, rarely, death [Reference Nachamkin, Szymanski and Blaser2]. About 50% of human cases are attributable to the consumption and handling of contaminated poultry; an additional 10–25% of cases are attributable to poultry through direct (e.g. animal contact) and indirect (e.g. drinking water) means [Reference Mullner3–Reference Mughini Gras6]. Recent models propose various complex transmission pathways of Campylobacter from animals and the environment to humans [1, Reference Arsenault7].
In Canada, the province of British Columbia (BC) (population 4·5 million) reported an annual average of 38·0/100000 cases of campylobacteriosis between 2005 and 2009, among the highest in the country [8]. In BC, several epidemiological studies have considered the role of the environment in human Campylobacter infection. One study showed the risk of sporadic campylobacteriosis to be eight times higher in people living in areas serviced by private single-dwelling wells compared to those on municipal water systems [Reference Uhlmann9]. Another study found the risk of campylobacteriosis was higher in urban areas and areas with high poultry density [Reference Laberge, Michel and Galanis10]. Finally, comparing locally acquired and travel-associated Campylobacter infections found higher odds of local infection in males, older people, those living in rural areas or close to Agricultural Land Reserves (ALRs), and in the warmer months [Reference Kuo11]. Some of these results are contradictory, which may be due to the fact that exposures have not yet been considered in combination or sequentially.
We conducted a study to improve our understanding of the non-foodborne transmission pathways of Campylobacter. We studied the combined effects of geographical factors, including drinking water and proximity to agriculture, associated with sporadic human campylobacteriosis in one region of BC, taking into account sociodemographic factors that may bias such associations. We hypothesized that campylobacteriosis is associated with exposure to undertreated drinking water contaminated through agricultural animal sources.
METHODS
A case-case design was used to compare reported cases of campylobacteriosis to other cases of reportable enteric diseases, referred to as enteric disease controls. All laboratory-confirmed cases of reportable enteric disease are reported from local BC laboratories to local health authorities and then to the BC Centre for Disease Control (BCCDC). Epidemiological data transferred on each case to the provincial level include age, sex, address of residence and ‘episode date’ (the earliest of onset, collection or reporting date). We did not have access to interview-based exposure information such as travel status, drinking water consumed or whether cases were part of a known outbreak.
We extracted all cases reported between 2005 and 2009 in the Fraser Health Authority (FHA) from the provincial database. This region encompasses 1·6 million people in an area east of the City of Vancouver and includes urban, semi-rural and rural settings [12]. It is the major poultry-growing region of the province. The drinking water supply includes 19 municipal systems and 450 smaller systems using surface, ground or mixed sources. These municipal and small water systems serve 95% and 1% of the population, respectively. The remaining 4% of the population obtain their drinking water primarily from ground water sources (private wells), for which owners are responsible for testing and treating [13].
In order to determine drinking water at patient's address of residence, we geocoded (assigned geographical coordinates to) cases and drinking water systems. All geocoding, spatial data analysis, and spatial database management were performed using ArcGIS v. 10·0 (ESRI, USA). Cases were geocoded by exact address of residence at the time of illness using a reference land parcel database [Reference Strickland14]. If exact address matching could not be accomplished, street address range interpolation or six-digit postal code geocoding was performed.
A spatial database was created to assign drinking water data to each case, via their address of residence, using spatial overlay and join functions, according to a previous study [Reference Uhlmann9]. Drinking water system management data (i.e. whether the water was on a municipal or a private system) were acquired for every land parcel from the 19 municipalities and two regional districts in the region. We excluded non-residential land parcels. Households that were not on a municipal or specific private system were deemed to have a private well. This was verified using the BC WELLS database [15]. We obtained information on the source of water for each system from FHA. A drinking water index was created taking into account system management (municipal or private) and drinking water source (ground, surface or mixed). Drinking water treatment data were assessed but deemed incomplete and not included in the analysis.
Geographical information systems (GIS)-based analyses were used to calculate proximity to the nearest ALR and nearest abattoir. ALR data were obtained from the Provincial Agricultural Land Commission, and abattoir data from the Canadian Food Inspection Agency and the BCCDC [16, 17]. Cases and controls were classified as living within or further than 2 km from an ALR, as previously validated [Reference Kuo11, Reference Turgeon18]. Distance to the nearest poultry abattoir and the nearest red meat abattoir was calculated for each case and control. An urban environment was defined as having a minimum population concentration of 1000 and a population density of > 400/km2; all other areas were defined as rural [19]. The material component of the 2006 Deprivation Index (DI) for Health Planning, based on education, employment and income indicators, was divided into tertiles and used as the socioeconomic status (SES) measure at the dissemination area census geographical unit [Reference Pampalon20, 21]. We assigned rural/urban environment and SES to each household using GIS-based spatial join functions.
Univariate analyses and multivariate logistic regression were conducted using SAS v. 9·3 (SAS Institute Inc., USA). Univariate comparisons between cases and controls were based on likelihood ratio χ 2 tests (comparing proportions) and t tests allowing unequal variance (comparing continuous measures). Multivariate analyses treated ‘status’ (campylobacteriosis case or control enteric disease) as the binary dependent variable and included water index (surface/municipal, mixed/municipal, ground-well/municipal, ground-well/private), age (0–14 or ⩾15 years), sex (male/female), urban area (yes/no), DI, average distance to the nearest abattoir, distance to nearest ALR (0–2 km or > 2 km) and month as categorical explanatory variables. The best fit model was based on backwards stepwise elimination of non-significant (alpha = 0·05) effects from the full model, which included all main effects and two-way interactions. At each step, goodness-of-fit statistics (Akaike's Information Criterion) were examined to confirm adequacy of model fit. This approach analyses the odds of campylobacteriosis infection relative to all other enteric diseases and assesses statistical significance using the likelihood ratio χ 2 statistic.
Ethics approval was received from the University of British Columbia (certificate no. H11-00430).
RESULTS
Between 2005 and 2009, 2992 campylobacteriosis cases and 4816 enteric disease controls were reported in FHA and included in the study. Enteric disease controls included cases of amoebiasis (n = 603), cryptosporidiosis (n = 209), cyclosporiasis (n = 76), shiga-toxigenic E. coli infection (n = 261), giardiasis (n = 1244), listeriosis (n = 16), salmonellosis (n = 1410), shigellosis (n = 372), Vibrio infection (n = 38), and yersiniosis (n = 587). During this period, 1800 Campylobacter isolates from across BC were speciated by the BC Public Health Microbiology and Reference Laboratory (BCPHMRL): 91·3% were C. jejuni, 5·8% were C. coli and 2·9% were other species (personal communication, 20 Dec. 2013: Dr Linda Hoang, Program Head, Public Health Advanced Bacteriology/Mycology Laboratory, BCPHMRL). There were no known outbreaks of waterborne disease or Campylobacter infection in this region during the study period.
Univariate analysis showed that campylobacteriosis patients were on average older (37·2 vs. 32·6 years) than enteric disease controls, with a greater proportion aged ⩾15 years (Table 1). Further, the proportion and odds of campylobacteriosis were noticeably higher for the ⩾15 years age group compared to enteric disease controls; thus, we dichotomized the age variable into 0–14 and ⩾15 years. The sex distribution was similar, but a greater proportion of campylobacteriosis cases had a lower DI (or higher SES) (40·0% vs. 30·6% of controls). Campylobacteriosis patients were more likely to live in rural settings (8·2% vs. 6·8% of controls) and closer to an ALR (57·3% vs. 53·3% of controls) as well as to live on land parcels serviced by private wells (6·8% vs. 4·8% of controls) (Fig. 1). Campylobacteriosis cases did not seem to live closer to abattoirs than controls. Finally, campylobacteriosis cases showed a peak in June-August whereas control enteric disease cases showed a less pronounced peak in August and a decrease in autumn (Fig. 2). Individual control enteric diseases showed a similar pattern or stable numbers throughout the year.
Spatial distribution of campylobacteriosis cases, enteric disease controls and water system characteristics, Fraser Health Authority, BC, Canada, 2005–009. Cases have been randomly offset by up to 250 m to prevent re-identification by location. The eastern and northern regions of Fraser Health Authority are not shown in the map extent.
Campylobacteriosis cases and control enteric disease cases by month, Fraser Health Authority, BC, Canada, 2005–2009.
Characteristics of campylobacteriosis and control enteric disease cases, Fraser Health Authority, BC, Canada, 2005–2009
ALR, Agricultural Land Reserve.
* Deprivation index data were suppressed from 186 cases because of small population size in certain areas.
One-hundred eighty-six individuals were excluded from further analyses because their SES data were suppressed due to a small population size. Private surface and private mixed water categories were also excluded from further analyses since the number of cases and controls on such systems was small (four individuals in each of the two categories).
Preliminary multivariate analyses showed that proximity to abattoirs did not have a consistent effect on the odds of campylobacteriosis and that excluding these variables improved model fit. Further, distance to the nearest ALR border was excluded from the multivariate analyses due to its collinearity with the urban/rural variable, which better explained the odds of campylobacteriosis.
The final multivariate model showed a significant effect of water index on the odds of campylobacteriosis compared to enteric disease control (P = 0·025), as well as significant interactions for rural area by month (P = 0·004) and age by DI (P = 0·033) (Table 2). Specifically, the odds of campylobacteriosis compared to enteric disease controls was higher for individuals serviced by private wells than municipal surface water systems (OR 1·4, 95% CI 1·1–1·8). In rural settings, the odds of campylobacteriosis were lower in October (P = 0·007) and higher in November (P = 0·014) compared to urban settings (Fig. 3). After peaking in July, the number of campylobacteriosis cases in rural areas dropped considerably in October and increased again for a single month (November), whereas the number of cases of other enteric diseases stayed relatively constant between June and October in rural areas and then dropped to a stable winter level starting in November (data not shown). The odds of campylobacteriosis were higher in individuals aged ⩾15 years, especially in those with a low and medium DI (Fig. 4).
Odds of campylobacteriosis as compared to enteric disease controls by month and by residence in rural or urban settings, Fraser Health Authority, BC, Canada, 2005–2009.
Odds of campylobacteriosis as compared to enteric disease controls by age and socioeconomic status, Fraser Health Authority, BC, Canada, 2005–2009. Low socioeconomic status refers to a high deprivation index and vice versa.
Multivariate regression model
DISCUSSION
We built a model to explain the possible transmission of Campylobacter from the environment to humans in a region of BC, Canada, comparing campylobacteriosis cases to cases of other enteric diseases. We found that private well water, living in a rural environment at certain times of the year and older age in conjunction with higher SES (lower DI) increased the odds of campylobacteriosis.
The main advantage of the case-case design is its ability to limit selection and information biases: controls have similar symptoms, are identified through the same system and are subject to the same recall biases as cases [Reference McCarthy and Giesecke22–Reference Wilson24]. In addition, this design is inexpensive and efficient, especially given that our hypothesis did not require interviewing cases and that we had access to enteric disease case data. A disadvantage of the case-case design is that controls differ from the general population in at least one significant way (their source of illness) and perhaps in general ways (host factors which increase their risk of illness) which limits the generalizability of results [Reference McCarthy and Giesecke22]. As suggested by McCarthy & Giesecke we used a variety of cases as controls to limit the impact that any single disease's exposures would have on the results [Reference McCarthy and Giesecke22]. We included as controls enteric diseases that are travel-associated, locally acquired or both because they all represent the population at risk of campylobacteriosis. Another disadvantage is that the case-case design underestimates the odds of disease when control diseases have similar exposures (e.g. both campylobacteriosis and giardiasis are associated with untreated drinking water) [Reference Gillespie23, Reference Wilson24]. However, 48·0% of controls were cases of salmonellosis, yersiniosis, shiga-toxigenic E. coli infection, Vibrio infection and listeriosis which are infections considered to be mainly foodborne [Reference Scallan25, Reference Thomas26]. The fact that we nonetheless found significant associations supports the hypothesis that transmission of Campylobacter is indeed modulated by the environment, including drinking water.
The odds of Campylobacter infection were higher than enteric disease controls for individuals living on land parcels serviced by a single-user private well as compared to municipal surface water. This has been reported previously for one of the FHA communities and by others [Reference Uhlmann9, Reference Kapperud27–Reference Schonberg-Norio29]. Private wells in our study region and elsewhere are often untreated (personal communication: M. Zubel, FHA, 2013) and can be at risk of contamination from sewage or agricultural runoff [Reference Clark30]. A 2004 survey conducted in our study region found up to 8·4% of private wells contained faecal coliforms [Reference Hii31]. Contamination was localized to a few specific wells associated with shallow construction, disrepair and lack of proper seal and was deemed to be due to local contamination from agricultural runoff. The proportion of wells contaminated was higher in August (7·9%) and October (8·4%) than in February (1·4%). In addition, the Sumas-Abbotsford aquifer which covers 91 km2 (0·6%) of our study region contains high levels of nitrates, which researchers attribute to manure spreading on agricultural land [Reference Hii31–Reference Carey and Cummings33].
We did not find that the higher odds of campylobacteriosis associated with private wells as opposed to enteric disease controls was mediated by proximity to ALR, residence in a rural area, or by season, which would be expected if agricultural runoff is the transmission pathway between the animal reservoir and human infection [Reference Febriani34]. This may be because the enteric disease controls we compared to have similar risk factors, had insufficient power to detect complex interactions (e.g. the role of drinking water may be influenced simultaneously by season and residence in a rural area), the association is explained by other factors (e.g. private wells are contaminated by human sewage rather than agricultural runoff), the presence of a private well is a proxy for another, unmeasured factor, such as animal contact. Although farm animal count data are available, the data are not sufficiently detailed for a valid assessment (e.g. urban areas are included in certain animal dense consolidated census subdivisions leading to misclassification of exposure) [Reference Kuo11].
The odds of Campylobacter infection compared to enteric disease controls differed by rural/urban residence depending on month of the year. In rural areas, the odds of campylobacteriosis were lower in October and higher in November than in urban settings. Many studies have found an increased risk of campylobacteriosis in rural areas and have attributed this to direct animal contact and the impact of animal waste on drinking water and the environment [Reference Spencer35–Reference Arsenault39]. Many other studies have shown an increased risk of campylobacteriosis in the warmer months [Reference Fleury40–Reference Nylen42]. The increased rate and odds of campylobacteriosis in rural areas in November is a novel finding; this was consistent throughout all 5 years of our study and cannot be explained by a known outbreak. When data from the rest of the province were reviewed, this trend was not observed (data not shown).
The reason for this seasonal pattern is unclear. Although we did not have travel status for most cases in this study, previous research has shown that most travel-related enteric disease in BC, including campylobacteriosis and salmonellosis, occurs in the winter and early spring months [Reference Taylor43]. We hypothesize that the increased odds in November may be associated with agricultural runoff in an animal dense rural environment. The start of the rainy season in October and the end of the period when manure spreading is allowed in late September may lead to an increased exposure to Campylobacter via the environment, including drinking water especially in poorly maintained drinking water systems [44, 45] (personal communication: O. Schmidt, BC Ministry of Agriculture, 2013). Febriani et al. also found an association between precipitation in autumn and an increased risk of acute gastrointestinal illness 3 weeks later [Reference Febriani46]. It is not clear why the odds of campylobacteriosis were lower in rural areas in October.
Rural residence was a stronger predictor of campylobacteriosis than proximity to an ALR. Residing in a rural area may represent factors other than proximity that are important, such as direct farm or wild animal contact, potentially contaminated drinking water and food choices. In addition, ALR includes animal growing, produce growing, and forested regions and therefore may be less specific to the hypothesized route of transmission between farm animals and humans [47].
In a preliminary analysis, we included distance to the nearest poultry abattoir because we hypothesized such facilities could increase the risk of campylobacteriosis to nearby residents (through environmental contamination by animal waste) or to abattoir workers [Reference Arsenault39, Reference Jones, Betaieb and Telford48–Reference Wilson50]. However, distance to the nearest abattoir was not a significant factor and did not improve the model fit. In fact, most FHA abattoirs were located in urban settings and probably discharged their waste into the municipal sewage system.
We found that the odds of campylobacteriosis were higher in individuals aged ⩾15 years with medium and high SES (middle and low DI) as opposed to enteric disease controls. In BC, as elsewhere, the incidence of campylobacteriosis is highest in children aged 1–4 years and in adults aged 20–29 years [8]. The risk in children is attributed to behaviours increasing exposure as well as to a lack of immunity; the specific risk factors for young adults remain unclear [Reference Nelson and Harris51]. Many authors have found an increased risk of campylobacteriosis in individuals with higher SES [Reference Spencer35, Reference Bessell36, Reference Arsenault52, Reference Simonsen, Frisch and Ethelberg53]. This could be explained by greater health-seeking behaviour, more frequent restaurant meals, a greater exposure to fresh (as opposed to frozen) chicken meat, travel to endemic countries or other factors. However, a review of travel exposure on a subset of our data (cases occurring in 2008) did not support increased travel in campylobacteriosis cases with higher SES (data not shown).
We could not account for exposures that may have occurred in locations other than the place of residence, but this applied to cases and controls. We were not able to quantify specific exposures such as water ingested and contact with animals. However, other data support our explanatory variables as indicators of exposure. For example, sources show that most (79%) residents in a FHA town drink their tap water (regardless of source) [Reference Teschke54]. Few (15%) residents in three urban and semi-rural FHA towns have farm animal contact (unpublished report: Fraser Health, Campylobacter risk factor and exposure analysis, 2012).
We did not have access to drinking water quality data to validate our water index. However, many authors have found private well water to be at higher risk [Reference Uhlmann9, Reference Kapperud27–Reference Schonberg-Norio29]. We collected the most recently available drinking water distribution data; it is possible that residences connected to municipal water systems after illness onset would have been misclassified. However, we do not believe this would affect cases and controls differently. Most (85%) individuals in our study are served by a single municipal surface water system, whereas only 3–6% of individuals used each of the other water systems. Thus, our estimates of the odds of campylobacteriosis are less precise for non-municipal surface water systems. Nevertheless, we were still able to find statistically significant differences in the odds of campylobacteriosis by water index. For the small water systems (n = 450), we relied on the x,y location coordinates of the drinking water system intake and number of users but lacked a detailed list of land parcels connected to it (as opposed to the parcels connected to the municipal water systems). This may have led to non-differential misclassification of cases and controls.
We conducted a case-case epidemiological study using individually assigned characteristics (e.g. age) as well as aggregate-level characteristics assigned to individuals (e.g. SES). The latter may be misclassified at the individual level, thereby introducing a bias in the estimate. However, this method allowed an efficient assessment of multiple environmental exposures (and possible confounders) associated with a case's place of residence. Some of the findings from this study, such as the increased odds of campylobacteriosis associated with private well water, are probably generalizable to other areas because they are widely supported in the literature. Others, such as the finding of varied risk associated with rural residence in autumn, need to be reproduced and explained further. All findings are relative to other enteric pathogens and may not represent the true population risk.
CONCLUSIONS
In conclusion, we found that the odds of campylobacteriosis were higher than other enteric diseases for those residing on land parcels serviced by private well water, independent of known confounders, but we were not able to show that this was mediated through food animal sources. Nonetheless, the odds of campylobacteriosis were lower in October and higher in November for rural residents, perhaps through exposure to agricultural runoff at the onset of the rainy season. Consideration should be given to further educate and support private well-water users against the risks of enteric diseases. Further work needs to be done to assess the role of agriculture on ground water and the ensuing risk of campylobacteriosis and other enteric diseases. Such studies and real-time surveillance of risk would greatly benefit from a central provincial repository of drinking water data so that province-wide data are accessible for study. In addition, further work is needed to understand the difference in Campylobacter transmission pathways in rural when compared to urban settings.
ACKNOWLEDGEMENTS
The authors acknowledge the financial support provided by the Public Health Agency of Canada. We thank the participating municipalities for providing access to drinking water data and Loretta Bogert O'Brien for obtaining and managing the drinking water data. Thanks also to BC laboratories for the diagnosis and reporting of cases and to the BC Enteric Policy Working Group for its support.
DECLARATION OF INTEREST
None.
