Optical maps distinguish individual strains of Escherichia coli O157 : H7

Rates of evolutionary change of resident Escherichia coli O157:H7 differ within the same ecological niche

Background

Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a pathogen known to reside in cattle feedlots. This retrospective study examined 181 STEC O157:H7 strains collected over 23 years from a closed-system feedlot. All strains were subjected to short-read sequencing, with a subset of 36 also subjected to long-read sequencing.

Results

Over 96% of the strains fell into four phylogenetically distinct clades. Clade membership was associated with multiple factors including stx composition and the alleles of a well-characterized polymorphism (tir 255 T > A). Small plasmids (2.7 to 40 kb) were found to be primarily clade specific. Within each clade, chromosomal rearrangements were observed along with a core phageome and clade specific phages. Across both core and mobile elements of the genome, multiple SNP alleles were in complete linkage disequilibrium across all strains within specific clades. Clade evolutionary rates varied between 0.9 and 2.8 SNP/genome/year with two tir A allele clades having the lowest evolutionary rates. Investigation into possible causes of the differing rates was not conclusive but revealed a synonymous based mutation in the DNA polymerase III of the fastest evolving clade. Phylogenetic trees generated through our bioinformatic pipeline versus the NCBI’s pathogen detection project were similar, with the two tir A allele clades matching individual NCBI SNP clusters, and the two tir T allele clades assigned to multiple closely-related SNP clusters.

Conclusions

In one ecological niche, a diverse STEC O157:H7 population exhibited different rates of evolution that associated with SNP alleles in linkage disequilibrium in the core genome and mobile elements, including tir 255 T > A.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08497-6.

Bistable Expression of a Toxin-Antitoxin System Located in a Cryptic Prophage of Escherichia coli O157:H7

Type II toxin-antitoxin (TA) systems are classically composed of two genes that encode a toxic protein and a cognate antitoxin protein. Both genes are organized in an operon whose expression is autoregulated at the level of transcription by the antitoxin-toxin complex, which binds operator DNA through the antitoxin’s DNA-binding domain. Here, we investigated the transcriptional regulation of a particular TA system located in the immunity region of a cryptic lambdoid prophage in the Escherichia coli O157:H7 EDL933 strain. This noncanonical paaA2-parE2 TA operon contains a third gene, paaR2, that encodes a transcriptional regulator that was previously shown to control expression of the TA. We provide direct evidence that the PaaR2 is a transcriptional regulator which shares functional similarities to the lambda CI repressor. Expression of the paaA2-parE2 TA operon is regulated by two other transcriptional regulators, YdaS and YdaT, encoded within the same region. We argue that YdaS and YdaT are analogous to lambda Cro and CII and that they do not constitute a TA system, as previously debated. We show that PaaR2 primarily represses the expression of YdaS and YdaT, which in turn controls the expression of paaR2-paaA2-parE2 operon. Overall, our results show that the paaA2-parE2 TA is embedded in an intricate lambdoid prophage-like regulation network. Using single-cell analysis, we observed that the entire locus exhibits bistability, which generates diversity of expression in the population. Moreover, we confirmed that paaA2-parE2 is addictive and propose that it could limit genomic rearrangements within the immunity region of the CP-933P cryptic prophage.

Genome structural variation in Escherichia coli O157:H7

The human zoonotic pathogen Escherichia coli O157:H7 is defined by its extensive prophage repertoire including those that encode Shiga toxin, the factor responsible for inducing life-threatening pathology in humans. As well as introducing genes that can contribute to the virulence of a strain, prophage can enable the generation of large-chromosomal rearrangements (LCRs) by homologous recombination. This work examines the types and frequencies of LCRs across the major lineages of the O157:H7 serotype. We demonstrate that LCRs are a major source of genomic variation across all lineages of E. coli O157:H7 and by using both optical mapping and Oxford Nanopore long-read sequencing prove that LCRs are generated in laboratory cultures started from a single colony and that these variants can be recovered from colonized cattle. LCRs are biased towards the terminus region of the genome and are bounded by specific prophages that share large regions of sequence homology associated with the recombinational activity. RNA transcriptional profiling and phenotyping of specific structural variants indicated that important virulence phenotypes such as Shiga-toxin production, type-3 secretion and motility can be affected by LCRs. In summary, E. coli O157:H7 has acquired multiple prophage regions over time that act to continually produce structural variants of the genome. These findings raise important questions about the significance of this prophage-mediated genome contingency to enhance adaptability between environments.

Chronological set of E. coli O157:H7 bovine strains establishes a role for repeat sequences and mobile genetic elements in genome diversification

Background

Enterohemorrhagic Escherichia coli O157:H7 (EHEC) is a significant foodborne pathogen that resides asymptomatically within cattle and other ruminants. The EHEC genome harbors an extensive collection of mobile genetic elements (MGE), including multiple prophage, prophage-like elements, plasmids, and insertion sequence (IS) elements.

Results

