Molecular evolution of the Escherichia coli chromosome. VI. Two regions of high effective recombination

A convergent evolutionary pathway attenuating cellulose production drives enhanced virulence of some bacteria

Bacteria adapt to selective pressure in their immediate environment in multiple ways. One mechanism involves the acquisition of independent mutations that disable or modify a key pathway, providing a signature of adaptation via convergent evolution. Extra-intestinal pathogenic Escherichia coli (ExPEC) belonging to sequence type 95 (ST95) represent a global clone frequently associated with severe human infections including acute pyelonephritis, sepsis, and neonatal meningitis. Here, we analysed a publicly available dataset of 613 ST95 genomes and identified a series of loss-of-function mutations that disrupt cellulose production or its modification in 55.3% of strains. We show the inability to produce cellulose significantly enhances ST95 invasive infection in a rat model of neonatal meningitis, leading to the disruption of intestinal barrier integrity in newborn pups and enhanced dissemination to the liver, spleen and brain. Consistent with these observations, disruption of cellulose production in ST95 augmented innate immune signalling and tissue neutrophil infiltration in a mouse model of urinary tract infection. Mutations that disrupt cellulose production were also identified in other virulent ExPEC STs, Shigella and Salmonella, suggesting a correlative association with many Enterobacteriaceae that cause severe human infection. Together, our findings provide an explanation for the emergence of hypervirulent Enterobacteriaceae clones.

High-resolution genomic analysis to investigate the impact of the invasive brushtail possum (Trichosurus vulpecula) and other wildlife on microbial water quality assessments

Escherichia coli are routine indicators of fecal contamination in water quality assessments. Contrary to livestock and human activities, brushtail possums (Trichosurus vulpecula), common invasive marsupials in Aotearoa/New Zealand, have not been thoroughly studied as a source of fecal contamination in freshwater. To investigate their potential role, Escherichia spp. isolates (n = 420) were recovered from possum gut contents and feces and were compared to those from water, soil, sediment, and periphyton samples, and from birds and other introduced mammals collected within the Mākirikiri Reserve, Dannevirke. Isolates were characterized using E. coli-specific real-time PCR targeting the uidA gene, Sanger sequencing of a partial gnd PCR product to generate a gnd sequence type (gST), and for 101 isolates, whole genome sequencing. Escherichia populations from 106 animal and environmental sample enrichments were analyzed using gnd metabarcoding. The alpha diversity of Escherichia gSTs was significantly lower in possums and animals compared with aquatic environmental samples, and some gSTs were shared between sample types, e.g., gST535 (in 85% of samples) and gST258 (71%). Forty percent of isolates gnd-typed and 75% of reads obtained by metabarcoding had gSTs shared between possums, other animals, and the environment. Core-genome single nucleotide polymorphism (SNP) analysis showed limited variation between several animal and environmental isolates (<10 SNPs). Our data show at an unprecedented scale that Escherichia clones are shared between possums, other wildlife, water, and the wider environment. These findings support the potential role of possums as contributors to fecal contamination in Aotearoa/New Zealand freshwater. Our study deepens the current knowledge of Escherichia populations in under-sampled wildlife. It presents a successful application of high-resolution genomic methods for fecal source tracking, thereby broadening the analytical toolbox available to water quality managers. Phylogenetic analysis of isolates and profiling of Escherichia populations provided useful information on the source(s) of fecal contamination and suggest that comprehensive invasive species management strategies may assist in restoring not only ecosystem health but also water health where microbial water quality is compromised.

O-Antigen Diversification Masks Identification of Highly Pathogenic Shiga Toxin-Producing Escherichia coli O104:H4-Like Strains

Shiga toxin-producing Escherichia coli (STEC) can give rise to a range of clinical outcomes from diarrhea to the life-threatening systemic condition hemolytic-uremic syndrome (HUS). Although STEC O157:H7 is the serotype most frequently associated with HUS, a major outbreak of HUS occurred in 2011 in Germany and was caused by a rare serotype, STEC O104:H4. Prior to 2011 and since the outbreak, STEC O104:H4 strains have only rarely been associated with human infections. From 2012 to 2020, intensified STEC surveillance was performed in Germany where the subtyping of ~8,000 clinical isolates by molecular methods, including whole-genome sequencing, was carried out. A rare STEC serotype, O181:H4, associated with HUS was identified, and like the STEC O104:H4 outbreak strain, this strain belongs to sequence type 678 (ST678). Genomic and virulence comparisons revealed that the two strains are phylogenetically related and differ principally in the gene cluster encoding their respective lipopolysaccharide O-antigens but exhibit similar virulence phenotypes. In addition, five other serotypes belonging to ST678 from human clinical infection, such as OX13:H4, O127:H4, OgN-RKI9:H4, O131:H4, and O69:H4, were identified from diverse locations worldwide. IMPORTANCE Our data suggest that the high-virulence ensemble of the STEC O104:H4 outbreak strain remains a global threat because genomically similar strains cause disease worldwide but that the horizontal acquisition of O-antigen gene clusters has diversified the O-antigens of strains belonging to ST678. Thus, the identification of these highly pathogenic strains is masked by diverse and rare O-antigens, thereby confounding the interpretation of their potential risk.

The impact of genetic recombination on pathogenic Leptospira

