Impact of homologous and non-homologous recombination in the genomic evolution of Escherichia coli

Recombination-prone bacterial strains form a reservoir from which epidemic clones emerge in agroecosystems

The acquisition of virulence-related genes through horizontal gene transfer can modify the pathogenic profiles of strains and lead to the emergence of new diseases. Xanthomonas arboricola is a bacterial species largely known for the damage it causes to stone and nut fruit trees worldwide. In addition to these host-specific populations called pathovars, many nonpathogenic strains have been identified in this species. Their evolutionary significance in the context of pathogen emergence is unknown. We looked at seven housekeeping genes amplified from 187 pathogenic and nonpathogenic strains isolated from various plants worldwide to analyze population genetics and recombination dynamics. We also examined the dynamics of the gains and losses of genes associated with life history traits (LHTs) during X. arboricola evolution. We discovered that X. arboricola presents an epidemic population structure. Successful pathovars of trees (i.e. pruni, corylina and juglandis) are epidemic clones whose emergence appears to be linked to the acquisition of eight genes coding for Type III effectors. The other strains of this species are part of a recombinant network, within which LHT-associated genes might have been lost. We suggest that nonpathogenic strains, because of their high genetic diversity and propensity for recombination, may promote the emergence of pathogenic strains.

Genomic content typifying a prevalent clade of bovine mastitis-associated Escherichia coli

E. coli represents a heterogeneous population with capabilities to cause disease in several anatomical sites. Among sites that can be colonised is the bovine mammary gland (udder) and a distinct class of mammary pathogenic E. coli (MPEC) has been proposed. MPEC are the principle causative agents of bovine mastitis in well-managed dairy farms, costing producers in the European Union an estimated €2 billion per year. Despite the economic impact, and the threat this disease presents to small and medium sized dairy farmers, the factors which mediate the ability for E. coli to thrive in bovine mammary tissue remain poorly elucidated. Strains belonging to E. coli phylogroup A are most frequently isolated from mastitis. In this paper, we apply a population level genomic analysis to this group of E. coli to uncover genomic signatures of mammary infectivity. Through a robust statistical analysis, we show that not all strains of E. coli are equally likely to cause mastitis, and those that do possess specific gene content that may promote their adaptation and survival in the bovine udder. Through a pan-genomic analysis, we identify just three genetic loci which are ubiquitous in MPEC, but appear dispensable for E. coli from other niches.

Benefit of transferred mutations is better predicted by the fitness of recipients than by their ecological or genetic relatedness

The effect of a mutation depends on its interaction with the genetic background in which it is assessed. Studies in experimental systems have demonstrated that such interactions are common among beneficial mutations and often follow a pattern consistent with declining evolvability of more fit genotypes. However, these studies generally examine the consequences of interactions between a small number of focal mutations. It is not clear, therefore, that findings can be extrapolated to natural populations, where new mutations may be transferred between genetically divergent backgrounds. We build on work that examined interactions between four beneficial mutations selected in a laboratory-evolved population of Escherichia coli to test how they interact with the genomes of diverse natural isolates of the same species. We find that the fitness effect of transferred mutations depends weakly on the genetic and ecological similarity of recipient strains relative to the donor strain in which the mutations were selected. By contrast, mutation effects were strongly inversely correlated to the initial fitness of the recipient strain. That is, there was a pattern of diminishing returns whereby fit strains benefited proportionally less from an added mutation. Our results strengthen the view that the fitness of a strain can be a major determinant of its ability to adapt. They also support a role for barriers of transmission, rather than differential selection of transferred DNA, as an explanation of observed phylogenetically determined patterns of restricted recombination among E. coli strains.

Rates of Lateral Gene Transfer in Prokaryotes: High but Why?

Lateral gene transfer is of fundamental importance to the evolution of prokaryote genomes and has important practical consequences, as evidenced by the rapid dissemination of antibiotic resistance and virulence determinants. Relatively little effort has so far been devoted to explicitly quantifying the rate at which accessory genes are taken up and lost, but it is possible that the combined rate of lateral gene transfer and gene loss is higher than that of point mutation. What evolutionary forces underlie the rate of lateral gene transfer are not well understood. We here use theory developed to explain the evolution of mutation rates to address this question and explore its consequences for the study of prokaryote evolution.

Evolutionary History of the Global Emergence of the Escherichia coli Epidemic Clone ST131

Escherichia colisequence type 131 (ST131) has emerged globally as the most predominant extraintestinal pathogenic lineage within this clinically important species, and its association with fluoroquinolone and extended-spectrum cephalosporin resistance impacts significantly on treatment. The evolutionary histories of this lineage, and of important antimicrobial resistance elements within it, remain unclearly defined. This study of the largest worldwide collection (n= 215) of sequenced ST131E. coliisolates to date demonstrates that the clonal expansion of two previously recognized antimicrobial-resistant clades, C1/H30R and C2/H30Rx, started around 25 years ago, consistent with the widespread introduction of fluoroquinolones and extended-spectrum cephalosporins in clinical medicine. These two clades appear to have emerged in the United States, with the expansion of the C2/H30Rx clade driven by the acquisition of ablaCTX-M-15-containing IncFII-like plasmid that has subsequently undergone extensive rearrangement. Several other evolutionary processes influencing the trajectory of this drug-resistant lineage are described, including sporadic acquisitions of CTX-M resistance plasmids and chromosomal integration ofblaCTX-Mwithin subclusters followed by vertical evolution. These processes are also occurring for another family of CTX-M gene variants more recently observed among ST131, theblaCTX-M-14/14-likegroup. The complexity of the evolutionary history of ST131 has important implications for antimicrobial resistance surveillance, epidemiological analysis, and control of emerging clinical lineages ofE. coli These data also highlight the global imperative to reduce specific antibiotic selection pressures and demonstrate the important and varied roles played by plasmids and other mobile genetic elements in the perpetuation of antimicrobial resistance within lineages.

