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.

Pathogenicity assessment of Shiga toxin-producing Escherichia coli strains isolated from wild birds in a major agricultural region in California.

Shiga toxin-producing Escherichia coli (STEC) consists of diverse strains differing in genetic make-up and virulence potential. To better understand the pathogenicity potential of STEC carried by the wildlife, three STEC and one E. coli strains isolated from wild birds near a major agricultural region in California were selected for comparative pathogenomic analyses. Three American crow (Corvus brachyrhynchos) strains, RM9088, RM9513, and RM10410, belonging to phylogroup A with serotypes O109:H48, O9:H30, and O113:H4, respectively, and a red-winged blackbird (Agelaius phoeniceus) strain RM14516 in phylogroup D with serotype O17:H18, were examined. Shiga toxin genes were identified in RM9088 (stx1a), RM10410 (stx1a + stx2d), and RM14516 (stx2a). Unlike STEC O157:H7 strain EDL933, none of the avian STEC strains harbored the pathogenicity islands OI-122, OI-57, and the locus of enterocyte effacement, therefore the type III secretion system biogenesis genes and related effector genes were absent in the three avian STEC genomes. Interestingly, all avian STEC strains exhibited greater (RM9088 and RM14516) or comparable (RM10410) cytotoxicity levels compared with EDL933. Comparative pathogenomic analyses revealed that RM9088 harbored numerous genes encoding toxins, toxins delivery systems, and adherence factors, including heat-labile enterotoxin, serine protease autotransporter toxin Pic, type VI secretion systems, protein adhesin Paa, fimbrial adhesin K88, and colonization factor antigen I. RM9088 also harbored a 36-Kb high pathogenicity island, which is related to iron acquisition and pathogenicity in Yersinia spp. Strain RM14516 carried an acid fitness island like the one in EDL933, containing a nine gene cluster involved in iron acquisition. Genes encoding extracellular serine protease EspP, subtilase cytotoxin, F1C fimbriae, and inverse autotransporter adhesin IatC were only detected in RM14516, and genes encoding serine protease autotransporter EspI and P fimbriae were only identified in RM10410. Although all curli genes were present in avian STEC strains, production of curli fimbriae was only detected for RM9088 and RM14516. Consistently, strong, moderate, and little biofilms were observed for RM9088, RM14516, and RM10410, respectively. Our study revealed novel combinations of virulence factors in two avian strains, which exhibited high level of cytotoxicity and strong biofilm formation. Comparative pathogenomics is powerful in assessing pathogenicity and health risk of STEC strains.

Comparative genomic and phenotypic analyses of the virulence potential in Shiga toxin-producing Escherichia coli O121:H7 and O121:H10.

Shiga toxin-producing Escherichia coli (STEC) O121 is among the top six non-O157 serogroups that are most frequently associated with severe disease in humans. While O121:H19 is predominant, other O121 serotypes have been frequently isolated from environmental samples, but their virulence repertoire is poorly characterized. Here, we sequenced the complete genomes of two animal isolates belonging to O121:H7 and O121:H10 and performed comparative genomic analysis with O121:H19 to assess their virulence potential. Both O121:H7 and O121:H10 strains carry a genome comparable in size with the O121:H19 genomes and belong to phylogroup B1. However, both strains appear to have evolved from a different lineage than the O121:H19 strains according to the core genes-based phylogeny and Multi Locus Sequence Typing. A systematic search of over 300 E. coli virulence genes listed in the Virulence Factor DataBase revealed a total of 73 and 71 in O121:H7 and O121:H10 strains, respectively, in comparison with an average of 135 in the O121:H19 strains. This variation in the virulence genes repertoire was mainly attributed to the reduction in the number of genes related to the Type III Secretion System in the O121:H7 and O121:H10 strains. Compared to the O121:H19 strains, the O121:H7 strain carries more adherence and toxin genes while the O121:H10 strain carries more genes related to the Type VI Secretion System. Although both O121:H7 and O121:H10 strains carry the large virulence plasmid pEHEC, they do not harbor all pEHEC virulence genes in O121:H19. Furthermore, unlike the O121:H19 strains, neither the O121:H7 nor O121:H10 strain carried the Locus of Enterocyte Effacement, OI-122, nor the tellurite resistance island. Although an incomplete Locus of Adhesion and Autoaggregation (LAA) was identified in the O121:H7 and O121:H10 strains, a limited number of virulence genes were present. Consistently, both O121:H7 and O121:H10 strains displayed significant reduced cytotoxicity than either the O157:H7 strain EDL933 or the O121:H19 strain RM8352. In fact, the O121:H7 strain RM8082 appeared to cause minimal cytotoxicity to Vero cells. Our study demonstrated distinct evolutionary lineages among the strains of serotypes O121:H19, O121:H10, and O121:H7 and suggested reduced virulence potentials in STEC strains of O121:H10 and O121:H7.

