Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA^w) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA^w and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in ‘carbohydrate transport and metabolism'' and ‘DNA replication, recombination and repair'' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.     

Specific Detection of Xanthomonas axonopodis pv. dieffenbachiae in Anthurium (Anthurium andreanum) Tissues by Nested PCR

Efficient control of Xanthomonas axonopodis pv. dieffenbachiae, the causal agent of anthurium bacterial blight, requires a sensitive and reliable diagnostic tool. A nested PCR test was developed from a sequence-characterized amplified region marker identified by randomly amplified polymorphic DNA PCR for the detection of X. axonopodis pv. dieffenbachiae. Serological and pathogenicity tests were performed concurrently with the nested PCR test with a large collection of X. axonopodis pv. dieffenbachiae strains that were isolated worldwide and are pathogenic to anthurium and/or other aroids. The internal primer pair directed amplification of the expected product (785 bp) for all 70 X. axonopodis pv. dieffenbachiae strains pathogenic to anthurium tested and for isolates originating from syngonium and not pathogenic to anthurium. This finding is consistent with previous studies which indicated that there is a high level of relatedness between strains from anthurium and strains from syngonium. Strains originating from the two host genera can be distinguished by restriction analysis of the amplification product. No amplification product was obtained with 98 strains of unrelated phytopathogenic bacteria or saprophytic bacteria from the anthurium phyllosphere, except for a weak signal obtained for one X. axonopodis pv. allii strain. Nevertheless, restriction enzyme analysis permitted the two pathovars to be distinguished. The detection threshold obtained with pure cultures or plant extracts (10³ CFU ml⁻¹) allowed detection of the pathogen from symptomless contaminated plants. This test could be a useful diagnostic tool for screening propagation stock plant material and for monitoring international movement of X. axonopodis pv. dieffenbachiae.     

Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA

We used multilocus sequence analysis (MLSA) on a worldwide collection of the plant pathogenic Ralstonia solanacearum (Betaproteobacteria) to retrace its complex evolutionary history. Using genetic imprints left during R. solanacearum evolution, we were able to delineate distinct evolutionary complex displaying contrasting dynamics. Among the phylotypes already described (I, IIA, IIB, III, IV), eight groups of strains with distinct evolutionary patterns, named clades, were identified. From our recombination analysis, we identified 21 recombination events that occurred within and across these lineages. Although appearing the most divergent and ancestral phylotype, phylotype IV was inferred as a gene donor for the majority of the recombination events that we detected. Whereas this phylotype apparently fuelled the species diversity, ongoing diversification was mainly detected within phylotype I, IIA and III. These three groups presented a recent expanding population structure, a high level of homologous recombination and evidences of long-distance migrations. Factors such as adaptation to a specific host or intense trading of infected crops may have promoted this diversification. Whether R. solanacearum lineages will eventually evolve in distinct species remains an open question. The intensification of cropping and increase of geographical dispersion may favour situations of phylotype sympatry and promote higher exchange of key factors for host adaptation from their common genetic pool.     

Genetic Structure and Population Dynamics of Xanthomonas axonopodis pv. manihotis in Colombia from 1995 to 1999

