Evolutionary implications of the frequent horizontal transfer of mismatch repair genes

Shifts in mutation spectra enhance access to beneficial mutations

Biased mutation spectra are pervasive, with wide variation in the magnitude of mutational biases that influence genome evolution and adaptation. How do such diverse biases evolve? Our experiments show that changing the mutation spectrum allows populations to sample previously undersampled mutational space, including beneficial mutations. The resulting shift in the distribution of fitness effects is advantageous: Beneficial mutation supply and beneficial pleiotropy both increase, while deleterious load reduces. More broadly, simulations indicate that reducing or reversing the direction of a long-term bias is always selectively favored. Such changes in mutation bias can occur easily via altered function of DNA repair genes. A phylogenetic analysis shows that these genes are repeatedly gained and lost in bacterial lineages, leading to frequent bias shifts in opposite directions. Thus, shifts in mutation spectra may evolve under selection and can directly alter the outcome of adaptive evolution by facilitating access to beneficial mutations.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

Antibiotic resistance in plant growth promoting bacteria: A comprehensive review and future perspectives to mitigate potential gene invasion risks

Plant growth-promoting bacteria (PGPB) are endowed with several attributes that can be beneficial for host plants. They opened myriad doors toward green technology approach to reduce the use of chemical inputs, improve soil fertility, and promote plants' health. However, many of these PGPB harbor antibiotic resistance genes (ARGs). Less attention has been given to multi-resistant bacterial bioinoculants which may transfer their ARGs to native soil microbial communities and other environmental reservoirs including animals, waters, and humans. Therefore, large-scale inoculation of crops by ARGs-harboring bacteria could worsen the evolution and dissemination of antibiotic resistance and aggravate the negative impacts on such ecosystem and ultimately public health. Their introduction into the soil could serve as ARGs invasion which may inter into the food chain. In this review, we underscore the antibiotic resistance of plant-associated bacteria, criticize the lack of consideration for this phenomenon in the screening and application processes, and provide some recommendations as well as a regulation framework relating to the development of bacteria-based biofertilizers to aid maximizing their value and applications in crop improvement while reducing the risks of ARGs invasion.

Speciation of Genes and Genomes: Conservation of DNA Polymorphism by Barriers to Recombination Raised by Mismatch Repair System

Some basic aspects of human and animal biology and evolution involve the establishment of biological uniqueness of species and individuals within their huge variety. The discrimination among closely related species occurs in their offspring at the level of chromosomal DNA sequence homology, which is required for fertility as the hallmark of species. Biological identification of individuals, i.e., of their biological “self”, occurs at the level of protein sequences presented by the MHC/HLA complex as part of the immune system that discriminates non-self from self. Here, a mechanistic molecular model is presented that can explain how DNA sequence divergence and the activity of key mismatch repair proteins, MutS and MutL, lead to 1) genetic separation of closely related species (sympatric speciation) (Fitch and Ayala, Proceedings of the National Academy of Sciences, 1994, 91, 6717–6720), 2) the stability of genomes riddled by diverged repeated sequences, and 3) conservation of highly polymorphic DNA sequence blocks that constitute the immunological self. All three phenomena involve suppression of recombination between diverged homologies, resulting in prevention of gene sharing between closely related genomes (evolution of new species) as well as sequence sharing between closely related genes within a genome (e.g., evolution of immunoglobulin, MHC, and other gene families bearing conserved polymorphisms).

Mutators Enhance Adaptive Micro-Evolution in Pathogenic Microbes

Adaptation to the changing environmental conditions experienced within a host requires genetic diversity within a microbial population. Genetic diversity arises from mutations which occur due to DNA damage from exposure to exogenous environmental stresses or generated endogenously through respiration or DNA replication errors. As mutations can be deleterious, a delicate balance must be obtained between generating enough mutations for micro-evolution to occur while maintaining fitness and genomic integrity. Pathogenic microorganisms can actively modify their mutation rate to enhance adaptive micro-evolution by increasing expression of error-prone DNA polymerases or by mutating or decreasing expression of genes required for DNA repair. Strains which exhibit an elevated mutation rate are termed mutators. Mutators are found in varying prevalence in clinical populations where large-effect beneficial mutations enhance survival and are predominately caused by defects in the DNA mismatch repair (MMR) pathway. Mutators can facilitate the emergence of antibiotic resistance, allow phenotypic modifications to prevent recognition and destruction by the host immune system and enable switching to metabolic and cellular morphologies better able to survive in the given environment. This review will focus on recent advances in understanding the phenotypic and genotypic changes occurring in MMR mutators in both prokaryotic and eukaryotic pathogens.

A Combination of Metagenomic and Cultivation Approaches Reveals Hypermutator Phenotypes within Vibrio cholerae-Infected Patients

Vibrio cholerae can cause a range of symptoms, from severe diarrhea to asymptomatic infection. Previous studies using whole-genome sequencing (WGS) of multiple bacterial isolates per patient showed that V. cholerae can evolve modest genetic diversity during symptomatic infection. To further explore the extent of V. cholerae within-host diversity, we applied culture-based WGS and metagenomics to a cohort of both symptomatic and asymptomatic cholera patients from Bangladesh. While metagenomics allowed us to detect more mutations in symptomatic patients, WGS of cultured isolates was necessary to detect V. cholerae diversity in asymptomatic carriers, likely due to their low V. cholerae load. Using both metagenomics and isolate WGS, we report three lines of evidence that V. cholerae hypermutators evolve within patients. First, we identified nonsynonymous mutations in V. cholerae DNA repair genes in 5 out of 11 patient metagenomes sequenced with sufficient coverage of the V. cholerae genome and in 1 of 3 patients with isolate genomes sequenced. Second, these mutations in DNA repair genes tended to be accompanied by an excess of intrahost single nucleotide variants (iSNVs). Third, these iSNVs were enriched in transversion mutations, a known hallmark of hypermutator phenotypes. While hypermutators appeared to generate mostly selectively neutral mutations, nonmutators showed signs of convergent mutation across multiple patients, suggesting V. cholerae adaptation within hosts. Our results highlight the power and limitations of metagenomics combined with isolate sequencing to characterize within-patient diversity in acute V. cholerae infections, while providing evidence for hypermutator phenotypes within cholera patients.

IMPORTANCE Pathogen evolution within patients can impact phenotypes such as drug resistance and virulence, potentially affecting clinical outcomes. V. cholerae infection can result in life-threatening diarrheal disease or asymptomatic infection. Here, we describe whole-genome sequencing of V. cholerae isolates and culture-free metagenomic sequencing from stool of symptomatic cholera patients and asymptomatic carriers. Despite the typically short duration of cholera, we found evidence for adaptive mutations in the V. cholerae genome that occur independently and repeatedly within multiple symptomatic patients. We also identified V. cholerae hypermutator phenotypes within several patients, which appear to generate mainly neutral or deleterious mutations. Our work sets the stage for future studies of the role of hypermutators and within-patient evolution in explaining the variation from asymptomatic carriage to symptomatic cholera.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