A chronological collection of EHEC strains (FRIK804, FRIK1275, and FRIK1625) isolated from a Wisconsin dairy farm (farm X) comprised a closely related clade genetically differentiated by structural alterations to the chromosome. Comparison of the FRIK804 genome with a reference EHEC strain Sakai found a unique prophage like element (PLE, indel 1) and an inversion (1.15 Mb) situated symmetrically with respect to the terminus region. Detailed analysis determined the inversion was due to homologous recombination between repeat sequences in prophage. The three farm X strains were distinguished by the presence or absence of indel 3 (61 kbp) and indel 4 (48 kbp); FRIK804 contained both of these regions, FRIK1275 lacked indel 4, and indels 3 and 4 were both absent in FRIK1625. Indel 3 was the stx2 prophage and indel 4 involved a deletion between two adjacent prophage with shared repeat sequences. Both FRIK804 and FRIK1275 produced functional phage while FRIK1625 did not, which is consistent with indel 3. Due to their involvement in recombination events, direct and inverted repeat sequences were identified, and their locations mapped to the chromosome. FRIK804 had a greater number and overall length of repeat sequences than E. coli K12 strain MG1655. Repeat sequences were most commonly associated with MGE.

Conclusions

This research demonstrated that three EHEC strains from a Wisconsin dairy farm were closely related and distinguished by variability within prophage regions and other MGE. Chromosome alterations were associated with recombination events between repeat sequences. An inventory of direct and inverted repeat sequences found a greater abundance and total length of repeat sequences in the EHEC strains compared to E. coli strain MG1655. The locations of the repeat sequences were biased towards MGE. The findings from this study expand our understanding of the precise molecular events and elements that contributed to genetic diversification of wild-type EHEC in the bovine and farm environments.

Cultivation-Free Typing of Bacteria Using Optical DNA Mapping

A variety of pathogenic bacteria can infect humans, and rapid species identification is crucial for the correct treatment. However, the identification process can often be time-consuming and depend on the cultivation of the bacterial pathogen(s). Here, we present a stand-alone, enzyme-free, optical DNA mapping assay capable of species identification by matching the intensity profiles of large DNA molecules to a database of fully assembled bacterial genomes (>10 000). The assay includes a new data analysis strategy as well as a general DNA extraction protocol for both Gram-negative and Gram-positive bacteria. We demonstrate that the assay is capable of identifying bacteria directly from uncultured clinical urine samples, as well as in mixtures, with the potential to be discriminative even at the subspecies level. We foresee that the assay has applications both within research laboratories and in clinical settings, where the time-consuming step of cultivation can be minimized or even completely avoided.

Analysis of genetic diversity of Xanthomonas oryzae pv. oryzae populations in Taiwan

Rice bacterial blight caused by Xanthomonas oryzae pv. oryzae (Xoo) is a major rice disease. In Taiwan, the tropical indica type of Oryza sativa originally grown in this area is mix-cultivated with the temperate japonica type of O. sativa, and this might have led to adaptive changes of both rice host and Xoo isolates. In order to better understand how Xoo adapts to this unique environment, we collected and analyzed fifty-one Xoo isolates in Taiwan. Three different genetic marker systems consistently identified five groups. Among these groups, two of them had unique sequences in the last acquired ten spacers in the clustered regularly interspaced short palindromic repeats (CRISPR) region, and the other two had sequences that were similar to the Japanese isolate MAFF311018 and the Philippines isolate PXO563, respectively. The genomes of two Taiwanese isolates with unique CRISPR sequence features, XF89b and XM9, were further completely sequenced. Comparison of the genome sequences suggested that XF89b is phylogenetically close to MAFF311018, and XM9 is close to PXO563. Here, documentation of the diversity of groups of Xoo in Taiwan provides evidence of the populations from different sources and hitherto missing information regarding distribution of Xoo populations in East Asia.

Applications of optical DNA mapping in microbiology

Optical mapping (OM) has been used in microbiology for the past 20 years, initially as a technique to facilitate DNA sequence-based studies; however, with decreases in DNA sequencing costs and increases in sequence output from automated sequencing platforms, OM has grown into an important auxiliary tool for genome assembly and comparison. Currently, there are a number of new and exciting applications for OM in the field of microbiology, including investigation of disease outbreaks, identification of specific genes of clinical and/or epidemiological relevance, and the possibility of single-cell analysis when combined with cell-sorting approaches. In addition, designing lab-on-a-chip systems based on OM is now feasible and will allow the integrated and automated microbiological analysis of biological fluids. Here, we review the basic technology of OM, detail the current state of the art of the field, and look ahead to possible future developments in OM technology for microbiological applications.

Advances in Molecular Serotyping and Subtyping of Escherichia coli

Escherichia coli plays an important role as a member of the gut microbiota; however, pathogenic strains also exist, including various diarrheagenic E. coli pathotypes and extraintestinal pathogenic E. coli that cause illness outside of the GI-tract. E. coli have traditionally been serotyped using antisera against the ca. 186 O-antigens and 53 H-flagellar antigens. Phenotypic methods, including bacteriophage typing and O- and H- serotyping for differentiating and characterizing E. coli have been used for many years; however, these methods are generally time consuming and not always accurate. Advances in next generation sequencing technologies have made it possible to develop genetic-based subtyping and molecular serotyping methods for E. coli, which are more discriminatory compared to phenotypic typing methods. Furthermore, whole genome sequencing (WGS) of E. coli is replacing established subtyping methods such as pulsed-field gel electrophoresis, providing a major advancement in the ability to investigate food-borne disease outbreaks and for trace-back to sources. A variety of sequence analysis tools and bioinformatic pipelines are being developed to analyze the vast amount of data generated by WGS and to obtain specific information such as O- and H-group determination and the presence of virulence genes and other genetic markers.

EHEC Genomics: Past, Present, and Future

