Commensal Escherichia coli Strains of Bovine Origin Competitively Mitigated Escherichia coli O157:H7 in a Gnotobiotic Murine Intestinal Colonization Model with or without Physiological Stress

Cattle are a primary reservoir of enterohemorrhagic Escherichia coli (EHEC) O157:H7. Currently, there are no effective methods of eliminating this important zoonotic pathogen from cattle, and colonization resistance in relation to EHEC O157:H7 in cattle is poorly understood. We developed a gnotobiotic EHEC O157:H7 murine model to examine aspects of the cattle pathogen-microbiota interaction, and to investigate competitive suppression of EHEC O157:H7 by 18 phylogenetically distinct commensal E. coli strains of bovine origin. As stress has been suggested to influence enteric colonization by EHEC O157:H7 in cattle, corticosterone administration (±) to incite a physiological stress response was included as an experimental variable. Colonization of the intestinal tract (IT) of mice by the bovine EHEC O157:H7 strain, FRIK-2001, mimicked characteristics of bovine IT colonization. In this regard, FRIK-2001 successfully colonized the IT and temporally incited minimal impacts on the host relative to other EHEC O157:H7 strains, including on the renal metabolome. The presence of the commensal E. coli strains decreased EHEC O157:H7 densities in the cecum, proximal colon, and distal colon. Moreover, histopathologic changes and inflammation markers were reduced in the distal colon of mice inoculated with commensal E. coli strains (both propagated separately and communally). Although stress induction affected the behavior of mice, it did not influence EHEC O157:H7 densities or disease. These findings support the use of a gnotobiotic murine model of enteric bovine EHEC O157:H7 colonization to better understand pathogen-host-microbiota interactions toward the development of effective on-farm mitigations for EHEC O157:H7 in cattle, including the identification of bacteria capable of competitively colonizing the IT.

Combining analytical epidemiology and genomic surveillance to identify risk factors associated with the spread of antimicrobial resistance in Salmonella enterica subsp. enterica serovar Heidelberg

Antimicrobial resistance (AMR) has become a critical threat to public health worldwide. The use of antimicrobials in food and livestock agriculture, including the production of poultry, is thought to contribute to the dissemination of antibiotic resistant bacteria (ARB) and the genes and plasmids that confer the resistant phenotype (ARG). However, the relative contribution of each of these processes to the emergence of resistant pathogens in poultry production and their potential role in the transmission of resistant pathogens in human infections, requires a deeper understanding of the dynamics of ARB and ARG in food production and the factors involved in the increased risk of transmission.

Systematic Evaluation of Whole-Genome Sequencing Based Prediction of Antimicrobial Resistance in Campylobacter jejuni and C. coli

The increasing prevalence of antimicrobial resistance (AMR) in Campylobacter spp. is a global concern. This study evaluated the use of whole-genome sequencing (WGS) to predict AMR in Campylobacter jejuni and C. coli. A panel of 271 isolates recovered from Canadian poultry was used to compare AMR genotype to antimicrobial susceptibility testing (AST) results (azithromycin, ciprofloxacin, erythromycin, gentamicin, tetracycline, florfenicol, nalidixic acid, telithromycin, and clindamycin). The presence of antibiotic resistance genes (ARGs) was determined for each isolate using five computational approaches to evaluate the effect of: ARG screening software, input data (i.e., raw reads, draft genome assemblies), genome coverage and genome assembly software. Overall, concordance between the genotype and phenotype was influenced by the computational pipelines, level of genome coverage and the type of ARG but not by input data. For example, three of the pipelines showed a 99% agreement between detection of a tet(O) gene and tetracycline resistance, whereas agreement between the detection of tet(O) and TET resistance was 98 and 93% for two pipelines. Overall, higher levels of genome coverage were needed to reliably detect some ARGs; for example, at 15X coverage a tet(O) gene was detected in >70% of the genomes, compared to <60% of the genomes for bla(OXA). No genes associated with florfenicol or gentamicin resistance were found in the set of strains included in this study, consistent with AST results. Macrolide and fluoroquinolone resistance was associated 100% with mutations in the 23S rRNA (A2075G) and gyrA (T86I) genes, respectively. A lower association between a A2075G 23S rRNA gene mutation and resistance to clindamycin and telithromycin (92.8 and 78.6%, respectively) was found. While WGS is an effective approach to predicting AMR in Campylobacter, this study demonstrated the impact that computational pipelines, genome coverage and the genes can have on the reliable identification of an AMR genotype.

Rapid and accurate SNP genotyping of clonal bacterial pathogens with BioHansel

Hierarchical genotyping approaches can provide insights into the source, geography and temporal distribution of bacterial pathogens. Multiple hierarchical SNP genotyping schemes have previously been developed so that new isolates can rapidly be placed within pre-computed population structures, without the need to rebuild phylogenetic trees for the entire dataset. This classification approach has, however, seen limited uptake in routine public health settings due to analytical complexity and the lack of standardized tools that provide clear and easy ways to interpret results. The BioHansel tool was developed to provide an organism-agnostic tool for hierarchical SNP-based genotyping. The tool identifies split k-mers that distinguish predefined lineages in whole genome sequencing (WGS) data using SNP-based genotyping schemes. BioHansel uses the Aho-Corasick algorithm to type isolates from assembled genomes or raw read sequence data in a matter of seconds, with limited computational resources. This makes BioHansel ideal for use by public health agencies that rely on WGS methods for surveillance of bacterial pathogens. Genotyping results are evaluated using a quality assurance module which identifies problematic samples, such as low-quality or contaminated datasets. Using existing hierarchical SNP schemes for Mycobacterium tuberculosis and Salmonella Typhi, we compare the genotyping results obtained with the k-mer-based tools BioHansel and SKA, with those of the organism-specific tools TBProfiler and genotyphi, which use gold-standard reference-mapping approaches. We show that the genotyping results are fully concordant across these different methods, and that the k-mer-based tools are significantly faster. We also test the ability of the BioHansel quality assurance module to detect intra-lineage contamination and demonstrate that it is effective, even in populations with low genetic diversity. We demonstrate the scalability of the tool using a dataset of ~8100 S. Typhi public genomes and provide the aggregated results of geographical distributions as part of the tool’s output. BioHansel is an open source Python 3 application available on PyPI and Conda repositories and as a Galaxy tool from the public Galaxy Toolshed. In a public health context, BioHansel enables rapid and high-resolution classification of bacterial pathogens with low genetic diversity.

Epidemiology of Campylobacter jejuni in raccoons (Procyon lotor) on swine farms and in conservation areas in southern Ontario

Campylobacter is a leading cause of foodborne illness in humans worldwide. Sources of infection are often difficult to identify, and are, generally, poorly understood. Recent work suggests that wildlife may represent a source of Campylobacter for human infections. Using a repeated cross-sectional study design, raccoons were trapped on five swine farms and five conservation areas in southern Ontario from 2011 to 2013. Our objectives were to: (a) assess the impact of seasonal, climatic, location, annual and raccoon demographic factors on the occurrence of Campylobacter jejuni in these animals; and (b) identify clusters of C. jejuni in space, time and space-time using spatial scan statistics. Multi-level multivariable logistic regression was used to examine the odds of isolating C. jejuni, with site and animal modelled as random intercepts. The following independent variables were examined: raccoon age and sex, year, location type, season, temperature and rainfall. A total of 1,096 samples were obtained from 627 raccoons; 46.3% were positive for C. jejuni. The following interactions and their main effects were significant (p < .05) and retained in the final model: season × temperature, year × rainfall, year × temperature. Based on the results from our multivariable model and spatial scan statistics, climatic variables (i.e. rainfall, temperature and season) were associated with the carriage of C. jejuni by raccoons, but the effects were not consistent, and varied by location and year. Although raccoons may pose a zoonotic risk due to their carriage of Campylobacter, further work is required to characterize the transmission and movement of this microorganism within the ecosystem.

Rates of fluoroquinolone resistance in domestically acquired Campylobacter jejuni are increasing in people living within a model study location in Canada

Antimicrobial resistance was evaluated in Campylobacter jejuni isolated from 1291 diarrheic people over a 15-year period (2004-2018) in southwestern Alberta, a model location in Canada with a high rate of campylobacteriosis. The prevalence of resistance to chloramphenicol, clindamycin, erythromycin, and gentamicin was low during the examination period (≤4.8%). Resistance to tetracycline remained consistently high (41.6%-65.1%), and resistance was primarily conferred by plasmid-borne tetO (96.2%). Resistance rates to ciprofloxacin and nalidixic acid increased substantially over the examination period, with a maximal fluoroquinolone resistance (FQR) prevalence of 28.9% in 2016. The majority of C. jejuni isolates resistant to ciprofloxacin (93.9%) contained a C257T single nucleotide polymorphism within the gyrA chromosomal gene. Follow up with infected people indicated that the observed increase in FQR was primarily due to domestically acquired infections. Moreover, the majority of FQ-resistant C. jejuni subtypes (82.6%) were endemic in Canada, primarily linked to cattle and chicken reservoirs; 18.4% of FQ-resistant isolates were assigned to three subtypes, predominantly associated with cattle. Study findings indicate the need to prioritize FQR monitoring in C. jejuni infections in Canada and to elucidate the dynamics of the emergence and transmission of resistant C. jejuni strains within and from cattle and chicken reservoirs.

