Utilizing a combination of molecular and spatial tools to assess the effect of a public health intervention

Prevalence and Genotyping of Campylobacter jejuni and Campylobacter coli from Ovine Carcasses in New Zealand

ABSTRACT

A pilot survey was performed to determine the prevalence of Campylobacter jejuni and Campylobacter coli on three age classes (lamb, hogget, and mutton) of ovine carcass trim postdressing and prechill. Sampling of hogget carcasses was undertaken 6 months before sampling of lamb and mutton carcasses. A total of 120 trim samples were collected from 11 processing plants across New Zealand. All samples were enriched and screened using PCR for the presence of C. jejuni and C. coli, and isolation was attempted for all screen-positive samples. Enumeration of Campylobacter from lamb trim samples showed that Campylobacter bacteria were present in very low numbers (<10 CFU/g). The overall prevalence of Campylobacter for ovine trim based on PCR detection was 33% (39 of 120 samples), with prevalences for hogget, lamb, and mutton carcass trim of 56% (28 of 50), 11% (4 of 35), and 20% (7 of 35), respectively. Whole genome sequencing was performed on a selection of C. jejuni and C. coli isolates, and the data were used to subtype using multilocus sequence typing (MLST) and whole genome MLST. Twenty-five MLST sequence types (STs) were identified among 44 isolates, including ST42, ST50, ST3222, and ST3072, which have been previously reported to be associated with ruminant sources. Four novel STs were also identified. Whole genome MLST analysis further discriminated isolates within a single ST type and demonstrated a genetic diversity among the ovine isolates collected. Genes associated with the oxacillinase class of β-lactamase enzymes were identified in 41 of 44 Campylobacter isolates. This study provides preliminary data that can be incorporated into existing source attribution models to assist in determining the potential contribution of ovine sources to the burden of campylobacteriosis in New Zealand.

HIGHLIGHTS

Comparison of the Pathogenic Potential of Campylobacter jejuni, C. upsaliensis and C. helveticus and Limitations of Using Larvae of Galleria mellonella as an Infection Model

Campylobacter enteritis in humans is primarily associated with C. jejuni/coli infection. Other species cause campylobacteriosis relatively infrequently; while this could be attributed to bias in diagnostic methods, the pathogenicity of non-jejuni/coli Campylobacter spp. such as C. upsaliensis and C. helveticus (isolated from dogs and cats) is uncertain. Galleria mellonella larvae are suitable models of the mammalian innate immune system and have been applied to C. jejuni studies. This study compared the pathogenicity of C. jejuni, C. upsaliensis, and C. helveticus isolates. Larvae inoculated with either C. upsaliensis or C. helveticus showed significantly higher survival than those inoculated with C. jejuni. All three Campylobacter species induced indistinguishable histopathological changes in the larvae. C. jejuni could be isolated from inoculated larvae up to eight days post-inoculation whereas C. upsaliensis and C. helveticus could only be isolated in the first two days. There was a significant variation in the hazard rate between batches of larvae, in Campylobacter strains, and in biological replicates as random effects, and in species and bacterial dose as fixed effects. The Galleria model is applicable to other Campylobacter spp. as well as C. jejuni, but may be subject to significant variation with all Campylobacter species. While C. upsaliensis and C. helveticus cannot be considered non-pathogenic, they are significantly less pathogenic than C. jejuni.

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

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.

Campylobacter jejuni Strains Associated with Wild Birds and Those Causing Human Disease in Six High-Use Recreational Waterways in New Zealand

Campylobacter jejuni, a leading cause of gastroenteritis worldwide, has been frequently isolated from recreational rivers and streams in New Zealand, yet the public health significance of this is unknown. This study uses molecular tools to improve our understanding of the epidemiology and sources of Campylobacter in recreational waterways, with a view to preventing human infection. Epidemiological and microbiological data were collected between 2005 and 2009 from six high-use recreational waterways in the Manawatu-Wanganui region of the North Island. Campylobacter spp. and C. jejuni were isolated from 33.2% and 20.4% of 509 samples, respectively. Isolation of Campylobacter was observed in both low and high river flows. After adjusting for the confounding effects of river flow, there was a significantly higher likelihood of isolating Campylobacter in the winter month of June compared to January. A high diversity of C. jejuni multilocus sequence types was seen, with the most commonly isolated being the water rail-associated ST-2381 (19/91 isolates [20.9%]), ST-1225 (8/91 isolates [8.8%]), and ST-45 (6/91 isolates [6.6%]). The ST-2381 was found in all rivers, while the most commonly isolated ST from human cases in New Zealand, the poultry-associated strain ST-474, was isolated only in one river. Although the majority of Campylobacter sequence types identified in river water were strains associated with wild birds that are rarely associated with human disease, poultry and ruminant-associated Campylobacter strains that are found in human infection were also identified and could present a public health risk.IMPORTANCE In 2016, there was a large-scale waterborne outbreak of campylobacteriosis in New Zealand, which was estimated to have affected over 5,000 people. This highlighted the need for a greater understanding of the sources of contamination of both surface and groundwater and risks associated with exposure to both drinking and recreational water. This study reports the prevalence and population structure of Campylobacter jejuni in six recreational waters of the Manawatu-Wanganui region of New Zealand and models the relationship between Campylobacter spp. and ruminant-associated Campylobacter and the parameters "sites," "months," and "river flow." Here, we demonstrate that both low and high river flows, month of the year, and recreational sites could influence the Campylobacter isolation from recreational waters. The presence of genotypes associated with human infection allowed us to describe potential risks associated with recreational waters.

