Haplotype diversity in "source-sink" dynamics of Escherichia coli urovirulence

Portrait of a generalist bacterium: pathoadaptation, metabolic specialization and extreme environments shape diversity of Staphylococcus saprophyticus

Staphylococcus saprophyticus is a Gram-positive, coagulase-negative staphylococcus found in diverse environments including soil and freshwater, meat, and dairy foods. S. saprophyticus is also an important cause of urinary tract infections (UTIs) in humans, and mastitis in cattle. However, the genetic determinants of virulence have not yet been identified, and it remains unclear whether there are distinct sub-populations adapted to human and animal hosts. Using a diverse sample of S. saprophyticus isolates from food, animals, environmental sources, and human infections, we characterized the population structure and diversity of global populations of S. saprophyticus . We found that divergence of the two major clades of S. saprophyticus is likely facilitated by barriers to horizontal gene transfer (HGT) and differences in metabolism. Using genome-wide association study (GWAS) tools we identified the first Type VII secretion system (T7SS) described in S. saprophyticus and its association with bovine mastitis. Finally, we found that in general, strains of S. saprophyticus from different niches are genetically similar with the exception of built environments, which function as a 'sink' for S. saprophyticus populations. This work increases our understanding of the ecology of S. saprophyticus and of the genomics of bacterial generalists.

Data summary

Raw sequencing data for newly sequenced S. saprophyticus isolates have been deposited to the NCBI SRA under the project accession PRJNA928770. A list of all genomes used in this work and their associated metadata are available in the supplementary material. Custom scripts used in the comparative genomics and GWAS analyses are available here: https://github.com/myoungblom/sapro_genomics .

Impact statement

It is not known whether human and cattle diseases caused by S. saprophyticus represent spillover events from a generalist adapted to survive in a range of environments, or whether the capacity to cause disease represents a specific adaptation. Seasonal cycles of S. saprophyticus UTIs and molecular epidemiological evidence suggest that these infections may be environmentally-acquired rather than via transmission from person to person. Using comparative genomics and genome wide association study tools, we found that S. saprophyticus appears adapted to inhabit a wide range of environments (generalist), with isolates from animals, food, natural environments and human infections being closely related. Bacteria that routinely switch environments, particularly between humans and animals, are of particular concern when it comes to the spread of antibiotic resistance from farm environments into human populations. This work provides a framework for comparative genomic analyses of bacterial generalists and furthers our understanding of how bacterial populations move between humans, animals, and the environment.

Persisting uropathogenic Escherichia coli lineages show signatures of niche-specific within-host adaptation mediated by mobile genetic elements

Large-scale genomic studies have identified within-host adaptation as a hallmark of bacterial infections. However, the impact of physiological, metabolic, and immunological differences between distinct niches on the pathoadaptation of opportunistic pathogens remains elusive. Here, we profile the within-host adaptation and evolutionary trajectories of 976 isolates representing 119 lineages of uropathogenic Escherichia coli (UPEC) sampled longitudinally from both the gastrointestinal and urinary tracts of 123 patients with urinary tract infections. We show that lineages persisting in both niches within a patient exhibit increased allelic diversity. Habitat-specific selection results in niche-specific adaptive mutations and genes, putatively mediating fitness in either environment. Within-lineage inter-habitat genomic plasticity mediated by mobile genetic elements (MGEs) provides the opportunistic pathogen with a mechanism to adapt to the physiological conditions of either habitat, and reduced MGE richness is associated with recurrence in gut-adapted UPEC lineages. Collectively, our results establish niche-specific adaptation as a driver of UPEC within-host evolution.

Recombinant FimH Adhesin Demonstrates How the Allosteric Catch Bond Mechanism Can Support Fast and Strong Bacterial Attachment in the Absence of Shear

The FimH protein of Escherichia coli is a model two-domain adhesin that is able to mediate an allosteric catch bond mechanism of bacterial cell attachment, where the mannose-binding lectin domain switches from an 'inactive' conformation with fast binding to mannose to an 'active' conformation with slow detachment from mannose. Because mechanical tensile force favors separation of the domains and, thus, FimH activation, it has been thought that the catch bonds can only be manifested in a fluidic shear-dependent mode of adhesion. Here, we used recombinant FimH variants with a weakened inter-domain interaction and show that a fast and sustained allosteric activation of FimH can also occur under static, non-shear conditions. Moreover, it appears that lectin domain conformational activation happens intrinsically at a constant rate, independently from its ability to interact with the pilin domain or mannose. However, the latter two factors control the rate of FimH deactivation. Thus, the allosteric catch bond mechanism can be a much broader phenomenon involved in both fast and strong cell-pathogen attachments under a broad range of hydrodynamic conditions. This concept that allostery can enable more effective receptor-ligand interactions is fundamentally different from the conventional wisdom that allostery provides a mechanism to turn binding off under specific conditions.

Adaptation to host-specific bacterial pathogens drive rapid evolution of novel PhoP/PhoQ regulation pathway modulating the virulence