IMPORTANCE

Escherichia coli, perennially a major bacterial pathogen, is becoming increasingly difficult to manage due to emerging resistance to all preferred antimicrobials. Resistance is concentrated within specificE. colilineages, such as sequence type 131 (ST131). Clarification of the genetic basis for clonally associated resistance is key to devising intervention strategies. We used high-resolution genomic analysis of a large global collection of ST131 isolates to define the evolutionary history of extended-spectrum beta-lactamase production in ST131. We documented diverse contributory genetic processes, including stable chromosomal integrations of resistance genes, persistence and evolution of mobile resistance elements within sublineages, and sporadic acquisition of different resistance elements. Both global distribution and regional segregation were evident. The diversity of resistance element acquisition and propagation within ST131 indicates a need for control and surveillance strategies that target both bacterial strains and mobile genetic elements.

Detection of homologous recombination in closely related strains

Detection of recombination events in a bacterial genome is both important from the evolutionary point of view, and of practical interest. Indeed, homologous recombination (HR) plays a major role in the exchange of antigenic determinants between strains. There exist statistical methods to detect recently recombined segments in whole-genome sequences that use a high local density of substitutions as a signal of HR events with a source outside considered strains. However, it is difficult to detect the HR events within a set of strains, which represent whole species diversity, due to a low number of substitutions in recombined segments and high level of diversity of strains. Here, we analyzed HR in 20 Escherichia coli (E. coli) strains to define what fraction of segments with a high substitution rate were introduced in a genome by HR. For detection of HR, we used the segmentation, performed by the adaptive weights smoothing (AWS) algorithm. It detects sharp changes in the structure of observed data analyzing only qualitative structural information. We validated the approach on simulated data, applied it to the analysis of E. coli strains, and determined the recombination rates between phylogroups.

The Impact of Selection, Gene Conversion, and Biased Sampling on the Assessment of Microbial Demography

Recent studies have linked demographic changes and epidemiological patterns in bacterial populations using coalescent-based approaches. We identified 26 studies using skyline plots and found that 21 inferred overall population expansion. This surprising result led us to analyze the impact of natural selection, recombination (gene conversion), and sampling biases on demographic inference using skyline plots and site frequency spectra (SFS). Forward simulations based on biologically relevant parameters from Escherichia coli populations showed that theoretical arguments on the detrimental impact of recombination and especially natural selection on the reconstructed genealogies cannot be ignored in practice. In fact, both processes systematically lead to spurious interpretations of population expansion in skyline plots (and in SFS for selection). Weak purifying selection, and especially positive selection, had important effects on skyline plots, showing patterns akin to those of population expansions. State-of-the-art techniques to remove recombination further amplified these biases. We simulated three common sampling biases in microbiological research: uniform, clustered, and mixed sampling. Alone, or together with recombination and selection, they further mislead demographic inferences producing almost any possible skyline shape or SFS. Interestingly, sampling sub-populations also affected skyline plots and SFS, because the coalescent rates of populations and their sub-populations had different distributions. This study suggests that extreme caution is needed to infer demographic changes solely based on reconstructed genealogies. We suggest that the development of novel sampling strategies and the joint analyzes of diverse population genetic methods are strictly necessary to estimate demographic changes in populations where selection, recombination, and biased sampling are present.

The Landscape of Realized Homologous Recombination in Pathogenic Bacteria

Recombination enhances the adaptive potential of organisms by allowing genetic variants to be tested on multiple genomic backgrounds. Its distribution in the genome can provide insight into the evolutionary forces that underlie traits, such as the emergence of pathogenicity. Here, we examined landscapes of realized homologous recombination of 500 genomes from ten bacterial species and found all species have "hot" regions with elevated rates relative to the genome average. We examined the size, gene content, and chromosomal features associated with these regions and the correlations between closely related species. The recombination landscape is variable and evolves rapidly. For example in Salmonella, only short regions of around 1 kb in length are hot whereas in the closely related species Escherichia coli, some hot regions exceed 100 kb, spanning many genes. Only Streptococcus pyogenes shows evidence for the positive correlation between GC content and recombination that has been reported for several eukaryotes. Genes with function related to the cell surface/membrane are often found in recombination hot regions but E. coli is the only species where genes annotated as "virulence associated" are consistently hotter. There is also evidence that some genes with "housekeeping" functions tend to be overrepresented in cold regions. For example, ribosomal proteins showed low recombination in all of the species. Among specific genes, transferrin-binding proteins are recombination hot in all three of the species in which they were found, and are subject to interspecies recombination.

Within-host evolution of bacterial pathogens

Whole-genome sequencing has opened the way for investigating the dynamics and genomic evolution of bacterial pathogens during the colonization and infection of humans. The application of this technology to the longitudinal study of adaptation in an infected host--in particular, the evolution of drug resistance and host adaptation in patients who are chronically infected with opportunistic pathogens--has revealed remarkable patterns of convergent evolution, suggestive of an inherent repeatability of evolution. In this Review, we describe how these studies have advanced our understanding of the mechanisms and principles of within-host genome evolution, and we consider the consequences of findings such as a potent adaptive potential for pathogenicity. Finally, we discuss the possibility that genomics may be used in the future to predict the clinical progression of bacterial infections and to suggest the best option for treatment.