"These Aren't the Strains You're Looking for": Recovery Bias of Common Campylobacter jejuni Subtypes in Mixed Cultures

Microbiological surveillance of the food chain plays a critical role in improving our understanding of the distribution and circulation of food-borne pathogens along the farm to fork continuum toward the development of interventions to reduce the burden of illness. The application of molecular subtyping to bacterial isolates collected through surveillance has led to the identification of strains posing the greatest risk to public health. Past evidence suggests that enrichment methods for Campylobacter jejuni, a leading bacterial foodborne pathogen worldwide, may lead to the differential recovery of subtypes, obscuring our ability to infer the composition of a mixed-strain sample and potentially biasing prevalence estimates in surveillance data. To assess the extent of potential selection bias resulting from enrichment-based isolation methods, we compared enrichment and non-enrichment isolation of mixed subtype cultures of C. jejuni, followed by subtype-specific enumeration using both colony plate-counts and digital droplet PCR. Results differed from the null hypothesis that similar proportions of C. jejuni subtypes are recovered from both methods. Our results also indicated a significant effect of subtype prevalence on isolation frequency post-recovery, with the recovery of more common subtypes being consistently favored. This bias was exacerbated when an enrichment step was included in the isolation procedure. Taken together, our results emphasize the importance of selecting multiple colonies per sample, and where possible, the use of both enrichment and non-enrichment isolation procedures to maximize the likelihood of recovering multiple subtypes present in a sample. Moreover, the effects of subtype-specific recovery bias should be considered in the interpretation of strain prevalence data toward improved risk assessment from microbiological surveillance data.

Campylobacter jejuni Strain Dynamics in a Raccoon (Procyon lotor) Population in Southern Ontario, Canada: High Prevalence and Rapid Subtype Turnover

Free-ranging wildlife are increasingly recognized as potential reservoirs of disease-causing Campylobacter species such as C. jejuni and C. coli. Raccoons (Procyon lotor), which live at the interface of rural, urban, and more natural environments, are ideal subjects for exploring the potential role that wildlife play in the epidemiology of campylobacteriosis. We studied the prevalence and genetic diversity of Campylobacter from live-captured raccoons on five swine farms and five conservation areas in southwest Ontario. From 2011 to 2013, we collected fecal swabs (n = 1,096) from raccoons, and (n = 50) manure pit samples from the swine farm environment. We subtyped the resulting Campylobacter isolates (n = 581) using Comparative Genomic Fingerprinting (CGF) and 114 distinct subtypes were observed, including 96 and 18 subtypes among raccoon and manure pit isolates, respectively. Campylobacter prevalence in raccoons was 46.3%, with 98.7% of isolates recovered identified as C. jejuni. Novel raccoon-specific CGF subtypes (n = 40/96) accounted for 24.6% (n = 143/581) of Campylobacter isolates collected in this study. Our results also show that C. jejuni is readily acquired and lost in this wild raccoon population and that a high Campylobacter prevalence is observed despite transient carriage typically lasting 30 days or fewer. Moreover, although raccoons appeared to be colonized by species-adapted subtypes, they also harbored agriculture-associated genotypes that accounted for the majority of isolates observed (66.4%) and that are strongly associated with human infections. This suggests that raccoons may act as vectors in the transmission of clinically-relevant C. jejuni subtypes at the interface of rural, urban, and more natural environments.

Integrating Whole-Genome Sequencing Data Into Quantitative Risk Assessment of Foodborne Antimicrobial Resistance: A Review of Opportunities and Challenges

Whole-genome sequencing (WGS) will soon replace traditional phenotypic methods for routine testing of foodborne antimicrobial resistance (AMR). WGS is expected to improve AMR surveillance by providing a greater understanding of the transmission of resistant bacteria and AMR genes throughout the food chain, and therefore support risk assessment activities. At this stage, it is unclear how WGS data can be integrated into quantitative microbial risk assessment (QMRA) models and whether their integration will impact final risk estimates or the assessment of risk mitigation measures. This review explores opportunities and challenges of integrating WGS data into QMRA models that follow the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR. We describe how WGS offers an opportunity to enhance the next-generation of foodborne AMR QMRA modeling. Instead of considering all hazard strains as equally likely to cause disease, WGS data can improve hazard identification by focusing on those strains of highest public health relevance. WGS results can be used to stratify hazards into strains with similar genetic profiles that are expected to behave similarly, e.g., in terms of growth, survival, virulence or response to antimicrobial treatment. The QMRA input distributions can be tailored to each strain accordingly, making it possible to capture the variability in the strains of interest while decreasing the uncertainty in the model. WGS also allows for a more meaningful approach to explore genetic similarity among bacterial populations found at successive stages of the food chain, improving the estimation of the probability and magnitude of exposure to AMR hazards at point of consumption. WGS therefore has the potential to substantially improve the utility of foodborne AMR QMRA models. However, some degree of uncertainty remains in relation to the thresholds of genetic similarity to be used, as well as the degree of correlation between genotypic and phenotypic profiles. The latter could be improved using a functional approach based on prediction of microbial behavior from a combination of 'omics' techniques (e.g., transcriptomics, proteomics and metabolomics). We strongly recommend that methodologies to incorporate WGS data in risk assessment be included in any future revision of the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR.

A strain comparison of Campylobacter isolated from retail poultry and human clinical cases in Atlantic Canada

Campylobacter is the leading cause of food-borne bacterial disease in Canada and many developed countries. One of the most common sources of human campylobacteriosis is considered to be the consumption or handling of raw or undercooked poultry. To date, few Canadian studies have investigated both the prevalence of Campylobacter on retail poultry and its potential impact on human clinical cases. The objective of this study was to evaluate the prevalence of Campylobacter spp. at the retail level and the correlation between subtypes recovered from chicken and those recovered from human clinical cases within the province of Nova Scotia, Canada. From this study 354 human clinical isolates were obtained from provincial hospital laboratories and a total of 480 packages of raw poultry cuts were sampled from retail outlets, yielding 312 isolates (65%), of all which were subtyped using comparative genomic fingerprinting (CGF). Of the 312 chicken isolates, the majority of isolates were C. jejuni (91.7%), followed by C. coli (7.7%) and C. lari (0.6%). Using CGF to subtype C. jejuni and C. coli isolates, 99 and 152 subtypes were recovered from chicken and clinical cases, respectively. The most prevalent human and chicken subtypes found in NS are similar to those observed nationally; indicating that the Campylobacter from this study appear to reflect of the profile of Campylobacter subtypes circulating nationally. Of the subtypes observed, only 36 subtypes were common between the two groups, however, these subtypes represented 48.3% of the clinical isolates collected. The findings from this study provides evidence that in Nova Scotia, retail poultry can act as a reservoir for Campylobacter subtypes that have been implicated in human illness.

Enhanced microbiological surveillance reveals that temporal case clusters contribute to the high rates of campylobacteriosis in a model agroecosystem

Infections by pathogenic Campylobacter species were determined in diarrheic (n = 2,217) and non-diarrheic control (n = 104) people in Southwestern Alberta (SWA), Canada over a 1-year period using specialized and conventional isolation, and direct PCR. Overall, 9.9% of diarrheic individuals were positive for C. jejuni (9.1%), C. upsaliensis (0.6%), and C. coli (0.5%). No C. lari was detected. Four diarrheic individuals were co-infected with C. jejuni and C. coli, and four different individuals were co-infected with C. jejuni and C. upsaliensis. Two control individuals were positive for C. jejuni. Approximately 50% of stools containing C. jejuni and/or C. coli were deemed negative by conventional isolation. Direct PCR for C. jejuni was less effective than culture-based detection. Most C. jejuni infections occurred in people living in the urban centers, but the prevalence of the bacterium was lower in females than males living in urban locations, and both males and females living in rural locations. Although C. jejuni was detected throughout the year, a trend for higher infection rates was observed in the late spring to early fall with a peak in August. Forty-six C. jejuni subtype clusters were identified, including 44 temporal case clusters attributed to 28 subtype groupings. The majority of infections (70.3%) were linked to subtypes associated with beef cattle. We conclude that many occurrences of pathogenic Campylobacter species were not detected by the conventional laboratory methodology, and temporal case clusters of C. jejuni subtypes associated with cattle contribute to the high rates of campylobacteriosis in SWA.

