Pathogen population structure can explain hospital outbreaks.

Hospitalized patients are at risk for increased length of stay, illness, or death due to hospital acquired infections. The majority of hospital transmission models describe dynamics on the level of the host rather than on the level of the pathogens themselves. Accordingly, epidemiologists often cannot complete transmission chains without direct evidence of either host-host contact or a large reservoir population. Here, we propose an ecology-based model to explain the transmission of pathogens in hospitals. The model is based upon metapopulation biology, which describes a group of interacting localized populations and island biogeography, which provides a basis for how pathogens may be moving between locales. Computational simulation trials are used to assess the applicability of the model. Results indicate that pathogens survive for extended periods without the need for large reservoirs by living in localized ephemeral populations while continuously transmitting pathogens to new seed populations. Computational simulations show small populations spending significant portions of time at sizes too small to be detected by most surveillance protocols and that the number and type of these ephemeral populations enable the overall pathogen population to be sustained. By modeling hospital pathogens as a metapopulation, many observations characteristic of hospital acquired infection outbreaks for which there has previously been no sufficient biological explanation, including how and why empirically successful interventions work, can now be accounted for using population dynamic hypotheses. Epidemiological links between temporally isolated outbreaks are explained via pathogen population dynamics and potential outbreak intervention targets are identified.

Different adaptive strategies in E. coli populations evolving under macronutrient limitation and metal ion limitation.

BACKGROUND: Adaptive responses to nutrient limitation involve mutations that increase the efficiency of usage or uptake of the limiting nutrient. However, starvation of different nutrients has contrasting effects on physiology, resulting in different evolutionary responses. Most studies performed to understand these evolutionary responses have focused only on macronutrient limitation. Hence our understanding of adaptation under limitation of other forms of nutrients is limited. In this study, we compared the evolutionary response in populations evolving under growth-limiting conditions for a macronutrient and a major cation. RESULTS: We evolved eight populations of E. coli in nutrient-limited chemostats for 400 generations to identify the genetic basis of the mechanisms involved in efficient usage of two nutrients: nitrogen and magnesium. Our population genomic sequencing work, based on this study and previous work, allowed us to identify targets of selection under these nutrient limiting conditions. Global transcriptional regulators glnGL were targets of selection under nitrogen starvation, while proteins involved in outer-membrane biogenesis (genes from the lpt operon) were targets of selection under magnesium starvation. The protein involved in cell-cycle arrest (yhaV) was a target of selection in both environments. We re-constructed specific mutants to analyze the effect of individual mutations on fitness in nutrient limiting conditions in chemostats and in batch cultures. We further demonstrated that adaptation to nitrogen starvation proceeds via a nutrient specific mechanism, while that to magnesium starvation involves a more general mechanism. CONCLUSIONS: Our results show two different forms of adaptive strategies under limitation of nutrients that effect cellular physiology in different ways. Adaptation to nitrogen starvation proceeds by upregulation of transcriptional regulator glnG and subsequently of transporter protein amtB, both of which results in increased nitrogen scavenging ability of the cell. On the other hand, adaptation to magnesium starvation proceeds via the restructuring of the cell outer-membrane, allowing magnesium to be redistributed to other biological processes. Also, adaptation to the chemostat environment involves selection for loss of function mutations in genes that under nutrient-limiting conditions interfere with continuous growth.

Evolutionary implications of Liebig's law of the minimum: Selection under low concentrations of two nonsubstitutable nutrients.

Interactions between different axes of an organism's niche determine the evolutionary trajectory of a population. An extreme case of these interactions is predicted from ecological theory in Liebig's law of the minimum. This law states that in environments where multiple nutrients are in relatively low concentrations, only one nutrient will affect the growth of the organism. This implies that the evolutionary response of the population would be dictated by the most growth-limiting nutrient. Alternatively, it is possible that an initial adaptation to the most limiting nutrient results in other nutrients present in low concentration affecting the evolutionary dynamics of the population. To test these hypotheses, we conducted twelve evolution experiments in chemostats using Escherichia coli populations: four under nitrogen limitation, four under magnesium limitation, and four in which both nitrogen and magnesium are in low concentrations. In the last environment, only magnesium seems to limit growth (Low Nitrogen Magnesium Limited environment, LNML). We observe a decrease in nitrogen concentration in the LNML environment over the course of our evolution experiment indicating that nitrogen might become limiting in these environments. Genetic reconstruction results show that clones adapted to magnesium limitation have genes involved in nitrogen starvation, that is, glnG (nitrogen starvation transcriptional regulator) and amtB (transport protein) to be upregulated only in the LNML environment as compared to magnesium-limiting environments. Together, our results highlights that in low-nutrient environments, adaptation to the growth-limiting nutrient results in other nutrients at low concentrations to play a role in the evolutionary dynamics of the population.