Leptospirosis is the most common zoonosis worldwide, and is increasingly common in poor urban communities, where there is inadequate sewage disposal and abundance of domestic and peridomestic animals. There are many risk factors associated with the disease, such as contaminated water exposure, close contact with animals, floods, recreational activities related to water, wet agriculture. The symptoms of leptospirosis are common to other infectious diseases and, if not treated, it can lead to meningitis, liver failure, kidney damage and death. Leptospirosis is caused by 38 pathogenic species of Leptospira, which are divided in almost 30 serogroups and more than 300 serovars. The serological classification (serogroups and serovars) is based on the expression of distinct lipopolysaccharide (LPS) antigens. These antigens are also associated to protective immunity; antibodies against a serovar protect from any member of the same serovar. Serologic and phylogenetic analyses are not congruent probably due to genetic recombination of LPS genes among different leptospiral species. To analyze the importance of recombination in leptospiral evolution, we performed a gene-by-gene tree topology comparison on closed genomes available in public databases at two levels: among core genomes of pathogenic species (34 recombinant among 1213 core genes), and among core genomes of L. interrogans isolates (178/798). We found that most recombinant genes code for proteins involved in translation, ribosomal structure and biogenesis, but also for cell wall, membrane and envelope biogenesis. Besides, our final results showed that half of LPS genes are recombinant (18/36). This is relevant because serovar classification and vaccine development are based on these epitopes. The frequent recombination of LPS-associated genes suggests that natural selection is promoting the survival of recombinant lineages. These results may help understanding the factors that make Leptospira a successful pathogen.

A role for ColV plasmids in the evolution of pathogenic Escherichia coli ST58

Escherichia coli ST58 has recently emerged as a globally disseminated uropathogen that often progresses to sepsis. Unlike most pandemic extra-intestinal pathogenic E. coli (ExPEC), which belong to pathogenic phylogroup B2, ST58 belongs to the environmental/commensal phylogroup B1. Here, we present a pan-genomic analysis of a global collection of 752 ST58 isolates from diverse sources. We identify a large ST58 sub-lineage characterized by near ubiquitous carriage of ColV plasmids, which carry genes encoding virulence factors, and by a distinct accessory genome including genes typical of the Yersiniabactin High Pathogenicity Island. This sub-lineage includes three-quarters of all ExPEC sequences in our study and has a broad host range, although poultry and porcine sources predominate. By contrast, strains isolated from cattle often lack ColV plasmids. Our data indicate that ColV plasmid acquisition contributed to the divergence of the major ST58 sub-lineage, and different sub-lineages inhabit poultry, swine and cattle.

Exopolysaccharides Producing Lactic Acid Bacteria in Wine and Other Fermented Beverages: For Better or for Worse?

Lactic acid bacteria (LAB) from fermented beverages such as wine, cider and beer produce a wide range of exopolysaccharides (EPS) through multiple biosynthetic pathways. These extracellular polysaccharides constitute key elements for bacterial species adaptation to such anthropic processes. In the food industry, LAB polysaccharides have been widely studied for their rheological, functional and nutritional properties; however, these have been poorly studied in wine, beer and cider until recently. In this review, we have gathered the information available on these specific polysaccharide structure and, biosynthetic pathways, as well as the physiology of their production. The genes associated with EPS synthesis are also presented and compared. Finally, the possible role of EPS for bacterial survival and spread, as well as the risks or possible benefits for the winemaker and the wine lover, are discussed.

The role of interspecies recombination in the evolution of antibiotic-resistant pneumococci

Multidrug-resistant Streptococcus pneumoniae emerge through the modification of core genome loci by interspecies homologous recombinations, and acquisition of gene cassettes. Both occurred in the otherwise contrasting histories of the antibiotic-resistant S. pneumoniae lineages PMEN3 and PMEN9. A single PMEN3 clade spread globally, evading vaccine-induced immunity through frequent serotype switching, whereas locally circulating PMEN9 clades independently gained resistance. Both lineages repeatedly integrated Tn916-type and Tn1207.1-type elements, conferring tetracycline and macrolide resistance, respectively, through homologous recombination importing sequences originating in other species. A species-wide dataset found over 100 instances of such interspecific acquisitions of resistance cassettes and flanking homologous arms. Phylodynamic analysis of the most commonly sampled Tn1207.1-type insertion in PMEN9, originating from a commensal and disrupting a competence gene, suggested its expansion across Germany was driven by a high ratio of macrolide-to-β-lactam consumption. Hence, selection from antibiotic consumption was sufficient for these atypically large recombinations to overcome species boundaries across the pneumococcal chromosome.

The population genetics of pathogenic Escherichia coli

Escherichia coli is a commensal of the vertebrate gut that is increasingly involved in various intestinal and extra-intestinal infections as an opportunistic pathogen. Numerous pathotypes that represent groups of strains with specific pathogenic characteristics have been described based on heterogeneous and complex criteria. The democratization of whole-genome sequencing has led to an accumulation of genomic data that render possible a population phylogenomic approach to the emergence of virulence. Few lineages are responsible for the pathologies compared with the diversity of commensal strains. These lineages emerged multiple times during E. coli evolution, mainly by acquiring virulence genes located on mobile elements, but in a specific chromosomal phylogenetic background. This repeated emergence of stable and cosmopolitan lineages argues for an optimization of strain fitness through epistatic interactions between the virulence determinants and the remaining genome.

gndDb, a Database of Partial gnd Sequences To Assist with Analysis of Escherichia coli Communities Using High-Throughput Sequencing

The use of culture methods to detect Escherichia coli diversity does not provide sufficient resolution to identify strains present at low levels. Here, we target the hypervariable gnd gene and describe a database containing 534 distinct partial gnd sequences and associated O groups for use with culture-independent E. coli community analysis.

The use of culture methods to detect Escherichia coli diversity does not provide sufficient resolution to identify strains present at low levels. Here, we target the hypervariable gnd gene and describe a database containing 534 distinct partial gnd sequences and associated O groups for use with culture-independent E. coli community analysis.

Diversity-Generating Machines: Genetics of Bacterial Sugar-Coating

Bacterial pathogens and commensals are surrounded by diverse surface polysaccharides which include capsules and lipopolysaccharides. These carbohydrates play a vital role in bacterial ecology and interactions with the environment. Here, we review recent rapid advancements in this field, which have improved our understanding of the roles, structures, and genetics of bacterial polysaccharide antigens. Genetic loci encoding the biosynthesis of these antigens may have evolved as bacterial diversity-generating machines, driven by selection from a variety of forces, including host immunity, bacteriophages, and cell-cell interactions. We argue that the high adaptive potential of polysaccharide antigens should be taken into account in the design of polysaccharide-targeting medical interventions like conjugate vaccines and phage-based therapies.

