| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
Division of Molecular Biology, Office of Applied Research and Safety Assessment, Center for Food Safety and Applied Nutrition, US Food and Drug Administration, Laurel, MD 20708, USA
Correspondence
Thomas A. Cebula
Thomas.Cebula{at}fda.hhs.gov
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). E. coli O157 : H7 infection causes severe abdominal cramping and a watery diarrhoea that progresses to severe bloody diarrhoea (Riley et al., 1983; Griffin et al., 1988). An estimated 15–20 % of people infected with E. coli O157 : H7 present with indications severe enough to require hospitalization and in about 2–7 % of patients, infections progress from haemorrhagic colitis into haemolytic-uraemic syndrome (HUS) with kidney failure complications (Su & Brandt, 1995).
E. coli O157 : H7 infections usually have a food-borne aetiology. Since cattle are a common reservoir of this organism, the vast majority of outbreaks since identification of the pathogen in 1982 have been traced to consumption of meat and dairy products contaminated with this organism (Mead & Griffin, 1998). Recent outbreaks of human disease, in contrast, have been linked to consumption of fresh produce tainted with E. coli O157 : H7. A recent multi-state outbreak in the USA, for instance, was associated with contaminated spinach and caused 199 cases of illness, including three deaths. Among the ill, 51 % were hospitalized and in 16 % of the cases, infection progressed to HUS and kidney failure (www.cdc.gov/foodborne/ecolispinach/100606.htm). The large number of patients hospitalized and high rate of kidney failure suggest that this outbreak was due to a more virulent strain of E. coli O157 : H7. Assessing genetic variation among temporally and geographically diverse strains of E. coli O157 : H7 thus becomes particularly important for understanding its public health consequences and adaptation to new environmental niches, as well as for the practical aspects of surveillance and epidemiological trace-back.
For the epidemiological tracing of clinical isolates to source strains, PFGE has been an effective tool for molecular analysis. Detecting variation among E. coli O157 : H7 strains at the gene and chromosome level, however, requires a more extensive analysis of strains. In the recent past, variation or diversity within bacterial species has been defined largely by genetic polymorphisms – early on, by analysing variants at the protein level as revealed by multilocus enzyme electrophoresis (MLEE, Whittam et al., 1993), and, more recently, by examining variation at the nucleotide level through the use of multilocus sequence typing (MLST, Feil & Enright, 2004; Urwin & Maiden, 2003). Studies have shown differences in the ability to differentiate individual strains by MLST and PFGE (Nemoy et al., 2005). Examination of a small number of genes has revealed only limited variation at the nucleotide level among E. coli O157 : H7 strains (Foley et al., 2004). Sampling larger portions of the genome with a variety of rapid microarray-based techniques, on the other hand, has provided ample evidence for significant allelic diversity among E. coli O157 : H7 strains. For example, when 4000 PCR-amplified genes (Fukiya et al., 2004) or 6000 50-mer oligonucleotide gene sequences (Wick et al., 2005) were targeted, and when 1199 specific genes (Zhang et al., 2006) or 1 % of the bacterial genome (Jackson et al., 2006) were examined, rich allelic diversity was revealed among a relatively homogeneous group of pathogenic E. coli strains.
Diversity among E. coli serotypes, at the chromosomal level, likewise was uncovered by direct genomic comparisons. For example, relative to the 4.639 Mbp genomic sequence of the prototypic E. coli K-12 strain, MG1655, the 5.528 Mbp genome of the pathogenic E. coli O157 : H7 strain, EDL933, was ‘peppered’ with insertions and deletions (Perna et al., 2001). That is, 177 segments of O157 : H7 DNA were defined as ‘O-islands’ because they were found in the O157 : H7 genome, but not in K-12. Correspondingly, 234 ‘K-islands’ were identified by their presence in the MG1655 genome and absence in the EDL933 genome. When the 5.231 Mbp sequence of E. coli uropathogenic strain CFT073 was added to this comparison, the complex mosaic structure of E. coli pathogens became even more evident (Welch et al., 2002). While about 39 % of the non-redundant proteins were common to all three strains, 46 % were unique to a single strain. The two E. coli pathogens differed from each other in 900 000 bp of unique sequence, and thus differed from each other by as much as each differed from the commensal strain. Despite the remarkable synteny of enteric bacterial chromosomes, as typified by E. coli and Salmonella, genomic comparisons have identified the movement and repositioning of large blocks or segments of genes as a major source of diversity. Comparisons of completed sequenced genomes using the NCBI gMap web site underscore the findings that related bacterial genomes contain multi-segmented arrays of genes inserted in different orientations and shuffled in chromosomal positions (http://www.ncbi.nlm.nih.gov/sutils/gmap.cgi; Kotewicz et al., 2003).
Optical mapping, which can scan and assess the architecture of complete bacterial genomes, thus should prove useful for comparative genomics as well as for epidemiology and microbial forensics (Cai et al., 1998). Optical mapping was pivotal in closing the genome sequence of E. coli O157 : H7 EDL933 (Lim et al., 2001), and its value was further underscored when optical maps were used to compare particular subtypes of Shigella with a sequenced strain (Zhou et al., 2004). In conjunction with subtractive hybridization to identify unique sequences, optical mapping also was used to map the sequences onto the enterotoxigenic E. coli H10407 chromosome (Chen et al., 2006). Herein, the chromosomal organization of each of 11 E. coli O157 : H7 strains was examined by BamHI optical mapping. A BamHI optical map for a typical 5–5.5 Mbp E. coli O157 : H7 chromosome presents as 500–700 ordered restriction fragments, with the number and arrangement of fragments varying with each strain examined. As restriction digestion is carried out on genomic DNA affixed to a glass substratum, the molecular size of contiguous BamHI fragments is determined and precisely mapped along the chromosome. In contrast, PFGE, used commonly in molecular epidemiology studies (Ribot et al., 2006), typically produces 50 genomic restriction fragments ordered only by fragment size within the gel, not by map position on the genome. Because of its greater resolution and specificity for strain identification, optical mapping holds promise for both the public health and forensic communities (Cebula et al., 2005). To this end, we present data to show that complex events, including inversions, deletions and insertions occurring within an individual chromosome, can be detected using optical mapping. Optical maps can serve as a unique DNA ‘bar code’ identifier for a particular microbe, not only at the species or subspecies level, but also at the individual strain or isolate level.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
. Reference strains that have been completely sequenced were EDL933 (GenBank AE005174), isolated from ground hamburger associated with a 1982 outbreak in the United States; and Sakai, RIMD 0509952 (GenBank BA000007), a clinical sample isolated during an outbreak in Japan in 1996. Selection of E. coli O157 : H7 strains for optical mapping was based on cladistic analysis of single nucleotide polymorphisms, which allowed grouping of the strains into different clades (Cebula et al., 2005). The FDA isolates of EDL933 and Sakai were assigned strain numbers EC1275 and EC1276, respectively. To clarify the comparisons of the optical maps of these isolates with the in silico maps derived from sequenced EDL933 and Sakai isolates, the EC1275 and EC1276 designations have been retained in this report.
|
). |
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). The accuracy of optical mapping was demonstrated by comparing the in silico BamHI map and the optical BamHI map of EDL933 (EC1275, Fig. 1b). With the exception of the loss of a number of small fragments (dropout fragments, discussed in the next section) the optical map recapitulates the sequence-based map for both EDL933 and Sakai strains.
|
). Of the 96 ‘dropout’ fragments, all were less than 2000 bp in size; a majority (60) was smaller than 600 bp. That is, of the fragments expected from the sequence of Sakai, 60/60 fragments in the 21–600 bp range, 22/30 fragments in the 601–1000 bp range, 10/39 fragments in the 1–1.5 kb range and 4/36 in the 1.5–2.0 kb range were not detected in the Sakai optical map (EC1276, Table 2). A similar set of results was found for EDL933. There were 112 dropout fragments in the optical map of EDL933 (EC1275, Table 2) compared to the in silico map. None of the 61 expected fragments in the 21–600 bp range were detected. Just as was the case for the optical map of Sakai, there was an increase in the mapping of expected fragments as they increased in size (Table 2). Optical mapping detected 100 % of expected fragments above 2 kb.
|
). This area of the chromosome allowed a reasonable comparison of the dropouts from strain to strain. Of the 110 restriction fragments in this portion of the genome, 25 showed dropouts among the nine strains; eight dropouts were below 600 bp in size, nine were 600–1000 bp and eight were from 1001 to 2000 bp. Thus, as was found for the EDL933 and Sakai optical maps relative to their in silico maps, no fragments below 600 bp were detected in the nine optical maps of unsequenced strains. In the optical maps of the nine unsequenced strains, five demonstrated all eight expected fragments in the 1–2 kb size range. Each of the remaining four maps was missing either one, two, three, or four of the expected eight fragments. Overall, for 25 expected fragments, the mean loss was 17.5; expected fragment loss ranged from 16 to 21 fragments among the nine strains. For the 1 Mbp portion of the genomes examined, this represented a mean loss of 10.9 kb, ranging from 7.5 to 15.8 kb among the nine strains. Based on fragments expected from the sequenced genomes, the optical maps showed an estimated loss of 44 kb for each genome. For the two sequenced reference strains, the optical maps showed 112 dropout fragments totalling 35.014 kb for EDL933 (EC1275) and 96 dropout fragments totalling 20.309 kb for Sakai (EC1276) (Table 2
).
|
There were a number of differences between the optical and the in silico maps of EDL933. In EDL933, if the 112 small drop-out fragments (21 bp to 2 kb) are subtracted from the sequence-based total of 645 fragments, there were 533 in silico fragments relative to the 525 optical fragments. The difference of 12 was a combination of BamHI polymorphisms (or failure to digest at six positions), which created a net loss of six fragments by fusion of fragments, and two and four fragment differences each in cryptic prophages O and V in runs of small fragments. These same types of similarities and differences held true for the Sakai sequence-based and optical maps. Determining whether sequencing errors or recombination between and among duplicated prophages would explain the seven differences awaits further investigation. For Sakai, the identity between the optical and in silico maps was stronger. Disregarding the 96 dropout fragments, there was only a single difference between the optical (EC1276) and in silico chromosome maps of Sakai. The difference represented a probable single restriction site change; the remaining 540 out of 541 restriction fragments were comparable in size and position in the optical and in silico maps of Sakai.
Chromosome size
Relative to the sequence-based fragment sizes, the standard deviation of measuring a 10 kb fragment in optical maps of 11 strains was about 300 bp or 3 %. The standard deviation for 20–30 kb fragments was typically 600 bp or 2 %. The variability of the total chromosome length derived from the sum of 530 measured fragments becomes a mean of means with a variability (v) that can be approximated by the formula v=
(530x0.22)=4.6 kb or about 0.1 % for a 5 Mbp genome. The estimated effect of fragment measurement accuracy (5 kb) was smaller than the calculated effect of dropout fragments (30–40 kb) on the overall length determination for a 5 Mbp genome. Whereas the reference sequence of the Sakai strain was 5 498 451 bp, the optical map size was 5 514 226 bp. For the EDL933 strain, the reference sequence was 5 528 446 bp, and the optical map was 5 535 378 bp. The optical map chromosome lengths for the two reference strains were larger than the sequence-based chromosome sizes by 16 and 7 kb. The chromosomal sizes measured for the nine additional strains of E. coli O157 : H7 ranged from 5.326 Mbp (EC536) to 5.579 Mbp (AB1), a size difference of 250 kb (Table 1).
These data, coupled with the comparison of each of the individual BamHI maps, demonstrated the rich diversity that exists among E. coli O157 : H7 strains.
The optical maps visually demonstrated chromosomal similarities and differences. Optical map alignment software semiquantitatively defines the similarities between genomes. In this study of 11 E. coli O157 : H7 strains, restriction fragment alignments were used to define similarities. Because of the extensive annotation of the sequenced EDL933 and Sakai sequences, similarities and differences identified in the sequences were verified in the optical maps of the two reference strains and, by inference, in the optical maps of the other isolates. All differences in the sequence-based maps were seen in the optical maps of the reference strains. Verification of chromosome similarities by Southern analyses was not performed; however, the presence of the normal genetic complement of E. coli O157 : H7 (some 5000 genes) was verified for each isolate by DNA microarray analysis (data not shown). Several deletions demonstrated in optical maps were verified by DNA microarray analysis (data not shown).
Fig. 2 shows sets of pairwise alignments of the BamHI optical maps of the strains examined in this study to reference in silico maps. The alignments demonstrate that although all of the strains showed large regions of similarity by contig analysis (green fragments), extensive regions of dissimilarity (white fragments) as well as five notable chromosomal inversions were found (yellow fragments, Fig. 2). Among the 11 E. coli O157 : H7 optical maps, there were 91 differences found at 28 chromosome sites or loci. The polymorphic sites were numbered consecutively from the chromosome origin; they were defined both by the closest numbered BamHI reference fragment and by the position of the fragment on the EDL933 chromosome (see Supplementary Table S1, available with the online version of this paper).
Simple fragment size polymorphisms
Comparison of the sequence-based EDL933 (645 fragments) and Sakai (639 fragments) in silico maps showed five restriction site polymorphisms and 13 simple fragment differences greater than 100 bp between equivalently positioned fragments. Most of these polymorphisms represented changes in fragment sizes below the limit of resolution and below the 2 kb limit of reliability. Only two of the differences between EDL933 and Sakai were useful for differentiating individual strains, polymorphic sites 21 and 25 (Supplementary Table S1). Polymorphic site 21 showed a difference between EDL933 restriction fragment e374 (2589 bp) and the comparable Sakai fragment s367 (3902 bp). The optical maps of EDL933 and Sakai showed fragments of 2111 bp and 4263 bp, respectively, in the two strains. Analysis of optical maps of the nine other strains revealed two classes of fragments at this chromosomal position; five strains harboured EDL933-like fragments of 2544±178 bp and four contained Sakai-like fragments of 4077±181 bp. The sequence-based polymorphic site 25 showed an EDL933 in silico fragment e600 of 20.654 kb and a comparable Sakai in silico fragment s593 of 59.412 kb. Optical maps of the nine remaining O157 : H7 strains showed a 20.652±0.377 kb EDL933-like fragment at this position. The number of useful simple fragment polymorphisms was too small to allow significant strain differentiation.
Two additional simple polymorphisms at sites 24 and 28 were found when the optical maps of the remaining nine strains were compared (Supplementary Table S1). These two simple discriminatory sites contained RFLPs. As detailed in the next section, although some RFLPs initially appeared simple, examination showed that some RFLPs were indicative of polymorphic sites where complex chromosomal changes and rearrangements had occurred. They were evident only when a greater number of strains was examined. Notably, many of the complex polymorphic sites were within prophages undergoing substitutions, inversions and possible truncations. The number and variability of these complex polymorphisms allowed strain differentiation.
Complex chromosomal polymorphisms; inversions
Complex alignment differences were frequent in the optical maps of the E. coli O157 : H7 strains and they were used successfully to discriminate individual strains. Inversions earmarked the first class of complex differences. The original genomic sequences and the in silico maps of EDL933 and Sakai strains demonstrated a chromosomal inversion spanning 430 kb, or 7.8 % of EDL933 relative to the Sakai and E. coli K-12 chromosomes. This inversion was prominent in the aligned in silico maps of EDL933 and Sakai (Fig. 2a). This inversion is flanked by prophages O and P at EDL933 O-islands 57 and 71.
Notably, when nine optical maps were compared with Sakai and EDL933, five showed chromosomal inversions (Fig. 2b). In each case, fragments corresponding to annotated prophages were found flanking the chromosomal inversions (Table 1). A 7.7 % inversion, distinctly different from, but similar in size and location to the 7.8 % inversion found in EDL933, was mapped in EC1225. Although differences between similar inversions might be due to assorted phages lysogenizing at the same insertion sites prior to or after inversion, they could also be explained by independent recombination events occurring at different sites within homologous regions of the flanking prophages, thus creating discrete configurations of restriction fragments at the junctions.
There were four other inversions. The chromosomes of the four strains appear complex when their optical maps are aligned with Sakai, but appear simple when aligned with EDL933. The chromosome structures of these strains are the result of double inversions. This is shown schematically for EC502 in Fig. 3(a). The order of inversion is not critical, but the structures observed in the optical maps can be most simply explained by an inversion occurring within another inversion (Fig. 3a). Fig. 3(b) illustrates this inversion within an inversion with the optical map alignments of EC502 with the Sakai and EDL933 chromosome maps. For EC502, the large inversion extended from the EDL933-equivalent position of O-island 45, which contains prophage W in EDL933, to the equivalent position of EDL933 O-island 93, the site for prophage V. For EC502, in the alignment with Sakai, within the large inverted segment (yellow), the internal inversion (green) has regained its orientation and alignment to Sakai. When shown aligned with EDL933, the entire inversion in EC502 is yellow, since EDL933 contains the internal inversion. For the other large inversion in strain EC533, the overall structure is very similar, except that the inversion extended from the equivalent positions of EDL933 prophage M (O-island 44) to prophage U (O-island 79). The two large inversions represented 1.656 kb (30 %) of the genome in EC502 and 1.482 kb (27 %) in EC533. These large inversions encompass the smaller 429 kb EDL933 inversion.
|
Complex chromosomal differences; substitutions, insertions and deletions
Optical mapping alignments of the various strains also showed that portions of the chromosome were not aligned over stretches up to 150 kb (strain EC1231, Fig. 2a, strain AB1, Fig. 2b). Many of the non-aligned fragments occurred at positions known in EDL933 to contain cryptic prophages. Supplementary Table S1 summarizes the prophage deletions and substitutions at novel chromosomal positions, including a number of novel insertions whose size and fragment distributions are suggestive of known prophages. Comparing each of the strains with the in silico EDL933 map, the extent of each change was measured and positioned relative to the EDL933 sequence. A total of 91 chromosomal mapping incongruities were identified at 28 unique loci within the 11 strains examined (Supplementary Table S1). Relative to strain EDL933, each strain exhibited chromosomal variation at seven to 18 sites, yielding a unique map with a different assortment of polymorphisms, including the five distinct chromosomal inversions. Two of the most variable sites are those occupied by the stx1 and stx2 prophages in EDL933, prophage W (stx2) at O-island 45 and prophage V (stx1) at O-island 93. Examples of each of the major types of variation are illustrated and detailed in the next section for strain EC536.
Complex chromosomal differences in strain EC536
The major chromosomal polymorphisms identified by the optical map of EC536 are shown in Fig. 4. An insertion containing five BamHI restriction fragments, totalling 26.9 kb, occurs in EC536 within a set of fragments corresponding to O-island 8 (Fig. 4a). In EDL933, O-island 8 contains two cryptic prophages H and I, spanning bp 300 060 to 323 540. Two chromosomal polymorphisms in EC536 involved loss of prophages relative to their chromosomal positions in EDL933. Whereas EDL933 harbours a P4-like prophage sequence at O-island 43 (bp 1 060 000) and prophage W at O-island 45 (bp 1 310 000), optical mapping demonstrated that EC536 was devoid of prophages at these sites (Fig. 4a). The absence of a prophage at one chromosomal location does not mean, however, that the prophage or a related prophage is not located elsewhere on the chromosome. An example is EDL933 prophage W, which carries the Shiga-like toxin II gene, stx2, at the chromosomal position annotated O-island 45 in EDL933. Although strain EC536 contains the stx2 gene as determined by PCR using stx2-specific primers, as well as by DNA microarray analysis, microarray results suggest that a number of prophage W genes in strain EC536 are absent (S. Jackson, I. Patel & M. L. Kotewicz, unpublished results). Since the O-island 45 site of prophage W is unoccupied in EC536 (Fig. 4a), the stx2 gene is probably located at a different chromosomal location. Five candidate sites were identified from the optical map: the insertion in the H/I prophage described above; the substitution within O-island complex 50, 51 and 52; both prophages, O and R, flanking the chromosomal inversion found in EC536; the phage-like insertion near O-island 79; and the altered prophage at O-island 93 (Fig. 4b). The latter two sites, chromosomal polymorphisms 18 and 20, represent prime candidates for the new location of stx2 in EC536. At site 18, a 55.7 kb insertion of six additional fragments, located 39.3 kb upstream of the O-island 79 site in EDL933, resembles the stx1 prophage V found in EDL933. At site 20, the BamHI polymorphism is found within fragments representing O-island 93, prophage CP-933 V containing the stx1 gene in EDL933. The stx2 gene found in EC536 may have replaced the stx1 gene.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
).
Many of the chromosomal differences found in the optical maps of nine E. coli O157 : H7 strains relative to two reference strains were associated with known prophages. There are 22 entries for integrases in the annotated EDL933 genome (Perna et al., 2001; see GenBank AE005174). Only one fully competent prophage, BP-933W, is found in EDL933 (Plunkett et al., 1999), indicating, by the measure of integrase genes, that 21 defective prophages are present in the chromosome. The other integrase genes are present within truncated phage genomes, presumed to be defective because of missing blocks of genes and their inability to form complete phage particles when induced. The number of defective prophages annotated within the Sakai and EDL933 genomes suggests that since the time their intact progenitors took up residence in each pathogen, extensive editing by deletion has occurred to ‘lock in’ prophage remnants in different E. coli O157 : H7 strains. This is consistent with hypotheses about mechanisms to prevent the continual expansion of bacterial genomes (Ochman, 2005).
Marked heterogeneity among prophages and their integration sites in different E. coli O157 : H7 strains has been observed by genomic PCR scanning (Ohnishi et al., 2002), and further refined by combining PCR scanning with microarray analysis for adjudging presence or absence of specific gene sequences (Ogura et al., 2006). Although both optical mapping and DNA microarray interrogation of novel strains rely heavily on knowledge of previously sequenced genomes, the latter technique is incapable of detecting novel insertions. In contrast, as evidenced from the present study, this is a major strength of optical mapping. Optical mapping complements PCR- and cloning-based strategies for genome sequencing, facilitating closure of bacterial genomes like E. coli O157 : H7 when multiple copies of homologous phages might hinder sequence assembly (Perna et al., 2001).
The optical maps of nine E. coli O157 : H7 strains presented in the current study contained 91 differences at 28 polymorphic chromosomal loci relative to the EDL933 and Sakai sequence-based maps, and relative to each other. Three of the 28 polymorphic loci showed simple changes, such as the gain or loss of a BamHI site, but 25 sites contained more complex differences. Fifteen of the loci demonstrating complex changes in the optical maps were O-islands previously annotated as prophage or cryptic prophage sequences, and seven others are candidates for novel prophage insertion sites. Moreover, at several of the latter sites, the optical maps showed positions where a large insertion had occurred in one strain, whereas smaller insertions were found in other strains. For example, at polymorphic site 18, EC1231 contained a 97 kb insertion relative to EDL933; other strains, however, contained smaller insertions at this position, such as EC869 (61 kb), EC536 (56 kb) and a number of strains with 6 kb insertions. This suggests that phage acquisition and truncation occurred at these sites (Supplementary Table S1, polymorphic loci 18, 22 and 26). No other means has yet been used to demonstrate that the shorter insertions are truncated phages, and it is possible that smaller insertions have occurred at these sites rather than the truncation of large prophages. Whether these sites result from strain-specific scars of phage insertion and truncation, or other insertions, determining the likely origins of these changes through sequence analysis will be of interest.
In addition to altering the primary structure, multiple homologous prophages have another effect on chromosome structure. Five chromosomal inversions were found among the nine strains examined. The optical maps revealed the sizes and locations of each inversion, and in each case the inversion was flanked by prophages. The positions of the ends of the inversions suggested that prophages served as targets for recombination, as has been found for inversions between artificially constructed chromosomal repeats (Miesel et al., 1994) or inversions found between copies of rRNA operons (Liu & Sanderson, 1995) and insertion elements in Salmonella (Alokam et al., 2002). The large number of cryptic prophages found in EDL933 and Sakai supply extensive 10–90 kb regions for homologous recombination. The inversion mapped for EC1225 was similar to the EDL933 inversion, and was flanked by prophages O and P (EDL933 O-islands 57 and 71). Careful examination of the prophage restriction fragments flanking the inversion in EC1225 suggested that either this inversion occurred within different restriction fragments than occurred in EDL933, i.e. it was a different inversion, or that different but related prophages were resident prior or subsequent to the inversion.
In the other four strains carrying inversions, AB1, EC536, EC533 and EC502, alignment of each optical map with the in silico chromosomal map of Sakai suggested a complex set of events over hundreds of kb. On the other hand, the alignments of these strains with the in silico EDL933 map suggested that, for all four strains, two simple chromosomal inversions in each strain might explain the chromosomal structures. For EC533 and EC502, two large inversions (26.8 and 29.9 % of the chromosome) flanking an EDL933-like inversion (7.8 % of the genome) explained their chromosomal maps. This does not exclude the reverse order of occurrence of a double inversion, that is, a novel large inversion followed by an EDL933-like inversion. The double inversions were flanked by prophage pairs, M/U or W/V (Table 1). For AB1 and EC536, a small 4.7 % inversion and a 4.2 % inversion within the larger 7.8 % EDL933-like inversion most simply explained their chromosomes. The prophage pair O and R was involved in these smaller inversions. As the flanking restriction fragments were different in EC536 and AB1, it is possible that they were independent inversions. Each of the four inversions is described most simply as an inversion within a chromosome containing an EDL933-like inversion, i.e. a double inversion.
There was concern about the preponderance of inversions found in the optical maps of these strains. Was it possible that the optical mapping assembly software created an artefact, inversions based on the assembly of similar flanking prophages? Two points refute this explanation. The prophages range in size from 10 to 90 kb. In contrast, the single molecule ‘reads’ performed in optical mapping are much longer. The minimum size of each molecule read is 200 kb, the mean size of the assembled molecules is 350 kb, and it is common to detect single molecules that extend over 600 kb. Additionally, optical maps are assembled to very high coverage; the minimal coverage was 30x and mean coverage was 50–100x. As a consequence, the overlap between molecules in the assembly is also very high, so that when two molecules are aligned, there is at least 150 kb of shared regions. In this context, even a relatively large prophage of 90 kb would be contained entirely within a single molecule. In order to overlap two molecules, the prophage would have to be flanked with at least 50 kb of additional information. Because of the length of the single mapped molecules, the accuracy of optical mapping is extremely robust in the presence of prophages, even when a similar prophage is found in multiple locations within a single genome.
The second point arguing against the chromosomal inversions being an artefact was independent PFGE and Southern experiments published during the preparation of this manuscript. Carefully detailed PFGE analysis of inversions occurring during passage of E. coli O157 : H7 EDL933 in the laboratory (Iguchi et al., 2006) found most of the inversions identified by optical mapping. An additional inversion between an EDL933 P4-like prophage M located at O-island 43 and another P4-like prophage at O-island 48 was also found after laboratory passage of EDL933. The large inversion in the optical map of strain EC533, representing 27 % of the EDL933 chromosome, was not found in the laboratory passage studies. Both optical mapping and PFGE results suggest that as additional strains are examined, additional inversions between other pairs of homologous phages will probably be found. There are, of course, chromosomal constraints that will modulate inversion formation, e.g. the positions of the origin and termination of chromosome replication (Iguchi et al., 2006).
In addition to showing that a single isolate of EDL933, during subculture in the laboratory, produced a number of different PFGE profiles derived from chromosomal inversions (Iguchi et al., 2006), three other E. coli O157 : H7 strains, including Sakai, showed changes in PFGE patterns after repeated subculture and prolonged storage (Shima et al., 2006). Chromosomal inversions in these strains were implicated in some of the changed PFGE profiles. Notably, a strain missing the stx1 prophage exhibited fewer changes in its PFGE profile than did strains carrying both stx1 and stx2 phages.
Overall, the optical map data presented here and other evidence (Iguchi et al., 2006; Shima et al., 2006) suggest that chromosomal inversions are ongoing and frequent in O157 : H7 strains. Although the frequency of inversions cannot be estimated from optical maps of temporally and geographically dispersed strains, the five sets of inversions in EDL933 identified in the laboratory studies arose during 20, 30 and 50 subcultures (Iguchi et al., 2006). Regarding the frequency of inversions and prophage changes in chromosomes, dissimilarities in the stx1 and stx2 gene profiles of particular O157 : H7 strains have been demonstrated to be not only the result of the heterogeneity in phages containing the toxin genes but also due to the unexpected variability in their integration sites (Shaikh & Tarr, 2003). Moreover, polymorphic PCR amplification patterns were an early indicator that insertions and deletions dominated the differences between O157 : H7 strains (Kudva et al., 2002a, b); many of these polymorphisms were found within prophages. The optical map data confirm earlier suggestions that polymorphisms in O157 are not sequence-based changes in XbaI or BamHI sites, but rather that PFGE polymorphisms reside within changing prophage and chromosome profiles.
Thus, chromosomal inversions and phage trafficking are major contributors to outbreak-specific PFGE profiles. The real-world exposure of an O157 strain to environmental stress, to host selection and to teeming populations of competing phages alters the chromosomes of this pathogen in a feedlot-by-feedlot, patient-by-patient and temporally and geographically specific manner. The scars of prophage deletions represent locked-in variants that accumulate and which will be documented over the next few years by sequencing of E. coli O157 : H7 strains. The first of these chromosomal scars was found by comparing EDL933 and Sakai genomic sequences.
It is thus not surprising that different sets of phages in different insertion sites create a backdrop where, in the course of an outbreak, recombinational scrambling of phages, exacerbated by antibiotic treatments causing induction of prophages (Zhang et al., 2000; Kohler et al., 2000) or stress and SOS repair (Kimmitt et al., 2000), would further increase the variety of phages found (Schmidt, 2001; Herold et al., 2004; Bielaszewska et al., 2006). The occurrence of multiple inversion events demonstrated by the optical mapping in four strains indicates that one inversion might shuttle phage genes between a defective and infective prophage, thus altering the configuration of virulence genes typified by stx1 and stx2 genes. The initial recombination that translocated phage and virulence genes could be masked by the occurrence of a second inversion within the same pair of phages, or exposed by a second inversion involving a different phage pair. While changes in the toxin structural genes themselves can moderate toxin production, so too can changes in transcriptional regulation of genes (Zhang et al., 2005; Koitabashi et al., 2006), further emphasizing how recombinational events might potentiate or attenuate virulence.
Genomic comparisons of an assortment of microbes have not only revealed the extensive diversity that exists within a bacterial species, but they also are helping to challenge the very concept of a bacterial ‘species’ itself. For example, although there are only 75 000 single nucleotide polymorphisms (SNPs) within the homologues of E. coli O157 : H7 strain EDL933 and E. coli K-12 strain MG1665, nearly a million base pairs of unique DNA distinguish these two strains. Moreover, as individual strains of O157 : H7 show insertions and deletions of hundreds of thousands of base pairs, while exhibiting only a few hundred SNPs, methods that assess both rates of mutation and recombination must be used to assess the genesis, evolution and lineage of a particular strain. Optical mapping offers a relatively economical and rapid method for analysing insertions, deletions, inversions and rearrangements – those anomalies that ultimately shape the organization and circumscribe the gene order of an individual genome or strain.
The resolution and detail that optical maps provide underscore their usefulness for molecular epidemiological studies. The method should also prove useful for microbial forensics investigations, which require technologies that discriminate between closely related strains or that detect insertions in engineered strains (Budowle et al., 2005). In this respect, optical mapping offers a tractable solution for triaging microbial strains for complete genome sequencing.
|
ACKNOWLEDGEMENTS |
|---|
Edited by: J. Parkhill
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Bielaszewska, M., Prager, R., Zhang, W., Friedrich, A. W., Mellmann, A., Tschape, H. & Karch, H. (2006). Chromosomal dynamism in progeny of outbreak-related sorbitol-fermenting enterohemorrhagic Escherichia coli O157 : NM. Appl Environ Microbiol 72, 1900–1909.
Brussow, H., Canchaya, C. & Hardt, W. D. (2004). Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev 68, 560–602.
Budowle, B., Johnson, M. D., Fraser, C. M., Leighton, T. J., Murch, R. S. & Chakraborty, R. (2005). Genetic analysis and attribution of microbial forensics evidence. Crit Rev Microbiol 31, 233–254.[CrossRef][Medline]
Cai, W., Jing, J., Irvin, B., Ohler, L., Rose, E., Shizuya, H., Kim, U. J., Simon, M., Anantharaman, T. & other authors (1998). High-resolution restriction maps of bacterial artificial chromosomes constructed by optical mapping. Proc Natl Acad Sci U S A 95, 3390–3395.
Cebula, T. A., Brown, E. W., Jackson, S. A., Mammel, M. K., Mukherjee, A. & LeClerc, J. E. (2005). Molecular applications for identifying microbial pathogens in the post-9/11 era. Expert Rev Mol Diagn 5, 431–445.[CrossRef][Medline]
Chen, Q., Savarino, S. J. & Venkatesan, M. M. (2006). Subtractive hybridization and optical mapping of the enterotoxigenic Escherichia coli H10407 chromosome: isolation of unique sequences and demonstration of significant similarity to the chromosome of E. coli K-12. Microbiology 152, 1041–1054.
Feil, E. J. & Enright, M. C. (2004). Analyses of clonality and the evolution of bacterial pathogens. Curr Opin Microbiol 7, 308–313.[CrossRef][Medline]
Foley, S. L., Simjee, S., Meng, J., White, D. G., McDermott, P. F. & Zhao, S. (2004). Evaluation of molecular typing methods for Escherichia coli O157 : H7 isolates from cattle, food, and humans. J Food Prot 67, 651–657.[Medline]
Fukiya, S., Mizoguchi, H., Tobe, T. & Mori, H. (2004). Extensive genomic diversity in pathogenic Escherichia coli and Shigella strains revealed by comparative genomic hybridization microarray. J Bacteriol 186, 3911–3921.
Griffin, P. M. & Tauxe, R. V. (1991). The epidemiology of infections caused by Escherichia coli O157 : H7, other enterohemorrhagic E. coli, and the associated hemolytic uremic syndrome. Epidemiol Rev 13, 60–98.[Medline]
Griffin, P. M., Ostroff, S. M., Tauxe, R. V., Greene, K. D., Wells, J. G., Lewis, J. H. & Blake, P. A. (1988). Illnesses associated with Escherichia coli O157 : H7 infections. A broad clinical spectrum. Ann Intern Med 109, 705–712.
Herold, S., Karch, H. & Schmidt, H. (2004). Shiga toxin-encoding bacteriophages–genomes in motion. Int J Med Microbiol 294, 115–121.[CrossRef][Medline]
Iguchi, A., Iyoda, S., Terajima, J., Watanabe, H. & Osawa, R. (2006). Spontaneous recombination between homologous prophage regions causes large-scale inversions within the Escherichia coli O157 : H7 chromosome. Gene 372, 199–207.[CrossRef][Medline]
Jackson, S. A., Mammel, M. K., Patel, I. R., Mays, T., Albert, T. J., LeClerc, J. E. & Cebula, T. A. (2006). Interrogating genomic diversity of Escherichia coli O157 : H7 using DNA tiling arrays. Forensic Sci Int 2006, 23 (Epub ahead of print)
Kimmitt, P. T., Harwood, C. R. & Barer, M. R. (2000). Toxin gene expression by shiga toxin-producing Escherichia coli: the role of antibiotics and the bacterial SOS response. Emerg Infect Dis 6, 458–465.[Medline]
Kohler, B., Karch, H. & Schmidt, H. (2000). Antibacterials that are used as growth promoters in animal husbandry can affect the release of Shiga-toxin-2-converting bacteriophages and Shiga toxin 2 from Escherichia coli strains. Microbiology 146, 1085–1090.
Koitabashi, T., Vuddhakul, V., Radu, S., Morigaki, T., Asai, N., Nakaguchi, Y. & Nishibuchi, M. (2006). Genetic characterization of Escherichia coli O157 : H7 strains carrying the stx2 gene but not producing Shiga toxin 2. Microbiol Immunol 50, 135–148.[Medline]
Kotewicz, M. L., Brown, E. W., LeClerc, J. E. & Cebula, T. A. (2003). Genomic variability among enteric pathogens: the case of the mutS-rpoS intergenic region. Trends Microbiol 11, 2–6.[CrossRef][Medline]
Kudva, I. T., Evans, P. S., Perna, N. T., Barrett, T. J., DeCastro, G. J., Ausubel, F. M., Blattner, F. R. & Calderwood, S. B. (2002a). Polymorphic amplified typing sequences provide a novel approach to Escherichia coli O157 : H7 strain typing. J Clin Microbiol 40, 1152–1159.
Kudva, I. T., Evans, P. S., Perna, N. T., Barrett, T. J., Ausubel, F. M., Blattner, F. R. & Calderwood, S. B. (2002b). Strains of Escherichia coli O157 : H7 differ primarily by insertions or deletions, not single-nucleotide polymorphisms. J Bacteriol 184, 1873–1879.
Lim, A., Dimalanta, E. T., Potamousis, K. D., Yen, G., Apodoca, J., Tao, C., Lin, J., Qi, R., Skiadas, J. & other authors (2001). Shotgun optical maps of the whole Escherichia coli O157 : H7 genome. Genome Res 11, 1584–1593.
Liu, S. L. & Sanderson, K. E. (1995). The chromosome of Salmonella paratyphi A is inverted by recombination between rrnH and rrnG. J Bacteriol 177, 6585–6592.
Mead, P. S. & Griffin, P. M. (1998). Escherichia coli O157 : H7. Lancet 352, 1207–1212.[CrossRef][Medline]
Miesel, L., Segall, A. & Roth, J. R. (1994). Construction of chromosomal rearrangements in Salmonella by transduction: inversions of non-permissive segments are not lethal. Genetics 137, 919–932.[Abstract]
Nemoy, L. L., Kotetishvili, M., Tigno, J., Keefer-Norris, A., Harris, A. D., Perencevich, E. N., Johnson, J. A., Torpey, D., Sulakvelidze, A. & other authors (2005). Multilocus sequence typing versus pulsed-field gel electrophoresis for characterization of extended-spectrum beta-lactamase-producing Escherichia coli isolates. J Clin Microbiol 43, 1776–1781.
Ochman, H. (2005). Genomes on the shrink. Proc Natl Acad Sci U S A 102, 11959–11960.
Ogura, Y., Kurokawa, K., Ooka, T., Tashiro, K., Tobe, T., Ohnishi, M., Nakayama, K., Morimoto, T., Terajima, J. & other authors (2006). Complexity of the genomic diversity in enterohemorrhagic Escherichia coli O157 revealed by the combinational use of the O157 Sakai OligoDNA microarray and the Whole Genome PCR scanning. DNA Res 13, 3–14.
Ohnishi, M., Terajima, J., Kurokawa, K., Nakayama, K., Murata, T., Tamura, K., Ogura, Y., Watanabe, H. & Hayashi, T. (2002). Genomic diversity of enterohemorrhagic Escherichia coli O157 revealed by whole genome PCR scanning. Proc Natl Acad Sci U S A 99, 17043–17048.
Perna, N. T., Plunkett, G., III, Burland, V., Mau, B., Glasner, J. D., Rose, D. J., Mayhew, G. F., Evans, P. S., Gregor, J. & other authors (2001). Genome sequence of enterohaemorrhagic Escherichia coli O157 : H7. Nature 409, 529–533 (Erratum in Nature 2001 410, 240)[CrossRef][Medline]
Plunkett, G., III, Rose, D. J., Durfee, T. J. & Blattner, F. R. (1999). Sequence of Shiga toxin 2 phage 933W from Escherichia coli O157 : H7: shiga toxin as a phage late-gene product. J Bacteriol 181, 1767–1778.
Ribot, E. M., Fair, M. A., Gautom, R., Cameron, D. N., Hunter, S. B., Swaminathan, B. & Barrett, T. J. (2006). Standardization of pulsed-field gel electrophoresis protocols for the subtyping of Escherichia coli O157 : H7, Salmonella, and Shigella for PulseNet. Foodborne Pathog Dis 3, 59–67.[CrossRef][Medline]
Riley, L. W., Remis, R. S., Helgerson, S. D., McGee, H. B., Wells, J. G., Davis, B. R., Hebert, R. J., Olcott, E. S., Johnson, L. M. & other authors (1983). Hemorrhagic colitis associated with a rare Escherichia coli serotype. N Engl J Med 308, 681–685.[Abstract]
Schmidt, H. (2001). Shiga-toxin-converting bacteriophages. Res Microbiol 152, 687–695.[Medline]
Shaikh, N. & Tarr, P. I. (2003). Escherichia coli O157 : H7 Shiga toxin-encoding bacteriophages: integrations, excisions, truncations, and evolutionary implications. J Bacteriol 185, 3596–3605.
Shima, K., Wu, Y., Sugimoto, N., Asakura, M., Nishimura, K. & Yamasaki, S. (2006). Comparison of a PCR-restriction fragment length polymorphism (PCR-RFLP) assay to pulsed-field gel electrophoresis to determine the effect of repeated subculture and prolonged storage on RFLP patterns of Shiga toxin-producing Escherichia coli O157 : H7. J Clin Microbiol 44, 3963–3968.
Su, C. & Brandt, L. J. (1995). Escherichia coli O157 : H7 infection in humans. Ann Intern Med 123, 698–714.
Urwin, R. & Maiden, M. C. (2003). Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol 11, 479–487.[CrossRef][Medline]
Welch, R. A., Burland, V., Plunkett, G., III, Redford, P., Roesch, P., Rasko, D., Buckles, E. L., Liou, S. R., Boutin, A. & other authors (2002). Extensive mosaic structure revealed by the complete genome sequence of uropathogenic Escherichia coli. Proc Natl Acad Sci U S A 99, 17020–17024.
Whittam, T. S., Wolfe, M. L., Wachsmuth, I. K., Orskov, F., Orskov, I. & Wilson, R. A. (1993). Clonal relationships among Escherichia coli strains that cause hemorrhagic colitis and infantile diarrhea. Infect Immun 61, 1619–1629.
Wick, L. M., Qi, W., Lacher, D. W. & Whittam, T. S. (2005). Evolution of genomic content in the stepwise emergence of Escherichia coli O157 : H7. J Bacteriol 187, 1783–1791.
Zhang, W., Bielaszewska, M., Friedrich, A. W., Kuczius, T. & Karch, H. (2005). Transcriptional analysis of genes encoding Shiga toxin 2 and its variants in Escherichia coli. Appl Environ Microbiol 71, 558–561.
Zhang, W., Qi, W., Albert, T. J., Motiwala, A. S., Alland, D., Hyytia-Trees, E. K., Ribot, E. M., Fields, P. I., Whittam, T. S. & Swaminathan, B. (2006). Probing genomic diversity and evolution of Escherichia coli O157 by single-nucleotide polymorphisms. Genome Res 16, 757–767.
Zhang, X., McDaniel, A. D., Wolf, L. E., Keusch, G. T., Waldor, M. K. & Acheson, D. W. (2000). Quinolone antibiotics induce Shiga toxin-encoding bacteriophages, toxin production, and death in mice. J Infect Dis 181, 664–670.[CrossRef][Medline]
Zhou, S., Kile, A., Bechner, M., Place, M., Kvikstad, E., Deng, W., Wei, J., Severin, J., Runnheim, R. & other authors (2004). Single-molecule approach to bacterial genomic comparisons via optical mapping. J Bacteriol 186, 7773–7782.
Received 20 November 2006;
revised 8 February 2007;
accepted 12 February 2007.
This article has been cited by other articles:
|
|
 |
|
  Z. Yang, J. Kim, C. Zhang, M. Zhang, J. Nietfeldt, C. M. Southward, M. G. Surette, S. D. Kachman, and A. K. Benson 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 J. Bacteriol., June 1, 2009; 191(11): 3553 - 3568. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. L. Kotewicz, M. K. Mammel, J. E. LeClerc, and T. A. Cebula Optical mapping and 454 sequencing of Escherichia coli O157 : H7 isolates linked to the US 2006 spinach-associated outbreak Microbiology, November 1, 2008; 154(11): 3518 - 3528. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. A. Rasko, M. J. Rosovitz, G. S. A. Myers, E. F. Mongodin, W. F. Fricke, P. Gajer, J. Crabtree, M. Sebaihia, N. R. Thomson, R. Chaudhuri, et al. The Pangenome Structure of Escherichia coli: Comparative Genomic Analysis of E. coli Commensal and Pathogenic Isolates J. Bacteriol., October 15, 2008; 190(20): 6881 - 6893. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |
