Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2004, p. 6501-6511, Vol. 70, No. 11
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.11.6501-6511.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Mathematics and Statistics, Lancaster University, Lancaster,1 DEFRA Epidemiology Fellowship Unit, Department of Veterinary Clinical Science, University of Liverpool, Liverpool, United Kingdom,3 Datalogisk Institut, Aarhus Universitet, Aarhus, Denmark2
Received 10 December 2003/ Accepted 31 May 2004
|
|
|---|
|
|
|---|
90%) are attributed to infections with Campylobacter jejuni (19), although Campylobacter coli is increasingly recognized as an important pathogen and is estimated to account for 25,000 cases of infectious intestinal disease, at a cost of four million pounds, by the UK National Health Service (41). Although most infections are self-limiting, some are associated with chronic, debilitating sequelae such as Guillain-Barré syndrome and Miller Fischer syndrome (30). Most infections are believed to result from the ingestion of contaminated food, although the role of other, nonfood exposures in the epidemiology of sporadic campylobacteriosis is still unknown. C. jejuni has been isolated from a range of food sources, including poultry (33), red meat, and milk (12, 18). Point source outbreaks are thought to be relatively uncommon compared to those by other major enteric pathogens, although there is increasing evidence for localized transmission (6). The primary source of contamination is believed to be animal feces. This is consistent with high carriage rates in poultry, pigs, and cattle (23) and with molecular evidence showing similar genotypes in farm animals and humans (9, 13). Contamination of the environment by domestic and wild animal feces (and that of humans) presents an alternative exposure pathway for human infection, for example, via drinking (10, 16, 37) and recreational water use (1). Humans may also be exposed to voided animal feces in the environment through outdoor activities such as camping, walking, and picnicking. The work described here forms part of a larger study of the epidemiology and genetic diversity of thermophilic Campylobacter spp. isolated by structured sampling of an area of the United Kingdom that is intensely farmed and also widely used for recreation. This paper describes the sampling protocol and the distribution of campylobacters in the environment to the level of the species, with particular emphasis on C. jejuni isolated from cattle feces. We examined the spatial distribution of C. jejuni, C. coli, Campylobacter lari, and Campylobacter hyointestinalis in cattle feces in the environment and tested the null hypothesis that Campylobacter spp. are distributed randomly, with no evidence of spatial clustering. We used model-based spatial statistical methods to describe the nature and extent of spatial and nonspatial clustering and to explore the determinants of the observed pattern of C. jejuni in cattle feces in the environment. In particular, we looked for evidence of spatial effects beyond the within-field effect that could be attributed to larger-scale spatial processes. Further species typing based on restriction fragment length polymorphism analysis (flaA and pulsed-field gel electrophoresis) and multilocus sequence typing will be reported elsewhere (25).
|
|
|---|
Figure 1a shows, for a subregion of the study area, the locations from which fecal material and other samples were collected. For the most part, the sampling locations were situated in clusters of four points forming the corners of squares with 50-m-long sides. These squares are referred to as "primary squares," the centers of which form a 15- by 15-km grid and are spaced 667 m apart. Each of the secondary points in the primary squares was in turn surrounded by four tertiary points (not shown) that were each 25 m apart. Soil samples were taken from the tertiary points, and fecal pats were counted and sampled in circles with a 5-m radius centered on the tertiary points. In three areas, a primary square was surrounded by 15 more primary squares, thus creating "fill-in" areas of 400 by 400 m. A small number of secondary points lay in areas which were not suitable for sampling (e.g., areas covered by buildings or roads) or on land for which sampling permission had not been given.
![]() View larger version (10K): [in a new window] |
FIG. 1. Locations of environmental samples for isolation of Campylobacter spp. in a subregion of the study area. (a) Secondary points where feces were sought. (b) Secondary points where cattle feces were present.
|
The sampling points were overlaid onto a 1:2,500-scale digitized map of the study area and were identified in the field by use of a combination of geographical features (e.g., field boundaries and buildings), a compass, and a tape measure. All sampling was carried out between 22 May and 26 July 2000. Figure 1b shows the secondary points in the subregion where cattle feces were present. Notice that no feces were collected from many primary squares and that other primary squares contained fewer than four secondary points where cattle feces were present. Only two of the fill-in squares in the study area contained cattle feces.
Sample collection.
Up to two cattle fecal samples were taken from the vicinity of each tertiary point: the nearest sample, defined as the sample nearest to the sampling point, and the freshest sample, defined as the sample which, on visual inspection, was considered to be the most recently voided. The nearest samples provided an indication of a typical pat in the sample region, whereas the freshest samples were expected to be the most reliable indicator of the presence of Campylobacter spp. in the area. When there was only one fecal pat, it was only sampled once and was considered the nearest sample, with the freshest pat being considered missing. Samples from the tertiary points were merged into two pooled samples for each secondary point, with one pool of the nearest samples and another of the freshest samples. Pooling was performed because the number of pats from tertiary points was too large for each one to be tested.
Soil samples were taken from the surface down to a depth of approximately 2 cm. If the sampling point was covered with feces, soil was taken from an uncontaminated location adjacent to the sampling point. Wildlife and sheep feces were identified by scanning the entire area of the primary square and were pooled in the field at the level of the primary square for wildlife samples and the secondary point for sheep samples. Wild mammal feces (mainly from rabbits and badgers) were easily identified by their size, shape, and in the case of badger samples, location in well-marked pits (latrines). Wild bird fecal samples were identified by their color, consistency, and location. Water samples were identified when we scanned the area of the primary square, but the samples were not pooled.
Isolation and species identification.
One milliliter of homogenized pooled fecal sample was added to 9 ml of campylobacter enrichment broth (Lab M, Bury, United Kingdom). The broth was incubated at 42°C for 24 h under microaerophilic conditions in a variable atmosphere incubator, inoculated onto campylobacter blood-free agar (Lab M) containing a CA antibiotic supplement (X112), and incubated as described above for a further 48 h. Three suspect colonies, representing the variability in colony morphology, were subcultured onto Columbia agar (Lab M) supplemented with 5% defibrinated horse blood (Southern Group Labs, Corby, United Kingdom) and incubated as described above for 24 to 48 h. Each isolate was presumptively identified by standard methods (no growth in O, presence of catalase and oxidase, Gram stain, cell size, and morphology) and then frozen in Microbank tubes (Pro-Lab Diagnostics, Neston, United Kingdom) at 70°C.
Presumptive positive isolates were identified as being either C. jejuni, C. coli, C. lari, C. hyointestinalis, or genus-specific Campylobacter by a series of single-reaction PCRs. Isolates were inoculated onto Columbia blood agar and grown microaerophilically for 48 h at 37°C. One colony was arbitrarily selected and suspended in 500 µl of sterile distilled water. The DNA was extracted by boiling the suspension for 20 min. One microliter of the boiled cell suspension was used in each 25-µl reaction volume.
The primers for C. jejuni were JEJ1 (5' CCTGCTACGGTGAAAGTTTTGC 3') and JEJ2 (5' GATCTTTTTGTTTTGTGCTGC 3'), resulting in an amplicon of 793 bp, and the primers for C. coli were COL1 (5' ATGAAAAAATATTTAGTTTTTGCA 3') and COL2 (5' ATTTTATTATTTGTAGCAGCG 3'), resulting in an amplicon of 894 bp (20), and the target for both was the putative virulence marker ceuE. The reaction mixtures for C. jejuni and C. coli contained the following: 0.2 mM (each) dATP, dCTP, dGTP, and dTTP, a 1 µM concentration of each primer, 20 mM (NH4)2SO4, 75 mM Tris-HCl (pH 8.8), 0.01% (wt/vol) Tween 20, 3.5 mM MgCl2, and 0.5 U of Taq DNA polymerase (Abgene). The reaction mixtures were initially held at 94°C for 5 min. The amplification cycle was performed 30 times and consisted of denaturing at 94°C for 30 s, annealing at 57°C for 30 s, and extension at 72°C for 1 min. At the end of cycling, the samples were held at 72°C for 5 min. The primers for C. lari were CL594F (5' CAAGTCTCTTGTGAAATCCAAC 3') and CL1155R (5'ATTTAGAGTGCTCACCCGAAG 3'), resulting in an amplicon of 561 bp (26). The primers for C. hyointestinalis were CFCH57F (5' GCAAGTCGAACGGAGTATTA 3') and CH1344R (5' GCGATTCCGGCTTCATGCTC 3'), resulting in an amplicon of 1,287 bp (26). The primers for genus-specific Campylobacter PCRs were C412F (5' GGATGACACTTTTCGGAGC 3') and C1288R (5' CATTGTAGCACGTGTGTC 3'), resulting in an amplicon of 816 bp (26). The reaction conditions for C. lari, C. hyointestinalis, and genus-specific Campylobacter were the same as those for C. jejuni and C. coli except that 2.5 mM MgCl2 was used. The target for genus-specific Campylobacter, C. lari, and C. hyointestinalis PCRs was the 16S ribosomal DNA region of the genome. Samples were initially held at 94°C for 4 min. The amplification consisted of 25 cycles of denaturation at 94°C for 1 min, annealing for 1 min at 64°C for C. lari, 65°C for C. hyointestinalis, and 55°C for genus-specific Campylobacter, and extension at 72°C for 1 min. At the end of cycling, the samples were held at 72°C for 5 min. The amplified DNAs were analyzed by electrophoresis in 1.5% agarose gels run at 120 V for 75 min, and the gels were stained with ethidium bromide.
Data recorded.
The age of each pat was scored by use of a four-point scale as follows: 1, freshly voided with no signs of aging (animal may be observed defecating, no crust); 2, recently voided with some early signs of aging (thin skin-like crust on surface); 3, clear signs of aging but still retaining integrity (thicker rigid crust, mold on surface, obvious insect activity); 4, aged (degenerated, fibrous, and dehydrated, with grass growing into fecal pat). The scoring system was designed and piloted with all observers to ensure consistency of scoring. The degree of consistency was continuously monitored throughout the period of study. Although the scoring is categorical and not continuous, intermediate scores, such as 2.5, were also permitted. A small number of missing values for pat age corresponded to pats for which the age was not recorded.
The weather conditions at the time of sampling were observed as a means to adjust for the potential effect of UV light on the isolation of Campylobacter spp. and the effect of rain on age scoring of cattle feces. Weather conditions were given simple scores (scores for rain, none, light, heavy; scores for sun, cloudy, overcast, sunny). The principal habitat at each secondary point was also known. The number of pats within 5 m of each tertiary point was recorded, and the average of these counts within a secondary point was considered to be a potential influence on the presence of C. jejuni in a pooled sample.
Statistical methods.
As discussed in Results, the spatial structure of each species was examined by the production of semivariograms (7). To construct a variogram, we calculated the spatial separations (Uij = ||Si Sj||) and the squared differences [Vij = (Yi Yj)2/2] for each pair of observations, Yi and Yj, located at points Si and Sj, respectively. Here Si refers to the location of the secondary point from around which the fecal pats comprising the pooled sample were taken. The collection of squared differences (Vij) and distances (Uij) is the variogram cloud, consisting of a large number of points which must be smoothed or averaged before they yield useful information. Typically, the spatial separations (Uij) are grouped into n equally spaced bins, with the bin width being termed the step size and with the variogram consisting of the bin centers (Uk, where k = 1...n) and the averages (
k) of all variogram ordinates (Vij) whose spatial separation (Uij) is in bin k.
Because of the regular distribution of sampling points in this study, the bins were given various sizes, which allowed for the fact that there were many points separated by distances of multiples of 666 m (the distance between primary squares) and very few with a separation of, say, 300 m. The first bin contained all pairs with a spatial separation of zero, i.e., the freshest and nearest samples from each secondary point. The second and third bins corresponded to distances within a primary square, namely, a 50-m separation for horizontal and vertical distances and a 71-m separation for diagonal distances. The seven subsequent bins corresponded to distances within fill-in squares, which were larger than the 50 m within primary squares but smaller than the 666 m between primary squares. The final bins related to pairs in neighboring primary squares and pairs whose primary squares were second neighbors. The variograms computed were omnidirectional variograms, meaning that only the separation distance and not the angle of pairs was taken into account.
In order to obtain a measure of uncertainty for variograms computed from the data, we compared them to variograms obtained from permuted data, for which the spatial locations were randomly permuted 99 times. Permuted data should not exhibit spatial dependence, and a species can be judged to have spatial dependence if its variogram shows more spatial structure than any of the permuted data sets. The analysis was done with the R package geoR (36).
A statistical model for the prevalence of C. jejuni was used to judge the contribution of various factors to the presence of C. jejuni. The relatively small numbers of positive isolates for other Campylobacter species and sample types precluded a similar analysis of these data. In this model, the subscripts i, j, and k are used to denote the ith primary square, the jth secondary point, and the kth tertiary point. We use
to refer to the type of sample, either the freshest or nearest pooled sample. Yij
refers to the presence or absence of C. jejuni in the pooled sample made from the pats near the tertiary points ijk of type
surrounding secondary point ij, taking the values 1 for the presence and 0 for the absence of C. jejuni.
In the model, each pooled sample Yij
is assumed to test positive for C. jejuni with the probability pij
, according to the formula Yij
Bernoulli(pij
).
The probability pij
varied from site to site and depended on the covariate Xij
. The environmental variables listed in Table 3 were considered possible covariates. Two weather-related covariates were the presence or absence of rain and a categorical variable denoting the amount of sun on the day of sampling. Other covariates were the predominant habitat within the square, the number of cattle pats in the sampling area, the presence or absence of bird or rabbit feces in the square, and whether the bird and rabbit feces present were positive for C. jejuni.
|
View this table: [in a new window] |
TABLE 3. Percentages of C. jejuni-positive cattle fecal samples stratified by various environment-related categories
|
was denoted Aijk
, and the age score used to describe each pooled sample was Ãij
= mink(Aijk
), or the age of the youngest pat of type
at point ij. The values assigned to the age categories were arbitrary, and there was no reason to expect that a pat of age 4 was twice as old as a pat of age 2. The average age score of the pats in the pooled sample was therefore unlikely to be meaningful. Furthermore, fitting the minimum age as a categorical variable, in which each age category has its own parameter, is likely to fit the data better than fitting age as a continuous variable, with the age effect being Ãij
ß.
Since pij
must be between 0 and 1, we modeled the inverse logit transform of the probabilities by using the following equation:
![]() | (1) |
A random effects model was used to ascertain the extent of spatial dependence in the data. The question of interest was as follows: after the effects of pat age and other covariates have been taken into account, are cattle fecal pats that are closer to each other more likely to have similar degrees of C. jejuni prevalence than pats which are further apart? Since the primary squares formed a regular lattice, a convenient model for these data was a Markov random field model (7).
In this spatial model, we allowed pij
to depend on the fixed effects, a latent spatial process (Ai), and a spatially uncorrelated random effect (Bi). Thus, all of the pats within a primary square had the same probability of being positive for C. jejuni (subject to the covariate Xij
values being equal), and neighboring primary squares were potentially more similar than primary squares that were further apart. The full model is Sij
= Xij
ß + Ai + Bi. Ai follows a Markov random field with the following formula:
![]() | (2) |
N(0,
B2). Here Ai = {Aj; j
i},
i refers to all of the primary squares which are neighbors of i, with the neighborhood being the four squares to the immediate north, south, east, and west, and ni is the number of points for this neighborhood in the study region. The fill-in areas were deemed to correspond to one primary square located at the center of the group to preserve the lattice structure.
The nature of the spatial transmission of C. jejuni can be inferred by examining the marginal variances of Ai and Bi in equation 2. Small values for the variances of both Ai and Bi would indicate no spatial dependence in the data. A large variance of Ai would be expected if there was a spatial correlation at distances of 600 m or more. If there was a spatial structure at short scales but not at 600 m or more,
B2 would be large, but var(Ai) would be negligible.
Markov random field models of this type are usually implemented within a Bayesian framework and fitted by use of the Gibbs sampler (17). For this study, we used BayesX (24) software, which has by default gamma(1, 0.005) distributions as priors on variance parameters, as recommended by Besag et al. (3), and uniform priors on the coefficients of the fixed effects. These priors are used throughout this paper. The chain was thinned by taking only every 300th sample and was run until 6,000 samples were obtained after a burn-in. Autocorrelation plots and time traces were used to monitor chain convergence.
Field-level variation.
Anticipating the results shown below, we explored random effects at the field level rather than the level of primary squares. One possible reason for the strong square effect shown below could be that the spatial correlation decays to zero somewhere between 50 and 600 m. Another possibility is that pats from the same primary square were likely to be located in the same field and therefore to come from the same herd or even from the same animal. For this study, detailed information on the herds and animals that had grazed each pasture was not available, although the field that each sample belonged to could be inferred from the spatial locations and habitat types recorded.
Recall that a sample tested for C. jejuni was derived from between one and four pats located at tertiary points. A sample was deemed to belong to a given field if all of the pats comprising the sample were located in that field. Samples consisting of pats from more than one field were not included in the analysis.
Using Yij to refer to pooled sample j in field i, we used the following model to incorporate field effects:
![]() | (3) |
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Percentages of samples that tested positive for Campylobacter spp., listed by species and sourcea
|
![]() View larger version (7K): [in a new window] |
FIG. 2. Locations of secondary sampling points that were negative (·) and positive ( ) for C. jejuni in cattle feces. Samples were identified as being positive by standard culture methods, and isolates were confirmed by PCR.
|
![]() View larger version (6K): [in a new window] |
FIG. 3. Centers of primary squares from which at least one sample of cattle feces tested positive ( ) and from which no samples tested positive (·) for C. lari and C. hyointestinalis. (a) C. lari. (b) C. hyointestinalis.
|
![]() View larger version (19K): [in a new window] |
FIG. 4. Variograms with observations at the same secondary point ( ) and at distances within primary squares (), within fill-in squares ( ), and between primary squares (+), along with permutation envelopes. (a) C. jejuni. (b) C. coli. (c) C. lari. (d) C. hyointestinalis. (e) Other Campylobacter spp.
|
|
View this table: [in a new window] |
TABLE 2. Percentages of positive samples for each Campylobacter spp. according to the age of the youngest pat comprising the pooled samplea
|
Spatial exploratory analysis.
Figure 4 shows variograms for the Campylobacter sp. data along with permutation envelopes. If the dissemination of a Campylobacter sp. were spatial, we would expect an upward trend on the variogram, with variogram ordinates to the left of the plot lying below the gray confidence region. If spatial effects were not important, the variogram ordinates should be flat and lie within the confidence region.
The various widths of the confidence region were due to the different numbers of points for different spatial separations. A primary square with two samples at each of the four secondary points yielded four pairs at coincident points and eight pairs each at separations of 50 and 71 m, causing the confidence regions to all be slightly wider at the first point than at the subsequent two points. The permutation envelopes were widest at distances between 75 and 600 m due to the fact that there were only small numbers of pairs at these separations, which are in the fill-in squares.
For the most part, the variograms lay inside the permutation envelopes, meaning that there is little evidence of spatial dependence. C. jejuni seems to have less variance within primary squares than do the permutations, which suggests that the data exhibit dependence at short ranges.
Modeling the presence of C. jejuni in cattle feces.
A statistical model was used to test the influence of several possible covariates, specifically the pat age, the presence of rain and sun, the presence of bird or rabbit feces and whether they were contaminated with C. jejuni, the number of pats in the area, and the habitat type of the sampling region. The first covariate considered was the ages of the fecal pats in order to test the hypothesis that older pats were less likely to test positive than fresher fecal material. As expected, the minimum age score proved a better predictor than the mean age score, and treating age scores as categorical was an improvement over fitting the age scores as continuous. The estimates of the age effect for young pats were broadly similar and had large standard errors. Combining age categories 1, 1.5, and 2 into a single category had a negligible effect on the deviance and reduced the number of parameters in the model. The presence of C. jejuni in bird feces and the presence of rain were significantly related to the isolation of C. jejuni, whereas the other variables were not statistically significant.
Table 4 shows parameter estimates and standard errors for the model incorporating the significant covariates. P values refer to the test that the true parameter value is zero, which along with the standard errors were based on a normal approximation of the model.
|
View this table: [in a new window] |
TABLE 4. Maximum likelihood estimates and standard errors for factors related to the isolation of C. jejuni from cattle feces
|
Modeling of spatial effects.
The nature of the spatial variation of C. jejuni could be inferred by examining the estimated variances of the independent square effect Bi and the spatially correlated square effect Ai. The posterior mean of the independent effect
B2 was 1.29, with a posterior standard deviation of 0.38. The marginal variance of the spatially correlated term Ai, calculated as the sample variance of the simulated Ai at each iteration, was an order of magnitude smaller, with a posterior mean of 0.040 and a posterior standard deviation of 0.084. The posterior means of the coefficients on the fixed effects were all within half of a standard error of those obtained above, and the posterior standard deviations were similar to the standard errors for the data described above.
The estimated variance of the spatial process suggests that, although feces sampled within the same primary square tended to have a similar likelihood of being positive for C. jejuni (as confirmed by the large
B2), neighboring primary squares were mostly independent of one another. Figure 5 shows the posterior means of Ai and Bi, with the shading of the squares indicating the values of the primary square effects. The irregularly shaped regions in Fig. 5 are explained below. Notice the different scales for the two figure panels, with the Bi values spanning a larger range than the Ai values. This larger variability in Bi is more evidence that the independent random effect term Bi has more influence on C. jejuni prevalence than the spatially correlated term Ai.
![]() View larger version (50K): [in a new window] |
FIG. 5. Posterior means of the spatial random effects, square random effects, and field random effects. White areas correspond to locations without observations. (a) Spatial effect Ai. (b) Independent random effect Bi.
|
The relative importance of the two effects can be affected by the priors used for the variance parameters. As a diagnostic tool, we fitted the model by using only the spatial Ai process, in effect forcing the algorithm to find a spatial pattern. The result of fitting this model was that the estimated spatial effects were similar to the Bi values fitted before. Thus, the conditionally independent square effect Bi proves to be the more important effect, persisting even when the model forces the latent process to be spatially correlated.
The two fill-in squares provided information about how the isolation of C. jejuni from cattle feces varies over short distances. Since the secondary points within the two squares formed a lattice, a Markov random field model could be fitted for these regions, with the spatial effect Aj and the independent effect Bj now representing the risk of C. jejuni presence at a secondary point j in a given fill-in square. The posterior means of the Aj values, with the mean corrected, were all <0.02 in magnitude, and the 90% credible intervals all contained zero. While the two fill-in squares had different risks for C. jejuni, there was no evidence that within a fill-in square neighboring points were more similar than squares that were further apart.
Modeling of field effects.
Table 5 shows the parameter estimates and credible intervals for the coefficient ß and the vari-ance of the field effect Ai in equation 3. These parameters are similar to those found above, as for the most part the fields and primary square groups were coincident.
|
View this table: [in a new window] |
TABLE 5. Parameter estimates for a model in which C. jejuni prevalence depends on the covariates listed and a random effect at the field level
|
Comparing the fit of the field model to the square model was difficult because the models are not nested. The variance of the field effect was only slightly higher than the square effect variance. The models could be fitted by using penalized pseudo-likelihood and then could be compared by examining their approximate AICs, although this showed an AIC that was only marginally in favor of the field model. We conclude here that the field effect is more appropriate for these data because of the more satisfactory interpretation of boundaries based on fields rather than primary squares.
A further question of interest was whether fields that are close together tend to have risk factors which are more similar than fields which are further apart. Figure 6 shows variograms of the posterior means of the field effects Ai, in which the location of each field is defined as the average of the x-axis and y-axis coordinates of the boundary vertices. The gray region is a 95% permutation-based confidence region, which is wide at short distances due to the small number of fields at a spatial separation of <600 m. The variogram ordinates at short ranges are close to the edge of the confidence region, and although the variogram does not indicate a statistically significant amount of spatial correlation, the issue warrants further investigation.
![]() View larger version (32K): [in a new window] |
FIG. 6. Variogram of the posterior means of the field risk factors for C. jejuni ( ) and of the permutation-based confidence region (shaded area).
|
The model described above was applied to the data for C. hyointestinalis, C. lari, and C. coli. Positive results for C. lari and C. coli were very sparse, which led to numerical problems when fitting the model both in the Bayesian framework and with pseudo-likelihood. The prevalence of these two species was therefore too low to permit the modeling of field effects, and we restricted our attention to C. hyointestinalis and C. jejuni.
Due the sparsity of C. hyointestinalis-positive results, the parameter estimates for C. hyointestinalis are less precise than those for C. jejuni. Figure 7 shows the field effects for C. jejuni plotted against the field effects for C. hyointestinalis. If the two species were likely to occur together, we would expect fields at a high risk of C. hyointestinalis to also be at a high risk for C. jejuni. Figure 7 shows no such relationship, suggesting that the species are independent of one another.
![]() View larger version (16K): [in a new window] |
FIG. 7. Posterior means of field risk factors for C. jejuni plotted against the field risk factors for C. hyointestinalis.
|
|
|
|---|
For this study, thermophilic Campylobacter spp. were commonly isolated from a wide range of samples. In addition to producing food for both local and national consumption, the intensively farmed area contains two small towns (with a population of approximately 2,000 in each) and is visited by a larger population for outdoor leisure activities. The presence of C. jejuni in environmental samples of domestic and wild animal feces and water found in this study suggests that humans are likely to be exposed to this microorganism through recreational and occupational activities in addition to food pathways. We have demonstrated a moderately high prevalence of C. jejuni in environmental cattle feces throughout the study area. The spatial pattern is consistent with that of a ubiquitous organism, with a short-scale variation in infection intensity that cannot be explained by variations in the ages of the fecal material.
The high prevalence of C. jejuni in fresh cattle feces is consistent with findings from other studies (21, 39). The role of ruminants as reservoirs of Campylobacter spp. and their potential importance as sources of human infection were reviewed by Stanley and Jones (38). In addition to direct contact through occupational and recreational exposure, C. jejuni-positive cattle feces may contaminate food products (milk and meat), public and private water supplies (10, 37), and recreational waters (1). There are more than 50 dairy farms in the study area that produce milk products for local and national consumption, including an "open" farm on which visitors may have direct contact with livestock. Milk and other dairy products have been associated with Campylobacter outbreaks (2, 16, 18), and there is evidence of an association between particular genotypes of C. jejuni and cattle (27).
In contrast to food pathways, the role of environmental contamination with cattle feces in human disease is less clear. It is possible that direct and indirect exposures to cattle feces through these multiple, nonfood pathways account for many of the apparently sporadic cases of human Campylobacter infection. There was a large amount of C. jejuni-positive fecal material in the environment of our study area. The recovery of C. jejuni from cattle feces was highest for freshly voided feces and declined with the apparent age of the fecal pat, but C. jejuni was still isolated from highly degraded fecal material. The C. jejuni prevalence was lower in samples collected during rainy periods, although this was most likely due to the effect of moisture on the age scoring of fecal material, but there may also have been a dilution effect due to the higher water content of samples collected during periods of wet weather. The decline in recovery with age was likely due to the death of the organism over time but may also have been due to the organism entering a viable but nonculturable state (29). Thermophilic Campylobacter spp. have been shown to survive well in dairy cattle slurry (40). Surface runoff and slurry spreading are therefore likely to disseminate Campylobacter spp., thereby increasing the potential for human exposure.
C. lari was isolated mainly from water and from cattle and wild bird feces in the study area. This species was identified as an opportunistic pathogen of an immunocompromised patient in 1984 (31) and as a human enteric pathogen in 1985 (42), although the origin of human infections is unclear. C. lari occurs at a much lower prevalence than C. jejuni and C. coli, being responsible for 0.1% of reported cases in the United Kingdom in 2000 (5). C. lari has been isolated from a range of environmental sources in other studies and has been associated with birds (43), shellfish (11), and water environments (22).
After we adjusted the data for the age of the fecal material, cattle feces were more likely to be positive for C. jejuni if they were in an area in which wild bird feces were present and were positive for C. jejuni. There are several possible explanations for this finding. Firstly, there may be a direct transmission between cattle and wild birds. Wild birds may be a source of transmission to cattle, with C. jejuni-positive wild bird feces presenting a fresh pool of infected fecal material for direct transmission, or alternatively, wild birds may be infected by the ingestion of infected cattle feces. This explanation was supported by a preliminary analysis of the genotypes of C. jejuni found in wild birds and cattle (flaA, pulsed-field gel electrophoresis, and multilocus sequence typing) that showed that many indistinguishable genotypes were present in cattle and wild bird feces in the study area (15). However, this finding may be attributable to factors unrelated to direct transmission that are associated with the isolation of C. jejuni-positive cattle and wild bird feces (e.g., another common source of infection or other common risk factors for carriage and shedding). Other studies have shown that wild birds are carriers of C. jejuni (34, 43), including strains that have been isolated from human infections (4).
Our analysis found no evidence of spatial dependence for C. jejuni in cattle feces at distances of >600 m. This observed pattern was therefore not consistent with a large-scale transmission attributable to watercourses, wildlife territories, or other geographical features that transcended field and farm boundaries. Similarly, although the data were sparse, there was no evidence for spatial dependence for any of the other Campylobacter spp. isolated, and the presence of one species of Campylobacter was not associated with the presence of another. After we allowed for the ages of the feces in the environment, samples taken from the same field proved to be correlated, most likely because these pats would have come from the same groups of animals within the same herd. In this study, we focused entirely on feces found in the environment, with the aim of systematically sampling the entire area with minimal disruption to the farmers and landowners. Therefore, the type and identity of the animals that had grazed the fields were not always available. Groups of cattle grazing the same field were likely to belong to the same age and management group, and these factors have been shown to be associated with the prevalence and level of fecal carriage (32, 39). However, since the spatial resolution was not sufficiently fine and since we could not obtain reliable information on the identities of farms using each field, we were unable to examine whether different fields within the same farm tended to be similar and whether neighboring fields from different farms were correlated. A further study of short-range spatial correlation is therefore warranted, with a focus on the effects of fields and herds. A Markov random field model could then be used to assess the transmission of C. jejuni between neighboring fields and herds.
Assessing the dependence of different species was difficult because of the scarcity of positive results for species other than C. jejuni. This may be due to the dominant species masking the presence of other, less common species. This phenomenon would weaken the dependence between the species, as fields with a high prevalence of C. jejuni would hide the C. hyointestinalis even if C. hyointestinalis was also present in large numbers. Since the intensities of the other species were low, many more samples per field would have to be tested and more colonies would need to be examined for inferences about field-level effects to be possible.
This work was funded by DEFRA UK.
|
|
|---|
lang/bayesx/bayesx.html.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»