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.

Clostridium difficile is an emergent pathogen, and the most common cause of nosocomial diarrhea. In an effort to understand the role of small noncoding RNAs (sRNAs) in C. difficile physiology and pathogenesis, we used an in silico approach to identify 511 sRNA candidates in both intergenic and coding regions. In parallel, RNA–seq and differential 5′-end RNA–seq were used for global identification of C. difficile sRNAs and their transcriptional start sites at three different growth conditions (exponential growth phase, stationary phase, and starvation). This global experimental approach identified 251 putative regulatory sRNAs including 94 potential trans riboregulators located in intergenic regions, 91 cis-antisense RNAs, and 66 riboswitches. Expression of 35 sRNAs was confirmed by gene-specific experimental approaches. Some sRNAs, including an antisense RNA that may be involved in control of C. difficile autolytic activity, showed growth phase-dependent expression profiles. Expression of each of 16 predicted c-di-GMP-responsive riboswitches was observed, and experimental evidence for their regulatory role in coordinated control of motility and biofilm formation was obtained. Finally, we detected abundant sRNAs encoded by multiple C. difficile CRISPR loci. These RNAs may be important for C. difficile survival in bacteriophage-rich gut communities. Altogether, this first experimental genome-wide identification of C. difficile sRNAs provides a firm basis for future RNome characterization and identification of molecular mechanisms of sRNA–based regulation of gene expression in this emergent enteropathogen.

In bacteria, the evolution of pathogenicity seems to be the result of the constant arrival of virulence factors (VFs) into the bacterial genome. However, the integration, retention, and/or expression of these factors may be the result of the interaction between the new arriving genes and the bacterial genomic background. To test this hypothesis, a phylogenetic analysis was done on a collection of 98 Escherichia coli/Shigella strains representing the pathogenic and commensal diversity of the species. The distribution of 17 VFs associated to the different E. coli pathovars was superimposed on the phylogenetic tree. Three major types of VFs can be recognized: (1) VFs that arrive and are expressed in different genetic backgrounds (such as VFs associated with the pathovars of mild chronic diarrhea: enteroaggregative, enteropathogenic, and diffusely-adhering E. coli), (2) VFs that arrive in different genetic backgrounds but are preferentially found, associated with a specific pathology, in only one particular background (such as VFs associated with extraintestinal diseases), and (3) VFs that require a particular genetic background for the arrival and expression of their virulence potential (such as VFs associated with pathovars typical of severe acute diarrhea: enterohemorragic, enterotoxigenic, and enteroinvasive E. coli strains). The possibility of a single arrival of VFs by chance, followed by a vertical transmission, was ruled out by comparing the evolutionary histories of some of these VFs to the strain phylogeny. These evidences suggest that important changes in the genome of E. coli have occurred during the diversification of the species, allowing the virulence factors associated with severe acute diarrhea to arrive in the population. Thus, the E. coli genome seems to be formed by an "ancestral" and a "derived" background, each one responsible for the acquisition and expression of different virulence factors.

Enteroaggregative Escherichia coli (EAEC) is recognized as an emerging cause of diarrhea in children and adults worldwide, and recent studies have implicated EAEC in persistent diarrhea in patients infected with human immunodeficiency virus (HIV). In this study, we identified aggregative adhesion fimbria type III (AAF-III) in isolate 55989, a typical EAEC strain. Analysis of the sequence of the plasmid-borne agg-3 gene cluster encoding AAF-III showed this cluster to be closely related to the agg and aaf operons and to the afa operons carried by diffusely adherent pathogenic E. coli. We investigated the adhesion properties of a collection of 25 EAEC strains isolated from HIV-infected patients presenting with persistent diarrhea. We found that a minority of strains (36%) carried sequences similar to those of the agg and aaf operons, which encode AAF-I and AAF-II, respectively. We developed PCR assays specific for the agg-3 operon. In our collection, the frequency of AAF-III strains was similar (12%) to that of AAF-I strains (16%) but higher than that of AAF-II isolates (0%). Differences between EAEC strains in terms of the virulence factors present render detection of these strains difficult with the available DNA probes. Based on comparison of the agg, aaf, and agg-3 operons, we defined an AAF probe internal to the adhesion gene clusters and demonstrated that it was efficient for the identification of EAEC strains. We investigated 32 EAEC isolates, of which only 34.4% were detected with the classical CVD432 probe (detecting pAA virulence plasmids) whereas 65.6% were detected with the AAF probe.

Pathogenic bacteria possess adhesion protein complexes that play essential roles in targeting host cells and in propagating infection. Although each family of adhesion proteins is generally associated with a specific human disease, the Dr family from Escherichia coli is a notable exception, as its members are associated with both diarrheal and urinary tract infections. These proteins are reported to form both fimbrial and afimbrial structures at the bacterial cell surface and target a common host cell receptor, the decay-accelerating factor (DAF or CD55). Using the newly solved three-dimensional structure of AfaE, we have constructed a robust atomic resolution model that reveals the structural basis for assembly by donor strand complementation and for the architecture of capped surface fibers.