This article examines the role of genomics in the understanding and identification of O157:H7 enterohemorrhagic Escherichia coli (EHEC). We highlight the development of novel molecular typing systems that are based on the genomic sequence that has been generated for this pathotype. The genomic comparisons of EHEC to other E. coli strains highlight the close relatedness of the O157 and O55 isolates and also identify other non-O157 clades of isolates that appear to have a different genomic history. Analysis within the EHEC isolates must be completed on a fine scale using whole-genome sequence-based approaches to assess both the conserved and lateral acquired gene content. The plethora of genomic data for EHEC isolates has provided the ability to examine this pathotype in detail, which has provided opportunities for novel surveillance, detection, and diagnostics.

Whole-Genome Mapping as a Novel High-Resolution Typing Tool for Legionella pneumophila

Legionella is the causative agent for Legionnaires' disease (LD) and is responsible for several large outbreaks in the world. More than 90% of LD cases are caused by Legionella pneumophila, and studies on the origin and transmission routes of this pathogen rely on adequate molecular characterization of isolates. Current typing of L. pneumophila mainly depends on sequence-based typing (SBT). However, studies have shown that in some outbreak situations, SBT does not have sufficient discriminatory power to distinguish between related and nonrelated L. pneumophila isolates. In this study, we used a novel high-resolution typing technique, called whole-genome mapping (WGM), to differentiate between epidemiologically related and nonrelated L. pneumophila isolates. Assessment of the method by various validation experiments showed highly reproducible results, and WGM was able to confirm two well-documented Dutch L. pneumophila outbreaks. Comparison of whole-genome maps of the two outbreaks together with WGMs of epidemiologically nonrelated L. pneumophila isolates showed major differences between the maps, and WGM yielded a higher discriminatory power than SBT. In conclusion, WGM can be a valuable alternative to perform outbreak investigations of L. pneumophila in real time since the turnaround time from culture to comparison of the L. pneumophila maps is less than 24 h.

A fast and scalable kymograph alignment algorithm for nanochannel-based optical DNA mappings

Optical mapping by direct visualization of individual DNA molecules, stretched in nanochannels with sequence-specific fluorescent labeling, represents a promising tool for disease diagnostics and genomics. An important challenge for this technique is thermal motion of the DNA as it undergoes imaging; this blurs fluorescent patterns along the DNA and results in information loss. Correcting for this effect (a process referred to as kymograph alignment) is a common preprocessing step in nanochannel-based optical mapping workflows, and we present here a highly efficient algorithm to accomplish this via pattern recognition. We compare our method with the one previous approach, and we find that our method is orders of magnitude faster while producing data of similar quality. We demonstrate proof of principle of our approach on experimental data consisting of melt mapped bacteriophage DNA.

Whole-genome optical mapping and finished genome sequence of Sphingobacterium deserti sp. nov., a new species isolated from the Western Desert of China

A novel Gram-negative bacterium, designated ZWT, was isolated from a soil sample of the Western Desert of China, and its phenotypic properties and phylogenetic position were investigated using a polyphasic approach. Growth occurred on TGY medium at 5-42°C with an optimum of 30°C, and at pH 7.0-11.0 with an optimum of pH 9.0. The predominant cellular fatty acids were summed feature 3 (C16:1ω7c/C16:1ω6c or C16:1ω6c/C16:1ω7c) (39.22%), iso-C15:0 (27.91%), iso-C17:0 3OH (15.21%), C16:0 (4.98%), iso-C15:0 3OH (3.03%), C16:0 3OH (5.39%) and C14:0 (1.74%). The major polar lipid of strain ZWT is phosphatidylethanolamine. The only menaquinone observed was MK-7. The GC content of the DNA of strain ZWT is 44.9 mol%. rDNA phylogeny, genome relatedness and chemotaxonomic characteristics all indicate that strain ZWT represents a novel species of the genus Sphingobacterium. We propose the name S. deserti sp. nov., with ZWT (= KCTC 32092T = ACCC 05744T) as the type strain. Whole genome optical mapping and next-generation sequencing was used to derive a finished genome sequence for strain ZWT, consisting of a circular chromosome of 4,615,818 bp in size. The genome of strain ZWT features 3,391 protein-encoding and 48 tRNA-encoding genes. Comparison of the predicted proteome of ZWT with those of other sphingobacteria identified 925 species-unique proteins that may contribute to the adaptation of ZWT to its native, extremely arid and inhospitable environment. As the first finished genome sequence for any Sphingobacterium, our work will serve as a useful reference for subsequent sequencing and mapping efforts for additional strains and species within this genus.

Large chromosomal rearrangements during a long-term evolution experiment with Escherichia coli

Large-scale rearrangements may be important in evolution because they can alter chromosome organization and gene expression in ways not possible through point mutations. In a long-term evolution experiment, twelve Escherichia coli populations have been propagated in a glucose-limited environment for over 25 years. We used whole-genome mapping (optical mapping) combined with genome sequencing and PCR analysis to identify the large-scale chromosomal rearrangements in clones from each population after 40,000 generations. A total of 110 rearrangement events were detected, including 82 deletions, 19 inversions, and 9 duplications, with lineages having between 5 and 20 events. In three populations, successive rearrangements impacted particular regions. In five populations, rearrangements affected over a third of the chromosome. Most rearrangements involved recombination between insertion sequence (IS) elements, illustrating their importance in mediating genome plasticity. Two lines of evidence suggest that at least some of these rearrangements conferred higher fitness. First, parallel changes were observed across the independent populations, with ~65% of the rearrangements affecting the same loci in at least two populations. For example, the ribose-utilization operon and the manB-cpsG region were deleted in 12 and 10 populations, respectively, suggesting positive selection, and this inference was previously confirmed for the former case. Second, optical maps from clones sampled over time from one population showed that most rearrangements occurred early in the experiment, when fitness was increasing most rapidly. However, some rearrangements likely occur at high frequency and may have simply hitchhiked to fixation. In any case, large-scale rearrangements clearly influenced genomic evolution in these populations.