Evolution of atypical enteropathogenic E. coli by repeated acquisition of LEE pathogenicity island variants

Atypical enteropathogenic Escherichia coli (aEPEC) is an umbrella term given to E. coli that possess a type III secretion system encoded in the locus of enterocyte effacement (LEE), but lack the virulence factors (stx, bfpA) that characterize enterohaemorrhagic E. coli and typical EPEC, respectively. The burden of disease caused by aEPEC has recently increased in industrialized and developing nations, yet the population structure and virulence profile of this emerging pathogen are poorly understood. Here, we generated whole-genome sequences of 185 aEPEC isolates collected during the Global Enteric Multicenter Study from seven study sites in Asia and Africa, and compared them with publicly available E. coli genomes. Phylogenomic analysis revealed ten distinct widely distributed aEPEC clones. Analysis of genetic variation in the LEE pathogenicity island identified 30 distinct LEE subtypes divided into three major lineages. Each LEE lineage demonstrated a preferred chromosomal insertion site and different complements of non-LEE encoded effector genes, indicating distinct patterns of evolution of these lineages. This study provides the first detailed genomic framework for aEPEC in the context of the EPEC pathotype and will facilitate further studies into the epidemiology and pathogenicity of EPEC by enabling the detection and tracking of specific clones and LEE variants.

Weakly Deleterious Mutations and Low Rates of Recombination Limit the Impact of Natural Selection on Bacterial Genomes

Free-living bacteria are usually thought to have large effective population sizes, and so tiny selective differences can drive their evolution. However, because recombination is infrequent, “background selection” against slightly deleterious alleles should reduce the effective population size (N_e) by orders of magnitude. For example, for a well-mixed population with 10¹² individuals and a typical level of homologous recombination (r/m = 3, i.e., nucleotide changes due to recombination [r] occur at 3 times the mutation rate [m]), we predict that N_e is <10⁷. An argument for high N_e values for bacteria has been the high genetic diversity within many bacterial “species,” but this diversity may be due to population structure: diversity across subpopulations can be far higher than diversity within a subpopulation, which makes it difficult to estimate N_e correctly. Given an estimate of N_e, standard population genetics models imply that selection should be sufficient to drive evolution if N_e × s is >1, where s is the selection coefficient. We found that this remains approximately correct if background selection is occurring or when population structure is present. Overall, we predict that even for free-living bacteria with enormous populations, natural selection is only a significant force if s is above 10⁻⁷ or so.

Because bacteria form huge populations with trillions of individuals, the simplest theoretical prediction is that the better allele at a site would predominate even if its advantage was just 10⁻⁹ per generation. In other words, virtually every nucleotide would be at the local optimum in most individuals. A more sophisticated theory considers that bacterial genomes have millions of sites each and selection events on these many sites could interfere with each other, so that only larger effects would be important. However, bacteria can exchange genetic material, and in principle, this exchange could eliminate the interference between the evolution of the sites. We used simulations to confirm that during multisite evolution with realistic levels of recombination, only larger effects are important. We propose that advantages of less than 10⁻⁷ are effectively neutral.

Comparison of Xenorhabdus bovienii bacterial strain genomes reveals diversity in symbiotic functions

BACKGROUND

Xenorhabdus bacteria engage in a beneficial symbiosis with Steinernema nematodes, in part by providing activities that help kill and degrade insect hosts for nutrition. Xenorhabdus strains (members of a single species) can display wide variation in host-interaction phenotypes and genetic potential indicating that strains may differ in their encoded symbiosis factors, including secreted metabolites.

METHODS

To discern strain-level variation among symbiosis factors, and facilitate the identification of novel compounds, we performed a comparative analysis of the genomes of 10 Xenorhabdus bovienii bacterial strains.

RESULTS

The analyzed X. bovienii draft genomes are broadly similar in structure (e.g. size, GC content, number of coding sequences). Genome content analysis revealed that general classes of putative host-microbe interaction functions, such as secretion systems and toxin classes, were identified in all bacterial strains. In contrast, we observed diversity of individual genes within families (e.g. non-ribosomal peptide synthetase clusters and insecticidal toxin components), indicating the specific molecules secreted by each strain can vary. Additionally, phenotypic analysis indicates that regulation of activities (e.g. enzymes and motility) differs among strains.

CONCLUSIONS

The analyses presented here demonstrate that while general mechanisms by which X. bovienii bacterial strains interact with their invertebrate hosts are similar, the specific molecules mediating these interactions differ. Our data support that adaptation of individual bacterial strains to distinct hosts or niches has occurred. For example, diverse metabolic profiles among bacterial symbionts may have been selected by dissimilarities in nutritional requirements of their different nematode hosts. Similarly, factors involved in parasitism (e.g. immune suppression and microbial competition factors), likely differ based on evolution in response to naturally encountered organisms, such as insect hosts, competitors, predators or pathogens. This study provides insight into effectors of a symbiotic lifestyle, and also highlights that when mining Xenorhabdus species for novel natural products, including antibiotics and insecticidal toxins, analysis of multiple bacterial strains likely will increase the potential for the discovery of novel molecules.

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.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

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.

Microbial diversity. Fine-scale diversity and extensive recombination in a quasisexual bacterial population occupying a broad niche

