Effect of chromosome location on bacterial mutation rates

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Ciprofloxacin induced antibiotic resistance in Salmonella Typhimurium mutants and genome analysis

Antibiotic resistance of Salmonella species is well reported. Ciprofloxacin is the frontline antibiotic for salmonellosis. The repeated exposure to ciprofloxacin leads to resistant strains. After 20 cycles of antibiotic exposure, resistant bacterial clones were evaluated. The colony size of the mutants was small and had an extended lag phase compared to parent strain. The whole genome sequencing showed 40,513 mutations across the genome. Small percentage (5.2%) of mutations was non-synonymous. Four-fold more transitions were observed than transversions. Ratio of < 1 transition vs transversion showed a positive selection for antibiotic resistant trait. Mutation distribution across the genome was uniform. The native plasmid was an exception and 2 mutations were observed on 90 kb plasmid. The important genes like dnaE, gyrA, iroC, metH and rpoB involved in antibiotic resistance had point mutations. The genome analysis revealed most of the metabolic pathways were affected.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

The Location of Substitutions and Bacterial Genome Arrangements

Increasing evidence supports the notion that different regions of a genome have unique rates of molecular change. This variation is particularly evident in bacterial genomes where previous studies have reported gene expression and essentiality tend to decrease, whereas substitution rates usually increase with increasing distance from the origin of replication. Genomic reorganization such as rearrangements occur frequently in bacteria and allow for the introduction and restructuring of genetic content, creating gradients of molecular traits along genomes. Here, we explore the interplay of these phenomena by mapping substitutions to the genomes of Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti, quantifying how many substitutions have occurred at each position in the genome. Preceding work indicates that substitution rate significantly increases with distance from the origin. Using a larger sample size and accounting for genome rearrangements through ancestral reconstruction, our analysis demonstrates that the correlation between the number of substitutions and the distance from the origin of replication is significant but small and inconsistent in direction. Some replicons had a significantly decreasing trend (E. coli and the chromosome of S. meliloti), whereas others showed the opposite significant trend (B. subtilis, Streptomyces, pSymA and pSymB in S. meliloti). dN, dS, and ω were examined across all genes and there was no significant correlation between those values and distance from the origin. This study highlights the impact that genomic rearrangements and location have on molecular trends in some bacteria, illustrating the importance of considering spatial trends in molecular evolutionary analysis. Assuming that molecular trends are exclusively in one direction can be problematic.

Mutation and Recombination Rates Vary Across Bacterial Chromosome

Bacteria evolve as a result of mutations and acquisition of foreign DNA by recombination processes. A growing body of evidence suggests that mutation and recombination rates are not constant across the bacterial chromosome. Bacterial chromosomal DNA is organized into a compact nucleoid structure which is established by binding of the nucleoid-associated proteins (NAPs) and other proteins. This review gives an overview of recent findings indicating that the mutagenic and recombination processes in bacteria vary at different chromosomal positions. Involvement of NAPs and other possible mechanisms in these regional differences are discussed. Variations in mutation and recombination rates across the bacterial chromosome may have implications in the evolution of bacteria.

Comparative pan-genomic analyses of Orientia tsutsugamushi reveal an exceptional model of bacterial evolution driving genomic diversity

Orientia tsutsugamushi, formerly Rickettsia tsutsugamushi, is an obligate intracellular pathogen that causes scrub typhus, an underdiagnosed acute febrile disease with high morbidity. Scrub typhus is transmitted by the larval stage (chigger) of Leptotrombidium mites and is irregularly distributed across endemic regions of Asia, Australia and islands of the western Pacific Ocean. Previous work to understand population genetics in O. tsutsugamushi has been based on sub-genomic sampling methods and whole-genome characterization of two genomes. In this study, we compared 40 genomes from geographically dispersed areas and confirmed patterns of extensive homologous recombination likely driven by transposons, conjugative elements and repetitive sequences. High rates of lateral gene transfer (LGT) among O. tsutsugamushi genomes appear to have effectively eliminated a detectable clonal frame, but not our ability to infer evolutionary relationships and phylogeographical clustering. Pan-genomic comparisons using 31 082 high-quality bacterial genomes from 253 species suggests that genomic duplication in O. tsutsugamushi is almost unparalleled. Unlike other highly recombinant species where the uptake of exogenous DNA largely drives genomic diversity, the pan-genome of O. tsutsugamushi is driven by duplication and divergence. Extensive gene innovation by duplication is most commonly attributed to plants and animals and, in contrast with LGT, is thought to be only a minor evolutionary mechanism for bacteria. The near unprecedented evolutionary characteristics of O. tsutsugamushi, coupled with extensive intra-specific LGT, expand our present understanding of rapid bacterial evolutionary adaptive mechanisms.