Culture independent analysis using gnd as a target gene to assess Escherichia coli diversity and community structure

Current culture methods to investigate changes in Escherichia coli community structure are often slow and laborious. Genes such as gnd (6-phosphogluconate dehydrogenase) have a highly variable nucleotide sequence and may provide a target for E. coli microbiome analysis using culture-independent methods. Metabarcoded PCR primers were used to generate separate libraries from calf faecal samples for high throughput sequencing. Although a total of 348 separate gnd sequence types (gSTs) were identified, 188 were likely to be due to sequencing errors. Of the remaining 160 gSTs, 92 did not match those in a database of 319 separate gnd sequences. ‘Animal’ was the main determinant of E. coli diversity with limited impact of sample type or DNA extraction method on intra-host E. coli community variation from faeces and recto-anal mucosal swab samples. This culture-independent study has addressed the difficulties of quantifying bacterial intra-species diversity and revealed that, whilst individual animals may harbour >50 separate E. coli strains, communities are dominated by <10 strains alongside a large pool of subdominant strains present at low abundances. This method will be useful for characterising the diversity and population structure of E. coli in experimental studies designed to assess the impact of interventions on the gut microbiome.

Unveiling the Hybrid Genome Structure of Escherichia coli RR1 (HB101 RecA+)

There have been extensive genome sequencing studies for Escherichia coli strains, particularly for pathogenic isolates, because fast determination of pathogenic potential and/or drug resistance and their propagation routes is crucial. For laboratory E. coli strains, however, genome sequence information is limited except for several well-known strains. We determined the complete genome sequence of laboratory E. coli strain RR1 (HB101 RecA⁺), which has long been used as a general cloning host. A hybrid genome sequence of K-12 MG1655 and B BL21(DE3) was constructed based on the initial mapping of Illumina HiSeq reads to each reference, and iterative rounds of read mapping, variant detection, and consensus extraction were carried out. Finally, PCR and Sanger sequencing-based finishing were applied to resolve non-single nucleotide variant regions with aberrant read depths and breakpoints, most of them resulting from prophages and insertion sequence transpositions that are not present in the reference genome sequence. We found that 96.9% of the RR1 genome is derived from K-12, and identified exact crossover junctions between K-12 and B genomic fragments. However, because RR1 has experienced a series of genetic manipulations since branching from the common ancestor, it has a set of mutations different from those found in K-12 MG1655. As well as identifying all known genotypes of RR1 on the basis of genomic context, we found novel mutations. Our results extend current knowledge of the genotype of RR1 and its relatives, and provide insights into the pedigree, genomic background, and physiology of common laboratory strains.

Genetics and evolution of Yersinia pseudotuberculosis O-specific polysaccharides: a novel pattern of O-antigen diversity

O-antigen polysaccharide is a major immunogenic feature of the lipopolysaccharide of Gram-negative bacteria, and most species produce a large variety of forms that differ substantially from one another. There are 18 known O-antigen forms in the Yersinia pseudotuberculosis complex, which are typical in being composed of multiple copies of a short oligosaccharide called an O unit. The O-antigen gene clusters are located between the hemH and gsk genes, and are atypical as 15 of them are closely related, each having one of five downstream gene modules for alternative main-chain synthesis, and one of seven upstream modules for alternative side-branch sugar synthesis. As a result, many of the genes are in more than one gene cluster. The gene order in each module is such that, in general, the earlier a gene product functions in O-unit synthesis, the closer the gene is to the 5΄ end for side-branch modules or the 3΄ end for main-chain modules. We propose a model whereby natural selection could generate the observed pattern in gene order, a pattern that has also been observed in other species.

Population Phylogenomics of Extraintestinal Pathogenic Escherichia coli

The emergence of genomics over the last 10 years has provided new insights into the evolution and virulence of extraintestinal Escherichia coli. By combining population genetics and phylogenetic approaches to analyze whole-genome sequences, it became possible to link genomic features to specific phenotypes, such as the ability to cause urinary tract infections. An E. coli chromosome can vary extensively in length, ranging from 4.3 to 6.2 Mb, encoding 4,084 to 6,453 proteins. This huge diversity is structured as a set of less than 2,000 genes (core genome) that are conserved between all the strains and a set of variable genes. Based on the core genome, the history of the species can be reliably reconstructed, revealing the recent emergence of phylogenetic groups A and B1 and the more ancient groups B2, F, and D. Urovirulence is most often observed in B2/F/D group strains and is a multigenic process involving numerous combinations of genes and specific alleles with epistatic interactions, all leading down multiple evolutionary paths. The genes involved mainly code for adhesins, toxins, iron capture systems, and protectins, as well as metabolic pathways and mutation-rate-control systems. However, the barrier between commensal and uropathogenic E. coli strains is difficult to draw as the factors that are responsible for virulence have probably also been selected to allow survival of E. coli as a commensal in the intestinal tract. Genomic studies have also demonstrated that infections are not the result of a unique and stable isolate, but rather often involve several isolates with variable levels of diversity that dynamically changes over time.

Pangenome Evolution in the Marine Bacterium Alteromonas

We have examined a collection of the free-living marine bacterium Alteromonas genomes with cores diverging in average nucleotide identities ranging from 99.98% to 73.35%, i.e., from microbes that can be considered members of a natural clone (like in a clinical epidemiological outbreak) to borderline genus level. The genomes were largely syntenic allowing a precise delimitation of the core and flexible regions in each. The core was 1.4 Mb (ca. 30% of the typical strain genome size). Recombination rates along the core were high among strains belonging to the same species (37.7-83.7% of all nucleotide polymorphisms) but they decreased sharply between species (18.9-5.1%). Regarding the flexible genome, its main expansion occurred within the boundaries of the species, i.e., strains of the same species already have a large and diverse flexible genome. Flexible regions occupy mostly fixed genomic locations. Four large genomic islands are involved in the synthesis of strain-specific glycosydic receptors that we have called glycotypes. These genomic regions are exchanged by homologous recombination within and between species and there is evidence for their import from distant taxonomic units (other genera within the family). In addition, several hotspots for integration of gene cassettes by illegitimate recombination are distributed throughout the genome. They code for features that give each clone specific properties to interact with their ecological niche and must flow fast throughout the whole genus as they are found, with nearly identical sequences, in different species. Models for the generation of this genomic diversity involving phage predation are discussed.

