Virulence Gene Profiles and Clonal Relationships of Escherichia coli O26:H11 Isolates from Feedlot Cattle as Determined by Whole-Genome Sequencing

Developing the PIP-eco: An integrated genomic pipeline for identification and characterization of Escherichia coli pathotypes encompassing hybrid forms

Pathogenic Escherichia coli (E. coli) strains are distinguished by their diverse virulence factors, which contribute to a wide spectrum of diseases. These pathogens evolve through the horizontal transfer of virulence factors, resulting in the emergence of hybrid pathotypes with complex and heterogeneous characteristics. Recognizing their profound impact on public health, this study introduces the PIP-eco pipeline, a comprehensive analytical tool designed for the precise identification and characterization of E. coli pathotypes. This PIP-eco pipeline advances beyond traditional molecular techniques by facilitating detailed analysis of both single and hybrid pathotypes. It integrates targeted marker gene analysis, virulence factor-based phylogenetic analysis, and pathogenicity islands (PAIs) profiling to elucidate the genetic diversity of E. coli pathotypes and support their accurate classification. This integrative approach enables PIP-eco to uncover connections among various E. coli pathotypes, highlight shared virulence factors, and provide insights into their evolutionary trajectories. By utilizing experimentally validated marker genes, the pipeline ensures robust identification of pathotypes, particularly those of hybrid pathotypes. Additionally, PAI analysis offers comprehensive genetic investigations, revealing strain-specific variations and potential virulence mechanisms. As a result, the PIP-eco pipeline emerges as a useful tool for dissecting the evolutionary dynamics of E. coli and characterizing complex pathotypes, addressing the critical need for accurate detection and understanding of hybrid pathotypes.

Epidemiology and antimicrobial resistance profiles of pathogenic Escherichia coli from commercial swine and poultry abattoirs and farms in South Africa: A One Health approach

Pathogenic Escherichia coli (PEC) are important foodborne bacteria that can cause severe illness in humans. The PECs thrive within the intestines of humans as well as animals and may contaminate multiple ecosystems, including food and water, via faecal transmission. Abattoir and farm employees are at high risk of PEC exposure, which could translate to community risk through person-to-person contact. To determine the epidemiology and resistome of PECs in Gauteng and Limpopo provinces of South Africa, 198 swine faecal samples, 220 poultry cloacal swabs, 108 human hand swabs, 11 run-off water samples from abattoirs and farms were collected from four swine and five poultry commercial abattoirs and two swine farms. One effluent sample each was collected from four wastewater treatment plants (WWTP) and a tertiary hospital setting. Phenotypic and genotypic techniques were used including polymerase chain reaction, pulsed-field gel electrophoresis (PFGE) and whole genome sequencing (WGS). Results showed EHEC and EPEC prevalence was 4.1 % (22/542) and 20.8 % (113/542), respectively, with the O26 serogroup detected the most in PEC isolates. According to the PFGE dendrogram, isolates from poultry, human hand swabs and run-off water clustered together. Diverse virulence factors such as the novel stx2k subtype and eae genes were detected among the 36 representative PEC isolates according to WGS. The results showed that 66.7 % (24/36) of sequenced PECs presented with multi-drug resistance (MDR) to β-lactamase 13.9 % (5/36), aminoglycoside 61.1 % (22/36), tetracycline 41.7 % (15/36) and quinolones 38.9 % (14/36). No colistin nor carbapenem resistance was detected. Sequence types (STs) associated with MDR in this study were: ST752, ST189, ST206, ST10, ST48 and ST38. The findings highlight the threat of zoonotic pathogens to close human contacts and the need for enhanced surveillance to mitigate the spread of MDR foodborne PECs.

Phylogenetic relationship and virulence composition of Escherichia coli O26:H11 cattle and human strain collections in Scotland; 2002-2020

O26 is the commonest non-O157 Shiga toxin (stx)-producing Escherichia coli serogroup reported in human infections worldwide. Ruminants, particularly cattle, are the primary reservoir source for human infection. In this study, we compared the whole genomes and virulence profiles of O26:H11 strains (n = 99) isolated from Scottish cattle with strains from human infections (n = 96) held by the Scottish Escherichia coli O157/STEC Reference Laboratory, isolated between 2002 and 2020. Bovine strains were from two national cross-sectional cattle surveys conducted between 2002-2004 and 2014-2015. A maximum likelihood phylogeny was constructed from a core-genome alignment with the O26:H11 strain 11368 reference genome. Genomes were screened against a panel of 2,710 virulence genes using the Virulence Finder Database. All stx-positive bovine O26:H11 strains belonged to the ST21 lineage and were grouped into three main clades. Bovine and human source strains were interspersed, and the stx subtype was relatively clade-specific. Highly pathogenic stx2a-only ST21 strains were identified in two herds sampled in the second cattle survey and in human clinical infections from 2010 onwards. The closest pairwise distance was 9 single-nucleotide polymorphisms (SNPs) between Scottish bovine and human strains and 69 SNPs between the two cattle surveys. Bovine O26:H11 was compared to public EnteroBase ST29 complex genomes and found to have the greatest commonality with O26:H11 strains from the rest of the UK, followed by France, Italy, and Belgium. Virulence profiles of stx-positive bovine and human strains were similar but more conserved for the stx2a subtype. O26:H11 stx-negative ST29 (n = 17) and ST396 strains (n = 5) were isolated from 19 cattle herds; all were eae-positive, and 10 of these herds yielded strains positive for ehxA, espK, and Z2098, gene markers suggestive of enterohaemorrhagic potential. There was a significant association (p < 0.001) between nucleotide sequence percent identity and stx status for the bacteriophage insertion site genes yecE for stx2 and yehV for stx1. Acquired antimicrobial resistance genes were identified in silico in 12.1% of bovine and 17.7% of human O26:H11 strains, with sul2, tet, aph(3″), and aph(6″) being most common. This study describes the diversity among Scottish bovine O26:H11 strains and investigates their relationship to human STEC infections.