Evolution of Resistance to Continuously Increasing Streptomycin Concentrations in Populations of Escherichia coli.

The evolution of antibiotic resistance in bacteria has become one of the defining problems in modern biology. Bacterial resistance to antimicrobial therapy threatens to eliminate one of the pillars of the practice of modern medicine. Yet, in spite of the importance of this problem, only recently have the dynamics of the shift from antibiotic sensitivity to resistance in a bacterial population been studied. In this study, a novel chemostat method was used to observe the evolution of resistance to streptomycin in a sensitive population of Escherichia coli, which grew while the concentration of antibiotic was constantly increasing. The results indicate that resistant mutants remain at a low frequency for longer than expected and do not begin to rise to a high frequency until the antibiotic concentrations are above the measured MIC, creating a "lull period" in which there were few bacterial cells growing in the chemostats. Overall, mutants resistant to streptomycin were found in >60% of the experimental trial replicates. All of the mutants detected were found to have MICs far above the maximum levels of streptomycin to which they were exposed and reached a high frequency within 96 h.

Tracking recent adaptive evolution in microbial species using TimeZone.

An important goal of the analysis of sequenced genomes of microbial pathogens is to improve the therapy of infectious diseases. In this context, a major challenge is to detect genomic-level evolutionary changes that increase microbial virulence. TimeZone, a genome analysis software package, is designed to detect footprints of positive selection for functionally adaptive point mutations. The uniqueness of TimeZone lies in its ability to predict recent adaptive mutations that are overlooked by conventional microevolutionary tools. This protocol describes the use of TimeZone to analyze adaptive footprints in either individual genes or in sets of genomes. Three major workflows are described: (i) extraction of orthologous gene sets from multiple genomes; (ii) alignment and phylogenetic analysis of genes; and (iii) identification of candidate genes under positive selection for point mutations, taking into account the effect of recombination events. This software package can be downloaded free from http://sourceforge.net/projects/timezone1/. In the case, for example, of the analysis of 14 Escherichia coli genomes, the protocol described here can be completed in ∼32 h.

Waste and yet want not.

In this issue of Molecular Cell, Shachrai et al. (2010) demonstrate that the cost of wasteful protein expression in E. coli is specific to the transition from stationary phase to balanced exponential growth, probably because of a shortage of ribosomes during this growth phase.

Evolution of northeastern and midwestern Borrelia burgdorferi, United States.

The per capita incidence of human Lyme disease in the northeastern United States is more than twice that in the Midwest. However, the prevalence of Borrelia burgdorferi, the bacterium that causes Lyme disease, in the tick vector is nearly identical in the 2 regions. The disparity in human Lyme disease incidence may result from a disparity in the human invasiveness of the bacteria in the Northeast and Midwest caused by fundamentally different evolutionary histories. B. burgdorferi populations in the Northeast and Midwest are geographically isolated, enabling evolutionary divergence in human invasiveness. However, we found that B. burgdorferi populations in the Northeast and Midwest shared a recent common ancestor, which suggests that substantial evolutionary divergence in human invasiveness has not occurred. We propose that differences in either animal ecology or human behavior are the root cause of the differences in human incidence between the 2 regions.

Transcription, translation, and the evolution of specialists and generalists.

We used DNA microarrays to measure transcription and iTRAQ 2D liquid chromatography-mass spectrometry/mass spectrometry (a mass-tag labeling proteomic technique) to measure protein expression in 14 strains of Escherichia coli adapted for hundreds of generations to growth-limiting concentrations of either lactulose, methylgalactoside, or a 72:28 mixture of the two. The two ancestors, TD2 and TD10, differ only in their lac operons and have similar transcription and protein expression profiles. Changes in transcription and protein expression are observed at 30-250 genes depending on the evolved strain. Lactulose specialists carry duplications of the lac operon and show increased transcription and translation at lac. Methylgalactoside specialists are galS(-) and so constitutively transcribe and translate mgl, which encodes a transporter of methylgalactoside. However, there are two strains that carry lac duplications, are galS(-), and show increased transcription and translation at both operons. One is a generalist, the other a lactulose specialist. The generalist fails to sweep to fixation because its lac(+), galS(+) competitor expresses the csg adhesin and sticks to the chemostat wall, thereby preventing complete washout. Transcription and translation are sometimes decoupled. Lactulose-adapted strains show increased protein expression at fru, a fructose transporter, without evidence of increased transcription. This suggests that fructose, produced by the action of beta-galactosidase on lactulose, may leach from cells before being recouped. Reduced expression, at "late" flagella genes and the constitutive gat operon, is an adaptation to starvation. A comparison with two other long-term evolution experiments suggests that certain aspects of adaptation are predictable, some are characteristic of an experimental system, whereas others seem erratic.

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.