The presence of the PhoP-PhoQ system is usually different in various bacterial groups, suggesting that PhoP can control the expression of different genes in species. However, little is known about the evolution of the PhoP-PhoQ system among bacterial pathogens. Here, we study the evolution of PhoP and PhoQ regulation in 15 species of Enterobacteriaceae family. We have determined that the regulatory objectives adopted by PhoP and PhoQ are mainly different, due to the result of horizontal gene transfer events and even the change in the genetic content between closely related species. We have compared many possibilities tests (M1 vs. M2 and M7 with M8) to determine the positive selection. Estimating parameters at M1 and M2, with positive selection in M2 of the two proteins. The proportions of positive selection sites significant with ω = 4.53076 for PhoP and ω = 4.21041 PhQ. M8 was significant for PhoP and PhQ proteins. To further confirm the positive selection results, we used the Selecton server to confer positive selection on individual sites using the Mechanistic-Empirical Combination model, and we noticed that several sites had been identified under selection pressure during the evolution. There was a strong indication for the positive selection in bacterial genes of PhoP and PhoQ showed the results. By the use of REL and IFEL, the positive selection for PhoP was detected 14 and 11 sites respectively at different codon positions. The positively selected sites of amino acids such as Arginine, Alanine, Lysine, and Leucine are more important for the production of signals. Our results suggest that the positive selection of PhoP-PhoQ genes in host adaptation during evolution raises an intriguing possibility causes subtle variations in actions of PhoP-PhoQ and also increases the opportunities that cause modification in protein structure for the evolution of increasing pathogenicity in bacterial pathogens.

Mutational convergence acts as a major player in adaptive parallel evolution of Shigella spp

Shigella spp., emerging from multiple origins of Escherichia coli, poses a significant health threat as a causative agent of bacillary dysentery. While multiple serotypes of four different species have evolved via independent lineages, Shigella spp. are designated as a single pathotype, primarily because of their common mode of pathogenesis. Convergent horizontal transfer events have so far been attributed to the commonalities in the evolution of virulence across diverse lineages. However, the role of mutational convergence in such parallel evolution is not yet well understood. Here we have carried out a genome-wide analysis of Shigella strains from all four species to detect the core genes (i.e. the ones present in all analyzed strains) acquiring convergent mutations of evolutionarily recent origin. Simulation studies show non-neutral accumulation of these convergent mutations across species, suggesting their adaptive role in the evolution of Shigella virulence. S. dysenteriae strain 197, representing highly virulent type 1 (Sd1) clone, carries excessively high number of core genes with recent convergent mutations compared to other analyzed strains. We propose that this high frequency of adaptive convergence in S. dysenteriae strain 197 could be linked to recent re-emergence of the Sd1 clone and its increased resistance to antimicrobials.

Recombination-independent rapid convergent evolution of the gastric pathogen Helicobacter pylori

BACKGROUND

Helicobacter pylori is a human stomach pathogen, naturally-competent for DNA uptake, and prone to homologous recombination. Extensive homoplasy (i.e., phylogenetically-unlinked identical variations) observed in H. pylori genes is considered a hallmark of such recombination. However, H. pylori also exhibits a high mutation rate. The relative adaptive role of homologous recombination and mutation in species diversity is a highly-debated issue in biology. Recombination results in homoplasy. While convergent mutation can also account for homoplasy, its contribution is thought to be minor. We demonstrate here that, contrary to dogma, convergent mutation is a key contributor to Helicobacter pylori homoplasy, potentially driven by adaptive evolution of proteins.

RESULTS

Our present genome-wide analysis shows that homoplastic nonsynonymous (amino acid replacement) changes are not typically accompanied by homoplastic synonymous (silent) variations. Moreover, the majority of the codon positions with homoplastic nonsynonymous changes also contain different (i.e. non-homoplastic) nonsynonymous changes arising from mutation only. This indicates that, to a considerable extent, nonsynonymous homoplasy is due to convergent mutations. High mutation rate or limited availability of evolvable sites cannot explain this excessive convergence, as suggested by our simulation studies. Rather, the genes with convergent mutations are overrepresented in distinct functional categories, suggesting possible selective responses to conditions such as distinct micro-niches in single hosts, and to differences in host genotype, physiology, habitat and diet.

CONCLUSIONS

We propose that mutational convergence is a key player in H. pylori's adaptation and extraordinary persistence in human hosts. High frequency of mutational convergence could be due to saturation of evolvable sites capable of responding to selection pressures, while the number of mutable residues is far from saturation. We anticipate a similar scenario of mutational vs. recombinational genome dynamics or plasticity for other naturally competent microbes where strong positive selection could favor frequent convergent mutations in adaptive protein evolution.

Whole-genome comparison of urinary pathogenic Escherichia coli and faecal isolates of UTI patients and healthy controls

The faecal flora is a common reservoir for urinary tract infection (UTI), and Escherichia coli (E. coli) is frequently found in this reservoir without causing extraintestinal infection. We investigated these E. coli reservoirs by whole-genome sequencing a large collection of E. coli from healthy controls (faecal), who had never previously had UTI, and from UTI patients (faecal and urinary) sampled from the same geographical area. We compared MLST types, phylogenetic relationship, accessory genome content and FimH type between patient and control faecal isolates as well as between UTI and faecal-only isolates, respectively. Comparison of the accessory genome of UTI isolates to faecal isolates revealed 35 gene families which were significantly more prevalent in the UTI isolates compared to the faecal isolates, although none of these were unique to one of the two groups. Of these 35, 22 belonged to a genomic island and three putatively belonged to a type VI secretion system (T6SS). MLST types and SNP phylogeny indicated no clustering of the UTI or faecal E. coli from patients distinct from the control faecal isolates, although there was an overrepresentation of UTI isolates belonging to clonal lineages CC73 and CC12. One combination of mutations in FimH, N70S/S78N, was significantly associated to UTI, while phylogenetic analysis of FimH and fimH identified no signs of distinct adaptation of UTI isolates compared to faecal-only isolates not causing UTI. In summary, the results showed that (i) healthy women who had never previously had UTI carried faecal E. coli which were overall closely related to UTI and faecal isolates from UTI patients; (ii) UTI isolates do not cluster separately from faecal-only isolates based on SNP analysis; and (iii) 22 gene families of a genomic island, putative T6SS proteins as well as specific metabolism and virulence associated proteins were significantly more common in UTI isolates compared to faecal-only isolates and (iv) evolution of fimH for these isolates was not linked to the clinical source of the isolates, apart from the mutation combination N70S/S78N, which was correlated to UTI isolates of phylogroup B2. Combined, these findings illustrate that faecal and UTI isolates, as well as faecal-only and faecal-UTI isolates, are closely related and can only be distinguished, if at all, by their accessory genome.

