Genomic comparison of Escherichia coli serotype O103:H2 isolates with and without verotoxin genes: implications for risk assessment of strains commonly found in ruminant reservoirs

Abstract

Introduction

Escherichia coli O103:H2 occurs as verotoxigenic E. coli (VTEC) carrying only vtx ₁ or vtx ₂ or both variants, but also as vtx-negative atypical enteropathogenic E. coli (aEPEC). The majority of E. coli O103:H2 identified from cases of human disease are caused by the VTEC form. If aEPEC strains frequently acquire verotoxin genes and become VTEC, they must be considered a significant public health concern. In this study, we have characterized and compared aEPEC and VTEC isolates of E. coli O103:H2 from Swedish cattle.

Methods

Fourteen isolates of E. coli O103:H2 with and without verotoxin genes were collected from samples of cattle feces taken during a nationwide cattle prevalence study 2011–2012. Isolates were sequenced with a 2×100 bp setup on a HiSeq2500 instrument producing >100× coverage per isolate. Single-nucleotide polymorphism (SNP) typing was performed using the genome analysis tool kit (GATK). Virulence genes and other regions of interest were detected. Susceptibility to transduction by two verotoxin-encoding phages was investigated for one representative aEPEC O103:H2 isolate.

Results and Discussion

This study shows that aEPEC O103:H2 is more commonly found (64%) than VTEC O103:H2 (36%) in the Swedish cattle reservoir. The only verotoxin gene variant identified was vtx _1a . Phylogenetic comparison by SNP analysis indicates that while certain subgroups of aEPEC and VTEC are closely related and have otherwise near identical virulence gene repertoires, they belong to separate lineages. This indicates that the uptake or loss of verotoxin genes is a rare event in the natural cattle environment of these bacteria. However, a representative of a VTEC-like aEPEC O103:H2 subgroup could be stably lysogenized by a vtx-encoding phage in vitro.

• VTEC
• STEC
• EHEC
• aEPEC
• cattle
• zoonosis
• next-generation sequencing
• high-throughput sequencing

Pathogenic Escherichia coli have evolved by the acquisition of mobile genetic elements containing virulence and fitness factors through horizontal gene transfer (1). Enteropathogenic E. coli (EPEC) are characterized by variants of the locus of enterocyte effacement (LEE), a genetic element that encodes intimin and components for a type III secretion system as well as several other virulence factors. These factors allow EPEC to cause attaching and effacing (A/E) lesions in the human intestinal epithelium (2). Most human-pathogenic verotoxigenic E. coli (VTEC, also known as STEC) originate from populations of EPEC lysogenized by one or more phages carrying genes encoding verotoxin 1 and/or verotoxin 2 (vtx ₁ , vtx ₂ ) (1). Although typical human-pathogenic EPEC carry the EPEC adherence factor plasmid, this plasmid is frequently absent in VTEC-like EPEC. Therefore, these strains are considered ‘atypical’ EPEC (aEPEC) (2). Ruminants are the natural reservoir of many forms of aEPEC and VTEC (2) and are colonized with a high prevalence and great diversity of these bacteria, but rarely develop symptoms as a result. Both VTEC and the different types of EPEC can cause diarrhea in infected humans, but VTEC can also cause more severe forms of bloody diarrhea and occasionally fatal complications like hemolytic uremic syndrome (HUS) (3). The E. coli O157:H7 serotype is the most common cause of human VTEC infection in Sweden as well as in many countries worldwide, but non-O157 serogroups including O26, O103, and O121 also cause a significant number of cases every year (4). A major outbreak of enteroaggregative VTEC O104:H4 in 2011 caused more than 4,000 cases in several European countries, including Sweden (5),and led to an increased interest in the role of verotoxin-encoding phages in the emergence of new VTEC types.

