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

Accelerated gene evolution through replication-transcription conflicts

Several mechanisms that increase the rate of mutagenesis across the entire genome have been identified; however, how the rate of evolution might be promoted in individual genes is unclear. Most genes in bacteria are encoded on the leading strand of replication. This presumably avoids the potentially detrimental head-on collisions that occur between the replication and transcription machineries when genes are encoded on the lagging strand. Here we identify the ubiquitous (core) genes in Bacillus subtilis and determine that 17% of them are on the lagging strand. We find a higher rate of point mutations in the core genes on the lagging strand compared with those on the leading strand, with this difference being primarily in the amino-acid-changing (nonsynonymous) mutations. We determine that, overall, the genes under strong negative selection against amino-acid-changing mutations tend to be on the leading strand, co-oriented with replication. In contrast, on the basis of the rate of convergent mutations, genes under positive selection for amino-acid-changing mutations are more commonly found on the lagging strand, indicating faster adaptive evolution in many genes in the head-on orientation. Increased gene length and gene expression amounts are positively correlated with the rate of accumulation of nonsynonymous mutations in the head-on genes, suggesting that the conflict between replication and transcription could be a driving force behind these mutations. Indeed, using reversion assays, we show that the difference in the rate of mutagenesis of genes in the two orientations is transcription dependent. Altogether, our findings indicate that head-on replication-transcription conflicts are more mutagenic than co-directional conflicts and that these encounters can significantly increase adaptive structural variation in the coded proteins. We propose that bacteria, and potentially other organisms, promote faster evolution of specific genes through orientation-dependent encounters between DNA replication and transcription.

Microbial variome database: point mutations, adaptive or not, in bacterial core genomes

