Reoccurring Escherichia coli O157:H7 Strain Linked to Leafy Greens-Associated Outbreaks, 2016-2019.

Genomic characterization of an Escherichia coli O157:H7 strain linked to leafy greens-associated outbreaks dates its emergence to late 2015. One clade has notable accessory genomic content and a previously described mutation putatively associated with increased arsenic tolerance. This strain is a reoccurring, emerging, or persistent strain causing illness over an extended period.

Conjugal Transfer, Whole-Genome Sequencing, and Plasmid Analysis of Four mcr-1-Bearing Isolates from U.S. Patients.

Four Enterobacteriaceae clinical isolates bearing mcr-1 gene-harboring plasmids were characterized. All isolates demonstrated the ability to transfer colistin resistance to Escherichia coli; plasmids were stable in conjugants after multiple passages on nonselective media. mcr-1 was located on an IncX4 (n = 3) or IncN (n = 1) plasmid. The IncN plasmid harbored 13 additional antimicrobial resistance genes. Results indicate that the mcr-1-bearing plasmids in this study were highly transferable in vitro and stable in the recipients.

Evaluating the occurrence of Escherichia albertii in chicken carcass rinses by PCR, Vitek analysis, and sequencing of the rpoB gene.

Escherichia albertii is a recently described species that has been associated with gastroenteritis in humans and with healthy and ill birds. Most recently, it has been identified as the causative agent in a food-borne outbreak in Japan. The distribution and clinical importance of E. albertii are not well studied because its importance is unclear. Culture methods for clinical isolation frequently miss E. albertii or incorrectly identify it as Shigella spp., Escherichia coli, or Hafnia alvei. This study was designed to determine if E. albertii could be recovered from chicken carcass rinses collected at slaughter during a 1-year period from November 2009 until October 2010. Colonies were isolated from chicken carcass rinses and tested by PCR for the presence or absence of clpX, lysP, mdh, intimin (eae), Shiga toxins 1 and 2 (stx1, stx2, and stx2f), heat-stable enterotoxin A (staA), and cytolethal distending toxins 1 and 2 (cdtB) genes. Sixty-five isolates were analyzed by sequencing a section of the rpoB gene. Analysis of the rpoB gene sequences revealed 14 fixed differences between E. albertii and other, closely related organisms. The fixed differences found in the rpoB gene could aid in future discrimination of E. albertii from closely related bacteria.

Genetic Diversity in Salmonella enterica in Outbreaks of Foodborne and Zoonotic Origin in the USA in 2006-2017.

Whole genome sequencing is replacing traditional laboratory surveillance methods as the primary tool to track and characterize clusters and outbreaks of the foodborne and zoonotic pathogen Salmonella enterica (S. enterica). In this study, 438 S. enterica isolates representing 35 serovars and 13 broad vehicle categories from one hundred epidemiologically confirmed outbreaks were evaluated for genetic variation to develop epidemiologically relevant interpretation guidelines for Salmonella disease cluster detection. The Illumina sequences were analyzed by core genome multi-locus sequence typing (cgMLST) and screened for antimicrobial resistance (AR) determinants and plasmids. Ninety-three of the one hundred outbreaks exhibited a close allele range (less than 10 allele differences with a subset closer than 5). The remaining seven outbreaks showed increased variation, of which three were considered polyclonal. A total of 16 and 28 outbreaks, respectively, showed variations in the AR and plasmid profiles. The serovars Newport and I 4,[5],12:i:-, as well as the zoonotic and poultry product vehicles, were overrepresented among the outbreaks, showing increased variation. A close allele range in cgMLST profiles can be considered a reliable proxy for epidemiological relatedness for the vast majority of S. enterica outbreak investigations. Variations associated with mobile elements happen relatively frequently during outbreaks and could be reflective of changing selective pressures.

Identification and Characterization of ten Escherichia coli Strains Encoding Novel Shiga Toxin 2 Subtypes, Stx2n as Well as Stx2j, Stx2m, and Stx2o, in the United States.

