- 1Department of Population Medicine, Ontario Veterinary College, Guelph, ON, Canada
- 2Independent Researcher, Mississauga, ON, Canada
- 3Centre for Foodborne, Environmental and Zoonotic Infectious Diseases, Public Health Agency of Canada, Guelph, ON, Canada
- 4Quadram Institute Bioscience, Norwich, United Kingdom
- 5National Microbiology Laboratory, Public Health Agency of Canada, Winnipeg, MB, Canada
- 6Department of Medical Microbiology and Infectious Diseases, University of Manitoba, Winnipeg, MB, Canada
- 7National Microbiology Laboratory, Public Health Agency of Canada, Guelph, ON, Canada
- 8Public Health Ontario, Toronto, ON, Canada
- 9Public Health Ontario, Kingston, ON, Canada
- 10Department of Biomedical and Molecular Science, Queen's University, Kingston, ON, Canada
- 11Department of Pathobiology, Ontario Veterinary College, Guelph, ON, Canada
- 12Wyoming Game and Fish Department, Laramie, WY, United States
- 13Department of Veterinary Sciences, University of Wyoming, Laramie, WY, United States
- 14Canadian Wildlife Health Cooperative, Ontario Veterinary College, Guelph, ON, Canada
Non-typhoidal Salmonella infections represent a substantial burden of illness in humans, and the increasing prevalence of antimicrobial resistance among these infections is a growing concern. Using a combination of Salmonella isolate short-read whole-genome sequence data from select human cases, raccoons, livestock and environmental sources, and an epidemiological framework, our objective was to determine if there was evidence for potential transmission of Salmonella and associated antimicrobial resistance determinants between these different sources in the Grand River watershed in Ontario, Canada. Logistic regression models were used to assess the potential associations between source type and the presence of select resistance genes and plasmid incompatibility types. A total of 608 isolates were obtained from the following sources: humans (n = 58), raccoons (n = 92), livestock (n = 329), and environmental samples (n = 129). Resistance genes of public health importance, including blaCMY−2, were identified in humans, livestock, and environmental sources, but not in raccoons. Most resistance genes analyzed were significantly more likely to be identified in livestock and/or human isolates than in raccoon isolates. Based on a 3,002-loci core genome multi-locus sequence typing (cgMLST) scheme, human Salmonella isolates were often more similar to isolates from livestock and environmental sources, than with those from raccoons. Rare instances of serovars S. Heidelberg and S. Enteritidis in raccoons likely represent incidental infections and highlight possible acquisition and dissemination of predominantly poultry-associated Salmonella by raccoons within these ecosystems. Raccoon-predominant serovars were either not identified among human isolates (S. Agona, S. Thompson) or differed by more than 350 cgMLST loci (S. Newport). Collectively, our findings suggest that the rural population of raccoons on swine farms in the Grand River watershed are unlikely to be major contributors to antimicrobial resistant human Salmonella cases in this region.
Introduction
Non-typhoidal Salmonella infections represent a major threat to public health worldwide (1). It has been estimated that non-typhoidal Salmonella cause roughly 93.8 million infections every year, resulting in approximately 155,000 deaths worldwide (2). The majority of infections result in self-limiting gastroenteritis, but invasive bloodstream infections do occur, with children, the elderly, and immunocompromised individuals at the greatest risk of developing life-threatening illness (3). The medical treatment of severe Salmonella infections can also become complicated by the presence of antimicrobial resistance [AMR; (4)]; resistant infections contribute to an increased risk of complications and death, in addition to prolonged recovery times (5, 6). The widespread use of antimicrobials in agriculture and in humans is thought to drive the selection of resistance and virulence among pathogenic organisms such as Salmonella (7, 8). Although wildlife species are increasingly being examined for their potential role in zoonotic diseases (9), their role in the epidemiology of Salmonella and AMR is not entirely clear. To achieve a greater understanding of the drivers of zoonotic infections in humans, and the transmission of AMR within the ecosystem, research approaches using multiple sampling sources, including wildlife and environmental samples, are increasingly being employed (10–14).
Recent work using comparative genomics has provided some epidemiological evidence that wild birds may contribute to the transmission of Salmonella between different human, animal, and environmental sources (15–18). With the vast majority of zoonotic disease investigations of Salmonella in wildlife focused on wild avian species, there is comparatively little work examining terrestrial wild mammals such as raccoons (Procyon lotor). A number of cross-sectional surveys have determined that raccoons can be asymptomatic carriers of non-typhoidal Salmonella, with serovars that overlap with those commonly responsible for causing illness in humans (19–22). Two studies which included genotypic assessments of Salmonella demonstrated the presence of identical pulsed-field gel electrophoresis (PFGE) patterns for certain human and raccoon isolates (23, 24).
In our previous work using whole-genome sequence (WGS) data to assess the potential for transmission of Salmonella, E. coli and AMR determinants between raccoons, soil, and manure pits on swine farms in southern Ontario, the identification of similar or identical subtypes (based on core-genome multi-locus sequence typing) among isolates from these different sources and farms was consistent with potential on-farm and between-farm transmission (25). Our present work assesses this same subset of samples originating from raccoons, soil, and manure pits on swine farms within the context of a greater geographic region in southern Ontario, by comparing those farm isolates to other Salmonella isolates from humans, livestock, and water sources collected during the same time period as part of public health surveillance programs. Our broad research aim was to address current knowledge gaps in our understanding of the role of raccoons in the epidemiology of human Salmonella infections in this region, as well as the movement of associated AMR determinants in this ecosystem, by using a combination of epidemiological modeling and WGS data analysis approaches. More specifically, our objective was to characterize Salmonella isolates from raccoons, humans, livestock and environmental sources, and to assess for patterns in the distribution of AMR determinants (i.e., genes, predicted plasmids) in these different sources using statistical modeling.
Methods
Dataset
Isolates examined within this study were obtained from three different collections of Salmonella cultures that were previously sampled through (1) a wildlife research study (21) and public health surveillance programs led by (2) the Public Health Agency of Canada, FoodNet Canada (https://www.canada.ca/en/public-health/services/surveillance/foodnet-canada.html), and (3) Public Health Ontario. Salmonella isolates collected through public health surveillance of human, livestock, and environmental sources were identified and included if the samples originated from the same geographic region and time period as the previous wildlife study (Wellington-Dufferin and Region of Waterloo Counties in southern Ontario, 2011–2013). The study region, which includes the cities of Waterloo, Kitchener, Guelph, Cambridge and surrounding areas, has a population of ~1 million people and is situated within a region of intensive agriculture. It is also located within the Grand River watershed (6,800 km2), the largest watershed in southern Ontario (21).
Raccoon, Swine Manure Pit, and Soil Isolates
Salmonella isolates for this study component were obtained from a previous three-year repeated cross-sectional study of raccoons on swine farms and in conservation areas (21); only the subset of Salmonella isolates collected on swine farms were selected for sequencing. These samples included fecal swabs from live-trapped raccoons, soil samples collected adjacent to trapped raccoons, and swine manure pit samples. The sampling sites and methods used for live-trapping and processing of raccoons have been described previously (21).
Livestock and Water Isolates
The recruitment and selection of farms for FoodNet Canada has previously been described (26). Livestock samples included pooled samples from swine, chickens (broilers and layers), and cattle (dairy and beef). Samples included both fresh fecal samples and manure pit samples for all species except for layers, from which only fresh fecal samples were obtained. Additional details about the sampling approach used are available in Flockhart et al. (26).
River water samples were obtained from five core water sites upstream of regional drinking water intake in the Grand River watershed. Beach water samples were obtained from local swimming sites in three different recreational areas in southern Ontario (Shade's Mills, Laurel Creek, Elora Gorge). Details regarding the approach and methods used to sample water from rivers and beaches have previously been described by Kadykalo et al. (27).
Human Isolates
Salmonella isolates from human clinical cases collected by Public Health Ontario as part of the reportable disease surveillance program were cultured from blood, stool, and urine specimens (https://www.publichealthontario.ca/en/diseases-and-conditions/infectious-diseases/enteric-foodborne-diseases/salmonellosis).
Previous Culture and Phenotypic Susceptibility Testing
Isolation of Salmonella was previously performed according to methods described in 'National Integrated Enteric Disease Surveillance Program: Sample Collection, Preparation and Laboratory Methodologies, January 2010' (https://www.phac-aspc.gc.ca/foodnetcanada/pdf/lab_sop-eng.pdf). One Salmonella isolate from each sample was sub-cultured and tested further. Presumptive-positive colonies were confirmed via biochemical reactions on triple sugar iron slant and urea slant agars (Becton, Dickinson). Additional confirmation was performed using Salmonella O Poly A-I & Vi antiserum (Becton, Dickinson). All isolates were previously tested for susceptibility to the following 15 antimicrobials by broth microdilution (Sensititre, CMV3AGNF panel, NARMS; Thermo Scientific): gentamicin (GEN), kanamycin (KAN), streptomycin (STR), amoxicillin-clavulanic acid (AMC), cefoxitin (FOX), ceftiofur (TIO), ceftriaxone (CRO), ampicillin (AMP), chloramphenicol (CHL), sulfisoxazole (SOX), trimethoprim-sulfamethoxazole (SXT), tetracycline (TCY), nalidixic acid (NAL), ciprofloxacin (CIP), and azithromycin (AZM). Susceptibility testing of livestock, wildlife, and environmental isolates was previously performed by the Antimicrobial Resistance Reference Laboratory, at the National Microbiology Laboratory (NML) at Guelph (Public Health Agency of Canada, Guelph, ON, Canada) in accordance with methods outlined by The Canadian Integrated Program for Antimicrobial Resistance Surveillance [CIPARS; (28)]. Isolates with intermediate resistance were classified as susceptible.
Selection of Isolates for Whole-Genome Sequencing
All raccoon fecal, soil, and manure pit samples originating from swine farms in the previous wildlife study (21) were selected for sequencing. Previously sequenced human Salmonella isolates originating from the Wellington-Dufferin-Guelph or Region of Waterloo Public Health Counties in Ontario between 2011 and 2013 were included in the present study. Previous selection of human isolates for sequencing as part of CIPARS surveillance included routine surveillance of 11 serovars (4, [5], 12,i:-, Dublin, Enteritidis, Heidelberg, Infantis, Kentucky, Newport, Paratyphi A, Paratyphi B, Typhi, and Typhimurium), in addition to supplementary sequencing by request or for specific research projects investigating AMR [e.g., ciprofloxacin resistance, extended-spectrum beta-lactamases, blaCMY−2 in S. Heidelberg; (29–31)]. For further details regarding CIPARS methods, please see the Design and Methods section of the annual report (32). Previously sequenced Salmonella from water and livestock collected by FoodNet Canada were included in the present study if they had been collected within the Region of Waterloo sentinel site between 2011 and 2013. Previous selection for sequencing of livestock isolates collected by FoodNet Canada was done on the basis of metadata; up to five isolates with identical metadata were sequenced (i.e., same submitter, collection date, commodity, serotype, and phage-type, if applicable); contaminated isolates and isolates without a serotype were not eligible for sequencing.
DNA Extraction, Whole-Genome Sequencing and Genome Assembly
Genomic DNA extractions were performed at the University of Guelph, or at the NML in Winnipeg, Manitoba. Briefly, Salmonella cultures were grown on Mueller Hinton agar and incubated at 35°C overnight. Cultures were then sent to the NML in Winnipeg for DNA extraction and sequencing, or these steps were performed on site, at the University of Guelph and the NML at Guelph, Ontario, respectively. Cultures of 1 ml Salmonella were used as input to the Qiagen DNEasy plant and tissue 96 kit, following manufacturer protocols (Qiagen, Hilden, Germany). Sequencing was performed using the Nextera XT libraries and Illumina MiSeq version 3 (600-cycle kit) or NextSeq550 platforms, according to the manufacturer's protocols. Raw reads were assembled using SPAdes (33), as part of the Shovill pipeline (version 1.0.1; https://github.com/tseeman/shovill) with the following settings: “–minlen 200 –mincov 2; –assembler spades; –trim”.
Analysis of Whole-Genome Assemblies
Prediction of legacy multi-locus sequence types was performed using MLST (version 2.19.0; https://github.com/tseemann/mlst), which uses the Achtman 7-loci scheme for Salmonella enterica (https://pubmlst.org/mlst/). Core genome multi-locus sequence typing (cgMLST) of all isolates was performed using the “fairly simple allele calling” tool fsac (version 1.2.0; https://github.com/dorbarker/fsac) with the 3,002-loci Enterobase scheme (https://enterobase.warwick.ac.uk/). Isolates with 20 or more missing loci were considered poor quality and excluded from any further analyses. Isolates with greater than two missing genotypes for resistant phenotypic test results (by antimicrobial class), and a sequencing depth of <30 using Mash (34) were also excluded from further analyses. A threshold of n = 10 was used as the maximum number of allelic differences present for Salmonella to be considered as possibly belonging to the same strain, consistent with the threshold used by PulseNet (35). Population structure was visualized using cgMLST data and minimum spanning trees generated by the standalone GrapeTree software package [version 1.5; (36)], using the “MSTreeV2” algorithm, which accounts for missing data. For the minimum spanning tree visualizing the overall population of Salmonella, a lenient clustering threshold was used (k = 30) to provide a qualitative assessment of overlap between isolates from different sources, while minimizing unnecessary noise. An additional minimum spanning tree was constructed to examine isolate overlap for only serovars present in both human and raccoon isolates; a stricter clustering threshold (k = 10) was used to provide additional resolution and to minimize clustering of isolates from potentially different strains, consistent with the PulseNet threshold (35). Allelic differences between cgMLST profiles from different sources were calculated using R (version 3.6.3). R code is available at https://github.com/nadinevogt21/salmonella-human-raccoon-WGS-project. In silico identification of serovars was performed using the SISTR command-line tool (version 1.1.1; https://github.com/phac-nml/sistr_cmd), and default settings with the “centroid” allele database.
Acquired resistance genes were identified using CARD-RGI (database version 3.0.8; https://github.com/arpcard/rgi) which uses a curated database that includes optimized identity and coverage settings for each gene (37). The following setting was used: “–exclude nudge” (i.e., loose hits were not nudged to strict). Only acquired resistance genes identified with perfect and strict hits under the “protein homolog model” were reported; regulatory genes, genes requiring mutation for expression, and genes representing chromosomally encoded resistance genes were not reported. In silico identification of plasmid incompatibility groups was performed using Abricate (version 1.0.1; https://github.com/tseemann/abricate) and the PlasmidFinder database (updated May-17 2020) which uses the database from the Center for Genomic Epidemiology, Technical University of Denmark, DTU (https://cge.cbs.dtu.dk/services/PlasmidFinder/). Settings included nucleotide coverage of 70% and percent identity of 98%.
For raccoon and human isolates with fewer than 10 allelic differences based on the cgMLST scheme, Snippy (version 4.4.0; https://github.com/tseemann/snippy) was used to further distinguish genetic differences based on single-nucleotide polymorphisms (SNPs) using the entire genome, with default settings. Salmonella enterica subsp. enterica serovar Heidelberg str. SL476 (Accession No.: NC_011083.1) was used as a reference genome.
Statistical Analyses
Multi-level logistic regression was used to model the odds of the presence of select plasmid incompatibility groups and acquired resistance genes found in Salmonella from different animal, human, and environmental sources. All statistical tests were performed using STATA (STATA Intercooled 14.2; StataCorp, College Station, Texas, USA). Only incompatibility groups or resistance genes with an overall prevalence of at least 10% and at most 90% were analyzed. Incompatibility types that were present only or predominantly (>95%) in one host-restricted serovar were reported descriptively only, due to the potential for complete confounding by serovar. Human and raccoon sources were retained as categories for analyses, and additional sources were grouped together into the following categories: livestock (cattle, swine, chicken), environmental (water, soil). Logistic regression models with source type as an independent variable were initially constructed using the “melogit” command, with a random intercept to account for clustering of isolates obtained from the same raccoon or swine manure pit (applicable to isolates from the wildlife research study). Site-level data for isolates obtained through FoodNet Canada were unavailable, thus, isolates that may have originated from the same farm or water source were treated as independent observations. For multi-level models that did not converge using “melogit”, the model was fitted using QR decomposition of parameters, using the command “meqrlogit”. All models were adjusted for the potential confounding effect of sampling year, if deemed appropriate. Sampling year was retained in the model based on the following criteria: it was statistically significant (p ≤ 0.05), or its removal resulted in a >20% change in the coefficient of the source type variable (38). To examine the fit of multi-level models, Pearson's residuals were assessed for outliers, and the best linear unbiased predictions (BLUPs) were assessed for normality and homoscedasticity. The random intercept was removed if the variance components were very small (<1 ×10−3) and if model fit was not improved based on changes to the Bayesian Information Criterion (BIC) value (38). If the random intercept was not retained in the final model, the model was fit using ordinary logistic regression. For ordinary logistic regression models, model fit was assessed using the Pearson's goodness-of-fit test, in addition to visual assessment of leverage, delta-beta, Pearson, and deviance residuals to identify any potential outliers. Outliers were investigated for recording errors, but otherwise left in the model. Exact logistic regression was used if low effective sample sizes posed estimation issues or to obtain median unbiased estimates if certain categories had zero positive observations; the score method was used to calculate p-values for these models. All tests were two-tailed, and a significance level of α = 0.05 was used.
Results
Description of Dataset
A total of 628 isolates were identified using the study criteria. After excluding isolates with 20 or more missing loci based on the 3,002-loci cgMLST scheme, 608 isolates were available for subsequent analyses from the following sources: human (n = 58), raccoon (n = 92), chicken (n = 214), cattle (n = 60), swine (from FoodNet Canada, n = 34), swine manure pit (from previous wildlife study; n = 21), soil (n = 45), water (n = 76), and beach (n = 8). Among isolates derived from cattle samples, these were roughly equally represented by beef (n = 29/60) and dairy animals (n = 31/60). The majority of chicken isolates were obtained from broilers (n = 204/214), with a handful of isolates from laying hens (n = 10/214). The subset of isolates obtained as part of the previous wildlife study (i.e., raccoon, soil, and swine manure pit samples) totaled 158; accession numbers for these isolates are available in Supplementary File 1. Isolates from all sources were evenly distributed across different sampling years: 2011 (36%), 2012 (34%), 2013 (30%). Since all human isolates that were included in this study exhibited phenotypic resistance to at least one antimicrobial (as part of the selection criteria for previous research projects), the prevalence of phenotypic multiclass resistance (to at least three drug classes) was highest in this source (65.5%, 95%CI: 51.8–77.5%). The prevalence of multiclass resistant and non-resistant isolates in the remaining source categories (all of which included phenotypically resistant and susceptible isolates, and were not selected based on phenotypic resistance) was as follows: livestock (10.6% multiclass resistant, 95%CI: 7.5–14.5%, 57.4% non-resistant, 95%CI: 51.9–62.8%), raccoon (2.2% multiclass resistant, 95%CI: 0.2–7.6%, 96.7% non-resistant, 95%CI: 90.8–99.3%), and environmental (3.1% multiclass resistant, 95%CI: 0.9–7.7%, 85.3% non-resistant, 95%CI: 77.9–90.9%). Results from previous phenotypic susceptibility testing of the 158 isolates from the swine farm subset is available in Bondo et al. (21).
Serovar Distribution
A total of 45 different serovars were identified, and the most commonly identified serovars among all isolates were S. Kentucky (n = 132), S. Heidelberg (n = 99), S. Newport (n = 55), S. Typhimurium (n = 39), S. Agona (n = 38), S. Infantis (n = 33) and S. Enteritidis (n = 30; Figure 1A). Population structure based on cgMLST and the distribution of source types for all isolates is illustrated in Figure 1B. For population structure of all isolates by detailed sampling source, see Supplementary Figure 1. Among raccoon isolates, the most frequently identified serovars were S. Newport (n = 33/92, 36%), S. Agona (n = 16/92, 17%), S. Infantis (n = 9/92, 9.8%), and S. Paratyphi var. Java (n = 9/92, 9.8%). The most prevalent serovars in human isolates were S. Enteritidis (n = 14/58, 24%), S. Heidelberg (n = 14/58, 24%), and S. Kentucky (n = 11/58, 19%); the latter two serovars were infrequently represented, or not present, among environmental and raccoon isolates, but were commonly identified in samples from chickens and cattle (Figure 1B). A heatmap of the distribution of serovars identified in humans in this study that were also identified in other sources is provided in Supplementary Table 1. Salmonella Enteritidis was isolated from all sources except for swine isolates. Several serovars that were infrequently identified in human isolates (<5%) were not identified in any other sample in our study (i.e., S. Litchfield, S. Dublin, S. Chester; Supplementary Table 1). Serovars present in both humans and raccoons, in decreasing frequency relative to raccoons, included: S. Newport, S. Infantis, S. Typhimurium, S. Heidelberg, and S. Enteritidis. The genetic relationships among isolates of S. Newport, S. Infantis, and S. Typhimurium are provided in Figure 2; since only one isolate each of S. Heidelberg and S. Enteritidis was identified in raccoons, these serovars were reported descriptively only. For detailed information regarding the distribution of all serovars by source, see Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-g001.jpg)
Figure 1. Population structure of 608 Salmonella enterica isolates from raccoons, livestock, humans, and environmental sources in southern Ontario, Canada based on 3002-loci cgMLST scheme from Enterobase. Minimum spanning tree created using k = 30 clustering threshold in GrapeTree. (A) Distribution of 45 serovars determined using SISTR. (B) Distribution of source types. Frequency counts are in square brackets. Bubble size is proportional to the number of isolates in each cluster, and each cluster contains isolates differing at a maximum of 30 cgMLST loci.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-g002.jpg)
Figure 2. Population structure of 127 Salmonella enterica isolates from raccoons, livestock, humans, and environmental sources in southern Ontario, Canada based on 3002-loci cgMLST scheme from Enterobase, for serovars S. Newport, S. Typhimurium, and S. Infantis (only serovars identified both in human and raccoon isolates; serovars S. Heidelberg and S. Enteritidis are not depicted here since only one raccoon isolate was identified for each, and these are reported in text). Minimum spanning tree created using k = 10 clustering threshold in GrapeTree. (A) Population structure with serovars determined using SISTR. (B) Distribution of source types. Frequency counts are in square brackets. Bubble size is proportional to the number of isolates in each cluster, and each cluster contains isolates differing at a maximum of 10 cgMLST loci.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t001.jpg)
Table 1. Distribution of Salmonella enterica serovarsa from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n = 608).
MLST Results
A total of 54 sequence types were identified, and the distribution of the 11 most common sequence types by source is provided in Table 2. Seven isolates were not typeable by MLST (two raccoon, one soil, two water, one swine, one chicken) due to a missing allele, or a partial match. The following internationally recognized sequence types (39) were identified in our study population: S. Heidelberg ST15 (n = 95), S. Typhimurium ST19 (n = 33), S. Infantis ST32 (n = 32), ST198 S. Kentucky (n = 15), ST96 S. Schwarzengrund (n = 7), ST10 S. Dublin (n = 2), ST45 S. Newport (n = 2), and ST65 S. Brandenburg (n = 2). As chickens were the most common source type in this study, sequence types from poultry-adapted serovars such as S. Kentucky ST152, and S. Heidelberg ST15 were most common (Table 2). Salmonella Kentucky isolates identified in humans (exclusively ST198) were infrequently identified in chicken isolates (n = 4; Table 2). Among S. Kentucky isolates from chickens, the vast majority (n = 109/114) of these were ST152, and the remaining isolates were ST198 (n = 4/114); one sample could not be typed due to a partial or missing allele. Most S. Newport isolates from raccoons and environmental sources were ST350, whereas the two human S. Newport isolates identified in our study were both ST45. Certain sequence types, including S. Infantis ST32 and S. Typhimurium ST19, were identified in every source. Other sequence types were identified only or predominantly in certain animal sources, but could also be found in water samples: S. Cerro ST367 and S. Give in cattle. For detailed information regarding the distribution of sequence types by source, see Supplementary Table 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t002.jpg)
Table 2. Counts of most frequently identified sequence types of Salmonella enterica isolates from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n = 608).
Population Structure
Serovars identified among both human and raccoon isolates were: S. Enteritidis, S. Heidelberg, S. Infantis, S. Newport, and S. Typhimurium (Figure 2). Based on our previously specified cluster threshold of 10 loci, none of the human and raccoon isolates clustered together for any of the serovars present in both host species: S. Infantis, S. Newport, and S. Typhimurium isolates with a minimum of 73, 359 and 76 different loci, respectively (Figure 2). Single antimicrobial susceptible isolates of S. Enteritidis ST11 and S. Heidelberg ST15 were identified in raccoons and differed from human isolates by a minimum of 61 and four loci, respectively. Further analysis of the single raccoon S. Heidelberg ST15 isolate compared to the five closely related human isolates (<10 cgMLST loci different) revealed a minimum of 15 and a maximum of 21 SNPs that differed between these isolates.
Salmonella Infantis was identified in all source types; all 33 isolates were within 107 loci of one another, whereas all non-human sources containing S. Infantis could be grouped together at a lower threshold of 86 loci. The minimum number of allelic differences to the most closely related S. Infantis human isolate was similar for raccoons (73 loci) and livestock (72 loci) but was lower for environmental samples (56 loci). In contrast, the 30 S. Enteritidis isolates which were found in all source types (but predominantly humans and livestock) clustered together at higher threshold of 441 loci. Among Salmonella Enteritidis isolates (n = 30), livestock and environmental isolates displayed greater similarity to human isolates (a minimum of 12 and 16 allelic differences, respectively) compared to the single raccoon S. Enteritidis isolate (at least 61 allelic differences). Salmonella Typhimurium, also isolated from all sources (n = 39 isolates), displayed a high clustering threshold for all isolates, with all but two isolates grouped together at a threshold of 377 loci. The minimum numbers of allelic differences for S. Typhimurium isolated from raccoons, livestock, and environmental samples compared with human isolates were comparable (76, 79, and 69 loci, respectively). Serovars such as S. Heidelberg (n = 99) that were predominantly isolated from chickens displayed a considerably lower clustering threshold, together with isolates from humans, water, cattle, and the single raccoon S. Heidelberg isolate (n < 37 allelic differences for 85% of isolates). Identical S. Heidelberg were identified between livestock (chicken) and human samples (0 allelic differences), whereas a minimum of at least 16 and four allelic differences were noted for environmental isolates and the single raccoon isolate, respectively.
Two distantly related clusters of S. Kentucky were identified (>2,700 allelic differences). Salmonella Kentucky (n = 132) identified in environmental samples differed from the nearest human isolates by 2,758 loci, compared to the minimum of 139 loci that differed between livestock and human isolates (S. Kentucky was not identified among raccoon isolates). Finally, for serovars such as S. Newport (n = 55) that were mainly identified in environmental sources (e.g., soil) and raccoon isolates, many had identical cgMLST types (0 allelic differences). The two human S. Newport isolates differed from the most closely related raccoon S. Newport isolate at 359 and 2,592 loci, respectively. In contrast, the minimum number of allelic differences between S. Newport from environmental isolates and human isolates was lower (126 loci).
In silico Determination of Plasmid Incompatibility Groups and Acquired AMR Genes
A total of 65 acquired resistance genes and 43 incompatibility groups were identified among all isolates. Commonly identified resistance genes with an overall prevalence >10% were: fosA7, tet(B), aph(6)-Id, aph(3”)-Ib, aac(6')-Iaa, and blaCMY−2 (Table 3). Human isolates contained the greatest diversity of resistance genes, with 58/65 genes identified. Other sources contained fewer resistance genes, in decreasing order: livestock (n = 29/65 genes), environmental (n = 27/65 genes), and raccoon isolates (n = 9/65 genes). A number of human and livestock isolates, along with three water isolates, contained blaCMY−2, which confers the cephamycinase-producing phenotype. Class A and C beta-lactamase genes responsible for conferring the extended-spectrum beta-lactamase (ESBL) phenotype (e.g., blaCTX−M−14, blaCTX−M−65) were only identified among human isolates. One human isolate was contained a mcr-1 gene conferring resistance to colistin. Two human isolates with the ACSSuT (ampicillin, chloramphenicol, streptomycin, sulfamethoxazole, and tetracycline) resistance pattern were identified, one in a ciprofloxacin-resistant S. Typhimurium, and another in a S. Infantis resistant to third-generation cephalosporins.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t003.jpg)
Table 3. Frequencies of acquired antimicrobial resistance genes identified using whole-genome sequence data from Salmonella enterica isolates from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada, 2011–2013 (n = 608).
Few resistance genes were identified among raccoon isolates (Table 3), corresponding with the n = 3/92 isolates that demonstrated phenotypic resistance (one isolate with STR-TCY, and two isolates with SOX-STR-TCY). Most of the resistance genes identified in raccoon isolates conferred resistance to aminoglycosides (n = 11/35 genes) and fosfomycin (n = 19/35 genes); no beta-lactam, phenicol, or quinolone resistance genes were identified in isolates from any of these animals (Table 3). In contrast, environmental isolates contained resistance genes to several major classes of antimicrobials (including beta-lactamases, phenicols, and quinolones; Table 3). The most commonly identified plasmid incompatibility groups are presented in Table 4. Some incompatibility groups were only identified in certain sources (e.g., ColE10 in swine), or appeared predominantly in certain sources (e.g., IncX1-3 in chicken isolates, IncFiip96a in raccoon and environmental isolates). Other incompatibility groups were identified in all source types (e.g., IncI1-Igamma, IncFIIS). For a summary of test sensitivity and specificity of in silico AMR prediction using phenotypic susceptibility testing results as the gold standard, see Supplementary Table 3. The overall sensitivity and specificity of in silico identification of AMR genes were 96.3 and 97.0%, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t004.jpg)
Table 4. Distribution of predicted plasmids identified using whole-genome sequence data by source type for Salmonella enterica isolates from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n = 608).
Associations Between Source Type and Carriage of AMR Genes and Plasmid Incompatibility Groups
Except for IncI1-Igamma, all incompatibility groups and resistance genes examined had a significant association with source type. Certain incompatibility groups were significantly more likely to be identified in human or livestock than in raccoon or environmental isolates (i.e., IncX1-1, colRNAI, colpVC, IncI1-Igamma), whereas others were significantly more likely to be identified in raccoon and environmental isolates compared to human and livestock isolates (i.e., IncFIIS, IncFiip96a; Tables 5, 6). Incompatibility groups IncX1-3 and IncFIB(AP001918) were not modeled since they were exclusively, or almost exclusively (>95%), identified in one serovar (Salmonella Kentucky), and the vast majority (>85%) of these isolates originated from broiler chickens. Most resistance genes analyzed were significantly more likely to be identified in livestock and/or human isolates compared with environmental and/or raccoon isolates (i.e., tetB, ). In particular, the odds of identifying tet(B) and aph(3”)-Ib were highest in livestock isolates, and the odds of identifying blaCMY−2 were highest in human isolates (Tables 5, 6). The odds of identifying blaCMY−2 were also significantly greater in livestock isolates compared to raccoon or environmental isolates (Table 6). Random intercepts were not retained in several models since variance components were extremely small (i.e., ColpVC, tetB, aac6'-Iaa, aph-6-Id, fosA7) or exact logistic regression was used (i.e., blaCMY−2, aph3”-Ib, IncFiip96a).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t005.jpg)
Table 5. Logistic regression models assessing the association between source type and the occurrence of select antimicrobial resistance genes and predicted plasmids in Salmonella enterica isolates from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n = 608).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/840416/fvets-09-840416-HTML-r1/image_m/fvets-09-840416-t006.jpg)
Table 6. Contrasts from logistic regression models assessing the statistically significant associations between source type and the occurrence of select antimicrobial resistance genes and predicted plasmids in Salmonella enterica isolates collected from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n = 608).
Discussion
A diversity of serovars and sequence types were isolated from humans, livestock, raccoons and environmental sources in our study region. Our course-grained epidemiological investigation of Salmonella from various sources in southern Ontario has provided insights into potential transmission between these different sources, and provides further evidence that, although raccoons have the potential to disseminate Salmonella and AMR to humans, their contribution appears to be minimal. Our findings build on existing work suggesting that wildlife play a largely indirect role in the transmission of Salmonella, serving primarily as biological intermediaries between humans, livestock and the environment, rather than acting as a primary driver or major reservoir (40–44). The identification of highly similar or identical cgMLST types (<10 allelic differences) was a rare occurrence in this population of over 600 isolates and was very specific to certain sources and serovars, highlighting potential transmission only in certain contexts (i.e., S. Heidelberg between poultry and humans, S. Newport between raccoons and soil). Likely these instances represent the acquisition of Salmonella by humans through undercooked poultry [for S. Heidelberg; (45)] or as a result of frequent exchange of Salmonella between raccoons and their immediate environment, as has been previously documented (25, 46). Alternatively, overlap in cgMLST types between these different sources could indicate exposure to a common source that was not identified here. Overall, we found little evidence of overlap between Salmonella from raccoons and human cases based on cgMLST, whereas livestock and environmental isolates from soil and water samples showed greater similarity to human isolates. The low prevalence of antimicrobial resistance among these raccoon isolates, and the rare occurrence or absence of genes conferring resistance to high-priority antimicrobials further advances the impression that raccoons are unlikely to be a primary driver or source of multi-drug resistant Salmonella infections for humans or livestock in this region. The selection of human isolates for inclusion in this study, however, was not random, since we included only pre-existing sequence data consisting of antimicrobial resistant human salmonellosis cases reported to public health, which merits consideration in the interpretation of our findings.
Upon closer examination of the distribution of human serovars in our study in relation to previous human Salmonella data from this region, there was some overlap in the top five serovars identified in our study (S. Enteritidis, S. Heidelberg, S. Kentucky, S. Typhimurium, S. Infantis), and the top serovars in humans previously reported by Flockhart et al. (26) in this sentinel site (S. Enteritidis, S. Typhimurium, S. Heidelberg, S. Newport, S. Thompson), and also at the provincial level (S. Enteritidis, S. Typhimurium, S. Heidelberg, S. Typhi, S. Newport; 47). The top three serovars responsible for human cases in our study, S. Enteritidis, S. Heidelberg, and S. Kentucky, accounted for over 65% of all human isolates in this study. The presence of S. Kentucky among top serovars causing human illness likely resulted from biases in the selection of human isolates for sequencing (and thus, inclusion in our study) based on phenotypic resistance. Of note, S. Kentucky previously accounted for <1% of human cases in an earlier study of this same FoodNet Canada sentinel site (2006–2011; 26). All human S. Kentucky isolates in our study were exclusively ST198, and all demonstrated phenotypic resistance to ciprofloxacin. Globally, the emergence of multi-drug resistant travel-associated S. Kentucky ST198 in humans is a growing concern (48); unfortunately, potential travel-related exposures among our human isolates could not be confirmed since epidemiological data concerning travel were unavailable. Apart from human isolates, the only other source containing ST198 isolates (all susceptible to ciprofloxacin) were chickens (n = 4/214), consistent with previous reports (48). Both S. Kentucky and S. Heidelberg were common among chicken isolates, but less common or absent in other sources (indeed, S. Kentucky was not isolated from raccoons). These findings are similar to previous work in this region implicating poultry as the most likely source of infection for human cases of S. Enteritidis and S. Heidelberg based on phage-typing and PFGE (26). In the present study, highly similar or identical S. Enteritidis ST11 and S. Heidelberg ST15 were identified in humans and from livestock (minimum 12 and 0 cgMLST loci, respectively). Single isolates of S. Enteritidis ST11 and S. Heidelberg ST15 identified in two different raccoons most likely represent an infection/colonization in the raccoon which originated from a poultry source, due to the apparent association with poultry isolates, along with previous work suggesting links between broiler chickens and these particular serovars in this region, and elsewhere (26, 49). With the exception of the single S. Heidelberg ST15 raccoon isolate that differed only at four loci (15 SNPs) from the most closely related human S. Heidelberg, none of the Salmonella isolates from raccoons clustered together with isolates from humans (based on our 10 loci threshold), thus suggesting that raccoon-human transmission of antimicrobial resistant Salmonella in our study populations rarely occurred, if at all.
Salmonella Newport was by far the most common serovar identified among raccoons, accounting for just over a third of all raccoon isolates; the two human S. Newport isolates were a different sequence type from those identified in raccoons, with a large number of allelic differences (>350 cgMLST loci), which suggests a different source for these human cases (possibly travel-related). Serovars that were occasionally isolated from raccoons (e.g., S. Typhimurium, S. Infantis) were commonly identified in livestock and environmental sources, suggesting widespread, and frequent dissemination of these generalist Salmonella serovars. Conversely, other Salmonella serovars were restricted to certain livestock hosts (e.g., S. Cerro in cattle), with a handful of water isolates that suggests sporadic environmental contamination with principally livestock-associated serovars (50, 51).
The sampling bias related to the inclusion of human isolates on the basis of demonstrated phenotypic resistance was clearly reflected in the prevalence of phenotypic multiclass resistance in the various sources (i.e., 65% of human isolates, and <15% in other sources), as well as the diversity of resistance genes present. By comparison, the prevalence of multiclass resistance among human Salmonella isolates identified across all of Ontario by CIPARS during the same time period was much lower, ranging between 1.5 and 37.1%, depending on the serovar of interest (52, 53). Among the remaining non-human sources, which included both phenotypically non-resistant and resistant isolates, the prevalence of multidrug resistance (3+ classes) and the diversity of resistance genes identified were highest among livestock isolates, followed by environmental isolates, and, lastly, raccoons. Although raccoons on swine farms were frequently infected with Salmonella [29% prevalence, as previously demonstrated by Bondo et al. (21)], antimicrobial resistance was uncommon (3.2%), and was comparable or lower than other studies of Salmonella in urban and suburban raccoons in Costa Rica and Japan [9–27%; (19, 22)]. Antimicrobial resistance genes were generally highest in human (e.g., blaCMY−2) and livestock sources (e.g., tetB) compared with raccoon and environmental isolates, consistent with a lack of direct exposure of free-ranging wildlife to antimicrobials (other than in wildlife rehabilitation settings). The predominant types of genotypic resistance identified in the subset of swine farm isolates from the previous wildlife research study were to aminoglycosides and fosfomycin (mediated by fosA7, but not tested phenotypically). This unusual finding of fosA7 or fosA7-like genes in wildlife has recently been documented in a white-tailed eagle (Haliaeetus albicilla) in Poland (54), and in Andean condors (Vultur gryphus) in Chile (55), following its initial discovery among S. Heidelberg from broiler chickens in British Colombia, Canada (56). Fosfomycin, an old antibiotic that has once again regained popularity due to the emergence of multi-drug resistant (particularly ciprofloxacin-resistant) lower urinary tract infections in humans (57) is currently considered a critically important antimicrobial by the World Health Organization. Although the identification of a high prevalence (20%) of fosA7 among raccoons in our study represents a potential public health concern, this gene was found in a similar prevalence across human, livestock and environmental samples. Fosfomycin is not widely used in veterinary medicine or animal feed, but is indicated for the treatment of infectious diseases in poultry and swine (58). The revival of fosfomycin in human and in veterinary medicine in response to the emergence of multi-drug resistance may represent a potential source of contamination of fosA7 for wildlife and their environment, thus, continued monitoring and surveillance of this particular type of resistance is needed.
Along with the minimal overlap observed between Salmonella from raccoons and humans based on cgMLST typing, the absence of certain resistance genes of public health concern (e.g., blaCMY−2, blaCARB−3, qnrS1) in the subset of swine farm isolates from the previous wildlife study suggests that raccoons and the swine farms sampled within the Grand River watershed did not play a substantial role in the direct dissemination of high-priority antimicrobial resistant strains of Salmonella to humans in this region. In contrast, water isolates from this watershed, similar to livestock, contained resistance genes to numerous antimicrobial classes (e.g., blaCMY−2, sul1, tetA, qnrS1, cmlA1), which supports the notion that environmental exposures in this region (e.g., the use of recreational areas for swimming) may be important contributors to resistant human Salmonella cases, especially given the high number of conservation areas located within this watershed (n = 11). Salmonella has consistently been isolated from water sources across Canada, including those within this watershed (14, 26, 27, 59, 60); further work is needed, however, to clarify and assess the potential contribution of water sources to cases of human salmonellosis in this region.
The distribution of plasmid incompatibility groups in the dataset as a whole also varied by source for five of the six predicted plasmids assessed, providing evidence of limited transmission between different sources in this region (one would expect widespread transmission and sharing of plasmids to produce a lack of statistical association with the “source type” variable). For instance, the prevalence of colRNAI and colpVC were highest in Salmonella from livestock sources, whereas IncX1-1 was most prevalent in Salmonella from humans. In contrast, IncFiip96a was only identified in Salmonella from raccoon and soil isolates and these particular predicted plasmids do not appear to be associated with resistance genes, given its high prevalence (48%) along with the low prevalence of resistance in these sampling sources (3%). We did not characterize plasmid mobility, thus it is unknown whether these incompatibility groups representing plasmids have the potential to be an important factor in the circulation of antimicrobial resistance genes, genes conferring resistance to heavy metals or disinfectants, or relevant virulence factors (not assessed here).
Limitations
As previously indicated, a major limitation of our study relates to systematic bias in the selection of human isolates for sequencing as part of public health surveillance programs. Previous work documented the presence of S. Thompson and S. Newport within the top ten serovars causing human illness in Ontario and in our study region between 2006 and 2011 (26, 47). Without a comprehensive examination of all human Salmonella cases in the region (including those not reported to public health authorities), it remains to be determined whether milder, non-resistant Salmonella infections in humans could be related to wildlife and environmental sources. For instance, S. Thompson and S. Newport were both highly prevalent serovars in raccoon and soil isolates from swine farms, both have previously been documented in water isolates in this study region (9, 26, 27), and they have also consistently been documented as top serovars causing human illness in this region (26, 47). Our collection of resistant human isolates only contained two S. Newport and no S. Thompson isolates; it remains to be seen whether other human infections with these serovars—not captured here—may be originating from wildlife and related environmental sources. In addition, associations and prevalence measures involving livestock isolates included in this study may also be impacted by selection bias, since the inclusion and continuing enrollment of farms in the FoodNet Canada surveillance program was not random, as previously discussed by Flockhart et al. (26). Similarly, the population of raccoons in our study should not be considered representative of all raccoons in the Grand River watershed, as a result of trapping animals (61), and focusing on a rural population of raccoons in agricultural areas.
Additional limitations relate to plasmid determination using incompatibility groups, and lack of characterization of plasmid mobility. Previous work has demonstrated that plasmid distribution varies by serovar (15), however, the aim of our work was to characterize patterns in the occurrence of genes and incompatibility groups representing plasmids in different sources to address potential transmission and propose primary sources. Due to the sheer number of serovars with few isolates that prevented many models from converging when serovar was included as a random intercept (data not shown), we were unable to control for the potential confounding effect of serovar on the occurrence of different plasmid types or resistance genes in different sources, which may be a particularly important factor for plasmids that are non-mobile or that have limited mobility.
Finally, our use of “livestock” as a broad category in statistical models to capture chicken, cattle, and swine isolates prevented us from identifying species-specific differences. However, this study provides a preliminary examination of potential transmission of Salmonella and AMR determinants between broad groups of animals, and the lumping of livestock into one category avoided potential misclassification of isolates associated with the surveillance-collected fecal samples from multi-species farms. Among these surveillance-collected isolates, fecal samples may have been obtained fresh, or from standing pooled manure piles; in the future, efforts focused on the impact of these sampling methods will be important for the utilization and interpretation of heterogeneous sampling methods and sources.
Conclusions
This study offers new insights about the epidemiology of Salmonella and associated AMR from human, animal, and environmental sources in a populous region of southern Ontario. Overall, raccoons sampled on swine farms in this study were rarely infected with internationally recognized Salmonella lineages (39), or with antimicrobial resistant isolates (<5%). Conversely, resistance was more commonly identified in livestock and water isolates, and a number of resistance genes of public health importance (e.g., blaCMY−2, blaCARB−3, qnrS1) were identified in some of these sources. Although other studies have demonstrated strong links between Salmonella isolates from wildlife and humans (17, 62, 63), our assessment of microbial population structure based on cgMLST suggests that raccoons captured on swine farms are an unlikely source of resistant Salmonella for humans in this study region. However, the occurrence of certain livestock-associated Salmonella serovars and various resistance determinants in raccoons supports their role as sentinels and potential disseminators of these organisms within the environment. Further research, including a comprehensive examination of all human cases (including phenotypically susceptible ones), and examination of additional raccoon populations (both urban and rural) may help to provide additional insight into the potential role of raccoons in human Salmonella infections.
Data Availability Statement
Previous epidemiological data from the wildlife study are available from the Agri-environmental Research Data Repository: http://dx.doi.org/10.5887/AERDR/10864/12074. Sequence data for isolates from the wildlife study have been deposited to Genbank and are available under BioProject PRJNA745182; accession numbers for those isolates are available in Supplementary Information. Interested researchers who would like to access FoodNet Canada or Public Health Ontario data will need to meet their requirements and make a data sharing agreement with those organizations.
Ethics Statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. The animal study was reviewed and approved by Research Ethics Board, University of Guelph. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
NV, BH, DP, CJ, EP, RR-S, and MM were involved in the conceptualization of the study. Funding acquisition was performed by CJ, DP, MM, EP, RR-S, and VA. DP, MM, EP, RR-S, CJ, and NR supervised the project. Curation of previously collected data was performed by KB, KZ, SK, AB, and JR. KB and SA were involved in study design of and field data collection for the previous wildlife study. NJ oversaw culture of isolates, and microbiological methodology for previous wildlife study. NV, AB, MM, JR, KZ, SK, VA, AM, EP, CJ, and JN assisted with acquisition of sequence data from public health, and coordination of sequencing of wildlife study isolates. Curation, cleaning, and analysis of WGS data was performed by NV, BH, and AV. Data analyses were performed by NV, BH, and AV, with assistance from DP and MM. Writing of the original draft was performed by NV. All authors have read and approved the final manuscript.
Funding
Funding was provided by the Ontario Ministry of Agriculture, Food, and Rural Affairs (OMAFRA), through the Ontario Agri-Food Innovation Alliance (UofG2016-2642). NV received stipend funding from the Ontario Veterinary College, University of Guelph and a National Sciences and Engineering Research Council Postgraduate Scholarship-Doctoral.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
Erin Harkness, Samantha Kagan, and Mary Thompson assisted with collection of field data, and Bryan Bloomfield collected swine manure pit samples and engaged the participation of landowners for the previous wildlife study. Field samples from the wildlife study were submitted by Tami Harvey and Barbara Jefferson. Culture and submission of samples from the previous wildlife study for sequencing was performed and coordinated by the McEwen Lab at the Canadian Research Institute for Food Safety (Sarah Martz, laboratory technicians, and co-op students), University of Guelph. Sequencing was performed at the National Microbiology Laboratory at Guelph and the National Microbiology Laboratory at Winnipeg, Public Health Agency of Canada (Bob Holtslander, Irene Young, Matt Cook, Ketna Mistry from Guelph Reference Services, OIE Reference Lab for Salmonellosis, Andrea Desruisseau, Chad Gill). Laura Mataseje and Manon Fleury assisted with acquisition of sequence data. Collection of water samples was performed by the Ontario Ministry of the Environment.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.840416/full#supplementary-material
Supplementary File 1. Core-genome multi-locus sequence types, legacy multi-locus sequence types, cgMLST serovars, and gene accession numbers for Salmonella enterica isolates from previous wildlife study.
Supplementary Figure 1. Population structure of 608 Salmonella enterica isolates from raccoons, livestock, humans, and environmental sources in southern Ontario, Canada based on 3,002-loci cgMLST scheme from Enterobase, showing distribution of detailed sampling sources.
Supplementary Table 1. Heatmap of all serovars identified in humans and their proportional occurrence in Salmonella enterica isolates from raccoons, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n=608). Counts of isolates and serovars are indicated on the left. Heatmap indicates the proportion of isolates from each source type in common with human serovars, with darker shading representing a higher proportion. Except for human sources, row totals of proportions do not add up to 100% since only serovars identified in human isolates are presented here.
Supplementary Table 2. Distribution of Salmonella enterica legacy multi-locus sequence types for isolates obtained from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n=608).
Supplementary Table 3. Test sensitivity and specificity for in silico identification of acquired antimicrobial resistance genes in Salmonella enterica isolates from raccoons, humans, livestock, and environmental sources in southern Ontario, Canada 2011–2013 (n=608).
References
1. World Health Organization. WHO estimates of the global burden of foodborne diseases. Foodborne disease burden epidemiology reference group, 2007-2015. (2015). Available online at: https://apps.who.int/iris/bitstream/handle/10665/199350/9789241565165_eng.pdf?sequence=1 (accessed August 10, 2020).
2. Majowicz SE, Musto J, Scallan E, Angulo FJ, Kirk M, O'Brien SJ, et al. The global burden of nontyphoidal Salmonella gastroenteritis. Clinical Infect Dis. (2010) 50:882–9. doi: 10.1086/650733
3. GBD 2017 Non-Typhoidal Salmonella Invasive Disease Collaborators. The global burden of non-typhoidal salmonella invasive disease: a systematic analysis for the Global Burden of Disease Study 2017. Lancet Infect Dis. (2019) 19:1312–24. doi: 10.1016/S1473-3099(19)30418-9
4. Michael GB, Schwarz S. Antimicrobial resistance in zoonotic nontyphoidal Salmonella: an alarming trend? Clin Microbiol Infect. (2016) 22:968–74. doi: 10.1016/j.cmi.2016.07.033
5. Kern WV, Rieg S. Burden of bacterial bloodstream infection-a brief update on epidemiology and significance of multidrug-resistant pathogens. Clin Microbiol Infect. (2020) 26:151–7. doi: 10.1016/j.cmi.2019.10.031
6. Tanwar J, Das S, Fatima Z, Hameed S. Multidrug resistance: an emerging crisis. Interdiscip Perspect Infect Dis. (2014) 2014:541340. doi: 10.1155/2014/541340
7. Kumar SB, Arnipalli SR, Ziouzenkova O. Antibiotics in food chain: the consequences for antibiotic resistance. Antibiotics. (2020) 9:688. doi: 10.3390/antibiotics9100688
8. Bawn M, Alikhan NF, Thilliez G, Kirkwood M, Wheeler NE, Petrovska L, et al. Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation. PLoS Genet. (2020) 16:e1008850. doi: 10.1371/journal.pgen.1008850
9. Torres RT, Carvalho J, Cunha MV, Serrano E, Palmeira JD, Fonseca C. Temporal and geographical research trends of antimicrobial resistance in wildlife-a bibliometric analysis. One Health. (2020) 11:100198. doi: 10.1016/j.onehlt.2020.100198
10. Munck N, Smith J, Bates J, Glass K, Hald T, Kirk MD. Source attribution of Salmonella in macadamia nuts to animal and environmental reservoirs in Queensland, Australia. Foodborne Pathog Dis. (2020) 17:357–64. doi: 10.1089/fpd.2019.2706
11. Aung KT, Chen HJ, Chau ML, Yap G, Lim XF, Humaidi M, et al. Salmonella in retail food and wild birds in Singapore-Prevalence, antimicrobial resistance, and sequence types. Int J Environ Res Public Health. (2019) 16:4235. doi: 10.3390/ijerph16214235
12. Guyomard-Rabenirina S, Reynaud Y, Pot M, Albina E, Couvin D, Ducat C, et al. Antimicrobial resistance in wildlife in Guadeloupe (French West Indies): distribution of a single bla[[sb]]CTX-M-1[[/S]]/IncI1/ST3 plasmid among humans and wild animals. Front Microbiol. (2020) 11:1524. doi: 10.3389/fmicb.2020.01524
13. De Lucia A, Rabie A, Smith RP, Davies R, Ostanello F, Ajayi D, et al. Role of wild birds and environmental contamination in the epidemiology of Salmonella infection in an outdoor pig farm. Vet Microbiol. (2018) 227:148–54. doi: 10.1016/j.vetmic.2018.11.003
14. Jokinen CC, Koot J, Cole L, Desruisseau A, Edge TA, Khan IU, et al. The distribution of Salmonella enterica serovars and subtypes in surface water from five agricultural regions across Canada. Water Res. (2015) 76:120–31. doi: 10.1016/j.watres.2015.02.038
15. Pornsukarom S, van Vliet AHM, Thakur S. Whole genome sequencing analysis of multiple Salmonella serovars provides insights into phylogenetic relatedness, antimicrobial resistance, and virulence markers across humans, food animals and agriculture environmental sources. BMC Genomics. (2018) 19:801. doi: 10.1186/s12864-018-5137-4
16. Bloomfield SJ, Benschop J, Biggs PJ, Marshall JC, Hayman D, Carter PE, et al. Genomic analysis of Salmonella enterica serovar Typhimurium DT160 associated with a 14-year outbreak, New Zealand, 1998-2012. Emerg Infect Dis. (2017) 23:906–13. doi: 10.3201/eid2306.161934
17. Toro M, Retamal P, Ayers S, Barreto M, Allard M, Brown EW, et al. Whole-genome sequencing analysis of Salmonella enterica serovar Enteritidis isolates in Chile provides insights into possible transmission between gulls, poultry, and humans. Appl Environ Microbiol. (2016) 82:6223–32. doi: 10.1128/AEM.01760-16
18. Mather AE, Lawson B, de Pinna E, Wigley P, Parkhill J, Thomson NR, et al. Genomic analysis of Salmonella enterica serovar Typhimurium from wild passerines in England and Wales. Appl Environ Microbiol. (2016) 82:6728–35. doi: 10.1128/AEM.01660-16
19. Baldi M, Barquero Calvo E, Hutter SE, Walzer C. Salmonellosis detection and evidence of antibiotic resistance in an urban raccoon population in a highly populated area, Costa Rica. Zoonoses Public Hlth. (2019) 66:852–60. doi: 10.1111/zph.12635
20. Rainwater KL, Marchese K, Slavinski S, Humberg LA, Dubovi EJ, Jarvis JA, et al. Health survey of free-ranging raccoons (Procyon lotor) in Central Park, New York, USA: implications for human and domestic animal health. J Wildl Dis. (2017) 53:272–84. doi: 10.7589/2016-05-096
21. Bondo KJ, Pearl DL, Janecko N, Boerlin P, Reid-Smith RJ, Parmley J, et al. Impact of season, demographic and environmental factors on Salmonella occurrence in raccoons (Procyon lotor) from swine farms and conservation areas in southern Ontario. PLoS ONE. (2016) 11:e0161497. doi: 10.1371/journal.pone.0161497
22. Lee K, Iwata T, Nakadai A, Kato T, Hayama S, Taniguchi T, et al. Prevalence of Salmonella, Yersinia and Campylobacter spp. in feral raccoons (Procyon lotor) and masked palm civets (Paguma larvata) in Japan. Zoonoses Public Health. (2011) 58:424–31. doi: 10.1111/j.1863-2378.2010.01384.x
23. Very KJ, Kirchner MK, Shariat N, Cottrell W, Sandt CH, Dudley EG, et al. Prevalence and spatial distribution of Salmonella infections in the Pennsylvania raccoon (Procyon lotor). Zoonoses Public Health. (2016) 63:223–33. doi: 10.1111/zph.12222
24. Compton JA, Baney JA, Donaldson SC, Houser BA, San Julian GJ, Yahner RH, et al. Salmonella infections in the common raccoon (Procyon lotor) in western Pennsylvania. J Clin Microbiol. (2008) 46:3084–6. doi: 10.1128/JCM.00685-08
25. Vogt NA, Hetman BM, Pearl DL, Vogt AA, Reid-Smith RJ, Parmley EJ, et al. Using whole-genome sequence data to examine the epidemiology of Salmonella, Escherichia coli and associated antimicrobial resistance in raccoons (Procyon lotor), swine manure pits, and soil samples on swine farms in southern Ontario, Canada. PLoS ONE. (2021) 16:e0260234. doi: 10.1371/journal.pone.0260234
26. Flockhart L, Pintar K, Cook A, McEwen S, Friendship R, Kelton D, et al. Distribution of Salmonella in humans, production animal operations and a watershed in a FoodNet Canada sentinel site. Zoonoses Public Health. (2017) 64:41–52. doi: 10.1111/zph.12281
27. Kadykalo S, Thomas J, Parmley EJ, Pintar K, Fleury M. Antimicrobial resistance of Salmonella and generic Escherichia coli isolated from surface water samples used for recreation and a source of drinking water in southwestern Ontario, Canada. Zoonoses Public Health. (2020) 67:566–75. doi: 10.1111/zph.12720
28. Government of Canada. Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) 2015 annual report. (2015). Guelph, ON: Public Health Agency of Canada.
29. Mulvey MR, Boyd DA, Finley R, Fakharuddin K, Langner S, Allen V, et al. Ciprofloxacin-resistant Salmonella enterica serovar Kentucky in Canada. Emerg Infect Dis. (2013) 19:999–1001. doi: 10.3201/eid1906.121351
30. Mulvey MR, Bharat A, Boyd DA, Irwin RJ, Wylie J. Characterization of a colistin-resistant Salmonella enterica 4,[5],12:i:- harbouring mcr-32 on a variant IncHI-2 plasmid identified in Canada. J Med Microbiol. (2018) 67:1673–5. doi: 10.1099/jmm.0.000854
31. Mangat CS, Bekal S, Avery BP, Côté G, Daignault D, Doualla-Bell F, et al. Genomic investigation of the emergence of invasive multidrug-resistant Salmonella enterica serovar Dublin in humans and animals in Canada. Antimicrob Agents Chemother. (2019) 63:e00108–19. doi: 10.1128/AAC.00108-19
32. Government of Canada Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) 2018: Integrated Findings. Guelph, Ontario: Public Health Agency of Canada (2020).
33. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. (2012) 19:455–77. doi: 10.1089/cmb.2012.0021
34. Ondov BD, Starrett GJ, Sappington A, Kostic A, Koren S, Buck CB, et al. Mash screen: high-throughput sequence containment estimation for genome discovery. Genome Biol. (2019) 20:232. doi: 10.1186/s13059-019-1841-x
35. Tolar B, Joseph LA, Schroeder MN, Stroika S, Ribot EM, Hise KB, et al. An overview of PulseNet USA databases. Foodborne Pathog Dis. (2019) 16:457–62. doi: 10.1089/fpd.2019.2637
36. Zhou Z, Alikhan NF, Sergeant MJ, Luhmann N, Vaz C, Francisco AP, et al. GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens. Genome Res. (2018) 28:1395–404. doi: 10.1101/gr.232397.117
37. Alcock BP, Raphenya AR, Lau TTY, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. (2020) 48:D517–25. doi: 10.1093/nar/gkz935
38. Dohoo I, Martin W, Stryhn H. Veterinary Epidemiologic Research, 3rd ed. Charlottetown, Prince Edward Island (2014).
39. Monte DF, Lincopan N, Berman H, Cerdeira L, Keelara S, Thakur S, et al. Genomic features of high-priority Salmonella enterica serovars circulating in the food production chain, Brazil, 2000-2016. Sci Rep. (2019) 9:11058. doi: 10.1038/s41598-019-45838-0
40. Uelze L, Bloch A, Borowiak M, Grobbel M, Deneke C, Fischer M, et al. What WGS reveals about Salmonella enterica subsp. enterica in wildlife in Germany. Microorganisms. (2021) 9:1911. doi: 10.3390/microorganisms9091911
41. Ahlstrom CA, Frick A, Pongratz C, Spink K, Xavier C, Bonnedahl j, et al. Genomic comparison of carbapenem-resistant Enterobacteriaceae from humans and gulls in Alaska. J Glob Antimicrob Resist. (2021) 25:23–5. doi: 10.1016/j.jgar.2021.02.028
42. La Tela I, Peruzy MF, D'Alessio N, Di Nocera F, Casalinuovo F, Carullo MR, et al. Serotyping and evaluation of antimicrobial resistance of Salmonella strains detected in wildlife and natural environments in southern Italy. Antibiotics. (2021) 10:353. doi: 10.3390/antibiotics10040353
43. Espunyes J, Cabezón O, Dias-Alves A, Miralles P, Ayats T, Cerdà-Cuéllar M. Assessing the role of livestock and sympatric wild ruminants in spreading antimicrobial resistant Campylobacter and Salmonella in alpine ecosystems. BMC Vet Res. (2021) 17:79. doi: 10.1186/s12917-021-02784-2
44. Skov MN, Madsen JJ, Rahbek C, Lodal J, Jespersen JB, Jørgensen JC, et al. Transmission of Salmonella between wildlife and meat-production animals in Denmark. J Appl Microbiol. (2008) 105:1558–68. doi: 10.1111/j.1365-2672.2008.03914.x
45. Currie A, MacDougall L, Aramini J, Gaulin C, Ahmed R, Isaacs S. Frozen chicken nuggets and strips and eggs are leading risk factors for Salmonella Heidelberg infections in Canada. Epidemiol Infect. (2005)133:809–16. doi: 10.1017/S0950268805004383
46. Bondo KJ, Pearl DL, Janecko N, Boerlin P, Reid-Smith RJ, Parmley J, et al. Epidemiology of Salmonella on the paws and in the faeces of free-ranging raccoons (Procyon lotor) in southern Ontario, Canada. Zoonoses Public Hlth. (2016) 63:303–10. doi: 10.1111/zph.12232
47. Public Public Health Agency of Canada 2006–2011: National Enteric Surveillance Program (NESP) Annual Summaries 2006-2011. Ottawa, ON: Government of Canada.
48. Hawkey J, Le Hello S, Doublet B, Granier SA, Hendriksen RS, Fricke WF, et al. Global phylogenomics of multidrug-resistant Salmonella enterica serotype Kentucky ST198. Microb Genom. (2019) 5:e000269. doi: 10.1099/mgen.0.000269
49. EFSA EFSA Panel on Biological Hazards (EFSA BIOHAZ Panel), Koutsoumanis K, Allende A, Alvarez-Ordóñez A, Bolton D, Bover-Cid S, et al. Salmonella control in poultry flocks and its public health impact. EFSA J. (2019) 17:e05596. doi: 10.2903/j.efsa.2019.5596
50. Cummings KJ, Warnick LD, Elton M, Rodriguez-Rivera LD, Siler JD, Wright EM, et al. Salmonella enterica serotype Cerro among dairy cattle in New York: An emerging pathogen? Foodborne Pathog Dis. (2010) 7:659–65. doi: 10.1089/fpd.2009.0462
51. Davies PR, Morrow WE, Jones FT, Deen J, Fedorka-Cray PJ, Harris IT. Prevalence of Salmonella in finishing swine raised in different production systems in North Carolina, USA. Epidemiol Infect. (1997) 119:237–44. doi: 10.1017/S095026889700784X
52. Government of Canada (2012). Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) 2012 Annual Report, Chapter 2–Antimicrobial resistance. Guelph, Ontario: Public Health Agency of Canada (2014).
53. Government of Canada (2013). Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) 2013 Annual Report, Chapter 2–Antimicrobial resistance. Guelph, Ontario: Public Health Agency of Canada (2015).
54. Skarzyńska M, Zaja CM, Bomba A, Bocian Ł, Kozdruń W, Polak M, et al. Antimicrobial resistance glides in the sky-free-living birds as a reservoir of resistant Escherichia coli with zoonotic potential. Front Microbiol. (2021) 12:656223. doi: 10.3389/fmicb.2021.656223
55. Fuentes-Castillo D, Esposito F, Cardoso B, Dalazen G, Moura Q, Fuga B, et al. Genomic data reveal international lineages of critical priority Escherichia coli harbouring wide resistome in Andean condors (Vultur gryphus Linnaeus, 1758). Mol Ecol. (2020) 29:1919–35. doi: 10.1111/mec.15455
56. Rehman MA, Yin X, Persaud-Lachhman MG, Diarra MS. First detection of a Fosfomycin resistance gene, fosA7, in Salmonella enterica serovar Heidelberg isolated from broiler chickens. Antimicrob Agents Chemother. (2017) 61:e00410–17. doi: 10.1128/AAC.00410-17
57. Gardiner BJ, Stewardson AJ, Abbott IJ, Peleg AY. Nitrofurantoin and fosfomycin for resistant urinary tract infections: old drugs for emerging problems. Aust Prescr. (2019) 42:14–9. doi: 10.18773/austprescr.2019.002
58. Pérez DS, Tapia MO, Soraci AL. Fosfomycin: uses and potentialities in veterinary medicine. Open Vet J. (2014) 4:26–43.
59. Thomas JL, Slawson RM, Taylor WD. Salmonella serotype diversity and seasonality in urban and rural streams. J Appl Microbiol. (2013) 114:907–22. doi: 10.1111/jam.12079
60. Edge TA, El-Shaarawi A, Gannon V, Jokinen C, Kent R, Khan IU, et al. Investigation of an Escherichia coli environmental benchmark for waterborne pathogens in agricultural watersheds in Canada. J Environ Qual. (2012) 41:21–30. doi: 10.2134/jeq2010.0253
61. Vogt NA, Stevens CPG, Pearl DL, Taboada EN, Jardine CM. Generalizability and comparability of prevalence estimates in the wild bird literature: methodological and epidemiological considerations. Anim Health Res Rev. (2020) 21:89–95. doi: 10.1017/S1466252320000043
62. Ford L, Ingle D, Glass K, Veitch M, Williamson DA, Harlock M, et al. Whole-genome sequencing of Salmonella Mississippi and Typhimurium Definitive Type 160, Australia and New Zealand. Emerg Infect Dis. (2019) 25:1690–7. doi: 10.3201/eid2509.181811
Keywords: antimicrobial resistance, foodborne illness, Procyon lotor, raccoon, Salmonella, whole-genome sequencing, wildlife
Citation: Vogt NA, Hetman BM, Vogt AA, Pearl DL, Reid-Smith RJ, Parmley EJ, Kadykalo S, Janecko N, Bharat A, Mulvey MR, Ziebell K, Robertson J, Nash J, Allen V, Majury A, Ricker N, Bondo KJ, Allen SE and Jardine CM (2022) Rural Raccoons (Procyon lotor) Not Likely to Be a Major Driver of Antimicrobial Resistant Human Salmonella Cases in Southern Ontario, Canada: A One Health Epidemiological Assessment Using Whole-Genome Sequence Data. Front. Vet. Sci. 9:840416. doi: 10.3389/fvets.2022.840416
Received: 21 December 2021; Accepted: 24 January 2022;
Published: 25 February 2022.
Edited by:
Min Yue, Zhejiang University, ChinaReviewed by:
Chien-Shun Chiou, Taiwan Centers for Disease Control (CDC), TaiwanMadelaine Norström, Norwegian Veterinary Institute (NVI), Norway
Copyright © 2022 Vogt, Hetman, Vogt, Pearl, Reid-Smith, Parmley, Kadykalo, Janecko, Bharat, Mulvey, Ziebell, Robertson, Nash, Allen, Majury, Ricker, Bondo, Allen and Jardine. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nadine A. Vogt, nvogt@uoguelph.ca
†Present Address: Vanessa Allen, Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, ON, Canada
Kristin J. Bondo, Pennsylvania Cooperative Fish and Wildlife Unit, The Pennsylvania State University, University Park, TX, United States