Horizontal gene transfer and recombination analysis of SARS-CoV-2 genes helps discover its close relatives and shed light on its origin

BACKGROUND

The SARS-CoV-2 pandemic is one of  the greatest  global medical and social challenges that have emerged in recent history. Human coronavirus strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein.

RESULTS

We provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-CoV-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215-1425] of gene S and region [534-727] of gene N. Moreover, some statistically significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster.

CONCLUSIONS

The results of our horizontal gene transfer and recombination analysis suggest that SARS-CoV-2 could not only be a chimera virus resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of this dangerous virus.

An interbacterial DNA deaminase toxin directly mutagenizes surviving target populations

When bacterial cells come in contact, antagonism mediated by the delivery of toxins frequently ensues. The potential for such encounters to have long-term beneficial consequences in recipient cells has not been investigated. Here, we examined the effects of intoxication by DddA, a cytosine deaminase delivered via the type VI secretion system (T6SS) of Burkholderia cenocepacia. Despite its killing potential, we observed that several bacterial species resist DddA and instead accumulate mutations. These mutations can lead to the acquisition of antibiotic resistance, indicating that even in the absence of killing, interbacterial antagonism can have profound consequences on target populations. Investigation of additional toxins from the deaminase superfamily revealed that mutagenic activity is a common feature of these proteins, including a representative we show targets single-stranded DNA and displays a markedly divergent structure. Our findings suggest that a surprising consequence of antagonistic interactions between bacteria could be the promotion of adaptation via the action of directly mutagenic toxins.

Evolutionary epidemiology of Streptococcus iniae: Linking mutation rate dynamics with adaptation to novel immunological landscapes

Pathogens continuously adapt to changing host environments where variation in their virulence and antigenicity is critical to their long-term evolutionary success. The emergence of novel variants is accelerated in microbial mutator strains (mutators) deficient in DNA repair genes, most often from mismatch repair and oxidized-guanine repair systems (MMR and OG respectively). Bacterial MMR/OG mutants are abundant in clinical samples and show increased adaptive potential in experimental infection models, yet the role of mutators in the epidemiology and evolution of infectious disease is not well understood. Here we investigated the role of mutation rate dynamics in the evolution of a broad host range pathogen, Streptococcus iniae, using a set of 80 strains isolated globally over 40 years. We have resolved phylogenetic relationships using non-recombinant core genome variants, measured in vivo mutation rates by fluctuation analysis, identified variation in major MMR/OG genes and their regulatory regions, and phenotyped the major traits determining virulence in streptococci. We found that both mutation rate and MMR/OG genotype are remarkably conserved within phylogenetic clades but significantly differ between major phylogenetic lineages. Further, variation in MMR/OG loci correlates with occurrence of atypical virulence-associated phenotypes, infection in atypical hosts (mammals), and atypical (osseous) tissue of a vaccinated primary host. These findings suggest that mutators are likely to facilitate adaptations preceding major diversification events and may promote emergence of variation permitting colonization of a novel host tissue, novel host taxa (host jumps), and immune-escape in the vaccinated host.

Escherichia coli with a Tunable Point Mutation Rate for Evolution Experiments

The mutation rate and mutations' effects on fitness are crucial to evolution. Mutation rates are under selection due to linkage between mutation rate modifiers and mutations' effects on fitness. The linkage between a higher mutation rate and more beneficial mutations selects for higher mutation rates, while the linkage between a higher mutation rate and more deleterious mutations selects for lower mutation rates. The net direction of selection on mutations rates depends on the fitness landscape, and a great deal of work has elucidated the fitness landscapes of mutations. However, tests of the effect of varying a mutation rate on evolution in a single organism in a single environment have been difficult. This has been studied using strains of antimutators and mutators, but these strains may differ in additional ways and typically do not allow for continuous variation of the mutation rate. To help investigate the effects of the mutation rate on evolution, we have genetically engineered a strain of Escherichia coli with a point mutation rate that can be smoothly varied over two orders of magnitude. We did this by engineering a strain with inducible control of the mismatch repair proteins MutH and MutL. We used this strain in an approximately 350 generation evolution experiment with controlled variation of the mutation rate. We confirmed the construct and the mutation rate were stable over this time. Sequencing evolved strains revealed a higher number of single nucleotide polymorphisms at higher mutations rates, likely due to either the beneficial effects of these mutations or their linkage to beneficial mutations.

Ancestral niche separation and evolutionary rate differentiation between sister marine flavobacteria lineages

Marine flavobacteria are specialists for polysaccharide degradation. They dominate in habitats enriched with polysaccharides, but are also prevalent in pelagic environments where polysaccharides are less available. These niches are likely occupied by distinct lineages, but evolutionary processes underlying their niche differentiation remain elusive. Here, genomic analyses and physiological assays indicate that the sister flavobacteria lineages Leeuwenhoekiella and Nonlabens likely explore polysaccharide-rich macroalgae and polysaccharide-poor pelagic niches respectively. Phylogenomic analyses inferred that the niche separation likely occurred anciently and coincided with increased sequence evolutionary rate in Nonlabens compared with Leeuwenhoekiella. Further analyses ruled out the known mechanisms likely driving evolutionary rate acceleration, including reduced selection efficiency, decreased generation time and increased mutation rate. In particular, the mutation rates were determined using an unbiased experimental method, which measures the present-day populations and may not reflect ancestral populations. These data collectively lead to a new hypothesis that an ancestral and transient mutation rate increase resulted in evolutionary rate increase in Nonlabens. This hypothesis was supported by inferring that gains and losses of genes involved in SOS response, a mechanism known to drive transiently increased mutation rate, coincided with evolutionary rate acceleration. Our analyses highlight the evolutionary mechanisms underlying niche differentiation of flavobacteria lineages.

Salmonella Pathogenicity Island 1 (SPI-1): The Evolution and Stabilization of a Core Genomic Type Three Secretion System