Restriction fragment length polymorphisms (RFLPs) were used to study the population genetics and temporal dynamics of the cassava bacterial pathogen Xanthomonas axonopodis pv. manihotis. The population dynamics were addressed by comparing samples collected from 1995 to 1999 from six locations, spanning four different edaphoclimatic zones (ECZs). Forty-five different X. axonopodis pv. manihotis RFLP types or haplotypes were identified between 1995 and 1999. High genetic diversity of the X. axonopodis pv. manihotis strains was evident within most of the fields sampled. In all but one site, diversity decreased over time within fields. Haplotype frequencies significantly differed over the years in all but one location. Studies of the rate of change of X. axonopodis pv. manihotis populations during the cropping cycle in two sites showed significant changes in the haplotype frequencies but not composition. However, variations in pathotype composition were observed from one year to the next at a single site in ECZs 1 and 2 and new pathotypes were described after 1997 in these ECZs, thus revealing the dramatic change in the pathogen population structure of X. axonopodis pv. manihotis. Disease incidence was used to show the progress of cassava bacterial blight in Colombia during the 5-year period in different ecosystems. Low disease incidence values were correlated with low rainfall in 1997 in ECZ 1.     

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the ‘birthday problem’. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10⁻⁶ s/s/y) and Campylobacter jejuni (3.4 x 10⁻⁶ s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame–analogous to a shared birthday–and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.     

The Recent Recombinant Evolution of a Major Crop Pathogen,Potato virus Y

Potato virus Y (PVY) is a major agricultural disease that reduces crop yields worldwide. Different strains of PVY are associated with differing degrees of pathogenicity, of which the most common and economically important are known to be recombinant. We need to know the evolutionary origins of pathogens to prevent further escalations of diseases, but putatively reticulate genealogies are challenging to reconstruct with standard phylogenetic approaches. Currently available phylogenetic hypotheses for PVY are either limited to non-recombinant strains, represent only parts of the genome, and/or incorrectly assume a strictly bifurcating phylogenetic tree. Despite attempts to date potyviruses in general, no attempt has been made to date the origins of pathogenic PVY. We test whether diversification of the major strains of PVY and recombination between them occurred within the time frame of the domestication and modern cultivation of potatoes. In so doing, we demonstrate a novel extension of a phylogenetic approach for reconstructing reticulate evolutionary scenarios. We infer a well resolved phylogeny of 44 whole genome sequences of PVY viruses, representative of all known strains, using recombination detection and phylogenetic inference techniques. Using Bayesian molecular dating we show that the parental strains of PVY diverged around the time potatoes were first introduced to Europe, that recombination between them only occurred in the last century, and that the multiple recombination events that led to highly pathogenic PVY^NTN occurred within the last 50 years. Disease causing agents are often transported across the globe by humans, with disastrous effects for us, our livestock and crops. Our analytical approach is particularly pertinent for the often small recombinant genomes involved (e.g. HIV/influenza A). In the case of PVY, increased transport of diseased material is likely to blame for uniting the parents of recombinant pathogenic strains: this process needs to be minimised to prevent further such occurrences.     

A Novel Two-Component Response Regulator Links rpf with Biofilm Formation and Virulence of Xanthomonas axonopodis pv. citri

Citrus bacterial canker caused by Xanthomonas axonopodis pv. citri is a serious disease that impacts citrus production worldwide, and X. axonopodis pv. citri is listed as a quarantine pest in certain countries. Biofilm formation is important for the successful development of a pathogenic relationship between various bacteria and their host(s). To understand the mechanisms of biofilm formation by X. axonopodis pv. citri strain XW19, the strain was subjected to transposon mutagenesis. One mutant with a mutation in a two-component response regulator gene that was deficient in biofilm formation on a polystyrene microplate was selected for further study. The protein was designated as BfdR for biofilm formation defective regulator. BfdR from strain XW19 shares 100% amino acid sequence identity with XAC1284 of X. axonopodis pv. citri strain 306 and 30–100% identity with two-component response regulators in various pathogens and environmental microorganisms. The bfdR mutant strain exhibited significantly decreased biofilm formation on the leaf surfaces of Mexican lime compared with the wild type strain. The bfdR mutant was also compromised in its ability to cause canker lesions. The wild-type phenotype was restored by providing pbfdR in trans in the bfdR mutant. Our data indicated that BfdR did not regulate the production of virulence-related extracellular enzymes including amylase, lipase, protease, and lecithinase or the expression of hrpG, rfbC, and katE; however, BfdR controlled the expression of rpfF in XVM2 medium, which mimics cytoplasmic fluids in planta. In conclusion, biofilm formation on leaf surfaces of citrus is important for canker development in X. axonopodis pv. citri XW19. The process is controlled by the two-component response regulator BfdR via regulation of rpfF, which is required for the biosynthesis of a diffusible signal factor.     

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts

Background

Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants.

Results

Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting.

Conclusion

Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1369-8) contains supplementary material, which is available to authorized users.     

A dynamic global equilibrium in carnivoran diversification over 20 million years