Frequent Implication of Multistress-Tolerant Campylobacter jejuni in Human Infections

Campylobacter jejuni, a major cause of bacterial foodborne illnesses, is considered highly susceptible to environmental stresses. In this study, we extensively investigated the stress tolerance of 121 clinical strains of C. jejuni against 5 stress conditions (aerobic stress, disinfectant exposure, freeze-thaw, heat treatment, and osmotic stress) that this pathogenic bacterium might encounter during foodborne transmission to humans. In contrast to our current perception about high stress sensitivity of C. jejuni, a number of clinical strains of C. jejuni were highly tolerant to multiple stresses. We performed population genetics analysis by using comparative genomic fingerprinting and showed that multistress-tolerant strains of C. jejuni constituted distinct clades. The comparative genomic fingerprinting subtypes belonging to multistress-tolerant clades were more frequently implicated in human infections than those in stress-sensitive clades. We identified unique stress-tolerant C. jejuni clones and showed the role of stress tolerance in human campylobacteriosis.

Genome-wide identification of geographical segregated genetic markers in Salmonella enterica serovar Typhimurium variant 4,[5],12:i:

Salmonella enterica ser. Typhimurium monophasic variant 4,[5],12:i:- has been associated with food-borne epidemics worldwide and swine appeared to be the main reservoir in most of the countries of isolation. However, the monomorphic nature of this serovar has, so far, hindered identification of the source due to expansion of clonal lineages in multiple hosts and food producing systems. Since geographically structured genetic signals can shape bacterial populations, identification of biogeographical markers in S. 1,4,[5],12:i:- genomes can contribute to improving source attribution. In this study, the phylogeographical structure of 148 geographically and temporally related Italian S. 1,4,[5],12:i:- has been investigated. The Italian isolates belong to a large population of clonal S. Typhimurium/1,4,[5],12:i:- isolates collected worldwide in two decades showing up to 2.5% of allele differences. Phylogenetic reconstruction revealed that isolates from the same geographical origin form highly supported monophyletic groups, suggesting discrete geographical segregation. These monophyletic groups are characterized by the gene content of a large sopE-containing prophage. Within this prophage, genome-wide comparison identified several genes overrepresented in strains of Italian origin. This suggests that certain lineages may be characterized by the acquisition of specific accessory genetic markers useful for improving identification of the source in ongoing epidemics.

A repeated cross-sectional study of the epidemiology of Campylobacter and antimicrobial resistant Enterobacteriaceae in free-living Canada geese in Guelph, Ontario, Canada

From May through October 2016, we conducted a repeated cross-sectional study examining the effects of temporal, spatial, flock and demographic factors (i.e. juvenile vs. adult) on the prevalence of Campylobacter and antimicrobial resistant Enterobacteriaceae among 344 fresh faecal samples collected from Canada geese (Branta canadensis) from four locations where birds nested in Guelph, Ontario, Canada. The overall prevalence of Campylobacter among all fresh faecal samples was 9.3% and was greatest in the fall when these birds became more mobile following the nesting season. Based on 40 gene comparative genomic fingerprinting (CGF40), the increase in prevalence noted in the fall was matched by an increase in the number of unique CGF40 subtypes identified. Resistance to colistin was detected most commonly, in 6% of Escherichia coli isolates, and was highest in the late summer months. All colistin-resistant isolates were negative for the mcr-1 to mcr-5 genes; a chromosomal resistance mechanism (PmrB) was identified in all of these isolates. The prevalence of samples with E. coli exhibiting multi-class resistance or extended spectrum beta-lactamase was low (i.e. <2% of samples). The intra-class correlation coefficients, estimated from the variance components of multilevel logistic regression models, indicated that the shedding of Campylobacter and antimicrobial resistant E. coli among geese within a flock (i.e. birds collected from the same site on the same day) was moderately correlated. Spatial, temporal, and spatiotemporal clusters identified using the spatial scan statistic, largely supported the findings from our multi-level models. Salmonella was not isolated from any of the fresh faecal samples collected suggesting that its prevalence in this population of birds was very low.

Therapeutic administration of enrofloxacin in mice does not select for fluoroquinolone resistance inCampylobacter jejuni

Enrofloxacin is registered for therapeutic use in beef cattle to treat bovine respiratory disease in Canada. A murine model was used to experimentally examine the impact of therapeutic administration of enrofloxacin on fluoroquinolone resistance development in Campylobacter jejuni. Administration of enrofloxacin to mice via subcutaneous injection or per os routes resulted in equivalent levels of bioactive enrofloxacin within the intestine, but bioactivity was short-lived (<48 h after cessation). Enrofloxacin administration did not affect densities of total bacteria, Firmicutes, or Bacteroidetes in digesta and had modest impacts on densities of Enterobacteriaceae. All mice inoculated with C. jejuni NCTC 11168 became persistently colonized by the bacterium. Enrofloxacin reduced C. jejuni cell densities within the cecal and colonic digesta for all treatments, and densities shed in feces as a function of antibiotic duration. None of the C. jejuni isolates recovered from mice after administration of enrofloxacin (n = 260) developed resistance to ciprofloxacin regardless of method or duration of administration. Furthermore, only modest shifts in the minimum inhibitory concentration of the isolates by treatment were noted. The study findings indicate that the risk posed by short-term subcutaneous administration of enrofloxacin for the development of fluoroquinolone resistance in mammals is low.

Epidemiology of Campylobacter, Salmonella and antimicrobial resistant Escherichia coli in free-living Canada geese (Branta canadensis) from three sources in southern Ontario

Antimicrobial resistant bacteria and zoonotic pathogens have previously been isolated from Canada geese. We examined the prevalence of three enteric bacteria (i.e. Campylobacter, Salmonella, Escherichia coli) among Canada geese from three sampling sources in southern Ontario from 2013 through 2015. Samples were obtained by convenience from hunting groups, diagnostic birds submitted for post-mortem, and fresh faeces from live birds in parks. Escherichia coli isolates were isolated and tested for susceptibility to 15 antimicrobials using the Canadian Integrated Program for Antimicrobial Resistance Surveillance test panel. The prevalences of Salmonella, Campylobacter and E. coli were 0%, 11.2% and 72.6%, respectively. Among E. coli isolates, 7.9% were resistant to ≥1 class of antimicrobials and 5.6% were resistant to ≥2 classes of antimicrobials, with some including resistance to antimicrobials of highest importance in human medicine. A significant association between season and E. coli resistance among samples from live birds was noted; summer samples had no resistant E. coli isolates, whereas spring samples demonstrated the highest prevalence of E. coli resistant to ≥1 class of antimicrobials (20.0%) among all sources. In addition, Campylobacter coli were only isolated from the spring faecal samples. Flock-level clustering was an important statistical consideration, as flock was a significant random effect in all but two of our models. Detection of Campylobacter and antimicrobial resistant E. coli from Canada geese suggests that these birds may play a role in disseminating these organisms within the environment.

Comprehensive assessment of the quality of Salmonella whole genome sequence data available in public sequence databases using the Salmonella in silico Typing Resource (SISTR)

Public health and food safety institutions around the world are adopting whole genome sequencing (WGS) to replace conventional methods for characterizing Salmonella for use in surveillance and outbreak response. Falling costs and increased throughput of WGS have resulted in an explosion of data, but questions remain as to the reliability and robustness of the data. Due to the critical importance of serovar information to public health, it is essential to have reliable serovar assignments available for all of the Salmonella records. The current study used a systematic assessment and curation of all Salmonella in the sequence read archive (SRA) to assess the state of the data and their utility. A total of 67 758 genomes were assembled de novo and quality-assessed for their assembly metrics as well as species and serovar assignments. A total of 42 400 genomes passed all of the quality criteria but 30.16 % of genomes were deposited without serotype information. These data were used to compare the concordance of reported and predicted serovars for two in silico prediction tools, multi-locus sequence typing (MLST) and the Salmonella in silico Typing Resource (SISTR), which produced predictions that were fully concordant with 87.51 and 91.91 % of the tested isolates, respectively. Concordance of in silico predictions increased when serovar variants were grouped together, 89.25 % for MLST and 94.98 % for SISTR. This study represents the first large-scale validation of serovar information in public genomes and provides a large validated set of genomes, which can be used to benchmark new bioinformatics tools.

Source attribution of human campylobacteriosis at the point of exposure by combining comparative exposure assessment and subtype comparison based on comparative genomic fingerprinting