The evolutionary logic of sepsis

The recently proposed Microbiome Mutiny Hypothesis posits that members of the human microbiome obtain information about the host individuals' health status and, when host survival is compromised, switch to an intensive exploitation strategy to maximize residual transmission. In animals and humans, sepsis is an acute systemic reaction to microbes invading the normally sterile body compartments. When induced by formerly mutualistic or neutral microbes, possibly in response to declining host health, sepsis appears to fit the 'microbiome mutiny' scenario except for its apparent failure to enhance transmission of the causative organisms. We propose that the ability of certain species of the microbiome to induce sepsis is not a fortuitous side effect of within-host replication, but rather it might, in some cases, be the result of their adaptive evolution. Whenever host health declines, inducing sepsis can be adaptive for those members of the healthy human microbiome that are capable of colonizing the future cadaver and spread by cadaver-borne transmission. We hypothesize that such microbes might exhibit switches along the 'mutualist - lethal pathogen - decomposer - mutualist again' scenario, implicating a previously unsuspected, surprising level of phenotypic plasticity. This hypothesis predicts that those species of the healthy microbiome that are recurring causative agents of sepsis can participate in the decomposition of cadavers, and can be transmitted as soil-borne or water-borne infections. Furthermore, in individual sepsis cases, the same microbial clones that dominate the systemic infection that precipitates sepsis, should also be present in high concentration during decomposition following death: this prediction is testable by molecular fingerprinting in experimentally induced animal models. Sepsis is a leading cause of human death worldwide. If further research confirms that some cases of sepsis indeed involve the 'mutiny' (facultative phenotypic switching) of normal members of the microbiome, then new strategies could be devised to prevent or treat sepsis by interfering with this process.

Pathoadaptive Mutations in Uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) are opportunistic human pathogens that primarily circulate as part of commensal intestinal microbiota. Though they have the ability to survive and proliferate in various urinary tract compartments, the urinary tract is a transient, occasional habitat for UPEC. Because of this, most of the UPEC traits have originally evolved to serve in intestinal colonization and transmission. Some of these bacterial traits serve as virulence factors - they are critical to or assist in survival of UPEC as pathogens, and the structure and/or function may be specialized for the infection. Other traits could serve as anti-virulence factors - they represent liability in the urinary tract and are under selection to be lost or inactivated during the infection. Inactivation, variation, or other changes of the bacterial genes that increase the pathogen's fitness during the infection are called pathoadaptive mutations. This chapter describes examples of pathoadaptive mutations in UPEC and provides rationale for their further in-depth study.

The adaptation rate of a quantitative trait in an environmental gradient

The spatial range of a species habitat is generally determined by the ability of the species to cope with biotic and abiotic variables that vary in space. Therefore, the species range is itself an evolvable property. Indeed, environmental gradients permit a mode of evolution in which range expansion and adaptation go hand in hand. This process can contribute to rapid evolution of drug resistant bacteria and viruses, because drug concentrations in humans and livestock treated with antibiotics are far from uniform. Here, we use a minimal stochastic model of discrete, interacting organisms evolving in continuous space to study how the rate of adaptation of a quantitative trait depends on the steepness of the gradient and various population parameters. We discuss analytical results for the mean-field limit as well as extensive stochastic simulations. These simulations were performed using an exact, event-driven simulation scheme that can deal with continuous time-, density- and coordinate-dependent reaction rates and could be used for a wide variety of stochastic systems. The results reveal two qualitative regimes. If the gradient is shallow, the rate of adaptation is limited by dispersion and increases linearly with the gradient slope. If the gradient is steep, the adaptation rate is limited by mutation. In this regime, the mean-field result is highly misleading: it predicts that the adaptation rate continues to increase with the gradient slope, whereas stochastic simulations show that it in fact decreases with the square root of the slope. This discrepancy underscores the importance of discreteness and stochasticity even at high population densities; mean-field results, including those routinely used in quantitative genetics, should be interpreted with care.

Adaptation of Escherichia coli traversing from the faecal environment to the urinary tract

The majority of extraintestinal pathogenic Escherichia coli (ExPEC) causing urinary tract infections (UTI) are found in the patient's own gut flora, but only limited knowledge is available on the potential adaptation that may occur in the bacteria in order to traverse the perineum and successfully infect the urinary tract. Here, matching pairs of faecal and UTI isolates from 42 patients were compared pairwise using in-depth whole-genome sequencing to investigate whether genetic changes were evident for successful colonization in these two different environments. The identified non-synonymous mutations (0-12 substitutions in each pair) were primarily associated to genes encoding virulence factors and nutrient metabolism; and indications of parallel evolution were observed in genes encoding the major phase-variable protein antigen 43, a toxin/antitoxin locus and haemolysin B. No differences in virulence potential were observed in a mouse UTI model for five matching faecal and UTI isolates with or without mutations in antigen 43 and haemolysin B. Variations in plasmid content were observed in only four of the 42 pairs. Although, we observed mutations in known UTI virulence genes for a few pairs, the majority showed no detectable differences with respect to mutations or mobilome when compared to their faecal counterpart. The results show that UPECs are successful in colonizing both the bladder and gut without adaptation.

Antimicrobial Resistance Profiles and Diversity in Salmonella from Humans and Cattle, 2004-2011