Extensive fine-scale genetic diversity is found in many microbial species across varied environments, but for most, the evolutionary scenarios that generate the observed variation remain unclear. Deep sequencing of a thermophilic cyanobacterial population and analysis of the statistics of synonymous single-nucleotide polymorphisms revealed a high rate of homologous recombination and departures from neutral drift consistent with the effects of genetic hitchhiking. A sequenced isolate genome resembled an unlinked random mixture of the allelic diversity at the sampled loci. These observations suggested a quasisexual microbial population that occupies a broad ecological niche, with selection driving frequencies of alleles rather than whole genomes.

Tackling Drug Resistant Infection Outbreaks of Global Pandemic Escherichia coli ST131 Using Evolutionary and Epidemiological Genomics

High-throughput molecular screening is required to investigate the origin and diffusion of antimicrobial resistance in pathogen outbreaks. The most frequent cause of human infection is Escherichia coli, which is dominated by sequence type 131 (ST131)-a set of rapidly radiating pandemic clones. The highly infectious clades of ST131 originated firstly by a mutation enhancing conjugation and adhesion. Secondly, single-nucleotide polymorphisms occurred enabling fluoroquinolone-resistance, which is near-fixed in all ST131. Thirdly, broader resistance through beta-lactamases has been gained and lost frequently, symptomatic of conflicting environmental selective effects. This flexible approach to gene exchange is worrying and supports the proposition that ST131 will develop an even wider range of plasmid and chromosomal elements promoting antimicrobial resistance. To stop ST131, deep genome sequencing is required to understand the origin, evolution and spread of antimicrobial resistance genes. Phylogenetic methods that decipher past events can predict future patterns of virulence and transmission based on genetic signatures of adaptation and gene exchange. Both the effect of partial antimicrobial exposure and cell dormancy caused by variation in gene expression may accelerate the development of resistance. High-throughput sequencing can decode measurable evolution of cell populations within patients associated with systems-wide changes in gene expression during treatments. A multi-faceted approach can enhance assessment of antimicrobial resistance in E. coli ST131 by examining transmission dynamics between hosts to achieve a goal of pre-empting resistance before it emerges by optimising antimicrobial treatment protocols.

Epidemic Clones, Oceanic Gene Pools, and Eco-LD in the Free Living Marine Pathogen Vibrio parahaemolyticus

We investigated global patterns of variation in 157 whole-genome sequences of Vibrio parahaemolyticus, a free-living and seafood associated marine bacterium. Pandemic clones, responsible for recent outbreaks of gastroenteritis in humans, have spread globally. However, there are oceanic gene pools, one located in the oceans surrounding Asia and another in the Mexican Gulf. Frequent recombination means that most isolates have acquired the genetic profile of their current location. We investigated the genetic structure in the Asian gene pool by calculating the effective population size in two different ways. Under standard neutral models, the two estimates should give similar answers but we found a 27-fold difference. We propose that this discrepancy is caused by the subdivision of the species into a hundred or more ecotypes which are maintained stably in the population. To investigate the genetic factors involved, we used 51 unrelated isolates to conduct a genome-wide scan for epistatically interacting loci. We found a single example of strong epistasis between distant genome regions. A majority of strains had a type VI secretion system associated with bacterial killing. The remaining strains had genes associated with biofilm formation and regulated by cyclic dimeric GMP signaling. All strains had one or other of the two systems and none of isolate had complete complements of both systems, although several strains had remnants. Further "top down" analysis of patterns of linkage disequilibrium within frequently recombining species will allow a detailed understanding of how selection acts to structure the pattern of variation within natural bacterial populations.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Whole-genome sequencing reveals clonal expansion of multiresistant Staphylococcus haemolyticus in European hospitals

OBJECTIVES

Staphylococcus haemolyticus is an emerging cause of nosocomial infections, primarily affecting immunocompromised patients. A comparative genomic analysis was performed on clinical S. haemolyticus isolates to investigate their genetic relationship and explore the coding sequences with respect to antimicrobial resistance determinants and putative hospital adaptation.

METHODS

Whole-genome sequencing was performed on 134 isolates of S. haemolyticus from geographically diverse origins (Belgium, 2; Germany, 10; Japan, 13; Norway, 54; Spain, 2; Switzerland, 43; UK, 9; USA, 1). Each genome was individually assembled. Protein coding sequences (CDSs) were predicted and homologous genes were categorized into three types: Type I, core genes, homologues present in all strains; Type II, unique core genes, homologues shared by only a subgroup of strains; and Type III, unique genes, strain-specific CDSs. The phylogenetic relationship between the isolates was built from variable sites in the form of single nucleotide polymorphisms (SNPs) in the core genome and used to construct a maximum likelihood phylogeny.

RESULTS

SNPs in the genome core regions divided the isolates into one major group of 126 isolates and one minor group of isolates with highly diverse genomes. The major group was further subdivided into seven clades (A-G), of which four (A-D) encompassed isolates only from Europe. Antimicrobial multiresistance was observed in 77.7% of the collection. High levels of homologous recombination were detected in genes involved in adherence, staphylococcal host adaptation and bacterial cell communication.

CONCLUSIONS

The presence of several successful and highly resistant clones underlines the adaptive potential of this opportunistic pathogen.

Extensive genomic variation within clonal bacterial groups resulted from homologous recombination