Human campylobacteriosis is a common zoonosis with a significant burden in many countries. Its prevention is difficult because humans can be exposed to Campylobacter through various exposures: foodborne, waterborne or by contact with animals. This study aimed at attributing campylobacteriosis to sources at the point of exposure. It combined comparative exposure assessment and microbial subtype comparison with subtypes defined by comparative genomic fingerprinting (CGF). It used isolates from clinical cases and from eight potential exposure sources (chicken, cattle and pig manure, retail chicken, beef, pork and turkey meat, and surface water) collected within a single sentinel site of an integrated surveillance system for enteric pathogens in Canada. Overall, 1518 non-human isolates and 250 isolates from domestically-acquired human cases were subtyped and their subtype profiles analyzed for source attribution using two attribution models modified to include exposure. Exposure values were obtained from a concurrent comparative exposure assessment study undertaken in the same area. Based on CGF profiles, attribution was possible for 198 (79%) human cases. Both models provide comparable figures: chicken meat was the most important source (65-69% of attributable cases) whereas exposure to cattle (manure) ranked second (14-19% of attributable cases), the other sources being minor (including beef meat). In comparison with other attributions conducted at the point of production, the study highlights the fact that Campylobacter transmission from cattle to humans is rarely meat borne, calling for a closer look at local transmission from cattle to prevent campylobacteriosis, in addition to increasing safety along the chicken supply chain.

A Genome-Wide Association Study to Identify Diagnostic Markers for Human Pathogenic Campylobacter jejuni Strains

Campylobacter jejuni is a leading human enteric pathogen worldwide and despite an improved understanding of its biology, ecology, and epidemiology, limited tools exist for identifying strains that are likely to cause disease. In the current study, we used subtyping data in a database representing over 24,000 isolates collected through various surveillance projects in Canada to identify 166 representative genomes from prevalent C. jejuni subtypes for whole genome sequencing. The sequence data was used in a genome-wide association study (GWAS) aimed at identifying accessory gene markers associated with clinically related C. jejuni subtypes. Prospective markers (n = 28) were then validated against a large number (n = 3,902) of clinically associated and non-clinically associated genomes from a variety of sources. A total of 25 genes, including six sets of genetically linked genes, were identified as robust putative diagnostic markers for clinically related C. jejuni subtypes. Although some of the genes identified in this study have been previously shown to play a role in important processes such as iron acquisition and vitamin B₅ biosynthesis, others have unknown function or are unique to the current study and warrant further investigation. As few as four of these markers could be used in combination to detect up to 90% of clinically associated isolates in the validation dataset, and such markers could form the basis for a screening assay to rapidly identify strains that pose an increased risk to public health. The results of the current study are consistent with the notion that specific groups of C. jejuni strains of interest are defined by the presence of specific accessory genes.

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole‐genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non‐UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.

The Validation and Implications of Using Whole Genome Sequencing as a Replacement for Traditional Serotyping for a National Salmonella Reference Laboratory

Salmonella serotyping remains the gold-standard tool for the classification of Salmonella isolates and forms the basis of Canada’s national surveillance program for this priority foodborne pathogen. Public health officials have been increasingly looking toward whole genome sequencing (WGS) to provide a large set of data from which all the relevant information about an isolate can be mined. However, rigorous validation and careful consideration of potential implications in the replacement of traditional surveillance methodologies with WGS data analysis tools is needed. Two in silico tools for Salmonella serotyping have been developed, the Salmonella in silico Typing Resource (SISTR) and SeqSero, while seven gene MLST for serovar prediction can be adapted for in silico analysis. All three analysis methods were assessed and compared to traditional serotyping techniques using a set of 813 verified clinical and laboratory isolates, including 492 Canadian clinical isolates and 321 isolates of human and non-human sources. Successful results were obtained for 94.8, 88.2, and 88.3% of the isolates tested using SISTR, SeqSero, and MLST, respectively, indicating all would be suitable for maintaining historical records, surveillance systems, and communication structures currently in place and the choice of the platform used will ultimately depend on the users need. Results also pointed to the need to reframe serotyping in the genomic era as a test to understand the genes that are carried by an isolate, one which is not necessarily congruent with what is antigenically expressed. The adoption of WGS for serotyping will provide the simultaneous collection of information that can be used by multiple programs within the current surveillance paradigm; however, this does not negate the importance of the various programs or the role of serotyping going forward.

Food Safety in the Age of Next Generation Sequencing, Bioinformatics, and Open Data Access

Public health labs and food regulatory agencies globally are embracing whole genome sequencing (WGS) as a revolutionary new method that is positioned to replace numerous existing diagnostic and microbial typing technologies with a single new target: the microbial draft genome. The ability to cheaply generate large amounts of microbial genome sequence data, combined with emerging policies of food regulatory and public health institutions making their microbial sequences increasingly available and public, has served to open up the field to the general scientific community. This open data access policy shift has resulted in a proliferation of data being deposited into sequence repositories and of novel bioinformatics software designed to analyze these vast datasets. There also has been a more recent drive for improved data sharing to achieve more effective global surveillance, public health and food safety. Such developments have heightened the need for enhanced analytical systems in order to process and interpret this new type of data in a timely fashion. In this review we outline the emergence of genomics, bioinformatics and open data in the context of food safety. We also survey major efforts to translate genomics and bioinformatics technologies out of the research lab and into routine use in modern food safety labs. We conclude by discussing the challenges and opportunities that remain, including those expected to play a major role in the future of food safety science.

Epidemiology of antimicrobial resistant Campylobacter spp. isolated from retail meats in Canada

Campylobacter is an important zoonotic pathogen found in livestock and can cause illness in humans following consumption of raw and undercooked meat products. The objectives of this study were to determine the prevalence of Campylobacter spp. in retail meat (poultry, turkey, pork and beef) purchased in Alberta, Canada and to assess antimicrobial resistance and genetic relatedness of recovered Campylobacter strains with previously isolated strains from clinical and environmental sources. A Comparative Genomic Fingerprinting (CGF) method was used for assessing genetic relatedness of isolates. A total of 606 samples comprising 204, 110, 145 and 147 samples of retail chicken, turkey, ground beef and pork, respectively, were obtained. Campylobacter was isolated from 23.5% (48/204) of chicken samples and 14.2% (8/110) of turkey samples. Pork and beef samples were negative for Campylobacter. Campylobacter jejuni was the most common (94.6%) spp. found followed by C. coli (5.4%). Resistance to tetracycline was found in 48.1% of isolates, followed by resistance to ciprofloxacin (5.5%), nalidixic acid (5.5%), azithromycin (1.78%), and erythromycin (1.78%). All isolates were susceptible to clindamycin, florfenicol, gentamicin and telithromycin. Tetracycline resistance was attributable to the presence of the tetO gene. CGF analysis showed that Campylobacter isolated from poultry meat in this study were genetically related to clinical isolates recovered from human infections and to those isolated from animals and the environment.

Prevalence and diversity of waterborne Arcobacter butzleri in southwestern Alberta, Canada

Arcobacter butzleri is a potential enteric pathogen to human beings, but its reservoirs and modes of transmission are largely unverified. Microbiological and molecular detection and subtyping techniques can facilitate surveillance of A. butzleri in hosts and environmental reservoirs. We isolated A. butzleri from 173 surface water samples (25.6%) and 81 treated wastewater samples (77.9%) collected in southwestern Alberta over a 1-year period. Arcobacter butzleri isolates (n = 500) were genotyped and compared to determine diversity of A. butzleri in southwestern Alberta. Culture methods affected the frequency of detection and genotype diversity of A. butzleri, and isolation comprehensiveness was different for surface waters and treated wastewaters. Detection of A. butzleri in the Oldman River Watershed corresponded with season, river flow rates, and fecal coliform densities. Arcobacter butzleri was detected most frequently in treated wastewater, in the Oldman River downstream from treated wastewater outfalls, and in tributaries near areas of intensive confined feeding operations. All sample sources possessed high genotype diversity, and A. butzleri isolates from treated wastewaters were genetically similar to isolates from the Oldman River downriver from treated wastewater outfall sites. In southwestern Alberta, municipal and agricultural activities contribute to the density and genotype diversity of A. butzleri in surface waters.

Genomic insights from whole genome sequencing of four clonal outbreak Campylobacter jejuni assessed within the global C. jejuni population

BACKGROUND

Whole genome sequencing (WGS) is useful for determining clusters of human cases, investigating outbreaks, and defining the population genetics of bacteria. It also provides information about other aspects of bacterial biology, including classical typing results, virulence, and adaptive strategies of the organism. Cell culture invasion and protein expression patterns of four related multilocus sequence type 21 (ST21) C. jejuni isolates from a significant Canadian water-borne outbreak were previously associated with the presence of a CJIE1 prophage. Whole genome sequencing was used to examine the genetic diversity among these isolates and confirm that previous observations could be attributed to differential prophage carriage. Moreover, we sought to determine the presence of genome sequences that could be used as surrogate markers to delineate outbreak-associated isolates.

RESULTS