Using Multinomial and Space-Time Permutation Models to Understand the Epidemiology of Infectious Bronchitis in California Between 2008 and 2012

Although infectious bronchitis virus (IBV) has been described as one of the most economically important viral respiratory diseases in poultry, there are few analyses of outbreaks that use spatial statistics. In order to better understand how the different genotypes of IBV behave spatially and temporally, we used geographic information system-based mapping coupled with spatial and spatial-temporal statistics to identify statistically significant clustering of multiple strains of infectious bronchitis (IB) between 2008 and 2012 in California. Specifically, space-time permutation and multinomial models were used to identify spatial and spatial-temporal clusters of various genotypes of IBV. Using time permutations (i.e., windows) spanning days to years, we identified three statistically significant ( P < 0.05) clusters. In contrast, multinomial models identified two statistically significant spatial-temporal clusters and one statistically significant spatial cluster. When comparing the space-time permutation and multinomial models against each other, we identified spatial and temporal overlap in two of the three statistically significant clusters. From a practical perspective, multinomial clustering approaches may be advantageous for studying IB because the model allows the different genotypes of IB to be independent nominal variables, thereby allowing for a more detailed spatial analysis. To that point, based on their risk ratios, the genotypes classified as vaccine-related were identified as the most significant contributor to two of the three mutinomial clusters. Additionally, statistically significant clusters were mapped and layered on a hot-spot analysis of commercial poultry farm density in order to qualitatively assess the relationship between farm density and clusters of IBV. Results showed that one of the three space-time permutations and one of the three multinomial clusters were spatially centered near the highest density farm areas, as determined by the hot-spot analysis.

Long-term Colonization by Campylobacter jejuni Within a Human Host: Evolution, Antimicrobial Resistance, and Adaptation

Background

Campylobacteriosis is inflammation of the gastrointestinal tract as a result of Campylobacter infection. Most campylobacteriosis cases are acute and self-limiting, with Campylobacter excretion ceasing a few weeks after symptoms cease. We identified a patient with fecal specimens positive for Campylobacter jejuni (ST45) intermittently during a 10-year period.

Methods

Sixteen Campylobacter isolates were collected from the patient during 2006-2016. The isolates' genomes were sequenced to determine their relatedness, and their antimicrobial susceptibility patterns and motility were measured to determine the effects of antibiotic therapy and long-term excretion on the Campylobacter population.

Results

Phylogenetic analyses estimated that the isolates shared a date of common ancestor between 1998 and 2006, coinciding with the onset of symptoms for the patient. Genomic analysis identified selection for changes in motility, and antimicrobial susceptibility testing suggested that the Campylobacter population developed resistance to several antibiotics coinciding with periods of antibiotic therapy.

Conclusions

The patient was consistently colonized with organisms from a Campylobacter population that adapted to the internal environment of the patient. Genomic and phylogenetic analyses can give insight into a patient's infection history and the effect of antimicrobial treatment on Campylobacter populations in this unusual situation of long-term colonization of an individual.

Extending statistical models for source attribution of zoonotic diseases: a study of campylobacteriosis

Preventing and controlling zoonoses through the design and implementation of public health policies requires a thorough understanding of transmission pathways. Modelling jointly the epidemiological data and genetic information of microbial isolates derived from cases provides a methodology for tracing back the source of infection. In this paper, the attribution probability for human cases of campylobacteriosis for each source, conditional on the extent to which each case resides in a rural compared to urban environment, is estimated. A model that incorporates genetic data and evolutionary processes is applied alongside a newly developed genetic-free model. We show that inference from each model is comparable except for rare microbial genotypes. Further, the effect of 'rurality' may be modelled linearly on the logit scale, with increasing rurality leading to the increasing likelihood of ruminant-sourced campylobacteriosis.