Due to divergence, genetic variation is generally believed to be high among distantly related strains, low among closely related ones and little or none within the same classified clonal groups. Several recent genome-wide studies, however, revealed that significant genetic variation resides in a considerable number of genes among strains with identical MLST (Multilocus sequence typing) types and much of the variation was introduced by homologous recombination. Recognizing and understanding genomic variation within clonal bacterial groups could shed new light on the evolutionary path of infectious agents and the emergence of particularly pathogenic or virulent variants. This commentary presents our recent contributions to this line of work.

An evolutionary analysis of genome expansion and pathogenicity in Escherichia coli

BACKGROUND

There are several studies describing loss of genes through reductive evolution in microbes, but how selective forces are associated with genome expansion due to horizontal gene transfer (HGT) has not received similar attention. The aim of this study was therefore to examine how selective pressures influence genome expansion in 53 fully sequenced and assembled Escherichia coli strains. We also explored potential connections between genome expansion and the attainment of virulence factors. This was performed using estimations of several genomic parameters such as AT content, genomic drift (measured using relative entropy), genome size and estimated HGT size, which were subsequently compared to analogous parameters computed from the core genome consisting of 1729 genes common to the 53 E. coli strains. Moreover, we analyzed how selective pressures (quantified using relative entropy and dN/dS), acting on the E. coli core genome, influenced lineage and phylogroup formation.

RESULTS

Hierarchical clustering of dS and dN estimations from the E. coli core genome resulted in phylogenetic trees with topologies in agreement with known E. coli taxonomy and phylogroups. High values of dS, compared to dN, indicate that the E. coli core genome has been subjected to substantial purifying selection over time; significantly more than the non-core part of the genome (p<0.001). This is further supported by a linear association between strain-wise dS and dN values (β = 26.94 ± 0.44, R2~0.98, p<0.001). The non-core part of the genome was also significantly more AT-rich (p<0.001) than the core genome and E. coli genome size correlated with estimated HGT size (p<0.001). In addition, genome size (p<0.001), AT content (p<0.001) as well as estimated HGT size (p<0.005) were all associated with the presence of virulence factors, suggesting that pathogenicity traits in E. coli are largely attained through HGT. No associations were found between selective pressures operating on the E. coli core genome, as estimated using relative entropy, and genome size (p~0.98).

CONCLUSIONS

On a larger time frame, genome expansion in E. coli, which is significantly associated with the acquisition of virulence factors, appears to be independent of selective forces operating on the core genome.

Within-population distribution of trimethoprim resistance in Escherichia coli before and after a community-wide intervention on trimethoprim use

A 2-year prospective intervention on the prescription of trimethoprim reduced the use by 85% in a health care region with 178,000 inhabitants. Here, we performed before-and-after analyses of the within-population distribution of trimethoprim resistance in Escherichia coli. Phylogenetic and population genetic methods were applied to multilocus sequence typing data of 548 consecutively collected E. coli isolates from clinical urinary specimens. Results were analyzed in relation to antibiotic susceptibility and the presence and genomic location of different trimethoprim resistance gene classes. A total of 163 E. coli sequence types (STs) were identified, of which 68 were previously undescribed. The isolates fell into one of three distinct genetic clusters designated BAPS 1 (E. coli phylogroup B2), BAPS 2 (phylogroup A and B1), and BAPS 3 (phylogroup D), each with a similar frequency before and after the intervention. BAPS 2 and BAPS 3 were positively and BAPS 1 was negatively associated with trimethoprim resistance (odds ratios of 1.97, 3.17, and 0.26, respectively). In before-and-after analyses, trimethoprim resistance frequency increased in BAPS 1 and decreased in BAPS 2. Resistance to antibiotics other than trimethoprim increased in BAPS 2. Analysis of the genomic location of different trimethoprim resistance genes in isolates of ST69, ST58, and ST73 identified multiple independent acquisition events in isolates of the same ST. The results show that despite a stable overall resistance frequency in E. coli before and after the intervention, marked within-population changes occurred. A decrease of resistance in one major genetic cluster was masked by a reciprocal increase in another major cluster.

Recombination blurs phylogenetic groups routine assignment in Escherichia coli: setting the record straight

The characterization of population structures plays a main role for understanding outbreaks and the dynamics of bacterial spreading. In Escherichia coli, the widely used combination of multiplex-PCR scheme together with goeBURST has some limitations. The purpose of this study is to show that the combination of different phylogenetic approaches based on concatenated sequences of MLST genes results in a more precise assignment of E. coli phylogenetic groups, complete understanding of population structure and reconstruction of ancestral clones. A collection of 80 Escherichia coli strains of different origins was analyzed following the Clermont and Doumith's multiplex-PCR schemes. Doumith's multiplex-PCR showed only 1.7% of misassignment, whereas Clermont's-2000 protocol reached 14.0%, although the discrepancies reached 30% and 38.7% respectively when recombinant C, F and E phylogroups were considered. Therefore, correct phylogroup attribution is highly variable and depends on the clonal composition of the sample. As far as population structure of these E. coli strains, including 48 E. coli genomes from GenBank, goeBURST provides a quite dispersed population structure; whereas NeighborNet approach reveals a complex population structure. MLST-based eBURST can infer different founder genotypes, for instance ST23/ST88 could be detected as the founder genotypes for STC23; however, phylogenetic reconstructions might suggest ST410 as the ancestor clone and several evolutionary trajectories with different founders. To improve our routine understanding of E. coli molecular epidemiology, we propose a strategy based on three successive steps; first, to discriminate three main groups A/B1/C, D/F/E and B2 following Doumith's protocol; second, visualization of population structure based on MLST genes according to goeBURST, using NeighborNet to establish more complex relationships among STs; and third, to perform, a cost-free characterization of evolutionary trajectories in variants emerging along the clonal expansion using parsimony methods of phylogenetic analysis.