The ecological and evolutionary processes leading to present-day biological diversity can be inferred by reconstructing the phylogeny of living organisms, and then modelling potential processes that could have produced this genealogy. A more direct approach is to estimate past processes from the fossil record. The Carnivora (Mammalia) has both substantial extant species richness and a rich fossil record. We compiled species-level data for over 10 000 fossil occurrences of nearly 1400 carnivoran species. Using this compilation, we estimated extinction, speciation and net diversification for carnivorans through the Neogene (22–2 Ma), while simultaneously modelling sampling probability. Our analyses show that caniforms (dogs, bears and relatives) have higher speciation and extinction rates than feliforms (cats, hyenas and relatives), but lower rates of net diversification. We also find that despite continual species turnover, net carnivoran diversification through the Neogene is surprisingly stable, suggesting a saturated adaptive zone, despite restructuring of the physical environment. This result is strikingly different from analyses of carnivoran diversification estimated from extant species alone. Two intervals show elevated diversification rates (13–12 Ma and 4–3 Ma), although the precise causal factors behind the two peaks in carnivoran diversification remain open questions.     

Phylogenetic and Structural Diversity in the Feline Leukemia Virus Env Gene

Feline leukemia virus (FeLV) belongs to the genus Gammaretrovirus, and causes a variety of neoplastic and non-neoplastic diseases in cats. Alteration of viral env sequences is thought to be associated with disease specificity, but the way in which genetic diversity of FeLV contributes to the generation of such variants in nature is poorly understood. We isolated FeLV env genes from naturally infected cats in Japan and analyzed the evolutionary dynamics of these genes. Phylogenetic reconstructions separated our FeLV samples into three distinct genetic clusters, termed Genotypes I, II, and III. Genotype I is a major genetic cluster and can be further classified into Clades 1–7 in Japan. Genotypes were correlated with geographical distribution; Genotypes I and II were distributed within Japan, whilst FeLV samples from outside Japan belonged to Genotype III. These results may be due to geographical isolation of FeLVs in Japan. The observed structural diversity of the FeLV env gene appears to be caused primarily by mutation, deletion, insertion and recombination, and these variants may be generated de novo in individual cats. FeLV interference assay revealed that FeLV genotypes did not correlate with known FeLV receptor subgroups. We have identified the genotypes which we consider to be reliable for evaluating phylogenetic relationships of FeLV, which embrace the high structural diversity observed in our sample. Overall, these findings extend our understanding of Gammaretrovirus evolutionary patterns in the field, and may provide a useful basis for assessing the emergence of novel strains and understanding the molecular mechanisms of FeLV transmission in cats.     

Multilocus sequence analysis of homologous recombination and diversity in Arthrobacter sensu lato named species and glacier-inhabiting strains

Members of the bacterial genus Arthrobacter sensu lato are Gram-positive actinomycetes distributed worldwide and found in numerous environments including soil, water, glacier ice, and sewage. Homologous recombination is an important driving force in bacterial evolution, but its impact on Arthrobacter sensu lato evolution is poorly understood. We evaluated homologous recombination among 41 Arthrobacter sensu lato named species, using multilocus sequence analysis (MLSA). A high level of recombination was found, associated with strong diversification and a reticulate evolutionary pattern of Arthrobacter sensu lato. We also collected a total of 31 cold-adapted Arthrobacter sensu lato strains from two cold glaciers located in northwest China and two temperate glaciers in southwest China, and evaluated their diversity and population structure by MLSA. The glacier strains displayed high diversity, but rates of recombination among the four glacier groups were quite low, indicating that barriers to homologous recombination formed in the past among the populations on different glaciers. Our findings indicate that historical glaciation events shaped the contemporary distributions, taxonomic relationships, and phylogeographic patterns of Arthrobacter sensu lato species on glaciers.     

Phylogenetic relatedness of Xanthomonas axonopodis pv. punicae, the causal agent of bacterial blight of pomegranate based on two loci, 16S rRNA and gyrB