Bacterial chromosomes are dynamic structures shaped by long histories of evolution. Among genomic changes, large-scale DNA rearrangements can have important effects on the presence, order, and expression of genes. Whole-genome sequencing that relies on short DNA reads cannot identify all large-scale rearrangements. Therefore, deciphering changes in the overall organization of genomes requires alternative methods, such as optical mapping. We analyzed the longest-running microbial evolution experiment (more than 25 years of evolution in the laboratory) by optical mapping, genome sequencing, and PCR analyses. We found multiple large genome rearrangements in all 12 independently evolving populations. In most cases, it is unclear whether these changes were beneficial themselves or, alternatively, hitchhiked to fixation with other beneficial mutations. In any case, many genome rearrangements accumulated over decades of evolution, providing these populations with genetic plasticity reminiscent of that observed in some pathogenic bacteria.

Future perspectives, applications and challenges of genomic epidemiology studies for food-borne pathogens: A case study of Enterohemorrhagic Escherichia coli (EHEC) of the O157:H7 serotype

The shiga-toxin (Stx)-producing human pathogen Escherichia coli serotype O157:H7 is a highly pathogenic subgroup of Stx-producing E. coli (STEC) with food-borne etiology and bovine reservoir. Each year in the U. S., approximately 100,000 patients are infected with enterohemorrhagic E. coli (EHEC) of the O157:H7 serotype. This food-borne pathogen is a global public health threat responsible for widespread outbreaks of human disease. Since its initial discovery in 1982, O157:H7 has rapidly become the dominant EHEC serotype in North America. Hospitalization rates among patients as high as 50% have been reported for severe outbreaks of human disease. Symptoms of disease can rapidly deteriorate and progress to life-threatening complications such as Hemolytic Uremic Syndrome (HUS), the leading cause of kidney failure in children, or Hemorrhagic Colitis. In depth understanding of the genomic diversity that exists among currently circulating EHEC populations has broad applications for improved molecular-guided biosurveillance, outbreak preparedness, diagnostic risk assessment, and development of alternative toxin-suppressing therapeutics.

Precise dissection of an Escherichia coli O157:H7 outbreak by single nucleotide polymorphism analysis

The current pathogen-typing methods have suboptimal sensitivities and specificities. DNA sequencing offers an opportunity to type pathogens with greater degrees of discrimination using single nucleotide polymorphisms (SNPs) than with pulsed-field gel electrophoresis (PFGE) and other methodologies. In a recent cluster of Escherichia coli O157:H7 infections attributed to salad bar exposures and romaine lettuce, a subset of cases denied exposure to either source, although PFGE and multiple-locus variable-number tandem-repeat analysis (MLVA) suggested that all isolates had the same recent progenitor. Interrogation of a preselected set of 3,442,673 nucleotides in backbone open reading frames (ORFs) identified only 1 or 2 single nucleotide differences in 3 of 12 isolates from the cases who denied exposure. The backbone DNAs of 9 of 9 and 3 of 3 cases who reported or were unsure about exposure, respectively, were isogenic. Backbone ORF SNP set sequencing offers pathogen differentiation capabilities that exceed those of PFGE and MLVA.

Whole genome mapping and re-organization of the nuclear and mitochondrial genomes of Babesia microti isolates

Babesia microti is the primary causative agent of human babesiosis, an emerging pathogen that causes a malaria-like illness with possible fatal outcome in immunocompromised patients. The genome sequence of the B. microti R1 strain was reported in 2012 and revealed a distinct evolutionary path for this pathogen relative to that of other apicomplexa. Lacking from the first genome assembly and initial molecular analyses was information about the terminal ends of each chromosome, and both the exact number of chromosomes in the nuclear genome and the organization of the mitochondrial genome remained ambiguous. We have now performed various molecular analyses to characterize the nuclear and mitochondrial genomes of the B. microti R1 and Gray strains and generated high-resolution Whole Genome maps. These analyses show that the genome of B. microti consists of four nuclear chromosomes and a linear mitochondrial genome present in four different structural types. Furthermore, Whole Genome mapping allowed resolution of the chromosomal ends, identification of areas of misassembly in the R1 genome, and genomic differences between the R1 and Gray strains, which occur primarily in the telomeric regions. These studies set the stage for a better understanding of the evolution and diversity of this important human pathogen.

High Resolution Typing by Whole Genome Mapping Enables Discrimination of LA-MRSA (CC398) Strains and Identification of Transmission Events