Analysis of long-term anti-microbial resistance (AMR) data is useful to understand source and transmission dynamics of AMR. We analysed 5124 human clinical isolates from Washington State Department of Health, 391 cattle clinical isolates from the Washington Animal Disease Diagnostic Laboratory and 1864 non-clinical isolates from foodborne disease research on dairies in the Pacific Northwest. Isolates were assigned profiles based on phenotypic resistance to 11 anti-microbials belonging to eight classes. Salmonella Typhimurium (ST), Salmonella Newport (SN) and Salmonella Montevideo (SM) were the most common serovars in both humans and cattle. Multinomial logistic regression showed ST and SN from cattle had greater probability of resistance to multiple classes of anti-microbials than ST and SN from humans (P < 0.0001). While these findings could be consistent with the belief that cattle are a source of resistant ST and SN for people, occurrence of profiles unique to cattle and not observed in temporally related human isolates indicates these profiles are circulating in cattle only. We used various measures to assess AMR diversity, conditional on the weighting of rare versus abundant profiles. AMR profile richness was greater in the common serovars from humans, although both source data sets were dominated by relatively few profiles. The greater profile richness in human Salmonella may be due to greater diversity of sources entering the human population compared to cattle or due to continuous evolution in the human environment. Also, AMR diversity was greater in clinical compared to non-clinical cattle Salmonella, and this could be due to anti-microbial selection pressure in diseased cattle that received treatment. The use of bootstrapping techniques showed that although there were shared profiles between humans and cattle, the expected and observed number of profiles was different, suggesting Salmonella and associated resistance from humans and cattle may not be wholly derived from a common population.

Frequency and fitness consequences of bacteriophage φ6 host range mutations

Viruses readily mutate and gain the ability to infect novel hosts, but few data are available regarding the number of possible host range-expanding mutations allowing infection of any given novel host, and the fitness consequences of these mutations on original and novel hosts. To gain insight into the process of host range expansion, we isolated and sequenced 69 independent mutants of the dsRNA bacteriophage Φ6 able to infect the novel host, Pseudomonas pseudoalcaligenes. In total, we found at least 17 unique suites of mutations among these 69 mutants. We assayed fitness for 13 of 17 mutant genotypes on P. pseudoalcaligenes and the standard laboratory host, P. phaseolicola. Mutants exhibited significantly lower fitnesses on P. pseudoalcaligenes compared to P. phaseolicola. Furthermore, 12 of the 13 assayed mutants showed reduced fitness on P. phaseolicola compared to wildtype Φ6, confirming the prevalence of antagonistic pleiotropy during host range expansion. Further experiments revealed that the mechanistic basis of these fitness differences was likely variation in host attachment ability. In addition, using computational protein modeling, we show that host-range expanding mutations occurred in hotspots on the surface of the phage's host attachment protein opposite a putative hydrophobic anchoring domain.

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

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

Evolution of enterohemorrhagic escherichia coli O26 based on single-nucleotide polymorphisms

Enterohemorrhagic Escherichia coli (EHEC) O26:H11/H⁻ is the predominant non-O157 EHEC serotype among patients with diarrhea, bloody diarrhea, and hemolytic uremic syndrome (HUS) worldwide. To elucidate their phylogeny and association between their phylogenetic background and clinical outcome of the infection, we investigated 120 EHEC O26:H11/H⁻ strains isolated between 1965 and 2012 from asymptomatic carriers and patients with diarrhea or HUS. Whole-genome shotgun sequencing (WGS) was applied to ten representative EHEC O26 isolates to determine single nucleotide polymorphism (SNP) localizations within a predefined set of core genes. A multiplex SNP assay, comprising a randomly distributed subset of 48 SNPs, was established to detect SNPs in 110 additional EHEC O26 strains. Within approximately 1 Mb of core genes, WGS resulted in 476 high-quality bi-allelic SNP localizations. Forty-eight of these were subsequently investigated in 110 EHEC O26 and four different SNP clonal complexes (SNP-CC) were identified. SNP-CC2 was significantly associated with the development of HUS. Within the subsequently established evolutionary model of EHEC O26, we dated the emergence of human EHEC O26 to approximately 19,700 years ago and demonstrated a recent evolution within humans into the 4 SNP-CCs over the past 1,650 years. WGS and subsequent SNP typing enabled us to gain new insights into the evolution of EHEC O26 suggesting a common theme in this EHEC group with analogies to EHEC O157. In addition, the SNP-CC analysis may help to assess a risk in infected individuals for the progression to HUS and to implement more specific infection control measures.

Genomic diversity and fitness of E. coli strains recovered from the intestinal and urinary tracts of women with recurrent urinary tract infection

Urinary tract infections (UTIs) are common in women, and recurrence is a major clinical problem. Most UTIs are caused by uropathogenic Escherichia coli (UPEC). UPEC are generally thought to migrate from the gut to the bladder to cause UTI. UPEC form specialized intracellular bacterial communities in the bladder urothelium as part of a pathogenic mechanism to establish a foothold during acute stages of infection. Evolutionarily, such a specific adaptation to the bladder environment would be predicted to result in decreased fitness in other habitats, such as the gut. To examine this prediction, we characterized 45 E. coli strains isolated from the feces and urine of four otherwise healthy women with recurrent UTI. Multilocus sequence typing and whole genome sequencing revealed that two patients maintained a clonal population in both these body habitats throughout their recurrent UTIs, whereas the other two exhibited a wholesale shift in the dominant UPEC strain colonizing both sites. In vivo competition studies in mouse models, using isolates taken from one of the patients with a wholesale population shift, revealed that the strain that dominated her last UTI episode had increased fitness in both the gut and the bladder relative to the strain that dominated in preceding episodes. Increased fitness correlated with differences in the strains' gene repertoires and carbohydrate and amino acid utilization profiles. Thus, UPEC appear capable of persisting in both the gut and urinary tract without a fitness trade-off, emphasizing the need to widen our consideration of potential reservoirs for strains causing recurrent UTI.

