Research Article

18 February 2020

Shifts in the Molecular Epidemiology of Campylobacter jejuni Infections in a Sentinel Region of New Zealand following Implementation of Food Safety Interventions by the Poultry Industry

Authors: Antoine Nohra, Alex Grinberg, Jonathan C. Marshall, Anne C. Midwinter, Julie M. Collins-Emerson, Nigel P. French https://orcid.org/0000-0002-6334-0657Authors Info & Affiliations

DOI: https://doi.org/10.1128/AEM.01753-19

Cite

PDF/EPUB

ABSTRACT

In 2006, New Zealand had the highest notification rate of campylobacteriosis in the world, and poultry was considered the leading source of campylobacteriosis. Implementation of food safety interventions by the poultry industry led to a decrease in the campylobacteriosis notification rate. The aim is to examine the impact of targeted food safety interventions implemented by the New Zealand poultry industry on the source attribution of Campylobacter jejuni infections in a sentinel region. Campylobacter jejuni isolates collected from the Manawatu region of New Zealand between 2005 and 2007 (“before intervention”) and 2008 and 2015 (“after intervention”) from human clinical cases, chicken meat, ruminant feces, environmental water, and wild bird sources were subtyped by multilocus sequence typing. Viable counts of Campylobacter spp. from carcasses were analyzed using a zero-inflated Poisson regression model. In the period before intervention, sequence type 474 (ST-474) was the most common sequence type (ST) recovered from human cases, accounting for 28.2% of the isolates. After intervention, the proportion of human cases positive for ST-474 reduced to 9.3%. Modeling indicated that chicken meat, primarily from one supplier, was the main source of C. jejuni infection in the Manawatu region before intervention. However, after intervention poultry collectively had a similar attribution to ruminants, but more human cases were attributed to ruminants than any single chicken supplier. Viable counts on carcasses were lower in all poultry suppliers after intervention. This study provides evidence of changes in the source attribution of campylobacteriosis following targeted food safety interventions in one sector of the food supply chain.

IMPORTANCE This study provides a unique insight into the effects of food safety interventions implemented in one sector of the food industry on the transmission routes of a major foodborne agent. Following the implementation of food safety interventions by the poultry industry, shifts in the molecular epidemiology of Campylobacter jejuni infections in a sentinel region of New Zealand were observed. Targeted interventions to reduce disease incidence are effective but require continued surveillance and analysis to indicate where further interventions may be beneficial.

INTRODUCTION

Farm-to-fork monitoring of foodborne pathogens helps to improve and maintain the safety of food supply chains. Campylobacteriosis is a significant burden on health care systems due to its high incidence. Campylobacter jejuni is the predominant species isolated from human campylobacteriosis cases in New Zealand and worldwide (1 –3). Most research on campylobacteriosis has focused on C. jejuni not only due to its high prevalence in cases of acute diarrhea but also for its association with serious sequelae, such as Guillain-Barré and Miller Fisher syndromes (4 –6).

Contaminated chicken and red meat, environmental water, milk, and contact with pets and farm animals are the most studied sources of human campylobacteriosis (7, 8). To assess the contribution of different infection sources to the burden of C. jejuni infection, DNA subtyping by multilocus sequence typing (MLST) is widely applied (9). Several studies, worldwide (10 –13), indicated the consumption of chicken meat as the main risk factor for C. jejuni infections. This finding is supported by the high level of contamination with C. jejuni often found in chicken meat (14, 15). However, other infection sources cannot be discounted. For instance, in Finland, bovines and poultry are identified as equally important sources of C. jejuni infections (16).