Large-Scale Phylogenetic Analysis Reveals a New Genetic Clade among Escherichia coli O26 Strains

Shiga toxin-producing Escherichia coli (STEC) O26 is the predominant non-O157 serogroup causing hemolytic uremic syndrome worldwide. Moreover, the serogroup is highly dynamic and harbors several pathogenic clones. Here, we investigated the phylogenetic relationship of STEC O26 at a global level based on 1,367 strains from 20 countries deposited in NCBI and Enterobase databases. The whole-genome-based analysis identified a new genetic clade, called ST29C4. The new clade was unique in terms of multilocus sequence type (ST29), CRISPR (group Ia), and dominant plasmid gene profile (ehxA+/katP-/espP-/etpD-). Moreover, the combination of multiple typing methods (core genome single nucleotide polymorphism [SNP] typing, CRISPR typing, and virulence genes analysis) demonstrated that this new lineage ST29C4 was in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains. Besides, we observed that ST29C4 harbored extraintestinal pathogenic E. coli (ExPEC)-related virulence gene (VG), tsh, and STEC-associated VG, stx2a, suggesting the emergence of a hybrid pathogen. The ST29C4 strains also exhibited high similarity in stx2a-prophage and integrase with the O104:H4 strain, further demonstrating its potential risk to human health. Collectively, the large-scale phylogenetic analysis extends the understanding of the clonal structure of O26 strains and provides new insights for O26 strain microevolution.

IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) O26 is the second prevalent STEC serogroup only to O157, which can cause a series of diseases ranging from mild diarrhea to life-threatening hemolytic uremic syndrome (HUS). The serogroup is highly diverse and multiple clones are characterized, including ST29C1-C3 and ST21C1-C2. However, the phylogenetic relationship of these clones remains fully unclear. In this study, we revealed a new genetic clade among O26 strains, ST29C4, which was unique in terms of CRISPR, multilocus sequence type (MLST), and plasmid gene profile (PGP). Moreover, the combination of multiple typing methods demonstrated that this new clone was located in the intermediate phylogenetic position between ST29C3 and other non-ST29C3 strains (i.e., ST29C1-C2 and ST21C1-C2). Overall, the large-scale phylogenetic analysis extends our current understanding of O26 microevolution.

SCRAMBLER: A Tool for De Novo CRISPR Array Reconstruction and Its Application for Analysis of the Structure of Prokaryotic Populations

CRISPR arrays are prokaryotic genomic loci consisting of repeat sequences alternating with unique spacers acquired from foreign nucleic acids. As one of the fastest-evolving parts of the genome, CRISPR arrays can be used to differentiate closely related prokaryotic lineages and track individual strains in prokaryotic communities. However, the assembly of full-length CRISPR arrays sequences remains a problem. Here, we developed SCRAMBLER, a tool that includes several pipelines for assembling CRISPR arrays from high-throughput short-read sequencing data. We assessed its performance with model data sets (Escherichia coli strains containing different CRISPR arrays and imitating prokaryotic communities of different complexities) and intestinal microbiomes of extant and extinct pachyderms. Evaluation of SCRAMBLER's performance using model data sets demonstrated its ability to assemble CRISPR arrays correctly from reads containing pairs of spacers, yielding a precision rate of >80% and a recall rate of 60-85% when checked against ground-truth data. Likewise, SCRAMBLER successfully assembled CRISPR arrays from the environmental samples, as attested by their matching with database entries. SCRAMBLER, an open-source software (github.com/biolab-tools/SCRAMBLER), can facilitate analysis of the composition and dynamics of CRISPR arrays in complex communities.

Inconsistent PCR detection of Shiga toxin-producing Escherichia coli: Insights from whole genome sequence analyses

Shiga toxin-producing Escherichia coli (STEC) have been linked to food-borne disease outbreaks. As PCR is routinely used to screen foods for STEC, it is important that factors leading to inconsistent detection of STEC by PCR are understood. This study used whole genome sequencing (WGS) to investigate causes of inconsistent PCR detection of stx1, stx2, and serogroup-specific genes. Fifty strains isolated from Alberta feedlot cattle from three different studies were selected with inconsistent or consistent detection of stx and serogroup by PCR. All isolates were initially classified as STEC by PCR. Sequencing was performed using Illumina MiSeq® with sample library by Nextera XT. Virtual PCRs were performed using Geneious and bacteriophage content was determined using PHASTER. Sequencing coverage ranged from 47 to 102x, averaging 74x, with sequences deposited in the NCBI database. Eleven strains were confirmed by WGS as STEC having complete stxA and stxB subunits. However, truncated stx fragments occurred in twenty-two other isolates, some having multiple stx fragments in the genome. Isolates with complete stx by WGS had consistent stx1 and stx2 detection by PCR, although one also having a stx2 fragment had inconsistent stx2 PCR. For all STEC and 18/39 non-STEC, serogroups determined by PCR agreed with those determined by WGS. An additional three WGS serotypes were inconclusive and two isolates were Citrobacter spp. Results demonstrate that stx fragments associated with stx-carrying bacteriophages in the E. coli genome may contribute to inconsistent detection of stx1 and stx2 by PCR. Fourteen isolates had integrated stx bacteriophage but lacked complete or fragmentary stx possibly due to partial bacteriophage excision after sub-cultivation or other unclear mechanisms. The majority of STEC isolates (7/11) did not have identifiable bacteriophage DNA in the contig(s) where stx was located, likely increasing the stability of stx in the bacterial genome and its detection by PCR.

Diversity of Non-O157 Shiga Toxin-Producing Escherichia coli Isolated from Cattle from Central and Southern Chile

Simple Summary