DNA adenine methylase, not the PstI restriction-modification system, regulates virulence gene expression in Shiga toxin-producing Escherichia coli.

We previously reported a distinct methylome between the two Shiga toxin-producing Escherichia coli (STEC) O145:H28 strains linked to the 2010 U.S. lettuce-associated outbreak (RM13514) and the 2007 Belgium ice cream-associated outbreak (RM13516), respectively. This difference was thought to be attributed to a prophage encoded type II restriction-modification system (PstI R-M) in RM13514. Here, we characterized this PstI R-M system in comparison to DNA adenine methylase (Dam), a highly conserved enzyme in γ proteobacteria, by functional genomics. Deficiency in Dam led to a differential expression of over 1000 genes in RM13514, whereas deficiency in PstI R-M only impacted a few genes transcriptionally. Dam regulated genes involved in diverse functions, whereas PstI R-M regulated genes mostly encoding transporters and adhesins. Dam regulated a large number of genes located on prophages, pathogenicity islands, and plasmids, including Shiga toxin genes, type III secretion system (TTSS) genes, and enterohemolysin genes. Production of Stx2 in dam mutant was significantly higher than in RM13514, supporting a role of Dam in maintaining lysogeny of Stx2-prophage. However, following mitomycin C treatment, Stx2 in RM13514 was significantly higher than that of dam or PstI R-M deletion mutant, implying that both Dam and PstI R-M contributed to maximum Stx2 production.

Differential induction of Shiga toxin in environmental Escherichia coli O145:H28 strains carrying the same genotype as the outbreak strains.

Shiga toxin-producing Escherichia coli (STEC) O145 is a major serotype associated with severe human disease. Production of Shiga toxins (Stxs), especially Stx2a, is thought to be correlated with STEC virulence. Since stx genes are located in prophages genomes, induction of prophages is required for effective Stxs production. Here, we investigated the production of Stxs in 12 environmental STEC O145:H28 strains under stresses STEC encounter in natural habitats and performed comparative analysis with two O145:H28 clinical strains, one linked to a 2010 U.S. lettuce-associated outbreak (RM13514) and the other linked to a 2007 Belgium ice cream-associated outbreak (RM13516). Similar to the outbreak strains, all environmental strains belong to Sequence Type (ST)-78 using the EcMLST typing scheme. Although all Stx1a-prophages were grouped together, variations in Stx1a production were observed prior to or following the inductions. Among all stx2a positive environmental strains, only the Stx2a-prophage in cattle isolate RM9154-C1 was clustered with the Stx2a-prophages in RM13514, the Stx2a-phage induced from a STEC O104:H4 strain linked to the 2011 outbreak of enterohemorrhagic infection in Germany, and the Stx2a-prophage in STEC O157:H7 strain EDL933, a prototype of enterohemorrhagic E. coli. Furthermore, the Stx2a-prophage in RM9154-C1 shared the same chromosomal insertion site and carried the same antiterminator Q gene and the late promoter PR' as the Stx2a-prophage in RM13514. Following mitomycin C or enrofloxacin treatment, the production of Stx2a in RM9154-C1 was the highest among all environmental strains tested. In contrast, following acid challenge and recovery, the production of Stx2a in RM9154-C1 was the lowest among all the environmental strains tested, at a level comparable to the clinical strains. A significant increase in Stx2a production was detected in all strains when exposed to H2O2, although the induction fold was much lower than those by other inducers. This low-efficiency induction of Stx-prophages by H2O2, a natural inducer of Stx-prophages, supports the hypothesis of bacterial altruism in controlling Stxs production, a strategy that assures the survival of the STEC population as a whole by sacrificing a small fraction of cells for Stxs production and release. Differential induction of Stxs among strains carrying nearly identical Stx-prophages suggests a role of host bacteria in regulating Stxs production. Our study revealed diverse Stx-prophages in STEC O145:H28 strains that were genotypically indistinguishable. Identification of a cattle isolate harboring a Stx2a-prophage associated with high virulence supports the premise that cattle, a natural reservoir of STEC, serve as a source of hypervirulent STEC strains.

Genomic Insight into Natural Inactivation of Shiga Toxin 2 Production in an Environmental Escherichia coli Strain Producing Shiga Toxin 1.