After its emergence in 2003, a livestock-associated (LA-)MRSA clade (CC398) has caused an impressive increase in the number of isolates submitted for the Dutch national MRSA surveillance and now comprises 40% of all isolates. The currently used molecular typing techniques have limited discriminatory power for this MRSA clade, which hampers studies on the origin and transmission routes. Recently, a new molecular analysis technique named whole genome mapping was introduced. This method creates high-resolution, ordered whole genome restriction maps that may have potential for strain typing. In this study, we assessed and validated the capability of whole genome mapping to differentiate LA-MRSA isolates. Multiple validation experiments showed that whole genome mapping produced highly reproducible results. Assessment of the technique on two well-documented MRSA outbreaks showed that whole genome mapping was able to confirm one outbreak, but revealed major differences between the maps of a second, indicating that not all isolates belonged to this outbreak. Whole genome mapping of LA-MRSA isolates that were epidemiologically unlinked provided a much higher discriminatory power than spa-typing or MLVA. In contrast, maps created from LA-MRSA isolates obtained during a proven LA-MRSA outbreak were nearly indistinguishable showing that transmission of LA-MRSA can be detected by whole genome mapping. Finally, whole genome maps of LA-MRSA isolates originating from two unrelated veterinarians and their household members showed that veterinarians may carry and transmit different LA-MRSA strains at the same time. No such conclusions could be drawn based spa-typing and MLVA. Although PFGE seems to be suitable for molecular typing of LA-MRSA, WGM provides a much higher discriminatory power. Furthermore, whole genome mapping can provide a comparison with other maps within 2 days after the bacterial culture is received, making it suitable to investigate transmission events and outbreaks caused by LA-MRSA.

Whole-genome mapping: a new paradigm in strain-typing technology

Strain-typing technology in support of outbreak identification and resolution has evolved from phenotypic analysis, such as serology and biotypes, to much-more-robust molecular genetic approaches, such as pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing. Whole-genome mapping (WGM) has been recently applied to subtyping analysis, and it bridges the gap between PFGE (∼20 bands sorted by size) and whole-genome sequencing. WGM utilizes restriction site analysis but arranges 200 to 500 bands in the order they appear on the chromosome. WGM is able to quickly and cost-effectively generate high-resolution, ordered whole-genome maps of bacteria.

Identification and characterization of spontaneous deletions within the Sp11-Sp12 prophage region of Escherichia coli O157:H7 Sakai

Prophages make up 12% of the enterohemorrhagic Escherichia coli genome and play prominent roles in the evolution and virulence of this food-borne pathogen. Acquisition and loss of and rearrangements within prophage regions are the primary causes of differences in pulsed-field gel electrophoresis (PFGE) patterns among strains of E. coli O157:H7. Sp11 and Sp12 are two tandemly integrated and putatively defective prophages carried by E. coli O157:H7 strain Sakai. In this study, we identified 3 classes of deletions that occur within the Sp11-Sp12 region, at a frequency of ca. 7.74 × 10(-4). One deletion resulted in a precise excision of Sp11, and the other two spanned the junction of Sp11 and Sp12. All deletions resulted in shifts in the XbaI fragment pattern observed by PFGE. We sequenced the inducible prophage pool of Sakai but did not identify any mature phage particles corresponding to either Sp11 or Sp12. Deletions containing pchB and psrC, which are Sp11-carried genes encoding proteins known or suspected to regulate type III secretion, did not affect the secretion levels of the EspA or EspB effector. Alignment of the Sp11-Sp12 DNA sequence with its corresponding regions in other E. coli O157:H7 and O55:H7 strains suggested that homologous recombination rather than integrase-mediated excision is the mechanism behind these deletions. Therefore, this study provides a mechanism behind the previously observed genetic instability of this genomic region of E. coli O157:H7.

Carriage of stx2a differentiates clinical and bovine-biased strains of Escherichia coli O157

Background

Shiga toxin (Stx) are cardinal virulence factors of enterohemorrhagic E. coli O157:H7 (EHEC O157). The gene content and genomic insertion sites of Stx-associated bacteriophages differentiate clinical genotypes of EHEC O157 (CG, typical of clinical isolates) from bovine-biased genotypes (BBG, rarely identified among clinical isolates). This project was designed to identify bacteriophage-mediated differences that may affect the virulence of CG and BBG.

Methods

Stx-associated bacteriophage differences were identified by whole genome optical scans and characterized among >400 EHEC O157 clinical and cattle isolates by PCR.

Results

Optical restriction maps of BBG strains consistently differed from those of CG strains only in the chromosomal insertion sites of Stx2-associated bacteriophages. Multiplex PCRs (stx1, stx2a, and stx2c as well as Stx-associated bacteriophage - chromosomal insertion site junctions) revealed four CG and three BBG that accounted for >90% of isolates. All BBG contained stx2c and Stx2c-associated bacteriophage – sbcB junctions. All CG contained stx2a and Stx2a-associated bacteriophage junctions in wrbA or argW.

Conclusions

Presence or absence of stx2a (or another product encoded by the Stx2a-associated bacteriophage) is a parsimonious explanation for differential virulence of BBG and CG, as reflected in the distributions of these genotypes in humans and in the cattle reservoir.

Comparative whole-genome mapping to determine Staphylococcus aureus genome size, virulence motifs, and clonality