Structural and population characterization of MrkD, the adhesive subunit of type 3 fimbriae

Type 3 fimbriae are adhesive organelles found in enterobacterial pathogens. The fimbriae promote biofilm formation on biotic and abiotic surfaces; however, the exact identity of the receptor for the type 3 fimbriae adhesin, MrkD, remains elusive. We analyzed naturally occurring structural and functional variabilities of the MrkD adhesin from Klebsiella pneumoniae and Escherichia coli isolates of diverse origins. We identified a total of 33 allelic variants of mrkD among 90 K. pneumoniae isolates and 10 allelic variants among 608 E. coli isolates, encoding 11 and 9 protein variants, respectively. Based on the level of accumulated silent variability between the alleles, mrkD was acquired a relatively long time ago in K. pneumoniae but recently in E. coli. However, unlike K. pneumoniae, mrkD in E. coli is actively evolving under a strong positive selection by accumulation of mutations, often targeting the same positions in the protein. Several naturally occurring MrkD protein variants from E. coli were found to be significantly less adherent when tested in a mannan-binding assay and showed reduced biofilm-forming capacity. Functional examination of the MrkD adhesin in flow chamber experiments determined that it interacts with Saccharomyces cerevisiae cells in a shear-dependent manner, i.e., the binding is catch-bond-like and enhanced under increasing shear conditions. Homology modeling strongly suggested that MrkD has a two-domain structure, comprising a pilin domain anchoring the adhesin to the fimbrial shaft and a lectin domain containing the binding pocket; this is similar to structures found in other catch-bond-forming fimbrial adhesins in enterobacteria.

The clonal distribution and diversity of extraintestinal Escherichia coli isolates vary according to patient characteristics

The clonal distribution of Escherichia coli across an unselected population in the current era of widespread antimicrobial resistance is incompletely defined. In this study, we used a newly described clonal typing strategy based on sequencing of fumC and fimH (i.e., CH typing) to infer multilocus sequence types (STs) for 299 consecutive, nonduplicate extraintestinal E. coli isolates from all cultures submitted to Olmsted County, MN, laboratories in February and March 2011 and then compared STs with epidemiological data. Forty-seven different STs were identified, most commonly ST131 (27%), ST95 (11%), ST73 (8%), ST127 (6%), and ST69 (5%). Isolates from these five STs comprised two-thirds of health care-associated (HA) isolates but only half of community-associated (CA) isolates. ST131 was represented overwhelmingly (88%) by a single recently expanded H30 subclone, which was the most extensively antimicrobial-resistant subclone overall and was especially predominant in HA infections and among adults >50 years old. In contrast, among patients 11 to 50 years old, ST69, -95, and -73 were more common. Because of the preponderance of the H30 subclone of ST131, ST diversity was lower among HA than CA isolates, and among antimicrobial-resistant than antimicrobial-susceptible isolates, which otherwise had similar ST distributions. In conclusion, in this U.S. Midwest region, the distribution and diversity of STs among extraintestinal E. coli clinical isolates vary by patient age, type of infection, and resistance phenotype. ST131 predominates among young children and the elderly, HA infections, and antimicrobial-resistant isolates, whereas other well-known pathogenic lineages are more common among adolescents and young adults, CA infections, and antimicrobial-susceptible isolates.

Predictive diagnostics for Escherichia coli infections based on the clonal association of antimicrobial resistance and clinical outcome

The ability to identify bacterial pathogens at the subspecies level in clinical diagnostics is currently limited. We investigated whether splitting Escherichia coli species into clonal groups (clonotypes) predicts antimicrobial susceptibility or clinical outcome. A total of 1,679 extraintestinal E. coli isolates (collected from 2010 to 2012) were collected from one German and 5 U.S. clinical microbiology laboratories. Clonotype identity was determined by fumC and fimH (CH) sequencing. The associations of clonotype with antimicrobial susceptibility and clinical variables were evaluated. CH typing divided the isolates into >200 CH clonotypes, with 93% of the isolates belonging to clonotypes with ≥ 2 isolates. Antimicrobial susceptibility varied substantially among clonotypes but was consistent across different locations. Clonotype-guided antimicrobial selection significantly reduced "drug-bug" mismatch compared to that which occurs with the use of conventional empirical therapy. With trimethoprim-sulfamethoxazole and fluoroquinolones, the drug-bug mismatch was predicted to decrease 62% and 78%, respectively. Recurrent or persistent urinary tract infection and clinical sepsis were significantly correlated with specific clonotypes, especially with CH40-30 (also known as H30), a recently described clonotype within sequence type 131 (ST131). We were able to clonotype directly from patient urine samples within 1 to 3 h of obtaining the specimen. In E. coli, subspecies-level identification by clonotyping can be used to significantly improve empirical predictions of antimicrobial susceptibility and clinical outcomes in a timely manner.

Convergent molecular evolution of genomic cores in Salmonella enterica and Escherichia coli

One of the strongest signals of adaptive molecular evolution of proteins is the occurrence of convergent hot spot mutations: repeated changes in the same amino acid positions. We performed a comparative genome-wide analysis of mutation-driven evolution of core (omnipresent) genes in 17 strains of Salmonella enterica subspecies I and 22 strains of Escherichia coli. More than 20% of core genes in both Salmonella and E. coli accumulated hot spot mutations, with a predominance of identical changes having recent evolutionary origin. There is a significant overlap in the functional categories of the adaptively evolving genes in both species, although mostly via separate molecular mechanisms. As a strong evidence of the link between adaptive mutations and virulence in Salmonella, two human-restricted serovars, Typhi and Paratyphi A, shared the highest number of genes with serovar-specific hot spot mutations. Many of the core genes affected by Typhi/Paratyphi A-specific mutations have known virulence functions. For each species, a list of nonrecombinant core genes (and the hot spot mutations therein) under positive selection is provided.