A case-control study performed in the 1990s in New Zealand indicated the odds of campylobacteriosis were directly and strongly correlated with recent consumption of raw or undercooked chicken and with chicken eaten in restaurants (17). In 2006, New Zealand had the highest rate of campylobacteriosis notification in the world, with >380 cases per 100,000 population (18). Poultry was estimated to be the leading source of campylobacteriosis, associated with 58% to 76% of cases, followed by the ruminant source (20% to 30%) (19). This triggered the New Zealand Food Safety Authority to announce the implementation of a Campylobacter Risk Management Strategy in collaboration with the poultry industry, starting in April 2008 (20). These interventions (21) eventually led to a 60% decrease in the campylobacteriosis notification rate (22). Preliminary data in the Manawatu region of New Zealand, the sentinel site for campylobacteriosis (23), showed that the decline in cases in the first 2 years of the postintervention period was associated with a decline in poultry-associated cases, which was accompanied by a relative increase in the contribution of other sources (24).

The aim of the current study was to compare the source attributions of campylobacteriosis caused by C. jejuni during two time periods, namely, before the implementation of the Campylobacter Risk Management Strategy (“before intervention” [2005 to 2007]) and an extended period “after intervention” (2008 to 2015). The Campylobacter microbial loads on chicken carcasses in the two time periods were also compared.

RESULTS

Description of the samples.

A total of 7,951 samples were collected from humans, chickens, ruminants, environmental water, and wild bird sources between 2005 and 2015 and were analyzed for Campylobacter sp. by culture. A comparison between samples collected before intervention and after intervention is shown in Table 1. The predominant species identified in both periods, among all sources, was C. jejuni, which accounted for >95% of human clinical isolates, followed by Campylobacter coli, accounting for a maximum of 12% in chicken samples.

TABLE 1

TABLE 1 Characteristics of samples collected before intervention and after intervention^a

Source	No. of collected samples	No. of Campylobacter sp.-positive samples	No. of C. jejuni isolates	No. of C. coli isolates	No. of isolates not further processed^b	No. of C. jejuni isolates typed by MLST
Human	774/1,525	670/1,185	652/1,128	17/38	1/19	652/1,040
Chicken	480/963	376/750	331/597	37/90	8/63	330/507
Ruminant	1,058/1,135	347/503	236/319	26/41	85/143	221/244
Environmental water	308/270	140/77	77/40	6/6	57/31	76/29
Wild bird	77/1,361	24/682	21/354	0/1	3/327	21/193
Total	2,697/5,254	1,557/3,197	1,317/2,438	86/176	154/583	1,300/2,013

Before intervention and after intervention are separated by a slash throughout the table.

Campylobacter species culture-positive isolates were not further identified if the isolates were PCR negative for C. jejuni and C. coli.

C. jejuni genetic relatedness.

A total of 3,313 C. jejuni isolates from all sources were typed by MLST (see Tables S1 and S2 in the supplemental material). Before intervention, sequence type 474 (ST-474) was the predominant ST, accounting for 17.2% of the 112 STs, followed by ST-45, ST-50, and ST-48, accounting for 11.9%, 8.0%, and 6.5% of the STs, respectively. After intervention, the proportions of the major STs (ST-45, ST-50, and ST-48) did not change significantly, except for ST-474, which dropped from 17.2% to 6.6%. ST-474 declined among the human clinical isolates from 28.2% (184/652) to 9.3% (97/1,040) and among the chicken meat isolates from 9.7% (32/330) to 5.3% (27/507). There was no strong variation in the frequency of most of the clonal complexes (CCs) between the two periods, except for CC ST-48, which decreased from 27.0% to 14.9%, and CC ST-45, which increased from 15.0% to 18.3%. The decrease in frequency of CC ST-48 and the increase of CC ST-45 were significant (Fisher’s two-tailed exact test, P < 0.05). Tables S1 and S2 in the supplemental material show the frequencies of all the C. jejuni CCs and STs in the different sources.

Figure 1 shows the rarefaction curves of C. jejuni STs from different sources before intervention (Fig. 1a) and after intervention (Fig. 1b). In both periods, human C. jejuni had the greatest number of STs. Rarefaction curves for the environmental water source before intervention and wild bird STs after intervention are steep, but the curves for chicken and ruminant STs appear to approach the asymptote in both periods, indicating a more exhaustive sampling of STs than that in the environmental water and wild bird populations.