Genomic analysis of the emergence of 20th century epidemic dysentery

Background

Shigella dysenteriae type 1 (Sd1) causes recurrent epidemics of dysentery associated with high mortality in many regions of the world. Sd1 infects humans at very low infectious doses (10 CFU), and treatment is complicated by the rapid emergence of antibiotic resistant Sd1 strains. Sd1 is only detected in the context of human infections, and the circumstances under which epidemics emerge and regress remain unknown.

Results

Phylogenomic analyses of 56 isolates collected worldwide over the past 60 years indicate that the Sd1 clone responsible for the recent pandemics emerged at the turn of the 20^th century, and that the two world wars likely played a pivotal role for its dissemination. Several lineages remain ubiquitous and their phylogeny indicates several recent intercontinental transfers. Our comparative genomics analysis reveals that isolates responsible for separate outbreaks, though closely related to one another, have independently accumulated antibiotic resistance genes, suggesting that there is little or no selection to retain these genes in-between outbreaks. The genomes appear to be subjected to genetic drift that affects a number of functions currently used by diagnostic tools to identify Sd1, which could lead to the potential failure of such tools.

Conclusions

Taken together, the Sd1 population structure and pattern of evolution suggest a recent emergence and a possible human carrier state that could play an important role in the epidemic pattern of infections of this human-specific pathogen. This analysis highlights the important role of whole-genome sequencing in studying pathogens for which epidemiological or laboratory investigations are particularly challenging.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-355) contains supplementary material, which is available to authorized users.

Global dissemination of a multidrug resistant Escherichia coli clone

Escherichia coli sequence type 131 (ST131) is a globally disseminated, multidrug resistant (MDR) clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with several factors, including resistance to fluoroquinolones, high virulence gene content, the possession of the type 1 fimbriae FimH30 allele, and the production of the CTX-M-15 extended spectrum β-lactamase (ESBL). Here, we used genome sequencing to examine the molecular epidemiology of a collection of E. coli ST131 strains isolated from six distinct geographical locations across the world spanning 2000-2011. The global phylogeny of E. coli ST131, determined from whole-genome sequence data, revealed a single lineage of E. coli ST131 distinct from other extraintestinal E. coli strains within the B2 phylogroup. Three closely related E. coli ST131 sublineages were identified, with little association to geographic origin. The majority of single-nucleotide variants associated with each of the sublineages were due to recombination in regions adjacent to mobile genetic elements (MGEs). The most prevalent sublineage of ST131 strains was characterized by fluoroquinolone resistance, and a distinct virulence factor and MGE profile. Four different variants of the CTX-M ESBL-resistance gene were identified in our ST131 strains, with acquisition of CTX-M-15 representing a defining feature of a discrete but geographically dispersed ST131 sublineage. This study confirms the global dispersal of a single E. coli ST131 clone and demonstrates the role of MGEs and recombination in the evolution of this important MDR pathogen.

Efficient inference of recombination hot regions in bacterial genomes

In eukaryotes, detailed surveys of recombination rates have shown variation at multiple genomic scales and the presence of “hotspots” of highly elevated recombination. In bacteria, studies of recombination rate variation are less developed, in part because there are few analysis methods that take into account the clonal context within which bacterial evolution occurs. Here, we focus in particular on identifying “hot regions” of the genome where DNA is transferred frequently between isolates. We present a computationally efficient algorithm based on the recently developed “chromosome painting” algorithm, which characterizes patterns of haplotype sharing across a genome. We compare the average genome wide painting, which principally reflects clonal descent, with the painting for each site which additionally reflects the specific deviations at the site due to recombination. Using simulated data, we show that hot regions have consistently higher deviations from the genome wide average than normal regions. We applied our approach to previously analyzed Escherichia coli genomes and revealed that the new method is highly correlated with the number of recombination events affecting each site inferred by ClonalOrigin, a method that is only applicable to small numbers of genomes. Furthermore, we analyzed recombination hot regions in Campylobacter jejuni by using 200 genomes. We identified three recombination hot regions, which are enriched for genes related to membrane proteins. Our approach and its implementation, which is downloadable from https://github.com/bioprojects/orderedPainting, will help to develop a new phase of population genomic studies of recombination in prokaryotes.

Shiga toxin-producing Escherichia coli

In the United States, it is estimated that non-O157 Shiga toxin-producing Escherichia coli (STEC) cause more illnesses than STEC O157:H7, and the majority of cases of non-O157 STEC infections are due to serogroups O26, O45, O103, O111, O121, and O145, referred to as the top six non-O157 STEC. The diseases caused by non-O157 STEC are generally milder than those induced by O157 STEC; nonetheless, non-O157 STEC strains have also been associated with serious illnesses such as hemorrhagic colitis and hemolytic uremic syndrome, as well as death. Ruminants, particularly cattle, are reservoirs for both O157 and non-O157 STEC, which are transmitted to humans by person-to-person or animal contact and by ingestion of food or water contaminated with animal feces. Improved strategies to control STEC colonization and shedding in cattle and contamination of meat and produce are needed. In general, non-O157 STEC respond to stresses such as acid, heat, and other stresses induced during food preparation similar to O157 STEC. Similar to O157:H7, the top six non-O157 STEC are classified as adulterants in beef by the USDA Food Safety and Inspection Service, and regulatory testing for these pathogens began in June 2012. Due to the genetic and phenotypic variability of non-O157 STEC strains, the development of accurate and reliable methods for detection and isolation of these pathogens has been challenging. Since the non-O157 STEC are responsible for a large portion of STEC-related illnesses, more extensive studies on their physiology, genetics, pathogenicity, and evolution are needed in order to develop more effective control strategies.