Shiga toxin-producing Escherichia coli (STEC) consists of a group of diverse strains differing greatly in genetic make-up and pathogenicity potential. Here, we investigated production of Shiga toxins (Stxs) in a bovine isolate carrying multiple Shiga toxin genes (stxs) after exposure to several antibiotics commonly used in food animals. Strain RM10809-C3 was co-isolated with a STEC O145:H28 strain from cattle feces near a leafy greens-growing region in California. The genome of RM10809-C3 is composed of a 5,128,479-bp chromosome and a 122,641-bp plasmid, encoding 5108 coding sequences. Strain RM10809-C3 belongs to serotype O22:H8 and is clustered together with two STEC O168:H8 food isolates using either multilocus sequence type or core genome-based phylogenetic analysis. Six intact prophages were identified in the genome of RM10809-C3, among which prophage 4 contained two sets of stx2d; whereas prophage 9 carried one set of stx1a. Increased production of Stx1 was detected in RM10809-C3 after exposure to mitomycin C and enrofloxacin, but not in cells exposed to tetracycline. In contrast, Stx2 remained undetectable in cells treated with any of the antibiotics examined. Comparison of Stx-converting prophages in strain RM10809-C3 with those in strain EDL933 revealed altered stx2 promoters in RM10809-C3, including deletion of the late promoter PR' and the mutations in qut, the binding site of antitermination protein Q. In contrast, both PR' and qut within the promoter of stx1 in RM10809-C3 were identical to the corresponding one in EDL933. Further, the protein Q encoded by Stx1-prophage in RM10809-C3 exhibited >94% identity with either of the two EDL933 protein Q; whereas both protein Q encoded by Stx2-prophage in RM10809-C3 were distantly related to any of the EDL933 protein Q. Natural silence of Stx2 production in strain RM10809-C3 emphasizes that not only the stx coding regions but also their regulatory factors are important in STEC risk assessment.

Complete Genome Sequences of Two Shiga Toxin-Producing Escherichia coli Strains Isolated from Crows.

Escherichia coli strains RM9088 and RM10410 were isolated from crows near a leafy greens-growing region in California in April and July 2009, respectively. Both strains carry genes encoding Shiga toxins and other virulence factors in enteric pathogens. Here, we report the complete genome sequences of RM9088 and RM10410.

A Clonal Shiga Toxin-Producing Escherichia coli O121:H19 Population Exhibits Diverse Carbon Utilization Patterns.

Shiga toxin-producing Escherichia coli (STEC) serotype O121:H19 is one of the major non-O157:H7 serotypes associated with severe human disease. Here we examined population structure, virulence potential, and metabolic profile of environmental STEC O121 strains recovered from a major produce production region in California and performed comparative analyses with STEC O121 clinical isolates. Multilocus sequence typing revealed that sequence type (ST)-655, a common ST in clinical strains, was the predominant genotype among the environmental strains. Phylotyping placed all STEC O121 strains in B1 group, a lineage containing other major non-O157 serogroups of STEC. Genes encoding different subtypes of Shiga toxin 1 and 2 were detected in O121, including stx1a, stx1d, stx2a, and stx2e. Furthermore, genes encoding intimin (eae) and enterohemolysin (ehxA) were detected in a majority of environmental strains (83.3%), suggesting that the majority of environmental STEC O121 strains are enterohemorrhagic E. coli. The STEC O121 strains with the same genotype were clustered together based on the carbon utilization pattern. Among the 122 carbon substrates that supported the growth of STEC O121 strains, 44 and 35 exhibited lineage (ST) and strain-specific metabolic profiles, respectively. Although clinical ST-655 strains displayed higher metabolic activity than environmental ST-655 strains for several carbon substrates, including l-alaninamide, 5-keto-d-gluconic acid, 3-O-ß-d-galactopyranosyl-d-arabinose, α-ketoglutaric acid, and lactulose, a few environmental strains with the enhanced metabolic potential for the above substrates were detected. Variations in curli biogenesis and swimming motility were also observed in ST-655 strains, suggesting that phenotypic variants are widespread in STEC. Considering the ecological niches that STEC colonizes, increased metabolic potential for plant-derived carbohydrates, mucus-derived substrates, or secondary metabolites produced by the indigenous microorganisms might have been selected. Such traits would confer STEC competitive advantages and facilitate survival and adaptation of STEC population to a given niche, including infected humans.

Complete Genome Sequences of Two Atypical Enteropathogenic Escherichia coli O145 Environmental Strains.

Escherichia coli O145 strains RM14715 and RM14723 were isolated from wildlife feces near a leafy greens-growing region in Yuma, Arizona. Both strains carry a distinct genotype compared with the E. coli O145 strains isolated from Salinas Valley, California. Here we report complete genome sequences and annotations of RM14715 and RM14723.

Complete Genome Sequence of a Natural Escherichia coli O145:H11 Isolate That Belongs to Phylogroup A.

Escherichia coli O145:H11 strain RM14721 was originally isolated from wildlife feces near a leafy greens-growing region in Yuma, AZ. This strain was initially positive for stx1; however, in subsequent cultures, stx1 was not detected by PCR. Here, we report the complete genome sequence and annotation of RM14721.

Complete Genome Sequence of a Shiga Toxin-Producing Enterobacter cloacae Clinical Isolate.

Enterobacter cloacae strain M12X01451 was isolated from a patient with mild diarrhea. This strain produces a novel subtype of Shiga toxin 1, Stx1e. The Stx1e-converting prophage in strain M12X01451 is stable and can infect other bacteria following induction. Here we report the complete genome sequence and annotation of strain M12X01451.