Verotoxin-encoding prophages can be stably maintained in the bacterial genome, but at least some strains of VTEC can lose the prophage and thereby the ability to produce the verotoxin, reverting to the aEPEC state (6, 7). During VTEC outbreaks, strains that are highly similar to the outbreak strain but lack verotoxin genes can sometimes be isolated from patients or suspected sources (6–8). Conversely, several studies have described successful experimental infection by verotoxin-encoding phages of a broad range of E. coli strains as well as of Shigella sonnei, while in vitro production of stable lysogens appears to be more challenging (9–12). The E. coli serotype O103:H2 occurs as verotoxin-negative aEPEC as well as VTEC carrying only vtx ₁ or vtx ₂ or both variants (13, 14). The verotoxin-negative variants have been shown to be the dominant O103:H2 population in ruminant reservoirs (15, 16), but E. coli O103:H2 with vtx ₁ is the most common cause of human VTEC illness (17–19).

The high prevalence of aEPEC O103:H2 in animals poses a potential threat for the human population. Assuming that frequent conversions by lysogeny of prophages occur in the animal host, in contaminated foodstuffs or even in the gastrointestinal tract of an infected human, the aEPEC strains must be considered a significant risk. In contrast, if the aEPEC and VTEC strains we find in ruminants represent distinct lineages, the conversions are likely to be rare events and the aEPEC strains less of a public health concern. In this study, we have compared aEPEC and VTEC isolates of E. coli O103:H2 from the Swedish cattle reservoir using high-throughput sequencing. We have also investigated the in vitro ability of strains of VTEC O103:H2 to revert to the aEPEC state and the ability of an aEPEC O103:H2 to acquire verotoxingenes by phage transduction.

Materials and methods

A nationwide VTEC O157:H7 prevalence study was performed on cattle at slaughter in Sweden during a one-year period, 2011–2012. One thousand fecal samples and 500 skin (ear) samples were collected from major slaughterhouses, with the number of samples proportional to the number of animals slaughtered at each slaughterhouse. From each collected sample, 10 g of feces or approximately 25 g of skin was pre-enriched in 90 ml of buffered peptone water (BPW) for 6–8 h at 37°C. Aliquots of 1.5 ml were stored at −70°C with 10–15% glycerol added. All samples were screened prior to freezing by real-time polymerase chain reaction (PCR) for the O103 wzx gene as follows. After enrichment, 10 µl of BPW was streaked on MacConkey agar and incubated overnight at 37°C. From the MacConkey plates, DNA extractions were prepared by the suspension of colony material in 500 µl nuclease-free water and heat lysis at 100°C for 10 min. The suspensions were then centrifuged at 13,000×g for 10 min and the supernatant was used as PCR template. PCR was performed according to Perelle et al. (20), with the following modifications: PCR was performed in 15-µl reaction volumes containing 2X PerfeCTa qPCR Toughmix with Low ROX (Quanta Biosciences, Gaithersburg, MD), 500 nM of each primer, 100 nM of probe, and 2 µl template. The PCR was run in an ABI 7500 Fast thermal cycler (Applied Biosystems, Foster City, CA) with the following temperature profile: 95°C for 3 min followed by 40 cycles of 95°C for 3 s and 60°C for 30 s. For this study, E. coli O103:H2 was isolated from a subset of the stored samples based on low C _t value when screening for the O103 wzx gene. For the isolation, the frozen BPW was thawed and mixed with 9 ml fresh BPW and incubated for 4 h at 37°C. One milliliter of each sample was then transferred to 90 ml modified tryptic soy broth with 8 mg novobiocin added and incubated at 41.5°C for 16–20 h. Immunomagnetic separation of these enriched samples was performed using the Dynabeads EPEC/VTEC O103 kit (Life Technologies, Carlsbad, USA), according to the manufacturer's instructions. The beads were distributed on three washed sheep blood agar plates. The plates were supplemented with 8 mg/l vancomycin and 0.05 mg/l cefixime to inhibit unspecific growth, in particular, gram-positive bacteria and swarming Proteus spp., respectively. The plates were incubated at 37°C for 18–24 h. Hemolytic colonies were tested by slide agglutination using sera for the O103 antigen (SSI, Copenhagen, Denmark). Agglutination-positive colonies were confirmed by PCR for the presence of the intimin eae gene (21), the O103 wzx gene (20), and the H2 fliC gene (22) as well as by E. coli O103:H2-specific CRISPR PCR (23). A single isolate from each positive sample was included in the study. A total of 14 isolates were collected and verified in this way and available for sequencing. Isolates that did not fulfill the criteria for aEPEC (intimin) or VTEC (any verotoxin) of the serotype O103:H2 were excluded from further analysis.