Cattle are the main reservoir of Shiga toxin-producing E. coli (STEC), foodborne pathogens that cause severe disease and outbreaks. However, not all STEC cause human illnesses or have the same virulence potential. Characterizing strains isolated worldwide allows insights into how strains spread and which isolates have a more significant risk potential. This study described STEC isolation rates from cattle in Chile and characterized 30 isolates. We obtained 93 STEC isolates from 56/446 (12.6%) fecal cattle samples. Then, 30 non-O157 STEC isolates were selected for complete characterization; we found isolates of 16 different sequence types and 17 serotypes. One isolate was resistant to tetracycline and carried resistance genes against the drug. Surveyed virulence genes (n = 31) were present from 13% to 100% of isolates, and one isolate carried 26/31 virulence genes. Most isolates (90%; 27/30) carried the stx2 gene, which is frequently linked to strains causing severe disease. A phylogenetic reconstruction demonstrated that isolates clustered based on serotypes, independent of their geographical origin (Central or Southern Chile). These results indicate that cattle in Chile carry a wide diversity of STEC potentially pathogenic for humans based on the presence of virulence genes.

Abstract

Cattle are the main reservoir of Shiga toxin-producing Escherichia coli (STEC), one of the world’s most important foodborne pathogens. The pathogen causes severe human diseases and outbreaks. This study aimed to identify and characterize non-O157 STEC isolated from cattle feces from central and southern Chile. We analyzed 446 cattle fecal samples and isolated non-O157 STEC from 12.6% (56/446); a total of 93 different isolates were recovered. Most isolates displayed β-glucuronidase activity (96.8%; 90/93) and fermented sorbitol (86.0%; 80/93), whereas only 39.8% (37/93) were resistant to tellurite. A subgroup of 30 representative non-O157 STEC isolates was selected for whole-genome sequencing and bioinformatics analysis. In silico analysis showed that they grouped into 16 different sequence types and 17 serotypes; the serotypes most frequently identified were O116:H21 and O168:H8 (13% each). A single isolate of serotype O26:H11 was recovered. One isolate was resistant to tetracycline and carried resistance genes tet(A) and tet(R); no other isolate displayed antimicrobial resistance or carried antimicrobial resistance genes. The intimin gene (eae) was identified in 13.3% (4/30) of the genomes and 90% (27/30) carried the stx2 gene. A phylogenetic reconstruction demonstrated that the isolates clustered based on serotypes, independent of geographical origin. These results indicate that cattle in Chile carry a wide diversity of STEC potentially pathogenic for humans based on the presence of critical virulence genes.

Characterization of Shiga toxin-producing Escherichia coli isolated from Cattle and Sheep in Xinjiang province, China, using whole-genome sequencing

Shiga toxin-producing Escherichia coli (STEC) is an important food-borne pathogen capable of causing severe gastrointestinal diseases in humans. Cattle and sheep are the natural reservoir hosts of STEC strains. Previously, we isolated 56 STEC strains from anal and carcass swab samples of cattle and sheep in farms and slaughterhouses. In this study, we performed whole-genome sequencing of these isolates and determined their serotypes, virulence profiles, sequence types (STs) and genetic relationships. Our results showed that the 56 isolates belong to 20 different STs, 29 O:H serotypes and 8 stx subtype combinations. The highly prevalent serotypes for bovine and ovine isolates were O8:H25 and O87:H16, respectively. Five serotypes of cattle or sheep isolates are novel. The majority (63%) of cattle isolates contain stx1 + stx2, subtyped into stx1a, stx2a and stx2c. In contrast, most of the sheep isolates contain stx1 only, primarily subtyped into stx1a and stx1c. None of the isolates tested eae-positive, but virulence factors such as ehxA and espP were present with variable prevalence rates. The prevalence of saa (19.6%) and espP (12.5%) in cattle isolates is much higher than that in sheep isolates, whereas that of subA (34%), katP (14.3%) and ireA (28.6%) in sheep isolates is considerably higher than that in cattle isolates. Core-genome SNP analysis revealed that the majority of isolates could be clustered based on their serotypes or STs, whereas some clustering is associated with more than one ST or serotype. Five sheep isolates (4 belonging to ST675 and serotype O76:H19 and 1 belonging to ST25 and serotype O128:H2) share STs, serotypes and stx profiles with two hemolytic uremic syndrome-associated enterohemorrhagic E. coli (HUSEC) isolates; a cattle isolate belonging to the same ST as HUSEC isolate HUSEC001 contains all the nine virulence genes tested. These data suggest a potential of the six isolates for causing severe human infections. Collectively, we described the characteristics of cattle and sheep STEC isolates from Xinjiang, China, which may be utilized in comparative studies of other geographic regions and sources of isolation, and for surveillance as well.

Precision long-read metagenomics sequencing for food safety by detection and assembly of Shiga toxin-producing Escherichia coli in irrigation water

Shiga toxin-producing Escherichia coli (STEC) contamination of agricultural water might be an important factor to recent foodborne illness and outbreaks involving leafy greens. Closed bacterial genomes from whole genome sequencing play an important role in source tracking. We aimed to determine the limits of detection and classification of STECs by qPCR and nanopore sequencing using 24 hour enriched irrigation water artificially contaminated with E. coli O157:H7 (EDL933). We determined the limit of STEC detection by qPCR to be 30 CFU/reaction, which is equivalent to 105 CFU/ml in the enrichment. By using Oxford Nanopore's EPI2ME WIMP workflow and de novo assembly with Flye followed by taxon classification with a k-mer analysis software (Kraken2), E. coli O157:H7 could be detected at 103 CFU/ml (68 reads) and a complete fragmented E. coli O157:H7 metagenome-assembled genome (MAG) was obtained at 105-108 CFU/ml. Using a custom script to extract the E. coli reads, a completely closed MAG was obtained at 107-108 CFU/ml and a complete, fragmented MAG was obtained at 105-106 CFU/ml. In silico virulence detection for E. coli MAGs for 105-108 CFU/ml showed that the virulotype was indistinguishable from the spiked E. coli O157:H7 strain. We further identified the bacterial species in the un-spiked enrichment, including antimicrobial resistance genes, which could have important implications to food safety. We propose this workflow provides proof of concept for faster detection and complete genomic characterization of STECs from a complex microbial sample compared to current reporting protocols and could be applied to determine the limit of detection and assembly of other foodborne bacterial pathogens.