Despite being a clonal pathogen, Staphylococcus aureus continues to acquire virulence and antibiotic-resistant genes located on mobile genetic elements such as genomic islands, prophages, pathogenicity islands, and the staphylococcal chromosomal cassette mec (SCCmec) by horizontal gene transfer from other staphylococci. The potential virulence of a S. aureus strain is often determined by comparing its pulsed-field gel electrophoresis (PFGE) or multilocus sequence typing profiles to that of known epidemic or virulent clones and by PCR of the toxin genes. Whole-genome mapping (formerly optical mapping), which is a high-resolution ordered restriction mapping of a bacterial genome, is a relatively new genomic tool that allows comparative analysis across entire bacterial genomes to identify regions of genomic similarities and dissimilarities, including small and large insertions and deletions. We explored whether whole-genome maps (WGMs) of methicillin-resistant S. aureus (MRSA) could be used to predict the presence of methicillin resistance, SCCmec type, and Panton-Valentine leukocidin (PVL)-producing genes on an S. aureus genome. We determined the WGMs of 47 diverse clinical isolates of S. aureus, including well-characterized reference MRSA strains, and annotated the signature restriction pattern in SCCmec types, arginine catabolic mobile element (ACME), and PVL-carrying prophage, PhiSa2 or PhiSa2-like regions on the genome. WGMs of these isolates accurately characterized them as MRSA or methicillin-sensitive S. aureus based on the presence or absence of the SCCmec motif, ACME and the unique signature pattern for the prophage insertion that harbored the PVL genes. Susceptibility to methicillin resistance and the presence of mecA, SCCmec types, and PVL genes were confirmed by PCR. A WGM clustering approach was further able to discriminate isolates within the same PFGE clonal group. These results showed that WGMs could be used not only to genotype S. aureus but also to identify genetic motifs in MRSA that may predict virulence.

Assessment of whole-genome mapping in a well-defined outbreak of Salmonella enterica serotype Saintpaul

We investigated the use of whole-genome mapping and pulsed-field gel electrophoresis (PFGE) with isolates from an outbreak of Salmonella enterica serotype Saintpaul. PFGE and whole-genome mapping were concordant with 22 of 23 isolates. Whole-genome mapping is a viable alternative tool for the epidemiological analysis of Salmonella food-borne disease investigations.

Genomic anatomy of Escherichia coli O157:H7 outbreaks

The rapid emergence of Escherichia coli O157:H7 from an unknown strain in 1982 to the dominant hemorrhagic E. coli serotype in the United States and the cause of widespread outbreaks of human food-borne illness highlights a need to evaluate critically the extent to which genomic plasticity of this important enteric pathogen contributes to its pathogenic potential and its evolution as well as its adaptation in different ecological niches. Aimed at a better understanding of the evolution of the E. coli O157:H7 pathogenome, the present study presents the high-quality sequencing and comparative phylogenomic analysis of a comprehensive panel of 25 E. coli O157:H7 strains associated with three nearly simultaneous food-borne outbreaks of human disease in the United States. Here we present a population genetic analysis of more than 200 related strains recovered from patients, contaminated produce, and zoonotic sources. High-resolution phylogenomic approaches allow the dynamics of pathogenome evolution to be followed at a high level of phylogenetic accuracy and resolution. SNP discovery and study of genome architecture and prophage content identified numerous biomarkers to assess the extent of genetic diversity within a set of clinical and environmental strains. A total of 1,225 SNPs were identified in the present study and are now available for typing of the E. coli O157:H7 lineage. These data should prove useful for the development of a refined phylogenomic framework for forensic, diagnostic, and epidemiological studies to define better risk in response to novel and emerging E. coli O157:H7 resistance and virulence phenotypes.

Investigating the global genomic diversity of Escherichia coli using a multi-genome DNA microarray platform with novel gene prediction strategies

BACKGROUND

The gene content of a diverse group of 183 unique Escherichia coli and Shigella isolates was determined using the Affymetrix GeneChip® E. coli Genome 2.0 Array, originally designed for transcriptome analysis, as a genotyping tool. The probe set design utilized by this array provided the opportunity to determine the gene content of each strain very accurately and reliably. This array constitutes 10,112 independent genes representing four individual E. coli genomes, therefore providing the ability to survey genes of several different pathogen types. The entire ECOR collection, 80 EHEC-like isolates, and a diverse set of isolates from our FDA strain repository were included in our analysis.

RESULTS

From this study we were able to define sets of genes that correspond to, and therefore define, the EHEC pathogen type. Furthermore, our sampling of 63 unique strains of O157:H7 showed the ability of this array to discriminate between closely related strains. We found that individual strains of O157:H7 differed, on average, by 197 probe sets. Finally, we describe an analysis method that utilizes the power of the probe sets to determine accurately the presence/absence of each gene represented on this array.

CONCLUSIONS

These elements provide insights into understanding the microbial diversity that exists within extant E. coli populations. Moreover, these data demonstrate that this novel microarray-based analysis is a powerful tool in the field of molecular epidemiology and the newly emerging field of microbial forensics.

Host specific diversity in Lactobacillus johnsonii as evidenced by a major chromosomal inversion and phage resistance mechanisms

Genetic diversity and genomic rearrangements are a driving force in bacterial evolution and niche adaptation. We sequenced and annotated the genome of Lactobacillus johnsonii DPC6026, a strain isolated from the porcine intestinal tract. Although the genome of DPC6026 is similar in size (1.97mbp) and GC content (34.8%) to the sequenced human isolate L. johnsonii NCC 533, a large symmetrical inversion of approximately 750 kb differentiated the two strains. Comparative analysis among 12 other strains of L. johnsonii including 8 porcine, 3 human and 1 poultry isolate indicated that the genome architecture found in DPC6026 is more common within the species than that of NCC 533. Furthermore a number of unique features were annotated in DPC6026, some of which are likely to have been acquired by horizontal gene transfer (HGT) and contribute to protection against phage infection. A putative type III restriction-modification system was identified, as were novel Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) elements. Interestingly, these particular elements are not widely distributed among L. johnsonii strains. Taken together these data suggest intra-species genomic rearrangements and significant genetic diversity within the L. johnsonii species and indicate towards a host-specific divergence of L. johnsonii strains with respect to genome inversion and phage exposure.