The sharing of genome sequences in online data repositories allows for large scale analyses of specific genes or gene families. This can result in the detection of novel gene subtypes as well as the development of improved detection methods. Here, we used publicly available WGS data to detect a novel Stx subtype, Stx2n in two clinical E. coli strains isolated in the USA. During this process, additional Stx2 subtypes were detected; six Stx2j, one Stx2m strain, and one Stx2o, were all analyzed for variability from the originally described subtypes. Complete genome sequences were assembled from short- or long-read sequencing and analyzed for serotype, and ST types. The WGS data from Stx2n- and Stx2o-producing STEC strains were further analyzed for virulence genes pro-phage analysis and phage insertion sites. Nucleotide and amino acid maximum parsimony trees showed expected clustering of the previously described subtypes and a clear separation of the novel Stx2n subtype. WGS data were used to design OMNI PCR primers for the detection of all known stx1 (283 bp amplicon), stx2 (400 bp amplicon), intimin encoded by eae (221 bp amplicon), and stx2f (438 bp amplicon) subtypes. These primers were tested in three different laboratories, using standard reference strains. An analysis of the complete genome sequence showed variability in serogroup, virulence genes, and ST type, and Stx2 pro-phages showed variability in size, gene composition, and phage insertion sites. The strains with Stx2j, Stx2m, Stx2n, and Stx2o showed toxicity to Vero cells. Stx2j carrying strain, 2012C-4221, was induced when grown with sub-inhibitory concentrations of ciprofloxacin, and toxicity was detected. Taken together, these data highlight the need to reinforce genomic surveillance to identify the emergence of potential new Stx2 or Stx1 variants. The importance of this surveillance has a paramount impact on public health. Per our description in this study, we suggest that 2017C-4317 be designated as the Stx2n type-strain.

Genomic and Phenotypic Characterization of Shiga Toxin-Producing Escherichia albertii Strains Isolated from Wild Birds in a Major Agricultural Region in California.

Escherichia albertii is an emerging foodborne pathogen. To better understand the pathogenesis and health risk of this pathogen, comparative genomics and phenotypic characterization were applied to assess the pathogenicity potential of E. albertii strains isolated from wild birds in a major agricultural region in California. Shiga toxin genes stx2f were present in all avian strains. Pangenome analyses of 20 complete genomes revealed a total of 11,249 genes, of which nearly 80% were accessory genes. Both core gene-based phylogenetic and accessory gene-based relatedness analyses consistently grouped the three stx2f-positive clinical strains with the five avian strains carrying ST7971. Among the three Stx2f-converting prophage integration sites identified, ssrA was the most common one. Besides the locus of enterocyte effacement and type three secretion system, the high pathogenicity island, OI-122, and type six secretion systems were identified. Substantial strain variation in virulence gene repertoire, Shiga toxin production, and cytotoxicity were revealed. Six avian strains exhibited significantly higher cytotoxicity than that of stx2f-positive E. coli, and three of them exhibited a comparable level of cytotoxicity with that of enterohemorrhagic E. coli outbreak strains, suggesting that wild birds could serve as a reservoir of E. albertii strains with great potential to cause severe diseases in humans.

Rapid identification of enteric bacteria from whole genome sequences using average nucleotide identity metrics.

Identification of enteric bacteria species by whole genome sequence (WGS) analysis requires a rapid and an easily standardized approach. We leveraged the principles of average nucleotide identity using MUMmer (ANIm) software, which calculates the percent bases aligned between two bacterial genomes and their corresponding ANI values, to set threshold values for determining species consistent with the conventional identification methods of known species. The performance of species identification was evaluated using two datasets: the Reference Genome Dataset v2 (RGDv2), consisting of 43 enteric genome assemblies representing 32 species, and the Test Genome Dataset (TGDv1), comprising 454 genome assemblies which is designed to represent all species needed to query for identification, as well as rare and closely related species. The RGDv2 contains six Campylobacter spp., three Escherichia/Shigella spp., one Grimontia hollisae, six Listeria spp., one Photobacterium damselae, two Salmonella spp., and thirteen Vibrio spp., while the TGDv1 contains 454 enteric bacterial genomes representing 42 different species. The analysis showed that, when a standard minimum of 70% genome bases alignment existed, the ANI threshold values determined for these species were ≥95 for Escherichia/Shigella and Vibrio species, ≥93% for Salmonella species, and ≥92% for Campylobacter and Listeria species. Using these metrics, the RGDv2 accurately classified all validation strains in TGDv1 at the species level, which is consistent with the classification based on previous gold standard methods.

