Promiscuous origin of a chimeric sequence in the Escherichia coli O157:H7 genome

The rpoS gene is predominantly inactivated during laboratory storage and undergoes source-sink evolution in Escherichia coli species

The rpoS gene codes for an alternative RNA polymerase sigma factor, which acts as a general regulator of the stress response. Inactivating alleles of rpoS in collections of natural Escherichia coli isolates have been observed at very variable frequencies, from less than 1% to more than 70% of strains. rpoS is easily inactivated in nutrient-deprived environments such as stab storage, which makes it difficult to determine the true frequency of rpoS inactivation in nature. We studied the evolutionary history of rpoS and compared it to the phylogenetic history of bacteria in two collections of 82 human commensal and extraintestinal E. coli strains. These strains were representative of the phylogenetic diversity of the species and differed only by their storage conditions. In both collections, the phylogenetic histories of rpoS and of the strains were congruent, indicating that horizontal gene transfer had not occurred at the rpoS locus, and rpoS was under strong purifying selection, with a ratio of the nonsynonymous mutation rate (Ka) to the synonymous substitution rate (Ks) substantially smaller than 1. Stab storage was associated with a high frequency of inactivating alleles, whereas almost no amino acid sequence variation was observed in RpoS in the collection studied directly after isolation of the strains from the host. Furthermore, the accumulation of variations in rpoS was typical of source-sink dynamics. In conclusion, rpoS is rarely inactivated in natural E. coli isolates within their mammalian hosts, probably because such strains rapidly become evolutionary dead ends. Our data should encourage bacteriologists to freeze isolates immediately and to avoid the use of stab storage.

Correlation between the genomic o454-nlpD region polymorphisms, virulence gene equipment and phylogenetic group of extraintestinal Escherichia coli (ExPEC) enables pathotyping irrespective of host, disease and source of isolation

BACKGROUND

The mutS-rpoS intergenic region in E. coli displays a mosaic structure which revealed pathotype specific patterns. To assess the importance of this region as a surrogate marker for the identification of highly virulent extraintestinal pathogenic E. coli (ExPEC) strains we aimed to: (i) characterize the genetic diversity of the mutS gene and the o454-nlpD genomic region among 510 E. coli strains from animals and humans; (ii) delineate associations between the polymorphism of this region and features such as phylogenetic background of E. coli, pathotype, host species, clinical condition, serogroup and virulence associated genes (VAG)s; and (iii) identify the most important VAGs for classification of the o454-nlpD region.

METHODS

Size variation in the o454-nlpD region was investigated by PCR amplification and sequencing. Phylogenetic relationships were assessed by Ecor- and Multilocus sequence- typing (MLST), and a comparative analysis between mutS gene phylogenetic tree obtained with RAxML and the MLST grouping method was performed. Correlation between o454-nlpD patterns and the features described above were analysed. In addition, the importance of 47 PCR-amplified ExPEC-related VAGs for classification of o454-nlpD patterns was investigated by means of Random Forest algorithm.

RESULTS

Four main structures (patterns I-IV) of the o454-nlpD region among ExPEC and commensal E. coli strains were identified. Statistical analysis showed a positive and exclusive association between pattern III and the ExPEC strains. A strong association between pattern III and either the Ecor group B2 or the sequence type complexes known to represent the phylogenetic background of highly virulent ExPEC strains (such as STC95, STC73 and STC131) was found as well. RF analyses determined five genes (csgA, malX, chuA, sit, and vat) to be suitable to predict pattern III strains.

CONCLUSION

The significant association between pattern III and group B2 strains suggested the o454-nlpD region to be of great value in identifying highly virulent strains among the mixed population of E. coli promising to be the basis of a future typing tool for ExPEC and their gut reservoir. Furthermore, top-ranked VAGs for classification and prediction of pattern III were identified. These data are most valuable for defining ExPEC pathotype in future in vivo assays.

Molecular applications for identifying microbial pathogens in the post-9/11 era