Differential carriage of the CJIE1 prophage was identified as the major genetic difference among the four outbreak isolates. High quality single-nucleotide variant (hqSNV) and core genome multilocus sequence typing (cgMLST) clustered these isolates within expanded datasets consisting of additional C. jejuni strains. The number and location of homopolymeric tract regions was identical in all four outbreak isolates but differed from all other C. jejuni examined. Comparative genomics and PCR amplification enabled the identification of large chromosomal inversions of approximately 93 kb and 388 kb within the outbreak isolates associated with transducer-like proteins containing long nucleotide repeat sequences. The 93-kb inversion was characteristic of the outbreak-associated isolates, and the gene content of this inverted region displayed high synteny with the reference strain.

CONCLUSIONS

The four outbreak isolates were clonally derived and differed mainly in the presence of the CJIE1 prophage, validating earlier findings linking the prophage to phenotypic differences in virulence assays and protein expression. The identification of large, genetically syntenous chromosomal inversions in the genomes of outbreak-associated isolates provided a unique method for discriminating outbreak isolates from the background population. Transducer-like proteins appear to be associated with the chromosomal inversions. CgMLST and hqSNV analysis also effectively delineated the outbreak isolates within the larger C. jejuni population structure.

Efficacy of wastewater treatment on Arcobacter butzleri density and strain diversity

Arcobacter butzleri is a suspected waterborne enteric pathogen that is ubiquitous in the environment, but the degree to which wastewater treatment prevents entry of A. butzleri into environmental waters and the risks posed are not well established. Untreated and treated wastewater samples (n = 260) were collected weekly from the Lethbridge and Fort Macleod wastewater treatment facilities (the two major municipal inputs in southwestern Alberta, Canada) from May 2008 to April 2009. Untreated wastewaters contained high densities of A. butzleri and fecal coliform indicators, and densities at Lethbridge were typically higher than at Fort Macleod. Data indicated that A. butzleri and fecal coliform densities in wastewater were greatest in autumn and lowest in winter. Mechanical and biological treatment of wastewaters reduced but did not eliminate fecal coliform indicators or A. butzleri. At Lethbridge, UVB irradiation of mechanically and biologically treated wastewater further reduced densities of fecal coliform indicators. There was high A. butzleri genotype diversity in all sample sources, and survival during treatment was not strain-dependent. No genotype was dominant in any sample source, but 8.9% of genotypes were recurrent over time, and 4.4% of genotypes were detected at both wastewater treatment facilities. The current study demonstrates that viable A. butzleri are able to survive wastewater treatment, including UVB irradiation, which may lead to increased density and genetic diversity of this suspected pathogen in environmental waters via wastewater effluent discharge.

Cluster Analysis of Campylobacter jejuni Genotypes Isolated from Small and Medium-Sized Mammalian Wildlife and Bovine Livestock from Ontario Farms