The O-antigen mediates differential survival of Salmonella against communities of natural predators

Antigenically distinct members of bacterial species can be differentially distributed in the environment. Predators known to consume antigenically distinct prey with different efficiencies are also differentially distributed. Here we show that antigenically distinct, but otherwise isogenic and physiologically indistinct, strains of Salmonella enterica show differential survival in natural soil, sediment and intestinal environments, where they would face a community of predators. Decline in overall cell numbers is attenuated by factors that inhibit the action of predators, including heat and antiprotozoal and antihelminthic drugs. Moreover, the fitness of strains facing these predators - calculated by comparing survival with and without treatments attenuating predator activity - varies between environments. These results suggest that relative survival in natural environments is arbitrated by communities of natural predators whose feeding preferences, if not species composition, vary between environments. These data support the hypothesis that survival against natural predators may drive the differential distribution of bacteria among microenvironments.

Molecular Characterization of Human Atypical Sorbitol-Fermenting Enteropathogenic Escherichia coli O157 Reveals High Diversity

Alongside the well-characterized enterohemorrhagic Escherichia coli (EHEC) O157:H7, serogroup O157 comprises sorbitol-fermenting typical and atypical enteropathogenic E. coli (EPEC/aEPEC) strains that carry the intimin-encoding gene eae but not Shiga toxin-encoding genes (stx). Since little is known about these pathogens, we characterized 30 clinical isolates from patients with hemolytic uremic syndrome (HUS) or uncomplicated diarrhea with respect to their flagellin gene (fliC) type and multilocus sequence type (MLST). Moreover, we applied whole-genome sequencing (WGS) to determine the phylogenetic relationship with other eae-positive EHEC serotypes and the composition of the rfbO157 region. fliC typing resulted in five fliC types (H7, H16, H34, H39, and H45). Isolates of each fliC type shared a unique ST. In comparison to the 42 HUS-associated E. coli (HUSEC) strains, only the stx-negative isolates with fliCH7 shared their ST with EHEC O157:H7/H(-) strains. With the exception of one O157:H(-) fliCH16 isolate, HUS was exclusively associated with fliCH7. WGS corroborated the separation of the fliCH7 isolates, which were closely related to the EHEC O157:H7/H(-) isolates, and the diverse group of isolates exhibiting different fliC types, indicating independent evolution of the different serotypes. This was also supported by the heterogeneity within the rfbO157 region that exhibited extensive recombinations. The genotypic subtypes and distribution of clinical symptoms suggested that the stx-negative O157 strains with fliCH7 were originally EHEC strains that lost stx The remaining isolates form a distinct and diverse group of atypical EPEC isolates that do not possess the full spectrum of virulence genes, underlining the importance of identifying the H antigen for clinical risk assessment.

Impermanence of bacterial clones

Bacteria reproduce asexually and pass on a single genome copied from the parent, a reproductive mode that assures the clonal descent of progeny; however, a truly clonal bacterial species is extremely rare. The signal of clonality can be interrupted by gene uptake and exchange, initiating homologous recombination that results in the unique sequence of one clone being incorporated into another. Because recombination occurs sporadically and on local scales, these events are often difficult to recognize, even when considering large samples of completely sequenced genomes. Moreover, several processes can produce the appearance of clonality in populations that undergo frequent recombination. The rates and consequences of recombination have been studied in Escherichia coli for over 40 y, and, during this time, there have been several shifting views of its clonal status, population structure, and rates of gene exchange. We reexamine the studies and retrace the evolution of the methods that have assessed the extent of DNA flux, largely focusing on its impact on the E. coli genome.

Novel recA-Independent Horizontal Gene Transfer in Escherichia coli K-12

In bacteria, mechanisms that incorporate DNA into a genome without strand-transfer proteins such as RecA play a major role in generating novelty by horizontal gene transfer. We describe a new illegitimate recombination event in Escherichia coli K-12: RecA-independent homologous replacements, with very large (megabase-length) donor patches replacing recipient DNA. A previously uncharacterized gene (yjiP) increases the frequency of RecA-independent replacement recombination. To show this, we used conjugal DNA transfer, combining a classical conjugation donor, HfrH, with modern genome engineering methods and whole genome sequencing analysis to enable interrogation of genetic dependence of integration mechanisms and characterization of recombination products. As in classical experiments, genomic DNA transfer begins at a unique position in the donor, entering the recipient via conjugation; antibiotic resistance markers are then used to select recombinant progeny. Different configurations of this system were used to compare known mechanisms for stable DNA incorporation, including homologous recombination, F'-plasmid formation, and genome duplication. A genome island of interest known as the immigration control region was specifically replaced in a minority of recombinants, at a frequency of 3 X 10(-12) CFU/recipient per hour.

Recombinant transfer in the basic genome of Escherichia coli

An approximation to the ∼4-Mbp basic genome shared by 32 strains of Escherichia coli representing six evolutionary groups has been derived and analyzed computationally. A multiple alignment of the 32 complete genome sequences was filtered to remove mobile elements and identify the most reliable ∼90% of the aligned length of each of the resulting 496 basic-genome pairs. Patterns of single base-pair mutations (SNPs) in aligned pairs distinguish clonally inherited regions from regions where either genome has acquired DNA fragments from diverged genomes by homologous recombination since their last common ancestor. Such recombinant transfer is pervasive across the basic genome, mostly between genomes in the same evolutionary group, and generates many unique mosaic patterns. The six least-diverged genome pairs have one or two recombinant transfers of length ∼40-115 kbp (and few if any other transfers), each containing one or more gene clusters known to confer strong selective advantage in some environments. Moderately diverged genome pairs (0.4-1% SNPs) show mosaic patterns of interspersed clonal and recombinant regions of varying lengths throughout the basic genome, whereas more highly diverged pairs within an evolutionary group or pairs between evolutionary groups having >1.3% SNPs have few clonal matches longer than a few kilobase pairs. Many recombinant transfers appear to incorporate fragments of the entering DNA produced by restriction systems of the recipient cell. A simple computational model can closely fit the data. Most recombinant transfers seem likely to be due to generalized transduction by coevolving populations of phages, which could efficiently distribute variability throughout bacterial genomes.