Rapid advances in molecular and optical technologies over the past 10 years have dramatically impacted the way biologic research is conducted today. Examples include microarrays, capillary sequencing, optical mapping and real-time sequencing (Pyrosequencing). These technologies are capable of rapidly delivering massive amounts of genetic information and are becoming routine mainstays of many laboratories. Fortunately, advances in scientific computing have provided the enormous computing power necessary to analyze these enormous data sets. The application of molecular technologies should prove useful to the burgeoning field of microbial forensics. In the post-9/11 era, when securing America's food supply is a major endeavor, the need for rapid identification of microbes that accidentally or intentionally find their way into foods is apparent. The principle that distinguishes a microbial forensic investigation from a molecular epidemiology study is that a biocrime has been committed. If proper attribution is to be attained, a link must be made between a particular microbe in the food and the perpetrator who placed it there. Therefore, the techniques used must be able to discriminate individual isolates of a particular microbe. A battery of techniques in development for distinguishing individual isolates of particular foodborne pathogens is discussed.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background

Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor σ^S (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain.

Results

Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms.

Conclusions

The ppGpp and σ^S variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Structure and distribution of the phosphoprotein phosphatase genes, prpA and prpB, among Shigella subgroups

Phosphoprotein phosphatases encoded by theprpAandprpBgenes function in signal transduction pathways for degradation of misfolded proteins in the extracytoplasmic compartments ofEscherichia coli. In order to trace the evolution ofprpgenes and assess their roles in other enteric pathogens, the structure and distribution of these genes among closely related Shigella subgroups were studied. PCR amplification, probe hybridization studies and DNA sequencing were used to determine theprpgenotypes of 58 strains from the four Shigella subgroups, Dysenteriae, Boydii, Sonnei and Flexneri. It was found that theprpalleles among Shigella subgroups were extremely susceptible to gene inactivation and that the mutations involved inprpallele inactivation were varied. They included IS insertions, gene replacement by an IS element, a small deletion within the gene or large deletion engulfing the entire gene region, and base substitutions that generated premature termination codons. As a result, of 58 strains studied, only eight (14 %) possessed intactprpAandprpBgenes. Of the Shigella strains examined, 76 % (44/58) showed at least one of theprpalleles inactivated by one or more IS elements, including IS1, IS4, IS600and IS629. Phylogenetic analysis revealed that IS elements have been independently acquired in multiple lineages of Shigella, suggesting that loss of functional alleles has been advantageous during Shigella strain evolution.

Insertions, deletions, and single-nucleotide polymorphisms at rare restriction enzyme sites enhance discriminatory power of polymorphic amplified typing sequences, a novel strain typing system for Escherichia coli O157:H7

Polymorphic amplified typing sequences (PATS) for Escherichia coli O157:H7 (O157) was previously based on indels containing XbaI restriction enzyme sites occurring in O-island sequences of the O157 genome. This strain-typing system, referred to as XbaI-based PATS, typed every O157 isolate tested in a reproducible, rapid, straightforward, and easy-to-interpret manner and had technical advantages over pulsed-field gel electrophoresis (PFGE). However, the system was less discriminatory than PFGE and was unable to differentiate fully between unrelated isolates. To overcome this drawback, we enhanced PATS by using another infrequently cutting restriction enzyme, AvrII (also known as BlnI), to identify additional polymorphic regions that could increase the discriminatory ability of PATS typing. Referred to as AvrII-based PATS, the system identified seven new polymorphic regions in the O157 genome. Unlike XbaI, polymorphisms involving AvrII sites were caused by both indels and single-nucleotide polymorphisms occurring in O-island and backbone sequences of the O157 genome. AvrII-based PATS by itself provided poor discrimination of the O157 isolates tested. However, when primer pairs amplifying the seven polymorphic AvrII sites were combined with those amplifying the eight polymorphic XbaI sites (combined PATS), the discriminatory power of PATS was enhanced. Combined PATS matched related O157 isolates better than PFGE while differentiating between unrelated isolates. PATS typed every O157 isolate tested and directly targeted polymorphic sequences responsible for differences in the restriction digest patterns of O157 genomic DNA, utilizing PCR rather than relying on gel electrophoresis. This enabled PATS to resolve the ambiguity in PFGE typing, including that arising from the "more distantly related" and "untypeable" profiles.

Limited boundaries for extensive horizontal gene transfer among Salmonella pathogens

Recombination is thought to be rare within Salmonella, as evidenced by absence of gene transfer among SARC strains that represent the broad genetic diversity of the eight primary subspecies of this common facultative intracellular pathogen. We adopted a phylogenetic approach to assess recombination within the mutS gene of 70 SARB strains, a genetically homogeneous population of Salmonella enterica subspecies I strains, which have in common the ability to infect warm-blooded animals. We report here that SARB strains show evidence for widespread recombinational exchange in contrast to results obtained with strains exhibiting species-level genetic variation. Besides extensive allele shuffling, SARB strains showed notably larger recombinagenic patch sizes for mutS (at least approximately 1.1 kb) than previously reported for S. enterica SARC strains. Explaining these experimental dichotomies provides important insight for understanding microbial evolution, because they suggest likely ecologic and genetic barriers that limit extensive gene transfer in the feral setting.

Molecular analysis of mutS expression and mutation in natural isolates of pathogenic Escherichia coli

Deficiencies in the MutS protein disrupt methyl-directed mismatch repair (MMR), generating a mutator phenotype typified by high mutation rates and promiscuous recombination. How such deficiencies might arise in the natural environment was determined by analysing pathogenic strains of Escherichia coli. Quantitative Western immunoblotting showed that the amount of MutS in a wild-type strain of the enterohaemorrhagic pathogen E. coli O157 : H7 decreased about 26-fold in stationary-phase cells as compared with the amount present during exponential-phase growth. The depletion of MutS in O157 : H7 is significantly greater than that observed for a laboratory-attenuated E. coli K-12 strain. In the case of stable mutators, mutS defects in strains identified among natural isolates were analysed, including two E. coli O157 : H7 strains, a diarrhoeagenic E. coli O55 : H7 strain, and a uropathogenic strain from the E. coli reference (ECOR) collection. No MutS could be detected in the four strains by Western immunoblot analyses. RNase T2 protection assays showed that the strains were either deficient in mutS transcripts or produced transcripts truncated at the 3' end. Nucleotide sequence analysis revealed extensive deletions in the mutS region of three strains, ranging from 7.5 to 17.3 kb relative to E. coli K-12 sequence, while the ECOR mutator contained a premature stop codon in addition to other nucleotide changes in the mutS coding sequence. These results provide insights into the status of the mutS gene and its product in pathogenic strains of E. coli.

Genomic variability among enteric pathogens: the case of the mutS-rpoS intergenic region

The mutS-rpoS intergenic region of enteric bacteria ranges in size from 88 bp in Yersinia to > 12000 bp in Salmonella. We interpret this expansion as the result of the horizontal transfer of segments of DNA from diverse origins. Both comparative genomic analysis and selective sequencing of a variety of Escherichia coli pathogens have provided additional evidence for reassortment of segments within this region.

Extended virulence genotype of pathogenic Escherichia coli isolates carrying the afa-8 operon: evidence of similarities between isolates from humans and animals with extraintestinal infections

The afimbrial AfaE-VIII adhesin is common among Escherichia coli isolates from calves with intestinal and/or extraintestinal infections and from humans with sepsis or pyelonephritis. The virulence genotypes of 77 Escherichia coli afa-8 isolates from farm animals and humans were compared to determine whether any trait of commonality exists between isolates of the different host species. Over half of the extraintestinal afa-8 isolates were associated with pap and f17Ac adhesin genes and contained virulence genes (pap, hly, and cnf1) which are characteristic of human extraintestinal pathogenic E. coli (ExPEC). PapG, which occurs as three known variants (variants I to III), is encoded by the corresponding three alleles of papG. Among the pap-positive strains, new papG variants (papGrs) that differed from the isolates with genes for the three adhesin classes predominated over isolates with papG allele III, which in turn were more prevalent than those with allele II. The data showed the substantial prevalence of the enteroaggregative E. coli heat-stable enterotoxin gene (east1) among afa-8 isolates. Most of the afa-8 isolates harbored the high-pathogenicity island (HPI) present in pathogenic Yersinia; however, two-thirds of the HPI-positive strains shared a truncated HPI integrase gene. The presence of ExPEC-associated virulence factors (VFs) in extraintestinal isolates that carry genes typical of enteric strains and that express O antigens associated with intestinal E. coli is consistent with transfer of VFs and O-antigen determinants between ExPEC and enteric strains. The similarities between animal and human ExPEC strains support the hypothesis of overlapping populations, with members of certain clones or clonal groups including animal and human strains. The presence of multiple-antibiotic-resistant bovine afa-8 strains among such clones may represent a potential public health risk.

Evolution of multi-gene segments in the mutS-rpoS intergenic region of Salmonella enterica serovar Typhimurium LT2

The nucleotide sequence of the 12.6 kb region between the mutS and rpoS genes of Salmonella enterica serovar Typhimurium LT2 (S. typhimurium) was compared to other enteric bacterial mutS-rpoS intergenic regions. The mutS-rpoS region is composed of three distinct segments, designated HK, O and S, as defined by sequence similarities to contiguous ORFs in other bacteria. Inverted chromosomal orientations of each of these segments are found between the mutS and rpoS genes in related ENTEROBACTERIACEAE: The HK segment is distantly related to a cluster of seven ORFs found in Haemophilus influenzae and a cluster of five ORFs found between the mutS and rpoS genes in Escherichia coli K-12. The O segment is related to the mutS-rpoS intergenic region found in E. coli O157:H7 and Shigella dysenteriae type 1. The third segment, S, is common to diverse Salmonella species, but is absent from E. coli. Despite the extensive collinearity and conservation of the overall genetic maps of S. typhimurium and E. coli K-12, the insertions, deletions and inversions in the mutS-rpoS region provide evidence that this region of the chromosome is an active site for horizontal gene transfer and rearrangement.

Detection of recombination among Salmonella enterica strains using the incongruence length difference test

Particular serovars of Salmonella enterica have emerged as significant foodborne pathogens in humans. At the chromosomal level, discrete regions in the Salmonella genome have been identified that are known to play important roles in the maintenance, survival, and virulence of S. enterica within the host. Interestingly, several of these loci appear to have been acquired by horizontal transfer of DNA among and between bacterial species. The profound importance of recombination in pathogen emergence is just now being realized, perhaps explaining the sudden interest in developing novel and facile ways for detecting putative horizontal transfer events in bacteria. The incongruence length difference (ILD) test offers one such means. ILD uses phylogeny to trace sequences that may have been acquired promiscuously by exchange of DNA during chromosome evolution. We show here that the ILD test readily detects recombinations that have taken place in several housekeeping genes in Salmonella as well as genes composing the type 1 pilin complex (14 min) and the inv-spa invasion gene complex (63 min). Moreover, the ILD test indicated that the mutS gene (64 min), whose product helps protect the bacterial genome from invasion by foreign DNA, appears to have undergone intragenic recombination within S. enterica subspecies I. ILD findings were supported using additional tests known to be independent of the ILD approach (e.g., split decomposition analysis and compatibility of sites). Taken together, these data affirm the application of the ILD test as one approach for identifying recombined sequences in the Salmonella chromosome. Furthermore, horizontally acquired sequences within mutS support a model whereby evolutionarily important recombinants of S. enterica are rescued from strains carrying defective mutS alleles via horizontal transfer.

Genomic subtraction to identify and characterize sequences of Shiga toxin-producing Escherichia coli O91:H21

To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.

Regulatory RNA as mediator in GacA/RsmA-dependent global control of exoproduct formation in Pseudomonas fluorescens CHA0

In Pseudomonas fluorescens CHA0, an antagonist of root-pathogenic fungi, the GacS/GacA two-component system tightly controls the expression of antifungal secondary metabolites and exoenzymes at a posttranscriptional level, involving the RNA-binding protein and global regulator of secondary metabolism RsmA. This protein was purified from P. fluorescens, and RNA bound to it was converted to cDNA, which served as a probe to isolate the corresponding chromosomal locus, rsmZ. This gene encoded a regulatory RNA of 127 nucleotides and a truncated form lacking 35 nucleotides at the 3' end. Expression of rsmZ depended on GacA, increased with increasing population density, and was stimulated by the addition of a solvent-extractable extracellular signal produced by strain CHA0 at the end of exponential growth. This signal appeared to be unrelated to N-acyl-homoserine lactones. A conserved upstream element in the rsmZ promoter, but not the stress sigma factor RpoS, was involved in rsmZ expression. Overexpression of rsmZ effectively suppressed the negative effect of gacS and gacA mutations on target genes, i.e., hcnA (for hydrogen cyanide synthase) and aprA (for the major exoprotease). Mutational inactivation of rsmZ resulted in reduced expression of these target genes in the presence of added signal. Overexpression of rsmA had a similar, albeit stronger negative effect. These results support a model in which GacA upregulates the expression of regulatory RNAs, such as RsmZ of strain CHA0, in response to a bacterial signal. By a titration effect, RsmZ may then alleviate the repressing activity of RsmA on the expression of target mRNAs.

Phylogenetic evidence for horizontal transfer of mutS alleles among naturally occurring Escherichia coli strains

mutS mutators accelerate the bacterial mutation rate 100- to 1,000-fold and relax the barriers that normally restrict homeologous recombination. These mutators thus afford the opportunity for horizontal exchange of DNA between disparate strains. While much is known regarding the mutS phenotype, the evolutionary structure of the mutS(+) gene in Escherichia coli remains unclear. The physical proximity of mutS to an adjacent polymorphic region of the chromosome suggests that this gene itself may be subject to horizontal transfer and recombination events. To test this notion, a phylogenetic approach was employed that compared gene phylogeny to strain phylogeny, making it possible to identify E. coli strains in which mutS alleles have recombined. Comparison of mutS phylogeny against predicted E. coli "whole-chromosome" phylogenies (derived from multilocus enzyme electrophoresis and mdh sequences) revealed striking levels of phylogenetic discordance among mutS alleles and their respective strains. We interpret these incongruences as signatures of horizontal exchange among mutS alleles. Examination of additional sites surrounding mutS also revealed incongruous distributions compared to E. coli strain phylogeny. This suggests that other regional sequences are equally subject to horizontal transfer, supporting the hypothesis that the 61.5-min mutS-rpoS region is a recombinational hot spot within the E. coli chromosome. Furthermore, these data are consistent with a mechanism for stabilizing adaptive changes promoted by mutS mutators through rescue of defective mutS alleles with wild-type sequences.

The Escherichia coli gene pool

Wild Escherichia coli are superbly adapted to survive in the intestines of their mammalian hosts and in the environment. E. coli K12 derivative (MG1655) encodes 4288 potential genes that provide the background genetic framework of this species. Particular E. coli clonal types encode additional chromosomal and extrachromosomal genes that facilitate the ability of E. coli to adapt to new environments. These additional genes are often clustered, have related functions (for example, virulence-associated genes in pathogenicity islands) and may be integrated at specific sites on the E. coli chromosome.

Three R's of bacterial evolution: how replication, repair, and recombination frame the origin of species

The genetic diversity of bacteria results not only from errors in DNA replication and repair but from horizontal exchange and recombination of DNA sequences from similar and disparate species as well. New individuals carrying adaptive changes are thus being spawned constantly among the population at large. When new selection pressures appear, these are the individuals that survive, at the expense of the general population, to forge new populations. Depending on the severity and uniqueness of the selection pressure, this could lead to new speciation. It is becoming more and more evident that, as nucleotide sequences of numerous loci from many bacterial strains continue to amass, horizontal transfer has played a key role in configuring the Escherichia coli chromosome. Here, we examine views, both old and new, for the role of recombination in the evolution of bacterial chromosomes. We present novel phylogenetic evidence for horizontal transfer of three genes involved in DNA replication and repair (mutS, uvrD, and polA). These data reveal a prominent role for horizontal transfer in the evolution of genes known to play a key role in the fidelity of DNA replication and, thus, ultimate survival of the organism. Our data underscore that recombination plays both a diversifying and a homogenizing role in defining the structure of the E. coli genome.

An Escherichia coli reference collection group B2- and uropathogen-associated polymorphism in the rpoS-mutS region of the E. coli chromosome

Chromosomal DNAs of enterohemorrhagic, uropathogenic, and laboratory attenuated Escherichia coli strains differ in the rpoS-mutS region. Many uropathogens lack a deletion and an insertion characteristic of enterohemorrhagic strains. At the same chromosomal position, they harbor a 2.1-kb insertion of unknown origin with a base composition suggestive of horizontal gene transfer. Unlike virulence determinants associated with urinary tract infection and/or neonatal meningitis (pap or prs, sfa, kps, and hly), the 2.1-kb insertion is shared by all group B2 strains of the E. coli Reference Collection.

Gene conservation and loss in the mutS-rpoS genomic region of pathogenic Escherichia coli

The extent and nature of DNA polymorphism in the mutS-rpoS region of the Escherichia coli genome were assessed in 21 strains of enteropathogenic E. coli (EPEC) and enterohemorrhagic E. coli (EHEC) and in 6 strains originally isolated from natural populations. The intervening region between mutS and rpoS was amplified by long-range PCR, and the resulting amplicons varied substantially in length (7.8 to 14.2 kb) among pathogenic groups. Restriction maps based on five enzymes and sequence analysis showed that strains of the EPEC 1, EPEC 2, and EHEC 2 groups have a long mutS-rpoS region composed of a approximately 6.0-kb DNA segment found in strain K-12 and a novel DNA segment ( approximately 2.9 kb) located at the 3' end of rpoS. The novel segment contains three genes (yclC, pad1, and slyA) that occur in E. coli O157:H7 and related strains but are not found in K-12 or members of the ECOR group A. Phylogenetic analysis of the common sequences indicates that the long intergenic region is ancestral and at least two separate deletion events gave rise to the shorter regions characteristic of the E. coli O157:H7 and K-12 lineages.