Bacterial blight caused by Xanthomonas axonopodis pv. punicae (Xap) is a major disease in pomegranate (Punica granatum) cultivation in India. The Xap strains from three distinct geographical origins, Delhi, Maharashtra and Andhra Pradesh were studied for their genetic variability and phylogenetic relationship with other Xanthomonads targeting two important loci 16S rRNA and gyrB. All Xap strains showed 100 % sequence conservation in both the loci, suggesting that geographical origin does not necessarily reflect variation to genetic make-up of the Xap. Phylogeny derived from 16S rRNA gene revealed that two Xanthomonas species, Xanthomonas citri subsp. malvacearum DSM 3849 T and X. axonopodis pv. manihotis NCPPB1834 formed a single cluster along with Xap. Further, analysis in the gyrB locus indicated that X. citri subsp. malvacearum shared 99.4 % identity while pathovars X. axonopodis pv. manihotis shared only 95 % identity with the Xap strains. Thus, we established that gyrB was the preferred locus over 16S rRNA gene to discriminate the Xap strains from closely related Xanthomonas species type strains. Nevertheless, our study demonstrated for the first time that pomegranate bacterial blight pathogen is phylogenetically very close to Xanthomonas citri subsp. malvacearum infecting cotton.     

Microevolution of Helicobacter pylori during Prolonged Infection of Single Hosts and within Families

Our understanding of basic evolutionary processes in bacteria is still very limited. For example, multiple recent dating estimates are based on a universal inter-species molecular clock rate, but that rate was calibrated using estimates of geological dates that are no longer accepted. We therefore estimated the short-term rates of mutation and recombination in Helicobacter pylori by sequencing an average of 39,300 bp in 78 gene fragments from 97 isolates. These isolates included 34 pairs of sequential samples, which were sampled at intervals of 0.25 to 10.2 years. They also included single isolates from 29 individuals (average age: 45 years) from 10 families. The accumulation of sequence diversity increased with time of separation in a clock-like manner in the sequential isolates. We used Approximate Bayesian Computation to estimate the rates of mutation, recombination, mean length of recombination tracts, and average diversity in those tracts. The estimates indicate that the short-term mutation rate is 1.4×10⁻⁶ (serial isolates) to 4.5×10⁻⁶ (family isolates) per nucleotide per year and that three times as many substitutions are introduced by recombination as by mutation. The long-term mutation rate over millennia is 5–17-fold lower, partly due to the removal of non-synonymous mutations due to purifying selection. Comparisons with the recent literature show that short-term mutation rates vary dramatically in different bacterial species and can span a range of several orders of magnitude.     

Serial evolutionary networks of within-patient HIV-1 sequences reveal patterns of evolution of X4 strains

Background  

The HIV virus is known for its ability to exploit numerous genetic and evolutionary mechanisms to ensure its proliferation, among them, high replication, mutation and recombination rates. Sliding MinPD, a recently introduced computational method [1], was used to investigate the patterns of evolution of serially-sampled HIV-1 sequence data from eight patients with a special focus on the emergence of X4 strains. Unlike other phylogenetic methods, Sliding MinPD combines distance-based inference with a nonparametric bootstrap procedure and automated recombination detection to reconstruct the evolutionary history of longitudinal sequence data. We present serial evolutionary networks as a longitudinal representation of the mutational pathways of a viral population in a within-host environment. The longitudinal representation of the evolutionary networks was complemented with charts of clinical markers to facilitate correlation analysis between pertinent clinical information and the evolutionary relationships.     

Disentangling the Population Structure and Evolution of the Clam Pathogen Vibrio tapetis