Use of Large-Scale Genomics to Identify the Role of Animals and Foods as Potential Sources of Extraintestinal Pathogenic Escherichia coli That Cause Human Illness

Extraintestinal pathogenic Escherichia coli (ExPEC) cause urinary tract and potentially life-threatening invasive infections. Unfortunately, the origins of ExPEC are not always clear. We used genomic data of E. coli isolates from five U.S. government organizations to evaluate potential sources of ExPEC infections. Virulence gene analysis of 38,032 isolates from human, food animal, retail meat, and companion animals classified the subset of 8142 non-diarrheagenic isolates into 40 virulence groups. Groups were identified as low, medium, and high relative risk of containing ExPEC strains, based on the proportion of isolates recovered from humans. Medium and high relative risk groups showed a greater representation of sequence types associated with human disease, including ST-131. Over 90% of food source isolates belonged to low relative risk groups, while >60% of companion animal isolates belonged to medium or high relative risk groups. Additionally, 18 of the 26 most prevalent antimicrobial resistance determinants were more common in high relative risk groups. The associations between antimicrobial resistance and virulence potentially limit treatment options for human ExPEC infections. This study demonstrates the power of large-scale genomics to assess potential sources of ExPEC strains and highlights the importance of a One Health approach to identify and manage these human pathogens.

Longitudinal surveillance and comparative characterization of Escherichia albertii in wild raccoons in the United States.

Escherichia albertii is an emerging enteric bacterial pathogen causing watery diarrhea, abdominal distension, vomiting and fever in humans. E. albertii has caused many foodborne outbreaks in Japan and was also reported in other countries worldwide. However, the important animal reservoirs of this pathogen are still largely unknown, impeding us to combat this emerging pathogen. Recently, we reported that wild raccoons (Procyon lotor) and broiler chickens are significant reservoirs of E. albertii in Japan and the U.S., respectively. Here, we performed a longitudinal surveillance to monitor prevalence of E. albertii in wild raccoons in the U.S. and conducted comprehensive comparative analyses of the E. albertii of different origins. A total of 289 fecal swab samples were collected from wild raccoons in Tennessee and Kentucky in the U.S. (2018-2020). Approximately 26% (74/289) of the raccoons examined were PCR-positive for E. albertii and eventually 22 E. albertii isolates were obtained. PFGE analysis showed the U.S. raccoon E. albertii were phylogenetically distant even though the corresponding raccoons were captured from a small area. Unlike the high prevalence of multidrug resistance (83%) observed in previous chicken E. albertii survey, antibiotic resistance was rarely observed in all the U.S. raccoon and 22 Japan raccoon strains with only one Japan strain displaying multidrug resistance (2%). Whole genome sequencing of 54 diverse E. albertii strains and subsequent comparative genomics analysis revealed unique clusters that displayed close evolutionary relationships and similar virulence gene profiles among the strains of different origins in terms of geographical locations (e.g., U.S. and Japan) and hosts (raccoon, chicken, swine, and human). Challenge experiment demonstrated raccoon E. albertii strains could successfully colonize in the chicken intestine at 3 and 8 days postinfection. A pilot environmental survey further showed all the four tested water samples from Tennessee river were E. albertii-positive; two different E. albertii strains, isolated from a single water sample, showed close relationships to those of human origin. Together, the findings from this study provide new insights into the ecology, evolution, and pathobiology of E. albertii, and underscore the need to control the emerging E. albertii in a complex ecosystem using One Health approach.

Use of Whole Genome Sequencing by the Federal Interagency Collaboration for Genomics for Food and Feed Safety in the United States.