For next-generation sequencing, DNA was extracted from colonies using the MasterPure DNA Purification kit (Epicentre, Madison, USA). Sequencing was performed using TruSeq libraries analyzed on a HiSeq2500 instrument (Illumina, San Diego, USA) with a 2×100 bp paired-end setup producing >100× average coverage per isolate. Read data sets were uploaded to the European Nucleotide Archive under project accession number PRJEB11451. For single-nucleotide polymorphism (SNP) typing, raw read data were trimmed and cleaned from duplicates using PrinSeq (24). The trimmed reads were mapped against the complete genome of the E. coli O103:H2 12009 strain (GenBank Acc. AP010958) (25) using Bowtie 2 (26). Local realignment and SNP discovery were performed using the Genome analysis toolkit (GATK) (27), with variants called using the Unified Genotyper algorithm. SNPs were filtered using the following parameters: QUAL (overall probability of polymorphism at position) > 1,000, QD (quality normalized by depth) > 2, FS (strand bias) < 60, ReadPosRankSumTest (alternative/reference allele positional bias) > − 8, MQ (overall mapping quality) > 40, and MQRankSum (alt/ref allele mapping quality bias) > − 12.5. Only SNPs in the chromosomal backbone of the reference were used in the SNP analysis, that is, prophages and integrative elements were excluded. The aligned SNP states in the different isolates were compared using the NeighborNet algorithm in SplitsTree4 (28). Virulence genes were detected in the read data using the CGE VirulenceFinder web tool (29), and the results were verified and manually corrected when necessary, using Bowtie 2 mapped reads against the virulence gene reference sequences.

To test the stability of the verotoxin-encoding prophages in vitro, aliquots of a fecal sample that had been tested and found to be free of any E. coli O103 were spiked with eight vtx ₁ -positive E. coli O103:H2 strains. These fecal samples were then treated like the other project samples as previously described, that is, enriched, frozen, thawed, and enriched again for the re-isolation of the spiked-in isolates. Experimental transduction by verotoxin-encoding phages φ731 (Δvtx2::cat), originally derived from VTEC O103:H25 (30, 31), and 933 W (Δvtx2::cat) (32), originally derived from VTEC O157:H7, was attempted for a single isolate of verotoxin-negative E. coli O103:H2 as previously described (11) with some modifications. Phage lysates were prepared from donor strains (E. coli C600 with phages 933 W, vtx2::cat or φ731, vtx2::cat). Donor strains were grown in LB-broth until the OD measured 0.3–0.5. Mitomycin C (Sigma-Aldrich, Saint Louis, USA) was added to a final concentration of 0.5 µg/ml and incubation was continued at 37°C overnight. Debris was removed by centrifugation at 3,000×g for 15 min and the supernatant was filtered (pore size 0.22 µM) to make a sterile phage lysate. The recipient strain was grown to OD 0.3–0.5 in LB-broth to attain a log-phase growing culture. To 4 ml of LB-broth, 500 µl phage lysate and 500 µl culture of recipient were added. After incubation for 4 h at 37°C, 100 µl was plated on LB-agar containing 25 µg/ml chloramphenicol. The plates were incubated overnight at 37°C before they were checked for resistant strains. Any growth was subsequently re-streaked and incubated on new chloramphenicol plates twice to isolate stably infected strains and remove any free phages that might interfere with PCR confirmation. To confirm that the resistant strains had been infected with the cat-carrying phage as opposed to becoming resistant to chloramphenicol by other means, the cat-gene was amplified using primers Cm-3 and Rho for φ731, and primers VCJ-F2 and VCJ-R2 for 933 W (Table 1). PCR was performed with a T100 thermal cycler (Bio-Rad, Hercules, USA) in 20 µl reactions containing 1 µl template DNA, 400 nM of each primer and 2X Multiplex PCR mastermix (Qiagen, Hilden, Germany). As expected, donor strains were positive and the recipient strain was negative for the presence of vtx2::cat. Four presumptive transductants were tested with PCR. The two isolates infected with φ731 (450 and 451) were PCR positive, confirming successful phage transduction. In contrast, the two isolates apparently infected with φ933 W (444 and 445) were PCR negative (Table 2), indicating that these isolates might have decreased sensitivity towards chloramphenicol resulting from mutations rather than from transduction with the 933 W phage.