Vibrio tapetis is a fastidious slow-growing microorganism that causes the Brown Ring Disease in clams. Recently, two subspecies for this bacterial pathogen have been proposed. We have developed a multilocus sequence typing scheme and performed evolutionary studies of V. tapetis population using the great majority of isolates of V. tapetis obtained worldwide until now (30 isolates). V. tapetis constitutes a high polymorphic population, showing low diversity indexes and some genetic discontinuity among the isolates. Mutation events are more frequent than recombination, although both are approximately equally important for genetic diversification. In fact, the divergence between subspecies occurred exclusively by mutation but the diversity observed among isolates of the same subspecies appeared to be generated mostly by recombination. Between the subspecies, genetic distance is very high and almost no recurrent gene flow exists. This pathogen displays a non-clonal population structure with an ancient spatial segregation population and some degree of geographical isolation, followed by a population expansion, at least for V. tapetis subsp. tapetis. A database from this study was created and hosted on publmlst.org (http://pubmlst.org/vtapetis/).     

Characterization of ISXax1, a Novel Insertion Sequence Restricted to Xanthomonas axonopodis pv. phaseoli (Variants fuscans and non-fuscans) and Xanthomonas axonopodis pv. vesicatoria

ISXax1 is a novel insertion sequence belonging to the IS256 and Mutator families. Dot blot, Southern blot, and PCR analyses revealed that ISXax1 is restricted to Xanthomonas axonopodis pv. phaseoli (variants fuscans and non-fuscans) and X. axonopodis pv. vesicatoria strains. Directed AFLP also showed that a high degree of polymorphism is associated with ISXax1 insertion in these strains.     

CRISPR Content Correlates with the Pathogenic Potential of Escherichia coli

Guide RNA molecules (crRNA) produced from clustered regularly interspaced short palindromic repeat (CRISPR) arrays, altogether with effector proteins (Cas) encoded by cognate cas (CRISPR associated) genes, mount an interference mechanism (CRISPR-Cas) that limits acquisition of foreign DNA in Bacteria and Archaea. The specificity of this action is provided by the repeat intervening spacer carried in the crRNA, which upon hybridization with complementary sequences enables their degradation by a Cas endonuclease. Moreover, CRISPR arrays are dynamic landscapes that may gain new spacers from infecting elements or lose them for example during genome replication. Thus, the spacer content of a strain determines the diversity of sequences that can be targeted by the corresponding CRISPR-Cas system reflecting its functionality. Most Escherichia coli strains possess either type I-E or I-F CRISPR-Cas systems. To evaluate their impact on the pathogenicity of the species, we inferred the pathotype and pathogenic potential of 126 strains of this and other closely related species and analyzed their repeat content. Our results revealed a negative correlation between the number of I-E CRISPR units in this system and the presence of pathogenicity traits: the median number of repeats was 2.5-fold higher for commensal isolates (with 29.5 units, range 0–53) than for pathogenic ones (12.0, range 0–42). Moreover, the higher the number of virulence factors within a strain, the lower the repeat content. Additionally, pathogenic strains of distinct ecological niches (i.e., intestinal or extraintestinal) differ in repeat counts. Altogether, these findings support an evolutionary connection between CRISPR and pathogenicity in E. coli.     

Genotypic Characterization of the Common Bean Bacterial Blight Pathogens, Xanthomonas axonopodis pv. phaseoli and Xanthomonas axonopodis pv. phaseoli var. fuscans by rep-PCR and PCR–RFLP of the Ribosomal Genes

Common bacterial blight (CBB), caused by Xanthomonas axonopodis pv. phaseoli and X. axonopodis pv. phaseoli var. fuscans is one of the most destructive diseases of common bean worldwide. The interrelatedness, genetic diversity and geographical distribution of the CBB pathogens was assessed using restriction fragment length polymorphism (RFLP) analysis of polymerase chain reaction amplified 16S ribosomal gene, including the 16S–23S intergenic spacer region and repetitive element PCR (rep‐PCR). RFLP profiles generated by the restriction endonucleases MboI, RsaI and HaeIII differentiated X. axonopodis pv. phaseoli from X. axonopodis pv. phaseoli var. fuscans and non‐pathogenic Xanthomonas species associated with common bean. Cluster analysis of rep‐PCR profiles revealed a high level of genetic differentiation (G_ST = 0.56) between the two CBB pathogens, showing that they are genetically distinct. Significant levels of genetic diversity were observed within each strain, indicating that the two bacteria are not clonal. More genetic diversity was observed in X. axonopodis pv. phaseoli (H = 0.134; I = 0.223) than X. axonopodis pv. phaseoli var. fuscans (H = 0.108; I = 0.184). However, no geographical differentiation was evident for either X. axonopodis pv. phaseoli var. fuscans (GST = 0.013) or X. axonopodis pv. phaseoli (G_ST = 0.017). This lack of geographical differentiation has important practical implications, as available host resistance genes are likely to be effective in controlling the disease in diverse geographical areas.     

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.