ABSTRACT: This multiagency report developed by the Interagency Collaboration for Genomics for Food and Feed Safety provides an overview of the use of and transition to whole genome sequencing (WGS) technology for detection and characterization of pathogens transmitted commonly by food and for identification of their sources. We describe foodborne pathogen analysis, investigation, and harmonization efforts among the following federal agencies: National Institutes of Health; Department of Health and Human Services, Centers for Disease Control and Prevention (CDC) and U.S. Food and Drug Administration (FDA); and the U.S. Department of Agriculture, Food Safety and Inspection Service, Agricultural Research Service, and Animal and Plant Health Inspection Service. We describe single nucleotide polymorphism, core-genome, and whole genome multilocus sequence typing data analysis methods as used in the PulseNet (CDC) and GenomeTrakr (FDA) networks, underscoring the complementary nature of the results for linking genetically related foodborne pathogens during outbreak investigations while allowing flexibility to meet the specific needs of Interagency Collaboration partners. We highlight how we apply WGS to pathogen characterization (virulence and antimicrobial resistance profiles) and source attribution efforts and increase transparency by making the sequences and other data publicly available through the National Center for Biotechnology Information. We also highlight the impact of current trends in the use of culture-independent diagnostic tests for human diagnostic testing on analytical approaches related to food safety and what is next for the use of WGS in the area of food safety.

Microarray-based analysis of IncA/C plasmid-associated genes from multidrug-resistant Salmonella enterica.

In the family Enterobacteriaceae, plasmids have been classified according to 27 incompatibility (Inc) or replicon types that are based on the inability of different plasmids with the same replication mechanism to coexist in the same cell. Certain replicon types such as IncA/C are associated with multidrug resistance (MDR). We developed a microarray that contains 286 unique 70-mer oligonucleotide probes based on sequences from five IncA/C plasmids: pYR1 (Yersinia ruckeri), pPIP1202 (Yersinia pestis), pP99-018 (Photobacterium damselae), pSN254 (Salmonella enterica serovar Newport), and pP91278 (Photobacterium damselae). DNA from 59 Salmonella enterica isolates was hybridized to the microarray and analyzed for the presence or absence of genes. These isolates represented 17 serovars from 14 different animal hosts and from different geographical regions in the United States. Qualitative cluster analysis was performed using CLUSTER 3.0 to group microarray hybridization results. We found that IncA/C plasmids occurred in two lineages distinguished by a major insertion-deletion (indel) region that contains genes encoding mostly hypothetical proteins. The most variable genes were represented by transposon-associated genes as well as four antimicrobial resistance genes (aphA, merP, merA, and aadA). Sixteen mercury resistance genes were identified and highly conserved, suggesting that mercury ion-related exposure is a stronger pressure than anticipated. We used these data to construct a core IncA/C genome and an accessory genome. The results of our studies suggest that the transfer of antimicrobial resistance determinants by transfer of IncA/C plasmids is somewhat less common than exchange within the plasmids orchestrated by transposable elements, such as transposons, integrating and conjugative elements (ICEs), and insertion sequence common regions (ISCRs), and thus pose less opportunity for exchange of antimicrobial resistance.

Sequence analysis of a group of low molecular-weight plasmids carrying multiple IS903 elements flanking a kanamycin resistance aph gene in Salmonella enterica serovars.

A group of low molecular-weight ColE1-like plasmids carrying the aph sequence type aph(ii) from three different Salmonella serovars were sequenced. These plasmids carry two or more copies of IS903 elements, with up to 21bp sequence differences to one another, two of which flank the aph gene. This group of plasmids did not appear to carry any known mobilization genes and instead carry three open reading frames encoding hypothetical proteins of unknown function possibly organized in an operon. The plasmid replication region (RNA I/II--rom) of this plasmid group showed extensive homology to that of pKPN2 plasmid of Klebsiella pneumoniae and pCol-let plasmid of Escherichia coli. Three of the four plasmids had identical sequences, and the fourth had an extra copy of IS903 with target duplication, suggesting a recent divergence in the different Salmonella serovars from a common ancestor.

Related antimicrobial resistance genes detected in different bacterial species co-isolated from swine fecal samples.