On the rapidity of antibiotic resistance evolution facilitated by a concentration gradient

The rapid emergence of bacterial strains resistant to multiple antibiotics is posing a growing public health risk. The mechanisms underlying the rapid evolution of drug resistance are, however, poorly understood. The heterogeneity of the environments in which bacteria encounter antibiotic drugs could play an important role. E.g., in the highly compartmentalized human body, drug levels can vary substantially between different organs and tissues. It has been proposed that this could facilitate the selection of resistant mutants, and recent experiments support this. To study the role of spatial heterogeneity in the evolution of drug resistance, we present a quantitative model describing an environment subdivided into relatively isolated compartments with various antibiotic concentrations, in which bacteria evolve under the stochastic processes of proliferation, migration, mutation and death. Analytical and numerical results demonstrate that concentration gradients can foster a mode of adaptation that is impossible in uniform environments. It allows resistant mutants to evade competition and circumvent the slow process of fixation by invading compartments with higher drug concentrations, where less resistant strains cannot subsist. The speed of this process increases sharply with the sensitivity of the growth rate to the antibiotic concentration, which we argue to be generic. Comparable adaptation rates in uniform environments would require a high selection coefficient (s > 0.1) for each forward mutation. Similar processes can occur if the heterogeneity is more complex than just a linear gradient. The model may also be applicable to other adaptive processes involving environmental heterogeneity and range expansion.

Adaptive evolution of class 5 fimbrial genes in enterotoxigenic Escherichia coli and its functional consequences

Class 5 fimbriae of enterotoxigenic Escherichia coli (ETEC) comprise eight serologically discrete colonization factors that mediate small intestinal adhesion. Their differentiation has been attributed to the pressure imposed by host adaptive immunity. We sequenced the major pilin and minor adhesin subunit genes of a geographically diverse population of ETEC elaborating CFA/I (n = 31), CS17 (n = 20), and CS2 (n = 18) and elucidated the functional effect of microevolutionary processes. Between the fimbrial types, the pairwise nucleotide diversity for the pilin or adhesin genes ranged from 35-43%. Within each fimbrial type, there were 17 non-synonymous and 1 synonymous point mutations among all pilin or adhesin gene copies, implying that each fimbrial type was acquired by ETEC strains very recently, consistent with a recent origin of this E. coli pathotype. The 17 non-synonymous allelic differences occurred in the CFA/I pilin gene cfaB (two changes) and adhesin gene cfaE (three changes), and CS17 adhesin gene csbD (12 changes). All but one amino acid change in the adhesins clustered around the predicted ligand-binding pocket. Functionally, these changes conferred an increase in cell adhesion in a flow chamber assay. In contrast, the two mutations in the non-adhesive CfaB subunit localized to the intersubunit interface and significantly reduced fimbrial adhesion in this assay. In conclusion, naturally occurring mutations in the ETEC adhesive and non-adhesive subunits altered function, were acquired under positive selection, and are predicted to impact bacteria-host interactions.

Sources and sinks: a stochastic model of evolution in heterogeneous environments

We study evolution driven by spatial heterogeneity in a stochastic model of source-sink ecologies. A sink is a habitat where mortality exceeds reproduction so that a local population persists only due to immigration from a source. Immigrants can, however, adapt to conditions in the sink by mutation. To characterize the adaptation rate, we derive expressions for the first arrival time of adapted mutants. The joint effects of migration, mutation, birth, and death result in two distinct parameter regimes. These results may pertain to the rapid evolution of drug-resistant pathogens and insects.

Escherichia coli O157

Escherichia coli O157 is an uncommon but serious cause of gastroenteritis. This bacterium is noteworthy because a few, but significant, number of infected people develop the haemolytic uraemic syndrome, which is the most frequent cause of acute renal failure in children in the Americas and Europe. Many infections of E coli O157 could be prevented by the more effective application of evidence-based methods, which is especially important because once an infection has been established, no therapeutic interventions are available to lessen the risk of the development of the haemolytic uraemic syndrome. This Review takes into account the evolution and geographical distibution of E coli O157 (and its close pathogenic relatives); the many and varied routes of transmission from its major natural hosts, ruminant farm animals; and other aspects of its epidemiology, its virulence factors, the diagnosis and management of infection and their complications, the repercussions of infection including costs, and prevention.

Experimental evidence that source genetic variation drives pathogen emergence

A pathogen can readily mutate to infect new host types, but this does not guarantee successful establishment in the new habitat. What factors, then, dictate emergence success? One possibility is that the pathogen population cannot sustain itself on the new host type (i.e. host is a sink), but migration from a source population allows adaptive sustainability and eventual emergence by delivering beneficial mutations sampled from the source's standing genetic variation. This idea is relevant regardless of whether the sink host is truly novel (host shift) or whether the sink is an existing or related, similar host population thriving under conditions unfavourable to pathogen persistence (range expansion). We predicted that sink adaptation should occur faster under range expansion than during a host shift owing to the effects of source genetic variation on pathogen adaptability in the sink. Under range expansion, source migration should benefit emergence in the sink because selection acting on source and sink populations is likely to be congruent. By contrast, during host shifts, source migration is likely to disrupt emergence in the sink owing to uncorrelated selection or performance tradeoffs across host types. We tested this hypothesis by evolving bacteriophage populations on novel host bacteria under sink conditions, while manipulating emergence via host shift versus range expansion. Controls examined sink adaptation when unevolved founding genotypes served as migrants. As predicted, adaptability was fastest under range expansion, and controls did not adapt. Large, similar and similarly timed increases in fitness were observed in the host-shift populations, despite declines in mean fitness of immigrants through time. These results suggest that source populations are the origin of mutations that drive adaptive emergence at the edge of a pathogen's ecological or geographical range.