A six-year epidemiological surveillance study in Split-Dalmatia County, Croatia: urban versus rural differences in human campylobacteriosis incidence

Data from population-based laboratory surveillance were used to examine the epidemiological pattern of campylobacteriosis in a sentinel site, Split-Dalmatia County (SDC),Croatia, from 2007 to 2012, and to evaluate the association between disease incidence and demographic, geographical, climatic, agricultural, and microbiological factors. A total of 2658 laboratory-confirmed Campylobacter infections were recorded. Overall mean incidence was 96/100,000, ranging from 61/100,000 in rural to 131/100,000 in urban areas; rates were highest in the age group 0-4 years. Overall mean and age- and sex-specific incidences were significantly higher in urban versus rural areas (p < 0.01). The number of infections peaked in early summer, and was correlated with higher average monthly temperature (r = 0.58) and lower humidity (r = - 0.27), but not with precipitation. Incidence was not associated with agricultural activities. A distinct campylobacteriosis pattern with consistently higher urban versus rural incidence was observed, which may help formulate further preventive measures.

Shiga Toxin-Producing Escherichia coli O157 Shedding Dynamics in an Australian Beef Herd

Shiga toxin-producing Escherichia coli (STEC) O157 is an important foodborne pathogen that can be transmitted to humans both directly and indirectly from the feces of beef cattle, its primary reservoir. Numerous studies have investigated the shedding dynamics of E. coli O157 by beef cattle; however, the spatiotemporal trends of shedding are still not well understood. Molecular tools can increase the resolution through the use of strain typing to explore transmission dynamics within and between herds and identify strain-specific characteristics that may influence pathogenicity and spread. Previously, the shedding dynamics and molecular diversity, through the use of multilocus variable number of tandem repeat analysis (MLVA) of STEC O157, were separately investigated in an Australian beef herd over a 9-month study period. Variation in shedding was observed over time, and 33 MLVA types were identified. The study presented here combines the two datasets previously published with an aim to clarify the relationship between epidemiological variables and strain types. Three major genetic clusters (GCs) were identified that were significantly associated with the location of the cattle in different paddocks. No significant association between GCs and individual cow was observed. Results from this molecular epidemiological study provide evidence for herd-level clonal replacement over time that may have been triggered by movement to a new paddock. In conclusion, this study has provided further insight into STEC O157 shedding dynamics and pathogen transmission. Knowledge gaps remain regarding the relationship of strain types and the shedding dynamics of STEC O157 by beef cattle that could be further clarified through the use of whole-genome sequencing.

Inferring source attribution from a multiyear multisource data set of Salmonella in Minnesota

Salmonella enterica is a global health concern because of its widespread association with foodborne illness. Bayesian models have been developed to attribute the burden of human salmonellosis to specific sources with the ultimate objective of prioritizing intervention strategies. Important considerations of source attribution models include the evaluation of the quality of input data, assessment of whether attribution results logically reflect the data trends and identification of patterns within the data that might explain the detailed contribution of different sources to the disease burden. Here, more than 12,000 non-typhoidal Salmonella isolates from human, bovine, porcine, chicken and turkey sources that originated in Minnesota were analysed. A modified Bayesian source attribution model (available in a dedicated R package), accounting for non-sampled sources of infection, attributed 4,672 human cases to sources assessed here. Most (60%) cases were attributed to chicken, although there was a spike in cases attributed to a non-sampled source in the second half of the study period. Molecular epidemiological analysis methods were used to supplement risk modelling, and a visual attribution application was developed to facilitate data exploration and comprehension of the large multiyear data set assessed here. A large amount of within-source diversity and low similarity between sources was observed, and visual exploration of data provided clues into variations driving the attribution modelling results. Results from this pillared approach provided first attribution estimates for Salmonella in Minnesota and offer an understanding of current data gaps as well as key pathogen population features, such as serotype frequency, similarity and diversity across the sources. Results here will be used to inform policy and management strategies ultimately intended to prevent and control Salmonella infection in the state.

The EpiQuant Framework for Computing Epidemiological Concordance of Microbial Subtyping Data