A potential factor leading to the spread of antimicrobial resistance (AR) in bacteria is the horizontal transfer of resistance genes between bacteria in animals or their environment. To investigate this, swine fecal samples were collected on-farm and cultured for Escherichia coli, Salmonella enterica, Campylobacter spp., and Enterococcus spp. which are all commonly found in swine. Forty-nine of the samples from which all four bacteria were recovered were selected yielding a total of 196 isolates for analysis. Isolates were tested for antimicrobial susceptibility followed by hybridization to a DNA microarray designed to detect 775 AR-related genes. E. coli and Salmonella isolated from the same fecal sample had the most AR genes in common among the four bacteria. Genes detected encoded resistance to aminoglycosides (aac(3), aadA1, aadB, and strAB), ß-lactams (ampC, ampR, and bla(TEM)), chloramphenicols (cat and floR), sulfanillic acid (sul1/sulI), tetracyclines (tet(A), tet(D), tet(C), tet(G), and tet(R)), and trimethoprim (dfrA1 and dfh). Campylobacter coli and Enterococcus isolated from the same sample frequently had tet(O) and aphA-3 genes detected in common. Almost half (47%) of E. coli and Salmonella isolated from the same fecal sample shared resistance genes at a significant level (χ², p < 0.0000001). These data suggest that there may have been horizontal exchange of AR genes between these bacteria or there may be a common source of AR genes in the swine environment for E. coli and Salmonella.

Prevalence of ColE1-like plasmids and kanamycin resistance genes in Salmonella enterica serovars.