Salmonella Pathogenicity Island 1 (SPI-1) encodes a type three secretion system (T3SS), effector proteins, and associated transcription factors that together enable invasion of epithelial cells in animal intestines. The horizontal acquisition of SPI-1 by the common ancestor of all Salmonella is considered a prime example of how gene islands potentiate the emergence of new pathogens with expanded niche ranges. However, the evolutionary history of SPI-1 has attracted little attention. Here, we apply phylogenetic comparisons across the family Enterobacteriaceae to examine the history of SPI-1, improving the resolution of its boundaries and unique architecture by identifying its composite gene modules. SPI-1 is located between the core genes fhlA and mutS, a hotspot for the gain and loss of horizontally acquired genes. Despite the plasticity of this locus, SPI-1 demonstrates stable residency of many tens of millions of years in a host genome, unlike short-lived homologous T3SS and effector islands including Escherichia ETT2, Yersinia YSA, Pantoea PSI-2, Sodalis SSR2, and Chromobacterium CPI-1. SPI-1 employs a unique series of regulatory switches, starting with the dedicated transcription factors HilC and HilD, and flowing through the central SPI-1 regulator HilA. HilA is shared with other T3SS, but HilC and HilD may have their evolutionary origins in Salmonella. The hilA, hilC, and hilD gene promoters are the most AT-rich DNA in SPI-1, placing them under tight control by the transcriptional repressor H-NS. In all Salmonella lineages, these three promoters resist amelioration towards the genomic average, ensuring strong repression by H-NS. Hence, early development of a robust and well-integrated regulatory network may explain the evolutionary stability of SPI-1 compared to T3SS gene islands in other species.

Mutation Rate Heterogeneity Increases Odds of Survival in Unpredictable Environments

Mutation rates affect both a population's present fitness and its capacity to adapt to future environmental changes. When the available genetic variability limits adaptation to environmental change, natural selection favors high mutations rates. However, constitutively high mutation rates compromise the fitness of a population in stable environments. This problem may be resolved if an increase in mutation rates is limited to times of stress, restricted to some genomic regions, and occurs only in a subpopulation of cells. Such within-population heterogeneity of mutation rates can result from genetic, environmental, and stochastic effects. The presence of subpopulations of transient mutator cells does not jeopardize the overall fitness of a population under stable environmental conditions. However, they can increase the odds of survival in changing environments because they represent reservoirs of increased genetic variability. This article presents evidence that such heterogeneity of mutation rates is more the norm than the exception.

Baker's Yeast Clinical Isolates Provide a Model for How Pathogenic Yeasts Adapt to Stress

Global outbreaks of drug-resistant fungi such as Candida auris are thought to be due at least in part to excessive use of antifungal drugs. Baker's yeast Saccharomyces cerevisiae has gained importance as an emerging opportunistic fungal pathogen that can cause infections in immunocompromised patients. Analyses of over 1000 S. cerevisiae isolates are providing rich resources to better understand how fungi can grow in human environments. A large percentage of clinical S. cerevisiae isolates are heterozygous across many nucleotide sites, and a significant proportion are of mixed ancestry and/or are aneuploid or polyploid. Such features potentially facilitate adaptation to new environments. These observations provide strong impetus for expanding genomic and molecular studies on clinical and wild isolates to understand the prevalence of genetic diversity and instability-generating mechanisms, and how they are selected for and maintained. Such work can also lead to the identification of new targets for antifungal drugs.

The Mutator Phenotype: Adapting Microbial Evolution to Cancer Biology

The mutator phenotype hypothesis was postulated almost 40 years ago to reconcile the observation that while cancer cells display widespread mutational burden, acquisition of mutations in non-transformed cells is a rare event. Moreover, it also suggested that cancer evolution could be fostered by increased genome instability. Given the evolutionary conservation throughout the tree of life and the genetic tractability of model organisms, yeast and bacterial species pioneered studies to dissect the functions of genes required for genome maintenance (caretaker genes) or for cell growth control (gatekeeper genes). In this review, we first provide an overview of what we learned from model organisms about the roles of these genes and the genome instability that arises as a consequence of their dysregulation. We then discuss our current understanding of how mutator phenotypes shape the evolution of bacteria and yeast species. We end by bringing clinical evidence that lessons learned from single-cell organisms can be applied to tumor evolution.

Mutation bias and GC content shape antimutator invasions

Mutators represent a successful strategy in rapidly adapting asexual populations, but theory predicts their eventual extinction due to their unsustainably large deleterious load. While antimutator invasions have been documented experimentally, important discrepancies among studies remain currently unexplained. Here we show that a largely neglected factor, the mutational idiosyncrasy displayed by different mutators, can play a major role in this process. Analysing phylogenetically diverse bacteria, we find marked and systematic differences in the protein-disruptive effects of mutations caused by different mutators in species with different GC compositions. Computer simulations show that these differences can account for order-of-magnitude changes in antimutator fitness for a realistic range of parameters. Overall, our results suggest that antimutator dynamics may be highly dependent on the specific genetic, ecological and evolutionary history of a given population. This context-dependency further complicates our understanding of mutators in clinical settings, as well as their role in shaping bacterial genome size and composition.

Evolution of Stress-Induced Mutagenesis in the Presence of Horizontal Gene Transfer

Stress-induced mutagenesis has been observed in multiple species of bacteria and yeast. It has been suggested that in asexual populations, a mutator allele that increases the mutation rate during stress can sweep to fixation with the beneficial mutations it generates. However, even asexual microbes can undergo horizontal gene transfer and rare recombination, which typically interfere with the spread of mutator alleles. Here we examine the effect of horizontal gene transfer on the evolutionary advantage of stress-induced mutator alleles. Our results demonstrate that stress-induced mutator alleles are favored by selection even in the presence of horizontal gene transfer and more so when the mutator alleles also increase the rate of horizontal gene transfer. We suggest that when regulated by stress, mutation and horizontal gene transfer can be complementary rather than competing adaptive strategies and that stress-induced mutagenesis has important implications for evolutionary biology, ecology, and epidemiology, even in the presence of horizontal gene transfer and rare recombination.

Pourquoi le taux de mutation n’est-il jamais égal à zéro ?

Alfred H. Sturtevant fut le premier à s’en étonner : le taux de mutation est faible mais n’atteint jamais zéro. Pourtant, la plupart des mutations qui modifient le phénotype ont un effet délétère, les individus qui produisent le moins de mutants génèrent donc plus de descendants viables et fertiles. La sélection naturelle devrait ainsi progressivement faire tendre le taux de mutation vers zéro au cours des générations. Des analyses récentes suggèrent que ce taux dépend principalement de la taille efficace des génomes et de l’effectif efficace des populations. Le maintien de taux de mutation plus élevés que nécessaire illustrerait les limites de la sélection naturelle dans un monde vivant constitué de populations de taille finie.

Genomic Characteristics of Desulfonema ishimotonii Tokyo 01T Implying Horizontal Gene Transfer Among Phylogenetically Dispersed Filamentous Gliding Bacteria