The Applications of Nanopore Sequencing Technology in Pathogenic Microorganism Detection

Infectious diseases are major threats to human health and lead to a serious public health burden. The emergence of new pathogens and the mutation of known pathogens challenge our ability to diagnose and control infectious diseases. Nanopore sequencing technology exhibited versatile applications in pathogenic microorganism detection due to its flexible data throughput. This review article introduced the applications of nanopore sequencing in clinical microbiology and infectious diseases management, including the monitoring of emerging infectious diseases outbreak, identification of pathogen drug resistance, and disease-related microbial communities characterization.

High-Zinc Supplementation of Weaned Piglets Affects Frequencies of Virulence and Bacteriocin Associated Genes Among Intestinal Escherichia coli Populations

To prevent economic losses due to post-weaning diarrhea (PWD) in industrial pig production, zinc (Zn) feed additives have been widely used, especially since awareness has risen that the regular application of antibiotics promotes buildup of antimicrobial resistance in both commensal and pathogenic bacteria. In a previous study on 179 Escherichia coli collected from piglets sacrificed at the end of a Zn feeding trial, including isolates obtained from animals of a high-zinc fed group (HZG) and a corresponding control group (CG), we found that the isolate collection exhibited three different levels of tolerance toward zinc, i.e., the minimal inhibitory concentration (MIC) detected was 128, followed by 256 and 512 μg/ml ZnCl2. We further provided evidence that enhanced zinc tolerance in porcine intestinal E. coli populations is clearly linked to excessive zinc feeding. Here we provide insights about the genomic make-up and phylogenetic background of these 179 E. coli genomes. Bayesian analysis of the population structure (BAPS) revealed a lack of association between the actual zinc tolerance level and a particular phylogenetic E. coli cluster or even branch for both, isolates belonging to the HZG and CG. In addition, detection rates for genes and operons associated with virulence (VAG) and bacteriocins (BAG) were lower in isolates originating from the HZG (41 vs. 65% and 22 vs. 35%, p < 0.001 and p = 0.002, resp.). Strikingly, E. coli harboring genes defining distinct pathotypes associated with intestinal disease, i.e., enterotoxigenic, enteropathogenic, and Shiga toxin-producing E. coli (ETEC, EPEC, and STEC) constituted 1% of the isolates belonging to the HZG but 14% of those from the CG. Notably, these pathotypes were positively associated with enhanced zinc tolerance (512 μg/ml ZnCl2 MIC, p < 0.001). Taken together, zinc excess seems to influence carriage rates of VAGs and BAGs in porcine intestinal E. coli populations, and high-zinc feeding is negatively correlated with enteral pathotype occurrences, which might explain earlier observations concerning the relative increase of Enterobacterales considering the overall intestinal microbiota of piglets during zinc feeding trials while PWD rates have decreased.

Occurrence and genomic characterization of ESBL-producing Escherichia coli ST29 strains from swine with abundant virulence genes

Food-production animals were considered to be a major reservoir of antimicrobial-resistant bacteria and clinically relevant pathogens. The potential of commensal Escherichia coli from pigs as a source of opportunistic pathogens associated with extraintestinal infections in humans needs to be assessed. In this study, 13 E. coli isolates from an intensive pig farm in China were analyzed using whole genome sequencing followed by in-depth in silico analysis. Genomic analysis showed comprehensive antimicrobial resistance profiles, with each isolate carrying between 4 and 22 antimicrobial resistance genes. Although these E. coli isolates were assigned to low-virulence phylogroup A and B1, 31 different virulence genes were detected at least once in the 13 sequenced isolates. Extraintestinal pathogenic E. coli-associated virulence genes, including iss, iha, tsh and iroN, were found in commensal E. coli isolates in this study. Of note, a large number of virulence genes (n = 22) were identified in ESBL-producing E. coli sequence type (ST) 29 isolates. Our study revealed the presence of comprehensive antimicrobial resistance and virulence gene profiles in commensal E. coli isolates of pigs. The emerged ESBL-producing E. coli ST 29 isolates harboring a high abundance of VAGs highlighted that this new clonal linage may pose a threat to public health.

Proteomic analysis and optimized production of Alkalihalobacillus patagoniensis PAT 05T extracellular proteases

Extracellular proteolytic extracts from the haloalkalitolerant strain Alkalihalobacillus patagoniensis PAT 05^T have proved highly efficient to reduce wool felting, as part of an ecofriendly treatment suitable for organic wool. In the present study, we identified the extracellular proteases produced by PAT 05^T and we optimized its growth conditions for protease production through statistical methods. A total of 191 proteins were identified in PAT 05^T culture supernatants through mass spectrometry analysis. Three of the 6 detected extracellular proteases belonged to the serine-endopeptidase family S8 (EC 3.4.21); two of them showed 86.3 and 67.9% identity with an alkaline protease from Bacillus alcalophilus and another one showed 50.4% identity with Bacillopeptidase F. The other 3 proteases exhibited 55.3, 49.4 and 61.1% identity with D-alanyl-D-alanine carboxypeptidase DacF, D-alanyl-D-alanine carboxypeptidase DacC and endopeptidase LytE, respectively. Using a Fractional Factorial Design followed by a Central Composite Design optimization, a twofold increase in protease production was reached. NaCl concentration was the most influential factor on protease production. The usefulness of PAT 05^T extracellular proteolytic extracts to reduce wool felting was possible associated with the activity of the serine-endopeptidases closely related to highly alkaline keratinolytic proteases. The other identified proteases could cooperate, improving protein hydrolysis. This study provided valuable information for the exploitation of PAT 05^T proteases which have potential for the valorization of organic wool as well as for other industrial applications.