Genome-wide recombination drives diversification of epidemic strains of Acinetobacter baumannii

Acinetobacter baumannii is an emerging human pathogen and a significant cause of nosocomial infections among hospital patients worldwide. The enormous increase in multidrug resistance among hospital isolates and the recent emergence of pan-drug-resistant strains underscores the urgency to understand how A. baumannii evolves in hospital environments. To this end, we undertook a genomic study of a polyclonal outbreak of multidrug-resistant A. baumannii at the research-based National Institutes of Health Clinical Center. Comparing the complete genome sequences of the three dominant outbreak strain types enabled us to conclude that, despite all belonging to the same epidemic lineage, the three strains diverged before their arrival at the National Institutes of Health. The simultaneous presence of three divergent strains from this lineage supports its increasing prevalence in international hospitals and suggests an ongoing adaptation to the hospital environment. Further genomic comparisons uncovered that much of the diversification that occurred since the divergence of the three outbreak strains was mediated by homologous recombination across 20% of their genomes. Inspection of recombinant regions revealed that several regions were associated with either the loss or swapping out of genes encoding proteins that are exposed to the cell surface or that synthesize cell-surface molecules. Extending our analysis to a larger set of international clinical isolates revealed a previously unappreciated ability of A. baumannii to vary surface molecules through horizontal gene transfer, with subsequent intraspecies dissemination by homologous recombination. These findings have immediate implications in surveillance, prevention, and treatment of A. baumannii infections.

Screening of Escherichia coli species biodiversity reveals new biofilm-associated antiadhesion polysaccharides

Bacterial biofilms often form multispecies communities in which complex but ill-understood competition and cooperation interactions occur. In light of the profound physiological modifications associated with this lifestyle, we hypothesized that the biofilm environment might represent an untapped source of natural bioactive molecules interfering with bacterial adhesion or biofilm formation. We produced cell-free solutions extracted from in vitro mature biofilms formed by 122 natural Escherichia coli isolates, and we screened these biofilm extracts for antiadhesion molecules active on a panel of Gram-positive and Gram-negative bacteria. Using this approach, we showed that 20% of the tested biofilm extracts contained molecules that antagonize bacterial growth or adhesion. We characterized a compound, produced by a commensal animal E. coli strain, for which activity is detected only in biofilm extract. Biochemical and genetic analyses showed that this compound corresponds to a new type of released high-molecular-weight polysaccharide whose biofilm-associated production is regulated by the RfaH protein. We demonstrated that the antiadhesion activity of this polysaccharide was restricted to Gram-positive bacteria and that its production reduced susceptibility to invasion and provided rapid exclusion of Staphylococcus aureus from mixed E. coli and S. aureus biofilms. Our results therefore demonstrate that biofilms contain molecules that contribute to the dynamics of mixed bacterial communities and that are not or only poorly detected in unconcentrated planktonic supernatants. Systematic identification of these compounds could lead to strategies that limit pathogen surface colonization and reduce the burden associated with the development of bacterial biofilms on medical devices.

We sought to demonstrate that bacterial biofilms are reservoirs for unknown molecules that antagonize bacterial adhesion. The use of natural strains representative of Escherichia coli species biodiversity showed that nonbiocidal antiadhesion polysaccharides are frequently found in mature biofilm extracts (bacterium-free suspensions which contain soluble molecules produced within the biofilm). Release of an antiadhesion polysaccharide confers a competitive advantage upon the producing strain against clinically relevant pathogens such as Staphylococcus aureus. Hence, exploring the biofilm environment provides a better understanding of bacterial interactions within complex communities and could lead to improved control of pathogen colonization.

Population structure and evolution of pathogenicity of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is an enteric human pathogen but is widespread in the environment. Pathogenicity is determined by a number of virulence factors, including the virulence plasmid pYV, the high-pathogenicity island (HPI), and the Y. pseudotuberculosis-derived mitogen (YPM), a superantigen. The presence of the 3 virulence factors varies among Y. pseudotuberculosis isolates. We developed a multilocus sequence typing (MLST) scheme to address the population structure of Y. pseudotuberculosis and the evolution of its pathogenicity. The seven housekeeping genes selected for MLST were mdh, recA, sucA, fumC, aroC, pgi, and gyrB. An MLST analysis of 83 isolates of Y. pseudotuberculosis, representing 19 different serotypes and six different genetic groups, identified 61 sequence types (STs) and 12 clonal complexes. Out of 26 allelic changes that occurred in the 12 clonal complexes, 13 were mutational events while 13 were recombinational events, indicating that recombination and mutation contributed equally to the diversification of the clonal complexes. The isolates were separated into 2 distinctive clusters, A and B. Cluster A is the major cluster, with 53 STs (including Y. pestis strains), and is distributed worldwide, while cluster B is restricted to the Far East. The YPM gene is widely distributed on the phylogenetic tree, with ypmA in cluster A and ypmB in cluster B. pYV is present in cluster A only but is sporadically absent in some cluster A isolates. In contrast, an HPI is present only in a limited number of lineages and must be gained by lateral transfer. Three STs carry all 3 virulence factors and can be regarded as high-pathogenicity clones. Isolates from the same ST may not carry all 3 virulence factors, indicating frequent gain or loss of these factors. The differences in pathogenicity among Y. pseudotuberculosis strains are likely due to the variable presence and instability of the virulence factors.

The genome sequence of E. coli W (ATCC 9637): comparative genome analysis and an improved genome-scale reconstruction of E. coli

BACKGROUND