A fundamental assumption in the use and interpretation of microbial subtyping results for public health investigations is that isolates that appear to be related based on molecular subtyping data are expected to share commonalities with respect to their origin, history, and distribution. Critically, there is currently no approach for systematically assessing the underlying epidemiology of subtyping results. Our aim was to develop a method for directly quantifying the similarity between bacterial isolates using basic sampling metadata and to develop a framework for computing the epidemiological concordance of microbial typing results. We have developed an analytical model that summarizes the similarity of bacterial isolates using basic parameters typically provided in sampling records, using a novel framework (EpiQuant) developed in the R environment for statistical computing. We have applied the EpiQuant framework to a data set comprising 654 isolates of the enteric pathogen Campylobacter jejuni from Canadian surveillance data in order to examine the epidemiological concordance of clusters obtained by using two leading C. jejuni subtyping methods. The EpiQuant framework can be used to directly quantify the similarity of bacterial isolates based on basic sample metadata. These results can then be used to assess the concordance between microbial epidemiological and molecular data, facilitating the objective assessment of subtyping method performance and paving the way for the improved application of molecular subtyping data in investigations of infectious disease.

Serotypes and Antimicrobial Resistance in Salmonella enterica Recovered from Clinical Samples from Cattle and Swine in Minnesota, 2006 to 2015

Salmonellosis remains one of the leading causes of foodborne disease worldwide despite preventive efforts at various stages of the food production chain. The emergence of multi-drug resistant (MDR) non-typhoidal Salmonella enterica represents an additional challenge for public health authorities. Food animals are considered a major reservoir and potential source of foodborne salmonellosis; thus, monitoring of Salmonella strains in livestock may help to detect emergence of new serotypes/MDR phenotypes and to gain a better understanding of Salmonella epidemiology. For this reason, we analyzed trends over a nine-year period in serotypes, and antimicrobial resistance, of Salmonella isolates recovered at the Minnesota Veterinary Diagnostic Laboratory (MVDL) from swine (n = 2,537) and cattle (n = 1,028) samples. Prevalence of predominant serotypes changed over time; in swine, S. Typhimurium and S. Derby decreased and S. Agona and S. 4,5,12:i:- increased throughout the study period. In cattle, S. Dublin, S. Montevideo and S. Cerro increased and S. Muenster became less frequent. Median minimum inhibitory concentration (MIC) values and proportion of antibiotic resistant isolates were higher for those recovered from swine compared with cattle, and were particularly high for certain antibiotic-serotype combinations. The proportion of resistant swine isolates was also higher than observed in the NARMS data, probably due to the different cohort of animals represented in each dataset. Results provide insight into the dynamics of antimicrobial resistant Salmonella in livestock in Minnesota, and can help to monitor emerging trends in antimicrobial resistance.

Timely Reporting and Interactive Visualization of Animal Health and Slaughterhouse Surveillance Data in Switzerland

The reporting of outputs from health surveillance systems should be done in a near real-time and interactive manner in order to provide decision makers with powerful means to identify, assess, and manage health hazards as early and efficiently as possible. While this is currently rarely the case in veterinary public health surveillance, reporting tools do exist for the visual exploration and interactive interrogation of health data. In this work, we used tools freely available from the Google Maps and Charts library to develop a web application reporting health-related data derived from slaughterhouse surveillance and from a newly established web-based equine surveillance system in Switzerland. Both sets of tools allowed entry-level usage without or with minimal programing skills while being flexible enough to cater for more complex scenarios for users with greater programing skills. In particular, interfaces linking statistical softwares and Google tools provide additional analytical functionality (such as algorithms for the detection of unusually high case occurrences) for inclusion in the reporting process. We show that such powerful approaches could improve timely dissemination and communication of technical information to decision makers and other stakeholders and could foster the early-warning capacity of animal health surveillance systems.

Seasonality ofCampylobacter jejuniisolates associated with human campylobacteriosis in the Manawatu region, New Zealand

A 9-year time-series of genotyped human campylobacteriosis cases from the Manawatu region of New Zealand was used to investigate strain-type seasonality. The data were collected from 2005 to 2013 and the samples were multi-locus sequence-typed (MLST). The four most prevalent clonal complexes (CCs), consisting of 1215 isolates, were CC48, CC21, CC45 and CC61. Seasonal decomposition and Poisson regression with autocorrelated errors, were used to display and test for seasonality of the most prevalent CCs. Of the four examined CCs, only CC45 showed a marked seasonal (summer) peak. The association of CC45 with summer peaks has been observed in other temperate countries, but has previously not been identified in New Zealand. This is the first in-depth study over a long time period employing MLST data to examine strain-type-associated seasonal patterns ofC. jejuniinfection in New Zealand.

Evaluation of the representativeness of a sentinel surveillance site for campylobacteriosis