Inference of the properties of the recombination process from whole bacterial genomes

Patterns of linkage disequilibrium, homoplasy, and incompatibility are difficult to interpret because they depend on several factors, including the recombination process and the population structure. Here we introduce a novel model-based framework to infer recombination properties from such summary statistics in bacterial genomes. The underlying model is sequentially Markovian so that data can be simulated very efficiently, and we use approximate Bayesian computation techniques to infer parameters. As this does not require us to calculate the likelihood function, the model can be easily extended to investigate less probed aspects of recombination. In particular, we extend our model to account for the bias in the recombination process whereby closely related bacteria recombine more often with one another. We show that this model provides a good fit to a data set of Bacillus cereus genomes and estimate several recombination properties, including the rate of bias in recombination. All the methods described in this article are implemented in a software package that is freely available for download at http://code.google.com/p/clonalorigin/.

Prevalence of autotransporters in Escherichia coli: what is the impact of phylogeny and pathotype?

Autotransporter (AT) proteins are widespread surface-exposed or secreted factors in Escherichia coli. Several ATs have been correlated with pathogenesis or specific phylogenetic lineages. Therefore, an application as biomarkers for individual extraintestinal pathogenic E.coli (ExPEC) or intestinal pathogenic E.coli (IPEC) has been proposed. To put this assumption on a solid foundation, we analyzed 111 publicly available E. coli genome sequences and screened them bioinformatically for the presence of 18 ATs. We determined the highest AT prevalence per strain in phylogroup B2 isolates and showed that AT distribution correlates rather with phylogenetic lineages than with pathotypes. Although a strict dependence between AT prevalence and pathotype was not observed, EspP, EhaA, and EhaG cluster with IPEC of phylogroup B1 and E, respectively, whereas UpaH is predominantly present in ExPEC of phylogroup B2. Furthermore, PicU, SepA, UpaB, UpaI, and UpaJ were associated with phylogroup B2. We detected UpaI and its positional ortholog EhaC in 93% of the E.coli strains tested. This AT variant is thus the most prevalent in E.coli irrespective of pathotype or phylogenetic background. Compared with the ATs UpaB, UpaC, and UpaJ of uropathogenic E.coli strain 536, UpaI had redundant functions, contributing to autoaggregation, biofilm formation, and binding to extracellular matrix proteins. The functional redundancy and wide distribution of ATs among pathogenic and non-pathogenic E.coli indicates that ATs cannot generally be regarded as specific biomarkers and virulence factors. Our results demonstrate that phylogeny has a bigger impact on the distribution of AT variants in E.coli than initially thought, especially in ExPEC.

The evolutionary path to extraintestinal pathogenic, drug-resistant Escherichia coli is marked by drastic reduction in detectable recombination within the core genome

Escherichia coli is a highly diverse group of pathogens ranging from commensal of the intestinal tract, through to intestinal pathogen, and extraintestinal pathogen. Here, we present data on the population diversity of E. coli, using Bayesian analysis to identify 13 distinct clusters within the population from multilocus sequence typing data, which map onto a whole-genome-derived phylogeny based on 62 genome sequences. Bayesian analysis of recombination within the core genome identified reduction in detectable core genome recombination as one moves from the commensals, through the intestinal pathogens down to the multidrug-resistant extraintestinal pathogenic clone E. coli ST131. Our data show that the emergence of a multidrug-resistant, extraintestinal pathogenic lineage of E. coli is marked by substantial reduction in detectable core genome recombination, resulting in a lineage which is phylogenetically distinct and sexually isolated in terms of core genome recombination.

Genomic evolution and transmission of Helicobacter pylori in two South African families

Helicobacter pylori infects the stomachs of one in two humans and can cause sequelae that include ulcers and cancer. Here we sequenced the genomes of 97 H. pylori isolates from 52 members of two families living in rural conditions in South Africa. From each of 45 individuals, two H. pylori strains were isolated from the antrum and corpus parts of the stomach, and comparisons of their genomes enabled us to study within-host evolution. In 5 of these 45 hosts, the two genomes were too distantly related to be derived from each other and therefore represented evidence of multiple infections. From the remaining 40 genome pairs, we estimated that the synonymous mutation rate was 1.38 × 10(-5) per site per year, with a low effective population size within host probably reflecting population bottlenecks and immune selection. Some individuals showed very little evidence for recombination, whereas in others, recombination introduced up to 100-times more substitutions than mutation. These differences may reflect unequal opportunities for recombination depending on the presence or absence of multiple infections. Comparing the genomes carried by distinct individuals enabled us to establish probable transmission links. Transmission events were found significantly more frequently between close relatives, and between individuals living in the same house. We found, however, that a majority of individuals (27/52) were not linked by transmission to other individuals. Our results suggest that transmission does not always occur within families, and that coinfection with multiple strains is frequent and evolutionarily important despite a fast turnover of the infecting strains within-host.

Chromosome painting in silico in a bacterial species reveals fine population structure