Escherichia coli is a model prokaryote, an important pathogen, and a key organism for industrial biotechnology. E. coli W (ATCC 9637), one of four strains designated as safe for laboratory purposes, has not been sequenced. E. coli W is a fast-growing strain and is the only safe strain that can utilize sucrose as a carbon source. Lifecycle analysis has demonstrated that sucrose from sugarcane is a preferred carbon source for industrial bioprocesses.

RESULTS

We have sequenced and annotated the genome of E. coli W. The chromosome is 4,900,968 bp and encodes 4,764 ORFs. Two plasmids, pRK1 (102,536 bp) and pRK2 (5,360 bp), are also present. W has unique features relative to other sequenced laboratory strains (K-12, B and Crooks): it has a larger genome and belongs to phylogroup B1 rather than A. W also grows on a much broader range of carbon sources than does K-12. A genome-scale reconstruction was developed and validated in order to interrogate metabolic properties.

CONCLUSIONS

The genome of W is more similar to commensal and pathogenic B1 strains than phylogroup A strains, and therefore has greater utility for comparative analyses with these strains. W should therefore be the strain of choice, or 'type strain' for group B1 comparative analyses. The genome annotation and tools created here are expected to allow further utilization and development of E. coli W as an industrial organism for sucrose-based bioprocesses. Refinements in our E. coli metabolic reconstruction allow it to more accurately define E. coli metabolism relative to previous models.

Impact of recombination on bacterial evolution

Genetic exchange plays a defining role in the evolution of many bacteria. The recent accumulation of nucleotide sequence data from multiple members of diverse bacterial genera has facilitated comparative studies that have revealed many features of this process. Here we focus on genetic exchange that has involved homologous recombination and illustrate how nucleotide sequence data have furthered our understanding of: (i) the frequency of recombination; (ii) the impact of recombination in different parts of the genome; and (iii) patterns of gene flow within bacterial populations. Summarizing the results obtained for a range of bacteria, we survey evidence indicating that the extent and nature of recombination vary widely among microbiological species and often among lineages assigned to the same microbiological species. These results have important implications in studies ranging from epidemiological investigations to examination of the bacterial species problem.

The population genetics of commensal Escherichia coli

The primary habitat of Escherichia coli is the vertebrate gut, where it is the predominant aerobic organism, living in symbiosis with its host. Despite the occurrence of recombination events, the population structure is predominantly clonal, allowing the delineation of major phylogenetic groups. The genetic structure of commensal E. coli is shaped by multiple host and environmental factors, and the determinants involved in the virulence of the bacteria may in fact reflect adaptation to commensal habitats. A better characterization of the commensal niche is necessary to understand how a useful commensal can become a harmful pathogen. In this Review we describe the population structure of commensal E. coli, the factors involved in the spread of different strains, how the bacteria can adapt to different niches and how a commensal lifestyle can evolve into a pathogenic one.

Derivation of Escherichia coli O157:H7 from its O55:H7 precursor

There are 29 E. coli genome sequences available, mostly related to studies of species diversity or mode of pathogenicity, including two genomes of the well-known O157:H7 clone. However, there have been no genome studies of closely related clones aimed at exposing the details of evolutionary change. Here we sequenced the genome of an O55:H7 strain, closely related to the major pathogenic O157:H7 clone, with published genome sequences, and undertook comparative genomic and proteomic analysis. We were able to allocate most differences between the genomes to individual mutations, recombination events, or lateral gene transfer events, in specific lineages. Major differences include a type II secretion system present only in the O55:H7 chromosome, fewer type III secretion system effectors in O55:H7, and 19 phage genomes or phagelike elements in O55:H7 compared to 23 in O157:H7, with only three common to both. Many other changes were found in both O55:H7 and O157:H7 lineages, but in general there has been more change in the O157:H7 lineages. For example, we found 50% more synonymous mutational substitutions in O157:H7 compared to O55:H7. The two strains also diverged at the proteomic level. Mutational synonymous SNPs were used to estimate a divergence time of 400 years using a new clock rate, in contrast to 14,000 to 70,000 years using the traditional clock rates. The same approaches were applied to three closely related extraintestinal pathogenic E. coli genomes, and similar levels of mutation and recombination were found. This study revealed for the first time the full range of events involved in the evolution of the O157:H7 clone from its O55:H7 ancestor, and suggested that O157:H7 arose quite recently. Our findings also suggest that E. coli has a much lower frequency of recombination relative to mutation than was observed in a comparable study of a Vibrio cholerae lineage.

Genome sequences of Escherichia coli B strains REL606 and BL21(DE3)

Escherichia coli K-12 and B have been the subjects of classical experiments from which much of our understanding of molecular genetics has emerged. We present here complete genome sequences of two E. coli B strains, REL606, used in a long-term evolution experiment, and BL21(DE3), widely used to express recombinant proteins. The two genomes differ in length by 72,304 bp and have 426 single base pair differences, a seemingly large difference for laboratory strains having a common ancestor within the last 67 years. Transpositions by IS1 and IS150 have occurred in both lineages. Integration of the DE3 prophage in BL21(DE3) apparently displaced a defective prophage in the lambda attachment site of B. As might have been anticipated from the many genetic and biochemical experiments comparing B and K-12 over the years, the B genomes are similar in size and organization to the genome of E. coli K-12 MG1655 and have >99% sequence identity over approximately 92% of their genomes. E. coli B and K-12 differ considerably in distribution of IS elements and in location and composition of larger mobile elements. An unexpected difference is the absence of a large cluster of flagella genes in B, due to a 41 kbp IS1-mediated deletion. Gene clusters that specify the LPS core, O antigen, and restriction enzymes differ substantially, presumably because of horizontal transfer. Comparative analysis of 32 independently isolated E. coli and Shigella genomes, both commensals and pathogenic strains, identifies a minimal set of genes in common plus many strain-specific genes that constitute a large E. coli pan-genome.

Understanding the differences between genome sequences of Escherichia coli B strains REL606 and BL21(DE3) and comparison of the E. coli B and K-12 genomes