Broad-scale redistribution of mRNA abundance and transcriptional machinery in response to growth rate in Salmonella enterica serovar Typhimurium

We have investigated the connection between the four-dimensional architecture of the bacterial nucleoid and the organism's global gene expression programme. By localizing the transcription machinery and the transcriptional outputs across the genome of the model bacterium Salmonella enterica serovar Typhimurium at different stages of the growth cycle, a surprising disconnection between gene dosage and transcriptional output was revealed. During exponential growth, gene output occurred chiefly in the Ori (origin), Ter (terminus) and NSL (non-structured left) domains, whereas the Left macrodomain remained transcriptionally quiescent at all stages of growth. The apparently high transcriptional output in Ter was correlated with an enhanced stability of the RNA expressed there during exponential growth, suggesting that longer mRNA half-lives compensate for low gene dosage. During exponential growth, RNA polymerase (RNAP) was detected everywhere, whereas in stationary phase cells, RNAP was concentrated in the Ter macrodomain. The alternative sigma factors RpoE, RpoH and RpoN were not required to drive transcription in these growth conditions, consistent with their observed binding to regions away from RNAP and regions of active transcription. Specifically, these alternative sigma factors were found in the Ter macrodomain during exponential growth, whereas they were localized at the Ori macrodomain in stationary phase.

Complex chromosomal neighborhood effects determine the adaptive potential of a gene under selection

How the organization of genes on a chromosome shapes adaptation is essential for understanding evolutionary paths. Here, we investigate how adaptation to rapidly increasing levels of antibiotic depends on the chromosomal neighborhood of a drug-resistance gene inserted at different positions of the Escherichia coli chromosome. Using a dual-fluorescence reporter that allows us to distinguish gene amplifications from other up-mutations, we track in real-time adaptive changes in expression of the drug-resistance gene. We find that the relative contribution of several mutation types differs systematically between loci due to properties of neighboring genes: essentiality, expression, orientation, termination, and presence of duplicates. These properties determine rate and fitness effects of gene amplification, deletions, and mutations compromising transcriptional termination. Thus, the adaptive potential of a gene under selection is a system-property with a complex genetic basis that is specific for each chromosomal locus, and it can be inferred from detailed functional and genomic data.

DOI:http://dx.doi.org/10.7554/eLife.25100.001

The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa

Background

Natural genetic variation ultimately arises from the process of mutation. Knowledge of how the raw material for evolution is produced is necessary for a full understanding of several fundamental evolutionary concepts. We performed a mutation accumulation experiment with wild-type and mismatch-repair deficient, mutator lines of the pathogenic bacterium Pseudomonas aeruginosa, and used whole-genome sequencing to reveal the genome-wide rate, spectrum, distribution, leading/lagging bias, and context-dependency of spontaneous mutations.

Results

Wild-type base-pair mutation and indel rates were ~10⁻¹⁰ and ~10⁻¹¹ per nucleotide per generation, respectively, and deficiencies in the mismatch-repair system caused rates to increase by over two orders of magnitude. A universal bias towards AT was observed in wild-type lines, but was reversed in mutator lines to a bias towards GC. Biases for which types of mutations occurred during replication of the leading versus lagging strand were detected reciprocally in both replichores. The distribution of mutations along the chromosome was non-random, with peaks near the terminus of replication and at positions intermediate to the replication origin and terminus. A similar distribution bias was observed along the chromosome in natural populations of P. aeruginosa. Site-specific mutation rates were higher when the focal nucleotide was immediately flanked by C:G pairings.

Conclusions

Whole-genome sequencing of mutation accumulation lines allowed the comprehensive identification of mutations and revealed what factors of molecular and genomic architecture affect the mutational process. Our study provides a more complete view of how several mechanisms of mutation, mutation repair, and bias act simultaneously to produce the raw material for evolution.

Electronic supplementary material

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

Evolution of new functions de novo and from preexisting genes