Using data collected from a cross-sectional study of 25 farms (eight beef, eight swine and nine dairy) in 2010, we assessed clustering of molecular subtypes of C. jejuni based on a Campylobacter-specific 40 gene comparative genomic fingerprinting assay (CGF40) subtypes, using unweighted pair-group method with arithmetic mean (UPGMA) analysis, and multiple correspondence analysis. Exact logistic regression was used to determine which genes differentiate wildlife and livestock subtypes in our study population. A total of 33 bovine livestock (17 beef and 16 dairy), 26 wildlife (20 raccoon (Procyon lotor), five skunk (Mephitis mephitis) and one mouse (Peromyscus spp.) C. jejuni isolates were subtyped using CGF40. Dendrogram analysis, based on UPGMA, showed distinct branches separating bovine livestock and mammalian wildlife isolates. Furthermore, two-dimensional multiple correspondence analysis was highly concordant with dendrogram analysis showing clear differentiation between livestock and wildlife CGF40 subtypes. Based on multilevel logistic regression models with a random intercept for farm of origin, we found that isolates in general, and raccoons more specifically, were significantly more likely to be part of the wildlife branch. Exact logistic regression conducted gene by gene revealed 15 genes that were predictive of whether an isolate was of wildlife or bovine livestock isolate origin. Both multiple correspondence analysis and exact logistic regression revealed that in most cases, the presence of a particular gene (13 of 15) was associated with an isolate being of livestock rather than wildlife origin. In conclusion, the evidence gained from dendrogram analysis, multiple correspondence analysis and exact logistic regression indicates that mammalian wildlife carry CGF40 subtypes of C. jejuni distinct from those carried by bovine livestock. Future studies focused on source attribution of C. jejuni in human infections will help determine whether wildlife transmit Campylobacter jejuni directly to humans.

Molecular and Statistical Analysis of Campylobacter spp. and Antimicrobial-Resistant Campylobacter Carriage in Wildlife and Livestock from Ontario Farms

The objectives of this study were to (i) compare the carriage of Campylobacter and antimicrobial-resistant Campylobacter among livestock and mammalian wildlife on Ontario farms, and (ii) investigate the potential sharing of Campylobacter subtypes between livestock and wildlife. Using data collected from a cross-sectional study of 25 farms in 2010, we assessed associations, using mixed logistic regression models, between Campylobacter and antimicrobial-resistant Campylobacter carriage and the following explanatory variables: animal species (beef, dairy, swine, raccoon, other), farm type (swine, beef, dairy), type of sample (livestock or wildlife) and Campylobacter species (jejuni, coli, other). Models included a random effect to account for clustering by farm where samples were collected. Samples were subtyped using a Campylobacter-specific 40 gene comparative fingerprinting assay. A total of 92 livestock and 107 wildlife faecal samples were collected, and 72% and 27% tested positive for Campylobacter, respectively. Pooled faecal samples from livestock were significantly more likely to test positive for Campylobacter than wildlife samples. Relative to dairy cattle, pig samples were at significantly increased odds of testing positive for Campylobacter. The odds of isolating Campylobacter jejuni from beef cattle samples were significantly greater compared to dairy cattle and raccoon samples. Fifty unique subtypes of Campylobacter were identified, and only one subtype was found in both wildlife and livestock samples. Livestock Campylobacter isolates were significantly more likely to exhibit antimicrobial resistance (AMR) compared to wildlife Campylobacter isolates. Campylobacter jejuni was more likely to exhibit AMR when compared to C. coli. However, C. jejuni isolates were only resistant to tetracycline, and C.  coli isolates exhibited multidrug resistance patterns. Based on differences in prevalence of Campylobacter spp. and resistant Campylobacter between livestock and wildlife samples, and the lack of similarity in molecular subtypes and AMR patterns, we concluded that the sharing of Campylobacter species between livestock and mammalian wildlife was uncommon.

SuperPhy: predictive genomics for the bacterial pathogen Escherichia coli

BACKGROUND

Predictive genomics is the translation of raw genome sequence data into a phenotypic assessment of the organism. For bacterial pathogens, these phenotypes can range from environmental survivability, to the severity of human disease. Significant progress has been made in the development of generic tools for genomic analyses that are broadly applicable to all microorganisms; however, a fundamental missing component is the ability to analyze genomic data in the context of organism-specific phenotypic knowledge, which has been accumulated from decades of research and can provide a meaningful interpretation of genome sequence data.

RESULTS

In this study, we present SuperPhy, an online predictive genomics platform ( http://lfz.corefacility.ca/superphy/ ) for Escherichia coli. The platform integrates the analytical tools and genome sequence data for all publicly available E. coli genomes and facilitates the upload of new genome sequences from users under public or private settings. SuperPhy provides real-time analyses of thousands of genome sequences with results that are understandable and useful to a wide community, including those in the fields of clinical medicine, epidemiology, ecology, and evolution. SuperPhy includes identification of: 1) virulence and antimicrobial resistance determinants 2) statistical associations between genotypes, biomarkers, geospatial distribution, host, source, and phylogenetic clade; 3) the identification of biomarkers for groups of genomes on the based presence/absence of specific genomic regions and single-nucleotide polymorphisms and 4) in silico Shiga-toxin subtype.

CONCLUSIONS

SuperPhy is a predictive genomics platform that attempts to provide an essential link between the vast amounts of genome information currently being generated and phenotypic knowledge in an organism-specific context.

The Salmonella In Silico Typing Resource (SISTR): An Open Web-Accessible Tool for Rapidly Typing and Subtyping Draft Salmonella Genome Assemblies

For nearly 100 years serotyping has been the gold standard for the identification of Salmonella serovars. Despite the increasing adoption of DNA-based subtyping approaches, serotype information remains a cornerstone in food safety and public health activities aimed at reducing the burden of salmonellosis. At the same time, recent advances in whole-genome sequencing (WGS) promise to revolutionize our ability to perform advanced pathogen characterization in support of improved source attribution and outbreak analysis. We present the Salmonella In Silico Typing Resource (SISTR), a bioinformatics platform for rapidly performing simultaneous in silico analyses for several leading subtyping methods on draft Salmonella genome assemblies. In addition to performing serovar prediction by genoserotyping, this resource integrates sequence-based typing analyses for: Multi-Locus Sequence Typing (MLST), ribosomal MLST (rMLST), and core genome MLST (cgMLST). We show how phylogenetic context from cgMLST analysis can supplement the genoserotyping analysis and increase the accuracy of in silico serovar prediction to over 94.6% on a dataset comprised of 4,188 finished genomes and WGS draft assemblies. In addition to allowing analysis of user-uploaded whole-genome assemblies, the SISTR platform incorporates a database comprising over 4,000 publicly available genomes, allowing users to place their isolates in a broader phylogenetic and epidemiological context. The resource incorporates several metadata driven visualizations to examine the phylogenetic, geospatial and temporal distribution of genome-sequenced isolates. As sequencing of Salmonella isolates at public health laboratories around the world becomes increasingly common, rapid in silico analysis of minimally processed draft genome assemblies provides a powerful approach for molecular epidemiology in support of public health investigations. Moreover, this type of integrated analysis using multiple sequence-based methods of sub-typing allows for continuity with historical serotyping data as we transition towards the increasing adoption of genomic analyses in epidemiology. The SISTR platform is freely available on the web at https://lfz.corefacility.ca/sistr-app/.

Extensive characterization of Campylobacter jejuni chicken isolates to uncover genes involved in the ability to compete for gut colonization

BACKGROUND

Campylobacter jejuni is responsible for human foodborne enteritis. This bacterium is a remarkable colonizer of the chicken gut, with some strains outcompeting others for colonization. To better understand this phenomenon, the objective of this study was to extensively characterize the phenotypic performance of C. jejuni chicken strains and associate their gut colonizing ability with specific genes.

RESULTS

C. jejuni isolates (n = 45) previously analyzed for the presence of chicken colonization associated genes were further characterized for phenotypic properties influencing colonization: autoagglutination and chemotaxis as well as adhesion to and invasion of primary chicken caecal cells. This allowed strains to be ranked according to their in vitro performance. After their in vitro capacity to outcompete was demonstrated in vivo, strains were then typed by comparative genomic fingerprinting (CGF). In vitro phenotypical properties displayed a linear variability among the tested strains. Strains possessing higher scores for phenotypical properties were able to outcompete others during chicken colonization trials. When the gene content of strains was compared, some were associated with different phenotypical scores and thus with different outcompeting capacities. Use of CGF profiles showed an extensive genetic variability among the studied strains and suggested that the outcompeting capacity is not predictable by CGF profile.

CONCLUSION

This study revealed a wide array of phenotypes present in C. jejuni strains, even though they were all recovered from chicken caecum. Each strain was classified according to its in vitro competitive potential and its capacity to compete for chicken gut colonization was associated with specific genes. This study also exposed the disparity existing between genetic typing and phenotypical behavior of C. jejuni strains.

Development of a comparative genomic fingerprinting assay for rapid and high resolution genotyping of Arcobacter butzleri

Background

Molecular typing methods are critical for epidemiological investigations, facilitating disease outbreak detection and source identification. Study of the epidemiology of the emerging human pathogen Arcobacter butzleri is currently hampered by the lack of a subtyping method that is easily deployable in the context of routine epidemiological surveillance. In this study we describe a comparative genomic fingerprinting (CGF) method for high-resolution and high-throughput subtyping of A. butzleri. Comparative analysis of the genome sequences of eleven A. butzleri strains, including eight strains newly sequenced as part of this project, was employed to identify accessory genes suitable for generating unique genetic fingerprints for high-resolution subtyping based on gene presence or absence within a strain.

Results

A set of eighty-three accessory genes was used to examine the population structure of a dataset comprised of isolates from various sources, including human and non-human animals, sewage, and river water (n=156). A streamlined assay (CGF₄₀) based on a subset of 40 genes was subsequently developed through marker optimization. High levels of profile diversity (121 distinct profiles) were observed among the 156 isolates in the dataset, and a high Simpson’s Index of Diversity (ID) observed (ID > 0.969) indicate that the CGF₄₀ assay possesses high discriminatory power. At the same time, our observation that 115 isolates in this dataset could be assigned to 29 clades with a profile similarity of 90% or greater indicates that the method can be used to identify clades comprised of genetically similar isolates.

Conclusions

The CGF₄₀ assay described herein combines high resolution and repeatability with high throughput for the rapid characterization of A. butzleri strains. This assay will facilitate the study of the population structure and epidemiology of A. butzleri.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0426-4) contains supplementary material, which is available to authorized users.

Comparative variation within the genome of Campylobacter jejuni NCTC 11168 in human and murine hosts

Campylobacteriosis incited by C. jejuni is a significant enteric disease of human beings. A person working with two reference strains of C. jejuni National Collection of Type Cultures (NCTC) 11168 developed symptoms of severe enteritis including bloody diarrhea. The worker was determined to be infected by C. jejuni. In excess of 50 isolates were recovered from the worker's stool. All of the recovered isolates and the two reference strains were indistinguishable from each other based on comparative genomic fingerprint subtyping. Whole genome sequence analysis indicated that the worker was infected with a C. jejuni NCTC 11168 obtained from the American Type Culture Collection; this strain (NCTC 11168-GSv) is the genome sequence reference. After passage through the human host, major genetic changes including indel mutations within twelve contingency loci conferring phase variations were detected in the genome of C. jejuni. Specific and robust single nucleotide polymorphism (SNP) changes in the human host were also observed in two loci (Cj0144c, Cj1564). In mice inoculated with an isolate of C. jejuni NCTC 11168-GSv from the infected person, the isolate underwent further genetic variation. At nine loci, mutations specific to inoculated mice including five SNP changes were observed. The two predominant SNPs observed in the human host reverted in mice. Genetic variations occurring in the genome of C. jejuni in mice corresponded to increased densities of C. jejuni cells associated with cecal mucosa. In conclusion, C. jejuni NCTC 11168-GSv was found to be highly virulent in a human being inciting severe enteritis. Host-specific mutations in the person with enteritis occurred/were selected for in the genome of C. jejuni, and many were not maintained in mice. Information obtained in the current study provides new information on host-specific genetic adaptation by C. jejuni.

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.

Evaluation of MALDI-TOF mass spectroscopy methods for determination of Escherichia coli pathotypes

It is rapidly becoming apparent that many E. coli pathotypes cause a considerable burden of human disease. Surveillance of these organisms is difficult because there are few or no simple, rapid methods for detecting and differentiating the different pathotypes. MALDI-TOF mass spectroscopy has recently been rapidly and enthusiastically adopted by many clinical laboratories as a diagnostic method because of its high throughput, relatively low cost, and adaptability to the laboratory workflow. To determine whether the method could be adapted for E. coli pathotype differentiation the Bruker Biotyper methodology and a second methodology adapted from the scientific literature were tested on isolates representing eight distinct pathotypes and two other groups of E. coli. A total of 136 isolates was used for this study. Results confirmed that the Bruker Biotyper methodology that included extraction of proteins from bacterial cells was capable of identifying E. coli isolates from all pathotypes to the species level and, furthermore, that the Bruker extraction and MALDI-TOF MS with the evaluation criteria developed in this work was effective for differentiating most pathotypes.

Detection of genetic diversity in Campylobacter jejuni isolated from a commercial turkey flock using flaA typing, MLST analysis and microarray assay

Campylobacter is genetically highly diverse and undergoes frequent intraspecific recombination. Turkeys have been identified as an important reservoir for Campylobacter jejuni which is of public health significance. The assessment of the genetic diversity among Campylobacter population is critical for our understanding of the epidemiology of this bacterium. The genetic profiles were different according to the molecular typing methods used. The performance of established flaA genotyping, multilocus sequencing typing (MLST) and DNA microarray assay based on the ArrayTube™ technology was evaluated using 14 Campylobacter jejuni isolated from a commercial turkey flock. The flaA typing was performed using PCR-RFLP with restriction enzymes Sau3AI, AluI, a ‘composite’ flaA analysis of AluI and Sau3AI and DdeI. The 14 isolates were differentiated into 3, 5, 7 and 9 genotypes, respectively. Entire flaA gene and short variable region (SVR) sequences were analysed. Sequencing of the entire flaA provided 11 different genotypes. flaA-SVR sequence analysis detected 8 flaA alleles and 4 flaA peptides. One new flaA allele type (528) was identified. MLST analysis represented 10 different sequence types (STs) and 5 clonal complexes (CCs). The microarray assay recognised 14 different genotypes. The discriminatory indices were 0.560, 0.802, 0.857, and 0.912 for flaA-RFLP depending on the used enzymes, 0.890 for flaA-SVR, 0.967 for entire flaA sequencing, 0.945 for MLST and 1.00 for the DNA microarray assay. The flaA gene was genetically stable over 20 passages on blood agar. In conclusion, the different typing tools demonstrated a high level of genetic heterogeneity of Campylobacter jejuni in a turkey flock, indicating that a single flock can be infected by multiple genotypes within one rearing cycle. DNA microarray-based assays had the highest discriminatory power when compared with other genotyping tools.

A framework for assessing the concordance of molecular typing methods and the true strain phylogeny of Campylobacter jejuni and C. coli using draft genome sequence data

Tracking of sources of sporadic cases of campylobacteriosis remains challenging, as commonly used molecular typing methods have limited ability to unambiguously link genetically related strains. Genomics has become increasingly prominent in the public health response to enteric pathogens as methods enable characterization of pathogens at an unprecedented level of resolution. However, the cost of sequencing and expertise required for bioinformatic analyses remains prohibitive, and these comprehensive analyses are limited to a few priority strains. Although several molecular typing methods are currently widely used for epidemiological analysis of campylobacters, it is not clear how accurately these methods reflect true strain relationships. To address this, we have developed a framework and associated computational tools to rapidly analyze draft genome sequence data for the assessment of molecular typing methods against a "gold standard" based on the phylogenetic analysis of highly conserved core (HCC) genes with high sequence quality. We analyzed 104 publicly available whole genome sequences (WGS) of C. jejuni and C. coli. In addition to in silico determination of multi-locus sequence typing (MLST), flaA, and porA type, as well as comparative genomic fingerprinting (CGF) type, we inferred a "reference" phylogeny based on 389 HCC genes. Molecular typing data were compared to the reference phylogeny for concordance using the adjusted Wallace coefficient (AWC) with confidence intervals. Although MLST targets the sequence variability in core genes and CGF targets insertions/deletions of accessory genes, both methods are based on multi-locus analysis and provided better estimates of true phylogeny than methods based on single loci (porA, flaA). A more comprehensive WGS dataset including additional genetically related strains, both epidemiologically linked and unlinked, will be necessary to more comprehensively assess the performance of subtyping methods for outbreak investigations and surveillance activities. Analyses of the strengths and weaknesses of widely used typing methodologies in inferring true strain relationships will provide guidance in the interpretation of this data for epidemiological purposes.

Spatial and temporal drivers of zoonotic pathogen contamination of an agricultural watershed

In regions where animal agriculture is prominent, such as southern Alberta, higher rates of gastrointestinal illness have been reported when compared with nonagricultural regions. This difference in the rate of illness is thought to be a result of increased zoonotic pathogen exposure through environmental sources such as water. In this study, temporal and spatial factors associated with bacterial pathogen contamination of the Oldman River, which transverses this region, were analyzed using classification and regression tree analysis. Significantly higher levels of fecal indicators; more frequent isolations of Campylobacter spp., Escherichia coli O157:H7, and Salmonella enterica spp.; and higher rates of detection of pig-specific Bacteroides markers occurred at downstream sites than at upstream sites, suggesting additive stream inputs. Fecal indicator densities were also significantly higher when any one of these three bacterial pathogens was present and where there were higher total animal manure units; however, occasionally pathogens were present when fecal indicator levels were low or undetectable. Overall, Salmonella spp., Campylobacter spp., and E. coli O157:H7 presence was associated with season, animal manure units, and total rainfall on the day of sampling and 3 d in advance of sampling. Several of the environmental variables analyzed in this study appear to influence pathogen prevalence and therefore may be useful in predicting water quality and safety and in the improvement of watershed management practices in this and other agricultural regions.

Identification of Salmonella enterica species- and subgroup-specific genomic regions using Panseq 2.0

The pan-genome of a taxonomic group consists of evolutionarily conserved core genes shared by all members and accessory genes that are present only in some members of the group. Group- and subgroup-specific core genes are thought to contribute to shared phenotypes such as virulence and niche specificity. In this study we analyzed 39 Salmonella enterica genomes (16 closed, 23 draft), a species that contains two human-specific serovars that cause typhoid fever, as well as a large number of zoonotic serovars that cause gastroenteritis in humans. Panseq 2.0 was used to define the pan-genome by adjusting the threshold at which group-specific "core" loci are defined. We found the pan-genome to be 9.03 Mbp in size, and that the core genome size decreased, while the number of SNPs/100 bp increased, as the number of strains used to define the core genome increased, suggesting substantial divergence among S. enterica subgroups. Subgroup-specific "core" genes, in contrast, had fewer SNPs/100 bp, likely reflecting their more recent acquisition. Phylogenetic trees were created from the concatenated and aligned pan-genome, the core genome, and multi-locus-sequence typing (MLST) loci. Branch support increased among the trees, and strains of the same serovar grouped closer together as the number of loci used to create the tree increased. Further, high levels of discrimination were achieved even amongst the most closely related strains of S. enterica Typhi, suggesting that the data generated by Panseq may also be of value in short-term epidemiological studies. Panseq provides an easy and fast way of performing pan-genomic analyses, which can include the identification of group-dominant as well as group-specific loci and is available as a web-server and a standalone version at http://lfz.corefacility.ca/panseq/.

In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence

Background

Many approaches have been used to study the evolution, population structure and genetic diversity of Escherichia coli O157:H7; however, observations made with different genotyping systems are not easily relatable to each other. Three genetic lineages of E. coli O157:H7 designated I, II and I/II have been identified using octamer-based genome scanning and microarray comparative genomic hybridization (mCGH). Each lineage contains significant phenotypic differences, with lineage I strains being the most commonly associated with human infections. Similarly, a clade of hyper-virulent O157:H7 strains implicated in the 2006 spinach and lettuce outbreaks has been defined using single-nucleotide polymorphism (SNP) typing. In this study an in silico comparison of six different genotyping approaches was performed on 19 E. coli genome sequences from 17 O157:H7 strains and single O145:NM and K12 MG1655 strains to provide an overall picture of diversity of the E. coli O157:H7 population, and to compare genotyping methods for O157:H7 strains.

Results

In silico determination of lineage, Shiga-toxin bacteriophage integration site, comparative genomic fingerprint, mCGH profile, novel region distribution profile, SNP type and multi-locus variable number tandem repeat analysis type was performed and a supernetwork based on the combination of these methods was produced. This supernetwork showed three distinct clusters of strains that were O157:H7 lineage-specific, with the SNP-based hyper-virulent clade 8 synonymous with O157:H7 lineage I/II. Lineage I/II/clade 8 strains clustered closest on the supernetwork to E. coli K12 and E. coli O55:H7, O145:NM and sorbitol-fermenting O157 strains.

Conclusion

The results of this study highlight the similarities in relationships derived from multi-locus genome sampling methods and suggest a "common genotyping language" may be devised for population genetics and epidemiological studies. Future genotyping methods should provide data that can be stored centrally and accessed locally in an easily transferable, informative and extensible format based on comparative genomic analyses.

Genomics-based molecular epidemiology of Campylobacter jejuni isolates from feedlot cattle and from people in Alberta, Canada

Feedlot cattle in Alberta, Canada, have been identified as reservoirs for Campylobacter jejuni, an important human pathogen. Oligonucleotide DNA microarrays were used as a platform to compare C. jejuni isolates from feedlot cattle and human clinical cases from Alberta. Comparative genomic hybridization (CGH) analysis was performed on 87 isolates (46 bovine, 41 human) obtained within the same geographical regions and time frame. Thirteen CGH clusters were obtained based on overall comparative genomic profile similarity. Nine CGH clusters contained human and cattle isolates, three contained only human isolates, and one contained only cattle isolates. The study isolates clustered regardless of temporal or geographical frameworks. In addition, array genes (n = 1,399) were investigated on a gene-by-gene basis to see if any were unequally distributed between human and cattle sources or between clusters dominated by either human or cattle isolates ("human enriched" versus "cattle enriched"). Using Fisher's exact test with the Westfall and Young correction for these comparisons, a small number of differentially distributed genes were identified. Our findings suggest that feedlot cattle and human C. jejuni strains are very similar and may be endemic within Alberta. Further, the common distribution of human clinical and bovine C. jejuni isolates within the same genetically based clusters suggests that dynamic and important transmission routes between cattle and human populations may exist.

Comparative genomic assessment of Multi-Locus Sequence Typing: rapid accumulation of genomic heterogeneity among clonal isolates of Campylobacter jejuni

Background

Multi-Locus Sequence Typing (MLST) has emerged as a leading molecular typing method owing to its high ability to discriminate among bacterial isolates, the relative ease with which data acquisition and analysis can be standardized, and the high portability of the resulting sequence data. While MLST has been successfully applied to the study of the population structure for a number of different bacterial species, it has also provided compelling evidence for high rates of recombination in some species. We have analyzed a set of Campylobacter jejuni strains using MLST and Comparative Genomic Hybridization (CGH) on a full-genome microarray in order to determine whether recombination and high levels of genomic mosaicism adversely affect the inference of strain relationships based on the analysis of a restricted number of genetic loci.

Results

Our results indicate that, in general, there is significant concordance between strain relationships established by MLST and those based on shared gene content as established by CGH. While MLST has significant predictive power with respect to overall genome similarity of isolates, we also found evidence for significant differences in genomic content among strains that would otherwise appear to be highly related based on their MLST profiles.

Conclusion

The extensive genomic mosaicism between closely related strains has important implications in the context of establishing strain to strain relationships because it suggests that the exact gene content of strains, and by extension their phenotype, is less likely to be "predicted" based on a small number of typing loci. This in turn suggests that a greater emphasis should be placed on analyzing genes of clinical interest as we forge ahead with the next generation of molecular typing methods.

Comparative genomic analysis of Campylobacter jejuni associated with Guillain-Barré and Miller Fisher syndromes: neuropathogenic and enteritis-associated isolates can share high levels of genomic similarity

Background

Campylobacter jejuni infection represents the most frequent antecedent infection triggering the onset of the neuropathic disorders Guillain-Barré syndrome (GBS) and Miller Fisher syndrome (MFS). Although sialylated ganglioside-mimicking lipo-oligosaccharide (LOS) structures are the strongest neuropathogenic determinants in C. jejuni, they do not appear to be the only requirement for a neuropathic outcome since strains capable of their production have been isolated from patients with uncomplicated cases of enteritis. Consequently, other pathogen and/or host-related factors contribute to the onset of neurological complications. We have used comparative genomic hybridization to perform a detailed genomic comparison of strains isolated from GBS/MFS and enteritis-only patients. Our dataset, in which the gene conservation profile for 1712 genes was assayed in 102 strains, including 56 neuropathogenic isolates, represents the largest systematic search for C. jejuni factors associated with GBS/MFS to date and has allowed us to analyze the genetic background of neuropathogenic C. jejuni strains with an unprecedented level of resolution.

Results

The majority of GBS/MFS strains can be assigned to one of six major lineages, suggesting that several genetic backgrounds can result in a neuropathogenic phenotype. A statistical analysis of gene conservation rates revealed that although genes involved in the sialylation of LOS structures were significantly associated with neuropathogenic strains, still many enteritis-control strains both bear these genes and share remarkable levels of genomic similarity with their neuropathogenic counterparts. Two capsule biosynthesis genes (Cj1421c and Cj1428c) showed higher conservation rates among neuropathogenic strains compared to enteritis-control strains. Any potential involvement of these genes in neuropathogenesis must be assessed. A single gene (HS:3 Cj1135) had a higher conservation rate among enteritis-control strains. This gene encodes a glucosyltransferase that is found in some of the LOS classes that do not express ganglioside mimics.

Conclusion

Our findings corroborate that neuropathogenic factors may be transferred between unrelated strains of different genetic background. Our results would also suggest that the failure of some strains isolated from uncomplicated cases of enteritis to elicit a neuropathic clinical outcome may be due to subtle genetic differences that silence their neuropathogenic potential and/or due to host-related factors.

The microarray data has been deposited in NCBI's Gene Expression Omnibus under accession number GSE3579.

Studying bacterial genome dynamics using microarray-based comparative genomic hybridization

Genome sequencing has revealed the remarkable amount of genetic diversity that can be encountered in bacterial genomes. In particular, the comparison of genome sequences from closely related strains has uncovered significant differences in gene content, hinting at the dynamic nature of bacterial genomes. The study of these genome dynamics is crucial to leveraging genomic information because the genome sequence of a single bacterial strain may not accurately represent the genome of the species. The dynamic nature of bacterial genome content has required us to apply the concepts of comparative genomics (CG) at the species level. Although direct genome sequence comparisons are an ideal method of performing CG, one current constraint is the limited availability of multiple genome sequences from a given bacterial species. DNA microarray-based comparative genomic hybridization (MCGH), which can be used to determine the presence or absence of thousands of genes in a single hybridization experiment, provides a powerful alternative for determining genome content and has been successfully used to investigate the genome dynamics of a wide number of bacterial species. Although MCGH-based studies have already provided a new vista on bacterial genome diversity, original methods for MCGH have been limited by the absence of novel gene sequences included in the microarray. New applications of the MCGH platform not only promise to accelerate the pace of novel gene discovery but will also help provide an integrated microarray-based approach to the study of bacterial CG.

Comparative genomics profiling of clinical isolates of Aeromonas salmonicida using DNA microarrays

BACKGROUND

Aeromonas salmonicida has been isolated from numerous fish species and shows wide variation in virulence and pathogenicity. As part of a larger research program to identify virulence genes and candidates for vaccine development, a DNA microarray was constructed using a subset of 2024 genes from the draft genome sequence of A. salmonicida subsp. salmonicida strain A449. The microarray included genes encoding known virulence-associated factors in A. salmonicida and homologs of virulence genes of other pathogens. We used microarray-based comparative genomic hybridizations (M-CGH) to compare selected A. salmonicida sub-species and other Aeromonas species from different hosts and geographic locations.

RESULTS

Results showed variable carriage of virulence-associated genes and generally increased variation in gene content across sub-species and species boundaries. The greatest variation was observed among genes associated with plasmids and transposons. There was little correlation between geographic region and degree of variation for all isolates tested.

CONCLUSION

We have used the M-CGH technique to identify subsets of conserved genes from amongst this set of A. salmonicida virulence genes for further investigation as potential vaccine candidates. Unlike other bacterial characterization methods that use a small number of gene or DNA-based functions, M-CGH examines thousands of genes and/or whole genomes and thus is a more comprehensive analytical tool for veterinary or even human health research.

A new approach for the analysis of bacterial microarray-based Comparative Genomic Hybridization: insights from an empirical study

BACKGROUND

Microarray-based Comparative Genomic Hybridization (M-CGH) has been used to characterize the extensive intraspecies genetic diversity found in bacteria at the whole-genome level. Although conventional microarray analytical procedures have proved adequate in handling M-CGH data, data interpretation using these methods is based on a continuous character model in which gene divergence and gene absence form a spectrum of decreasing gene conservation levels. However, whereas gene divergence may yet be accompanied by retention in gene function, gene absence invariably leads to loss of function. This distinction, if ignored, leads to a loss in the information to be gained from M-CGH data. We present here results from experiments in which two genome-sequenced strains of C. jejuni were compared against each other using M-CGH. Because the gene content of both strains was known a priori, we were able to closely examine the effects of sequence divergence and gene absence on M-CGH data in order to define analytical parameters for M-CGH data interpretation. This would facilitate the examination of the relative effects of sequence divergence or gene absence in comparative genomics analyses of multiple strains of any species for which genome sequence data and a DNA microarray are available.

RESULTS

As a first step towards improving the analysis of M-CGH data, we estimated the degree of experimental error in a series of experiments in which identical samples were compared against each other by M-CGH. This variance estimate was used to validate a Log Ratio-based methodology for identification of outliers in M-CGH data. We compared two genome strains by M-CGH to examine the effect of probe/target identity on the Log Ratios of signal intensities using prior knowledge of gene divergence and gene absence to establish Log Ratio thresholds for the identification of absent and conserved genes.

CONCLUSION

The results from this empirical study validate the Log Ratio thresholds that have been used in other studies to establish gene divergence/absence. Moreover, the analytical framework presented here enhances the information content derived from M-CGH data by shifting the focus from divergent/absent gene detection to accurate detection of conserved and absent genes. This approach closely aligns the technical limitations of M-CGH analysis with practical limitations on the biological interpretation of comparative genomics data.

Large-scale comparative genomics meta-analysis of Campylobacter jejuni isolates reveals low level of genome plasticity

We have used comparative genomic hybridization (CGH) on a full-genome Campylobacter jejuni microarray to examine genome-wide gene conservation patterns among 51 strains isolated from food and clinical sources. These data have been integrated with data from three previous C. jejuni CGH studies to perform a meta-analysis that included 97 strains from the four separate data sets. Although many genes were found to be divergent across multiple strains (n = 350), many genes (n = 249) were uniquely variable in single strains. Thus, the strains in each data set comprise strains with a unique genetic diversity not found in the strains in the other data sets. Despite the large increase in the collective number of variable C. jejuni genes (n = 599) found in the meta-analysis data set, nearly half of these (n = 276) mapped to previously defined variable loci, and it therefore appears that large regions of the C. jejuni genome are genetically stable. A detailed analysis of the microarray data revealed that divergent genes could be differentiated on the basis of the amplitudes of their differential microarray signals. Of 599 variable genes, 122 could be classified as highly divergent on the basis of CGH data. Nearly all highly divergent genes (117 of 122) had divergent neighbors and showed high levels of intraspecies variability. The approach outlined here has enabled us to distinguish global trends of gene conservation in C. jejuni and has enabled us to define this group of genes as a robust set of variable markers that can become the cornerstone of a new generation of genotyping methods that use genome-wide C. jejuni gene variability data.