Analysis of genetic differences (gene presence/absence and nucleotide polymorphisms) among strains of a bacterial species is crucial to understanding molecular mechanisms of bacterial pathogenesis and selecting targets for novel antibacterial therapeutics. However, lack of genome-wide association studies on large and epidemiologically well-defined strain collections from the same species makes it difficult to identify the genes under positive selection and define adaptive polymorphisms in those genes. To address this need and to overcome existing limitations, we propose to create a "microbial variome"--a species-specific resource database of genomic variations based on molecular evolutionary analysis. Here, we present prototype variome databases of Escherichia coli and Salmonella enterica subspecies enterica (http://depts.washington.edu/sokurel/variome, last accessed March 26, 2013). The prototypes currently include the point mutations data of core protein-coding genes from completely sequenced genomes of 22 E. coli and 17 S. enterica strains. These publicly available databases allow for single- and multiple-field sorting, filtering, and searching of the gene variability data and the potential adaptive significance. Such resource databases would immensely help experimental research, clinical diagnostics, epidemiology, and environmental control of human pathogens.

Draft Genome sequence of Escherichia coli AI27, a porcine isolate belonging to phylogenetic group B1

Escherichia coli AI27 is a putatively commensal strain isolated from feces of a pig. Here we report the draft genome sequence of E. coli AI27. This is the first porcine strain in the phylogenetic group B1 whose genome sequence has been determined.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

Estimating variation within the genes and inferring the phylogeny of 186 sequenced diverse Escherichia coli genomes

Background

Escherichia coli exists in commensal and pathogenic forms. By measuring the variation of individual genes across more than a hundred sequenced genomes, gene variation can be studied in detail, including the number of mutations found for any given gene. This knowledge will be useful for creating better phylogenies, for determination of molecular clocks and for improved typing techniques.

Results

We find 3,051 gene clusters/families present in at least 95% of the genomes and 1,702 gene clusters present in 100% of the genomes. The former 'soft core' of about 3,000 gene families is perhaps more biologically relevant, especially considering that many of these genome sequences are draft quality. The E. coli pan-genome for this set of isolates contains 16,373 gene clusters.

A core-gene tree, based on alignment and a pan-genome tree based on gene presence/absence, maps the relatedness of the 186 sequenced E. coli genomes. The core-gene tree displays high confidence and divides the E. coli strains into the observed MLST type clades and also separates defined phylotypes.

Conclusion

The results of comparing a large and diverse E. coli dataset support the theory that reliable and good resolution phylogenies can be inferred from the core-genome. The results further suggest that the resolution at the isolate level may, subsequently be improved by targeting more variable genes. The use of whole genome sequencing will make it possible to eliminate, or at least reduce, the need for several typing steps used in traditional epidemiology.

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.

Gene transfer agents: phage-like elements of genetic exchange

Horizontal gene transfer is important in the evolution of bacterial and archaeal genomes. An interesting genetic exchange process is carried out by diverse phage-like gene transfer agents (GTAs) that are found in a wide range of prokaryotes. Although GTAs resemble phages, they lack the hallmark capabilities that define typical phages, and they package random pieces of the producing cell's genome. In this Review, we discuss the defining characteristics of the GTAs that have been identified to date, along with potential functions for these agents and the possible evolutionary forces that act on the genes involved in their production.

Identification of hot and cold spots in genome of Mycobacterium tuberculosis using Shewhart Control Charts

The organization of genomic sequences is dynamic and undergoes change during the process of evolution. Many of the variations arise spontaneously and the observed genomic changes can either be distributed uniformly throughout the genome or be preferentially localized to some regions (hot spots) compared to others. Conversely cold spots may tend to accumulate very few variations or none at all. In order to identify such regions statistically, we have developed a method based on Shewhart Control Chart. The method was used for identification of hot and cold spots of single-nucleotide variations (SNVs) in Mycobacterium tuberculosis genomes. The predictions have been validated by sequencing some of these regions derived from clinical isolates. This method can be used for analysis of other genome sequences particularly infectious microbes.

Microevolution in fimH gene of mucosa-associated Escherichia coli strains isolated from pediatric patients with inflammatory bowel disease

Several studies reported increased numbers of mucosa-associated Escherichia coli strains in patients with inflammatory bowel disease (IBD), encompassing Crohn's disease (CD) and ulcerative colitis (UC). The majority of E. coli strains possess type 1 fimbriae, whose tip fibrillum protein, FimH, naturally undergoes amino acid replacements, an important process in the adaptation of commensal E. coli strains to environmental changes, like those observed in IBD and urinary tract infections. In this study, we analyzed mutational patterns in the fimH gene of 52 mucosa-associated E. coli strains isolated from IBD and non-IBD pediatric patients, in order to investigate microevolution of this genetic trait. FimH-positive strains were also phylogenetically typed and tested for their adhesive ability on Caco-2 cells. Specific FimH alleles for each grouping feature were found. Mutations G66S and V27A were related to CD, while mutations A242V, V163A, and T74I were attributed to UC. Otherwise, the G66S, N70S, and S78N mutations were specifically attributed to B2/D phylogroups. The N70S and A119V mutations were related to adhesive E. coli strains. Phylogroup B2, adhesive, and IBD E. coli strains showed a higher site substitution rate (SSR) in the fimH gene, together with a higher number of mutations. The degree of naïve mucosal inflammation was related to specific FimH alleles. Moreover, we could suggest that the V27A mutation is pathoadaptive for the CD intestinal habitat, while we could also suggest that both the N70S and S78N mutations are related to the more virulent E. coli B2 phylogroup. In conclusion, we found some FimH variants that seem to be more involved than others in the evolution of IBD pathogenesis.

Toward an efficient method of identifying core genes for evolutionary and functional microbial phylogenies

Microbial community metagenomes and individual microbial genomes are becoming increasingly accessible by means of high-throughput sequencing. Assessing organismal membership within a community is typically performed using one or a few taxonomic marker genes such as the 16S rDNA, and these same genes are also employed to reconstruct molecular phylogenies. There is thus a growing need to bioinformatically catalog strongly conserved core genes that can serve as effective taxonomic markers, to assess the agreement among phylogenies generated from different core gene, and to characterize the biological functions enriched within core genes and thus conserved throughout large microbial clades. We present a method to recursively identify core genes (i.e. genes ubiquitous within a microbial clade) in high-throughput from a large number of complete input genomes. We analyzed over 1,100 genomes to produce core gene sets spanning 2,861 bacterial and archaeal clades, ranging in size from one to >2,000 genes in inverse correlation with the α-diversity (total phylogenetic branch length) spanned by each clade. These cores are enriched as expected for housekeeping functions including translation, transcription, and replication, in addition to significant representations of regulatory, chaperone, and conserved uncharacterized proteins. In agreement with previous manually curated core gene sets, phylogenies constructed from one or more of these core genes agree with those built using 16S rDNA sequence similarity, suggesting that systematic core gene selection can be used to optimize both comparative genomics and determination of microbial community structure. Finally, we examine functional phylogenies constructed by clustering genomes by the presence or absence of orthologous gene families and show that they provide an informative complement to standard sequence-based molecular phylogenies.

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

Single nucleotide polypmorphisms of fimH associated with adherence and biofilm formation by serovars of Salmonella enterica

Type 1 fimbriae produced by serovars of Salmonella are characterized by their ability to agglutinate guinea pig erythrocytes in the absence of d-mannose but not in its presence. The FimH protein is the adhesin that mediates this reaction; it is distinct from the major fimbrial protei.n (FimA) that composes the fimbrial shaft. Avian-adapted serovars of Salmonella produce non-haemagglutinating fimbriae that have been reported to mediate adherence to avian cells. A single amino acid substitution is present in the FimH adhesin of these strains compared to that of a Typhimurium isolate. Also, previous studies have shown that single nucleotide polymorphisms in two strains of the Typhimurium fimH alter the binding specificity. We therefore investigated the allelic variation of fimH from a range of serotypes (both host-adapted and non-host-adapted) and isolates of Salmonella. Most FimH adhesins mediated the mannose-sensitive haemagglutination of guinea pig erythrocytes, but many did not facilitate adherence to HEp-2 cells. A small number of isolates also produced fimbriae but did not mediate adherence to either cell type. Transformants possessing cloned fimH genes exhibited a number of different substitutions within the predicted amino acid sequence of the FimH polypeptide. No identical FimH amino sequence was found between strains that adhere to erythrocytes and/or HEp-2 cells and those produced by non-adherent strains. FimH-mediated adherence to HEp-2 cells was invariably associated with the ability to form biofilms on mannosylated bovine serum albumin.

Functional and metabolic effects of adaptive glycerol kinase (GLPK) mutants in Escherichia coli

Herein we measure the effect of four adaptive non-synonymous mutations to the glycerol kinase (glpK) gene on catalytic function and regulation, to identify changes that correlate to increased fitness in glycerol media. The mutations significantly reduce affinity for the allosteric inhibitor fructose-1,6-bisphosphate (FBP) and formation of the tetramer, which are structurally related, in a manner that correlates inversely with imparted fitness during growth on glycerol, which strongly suggests that these enzymatic parameters drive growth improvement. Counterintuitively, the glpK mutations also increase glycerol-induced auto-catabolite repression that reduces glpK transcription in a manner that correlates to fitness. This suggests that increased specific GlpK activity is attenuated by negative feedback on glpK expression via catabolite repression, possibly to prevent methylglyoxal toxicity. We additionally report that glpK mutations were fixed in 47 of 50 independent glycerol-adapted lineages. By far the most frequently mutated locus (nucleotide 218) was mutated in 20 lineages, strongly suggesting this position has an elevated mutation rate. This study demonstrates that fitness correlations can be used to interrogate adaptive processes at the protein level and to identify the regulatory constraints underlying selection and improved growth.

Dynamical phase transition in a model for evolution with migration

We study a simple quasispecies model for evolution in two different habitats, with different fitness landscapes, coupled through one-way migration. Our key finding is a dynamical phase transition at a critical value of the migration rate, at which the time to reach the steady state diverges. The genetic composition of the population is qualitatively different above and below the transition. Using results from localization theory, we show that the critical migration rate may be very small-demonstrating that evolutionary outcomes can be very sensitive to even a small amount of migration.

Microevolution of Helicobacter pylori during prolonged infection of single hosts and within families

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

Mutation rates in bacteria have generally been considered to be much slower than in viruses. This is partly because estimates of long-term mutation rates for the evolution of distinct species have been inappropriately used for dating divergence within species. Furthermore, the most commonly used long-term mutation rate is based on geological dates that are no longer accepted. In addition, only few short-term mutation rates have been calculated within bacterial species, and these differ with the species by several orders of magnitude. Here, we provide robust estimates for short-term mutation and recombination rates within Helicobacter pylori, a bacterium that commonly infects the human gastric mucosa, based on serial isolates from long-term infections and on differences between isolates from multiple family members. These short-term mutation rates are 5–17-fold faster than long-term mutation rates in H. pylori that have been calibrated by parallel ancient migrations of humans. Short-term mutation rates in bacteria, including those for H. pylori, can be quite fast, partially overlapping with those for viruses. Future calculations of ages of bacterial species will need to account for dramatic differences in mutation rate between species and for dramatic differences between short- and long-term mutation rates.

Small variable segments constitute a major type of diversity of bacterial genomes at the species level

BACKGROUND

Analysis of large scale diversity in bacterial genomes has mainly focused on elements such as pathogenicity islands, or more generally, genomic islands. These comprise numerous genes and confer important phenotypes, which are present or absent depending on strains. We report that despite this widely accepted notion, most diversity at the species level is composed of much smaller DNA segments, 20 to 500 bp in size, which we call microdiversity.

RESULTS

We performed a systematic analysis of the variable segments detected by multiple whole genome alignments at the DNA level on three species for which the greatest number of genomes have been sequenced: Escherichia coli, Staphylococcus aureus, and Streptococcus pyogenes. Among the numerous sites of variability, 62 to 73% were loci of microdiversity, many of which were located within genes. They contribute to phenotypic variations, as 3 to 6% of all genes harbor microdiversity, and 1 to 9% of total genes are located downstream from a microdiversity locus. Microdiversity loci are particularly abundant in genes encoding membrane proteins. In-depth analysis of the E. coli alignments shows that most of the diversity does not correspond to known mobile or repeated elements, and it is likely that they were generated by illegitimate recombination. An intriguing class of microdiversity includes small blocks of highly diverged sequences, whose origin is discussed.

CONCLUSIONS

This analysis uncovers the importance of this small-sized genome diversity, which we expect to be present in a wide range of bacteria, and possibly also in many eukaryotic genomes.

Genomewide screening for novel genetic variations associated with ciprofloxacin resistance in Bacillus anthracis

Fluoroquinolone (FQ) resistance of Bacillus anthracis is a serious concern in the fields of biodefense and bioterrorism since FQs are very effective antibiotics and are recommended as first-line treatment against this lethal bacterium. In this study, we obtained 2 strains of B. anthracis showing resistance or intermediate resistance to ciprofloxacin (CIP) by a stepwise selection procedure with increasing CIP concentrations. Fifteen genetic variations were identified between the parental and CIP-resistant strains by next-generation sequencing. Nonsynonymous mutations in the quinolone resistance-determining region (QRDR) of type II DNA topoisomerase were identified in the resistant strain but not in the intermediate-resistant strain. The GBAA0834 (TetR-type transcriptional regulator) locus was also revealed to be a novel "mutation hot spot" that leads to the increased expression of multidrug efflux systems for CIP resistance. As an initial step of CIP resistance in B. anthracis, such disruptive mutations of GBAA0834 appear to be more easily acquired than those in an essential gene, such as that encoding type II DNA topoisomerase. Such an intermediate-resistant phenotype could increase a cell population under CIP-selective pressure and might promote the emergence of highly resistant isolates. Our findings reveal, in addition to QRDR, crucial genetic targets for the investigation of intermediate resistance of B. anthracis to FQs.

Comparative structure-function analysis of mannose-specific FimH adhesins from Klebsiella pneumoniae and Escherichia coli

FimH, the adhesive subunit of type 1 fimbriae expressed by many enterobacteria, mediates mannose-sensitive binding to target host cells. At the same time, fine receptor-structural specificities of FimH from different species can be substantially different, affecting bacterial tissue tropism and, as a result, the role of the particular fimbriae in pathogenesis. In this study, we compared functional properties of the FimH proteins from Escherichia coli and Klebsiella pneumoniae, which are both 279 amino acids in length but differ by some approximately 15% of residues. We show that K. pneumoniae FimH is unable to mediate adhesion in a monomannose-specific manner via terminally exposed Manalpha(1-2) residues in N-linked oligosaccharides, which are the structural basis of the tropism of E. coli FimH for uroepithelial cells. However, K. pneumoniae FimH can bind to the terminally exposed Manalpha(1-3)Manbeta(1-4)GlcNAcbeta1 trisaccharide, though only in a shear-dependent manner, wherein the binding is marginal at low shear force but enhanced sevenfold under increased shear. A single mutation in the K. pneumoniae FimH, S62A, converts the mode of binding from shear dependent to shear independent. This mutation has occurred naturally in the course of endemic circulation of a nosocomial uropathogenic clone and is identical to a pathogenicity-adaptive mutation found in highly virulent uropathogenic strains of E. coli, in which it also eliminates the dependence of E. coli binding on shear. The shear-dependent binding properties of the K. pneumoniae and E. coli FimH proteins are mediated via an allosteric catch bond mechanism. Thus, despite differences in FimH structure and fine receptor specificity, the shear-dependent nature of FimH-mediated adhesion is highly conserved between bacterial species, supporting its remarkable physiological significance.