Desulfonema ishimotonii strain Tokyo 01^T is a filamentous sulfate-reducing bacterium isolated from a marine sediment. In this study, the genome of this strain was sequenced and analyzed with a focus on gene transfer from phylogenetically distant organisms. While the strain belongs to the class Deltaproteobacteria, hundreds of proteins encoded in the genome showed the highest sequence similarities to those of organisms outside of the class Deltaproteobacteria, suggesting that more than 20% of the genome is putatively of foreign origins. Many of these proteins had the highest sequence identities with proteins encoded in the genomes of filamentous bacteria, including giant sulfur oxidizers of the orders Thiotrichales, cyanobacteria of various genera, and uncultured bacteria of the candidate phylum KSB3. As mobile genetic elements transferred from phylogenetically distant organisms, putative inteins were identified in the GyrB and DnaE proteins encoded in the genome of strain Tokyo 01^T. Genes involved in DNA recombination and repair were enriched in comparison to the closest relatives in the same family. Some of these genes were also related to those of organisms outside of the class Deltaproteobacteria, suggesting that they were acquired by horizontal gene transfer from diverse bacteria. The genomic data suggested significant genetic transfer among filamentous gliding bacteria in phylogenetically dispersed lineages including filamentous sulfate reducers. This study provides insights into the genomic evolution of filamentous bacteria belonging to diverse lineages, characterized by various physiological functions and different ecological roles.

Incompatibilities in Mismatch Repair Genes MLH1-PMS1 Contribute to a Wide Range of Mutation Rates in Human Isolates of Baker's Yeast

Laboratory baker's yeast strains bearing an incompatible combination of MLH1 and PMS1 mismatch repair alleles are mutators that can adapt more rapidly to stress, but do so at the cost of long-term fitness. We identified 18 baker's yeast isolates from 1011 surveyed that contain the incompatible MLH1-PMS1 genotype in a heterozygous state. Surprisingly, the incompatible combination from two human clinical heterozygous diploid isolates, YJS5845 and YJS5885, contain the exact MLH1 (S288c-derived) and PMS1 (SK1-derived) open reading frames originally shown to confer incompatibility. While these isolates were nonmutators, their meiotic spore clone progeny displayed mutation rates in a DNA slippage assay that varied over a 340-fold range. This range was 30-fold higher than observed between compatible and incompatible combinations of laboratory strains. Genotyping analysis indicated that MLH1-PMS1 incompatibility was the major driver of mutation rate in the isolates. The variation in the mutation rate of incompatible spore clones could be due to background suppressors and enhancers, as well as aneuploidy seen in the spore clones. Our data are consistent with the observed variance in mutation rate contributing to adaptation to stress conditions (e.g., in a human host) through the acquisition of beneficial mutations, with high mutation rates leading to long-term fitness costs that are buffered by mating or eliminated through natural selection.

New Biological Insights Into How Deforestation in Amazonia Affects Soil Microbial Communities Using Metagenomics and Metagenome-Assembled Genomes

Deforestation in the Brazilian Amazon occurs at an alarming rate, which has broad effects on global greenhouse gas emissions, carbon storage, and biogeochemical cycles. In this study, soil metagenomes and metagenome-assembled genomes (MAGs) were analyzed for alterations to microbial community composition, functional groups, and putative physiology as it related to land-use change and tropical soil. A total of 28 MAGs were assembled encompassing 10 phyla, including both dominant and rare biosphere lineages. Amazon Acidobacteria subdivision 3, Melainabacteria, Microgenomates, and Parcubacteria were found exclusively in pasture soil samples, while Candidatus Rokubacteria was predominant in the adjacent rainforest soil. These shifts in relative abundance between land-use types were supported by the different putative physiologies and life strategies employed by the taxa. This research provides unique biological insights into candidate phyla in tropical soil and how deforestation may impact the carbon cycle and affect climate change.

Horizontal gene transfer of Chlamydia: Novel insights from tree reconciliation

Recent comparative genomics studies have suggested that horizontal gene transfer (HGT) is one of the major processes in bacterial evolution. In this study, HGT events of 64 Chlamydia strains were investigated based on the pipeline employed in HGTree database constructed in our recent study. Tree reconciliation method was applied in order to calculate feasible HGT events. Following initial detection and an evaluation procedure, evidence of the HGT was identified in 548 gene families including 42 gene families transferred from outside of Chlamydiae phylum with high reliability. The donor species of inter-phylum HGT consists of 12 different bacterial and archaeal phyla, suggesting that Chlamydia might have even more various host range than in previous reports. In addition, each species of Chlamydia showed varying preference towards HGT, and genes engaged in HGT within Chlamydia and between other species showed different functional distribution. Also, examination of individual gene flows of niche-specific genes suggested that many of such genes are transferred mainly within Chlamydia genus. Our results uncovered novel features of HGT acting on Chlamydia genome evolution, and it would be also strong evidence that HGT is an ongoing process for intracellular pathogens. We expect that the results provide more insight into lineage- and niche-specific adaptations regarding their infectivity and pathogenicity.

Mutator genomes decay, despite sustained fitness gains, in a long-term experiment with bacteria

Understanding the extreme variation among bacterial genomes remains an unsolved challenge in evolutionary biology, despite long-standing debate about the relative importance of natural selection, mutation, and random drift. A potentially important confounding factor is the variation in mutation rates between lineages and over evolutionary history, which has been documented in several species. Mutation accumulation experiments have shown that hypermutability can erode genomes over short timescales. These results, however, were obtained under conditions of extremely weak selection, casting doubt on their general relevance. Here, we circumvent this limitation by analyzing genomes from mutator populations that arose during a long-term experiment with Escherichia coli, in which populations have been adaptively evolving for >50,000 generations. We develop an analytical framework to quantify the relative contributions of mutation and selection in shaping genomic characteristics, and we validate it using genomes evolved under regimes of high mutation rates with weak selection (mutation accumulation experiments) and low mutation rates with strong selection (natural isolates). Our results show that, despite sustained adaptive evolution in the long-term experiment, the signature of selection is much weaker than that of mutational biases in mutator genomes. This finding suggests that relatively brief periods of hypermutability can play an outsized role in shaping extant bacterial genomes. Overall, these results highlight the importance of genomic draft, in which strong linkage limits the ability of selection to purge deleterious mutations. These insights are also relevant to other biological systems evolving under strong linkage and high mutation rates, including viruses and cancer cells.

Characterization of the DNA mismatch repair proteins MutS and MutL in a hypermutator Acinetobacter baumannii