How the enormous structural and functional diversity of new genes and proteins was generated (estimated to be 10(10)-10(12) different proteins in all organisms on earth [Choi I-G, Kim S-H. 2006. Evolution of protein structural classes and protein sequence families. Proc Natl Acad Sci 103: 14056-14061] is a central biological question that has a long and rich history. Extensive work during the last 80 years have shown that new genes that play important roles in lineage-specific phenotypes and adaptation can originate through a multitude of different mechanisms, including duplication, lateral gene transfer, gene fusion/fission, and de novo origination. In this review, we focus on two main processes as generators of new functions: evolution of new genes by duplication and divergence of pre-existing genes and de novo gene origination in which a whole protein-coding gene evolves from a noncoding sequence.

The Rate and Molecular Spectrum of Spontaneous Mutations in the GC-Rich Multichromosome Genome of Burkholderia cenocepacia

Spontaneous mutations are ultimately essential for evolutionary change and are also the root cause of many diseases. However, until recently, both biological and technical barriers have prevented detailed analyses of mutation profiles, constraining our understanding of the mutation process to a few model organisms and leaving major gaps in our understanding of the role of genome content and structure on mutation. Here, we present a genome-wide view of the molecular mutation spectrum in Burkholderia cenocepacia, a clinically relevant pathogen with high %GC content and multiple chromosomes. We find that B. cenocepacia has low genome-wide mutation rates with insertion-deletion mutations biased toward deletions, consistent with the idea that deletion pressure reduces prokaryotic genome sizes. Unlike prior studies of other organisms, mutations in B. cenocepacia are not AT biased, which suggests that at least some genomes with high %GC content experience unusual base-substitution mutation pressure. Importantly, we also observe variation in both the rates and spectra of mutations among chromosomes and elevated G:C > T:A transversions in late-replicating regions. Thus, although some patterns of mutation appear to be highly conserved across cellular life, others vary between species and even between chromosomes of the same species, potentially influencing the evolution of nucleotide composition and genome architecture.

Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins

The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates' recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins, implemented in Python and C and supported on Linux and Mac OS X.

Mutation frequency and spectrum of mutations vary at different chromosomal positions of Pseudomonas putida

It is still an open question whether mutation rate can vary across the bacterial chromosome. In this study, the occurrence of mutations within the same mutational target sequences at different chromosomal locations of Pseudomonas putida was monitored. For that purpose we constructed two mutation detection systems, one for monitoring the occurrence of a broad spectrum of mutations and transposition of IS element IS1411 inactivating LacI repressor, and another for detecting 1-bp deletions. Our results revealed that both the mutation frequency and the spectrum of mutations vary at different chromosomal positions. We observed higher mutation frequencies when the direction of transcription of the mutational target gene was opposite to the direction of replisome movement in the chromosome and vice versa, lower mutation frequency was accompanied with co-directional transcription and replication. Additionally, asymmetry of frameshift mutagenesis at homopolymeric and repetitive sequences during the leading and lagging-strand replication was found. The transposition frequency of IS1411 was also affected by the chromosomal location of the target site, which implies that regional differences in chromosomal topology may influence transposition of this mobile element. The occurrence of mutations in the P. putida chromosome was investigated both in growing and in stationary-phase bacteria. We found that the appearance of certain mutational hot spots is strongly affected by the chromosomal location of the mutational target sequence especially in growing bacteria. Also, artificial increasing transcription of the mutational target gene elevated the frequency of mutations in growing bacteria.

A specific mutation in the promoter region of the silent cel cluster accounts for the appearance of lactose-utilizing Lactococcus lactis MG1363

The Lactococcus lactis laboratory strain MG1363 has been described to be unable to utilize lactose. However, in a rich medium supplemented with lactose as the sole carbon source, it starts to grow after prolonged incubation periods. Transcriptome analyses showed that L. lactis MG1363 Lac(+) cells expressed celB, encoding a putative cellobiose-specific phosphotransferase system (PTS) IIC component, which is normally silent in MG1363 Lac(-) cells. Nucleotide sequence analysis of the cel cluster of a Lac(+) isolate revealed a change from one of the guanines to adenine in the promoter region. We showed here that one particular mutation, taking place at increased frequency, accounts for the lactose-utilizing phenotype occurring in MG1363 cultures. The G-to-A transition creates a -10 element at an optimal distance from the -35 element. Thus, a fully active promoter is created, allowing transcription of the otherwise cryptic cluster. Nuclear magnetic resonance (NMR) spectroscopy results show that MG1363 Lac(+) uses a novel pathway of lactose utilization.

Escherichia coli frameshift mutation rate depends on the chromosomal context but not on the GATC content near the mutation site

Different studies have suggested that mutation rate varies at different positions in the genome. In this work we analyzed if the chromosomal context and/or the presence of GATC sites can affect the frameshift mutation rate in the Escherichia coli genome. We show that in a mismatch repair deficient background, a condition where the mutation rate reflects the fidelity of the DNA polymerization process, the frameshift mutation rate could vary up to four times among different chromosomal contexts. Furthermore, the mismatch repair efficiency could vary up to eight times when compared at different chromosomal locations, indicating that detection and/or repair of frameshift events also depends on the chromosomal context. Also, GATC sequences have been proved to be essential for the correct functioning of the E. coli mismatch repair system. Using bacteriophage heteroduplexes molecules it has been shown that GATC influence the mismatch repair efficiency in a distance- and number-dependent manner, being almost nonfunctional when GATC sequences are located at 1 kb or more from the mutation site. Interestingly, we found that in E. coli genomic DNA the mismatch repair system can efficiently function even if the nearest GATC sequence is located more than 2 kb away from the mutation site. The results presented in this work show that even though frameshift mutations can be efficiently generated and/or repaired anywhere in the genome, these processes can be modulated by the chromosomal context that surrounds the mutation site.

Escape from growth restriction in small colony variants of Salmonella typhimurium by gene amplification and mutation

Antibiotic resistance in bacteria is generally associated with fitness costs that often can be reduced by second-site compensatory mutations. Here, we examined how a protamine-resistant small colony variant of Salmonella typhimurium adapts to the growth reduction conferred by a resistance mutation in hemC (encoding a haem-biosynthesis enzyme). We show that adaptation occurs in a multi-step process where fitness is successively increased. Thus, the initial adaptive response was selection for an unstable gene amplification of the mutant hemC gene that provided a small fitness increase. Fitness was increased further by a mutation that restored HemC function in one gene copy, relaxing selection for the amplification. Subsequently, the amplification segregated back to the haploid state and even higher fitness. The end result was in most cases mutant strains with a hemC sequence different from that of the wild-type strain. These findings suggest that gene amplification facilitates adaptive evolution. A higher gene dosage increases the target size for compensatory mutations and improves fitness of the cell, thereby allowing an increase in the population size, further increasing the probability of a subsequent stable mutation. Our results provide a novel genetic basis for growth compensation in small colony variants.

Short-term signatures of evolutionary change in the Salmonella enterica serovar typhimurium 14028 genome

Salmonella enterica serovar Typhimurium is a Gram-negative pathogen that causes gastroenteritis in humans and a typhoid-like disease in mice and is often used as a model for the disease promoted by the human-adapted S. enterica serovar Typhi. Despite its health importance, the only S. Typhimurium strain for which the complete genomic sequence has been determined is the avirulent LT2 strain, which is extensively used in genetic and physiologic studies. Here, we report the complete genomic sequence of the S. Typhimurium strain 14028s, as well as those of its progenitor and two additional derivatives. Comparison of these S. Typhimurium genomes revealed differences in the patterns of sequence evolution and the complete inventory of genetic alterations incurred in virulent and avirulent strains, as well as the sequence changes accumulated during laboratory passage of pathogenic organisms.

Adaptive divergence in experimental populations of Pseudomonas fluorescens. IV. Genetic constraints guide evolutionary trajectories in a parallel adaptive radiation

The capacity for phenotypic evolution is dependent upon complex webs of functional interactions that connect genotype and phenotype. Wrinkly spreader (WS) genotypes arise repeatedly during the course of a model Pseudomonas adaptive radiation. Previous work showed that the evolution of WS variation was explained in part by spontaneous mutations in wspF, a component of the Wsp-signaling module, but also drew attention to the existence of unknown mutational causes. Here, we identify two new mutational pathways (Aws and Mws) that allow realization of the WS phenotype: in common with the Wsp module these pathways contain a di-guanylate cyclase-encoding gene subject to negative regulation. Together, mutations in the Wsp, Aws, and Mws regulatory modules account for the spectrum of WS phenotype-generating mutations found among a collection of 26 spontaneously arising WS genotypes obtained from independent adaptive radiations. Despite a large number of potential mutational pathways, the repeated discovery of mutations in a small number of loci (parallel evolution) prompted the construction of an ancestral genotype devoid of known (Wsp, Aws, and Mws) regulatory modules to see whether the types derived from this genotype could converge upon the WS phenotype via a novel route. Such types-with equivalent fitness effects-did emerge, although they took significantly longer to do so. Together our data provide an explanation for why WS evolution follows a limited number of mutational pathways and show how genetic architecture can bias the molecular variation presented to selection.

Contribution of gene amplification to evolution of increased antibiotic resistance in Salmonella typhimurium

The use of beta-lactam antibiotics has led to the evolution and global spread of a variety of resistance mechanisms, including beta-lactamases, a group of enzymes that degrade the beta-lactam ring. The evolution of increased beta-lactam resistance was studied by exposing independent lineages of Salmonella typhimurium to progressive increases in cephalosporin concentration. Each lineage carried a beta-lactamase gene (bla(TEM-1)) that provided very low resistance. In most lineages, the initial response to selection was an amplification of the bla(TEM-1) gene copy number. Amplification was followed in some lineages by mutations (envZ, cpxA, or nmpC) that reduced expression of the uptake functions, the OmpC, OmpD, and OmpF porins. The initial resistance provided by bla(TEM-1) amplification allowed the population to expand sufficiently to realize rare secondary point mutations. Mathematical modeling showed that amplification often is likely to be the initial response because events that duplicate or further amplify a gene are much more frequent than point mutations. These models show the importance of the population size to appearance of later point mutations. Transient gene amplification is likely to be a common initial mechanism and an intermediate in stable adaptive improvement. If later point mutations (allowed by amplification) provide sufficient adaptive improvement, the amplification may be lost.

Horizontal gene transfer between Ralstonia solanacearum strains detected by comparative genomic hybridization on microarrays

The plant pathogenic Betaproteobacterium Ralstonia solanacearum is a complex species in that most of the strains share the common characteristic of being naturally transformable. In this study, we used a new approach based on comparative genomic hybridization (CGH) on microarrays to investigate the extent of horizontal gene transfers (HGTs) between different strains of R. solanacearum. Recipient strains from phylotypes I, II and III were naturally transformed in vitro by genomic DNA from the GMI1000 reference strain (phylotype I) and the resulting DNAs were hybridized on a microarray representative of the 5120 predicted genes from the GMI1000 strain. In addition to transfer of the antibiotic resistance marker, in 8 of the 16 tested transformants, CGH on microarrays detected other transferred GMI1000 genes and revealed their number, category, function and localization along the genome. We showed that DNA blocks up to 30 kb and 33 genes could be integrated during a single event. Most of these blocks flanked the marker gene DNA but, interestingly, multiple DNA acquisitions along the genome also occurred in a single recombinant clone in one transformation experiment. The results were confirmed by PCR amplification, cloning and sequencing and Southern blot hybridization. This represents the first comprehensive identification of gene acquisitions and losses along the genome of the recipient bacterial strain during natural transformation experiments. In future studies, this strategy should help to answer many questions related to HGT mechanisms.

Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the ∼18,000 families of orthologous genes, we found ∼2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genome's long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

Although abundant knowledge has been accumulated regarding the E. coli laboratory strain K-12, little is known about the evolutionary trajectories that have driven the high diversity observed among natural isolates of the species, which encompass both commensal and highly virulent intestinal and extraintestinal pathogenic strains. We have annotated or re-annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Although recombination rates are much higher than mutation rates, we were able to reconstruct a robust phylogeny based on the ∼2,000 genes common to all strains. Based on this phylogeny, we established the evolutionary scenario of gains and losses of thousands of specific genes, identifying functional classes under opposite selection pressures. This genome flux is confined to very few positions in the chromosome, which are the same for every genome. Notably, we identified few or no extraintestinal virulence-specific genes. We also defined a long-scale structure of recombination in the genome with lower recombination rates at the terminus of replication. These findings demonstrate that, despite a very high gene flow, genes can co-exist in an organised genome.

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.

Natural isolates of Salmonella enterica serovar Dublin carry a single nadA missense mutation

Nicotinic acid is required by most isolates of Salmonella enterica (serovar Dublin), a pathogen of cattle. A single nadA missense mutation causes the nutritional requirement of all serovar Dublin isolates tested. Models for persistence of this allele are tested and discussed.

Reduction of the fitness burden of quinolone resistance in Pseudomonas aeruginosa

OBJECTIVES

Quinolone resistance in the opportunistic pathogen Pseudomonas aeruginosa is commonly caused by mutations that alter the target molecules DNA gyrase/topoisomerase IV, or cause activation of various efflux systems. We have analysed the effect of quinolone resistance caused by DNA gyrase/topoisomerase IV mutations on bacterial fitness.

METHODS

Norfloxacin-resistant mutants were isolated and by DNA sequencing the mutations conferring resistance were identified. Mutant fitness was determined by measuring growth rates in vitro. Mutants with reduced growth rates were serially passaged to obtain growth-compensated mutants. The level of DNA supercoiling was determined by isolating plasmid DNA from the susceptible, resistant and compensated mutants and comparing the topoisomer distribution patterns by gel electrophoresis in the presence of chloroquine.

RESULTS

Low-level resistance (4-48 mg/L) was caused by single mutations in gyrA or gyrB. Among these strains, three out of eight mutants showed lower fitness, whereas high-level resistant (>256 mg/L) mutants with double mutations in gyrA and parC, parE, nfxB or unknown genes all showed a reduced fitness. Slow-growing resistant mutants with a gyrA mutation had decreased DNA supercoiling. After serial passage in laboratory medium, mutant fitness was increased by compensatory mutation(s) that restored supercoiling to normal levels. The compensatory mutation(s) was not located in any of the genes (gyrAB, topA, parCE, hupB, fis, hupN, himAD or PA5348) that were expected to affect supercoiling.

CONCLUSIONS

Our results show that 'no cost' and compensatory mutations are common in quinolone-resistant P. aeruginosa.

Experimental adaptation of Salmonella typhimurium to mice

Experimental evolution is a powerful approach to study the dynamics and mechanisms of bacterial niche specialization. By serial passage in mice, we evolved 18 independent lineages of Salmonella typhimurium LT2 and examined the rate and extent of adaptation to a mainly reticuloendothelial host environment. Bacterial mutation rates and population sizes were varied by using wild-type and DNA repair-defective mutator (mutS) strains with normal and high mutation rates, respectively, and by varying the number of bacteria intraperitoneally injected into mice. After <200 generations of adaptation all lineages showed an increased fitness as measured by a faster growth rate in mice (selection coefficients 0.11-0.58). Using a generally applicable mathematical model we calculated the adaptive mutation rate for the wild-type bacterium to be >10(-6)/cell/generation, suggesting that the majority of adaptive mutations are not simple point mutations. For the mutator lineages, adaptation to mice was associated with a loss of fitness in secondary environments as seen by a reduced metabolic capability. During adaptation there was no indication that a high mutation rate was counterselected. These data show that S. typhimurium can rapidly and extensively increase its fitness in mice but this niche specialization is, at least in mutators, associated with a cost.

The replication-related organization of bacterial genomes

The replication of the chromosome is among the most essential functions of the bacterial cell and influences many other cellular mechanisms, from gene expression to cell division. Yet the way it impacts on the bacterial chromosome was not fully acknowledged until the availability of complete genomes allowed one to look upon genomes as more than bags of genes. Chromosomal replication includes a set of asymmetric mechanisms, among which are a division in a lagging and a leading strand and a gradient between early and late replicating regions. These differences are the causes of many of the organizational features observed in bacterial genomes, in terms of both gene distribution and sequence composition along the chromosome. When asymmetries or gradients increase in some genomes, e.g. due to a different composition of the DNA polymerase or to a higher growth rate, so do the corresponding biases. As some of the features of the chromosome structure seem to be under strong selection, understanding such biases is important for the understanding of chromosome organization and adaptation. Inversely, understanding chromosome organization may shed further light on questions relating to replication and cell division. Ultimately, the understanding of the interplay between these different elements will allow a better understanding of bacterial genetics and evolution.

Transcription increases multiple spontaneous point mutations in Salmonella enterica

The spontaneous rate of G.C-->A.T mutations and a hotspot T.A-->G.C transversion are known to increase with the frequency of transcription-increases that have been ascribed primarily to processes that affect only these specific mutations. To investigate how transcription induces other spontaneous point mutations, we tested for its effects in repair-proficient Salmonella enterica using reversion assays of chromosomally inserted alleles. Our results indicate that transcription increases rates of all tested point mutations in the induced gene: induction significantly increased the individual rates of an A.T-->T.A transversion, an A.T-->G.C transition and the pooled rates of the three other point mutations assayed. Although the S.enterica genome is thought to have a mutational bias towards G.C base pairs, transitions creating A.T pairs were approximately 10 times more frequent than the reverse mutation, resulting in an overall mutation pressure to lower G+C contents. Transitions occurred at roughly twice the rate of transversions, similar to results from sequence comparisons; however, several individual transversions are more frequent than the least common transition.

The YggX protein of Salmonella enterica is involved in Fe(II) trafficking and minimizes the DNA damage caused by hydroxyl radicals: residue CYS-7 is essential for YggX function

Previous work from our laboratory identified YggX as a protein whose accumulation increased the resistance of Salmonella enterica to superoxide stress, reversed defects attributed to oxidized [Fe-S] clusters, and decreased the spontaneous mutation frequency of the cells. Here we present work aimed at determining why the accumulation of YggX correlates with reduced mutation frequency. Genetic and biochemical data showed that accumulation of YggX reduced the damage to DNA by hydroxyl radicals. The ability of purified YggX to protect DNA from Fenton chemistry mediated damage in vitro and to decrease the concentration of Fe(II) ions in solution available for chelation provided a framework for the interpretation of data obtained from in vivo experiments. The interpretation of in vitro assay results, within the context of the in vivo phenotypes, was validated by a mutant variant of YggX (C7S) that was unable to function in vivo or in vitro. We propose a model, based on data presented here and reported earlier, that suggests YggX is a player in Fe(II) trafficking in bacteria.

Compensatory adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium

Most chromosomal mutations that cause antibiotic resistance impose fitness costs on the bacteria. This biological cost can often be reduced by compensatory mutations. In Salmonella typhimurium, the nucleotide substitution AAA42 --> AAC in the rpsL gene confers resistance to streptomycin. The resulting amino acid substitution (K42N) in ribosomal protein S12 causes an increased rate of ribosomal proofreading and, as a result, the rate of protein synthesis, bacterial growth and virulence are decreased. Eighty-one independent lineages of the low-fitness, K42N mutant were evolved in the absence of antibiotic to ameliorate the costs. From the rate of fixation of compensated mutants and their fitness, the rate of compensatory mutations was estimated to be > or = 10-7 per cell per generation. The size of the population bottleneck during evolution affected fitness of the adapted mutants: a larger bottleneck resulted in higher average fitness. Only four of the evolved lineages contained streptomycin-sensitive revertants. The remaining 77 lineages contained mutants that were still fully streptomycin resistant, had retained the original resistance mutation and also acquired compensatory mutations. Most of the compensatory mutations, resulting in at least 35 different amino acid substitutions, were novel single-nucleotide substitutions in the rpsD, rpsE, rpsL or rplS genes encoding the ribosomal proteins S4, S5, S12 and L19 respectively. Our results show that the deleterious effects of a resistance mutation can be compensated by an unexpected variety of mutations.

Gene location and bacterial sequence divergence

Previous comparison of a relatively small set of homologous genes from Escherichia coli and Salmonella typhimurium revealed that genes nearer to the origin of replication had substitution rates lower than genes closer to the replication terminus. The recently completed sequences of numerous bacterial genomes have allowed us to test whether this effect of distance from the replication origin on substitution rates, as observed for the E. coli-S. typhimurium comparison, is a general feature of bacterial genomes. Extending the analysis to all 3,000 E. coli-S. typhimurium homologs confirmed the significant association between chromosomal position and synonymous site divergence. However, the effect, though still significant, is not as dramatic as originally thought. A similar association between relative chromosomal location and synonymous substitution rate was detected in the majority of other bacterial species comparisons within alpha- and gamma- Proteobacteria, and Firmicutes but was absent in Chlamydiales. The opposite trend, i.e., a decrease in synonymous divergence with distance from the replication origin, was detected in Mycobacteria. Analysis of the patterns of nucleotide substitutions revealed that the distance effect is not affected by gene orientation and is mainly caused by an increase in rates of transversions, suggesting that this effect may not be caused by recombinational repair or biased gene conversion, as originally suggested.