The propensity of different Borrelia burgdorferi sensu stricto genotypes to cause disseminated infections in humans.

Lineages of Borrelia burgdorferi, the bacterium that causes Lyme disease, can be characterized by distinct alleles at the outer surface protein C (ospC) locus. The lineages marked by ospC genotypes have been shown to be differentially invasive in different species of mammals, including humans; genotypes A, B, I, and K effectively disseminate to human blood and cerebrospinal fluid. In this report, we extend the sample of genotypes isolated from human blood to include genotypes N, H, C, M, and D, and rank each by their probability of disseminating from ticks to the blood of humans. Our results demonstrate that only some genotypes of B. burgdorferi present in ticks have a high propensity to disseminate in humans.

The cost of expression of Escherichia coli lac operon proteins is in the process, not in the products.

Transcriptional regulatory networks allow bacteria to express proteins only when they are needed. Adaptive hypotheses explaining the evolution of regulatory networks assume that unneeded expression is costly and therefore decreases fitness, but the proximate cause of this cost is not clear. We show that the cost in fitness to Escherichia coli strains constitutively expressing the lactose operon when lactose is absent is associated with the process of making the lac gene products, i.e., associated with the acts of transcription and/or translation. These results reject the hypotheses that regulation exists to prevent the waste of amino acids in useless protein or the detrimental activity of unnecessary proteins. While the cost of the process of protein expression occurs in all of the environments that we tested, the expression of the lactose permease could be costly or beneficial, depending on the environment. Our results identify the basis of a single selective pressure likely acting across the entire E. coli transcriptome.

Conspicuous impacts of inconspicuous hosts on the Lyme disease epidemic.

Emerging zoonotic pathogens are a constant threat to human health throughout the world. Control strategies to protect public health regularly fail, due in part to the tendency to focus on a single host species assumed to be the primary reservoir for a pathogen. Here, we present evidence that a diverse set of species can play an important role in determining disease risk to humans using Lyme disease as a model. Host-targeted public health strategies to control the Lyme disease epidemic in North America have focused on interrupting Borrelia burgdorferi sensu stricto (ss) transmission between blacklegged ticks and the putative dominant reservoir species, white-footed mice. However, B. burgdorferi ss infects more than a dozen vertebrate species, any of which could transmit the pathogen to feeding ticks and increase the density of infected ticks and Lyme disease risk. Using genetic and ecological data, we demonstrate that mice are neither the primary host for ticks nor the primary reservoir for B. burgdorferi ss, feeding 10% of all ticks and 25% of B. burgdorferi-infected ticks. Inconspicuous shrews feed 35% of all ticks and 55% of infected ticks. Because several important host species influence Lyme disease risk, interventions directed at a multiple host species will be required to control this epidemic.

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).

Selection for functional diversity drives accumulation of point mutations in Dr adhesins of Escherichia coli.

Immune escape is considered to be the driving force behind structural variability of major antigens on the surface of bacterial pathogens, such as fimbriae. In the Dr family of Escherichia coli adhesins, structural and adhesive functions are carried out by the same subunit. Dr adhesins have been shown to bind decay-accelerating factor (DAF), collagen IV, and carcinoembryonic antigen-related cell adhesion molecules (CEACAMs). We show that genes encoding Dr adhesins from 100 E. coli strains form eight structural groups with a high level of amino acid sequence diversity between them. However, genes comprising each group differ from each other by only a small number of point mutations. Out of 66 polymorphisms identified within the groups, only three were synonymous mutations, indicating strong positive selection for amino acid replacements. Functional analysis of intragroup variants comprising the Dr haemagglutinin (DraE) group revealed that the point mutations result in distinctly different binding phenotypes, with a tendency of increased affinity to DAF, decreased sensitivity of DAF binding to inhibition by chloramphenicol, and loss of binding capability to collagen, CEACAM3 and CEACAM6. Thus, variability by point mutation of major antigenic proteins on the bacterial surface can be a signature of selection for functional modification.

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

FimH, the mannose-specific, type 1 fimbrial adhesin of Escherichia coli, acquires amino acid replacements adaptive in extraintestinal niches (the genitourinary tract) but detrimental in the main habitat (the large intestine). This microevolutionary dynamics is reminiscent of an ecological "source-sink" model of continuous species spread from a stable primary habitat (source) into transient secondary niches (sink), with eventual extinction of the sink-evolved populations. Here, we have adapted two ecological analytical tools-diversity indexes DS and alpha--to compare size and frequency distributions of fimH haplotypes between evolutionarily conserved FimH variants ("source" haplotypes) and FimH variants with adaptive mutations (putative "sink" haplotypes). Both indexes show two- to threefold increased diversity of the sink fimH haplotypes relative to the source haplotypes, a pattern that ran opposite to those seen with nonstructural fimbrial genes (fimC and fimI) and housekeeping loci (adk and fumC) but similar to that seen with another fimbrial adhesin of E. coli, papG-II, also implicated in extraintestinal infections. The increased diversity of the sink pool of adhesin genes is due to the increased richness of the haplotypes (the number of unique haplotypes), rather than their evenness (the extent of similarity in relative abundances). Taken together, this pattern supports a continuous emergence and extinction of the gene alleles adaptive to virulence sink habitats of E. coli, rather than a one-time change in the habitat conditions. Thus, ecological methods of species diversity analysis can be successfully adapted to characterize the emergence of microbial virulence in bacterial pathogens subject to source-sink dynamics.