FIG 1

C. jejuni population genetic structure.

Minimum spanning trees visualizing the STs before intervention and after intervention are shown in Fig. S1 and S2, respectively, in the supplemental material. Most of the highly abundant STs in both periods occurred in at least three sources (humans, chickens, and ruminants), except for ST-48, which occurred only in human clinical cases and chicken. STs from environmental water isolates (Fig. S1) and wild bird isolates (Fig. S2) generally appeared to be dissociated from the other sources.

Campylobacteriosis source attribution.

Before intervention, the proportional similarity index (PSI) between supplier A and human C. jejuni was significantly higher than that between C. jejuni from other sources and humans (Table 2). Supplier B, supplier “others,” and ruminant sources shared almost similar PSI results, and their PSIs were significantly greater than the PSI between isolates from environmental water or wild bird sources and human isolates. However, after intervention, all isolates from chicken suppliers and ruminant sources had very similar PSIs with human clinical isolates, but these remained significantly greater than the PSI between isolates from human and environmental water or wild bird sources. There was a significant increase, from the before intervention to the after intervention period, in the PSI between isolates from human and supplier B, supplier “others,” and ruminant sources, with a decrease in the PSI between human and supplier A isolates from 0.61 to 0.49.

TABLE 2

TABLE 2 Multilocus sequence type proportional similarity indices between different C. jejuni sources and human clinical isolates before intervention and after intervention

Source	PSI (95% CrI)
Source	Before intervention	After intervention
Supplier A	0.61 (0.51–0.66)	0.49 (0.42–0.52)
Supplier B	0.34 (0.27–0.38)	0.48 (0.40–0.51)
Supplier “others”	0.35 (0.27–0.39)	0.49 (0.39–0.52)
Ruminants	0.37 (0.32–0.41)	0.52 (0.46–0.56)
Environmental water	0.16 (0.07–0.21)	0.18 (0.06–0.19)
Wild birds	0.10 (0.07–0.12)	0.22 (0.17–0.24)

The asymmetric island and Bayesian hierarchical models showed similar attribution estimates (Fig. 2), attributing 59% (95% credible interval [CrI], 52.5% to 67%) and 57% (95% CrI, 47.7% to 67.7%) of human infections, respectively, to supplier A, before intervention. In this period, ruminant sources were the second most common source of infection attributed by the asymmetric island and Bayesian hierarchical models, with attributable values of 24.1% (95% CrI, 18.5% to 29.6%) and 20.2% (95% CrI, 13.3% to 28.4%), respectively. However, both models defined ruminants as the main source for human infection compared with any single chicken supplier after intervention (the asymmetric island model attributed 45.8% [95% CrI, 40.0% to 51.7%] and the Bayesian hierarchical model 40.0% [95% CrI, 32.0% to 46.6%]). Both models showed a significant decrease in the infections attributable to supplier A in the second period.

FIG 2

Altogether, the above models attributed the majority of human campylobacteriosis infections to ruminant and chicken sources, whereas exposure to environmental water and wild bird sources appeared to contribute significantly less to the burden of disease.

Bacterial count results.

A total of 1,126/1,443 (78%) chicken meat samples were Campylobacter sp. positive by culture. The zero-inflated Poisson (ZIP) model (Table 3) showed that all suppliers had a drop in mean log₁₀ counts after the intervention. The estimated proportions of negatives (i.e., probability of zero) before and after intervention were similar among suppliers A and B; the point estimate increased for supplier A and decreased for supplier B, but the confidence intervals overlapped for both suppliers. However, supplier “others” had a higher probability of zero after intervention in addition to lower mean counts. Therefore, considering both the probability of zero counts and the mean counts, all suppliers showed reduced contamination after intervention.

TABLE 3

TABLE 3 ZIP model predictions^a