Table 1 Primers used for confirmation of phage-infected isolates

Name	Sequence 5′–3′	Reference
Cm-3	catatgaatatcctccttag	(32)
Rho	atatctgcgccgggtctg	(32)
VCJ-F2	catcttcgttgtgtaggctg	This study
VCJ-R2	atacttaactcctcatatgaatatcctc	This study

Table 2 Strains and isolates used in this study

		PCR results
ID	Properties	Cm-3 + Rho	VCJ-F2 + VCJ-R2	Reference
TJ10	C600 933 W, stx2::cat	+	nd	(33)
175	C600 Φ731, stx2::cat	nd	+	(32)
11-BKT072753	O103:H2, eaeA+, stx-	–	–
444/445	Infection of 11-BKT072753 with 933 W	nd	–
450/451	infection of 11-BKT072753 with Φ731	+	nd

nd, not done.

Results and discussion

Among cattle fecal samples originally collected for a nationwide O157:H7 prevalence study, we found vtx-negative (64%) and vtx _1a -positive (36%) isolates of intimin-positive E. coli O103:H2 in cattle, but no O103:H2 isolates with vtx ₂ . In agreement, Sekse et al. found that less than 20% of intimin-positive E. coli O103:H2 isolates from Norwegian sheep encoded vtx ₁ genes, while vtx ₂ was completely absent (15). Furthermore, a recent study in Australia identified only intimin-positive E. coli O103:H2 isolates from ruminants and none that were verotoxin positive (16). Thus, the most common VTEC-like O103:H2 population in ruminants appears to be intimin positive, but without verotoxin genes. A possible alternative explanation would be that the verotoxin genes are lost during the isolation procedure. To test this, we added eight vtx ₁ -positive isolates to portions of an O103-negative cattle fecal sample and re-isolated them, with all isolates remaining verotoxin positive.

The primary goal of this study was to characterize and compare E. coli O103:H2 isolates. Therefore, the isolation of E. coli O103:H2 was not attempted from all samples in the prevalence study, and no prevalence of VTEC/aEPEC O103:H2 in Swedish cattle can be reliably calculated from the data. In addition to the VTEC and aEPEC described here, we found an abundance of E. coli O103:H2 with a hemolysis phenotype that lacked both verotoxin and intimin. For many of these avirulent isolates, the O103:H2 CRISPR PCR (23) for specific detection of E. coli O103:H2 was positive, suggesting these isolates are phylogenetically related to aEPEC and VTEC O103:H2. Similar findings were reported from a prevalence study of E. coli O103:H2 in Norwegian sheep (15). These incidental findings can be a problem when trying to find a particular strain of O103:H2 in ruminants, for example, for investigating suspected cases of infection from animals to humans.

SNP typing of the 14 E. coli aEPEC and VTEC O103:H2 isolates separated them into five major groups, two of which consisted only of VTEC and the remaining three only of aEPEC (Fig. 1). The fact that VTEC and aEPEC form separate lineages in this limited set of isolates might indicate that the uptake or loss of verotoxin-encoding phages is a rare event in the natural cattle environment of E. coli O103:H2. aEPEC O103:H2 may consequently pose less of a risk for human health than we could have assumed if any lineage contained both verotoxin states. As more E. coli O103:H2 genome sequences become available in the next few years, it will be evident whether this observation also holds for lineages of E. coli O103:H2, which are absent in the Swedish cattle reservoir, for example, those carrying vtx ₂ .

Fig. 1 Virulence gene repertoire (left) and single-nucleotide polymorphism variation network (right), inferred from next-generation sequencing of Swedish cattle isolates of E. coli O103:H2. Black squares indicate present genes, white absent. For efa1, a ‘T’ indicates the presence of a truncated gene.