Identifying population structure forms an important basis for genetic and evolutionary studies. Most current methods to identify population structure have limitations in analyzing haplotypes and recombination across the genome. Recently, a method of chromosome painting in silico has been developed to overcome these shortcomings and has been applied to multiple human genome sequences. This method detects the genome-wide transfer of DNA sequence chunks through homologous recombination. Here, we apply it to the frequently recombining bacterial species Helicobacter pylori that has infected Homo sapiens since their birth in Africa and shows wide phylogeographic divergence. Multiple complete genome sequences were analyzed including sequences from Okinawa, Japan, that we recently sequenced. The newer method revealed a finer population structure than revealed by a previous method that examines only MLST housekeeping genes or a phylogenetic network analysis of the core genome. Novel subgroups were found in Europe, Amerind, and East Asia groups. Examination of genetic flux showed some singleton strains to be hybrids of subgroups and revealed evident signs of population admixture in Africa, Europe, and parts of Asia. We expect this approach to further our understanding of intraspecific bacterial evolution by revealing population structure at a finer scale.

Comparisons of infant Escherichia coli isolates link genomic profiles with adaptation to the ecological niche

BACKGROUND

Despite being one of the most intensely studied model organisms, many questions still remain about the evolutionary biology and ecology of Escherichia coli. An important step toward achieving a more complete understanding of E.coli biology entails elucidating relationships between gene content and adaptation to the ecological niche.

RESULTS

Here, we present genome comparisons of 16 E.coli strains that represent commensals and pathogens isolated from infants during a specific time period in Trondheim, Norway. Using differential gene content, we characterized enrichment profiles of the collection of strains relating to phylogeny, early vs. late colonization, pathogenicity and growth rate. We found clear gene content distinctions relating to the various grouping criteria. We also found that different categories of strains use different genetic elements for similar biological processes. The sequenced genomes included two pairs of strains where each pair was isolated from the same infant at different time points. One pair, in which the strains were isolated four months apart, showed maintenance of an early colonizer genome profile but also gene content and codon usage changes toward the late colonizer profile. Lastly, we placed our sequenced isolates into a broader genomic context by comparing them with 25 published E.coli genomes that represent a variety of pathotypes and commensal strains. This analysis demonstrated the importance of geography in shaping strain level gene content profiles.

CONCLUSIONS

Our results indicate a general pattern where alternative genetic pathways lead toward a consistent ecological role for E.coli as a species. Within this framework however, we saw selection shaping the coding repertoire of E.coli strains toward distinct ecotypes with different phenotypic properties.

Homologous recombination drives both sequence diversity and gene content variation in Neisseria meningitidis

The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis (Nm) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm. Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm. These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.

Inference of reticulate evolutionary histories by maximum likelihood: the performance of information criteria

BACKGROUND

Maximum likelihood has been widely used for over three decades to infer phylogenetic trees from molecular data. When reticulate evolutionary events occur, several genomic regions may have conflicting evolutionary histories, and a phylogenetic network may provide a more adequate model for representing the evolutionary history of the genomes or species. A maximum likelihood (ML) model has been proposed for this case and accounts for both mutation within a genomic region and reticulation across the regions. However, the performance of this model in terms of inferring information about reticulate evolution and properties that affect this performance have not been studied.

RESULTS

In this paper, we study the effect of the evolutionary diameter and height of a reticulation event on its identifiability under ML. We find both of them, particularly the diameter, have a significant effect. Further, we find that the number of genes (which can be generalized to the concept of "non-recombining genomic regions") that are transferred across a reticulation edge affects its detectability. Last but not least, a fundamental challenge with phylogenetic networks is that they allow an arbitrary level of complexity, giving rise to the model selection problem. We investigate the performance of two information criteria, the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), for addressing this problem. We find that BIC performs well in general for controlling the model complexity and preventing ML from grossly overestimating the number of reticulation events.

CONCLUSION

Our results demonstrate that BIC provides a good framework for inferring reticulate evolutionary histories. Nevertheless, the results call for caution when interpreting the accuracy of the inference particularly for data sets with particular evolutionary features.

Hierarchical clustering of genetic diversity associated to different levels of mutation and recombination in Escherichia coli: a study based on Mexican isolates

Escherichia coli occur as either free-living microorganisms, or within the colons of mammals and birds as pathogenic or commensal bacteria. Although the Mexican population of intestinal E. coli maintains high levels of genetic diversity, the exact mechanisms by which this occurs remain unknown. We therefore investigated the role of homologous recombination and point mutation in the genetic diversification and population structure of Mexican strains of E. coli. This was explored using a multi locus sequence typing (MLST) approach in a non-outbreak related, host-wide sample of 128 isolates. Overall, genetic diversification in this sample appears to be driven primarily by homologous recombination, and to a lesser extent, by point mutation. Since genetic diversity is hierarchically organized according to the MLST genealogy, we observed that there is not a homogeneous recombination rate, but that different rates emerge at different clustering levels such as phylogenetic group, lineage and clonal complex (CC). Moreover, we detected clear signature of substructure among the A+B1 phylogenetic group, where the majority of isolates were differentiated into four discrete lineages. Substructure pattern is revealed by the presence of several CCs associated to a particular life style and host as well as to different genetic diversification mechanisms. We propose these findings as an alternative explanation for the maintenance of the clear phylogenetic signal of this species despite the prevalence of homologous recombination. Finally, we corroborate using both phylogenetic and genetic population approaches as an effective mean to establish epidemiological surveillance tailored to the ecological specificities of each geographic region.