Supplier	Probability of negative sample (CI^b)	Mean log₁₀/ml (CI)	Data range (CFU/ml)
Before intervention
A	0.14 (0.08–0.23)	2.08 (2.08–2.08)	0–29,300
B	0.43 (0.35–0.54)	2.00 (2.00–2.00)	0–101,000
Others	0.26 (0.15–0.41)	1.52 (1.51–1.53)	0–37,500
After intervention
A	0.20 (0.15–0.25)	1.59 (1.59–1.59)	0–12,880
B	0.33 (0.27–0.38)	1.15 (1.15–1.16)	0–4,500
Others	0.45 (0.40–0.51)	1.16 (1.16–1.17)	0–3,160

The data range is from the raw data, while the other two values (probability of negative sample and mean log₁₀/ml) are from the ZIP model.

CI, confidence interval.

DISCUSSION

We compared the molecular epidemiology and source attribution of human C. jejuni infection in the Manawatu region of New Zealand before and after food safety interventions implemented by the poultry industry. C. jejuni was the dominant species among all sources in both periods, accounting for more than 95% of the clinical isolates, which is consistent with reports from other countries (1, 2, 25, 26). Both periods combined, the estimated percentage of Campylobacter sp.-positive chicken carcasses by culture remained high (∼78%) and was also similar to the percentage observed in other countries (14, 27, 28). However, our model results (Table 3) suggest that the implemented food safety interventions may have led to a significant decrease in the Campylobacter sp. CFU counts on chicken carcasses, which is important to reduce disease incidence.

Our data show a marked drop in the prevalence of the major C. jejuni clonal complex (CC) ST-48, from 27% before intervention, to 13.6% after intervention. This change in CC ST-48 prevalence was attributable to a large decrease in the isolation of the (internationally rare but common in New Zealand) genotype ST-474 among all sources, from 17.2% (224/1,300) to 6.6% (132/2,013). Moreover, ST-474 was a dominant ST isolated from clinical cases before intervention but dropped from 28.2% to 9.3% after intervention and dropped from 9.7% to 5.3% after intervention in all chicken sources, which was also shown in a 2013 study (29). The significant decrease in ST-474 prevalence observed in the sentinel Manawatu region might explain the national 60% decrease in campylobacteriosis notification rate (21). Higher resolution genotyping, such as whole-genome MLST or single nucleotide polymorphism (SNP)-based analyses of ST-474, may help elucidate whether this decline was associated with a change in the population structure of ST-474. Another significant change was the increase in the prevalence of CC ST-45 from 15% to 18.3%. This was associated with a significant increase in the proportion of ST-583 and the detection of 16 different STs belonging to CC ST-45 that were not detected between 2005 and 2007.

The overlap between human- and chicken-associated genotypes is consistent with poultry being important in disease transmission in New Zealand (15, 17, 24, 30), albeit at a lower level after intervention. The PSIs were calculated before intervention and after intervention to evaluate the similarity between ST frequency distributions. Before intervention, the largest similarity was between STs from isolates from humans and supplier A, mainly due to the high prevalence of ST-474, both among clinical cases and supplier A, and also to the low prevalence of ST-474 in samples from nonpoultry sources in New Zealand (31). On the other hand, the similarity between isolates from clinical cases and supplier B, supplier “others,” and ruminant sources increased significantly after intervention and dropped between clinical cases and those from supplier A (the 95% CrIs overlapped among chicken and ruminant sourced isolates). Before intervention, there was a good agreement between the results of the models, where supplier A was the largest significant contributor to human infections (Fig. 2). After intervention, the contribution of the ruminant source increased significantly to become the major contributor to human infections, unlike the supplier A source, which decreased significantly according to the Bayesian hierarchical and asymmetric island models. However, if we combine all the suppliers as one poultry source, the contributions to human infections by poultry and ruminant sources are similar (see Table S3 in the supplemental material), which was not the case before intervention, where chicken contributed between 58% and 76% and ruminants contributed between ∼20% and ∼30% (19).