Darwinian interventions: taming pathogens through evolutionary ecology

From the seemingly inevitable emergence of antibiotic resistance following drug use to the necessity of annual influenza vaccine updates, the ability of pathogens to evolve appears to preclude the development of effective, long-term interventions for many diseases. Despite this gloomy forecast, recent research explores how various principles and concepts from evolutionary ecology might inform practical attempts to bring these disease-causing agents under greater human control. By utilizing evolutionary and ecological information at various scales, these works investigate some promising avenues of disease intervention, from trapping pathogens in evolutionary dead ends, through slowing or inhibiting the process of pathogen adaptation, to more accurate forecasting of pathogen evolutionary trajectories. Thus, an evolutionary ecology perspective might ultimately provide powerful new tools in the pursuit of enduring solutions to the problem of treatment-driven pathogen evolution.

Adhesive organelles of Gram-negative pathogens assembled with the classical chaperone/usher machinery: structure and function from a clinical standpoint

This review summarizes current knowledge on the structure, function, assembly and biomedical applications of the superfamily of adhesive fimbrial organelles exposed on the surface of Gram-negative pathogens with the classical chaperone/usher machinery. High-resolution three-dimensional (3D) structure studies of the minifibers assembling with the FGL (having a long F1-G1 loop) and FGS (having a short F1-G1 loop) chaperones show that they exploit the same principle of donor-strand complementation for polymerization of subunits. The 3D structure of adhesive subunits bound to host-cell receptors and the final architecture of adhesive fimbrial organelles reveal two functional families of the organelles, respectively, possessing polyadhesive and monoadhesive binding. The FGL and FGS chaperone-assembled polyadhesins are encoded exclusively by the gene clusters of the γ3- and κ-monophyletic groups, respectively, while gene clusters belonging to the γ1-, γ2-, γ4-, and π-fimbrial clades exclusively encode FGS chaperone-assembled monoadhesins. Novel approaches are suggested for a rational design of antimicrobials inhibiting the organelle assembly or inhibiting their binding to host-cell receptors. Vaccines are currently under development based on the recombinant subunits of adhesins.

Bacteriophages as model organisms for virus emergence research

Viruses fully emerge by gaining the ability to sustainably infect new host populations. When the hosts are humans, emerging viruses can present major public health issues, as exemplified by the AIDS pandemic. Therefore, heuristic approaches to identify nascent diseases before they become pandemic would be valuable. Unfortunately, the current patient-based and epidemiological approaches are ill-suited in this regard because they are largely responsive and not predictive. Alternative approaches based on virus evolutionary ecology might have greater potential to predict virus emergence. However, given the difficulties encountered when studying metazoan viruses in this context, the development of new model systems is greatly desirable. Here, I highlight studies that show that bacteriophages are appropriate model organisms for virus emergence research because of the ease in which important population parameters can be manipulated. Ideally this research will permit identifying major factors determining the persistence or extinction of emerging viruses. If such viruses could be recognized in advance, patient-based and epidemiological strategies could be better mobilized to deal with them.

High frequency of hotspot mutations in core genes of Escherichia coli due to short-term positive selection

Core genes comprising the ubiquitous backbone of bacterial genomes are not subject to frequent horizontal transfer and generally are not thought to contribute to the adaptive evolution of bacterial pathogens. We determined, however, that at least one-third and possibly more than one-half of the core genes in Escherichia coli genomes are targeted by repeated replacement substitutions in the same amino acid positions-hotspot mutations. Occurrence of hotspot mutations is driven by positive selection, as their rate is significantly higher than expected by random chance alone, and neither intragenic recombination nor increased mutability can explain the observed patterns. Also, commensal E. coli strains have a significantly lower frequency of mutated genes and mutations per genome than pathogenic strains. E. coli strains causing extra-intestinal infections accumulate hotspot mutations at the highest rate, whereas the highest total number of mutated genes has been found among Shigella isolates, suggesting the pathoadaptive nature of such mutations. The vast majority of hotspot mutations are of recent evolutionary origin, implying short-term positive selection, where adaptive mutations emerge repeatedly but are not sustained in natural circulation for long. Such pattern of dynamics is consistent with source-sink model of virulence evolution.

Comparative evolutionary analysis of the major structural subunit of Vibrio vulnificus type IV pili

Type IV pili contribute to virulence in Vibrio vulnificus, the bacterium responsible for the majority of fatal seafood-related infections. Here, we performed within- and between-species evolutionary analysis of the gene that encodes the major structural subunit of the pilus, pilA, by comparing it with pilD and gyrB, the genes encoding the type IV prepilin peptidase and beta subunit of DNA gyrase, respectively. Although the diversity in pilD and gyrB is similar to each other and likely to have accumulated after speciation of V. vulnificus, pilA is several times more diverse at both nonsynonymous and synonymous levels. Also, in contrast to pilD and gyrB, there are virtually unrestricted and highly localized horizontal movements of pilA alleles between the major phylogenetic groups of V. vulnificus. The frequent movement of pilA involves homologous recombination of the entire gene with no evidence for intragenic recombination between the alleles. We propose that pilA allelic diversity and horizontal movement is maintained in the population by both diversifying and frequency-dependent selection most likely to escape shellfish innate immunity defense or lytic phages. Other possibilities leading to such selection dynamics of V. vulnificus pilA could involve adaptation to diverse host populations or within-host compartments, or natural DNA uptake and transformation. We show that the history of nucleotide diversification in pilA predates V. vulnificus speciation and this diversification started at or before the time of the last common ancestor for V. vulnificus, Vibrio parahaemolyticus, and Vibrio cholerae. At the same time, it appears that within the various pilA groups of V. vulnificus, there is no positive selection for structural mutations and consequently no evidence for source-sink selection. In contrast, pilD has accumulated a number of apparently adaptive mutations in the regions encoding the membrane-spanning portions of the prepilin peptidase, possibly affecting fimbrial expression and/or function, and is being subjected to source-sink selection dynamics.