Methods for genotyping verotoxin-producing Escherichia coli

Verotoxin-producing Escherichia coli (VTEC) is annually incriminated in more than 100,000 cases of enteric foodborne human disease and in losses amounting to $US 2.5 billion every year. A number of genotyping methods have been developed to track VTEC infections and determine diversity and evolutionary relationships among these microorganisms. These methods have facilitated monitoring and surveillance of foodborne VTEC outbreaks and early identification of outbreaks or clusters of outbreaks. Pulsed-field gel electrophoresis (PFGE) has been used extensively to track and differentiate VTEC because of its high discriminatory power, reproducibility and ease of standardization. Multiple-locus variable-number tandem-repeats analysis (MLVA) and microarrays are the latest genotyping methods that have been applied to discriminate VTEC. MLVA, a simpler and less expensive method, is proving to have a discriminatory power comparable to that of PFGE. Microarrays are successfully being applied to differentiate VTEC and make inferences on genome diversification. Novel methods that are being evaluated for subtyping VTEC include the detection of single nucleotide polymorphisms and optical mapping. This review discusses the principles, applications, advantages and disadvantages of genotyping methods that have been used to differentiate VTEC strains. These methods have been mainly used to differentiate strains of O157:H7 VTEC and to a lesser extent non-O157 VTEC.

Genomic signatures of strain selection and enhancement in Bacillus atrophaeus var. globigii, a historical biowarfare simulant

BACKGROUND

Despite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS).

RESULTS

Archival strains and two "present day" type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the "military" isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of "military" isolates on lactate-containing media, and showed that the "military" strains exhibited a hypersporulating phenotype.

CONCLUSIONS

Our analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.

Genome signatures of Escherichia coli O157:H7 isolates from the bovine host reservoir

Cattle comprise a main reservoir of Shiga toxin-producing Escherichia coli O157:H7 (STEC). The significant differences in host prevalence, transmissibility, and virulence phenotypes among strains from bovine and human sources are of major interest to the public health community and livestock industry. Genomic analysis revealed divergence into three lineages: lineage I and lineage I/II strains are commonly associated with human disease, while lineage II strains are overrepresented in the asymptomatic bovine host reservoir. Growing evidence suggests that genotypic differences between these lineages, such as polymorphisms in Shiga toxin subtypes and synergistically acting virulence factors, are correlated with phenotypic differences in virulence, host ecology, and epidemiology. To assess the genomic plasticity on a genome-wide scale, we have sequenced the whole genome of strain EC869, a bovine-associated E. coli O157:H7 isolate. Comparative phylogenomic analysis of this key isolate enabled us to place accurately bovine lineage II strains within the genetically homogenous E. coli O157:H7 clade. Identification of polymorphic loci that are anchored both in the chromosomal backbone and horizontally acquired regions allowed us to associate bovine genotypes with altered virulence phenotypes and host prevalence. This study catalogued numerous novel lineage II-specific genome signatures, some of which appear to be associated intimately with the altered pathogenic potential and niche adaptation within the bovine rumen. The presented extended list of polymorphic markers is valuable in the development of a robust typing system critical for forensic, diagnostic, and epidemiological studies of this emerging human pathogen.

Use of optical mapping to sort uropathogenic Escherichia coli strains into distinct subgroups

Optical maps were generated for 33 uropathogenic Escherichia coli (UPEC) isolates. For individual genomes, the NcoI restriction fragments aligned into a unique chromosome map for each individual isolate, which was then compared with the in silico restriction maps of all of the sequenced E. coli and Shigella strains. All of the UPEC isolates clustered separately from the Shigella strains as well as the laboratory and enterohaemorrhagic E. coli strains. Moreover, the individual strains appeared to cluster into distinct subgroups based on the dendrogram analyses. Phylogenetic grouping of these 33 strains showed that 32/33 were the B2 subgroup and 1/33 was subgroup A. To further characterize the similarities and differences among the 33 isolates, pathogenicity island (PAI), haemolysin and virulence gene comparisons were performed. A strong correlation was observed between individual subgroups and virulence factor genes as well as haemolysis activity. Furthermore, there was considerable conservation of sequenced-strain PAIs in the specific subgroups. Strains with different antibiotic-resistance patterns also appeared to sort into separate subgroups. Thus, the optical maps distinguished the UPEC strains from other E. coli strains and further subdivided the strains into distinct subgroups. This optical mapping procedure holds promise as an alternative way to subgroup all E. coli strains, including those involved in infections outside of the intestinal tract and epidemic strains with distinct patterns of antibiotic resistance.

Optical mapping reveals a large genetic inversion between two methicillin-resistant Staphylococcus aureus strains