In summary, our results indicate that >95% of human campylobacteriosis cases are due to C. jejuni in the Manawatu region sentinel site study in New Zealand. The majority of human campylobacteriosis infections could be attributed to ruminant and chicken sources, whereas exposure to environmental water and wild bird sources appeared to contribute significantly less to the burden of disease. Our source attribution comparative analysis after intervention estimated ruminant and chicken sources as equal contributors to the human disease burden, unlike other countries where chickens remain the dominant risk factor, followed by ruminants (32, 33). This apparent change in the source of C. jejuni infections following food safety interventions implemented by the poultry industry provides valuable information to inform decision making by industry and food safety regulators.

MATERIALS AND METHODS

C. jejuni isolates.

(i) Clinical human isolates. We used C. jejuni isolated from human feces that was submitted to clinical microbiology laboratories before intervention (February 2005 to December 2007) and after intervention (January 2008 to December 2015) in the Manawatu region (23). Officially, the implementation of the Campylobacter Risk Management Strategy started in April 2008, but the poultry companies started trial interventions months earlier (24).

(ii) Chicken meat isolates. Over the same period, fresh chicken carcasses from different commercial chicken suppliers were sampled each month from supermarkets in Palmerston North City (suppliers are designated supplier A, supplier B, and five smaller suppliers collectively designated “others”). Culture for Campylobacter sp. was performed as previously described (3).

(iii) Isolates from environmental water, wild birds, and farmed ruminants. (a) Environmental water isolates. C. jejuni bacteria isolated from recreational waterways and pretreatment drinking water samples from the Manawatu area between 2006 and 2014 were used.

(b) Isolates from wild birds and ruminants. A total of 214 C. jejuni isolates from feces of wild birds were obtained between 2005 and 2013. A total of 465 isolates from cattle and sheep feces from 30 different farms obtained between 2005 and 2015 as part of a monitoring program implemented in the area were also used.

Isolates from pigs were not included due to the low prevalence of Campylobacter sp. from this source in New Zealand (we obtained 5 C. jejuni isolates from 650 pig meat samples; data not shown).

C. jejuni isolation, identification, and subtyping.

Methods for Campylobacter sp. isolation from water, feces, and meat samples and for identification and subtyping by MLST have been previously described (3), except that between August 2014 and December 2015, human clinical Campylobacter spp. were identified by multiplex PCR to detect the ceuE gene associated with C. coli and the hipO gene associated with C. jejuni (34).

Campylobacter sp. viable bacterial counts from chicken meat samples.

Enumeration of Campylobacter colonies from the chicken samples was done using manual spread plating and a spiral plater (Wasp; Don Whitley, West Yorkshire, UK). During the rinsing of chicken carcasses, duplicate modified charcoal-cefoperazone-deoxycholate agar (mCCDA) (Fort Richard Laboratories, Auckland, New Zealand) plates were inoculated with 1 ml (manual spread plate) and 50 μl (spiral plater) of rinse and 100-μl (spiral plater) aliquots of resuspended rinse pellet. Plates were incubated microaerobically at 42°C for 48 h and Campylobacter sp.-like colonies were counted using a plate reader (aCOLyte; Synbiosis, England) or manually. The number of CFUs in 1 ml of 200 ml chicken rinse was then extrapolated.

Analysis of C. jejuni genetic relatedness.

The genetic diversities of C. jejuni from different sources were compared by rarefaction analysis using the package vegan in R version 3.1.3 (35).

Minimum spanning trees were used to visualize allelic differences between the MLST of C. jejuni isolates from different sources, using pairwise Hamming distances. Two trees were generated to compare before intervention and after intervention, using the Prim’s algorithm (36), as implemented in the Bionumerics software (Applied Maths).

C. jejuni source attribution.