Multi-antimicrobial-resistant Salmonella enterica strains frequently carry resistance genes on plasmids. Recent studies focus heavily on large conjugative plasmids, and the role that small plasmids play in resistance gene transfer is largely unknown. To expand our previous studies in assessing the prevalence of the isolates harboring ColE1-like plasmids carrying the aph gene responsible for kanamycin resistance (Kan(r)) phenotypes, 102 Kan(r) Salmonella isolates collected through the National Antimicrobial Resistance Monitoring System (NARMS) in 2005 were screened by PCR using ColE1 primer sets. Thirty isolates were found to be positive for ColE1-like replicon. Plasmids from 23 isolates were able to propagate in Escherichia coli and were subjected to further characterization. Restriction mapping revealed three major plasmid groups found in three or more isolates, with each group consisting of two to three subtypes. The aph genes from the Kan(r) Salmonella isolates were amplified by PCR, sequenced, and showed four different aph(3')-I genes. The distribution of the ColE1 plasmid groups in association with the aph gene, Salmonella serovar, and isolate source demonstrated a strong linkage of the plasmid with S. enterica serovar Typhimurium DT104. Due to their high copy number and mobility, the ColE1-like plasmids may play a critical role in transmission of antibiotic resistance genes among enteric pathogens, and these findings warrant a close monitoring of this plasmid incompatibility group.

The Virulence of Escherichia coli O157:H7 Isolates in Mice Depends on Shiga Toxin Type 2a (Stx2a)-Induction and High Levels of Stx2a in Stool.

In this study we compared nine Shiga toxin (Stx)-producing Escherichia coli O157:H7 patient isolates for Stx levels, stx-phage insertion site(s), and pathogenicity in a streptomycin (Str)-treated mouse model. The strains encoded stx2a, stx1a and stx2a, or stx2a and stx2c. All of the strains elaborated 105-106 cytotoxic doses 50% (CD50) into the supernatant after growth in vitro as measured on Vero cells, and showed variable levels of increased toxin production after growth with sub-inhibitory levels of ciprofloxacin (Cip). The stx2a+stx2c+ isolates were 90-100% lethal for Str-treated BALB/c mice, though one isolate, JH2013, had a delayed time-to-death. The stx2a+ isolate was avirulent. Both an stx2a and a recA deletion mutant of one of the stx2a+stx2c+ strains, JH2010, exhibited at least a three-log decrease in cytotoxicity in vitro and both were avirulent in the mice. Stool from Str-treated mice infected with the highly virulent isolates were 10- to 100-fold more cytotoxic than feces from mice infected with the clinical isolate, JH2012, that made only Stx2a. Taken together these findings demonstrate that the stx2a-phage from JH2010 induces to higher levels in vivo than does the phage from JH2012. The stx1a+stx2a+ clinical isolates were avirulent and neutralization of Stx1 in stool from mice infected with those strains indicated that the toxin produced in vivo was primarily Stx1a. Treatment of mice infected with Stx1a+Stx2a+ isolates with Cip resulted in an increase in Stx2a production in vivo and lethality in the mice. Our data suggest that high levels of Stx2a in stool are predictive of virulence in mice.

Inc A/C plasmids are prevalent in multidrug-resistant Salmonella enterica isolates.

Bacterial plasmids are fragments of extrachromosomal double-stranded DNA that can contain a variety of genes that are beneficial to the host organism, like those responsible for antimicrobial resistance. The objective of this study was to characterize a collection of 437 Salmonella enterica isolates from different animal sources for their antimicrobial resistance phenotypes and plasmid replicon types and, in some cases, by pulsed-field gel electrophoresis (PFGE) in an effort to learn more about the distribution of multidrug resistance in relation to replicon types. A PCR-based replicon typing assay consisting of three multiplex PCRs was used to detect 18 of the 26 known plasmid types in the Enterobacteriaceae based on their incompatibility (Inc) replicon types. Linkage analysis was completed with antibiograms, replicon types, serovars, and Inc A/C. Inc A/C plasmids were prevalent in multidrug-resistant isolates with the notable exception of Salmonella enterica serovar Typhimurium. Cluster analysis based on PFGE of a subset of 216 isolates showed 155 unique types, suggesting a variable population, but distinct clusters of isolates with Inc A/C plasmids were apparent. Significant linkage of serovar was also seen with Inc replicon types B/O, I1, Frep, and HI1. The present study showed that the combination of Salmonella, the Inc A/C plasmids, and multiple-drug-resistant genes is very old. Our results suggest that some strains, notably serovar Typhimurium and closely related types, may have once carried the plasmid but that the resistance genes were transferred to the chromosome with the subsequent loss of the plasmid.

PacBio Genome Sequences of Eight Escherichia albertii Strains Isolated from Humans in the United States.

Escherichia albertii is an emerging pathogen that is closely related to Escherichia coli and can carry some of the same virulence genes as E. coli. Here, we report the release of Illumina-corrected PacBio sequences for eight E. albertii genomes. Two of these strains carry Shiga toxin 2f.

Isolation and characterization of two novel groups of kanamycin-resistance ColE1-like plasmids in Salmonella enterica serotypes from food animals.

While antimicrobial resistance in Salmonella enterica is mainly attributed to large plasmids, small plasmids may also harbor antimicrobial resistance genes. Previously, three major groups of ColE1-like plasmids conferring kanamycin-resistance (KanR) in various S. enterica serotypes from diagnostic samples of human or animals were reported. In this study, over 200 KanR S. enterica isolates from slaughter samples, collected in 2010 and 2011 as a part of the animal arm of the National Antimicrobial Resistance Monitoring System, were screened for the presence of ColE1-like plasmids. Twenty-three KanR ColE1-like plasmids were successfully isolated. Restriction fragment mapping revealed five major plasmid groups with subgroups, including two new groups, X (n = 3) and Y/Y2/Y3 (n = 4), in addition to the previously identified groups A (n = 7), B (n = 6), and C/C3 (n = 3). Nearly 75% of the plasmid-carrying isolates were from turkey and included all the isolates carrying X and Y plasmids. All group X plasmids were from serotype Hadar. Serotype Senftenberg carried all the group Y plasmids and one group B plasmid. All Typhimurium isolates (n = 4) carried group A plasmids, while Newport isolates (n = 3) each carried a different plasmid group (A, B, or C). The presence of the selection bias in the NARMS strain collection prevents interpretation of findings at the population level. However, this study demonstrated that KanR ColE1-like plasmids are widely distributed among different S. enterica serotypes in the NARMS isolates and may play a role in dissemination of antimicrobial resistance genes.

High-Quality Whole-Genome Sequences for 59 Historical Shigella Strains Generated with PacBio Sequencing.

Shigella spp. are enteric pathogens that cause shigellosis. We report here the high-quality whole-genome sequences of 59 historical Shigella strains that represent the four species and a variety of serotypes.

High-Quality Whole-Genome Sequences for 77 Shiga Toxin-Producing Escherichia coli Strains Generated with PacBio Sequencing.

Shiga toxin-producing Escherichia coli (STEC) is an enteric foodborne pathogen that can cause mild to severe illness. Here, we report the availability of high-quality whole-genome sequences for 77 STEC strains generated using the PacBio sequencing platform.