Each difference between the genome sequences of Escherichia coli B strains REL606 and BL21(DE3) can be interpreted in light of known laboratory manipulations plus a gene conversion between ribosomal RNA operons. Two treatments with 1-methyl-3-nitro-1-nitrosoguanidine in the REL606 lineage produced at least 93 single-base-pair mutations ( approximately 90% GC-to-AT transitions) and 3 single-base-pair GC deletions. Two UV treatments in the BL21(DE3) lineage produced only 4 single-base-pair mutations but 16 large deletions. P1 transductions from K-12 into the two B lineages produced 317 single-base-pair differences and 9 insertions or deletions, reflecting differences between B DNA in BL21(DE3) and integrated restriction fragments of K-12 DNA inherited by REL606. Two sites showed selective enrichment of spontaneous mutations. No unselected spontaneous single-base-pair mutations were evident. The genome sequences revealed that a progenitor of REL606 had been misidentified, explaining initially perplexing differences. Limited sequencing of other B strains defined characteristic properties of B and allowed assembly of the inferred genome of the ancestral B of Delbrück and Luria. Comparison of the B and K-12 genomes shows that more than half of the 3793 proteins of their basic genomes are predicted to be identical, although approximately 310 appear to be functional in either B or K-12 but not in both. The ancestral basic genome appears to have had approximately 4039 coding sequences occupying approximately 4.0 Mbp. Repeated horizontal transfer from diverged Escherichia coli genomes and homologous recombination may explain the observed variable distribution of single-base-pair differences. Fifteen sites are occupied by phage-related elements, but only six by comparable elements at the same site. More than 50 sites are occupied by IS elements in both B and K, 16 in common, and likely founding IS elements are identified. A signature of widespread cryptic phage P4-type mobile elements was identified. Complex deletions (dense clusters of small deletions and substitutions) apparently removed nonessential genes from approximately 30 sites in the basic genomes.

Looking for Darwin's footprints in the microbial world

As we observe the 200th anniversary of Charles Darwin's birth, microbiologists interested in the application of Darwin's ideas to the microscopic world have a lot to celebrate: an emerging picture of the (mostly microbial) Tree of Life at ever-increasing resolution, an understanding of horizontal gene transfer as a driving force in the evolution of microbes, and thousands of complete genome sequences to help formulate and refine our theories. At the same time, quantitative models of the microevolutionary processes shaping microbial populations remain just out of reach, a point that is perhaps most dramatically illustrated by the lack of consensus on how (or even whether) to define bacterial species. Here, we summarize progress and prospects in bacterial population genetics, with an emphasis on detecting the footprint of positive Darwinian selection in microbial genomes.

Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the ∼18,000 families of orthologous genes, we found ∼2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

Although abundant knowledge has been accumulated regarding the E. coli laboratory strain K-12, little is known about the evolutionary trajectories that have driven the high diversity observed among natural isolates of the species, which encompass both commensal and highly virulent intestinal and extraintestinal pathogenic strains. We have annotated or re-annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Although recombination rates are much higher than mutation rates, we were able to reconstruct a robust phylogeny based on the ∼2,000 genes common to all strains. Based on this phylogeny, we established the evolutionary scenario of gains and losses of thousands of specific genes, identifying functional classes under opposite selection pressures. This genome flux is confined to very few positions in the chromosome, which are the same for every genome. Notably, we identified few or no extraintestinal virulence-specific genes. We also defined a long-scale structure of recombination in the genome with lower recombination rates at the terminus of replication. These findings demonstrate that, despite a very high gene flow, genes can co-exist in an organised genome.

Population genomics and the bacterial species concept

In recent years, the importance of horizontal gene transfer (HGT) in bacterial evolution has been elevated to such a degree that many bacteriologists now question the very existence of bacterial species. If gene transfer is as rampant as comparative genomic studies have suggested, how could bacterial species survive such genomic fluidity? And yet, most bacteriologists recognize, and name, as species, clusters of bacterial isolates that share complex phenotypic properties. The Core Genome Hypothesis (CGH) has been proposed to explain this apparent paradox of fluid bacterial genomes associated with stable phenotypic clusters. It posits that there is a core of genes responsible for maintaining the species-specific phenotypic clusters observed throughout bacterial diversity and argues that, even in the face of substantial genomic fluidity, bacterial species can be rationally identified and named.

Recombination and positive selection contribute to evolution of Listeria monocytogenes inlA

The surface molecule InlA interacts with E-cadherin to promote invasion of Listeria monocytogenes into selected host cells. DNA sequencing of inlA for 40 L. monocytogenes isolates revealed 107 synonymous and 45 nonsynonymous substitutions. A frameshift mutation in a homopolymeric tract encoding part of the InlA signal peptide was identified in three lineage II isolates, which also showed reduced ability to invade human intestinal epithelial cells. Phylogenies showed clear separation of inlA sequences into lineages I and II. Thirteen inlA recombination events, predominantly involving lineage II strains as recipients (12 events), were detected and a number of amino acid residues were shown to be under positive selection. Four of the 45 non-synonymous changes were found to be under positive selection with posterior probabilities >95 %. Mapping of polymorphic and positively selected amino acid sites on the partial crystal structure for InlA showed that the internalin surface of the leucine-rich repeat (LRR) region that faces the InlA receptor E-cadherin does not include any polymorphic sites; all polymorphic and positively selected amino acids mapped to the outer face of the LRR region or to other InlA regions. The data show that (i) inlA is highly polymorphic and evolution of inlA involved a considerable number of recombination events in lineage II isolates; (ii) positive selection at specific amino acid sites appears to contribute to evolution of inlA, including fixation of recombinant events; and (iii) single-nucleotide deletions in a lineage II-specific 3' homopolymeric tract in inlA lead to complete loss of InlA or to production of truncated InlA, which conveys reduced invasiveness. In conclusion, inlA has a complex evolutionary history, which is consistent with L. monocytogenes' natural history as an environmental pathogen with broad host-range, including its adaptation to environments and hosts where different inlA alleles may provide a selective advantage or where inlA may not be required.