The relative contribution of various C. jejuni sources to human disease was estimated over the two sampling periods. The similarity between the frequency distribution of human C. jejuni STs and those of the different sources was estimated using proportional similarity indices (PSIs) and their bootstrap credible intervals (CrIs), as previously described (37). The PSI measures the area of intersection between two frequency distributions (38) and ranges between 0 and 1, where 0 indicates no similarity and 1 indicates identical frequency distributions (39). Calculations were performed using R version 3.1.3 (35).

The asymmetric island model was used to probabilistically assign each human isolate to one of four sources (chickens, ruminants, environmental water, and wild birds) (13).

A Bayesian hierarchical model, based on the modified Hald model (40), was used to implement a nonparametric source attribution model to attribute human campylobacteriosis cases to sources in a Bayesian framework with source and sequence type effects (41).

Analysis of Campylobacter colony counts.

A zero-inflated Poisson (ZIP) model was used to assess differences in the Campylobacter viable counts between chicken meat samples collected in the two sampling periods. This model was used because the variance of the viable counts was higher than the mean due to the presence of samples with zero counts, which led to overdispersion. The R Studio version 3.1.3 software (“pscl” package) (35) was used to apply the following equation:

Y_{ij} \sim {\begin{cases} Poisson (V_{ij} λ_{ij}), & Z_{ij} = 0 \\ 0, & Z_{ij} = 1 \end{cases}

Z_{ij} \sim Bernoulli (P_{ij})

where

Y_{ij}

is the total number of colonies across all replicates in isolate i from supplier j,

P_{ij}

is the probability of a zero count,

V_{ij}

is the total volume of rinse that is plated across replicates, and

λ_{ij}

is the average colony count per ml, where

\log (λ_{ij}) = β_{0 j} + β_{1 j} Interventio n_{ij}

logit (P_{ij}) = α_{0 j} + α_{1 j} Interventio n_{ij}

These equations allow both the probability of a zero count (i.e., probability of a negative) and the expected colony count, given that it is positive, to differ between the two sampling periods and between suppliers. The use of total counts and volumes across replicates for each isolate removes potential clustering by isolate, such that each observation in the model may be regarded as independent given supplier and time period. This does not result in a loss of information, as the Poisson model is invariant to aggregation across data with a constant rate.

ACKNOWLEDGMENTS

The sentinel site, as well as microbiological and molecular analyses performed therein, was funded by the New Zealand Ministry for Primary Industries (formerly the New Zealand Food Safety Authority). N.P.F. is supported by the New Zealand Food Safety Science and Research Centre.

We thank Rukhshana Akhter and Lynn Rogers from the mEpiLab team for their laboratory work contribution, Ahmed Fayaz for the extraction of the data, and Phil Carter from the Institute of Environmental Science and Research for help with sequence typing.

Supplemental Material

File (aem.01753-19-s0001.pdf)

Download
622.91 KB

REFERENCES

Horrocks SM, Anderson RC, Nisbet DJ, Ricke SC. 2009. Incidence and ecology of Campylobacter jejuni and coli in animals. Anaerobe 15:18–25.

Request Username

Create a new account

Change Password

Your password must have 8 characters or more and contain 3 of the following:

Password Changed Successfully

Verify Phone

Congrats!

ABSTRACT

INTRODUCTION

RESULTS

Description of the samples.

C. jejuni genetic relatedness.

C. jejuni population genetic structure.

Campylobacteriosis source attribution.

Bacterial count results.

DISCUSSION

MATERIALS AND METHODS

C. jejuni isolates.

C. jejuni isolation, identification, and subtyping.

Campylobacter sp. viable bacterial counts from chicken meat samples.

Analysis of C. jejuni genetic relatedness.

C. jejuni source attribution.

Analysis of Campylobacter colony counts.

ACKNOWLEDGMENTS

Supplemental Material

REFERENCES

Information

Published In

Copyright

History

Permissions

Keywords

Contributors

Authors

Editor

Notes

Metrics

Citations

Figures

Media

Share

Share the article link

Share with email

Share on social media