It is important to assess the suitability of sentinel sites for human disease; however, there have been few publications documenting the process of formal evaluation. We describe an approach to examining the representativeness of a single sentinel site employed for campylobacteriosis surveillance and source attribution, utilizing a selection of data sources and statistical comparisons of demographic, epidemiological and pathogen genotyping data across selected regions of New Zealand. Our findings showed that while this region captured the national variability in many variables, for example by containing sizable urban and rural populations, the relative frequency of these features did vary from other regions of New Zealand. We discuss the value of choosing a sentinel site that represents the national distribution of key variables, compared to a site that captures the broad features of the wider population, but provides greater power for the monitoring of sub-populations.

Current methods for molecular typing of Campylobacter species

Campylobacter remains one of the most common bacterial causes of gastroenteritis worldwide. Tracking sources of this organism is challenging due to the large numbers of human cases, and the prevalence of this organism throughout the environment due to growth in a wide range of animal species. Many molecular subtyping methods have been developed to characterize Campylobacter species, but only a few are commonly used in molecular epidemiology studies. This review examines the applicability of these methods, as well as the role that emerging whole genome sequencing technologies will play in tracking sources of Campylobacter spp. infection.

Comparison of time series models for predicting campylobacteriosis risk in New Zealand

Predicting campylobacteriosis cases is a matter of considerable concern in New Zealand, after the number of the notified cases was the highest among the developed countries in 2006. Thus, there is a need to develop a model or a tool to predict accurately the number of campylobacteriosis cases as the Microbial Risk Assessment Model used to predict the number of campylobacteriosis cases failed to predict accurately the number of actual cases. We explore the appropriateness of classical time series modelling approaches for predicting campylobacteriosis. Finding the most appropriate time series model for New Zealand data has additional practical considerations given a possible structural change, that is, a specific and sudden change in response to the implemented interventions. A univariate methodological approach was used to predict monthly disease cases using New Zealand surveillance data of campylobacteriosis incidence from 1998 to 2009. The data from the years 1998 to 2008 were used to model the time series with the year 2009 held out of the data set for model validation. The best two models were then fitted to the full 1998-2009 data and used to predict for each month of 2010. The Holt-Winters (multiplicative) and ARIMA (additive) intervention models were considered the best models for predicting campylobacteriosis in New Zealand. It was noticed that the prediction by an additive ARIMA with intervention was slightly better than the prediction by a Holt-Winter multiplicative method for the annual total in year 2010, the former predicting only 23 cases less than the actual reported cases. It is confirmed that classical time series techniques such as ARIMA with intervention and Holt-Winters can provide a good prediction performance for campylobacteriosis risk in New Zealand. The results reported by this study are useful to the New Zealand Health and Safety Authority's efforts in addressing the problem of the campylobacteriosis epidemic.

Whole-genome comparison of two Campylobacter jejuni isolates of the same sequence type reveals multiple loci of different ancestral lineage

Campylobacter jejuni ST-474 is the most important human enteric pathogen in New Zealand, and yet this genotype is rarely found elsewhere in the world. Insight into the evolution of this organism was gained by a whole genome comparison of two ST-474, flaA SVR-14 isolates and other available C. jejuni isolates and genomes. The two isolates were collected from different sources, human (H22082) and retail poultry (P110b), at the same time and from the same geographical location. Solexa sequencing of each isolate resulted in 1.659 Mb (H22082) and 1.656 Mb (P110b) of assembled sequences within 28 (H22082) and 29 (P110b) contigs. We analysed 1502 genes for which we had sequences within both ST-474 isolates and within at least one of 11 C. jejuni reference genomes. Although 94.5% of genes were identical between the two ST-474 isolates, we identified 83 genes that differed by at least one nucleotide, including 55 genes with non-synonymous substitutions. These covered 101 kb and contained 672 point differences. We inferred that 22 (3.3%) of these differences were due to mutation and 650 (96.7%) were imported via recombination. Our analysis estimated 38 recombinant breakpoints within these 83 genes, which correspond to recombination events affecting at least 19 loci regions and gives a tract length estimate of 2 kb. This includes a 12 kb region displaying non-homologous recombination in one of the ST-474 genomes, with the insertion of two genes, including ykgC, a putative oxidoreductase, and a conserved hypothetical protein of unknown function. Furthermore, our analysis indicates that the source of this recombined DNA is more likely to have come from C. jejuni strains that are more closely related to ST-474. This suggests that the rates of recombination and mutation are similar in order of magnitude, but that recombination has been much more important for generating divergence between the two ST-474 isolates.