The results of the virulence gene typing are presented in Fig. 1. All isolates had the epsilon intimin variant. All isolates lacked the bfpA gene characteristic of typical EPEC, confirming their status as aEPEC or VTEC. Genes encoding verotoxin 1a were found in all isolates in SNP clusters I and II. Virulence genes in the LEE were conserved, apart from a truncated version of efa1 found in aEPEC SNP clusters IV and V. The same two clusters also appear to lack espJ and the prophage-associated O103:H2 PO2 region, which encodes several putative virulence determinants (34). The large virulence plasmid pO103 was present and intact in all isolates except the only representative of cluster IV, which lacked approximately half the E. coli O103:H2 12009 pO103 plasmid sequence, including the etpD gene. Apart from the presence of verotoxin genes, there were no consistent differences in virulence gene repertoires of the VTEC clusters and the most closely related aEPEC cluster III. However, katP was found only in VTEC cluster I and has been shown to be significantly associated with E. coli O103:H2 isolates from human infections compared with isolates from cattle (17). Thus, the possible existence of distinct lineages of VTEC E. coli O103:H2 that differ in their virulence for the human host warrants further investigation.

The actual benefit for EHEC strains carrying verotoxin genes is controversial (35), but the verotoxin-negative lineages of E. coli O103:H2 are obviously at least as successful as their verotoxin-positive counterparts in the ruminant reservoirs that are the natural habitat of these bacteria. As the local ruminant reservoir is probably the source for a significant portion of cases of E. coli O103:H2 infection in humans, the scarcity of E. coli O103:H2 lineages that carry verotoxin 2, in particular, seems fortunate. Most domestic cases of E. coli O103:H2 infection in Sweden are caused by strains carrying only vtx ₁ and produce comparatively mild symptoms (36). In general, the risk of HUS appears to be higher for those infected by VTEC strains carrying vtx ₂ , in particular, the subtype vtx _2a , while vtx ₁ rarely causes severe disease (37, 38). However, E. coli O103:H2 with only vtx ₁ is a known cause of HUS (18, 19).

A representative of a VTEC-like aEPEC O103:H2 subgroup III (11-BKT072753) was successfully lysogenized by φ731 (Δvtx2::cat) phage transduction in vitro, while transduction by the 933W (Δvtx2::cat) phage failed. Further studies will be needed to explain why only one of the phages was able to lysogenize this particular strain. The authors anticipate that the two phages are able to use the same insertion site in the host genome, so this might be due to other unknown host specificity factors or experimental differences such as different concentrations of the induced phages. The presence of the φ731 phage was persistent after the subculturing of the strain, suggesting stable insertion in the host genome. This shows the ability of Swedish aEPEC E. coli O103:H2 to convert to verotoxin 2-positive status, at least under laboratory conditions. As previous studies have found verotoxin-producing E. coli to be abundant and diverse in the gastrointestinal tract of ruminants (reviewed by Mauro et al. (39)), lack of compatible donor strains is unlikely to explain the rarity or absence of verotoxin 2-positive O103:H2 strains in the cattle reservoir. Historically, VTEC research has, for obvious reasons, mostly been focused on strains directly responsible for causing severe illness in humans, while the extensive diversity of VTEC and VTEC-like bacteria occupying niches in animal reservoirs or the environment has been neglected. This study highlights the importance of collecting and characterizing both VTEC and aEPEC in order to increase our understanding of the role of verotoxin-encoding phages in the evolution and ecology of human-pathogenic VTEC.

Conflict of interest and funding

The authors declare no conflict of interest. This study was financed by the Swedish Civil Contingencies Agency. EB-R was partially supported by the EU FP7 project ALLBIO, project reference: 289452, funded under FP7-KBBE.

Acknowledgements

The authors acknowledge support from Science for Life Laboratory, the National Genomics Infrastructure, NGI, and Uppmax/Uppnex for providing assistance with massive parallel sequencing and computational infrastructure. The authors also thank the diagnostics staff at the molecular biology and VTEC labs at the SVA for assistance with bacteriology and PCR assays.