Preliminary data on the antimicrobial effect of Cannabis sativa L. variety Futura 75 against food-borne pathogens in vitro as well as against naturally occurring microbial populations on minced meat during storage

In the present study, the antimicrobial effect of Cannabis sativa Futura 75 was evaluated both in vitro against foodborne bacterial pathogens, and on food against naturally occurring microbial groups of minced meat stored for 8 days at 4°C. Ethanol extraction was performed on the grind of the inflorescence. After extraction, ethanol was completely evaporated and substituted by water. Serial dilutions of the extract, the grind and cannabidiol 99% were added to Nutrient Agar and spotted with Listeria monocytogenes, Salmonella Typhimurium, Escherichia coli and Staphylococcus spp. Regarding the evaluation on food, 50 mL of extract, characterised by CBD at concentration of 322,70 μg/mL, were added to 2.5 kg of minced beef meat. Meat was divided into aliquots and stored for 8 days at 4°C. At 0, 1, 2, 3, 4, 7, and 8 days, aerobic bacteria, enterobacteria, coliforms and E. coli were enumerated. All tested products were efficient against Gram +. In particular, extract corresponding to CBD concentration of 0.017 and 0.3 mg/mL were effective against L. monocytogenes and Staphylococcus spp. respectively. After 8 days of storage at 4°C, treated minced meat showed a bright red colour in comparison to a brownish control meat. Moreover, Enterobacteriaceae and coliforms were significantly reduced of 2.3 log CFU/g and 1.6 log CFU/g respectively in treated meat in comparison to the control. Although preliminary, the present study suggests the antimicrobial properties of the extract of Cannabis sativa both in vitro and in minced meat.

Nanopore sequencing for fast determination of plasmids, phages, virulence markers, and antimicrobial resistance genes in Shiga toxin-producing Escherichia coli

Whole genome sequencing can provide essential public health information. However, it is now known that widely used short-read methods have the potential to miss some randomly-distributed segments of genomes. This can prevent phages, plasmids, and virulence factors from being detected or properly identified. Here, we compared assemblies of three complete Shiga toxin-producing Escherichia coli (STEC) O26:H11/H- genomes from two different sequence types (ST21 and 29), each acquired using the Nextera XT MiSeq, MinION nanopore-based sequencing, and Pacific Biosciences (PacBio) sequencing. Each closed genome consisted of a single chromosome, approximately 5.7 Mb for CFSAN027343, 5.6 Mb for CFSAN027346, and 5.4 MB for CFSAN027350. However, short-read whole genome sequencing (WGS) using Nextera XT MiSeq failed to identify some virulence genes in plasmids and on the chromosome, both of which were detected using the long-read platforms. Results from long-read MinION and PacBio allowed us to identify differences in plasmid content: a single 88 kb plasmid in CFSAN027343; a 157kb plasmid in CFSAN027350; and two plasmids in CFSAN027346 (one 95 Kb, one 72 Kb). These data enabled rapid characterization of the virulome, detection of antimicrobial genes, and composition/location of Stx phages. Taken together, positive correlations between the two long-read methods for determining plasmids, virulome, antimicrobial resistance genes, and phage composition support MinION sequencing as one accurate and economical option for closing STEC genomes and identifying specific virulence markers.

Whole genome sequencing to characterize shiga toxin-producing Escherichia coli O26 in a public health setting

BACKGROUND

Shiga-toxin producing Escherichia coli (STEC) O26:H11 is the second most common cause of severe diarrhea and hemolytic uremic syndrome worldwide. The implementation of whole genome sequencing (WGS) enhances the detection and in-depth characterization of these non-O157 STEC strains. The aim of this study was to compare WGS to phenotypic serotyping and pulse field gel electrophoresis (PFGE) for characterization of STECO26 strains following a zoonotic outbreak from cattle to humans.

METHODS AND RESULTS

This study evaluated seven E. coli strains; two strains isolated from two children with gastrointestinal symptoms and five strains from five calves suspected as the source of infection. Six of these isolates were serotyped phenotypically and by WGS as E. coli O26:H11 while one bovine isolate could be serotyped only by WGS as E. coli O182:H25. Stx1 was detected in two human- and two bovine-isolates using PCR and WGS. Using WGS, all four STECO26 isolates belong to sequence type (ST) 21 while the two stx1 negative E. coli O26 were ST29. All four STECO26 isolates were indistinguishable by PFGE. However, the data generated by WGS linked the two human STECO26 isolates to only one bovine STECO26 strain by having identical high-quality single nucleotide polymorphisms (hqSNPs) and identical virulence factor profiles while the remaining bovine STECO26 isolate differed by 7 hqSNPs and lacked virulence factor toxB.

CONCLUSIONS

These data demonstrated that WGS provided significant information beyond traditional epidemiological tools allowing for comprehensive characterization of the STEC. Using this approach, WGS was able to identify the specific source of infection in this study.

Virulence gene profiles and phylogeny of Shiga toxin-positive Escherichia coli strains isolated from FDA regulated foods during 2010-2017

Illnesses caused by Shiga toxin-producing Escherichia coli (STECs) can be life threatening, such as hemolytic uremic syndrome (HUS). The STECs most frequently identified by USDA's Microbiological Data Program (MDP) carried toxin gene subtypes stx1a and/or stx2a. Here we described the genome sequences of 331 STECs isolated from foods regulated by the FDA 2010-2017, and determined their genomic identity, serotype, sequence type, virulence potential, and prevalence of antimicrobial resistance. Isolates were selected from the MDP archive, routine food testing by FDA field labs (ORA), and food testing by a contract company. Only 276 (83%) strains were confirmed as STECs by in silico analysis. Foods from which STECs were recovered included cilantro (6%), spinach (25%), lettuce (11%), and flour (9%). Phylogenetic analysis using core genome MLST revealed these STEC genomes were highly variable, with some clustering associated with ST types and serotypes. We detected 95 different sequence types (ST); several ST were previously associated with HUS: ST21 and ST29 (O26:H11), ST11 (O157:H7), ST33 (O91:H14), ST17 (O103:H2), and ST16 (O111:H-). in silico virulome analyses showed ~ 51% of these strains were potentially pathogenic [besides stx gene they also carried eae (25%) or 26% saa (26%)]. Virulence gene prevalence was also determined: stx1 only (19%); stx2 only (66%); and stx1/sxt2 (15%). Our data form a new WGS dataset that can be used to support food safety investigations and monitor the recurrence/emergence of E. coli in foods.