Mutations of mutS and mutL genes have been linked with the emergence of hypermutator (HPM) phenotype in several bacteria. Nevertheless, there is scarce evidence that these mutations occurred in HPM Acinetobacter baumannii, therefore, it remains unknown whether the mutations located in domains mediating the functions of MutS and MutL. To address this information gap, the nucleotide sequences of mutS and mutL were characterized and their mutations were identified. Additionally, we proposed in silico models of mutated proteins and analyzed the secondary and tertiary structures, and the interaction interfaces of MutL and MutS. The HPM A. baumannii and a wild-type strain were subjected to PCR amplification of full length mutS and mutL, cloning, and sequencing. Following several reads of both strands of each gene and sequence assembly, the mutations were identified. Thereafter, the three-dimensional (3-D) structure of A. baumannii ATCC 19606 was developed and utilized as a template for homology modeling of the mutated amino acid sequences using the Phyre² and I-TASSER, VMD 1.9.3, LigPlus v.1.4.5, PyMOL v.0.99 software. Regardless of silent mutations (n = 43), 11 missense mutations were identified in the MutS domains of HPM strain such as A4T, T272S, D278N in N-terminus, connector, and core domains, respectively. Three mutations -I357T, A408S, N447S- and 16 silent mutations were observed in MutL. Secondary structure prediction of MutS revealed that the amount of alpha helices, beta sheets, and coils in HPM were 35, 29, and 63, respectively, while these values were 36, 28, and 63 for A. baumannii ATCC 19606 as non mutator. In the case of MutL, for both HPM and non-mutator, 20, 21, and 39 of complete protein were alpha helices, beta sheets, and coils, respectively. Superimposition of structures of MutS of HPM on non-mutator revealed that T272, D278, G457, S528, A533, Y715, and E747 are closely matched with S272, D278, A457, P528, V533, C715, and K747, respectively in non-mutator strain. When the structure of MutL model in HPM was superimposed on its counterpart in non-mutator, all but residues S447, S408, and T357 were identical. Many mutations along the mutS and mutL were noted, but most of the mutations were observed in the interaction interfaces of MutS and MutL. Other substitutions were predominantly detected in C-terminus of MutS that may lead to reduced ATP binding and hydrolysis. Three substitution mutations were adjacent to C-terminus of MutL and are raising the suggestion of reduction in MutL dimerization. It seems that a combination of these mutations is implicated in increased mutation frequency and accordingly emergence of HPM strain.

Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii

Pathogenic microbes confront an evolutionary conflict between the pressure to maintain genome stability and the need to adapt to mounting external stresses. Bacteria often respond with elevated mutation rates, but little evidence exists of stable eukaryotic hypermutators in nature. Whole genome resequencing of the human fungal pathogen Cryptococcus deuterogattii identified an outbreak lineage characterized by a nonsense mutation in the mismatch repair component MSH2. This defect results in a moderate mutation rate increase in typical genes, and a larger increase in genes containing homopolymer runs. This allows facile inactivation of genes with coding homopolymer runs including FRR1, which encodes the target of the immunosuppresive antifungal drugs FK506 and rapamycin. Our study identifies a eukaryotic hypermutator lineage spread over two continents and suggests that pathogenic eukaryotic microbes may experience similar selection pressures on mutation rate as bacterial pathogens, particularly during long periods of clonal growth or while expanding into new environments.

As humans, we often think of genetic mutations as being bad. Over the past several decades we have seen health warnings issued on a variety of environmental exposures, from cigarettes to tanning beds, and with good reason because they cause mutations. For multicellular organisms like humans, these mutations are strongly associated with cancer. But in bacteria, this is not true. In fact, the rate at which mutations occur sometimes increases to help bacteria cope with stressful environments.

Unlike bacteria, humans are eukaryotes – the name given to organisms whose cells contain different compartments separated by membranes, such as the nucleus of the cell. For years, we have assumed that eukaryotic microbes, like fungi and parasites, act more like humans than like bacteria because work in budding yeast (another eukaryote) has suggested this to be the case. However, recent work in disease-causing fungi has shown that, much like bacteria, elevated mutation rates may help them to respond to stress. This could also enable fungi to become resistant to drugs used to treat fungal infections.

Cryptococcus deuterogattii is a fungus that causes human diseases including meningoencephalitis and a lung infection called pulmonary cryptococcosis. An ongoing outbreak of the fungus began in the Pacific Northwest of Canada in the late 1990s and emerged in the United States in 2006/2007. Among isolates closely related to those fungi causing the outbreak, three were found that appear to have a specific mutation in their DNA mismatch repair pathway, meaning that they may also experience a higher mutation rate. These strains are also less able to cause disease than others.

Billmyre et al. now demonstrate experimentally that all three isolates have a specific DNA mismatch repair defect, and show that these fungi experience elevated mutation rates, resulting in what is known as a hypermutator state. Furthermore, whole genome sequencing and phylogenetic analysis showed that these hypermutator strains are derived from the outbreak-causing fungi, and that their reduced ability to cause disease is likely a result of accumulating mutations and the loss of the ability to grow at the higher temperatures found in the human body.

Fungal infections are difficult to treat, in part because there are a limited number of available drugs. Elevated mutation rates will likely increase how often and how rapidly fungi develop resistance to these drugs. Understanding how commonly fungi exhibit a hypermutator state that could impact the development of drug resistance will therefore be important for treating patients with fungal infections, which account for millions of infections and hundreds of thousands of deaths annually worldwide.

Constraint and Contingency Pervade the Emergence of Novel Phenotypes in Complex Metabolic Systems

An evolutionary constraint is a bias or limitation in phenotypic variation that a biological system produces. We know examples of such constraints, but we have no systematic understanding about their extent and causes for any one biological system. We here study metabolisms, genomically encoded complex networks of enzyme-catalyzed biochemical reactions, and the constraints they experience in bringing forth novel phenotypes that allow survival on novel carbon sources. Our computational approach does not limit us to analyzing constrained variation in any one organism, but allows us to quantify constraints experienced by any metabolism. Specifically, we study metabolisms that are viable on one of 50 different carbon sources, and quantify how readily alterations of their chemical reactions create the ability to survive on a novel carbon source. We find that some metabolic phenotypes are much less likely to originate than others. For example, metabolisms viable on D-glucose are 1835 times more likely to give rise to metabolisms viable on D-fructose than on acetate. Likewise, we observe that some novel metabolic phenotypes are more contingent on parental phenotypes than others. Biochemical similarities among carbon sources can help explain the causes of these constraints. In addition, we study metabolisms that can be produced by recombination among 55 metabolisms of different bacterial strains or species, and show that their novel phenotypes are also contingent on and constrained by parental genotypes. To our knowledge, our analysis is the first to systematically quantify the incidence of constrained evolution in a broad class of biological system that is central to life and its evolution.

Contribution of increased mutagenesis to the evolution of pollutants-degrading indigenous bacteria

Bacteria can rapidly evolve mechanisms allowing them to use toxic environmental pollutants as a carbon source. In the current study we examined whether the survival and evolution of indigenous bacteria with the capacity to degrade organic pollutants could be connected with increased mutation frequency. The presence of constitutive and transient mutators was monitored among 53 pollutants-degrading indigenous bacterial strains. Only two strains expressed a moderate mutator phenotype and six were hypomutators, which implies that constitutively increased mutability has not been prevalent in the evolution of pollutants degrading bacteria. At the same time, a large proportion of the studied indigenous strains exhibited UV-irradiation-induced mutagenesis, indicating that these strains possess error-prone DNA polymerases which could elevate mutation frequency transiently under the conditions of DNA damage. A closer inspection of two Pseudomonas fluorescens strains PC20 and PC24 revealed that they harbour genes for ImuC (DnaE2) and more than one copy of genes for Pol V. Our results also revealed that availability of other nutrients in addition to aromatic pollutants in the growth environment of bacteria affects mutagenic effects of aromatic compounds. These results also implied that mutagenicity might be affected by a factor of how long bacteria have evolved to use a particular pollutant as a carbon source.