Source-sink dynamics of virulence evolution.

To understand the evolution of genetic diversity within species--bacterial and others--we must dissect the first steps of genetic adaptation to novel habitats, particularly habitats that are suboptimal for sustained growth where there is strong selection for adaptive changes. Here, we present the view that bacterial human pathogens represent an excellent model for understanding the molecular mechanisms of the adaptation of a species to alternative habitats. In particular, bacterial pathogens allow us to develop analytical methods to detect genetic adaptation using an evolutionary 'source-sink' model, with which the evolution of bacterial pathogens can be seen from the angle of continuous switching between permanent (source) and transient (sink) habitats. The source-sink model provides a conceptual framework for understanding the population dynamics and molecular mechanisms of virulence evolution.

A modest model explains the distribution and abundance of Borrelia burgdorferi strains.

The distribution and abundance of Borrelia burgdorferi, including human Lyme disease strains, is a function of its interactions with vertebrate species. We present a mathematical model describing important ecologic interactions affecting the distribution and abundance of B. burgdorferi strains, marked by the allele at the outer surface protein C locus, in Ixodes scapularis ticks, the principal vector. The frequency of each strain in ticks can be explained by the vertebrate species composition, the density of each vertebrate species, the number of ticks that feed on individuals of each species, and the rate at which those ticks acquire different strains. The model results are consistent with empirical data collected in a major Lyme disease focus in New England. An applicable extension of these results would be to predict the proportion of ticks carrying human infectious strains of B. burgdorferi from disease host densities and thus predict the local risk of contracting Lyme disease.

ospC diversity in Borrelia burgdorferi: different hosts are different niches.

The outer surface protein C (ospC) locus of the Lyme disease bacterium, Borrelia burgdorferi, is at least an order of magnitude more variable than other genes in the species. This variation is classified into 22 ospC major groups, 15 of which are found in the northeastern United States. The frequency distributions of ospC within populations suggest that this locus is under balancing selection. In multiple-niche polymorphism, a type of balancing selection, diversity within a population can be maintained when the environment is heterogeneous and no one genotype has the highest fitness in all environments. Genetically different individuals within vertebrate species and different vertebrate species constitute diverse environments for B. burgdorferi. We examined four important host species of B. burgdorferi and found that the strains that infected each species had different sets of ospC major groups. We found no variation among conspecific hosts in the ospC major groups of their infecting strains. These results suggest multiple niches create balancing selection at the ospC locus.

Evolution of specialists in an experimental microcosm.

The impact of adaptation on the persistence of a balanced polymorphism was explored using the lactose operon of Escherichia coli as a model system. Competition in chemostats for two substitutable resources, methylgalactoside and lactulose, generates stabilizing frequency-dependent selection when two different naturally isolated lac operons (TD2 and TD10) are used. The fate of this balanced polymorphism was tracked over evolutionary time by monitoring the frequency of fhuA-, a linked neutral genetic marker that confers resistance to the bacteriophage T5. In four out of nine chemostats the lac polymorphism persisted for 400-600 generations when the experiments were terminated. In the other five chemostats the fhuA polymorphism, and consequently the lac operon polymorphism, was lost between 86 and 219 generations. Four of 13 chemostat cultures monomorphic for the lac operon retained the neutral fhuA polymorphism for 450-550 generations until they were terminated; the remainder became monomorphic at fhuA between 63 and 303 generations. Specialists on each galactoside were isolated from chemostats that maintained the fhuA polymorphism, whether polymorphic or monomorphic at the lac operon. Strains isolated from three of four chemostats in which the lac polymorphism was preserved had switched their galactoside preference. Most of the chemostats where the fhuA polymorphism was lost also contained specialists. These results demonstrate that the initial polymorphism at lac was of little consequence to the outcome of long-term adaptive evolution. Instead, the fitnesses of evolved strains were dominated by mutations arising elsewhere in the genome, a fact confirmed by showing that operons isolated from their evolved backgrounds were alone unable to explain the presence of both specialists. Our results suggest that, once stabilized, ecological specialization prevented selective sweeps through the entire population, thereby promoting the maintenance of linked neutral polymorphisms.