DNA Microarray-Based Genomic Characterization of the Pathotypes of Escherichia coli O26, O45, O103, O111, and O145 Isolated from Feces of Feedlot Cattle †

Shiga toxin-producing Escherichia coli (STEC) serogroups O26, O45, O103, O111, O121, and O145, referred to as the top six non-O157 serogroups, are responsible for more than 70% of human non-O157 STEC infections in North America. Cattle harbor non-O157 strains in the hindgut and shed them in the feces. The objective of this study was to use the U.S. Food and Drug Administration (FDA) E. coli identification (ECID) DNA microarray to identify the serotype, assess the virulence potential of each, and determine the phylogenetic relationships among five of the six non-O157 E. coli serogroups isolated from feedlot cattle feces. Forty-four strains of STEC, enterohemorrhagic E. coli (EHEC), enteropathogenic E. coli (EPEC), or putative nonpathotype E. coli (NPEC) of cattle origin and five human clinical strains of EHEC were assayed with the FDA-ECID DNA microarray. The cattle strains harbored diverse flagellar genes. The bovine and human strains belonging to serogroups O26, O45, and O103 carried stx₁ only, O111 carried both stx₁ and stx₂, and O145 carried either stx₁ or stx₂. The strains were also positive for various subtypes of intimin and other adhesins (IrgA homologue adhesin, long polar fimbriae, mannose-specific adhesin, and curli). Both human and cattle strains were positive for LEE-encoded type III secretory system genes and non-LEE-encoded effector genes. SplitsTree4, a program used to determine the phylogenetic relationship among the strains, revealed that the strains within each serogroup clustered according to their pathotype. In addition to genes encoding Shiga toxins, bovine non-O157 E. coli strains possessed other major virulence genes, including those for adhesins, type III secretory system proteins, and plasmid-borne virulence genes, similar to human clinical strains. Because virulence factors encoded by these genes are involved in the pathogenesis of various pathotypes of E. coli, the bovine non-O157 strains could cause human illness. The FDA-ECID DNA microarray assay rapidly provided a profile of the virulence genes for assessment of the virulence potential of each strain.

Whole genome sequencing based typing and characterisation of Shiga-toxin producing Escherichia coli strains belonging to O157 and O26 serotypes and isolated in dairy farms

In the present study, the genetic relationships as well as the virulome and resistome of newly sequenced O26 and O157 Shiga-toxin producing E. coli (STEC) isolates, collected from dairy farms in Italy, were investigated in comparison to publicly available genomes collected worldwide. The whole genome of Italian isolates was sequenced on Illumina MiSeq Platform. Reads quality control, de novo draft genome assembly, species confirmation and the 7- loci Multi-Locus Sequence Type assignment were performed using INNUca pipeline. Reference-based SNPs calling was performed on O157 and O26 genomes, separately, mapping contigs to high-quality finished genomes. Virulence and antimicrobial resistance determinants were detected in silico using the tool ABRicate. Phylogenetic reconstructions revealed that genomes clustered mainly based on their 7-loci MLST type. The virulome of tested genomes included 190 determinants. O157 genomes carried chu genes associated to heme mediated iron uptake, whereas O26 genomes harboured genes ybt associated to siderophore mediated iron uptake. Resistome analysis showed the presence of tet(34) on all but one O157 genomes and on only one O26 genomes. Only 4 genomes carried genes associated to multiresistance. In the present study, the genes chu and ybt were identified as potential biomarker for the differentiation of O157 and O26 serotypes.

The Present and Future of Whole Genome Sequencing (WGS) and Whole Metagenome Sequencing (WMS) for Surveillance of Antimicrobial Resistant Microorganisms and Antimicrobial Resistance Genes across the Food Chain

Antimicrobial resistance (AMR) surveillance is a critical step within risk assessment schemes, as it is the basis for informing global strategies, monitoring the effectiveness of public health interventions, and detecting new trends and emerging threats linked to food. Surveillance of AMR is currently based on the isolation of indicator microorganisms and the phenotypic characterization of clinical, environmental and food strains isolated. However, this approach provides very limited information on the mechanisms driving AMR or on the presence or spread of AMR genes throughout the food chain. Whole-genome sequencing (WGS) of bacterial pathogens has shown potential for epidemiological surveillance, outbreak detection, and infection control. In addition, whole metagenome sequencing (WMS) allows for the culture-independent analysis of complex microbial communities, providing useful information on AMR genes occurrence. Both technologies can assist the tracking of AMR genes and mobile genetic elements, providing the necessary information for the implementation of quantitative risk assessments and allowing for the identification of hotspots and routes of transmission of AMR across the food chain. This review article summarizes the information currently available on the use of WGS and WMS for surveillance of AMR in foodborne pathogenic bacteria and food-related samples and discusses future needs that will have to be considered for the routine implementation of these next-generation sequencing methodologies with this aim. In particular, methodological constraints that impede the use at a global scale of these high-throughput sequencing (HTS) technologies are identified, and the standardization of methods and protocols is suggested as a measure to upgrade HTS-based AMR surveillance schemes.

Comparative genomics reveals differences in mobile virulence genes of Escherichia coli O103 pathotypes of bovine fecal origin