A null model for microbial diversification

Whether prokaryotes (Bacteria and Archaea) are naturally organized into phenotypically and genetically cohesive units comparable to animal or plant species remains contested, frustrating attempts to estimate how many such units there might be, or to identify the ecological roles they play. Analyses of gene sequences in various closely related prokaryotic groups reveal that sequence diversity is typically organized into distinct clusters, and processes such as periodic selection and extensive recombination are understood to be drivers of cluster formation ("speciation"). However, observed patterns are rarely compared with those obtainable with simple null models of diversification under stochastic lineage birth and death and random genetic drift. Via a combination of simulations and analyses of core and phylogenetic marker genes, we show that patterns of diversity for the genera Escherichia, Neisseria, and Borrelia are generally indistinguishable from patterns arising under a null model. We suggest that caution should thus be taken in interpreting observed clustering as a result of selective evolutionary forces. Unknown forces do, however, appear to play a role in Helicobacter pylori, and some individual genes in all groups fail to conform to the null model. Taken together, we recommend the presented birth-death model as a null hypothesis in prokaryotic speciation studies. It is only when the real data are statistically different from the expectations under the null model that some speciation process should be invoked.

Perspectiva histórica del origen evolutivo de la resistencia a antibióticos

La resistencia a antimicrobianos representa un aspecto natural de evolución bacteriana, que puede resultar de mutaciones o por adquisición de genes foráneos. Hay diferentes posturas sobre el origen de ésta resistencia que explican la habilidad de estos microorganismos de adquirir nuevas características. Las teorías de la evolución de Lamarck y Darwin, han dado pie a experimentos diseñados para explorar el origen de la variación bacteriana y surgimiento de nuevas características. Estos estudios muestran que la resistencia está relacionada con mutaciones en genes cromosomales y/o la transferencia de elementos genéticos extracromosomales, que se expresan según la presión antibiótica ejercida. Está revisión recopila los principales experimentos y las conclusiones derivadas para explicar el fenómeno de resistencia a antibióticos.

Mismatch Repair Incompatibilities in Diverse Yeast Populations

An elevated mutation rate can provide cells with a source of mutations to adapt to changing environments. We identified a negative epistatic interaction involving naturally occurring variants in the MLH1 and PMS1 mismatch repair (MMR) genes of Saccharomyces cerevisiae We hypothesized that this MMR incompatibility, created through mating between divergent S. cerevisiae, yields mutator progeny that can rapidly but transiently adapt to an environmental stress. Here we analyzed the MLH1 and PMS1 genes across 1010 S. cerevisiae natural isolates spanning a wide range of ecological sources (tree exudates, Drosophila, fruits, and various fermentation and clinical isolates) and geographical sources (Europe, America, Africa, and Asia). We identified one homozygous clinical isolate and 18 heterozygous isolates containing the incompatible MMR genotype. The MLH1-PMS1 gene combination isolated from the homozygous clinical isolate conferred a mutator phenotype when expressed in the S288c laboratory background. Using a novel reporter to measure mutation rates, we showed that the overall mutation rate in the homozygous incompatible background was similar to that seen in compatible strains, indicating the presence of suppressor mutations in the clinical isolate that lowered its mutation rate. This observation and the identification of 18 heterozygous isolates, which can lead to MMR incompatible genotypes in the offspring, are consistent with an elevated mutation rate rapidly but transiently facilitating adaptation. To avoid long-term fitness costs, the incompatibility is apparently buffered by mating or by acquiring suppressors. These observations highlight effective strategies in eukaryotes to avoid long-term fitness costs associated with elevated mutation rates.

A Mobile Element in mutS Drives Hypermutation in a Marine Vibrio

Bacteria face a trade-off between genetic fidelity, which reduces deleterious mistakes in the genome, and genetic innovation, which allows organisms to adapt. Evidence suggests that many bacteria balance this trade-off by modulating their mutation rates, but few mechanisms have been described for such modulation. Following experimental evolution and whole-genome resequencing of the marine bacterium Vibrio splendidus 12B01, we discovered one such mechanism, which allows this bacterium to switch to an elevated mutation rate. This switch is driven by the excision of a mobile element residing in mutS, which encodes a DNA mismatch repair protein. When integrated within the bacterial genome, the mobile element provides independent promoter and translation start sequences for mutS—different from the bacterium’s original mutS promoter region—which allow the bacterium to make a functional mutS gene product. Excision of this mobile element rejoins the mutS gene with host promoter and translation start sequences but leaves a 2-bp deletion in the mutS sequence, resulting in a frameshift and a hypermutator phenotype. We further identified hundreds of clinical and environmental bacteria across Betaproteobacteria and Gammaproteobacteria that possess putative mobile elements within the same amino acid motif in mutS. In a subset of these bacteria, we detected excision of the element but not a frameshift mutation; the mobile elements leave an intact mutS coding sequence after excision. Our findings reveal a novel mechanism by which one bacterium alters its mutation rate and hint at a possible evolutionary role for mobile elements within mutS in other bacteria.

DNA mutations are a double-edged sword. Most mutations are harmful; they can scramble precise genetic sequences honed over thousands of generations. But in rare cases, mutations also produce beneficial new traits that allow populations to adapt to changing environments. Recent evidence suggests that some bacteria balance this trade-off by altering their mutation rates to suit their environment. To date, however, we know of few mechanisms that allow bacteria to change their mutation rates. We describe one such mechanism, driven by the action of a mobile element, in the marine bacterium Vibrio splendidus 12B01. We also found similar mobile genetic sequences in the mutS genes of many different bacteria, including clinical and agricultural pathogens. These mobile elements might play an as yet unknown role in the evolution of these important bacteria.

Mutability in Pseudomonas viridiflava as a programmed balance between antibiotic resistance and pathogenicity