Staphylococcus aureus is a highly versatile and evolving bacterium of great clinical importance. S. aureus can evolve by acquiring single nucleotide polymorphisms and mobile genetic elements and by recombination events. Identification and location of novel genomic elements in a bacterial genome are not straightforward, unless the whole genome is sequenced. Optical mapping is a new tool that creates a high-resolution, in situ ordered restriction map of a bacterial genome. These maps can be used to determine genomic organization and perform comparative genomics to identify genomic rearrangements, such as insertions, deletions, duplications, and inversions, compared to an in silico (virtual) restriction map of a known genome sequence. Using this technology, we report here the identification, approximate location, and characterization of a genetic inversion of approximately 500 kb of a DNA element between the NRS387 (USA800) and FPR3757 (USA300) strains. The presence of the inversion and location of its junction sites were confirmed by site-specific PCR and sequencing. At both the left and right junction sites in NRS387, an IS1181 element and a 73-bp sequence were identified as inverted repeats, which could explain the possible mechanism of the inversion event.

Genomic instability in regions adjacent to a highly conserved pch prophage in Escherichia coli O157:H7 generates diversity in expression patterns of the LEE pathogenicity island

The LEE pathogenicity island has been acquired on multiple occasions within the different lineages of enteropathogenic and enterohemorrhagic Escherichia coli. In each lineage, LEE expression is regulated by complex networks of pathways, including core pathways shared by all lineages and lineage-specific pathways. Within the O157:H7 lineage of enterohemorrhagic E. coli, strain-to-strain variation in LEE expression has been observed, implying that expression patterns can diversify even within highly related subpopulations. Using comparative genomics of E. coli O157:H7 subpopulations, we have identified one source of strain-level variation affecting LEE expression. The variation occurs in prophage-dense regions of the genome that lie immediately adjacent to the late regions of the pch prophage carrying pchA, pchB, pchC, and a newly identified pch gene, pchX. Genomic segments extending from the holin S region to the pchA, pchB, pchC, and pchX genes of their respective prophage are highly conserved but are nonetheless embedded within adjacent genomic segments that are extraordinarily variable, termed pch adjacent genomic regions (pch AGR). Despite the remarkable degree of variation, the pattern of variation in pch AGR is highly correlated with the distribution of phylogenetic markers on the backbone of the genome. Quantitative analysis of transcription from the LEE1 promoter further revealed that variation in the pch AGR has substantial effects on absolute levels and patterns of LEE1 transcription. Variation in the pch AGR therefore serves as a mechanism to diversify LEE expression patterns, and the lineage-specific pattern of pch AGR variation could ultimately influence ecological or virulence characteristics of subpopulations within each lineage.

Optical mapping and 454 sequencing of Escherichia coli O157 : H7 isolates linked to the US 2006 spinach-associated outbreak

Optical maps for five representative clinical, food-borne and bovine-derived isolates from the 2006 Escherichia coli O157 : H7 outbreak linked to fresh spinach in the United States showed a common set of 14 distinct chromosomal markers that define the outbreak strain. Partial 454 DNA sequencing was used to characterize the optically mapped chromosomal markers. The markers included insertions, deletions, substitutions and a simple single nucleotide polymorphism creating a BamHI site. The Shiga toxin gene profile of the spinach-associated outbreak isolates (stx1(-) stx2(+) stx2c(+)) correlated with prophage insertions different from those in the prototypical EDL933 and Sakai reference strains (stx1(+) stx2(+) stx2c(-)). The prophage occupying the yehV chromosomal position in the spinach-associated outbreak isolates was similar to the stx1(+) EDL933 cryptic prophage V, but it lacked the stx1 gene. In EDL933, the stx2 genes are within prophage BP933-W at the wrbA chromosomal locus; this locus was unoccupied in the spinach outbreak isolates. Instead, the stx2 genes were found within a chimeric BP933-W-like prophage with a different integrase, inserted at the argW locus in the outbreak isolates. An extra set of Shiga toxin genes, stx2c, was found in the outbreak isolates within a prophage integrated at the sbcB locus. The optical maps of two additional clinical isolates from the outbreak showed a single, different prophage variation in each, suggesting that changes occurred in the source strain during the course of this widespread, multi-state outbreak.

The pangenome structure of Escherichia coli: comparative genomic analysis of E. coli commensal and pathogenic isolates

Whole-genome sequencing has been skewed toward bacterial pathogens as a consequence of the prioritization of medical and veterinary diseases. However, it is becoming clear that in order to accurately measure genetic variation within and between pathogenic groups, multiple isolates, as well as commensal species, must be sequenced. This study examined the pangenomic content of Escherichia coli. Six distinct E. coli pathovars can be distinguished using molecular or phenotypic markers, but only two of the six pathovars have been subjected to any genome sequencing previously. Thus, this report provides a seminal description of the genomic contents and unique features of three unsequenced pathovars, enterotoxigenic E. coli, enteropathogenic E. coli, and enteroaggregative E. coli. We also determined the first genome sequence of a human commensal E. coli isolate, E. coli HS, which will undoubtedly provide a new baseline from which workers can examine the evolution of pathogenic E. coli. Comparison of 17 E. coli genomes, 8 of which are new, resulted in identification of approximately 2,200 genes conserved in all isolates. We were also able to identify genes that were isolate and pathovar specific. Fewer pathovar-specific genes were identified than anticipated, suggesting that each isolate may have independently developed virulence capabilities. Pangenome calculations indicate that E. coli genomic diversity represents an open pangenome model containing a reservoir of more than 13,000 genes, many of which may be uncharacterized but important virulence factors. This comparative study of the species E. coli, while descriptive, should provide the basis for future functional work on this important group of pathogens.