Escherichia coli O103, harbored in the hindgut and shed in the feces of cattle, can be enterohemorrhagic (EHEC), enteropathogenic (EPEC), or putative non-pathotype. The genetic diversity particularly that of virulence gene profiles within O103 serogroup is likely to be broad, considering the wide range in severity of illness. However, virulence descriptions of the E. coli O103 strains isolated from cattle feces have been primarily limited to major genes, such as Shiga toxin and intimin genes. Less is known about the frequency at which other virulence genes exist or about genes associated with the mobile genetic elements of E. coli O103 pathotypes. Our objective was to utilize whole genome sequencing (WGS) to identify and compare major and putative virulence genes of EHEC O103 (positive for Shiga toxin gene, stx1, and intimin gene, eae; n = 43), EPEC O103 (negative for stx1 and positive for eae; n = 13) and putative non-pathotype O103 strains (negative for stx and eae; n = 13) isolated from cattle feces. Six strains of EHEC O103 from human clinical cases were also included. All bovine EHEC strains (43/43) and a majority of EPEC (12/13) and putative non-pathotype strains (12/13) were O103:H2 serotype. Both bovine and human EHEC strains had significantly larger average genome sizes (P < 0.0001) and were positive for a higher number of adherence and toxin-based virulence genes and genes on mobile elements (prophages, transposable elements, and plasmids) than EPEC or putative non-pathotype strains. The genome size of the three pathotypes positively correlated (R² = 0.7) with the number of genes carried on mobile genetic elements. Bovine strains clustered phylogenetically by pathotypes, which differed in several key virulence genes. The diversity of E. coli O103 pathotypes shed in cattle feces is likely reflective of the acquisition or loss of virulence genes carried on mobile genetic elements.

Population structure of Escherichia coli O26 : H11 with recent and repeated stx2 acquisition in multiple lineages

A key virulence factor of enterohaemorrhagic Escherichia coli (EHEC) is the bacteriophage-encoded Shiga toxin (Stx). Stxs are classified into two types, Stx1 and Stx2, and Stx2-producing strains are thought to cause more severe infections than strains producing only Stx1. Although O26 : H11 is the second most prevalent EHEC following O157 : H7, the majority of O26 : H11 strains produce Stx1 alone. However, Stx2-producing O26 strains have increasingly been detected worldwide. Through a large-scale genome analysis, we present a global phylogenetic overview and evolutionary timescale for E. coli O26 : H11. The origin of O26 has been estimated to be 415 years ago. Sequence type 21C1 (ST21C1), one of the two sublineages of ST21, the most predominant O26 : H11 lineage worldwide, emerged 213 years ago from one of the three ST29 sublineages (ST29C2). The other ST21 lineage (ST21C2) emerged 95 years ago from ST21C1. Increases in population size occurred in the late 20th century for all of the O26 lineages, but most remarkably for ST21C2. Analysis of the distribution of stx2-positive strains revealed the recent and repeated acquisition of the stx2 gene in multiple lineages of O26, both in ST21 and ST29. Other major EHEC virulence genes, such as type III secretion system effector genes and plasmid-encoded virulence genes, were well conserved in ST21 compared to ST29. In addition, more antimicrobial-resistance genes have accumulated in the ST21C1 lineage. Although current attention is focused on several highly virulent ST29 clones that have acquired the stx2 gene, there is also a considerable risk that the ST21 lineage could yield highly virulent clones.

The Mobilome; A Major Contributor to Escherichia coli stx2-Positive O26:H11 Strains Intra-Serotype Diversity

Shiga toxin-producing Escherichia coli of serotype O26:H11/H- constitute a diverse group of strains and several clones with distinct genetic characteristics have been identified and characterized. Whole genome sequencing was performed using Illumina and PacBio technologies on eight stx2-positive O26:H11 strains circulating in France. Comparative analyses of the whole genome of the stx2-positive O26:H11 strains indicate that several clones of EHEC O26:H11 are co-circulating in France. Phylogenetic analysis of the French strains together with stx2-positive and stx-negative E. coli O26:H11 genomes obtained from Genbank indicates the existence of four clonal complexes (SNP-CCs) separated in two distinct lineages, one of which comprises the "new French clone" (SNP-CC1) that appears genetically closely related to stx-negative attaching and effacing E. coli (AEEC) strains. Interestingly, the whole genome SNP (wgSNP) phylogeny is summarized in the cas gene phylogeny, and a simple qPCR assay targeting the CRISPR array specific to SNP-CC1 (SP_O26-E) can distinguish between the two main lineages. The PacBio sequencing allowed a detailed analysis of the mobile genetic elements (MGEs) of the strains. Numerous MGEs were identified in each strain, including a large number of prophages and up to four large plasmids, representing overall 8.7-19.8% of the total genome size. Analysis of the prophage pool of the strains shows a considerable diversity with a complex history of recombination. Each clonal complex (SNP-CC) is characterized by a unique set of plasmids and phages, including stx-prophages, suggesting evolution through separate acquisition events. Overall, the MGEs appear to play a major role in O26:H11 intra-serotype clonal diversification.

Genome sequencing and comparative genomics of enterohemorrhagic Escherichia coli O145:H25 and O145:H28 reveal distinct evolutionary paths and marked variations in traits associated with virulence & colonization

Background

Enterohemorrhagic Escherichia coli (EHEC) O145 are among the top non-O157 serogroups associated with severe human disease worldwide. Two serotypes, O145:H25 and O145:H28 have been isolated from human patients but little information is available regarding the virulence repertoire, origin and evolutionary relatedness of O145:H25. Hence, we sequenced the complete genome of two O145:H25 strains associated with hemolytic uremic syndrome (HUS) and compared the genomes with those of previously sequenced O145:H28 and other EHEC strains.

Results