Mutable bacterial cells are defective in their DNA repair system and often have a phenotype different from that of their wild-type counterparts. In human bacterial pathogens, the mutable and hypermutable phenotypes are often associated with general antibiotic resistance. Here, we quantified the occurrence of mutable cells in Pseudomonas viridiflava, a phytopathogenic bacterium in the P. syringae complex with a broad host range and capacity to live as a saprophyte. Two phenotypic variants (transparent and mucoid) were produced by this bacterium. The transparent variant had a mutator phenotype, showed general antibiotic resistance and could not induce disease on the plant species tested (bean). In contrast, the mucoid variant did not display mutability or resistance to antibiotics and was capable of inducing disease on bean. Both the transparent and mucoid variants were less fit when grown in vitro, whereas, in planta, both of the variants and wild-types attained similar population densities. Given the importance of the methyl-directed mismatch repair system (MMR) in the occurrence of mutable and hypermutable cells in human bacterial pathogens, we investigated whether mutations in mut genes were associated with mutator transparent cells in P. viridiflava. Our results showed no mutations in MMR genes in any of the P. viridiflava cells tested. Here, we report that a high mutation rate and antibiotic resistance are inversely correlated with pathogenicity in P. viridiflava, but are not associated with mutations in MMR. In addition, P. viridiflava variants differ from variants produced by other phytopathogenic bacteria in the absence of reversion to the wild-type phenotype.

DNA Repair Is Associated with Information Content in Bacteria, Archaea, and DNA Viruses

The concept of a "proteomic constraint" proposes that DNA repair capacity is positively correlated with the information content of a genome, which can be approximated to the size of the proteome (P). This in turn implies that DNA repair genes are more likely to be present in genomes with larger values of P. This stands in contrast to the common assumption that informational genes have a core function and so are evenly distributed across organisms. We examined the presence/absence of 18 DNA repair genes in bacterial genomes. A positive relationship between gene presence and P was observed for 17 genes in the total dataset, and 16 genes when only nonintracellular bacteria were examined. A marked reduction of DNA repair genes was observed in intracellular bacteria, consistent with their reduced value of P. We also examined archaeal and DNA virus genomes, and show that the presence of DNA repair genes is likewise related to a larger value of P. In addition, the products of the bacterial genes mutY, vsr, and ndk, involved in the correction of GC/AT mutations, are strongly associated with reduced genome GC content. We therefore propose that a reduction in information content leads to a loss of DNA repair genes and indirectly to a reduction in genome GC content in bacteria by exposure to the underlying AT mutation bias. The reduction in P may also indirectly lead to the increase in substitution rates observed in intracellular bacteria via loss of DNA repair genes.

A Genetic Incompatibility Accelerates Adaptation in Yeast

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker’s yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

In nature, bacterial populations with high mutation rates can adapt faster to new environments by acquiring beneficial mutations. However, such populations also accumulate harmful mutations that reduce their fitness. We show that the model eukaryote baker’s yeast can use a similar mutator strategy to adapt to new environments. The mutator state that we observed resulted from an incompatibility involving two genes, MLH1 and PMS1, that work together to remove DNA replication errors through a spellchecking mismatch repair mechanism. This incompatibility can occur through mating between baker’s yeast from different genetic backgrounds, yielding mutator offspring containing an MLH1-PMS1 combination not present in either parent. Interestingly, these offspring adapted more rapidly to stress, compared to the parental strains, and did so without an overall loss in fitness. DNA sequencing analyses of baker’s yeast strains from across the globe support the presence of incompatible hybrid yeast strains in nature. These observations provide a powerful model to understand how the segregation of defects in DNA mismatch repair can serve as an effective strategy to enable eukaryotes to adapt to changing environments.

Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.

Development of a stress-induced mutagenesis module for autonomous adaptive evolution of Escherichia coli to improve its stress tolerance

BACKGROUND

Microbial tolerance to different environmental stresses is of importance for efficient production of biofuels and biochemical. Such traits are often improved by evolutionary engineering approaches including mutagen-induced mutagenesis and successive passage. In contrast to these approaches which generate mutations in rapidly growing cells, recent research showed that mutations could be generated in non-dividing cells under stressful but non-lethal conditions, leading to the birth of the theory of stress-induced mutagenesis (SIM). A molecular mechanism of SIM has been elucidated to be mutagenic repair of DNA breaks. This inspired us to develop a synthetic SIM module to simulate the mutagenic cellular response so as to accelerate microbial adaptive evolution for an improved stress tolerance.

RESULTS

A controllable SIM evolution module was devised based on a genetic toggle switch in Escherichia coli. The synthetic module enables expression and repression of the genes related to up- and down-regulation responses during SIM in a bistable way. Upon addition of different inducers, the module can be turned on or off, triggering transition to a mutagenic or a high-fidelity state and thus allowing periodic adaptive evolution. Six genes (recA, dinB, umuD, ropS, ropE, and nusA) in the up-regulation responses were evaluated for their potentials to enhance the SIM rate. Expression of recA, dinB, or ropS alone increased the SIM rate by 4.5- to 13.7-fold, whereas their combined expression improved the rate by 31.9-fold. Besides, deletion of mutL increased the SIM rate by 82-fold. Assembly of these genes into the SIM module in the mutL-deletion E. coli strain elevated the SIM rate by nearly 3000-fold. Accelerated adaptive evolution of E. coli equipped with this synthetic SIM module was demonstrated under n-butanol stress, with the minimal inhibitory concentration of n-butanol increasing by 56 % within 2.5 months.

CONCLUSIONS

A synthetic SIM module was constructed to simulate cellular mutagenic responses during SIM. Based on this, a novel evolutionary engineering approach-SIM-based adaptive evolution-was developed to improve the n-butanol tolerance of E. coli.

From passengers to drivers: Impact of bacterial transposable elements on evolvability

Microbes have several mechanisms that promote evolutionary adaptation in stressful environments. The corresponding molecular pathways promote diversity through modulating rates of recombination, mutation or influence the activity of transposable genetic elements. Recent experimental studies suggest an evolutionary conflict between these mechanisms. Specifically, presence of mismatch repair mutator alleles in a bacterial population dramatically reduced fixation of bacterial insertion sequence elements. When rare, these elements had only a limited impact on adaptive evolution compared with other mutation-generating pathways. IS elements may initially spread like molecular parasites, but once present in many copies in a given genome, they might become generators of novelty during bacterial evolution.

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

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

Correlation between the genomic o454-nlpD region polymorphisms, virulence gene equipment and phylogenetic group of extraintestinal Escherichia coli (ExPEC) enables pathotyping irrespective of host, disease and source of isolation

BACKGROUND

The mutS-rpoS intergenic region in E. coli displays a mosaic structure which revealed pathotype specific patterns. To assess the importance of this region as a surrogate marker for the identification of highly virulent extraintestinal pathogenic E. coli (ExPEC) strains we aimed to: (i) characterize the genetic diversity of the mutS gene and the o454-nlpD genomic region among 510 E. coli strains from animals and humans; (ii) delineate associations between the polymorphism of this region and features such as phylogenetic background of E. coli, pathotype, host species, clinical condition, serogroup and virulence associated genes (VAG)s; and (iii) identify the most important VAGs for classification of the o454-nlpD region.

METHODS