A precise reconstruction of the emergence and constrained radiations of Escherichia coli O157 portrayed by backbone concatenomic analysis

Single nucleotide polymorphisms (SNPs) in stable genome regions provide durable measurements of species evolution. We systematically identified each SNP in concatenations of all backbone ORFs in 7 newly or previously sequenced evolutionarily instructive pathogenic Escherichia coli O157:H7, O157:H(-), and O55:H7. The 1,113 synonymous SNPs demonstrate emergence of the largest cluster of this pathogen only in the last millennium. Unexpectedly, shared SNPs within circumscribed clusters of organisms suggest severely restricted survival and limited effective population sizes of pathogenic O157:H7, tenuous survival of these organisms in nature, source-sink evolutionary dynamics, or, possibly, a limited number of mutations that confer selective advantage. A single large segment spanning the rfb-gnd gene cluster is the only backbone region convincingly acquired by recombination as O157 emerged from O55. This concatenomic analysis also supports using SNPs to differentiate closely related pathogens for infection control and forensic purposes. However, constrained radiations raise the possibility of making false associations between isolates.

The rate of environmental change drives adaptation to an antibiotic sink

Recent accelerated trends of human-induced global changes are providing many examples of adaptation to novel environments. Although the rate of environmental change can vary dramatically, its effect on evolving populations is unknown. A crucial feature explaining the adaptation to harsh environments is the supply of beneficial mutations via immigration from a 'source' population. In this study, we tested the effect of immigration on adaptation to increasing concentrations of antibiotics. Using experimental population of Pseudomonas aeruginosa, a pathogenic bacterium, we show that higher immigration rates and a slow increase in antibiotic concentration result in a more rapid evolution of resistance; however, a high immigration rate combined with rapid increases in concentration resulted in higher maximal levels of resistance. These findings, which support recent theoretical work, have important implications for the control of antibiotic resistance because they show that rapid rates of change can produce variants with the ability to resist future treatments.

Source-sink dynamics shape the evolution of antibiotic resistance and its pleiotropic fitness cost

Understanding the conditions that favour the evolution and maintenance of antibiotic resistance is the central goal of epidemiology. A crucial feature explaining the adaptation to harsh, or 'sink', environments is the supply of beneficial mutations via migration from a 'source' population. Given that antibiotic resistance is frequently associated with antagonistic pleiotropic fitness costs, increased migration rate is predicted not only to increase the rate of resistance evolution but also to increase the probability of fixation of resistance mutations with minimal fitness costs. Here we report in vitro experiments using the nosocomial pathogenic bacterium Pseudomonas aeruginosa that support these predictions: increasing rate of migration into environments containing antibiotics increased the rate of resistance evolution and decreased the associated costs of resistance. Consistent with previous theoretical work, we found that resistance evolution arose more rapidly in the presence of a single antibiotic than two. Evolution of resistance was also more rapid when bacteria were subjected to sequential exposure with two antibiotics (cycling therapy) compared with simultaneous exposure (bi-therapy). Furthermore, pleiotropic fitness costs of resistance to two antibiotics were higher than for one antibiotic, and were also higher under bi-therapy than cycling therapy, although the cost of resistance depended on the order of the antibiotics through time. These results may be relevant to the clinical setting where immigration is known to be important between chemotherapeutically treated patients, and demonstrate the importance of ecological and evolutionary dynamics in the control of antibiotic resistance.

ZPS: visualization of recent adaptive evolution of proteins

Background

Detection of adaptive amino acid changes in proteins under recent short-term selection is of great interest for researchers studying microevolutionary processes in microbial pathogens or any other biological species. However, independent occurrence of such point mutations within genetically diverse haplotypes makes it difficult to detect the selection footprint by using traditional molecular evolutionary analyses. The recently developed Zonal Phylogeny (ZP) has been shown to be a useful analytic tool for identifying the footprints of short-term positive selection. ZP separates protein-encoding genes into evolutionarily long-term (with silent diversity) and short-term (without silent diversity) categories, or zones, followed by statistical analysis to detect signs of positive selection in the short-term zone. However, successful broad application of ZP for analysis of large haplotype datasets requires automation of the relatively labor-intensive computational process.

Results

Here we present Zonal Phylogeny Software (ZPS), an application that describes the distribution of single nucleotide polymorphisms (SNPs) of synonymous (silent) and non-synonymous (replacement) nature along branches of the DNA tree for any given protein-coding gene locus. Based on this information, ZPS separates the protein variant haplotypes with silent variability (Primary zone) from those that have recently evolved from the Primary zone variants by amino acid changes (External zone). Further comparative analysis of mutational hot-spot frequencies and haplotype diversity between the two zones allows determination of whether the External zone haplotypes emerged under positive selection.

Conclusions

As a visualization tool, ZPS depicts the protein tree in a DNA tree, indicating the most parsimonious numbers of synonymous and non-synonymous changes along the branches of a maximum-likelihood based DNA tree, along with information on homoplasy, reversion and structural mutation hot-spots. Through zonal differentiation, ZPS allows detection of recent adaptive evolution via selection of advantageous structural mutations, even when the advantage conferred by such mutations is relatively short-term (as in the case of "source-sink" evolutionary dynamics, which may represent a major mode of virulence evolution in microbes).