The genomes of the two O145:H25 strains were 5.3 Mbp in size; slightly smaller than those of O145:H28 and other EHEC strains. Both strains contained three nearly identical plasmids and several prophages and integrative elements, many of which differed significantly in size, gene content and organization as compared to those present in O145:H28 and other EHECs. Furthermore, notable variations were observed in several fimbrial gene cluster and intimin types possessed by O145:H25 and O145:H28 indicating potential adaptation to distinct areas of host colonization. Comparative genomics further revealed that O145:H25 are genetically more similar to other non-O157 EHEC strains than to O145:H28.

Conclusion

Phylogenetic analysis accompanied by comparative genomics revealed that O145:H25 and O145:H28 evolved from two separate clonal lineages and that horizontal gene transfer and gene loss played a major role in the divergence of these EHEC serotypes. The data provide further evidence that ruminants might be a possible reservoir for O145:H25 but that they might be impaired in their ability to establish a persistent colonization as compared to other EHEC strains.

Electronic supplementary material

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

A survey of methods used for antimicrobial susceptibility testing in veterinary diagnostic laboratories in the United States

Antimicrobial resistance is a serious threat to animal and human health worldwide, requiring a collaborative, holistic approach. The U.S. Government has developed a national strategy to address antimicrobial resistance, with one component being to monitor antimicrobial resistance in agricultural settings. We developed a survey to collect information about antimicrobial susceptibility testing (AST) from the veterinary diagnostic laboratory community in the United States, assessing current practices and technologies and determining how AST information is shared. Of the 132 surveys administered, 52 (39%) were returned. Overall, responding laboratories conducted susceptibility tests on 98,788 bacterial isolates in 2014, with Escherichia coli being the most common pathogen tested across all animal species. The 2 most common AST methods employed were the disk diffusion method (71%) and the Sensititre platform broth microdilution system (59%). Laboratories primarily used the Clinical Laboratory Standards Institute (CLSI) VET-01 standard (69%) and the automatically calculated interpretations provided by the commercial AST systems (61%) for interpreting their AST data. Only 22% of laboratories published AST data on a periodic basis, usually via annual reports published on the laboratory's website or through peer-reviewed journals for specific pathogens. Our results confirm that disk diffusion and broth microdilution remain the standard AST methods employed by U.S. veterinary diagnostic laboratories, and that CLSI standards are commonly used for interpreting AST results. This information will help determine the most efficient standardized methodology for future surveillance. Furthermore, the current infrastructure within laboratories, once harmonized, will help provide a mechanism for conducting national surveillance programs.

Defining a Core Genome Multilocus Sequence Typing Scheme for the Global Epidemiology of Vibrio parahaemolyticus

Vibrio parahaemolyticus is an important human foodborne pathogen whose transmission is associated with the consumption of contaminated seafood, with a growing number of infections reported over recent years worldwide. A multilocus sequence typing (MLST) database for V. parahaemolyticus was created in 2008, and a large number of clones have been identified, causing severe outbreaks worldwide (sequence type 3 [ST3]), recurrent outbreaks in certain regions (e.g., ST36), or spreading to other regions where they are nonendemic (e.g., ST88 or ST189). The current MLST scheme uses sequences of 7 genes to generate an ST, which results in a powerful tool for inferring the population structure of this pathogen, although with limited resolution, especially compared to pulsed-field gel electrophoresis (PFGE). The application of whole-genome sequencing (WGS) has become routine for trace back investigations, with core genome MLST (cgMLST) analysis as one of the most straightforward ways to explore complex genomic data in an epidemiological context. Therefore, there is a need to generate a new, portable, standardized, and more advanced system that provides higher resolution and discriminatory power among V. parahaemolyticus strains using WGS data. We sequenced 92 V. parahaemolyticus genomes and used the genome of strain RIMD 2210633 as a reference (with a total of 4,832 genes) to determine which genes were suitable for establishing a V. parahaemolyticus cgMLST scheme. This analysis resulted in the identification of 2,254 suitable core genes for use in the cgMLST scheme. To evaluate the performance of this scheme, we performed a cgMLST analysis of 92 newly sequenced genomes, plus an additional 142 strains with genomes available at NCBI. cgMLST analysis was able to distinguish related and unrelated strains, including those with the same ST, clearly showing its enhanced resolution over conventional MLST analysis. It also distinguished outbreak-related from non-outbreak-related strains within the same ST. The sequences obtained from this work were deposited and are available in the public database (http://pubmlst.org/vparahaemolyticus). The application of this cgMLST scheme to the characterization of V. parahaemolyticus strains provided by different laboratories from around the world will reveal the global picture of the epidemiology, spread, and evolution of this pathogen and will become a powerful tool for outbreak investigations, allowing for the unambiguous comparison of strains with global coverage.

Identification of a New Virulent Clade in Enterohemorrhagic Escherichia coli O26:H11/H- Sequence Type 29

Enterohemorrhagic Escherichia coli (EHEC) O26 infections cause severe human diseases such as hemolytic uremic syndrome and encephalopathy, and is the predominant serogroup among non-O157 EHEC in many countries. Shiga toxin (Stx), which consists of two distinct types (Stx1 and Stx2), plays a central role in EHEC pathogenesis. The major stx gene type in EHEC O26 strains is stx1, although isolates with only stx2 have emerged in Japan since 2012 and have been reported in Europe. In this study, we selected 27 EHEC O26 strains isolated in Japan and identified a distinct genetic clade within sequence type (ST) 29, designated ST29C1, that carried only stx2 and had the plasmid gene profile ehxA+/katP−/espP+/etpD−. We showed that ST29C1 strains produced higher Stx2a levels, and greater virulence in Vero cells and in germ-free mice than other lineages. We also showed that ST29C1 was a distinct phylogenetic clade by SNP analysis using whole genome sequences and clearly differed from the major European EHEC O26 virulent clone, which was designated ST29C2 in this study. The combination of toxin production analysis, virulence analysis in Vero cells and germ-free mice, and phylogenetic analysis identified a newly emerging virulent EHEC clade.