Size variation in the o454-nlpD region was investigated by PCR amplification and sequencing. Phylogenetic relationships were assessed by Ecor- and Multilocus sequence- typing (MLST), and a comparative analysis between mutS gene phylogenetic tree obtained with RAxML and the MLST grouping method was performed. Correlation between o454-nlpD patterns and the features described above were analysed. In addition, the importance of 47 PCR-amplified ExPEC-related VAGs for classification of o454-nlpD patterns was investigated by means of Random Forest algorithm.

RESULTS

Four main structures (patterns I-IV) of the o454-nlpD region among ExPEC and commensal E. coli strains were identified. Statistical analysis showed a positive and exclusive association between pattern III and the ExPEC strains. A strong association between pattern III and either the Ecor group B2 or the sequence type complexes known to represent the phylogenetic background of highly virulent ExPEC strains (such as STC95, STC73 and STC131) was found as well. RF analyses determined five genes (csgA, malX, chuA, sit, and vat) to be suitable to predict pattern III strains.

CONCLUSION

The significant association between pattern III and group B2 strains suggested the o454-nlpD region to be of great value in identifying highly virulent strains among the mixed population of E. coli promising to be the basis of a future typing tool for ExPEC and their gut reservoir. Furthermore, top-ranked VAGs for classification and prediction of pattern III were identified. These data are most valuable for defining ExPEC pathotype in future in vivo assays.

Microbial communities evolve faster in extreme environments

Evolutionary analysis of microbes at the community level represents a new research avenue linking ecological patterns to evolutionary processes, but remains insufficiently studied. Here we report a relative evolutionary rates (rERs) analysis of microbial communities from six diverse natural environments based on 40 metagenomic samples. We show that the rERs of microbial communities are mainly shaped by environmental conditions, and the microbes inhabiting extreme habitats (acid mine drainage, saline lake and hot spring) evolve faster than those populating benign environments (surface ocean, fresh water and soil). These findings were supported by the observation of more relaxed purifying selection and potentially frequent horizontal gene transfers in communities from extreme habitats. The mechanism of high rERs was proposed as high mutation rates imposed by stressful conditions during the evolutionary processes. This study brings us one stage closer to an understanding of the evolutionary mechanisms underlying the adaptation of microbes to extreme environments.

Experimental evolution and the dynamics of genomic mutation rate modifiers

Because genes that affect mutation rates are themselves subject to mutation, mutation rates can be influenced by natural selection and other evolutionary forces. The population genetics of mutation rate modifier alleles has been a subject of theoretical interest for many decades. Here, we review experimental contributions to our understanding of mutation rate modifier dynamics. Numerous evolution experiments have shown that mutator alleles (modifiers that elevate the genomic mutation rate) can readily rise to high frequencies via genetic hitchhiking in non-recombining microbial populations. Whereas these results certainly provide an explanatory framework for observations of sporadically high mutation rates in pathogenic microbes and in cancer lineages, it is nonetheless true that most natural populations have very low mutation rates. This raises the interesting question of how mutator hitchhiking is suppressed or its phenotypic effect reversed in natural populations. Very little experimental work has addressed this question; with this in mind, we identify some promising areas for future experimental investigation.

Mechanisms of ectopic gene conversion

Gene conversion (conversion), the unidirectional transfer of DNA sequence information, occurs as a byproduct of recombinational repair of broken or damaged DNA molecules. Whereas excision repair processes replace damaged DNA by copying the complementary sequence from the undamaged strand of duplex DNA, recombinational mechanisms copy similar sequence, usually in another molecule, to replace the damaged sequence. In mitotic cells the other molecule is usually a sister chromatid, and the repair does not lead to genetic change. Less often a homologous chromosome or homologous sequence in an ectopic position is used. Conversion results from repair in two ways. First, if there was a double-strand gap at the site of a break, homologous sequence will be used as the template for synthesis to fill the gap, thus transferring sequence information in both strands. Second, recombinational repair uses complementary base pairing, and the heteroduplex molecule so formed is a source of conversion, both as heteroduplex and when donor (undamaged template) information is retained after correction of mismatched bases in heteroduplex. There are mechanisms that favour the use of sister molecules that must fail before ectopic homology can be used. Meiotic recombination events lead to the formation of crossovers required in meiosis for orderly segregation of pairs of homologous chromosomes. These events result from recombinational repair of programmed double-strand breaks, but in contrast with mitotic recombination, meiotic recombinational events occur predominantly between homologous chromosomes, so that transfer of sequence differences by conversion is very frequent. Transient recombination events that do not form crossovers form both between homologous chromosomes and between regions of ectopic homology, and leave their mark in the occurrence of frequent non-crossover conversion, including ectopic conversion.

Frequency-dependent selection can lead to evolution of high mutation rates

Theoretical and experimental studies have shown that high mutation rates can be advantageous, especially in novel or fluctuating environments. Here we examine how frequency-dependent competition may lead to fluctuations in trait frequencies that exert upward selective pressure on mutation rates. We use a mathematical model to show that cyclical trait dynamics generated by "rock-paper-scissors" competition can cause the mutation rate in a population to converge to a high evolutionarily stable mutation rate, reflecting a trade-off between generating novelty and reproducing past success. Introducing recombination lowers the evolutionarily stable mutation rate but allows stable coexistence between mutation rates above and below the evolutionarily stable rate. Even considering strong mutational load and ignoring the costs of faithful replication, evolution favors positive mutation rates if the selective advantage of prevailing in competition exceeds the ratio of recombining to nonrecombining offspring. We discuss a number of genomic mechanisms that may meet our theoretical requirements for the adaptive evolution of mutation. Overall, our results suggest that local mutation rates may be higher on genes influencing cyclical competition and that global mutation rates in asexual species may be higher in populations subject to strong cyclical competition.

Molecular applications for identifying microbial pathogens in the post-9/11 era

Rapid advances in molecular and optical technologies over the past 10 years have dramatically impacted the way biologic research is conducted today. Examples include microarrays, capillary sequencing, optical mapping and real-time sequencing (Pyrosequencing). These technologies are capable of rapidly delivering massive amounts of genetic information and are becoming routine mainstays of many laboratories. Fortunately, advances in scientific computing have provided the enormous computing power necessary to analyze these enormous data sets. The application of molecular technologies should prove useful to the burgeoning field of microbial forensics. In the post-9/11 era, when securing America's food supply is a major endeavor, the need for rapid identification of microbes that accidentally or intentionally find their way into foods is apparent. The principle that distinguishes a microbial forensic investigation from a molecular epidemiology study is that a biocrime has been committed. If proper attribution is to be attained, a link must be made between a particular microbe in the food and the perpetrator who placed it there. Therefore, the techniques used must be able to discriminate individual isolates of a particular microbe. A battery of techniques in development for distinguishing individual isolates of particular foodborne pathogens is discussed.