Differential Salmonella survival against communities of intestinal amoebae

Predation from intestinal amoebae may provide selective pressure for the maintenance of high genetic diversity at the Salmonella enterica rfb locus, whereby serovars better escape predators in particular environments depending on the O-antigens they express. Here, the hypothesis that amoebae from a particular intestinal environment collectively prefer one serovar over another is tested. Collections of Acanthamoeba, Tetramitus, Naegleria and Hartmannella were isolated from the intestinal tracts of several vertebrate hosts, including bullfrog tadpoles, goldfish, turtles and bearded dragons, and their feeding preferences were determined. Congeneric amoebae from the same environment had significantly similar feeding preferences. Strikingly, even unrelated amoebae - such as Naegleria and Tetramitus from goldfish - also had significantly similar feeding preferences. Yet amoebae isolated from different environments showed no similarity in prey choice. Thus, feeding preferences of amoebae appear to reflect their environment, not their taxonomic relationships. A mechanism mediating this phenotypic convergence is discussed.

Temporal fragmentation of speciation in bacteria

Because bacterial recombination involves the occasional transfer of small DNA fragments between strains, different sets of niche-specific genes may be maintained in populations that freely recombine at other loci. Therefore, genetic isolation may be established at different times for different chromosomal regions during speciation as recombination at niche-specific genes is curtailed. To test this model, we separated sequence divergence into rate and time components, revealing that different regions of the Escherichia coli and Salmonella enterica chromosomes diverged over a approximately 70-million-year period. Genetic isolation first occurred at regions carrying species-specific genes, indicating that physiological distinctiveness between the nascent Escherichia and Salmonella lineages was maintained for tens of millions of years before the complete genetic isolation of their chromosomes.

Clonal analysis reveals high rate of structural mutations in fimbrial adhesins of extraintestinal pathogenic Escherichia coli

Type 1 fimbriae of Escherichia coli mediate mannose-specific adhesion to host epithelial surfaces and consist of a major, antigenically variable pilin subunit, FimA, and a minor, structurally conserved adhesive subunit, FimH, located on the fimbrial tip. We have analysed the variability of fimA and fimH in strains of vaginal and other origin that belong to one of the most prominent clonal groups of extraintestinal pathogenic E. coli, comprised of O1:K1-, O2:K1- and O18:K1-based serotypes. Multiple locus sequence typing (MLST) of this group revealed that the strains have identical (at all but one nucleotide position) eight housekeeping loci around the genome and belong to the ST95 complex defined by the publicly available E. coli MLST database. Multiple highly diverse fimA alleles have been introduced into the ST95 clonal complex via horizontal transfer, at a frequency comparable to that of genes defining the major O- and H-antigens. However, no further significant FimA diversification has occurred via point mutation after the transfers. In contrast, while fimH alleles also move horizontally (along with the fimA loci), they acquire point amino acid replacements at a higher rate than either housekeeping genes or fimA. These FimH mutations enhance binding to monomannose receptors and bacterial tropism for human vaginal epithelium. A similar pattern of rapid within-clonal structural evolution of the adhesive, but not pilin, subunit is also seen, respectively, in papG and papA alleles of the di-galactose-specific P-fimbriae. Thus, while structurally diverse pilin subunits of E. coli fimbriae are under selective pressure for frequent horizontal transfer between clones, the adhesive subunits of extraintestinal E. coli are under strong positive selection (Dn/Ds > 1 for fimH and papG) for functionally adaptive amino acid replacements.

Evolutionary changes of the flhDC flagellar master operon in Shigella strains

Shigella strains are nonmotile. The master operon of flagellar synthesis, flhDC, was analyzed for genetic damage in 46 Shigella strains representing all known serotypes. In 11 strains (B1, B3, B6, B8, B10, B18, D5, F1B, D10, F3A, and F3C) the flhDC operon was completely deleted. PCR and sequence analysis of the flhDC region of the remaining 35 strains revealed many insertions or deletions associated with insertion sequences, and the majority of the strains were found to be defective in their flhDC genes. As these genes also play a role in regulation of non-flagellar genes, the loss may have other consequences or be driven by selection pressures other than those against flagellar motility. It has been suggested that Shigella strains fall mostly into three clusters within Escherichia coli, with five outlier strains, four of which are also within E. coli (G. M. Pupo, R. Lan, and P. R. Reeves, Proc. Natl. Acad. Sci. USA 97:10567-10572, 2000). The distribution of genetic changes in the flhDC region correlated very well with the three clusters and outlier strains found using housekeeping gene DNA sequences, enabling us to follow the sequence of mutational change in the flhDC locus. Two cluster 2 strains were found to have unique flhDC sequences, which are most probably due to recombination during the exchange of the adjacent O-antigen gene clusters.

Codon usage comparison of novel genes in clinical isolates of Haemophilus influenzae

A similarity statistic for codon usage was developed and used to compare novel gene sequences found in clinical isolates of Haemophilus influenzae with a reference set of 80 prokaryotic, eukaryotic and viral genomes. These analyses were performed to obtain an indication as to whether individual genes were Haemophilus-like in nature, or if they probably had more recently entered the H.influenzae gene pool via horizontal gene transfer from other species. The average and SD values were calculated for the similarity statistics from a study of the set of all genes in the H.influenzae Rd reference genome that encoded proteins of 100 amino acids or longer. Approximately 80% of Rd genes gave a statistic indicating that they were most like other Rd genes. Genes displaying codon usage statistics >1 SD above this range were either considered part of the highly expressed group of H.influenzae genes, or were considered of foreign origin. An alternative determinant for identifying genes of foreign origin was when the similarity statistics produced a value that was much closer to a non-H.influenzae reference organism than to any of the Haemophilus species contained in the reference set. Approximately 65% of the novel sequences identified in the H.influenzae clinical isolates displayed codon usages most similar to Haemophilus sp. The remaining novel sequences produced similarity statistics closer to one of the other reference genomes thereby suggesting that these sequences may have entered the H.influenzae gene pool more recently via horizontal transfer.