Transcription increases multiple spontaneous point mutations in Salmonella enterica

Mutation bias interacts with composition bias to influence adaptive evolution

Mutation is a biased stochastic process, with some types of mutations occurring more frequently than others. Previous work has used synthetic genotype-phenotype landscapes to study how such mutation bias affects adaptive evolution. Here, we consider 746 empirical genotype-phenotype landscapes, each of which describes the binding affinity of target DNA sequences to a transcription factor, to study the influence of mutation bias on adaptive evolution of increased binding affinity. By using empirical genotype-phenotype landscapes, we need to make only few assumptions about landscape topography and about the DNA sequences that each landscape contains. The latter is particularly important because the set of sequences that a landscape contains determines the types of mutations that can occur along a mutational path to an adaptive peak. That is, landscapes can exhibit a composition bias—a statistical enrichment of a particular type of mutation relative to a null expectation, throughout an entire landscape or along particular mutational paths—that is independent of any bias in the mutation process. Our results reveal the way in which composition bias interacts with biases in the mutation process under different population genetic conditions, and how such interaction impacts fundamental properties of adaptive evolution, such as its predictability, as well as the evolution of genetic diversity and mutational robustness.

Mutation is often depicted as a random process due its unpredictable nature. However, such randomness does not imply uniformly distributed outcomes, because some DNA sequence changes happen more frequently than others. Mutation bias can be an orienting factor in adaptive evolution, influencing the mutational trajectories populations follow toward higher-fitness genotypes. Because these trajectories are typically just a small subset of all possible mutational trajectories, they can exhibit composition bias—an enrichment of a particular kind of DNA sequence change, such as transition or transversion mutations. Here, we use empirical data from eukaryotic transcriptional regulation to study how mutation bias and composition bias interact to influence adaptive evolution.

Genetic Variants and Phenotypic Characteristics of Salmonella Typhimurium-Resistant Mutants after Exposure to Carvacrol

The emergence of antimicrobial resistance has raised questions about the safety of essential oils and their individual constituents as food preservatives and as disinfection agents. Further research is required to understand how and under what conditions stable genotypic resistance might occur in food pathogens. Evolution experiments on Salmonella Typhimurium cyclically exposed to sublethal and lethal doses of carvacrol permitted the isolation of SeSCar and SeLCar strains, respectively. Both evolved strains showed a significant increase in carvacrol resistance, assessed by minimum inhibitory and bactericidal concentrations, the study of growth kinetics in the presence of carvacrol, and the evaluation of survival under lethal conditions. Moreover, antibiotic susceptibility tests revealed a development of SeLCar resistance to a wide range of antibiotics. Whole genome sequencing allowed the identification of single nucleotide variations in transcriptional regulators of oxidative stress-response: yfhP in SeSCar and soxR in SeLCar, which could be responsible for the increased resistance by improving the response to carvacrol and preventing its accumulation inside the cell. This study demonstrates the emergence of S. Typhimurium-resistant mutants against carvacrol, which might pose a risk to food safety and should therefore be considered in the design of food preservation strategies, or of cleaning and disinfection treatments.

Irregularities in genetic variation and mutation rates with environmental stresses

The appearance of new mutations is determined by the equilibrium between DNA error formation and repair. In bacteria like Escherichia coli, stresses are thought shift this balance towards increased mutagenesis. Recent findings, however, suggest a very uneven relationship between stress and mutations. Only a subset of stressful environments increase the net rate of mutation and different forms of nutritional stress (such as oxygen, carbon or phosphorus limitations) result in markedly different mutation rates after similar reductions in growth rate. Moreover, different stresses result in altered mutational spectra, with some increasing transposition and others increasing indel formation. Single-base substitution rates are lower with some stresses than in unstressed bacteria. Indeed, changes to the mix of mutations with stress are more widespread than a marked increase in net mutation rate. Much remains to be learned on how environments have unique mutational signatures and why some stresses are more mutagenic than others. Even beyond stress-induced genetic variation, the fundamental unresolved question in the stress-mutation relationship is the adaptive value of different types of mutations and mutation rates; is transposition, for example, more advantageous under anaerobic conditions? It remains to be investigated whether stress-specific genetic variation impacts on evolvability differentially in distinct environments.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis

Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen.

Some types of mutations occur more frequently than expected. This study shows that such bias —an excess of transitions over transversions—influences the evolution of antibiotic resistance in a key global pathogen, Mycobacterium tuberculosis.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Spatial Vulnerabilities of the Escherichia coli Genome to Spontaneous Mutations Revealed with Improved Duplex Sequencing

Investigation of spontaneous mutations by next-generation sequencing technology has attracted extensive attention lately due to the fundamental roles of spontaneous mutations in evolution and pathological processes. However, these studies only focused on the mutations accumulated through many generations during long-term (possibly be years of) culturing, but not the freshly generated mutations that occur at very low frequencies. In this study, we established a molecularly barcoded deep sequencing strategy to detect low abundant spontaneous mutations in genomes of bacteria cell cultures. Genome-wide spontaneous mutations in 15 Escherichia coli cell culture samples were defined with a high confidence (P < 0.01). We also developed a hotspot-calling approach based on the run-length encoding algorithm to find the genomic regions that are vulnerable to the spontaneous mutations. The hotspots for the mutations appeared to be highly conserved across the bacteria samples. Further biological annotation of these regions indicated that most of the spontaneous mutations were located at the repeat domains or nonfunctional domains of the genomes, suggesting the existence of mechanisms that could somehow prevent the occurrence of mutations in crucial genic areas. This study provides a more faithful picture of mutation occurrence and spectra in a single expansion process without long-term culturing.

Mutational signatures indicative of environmental stress in bacteria

Evolutionary innovations are dependent on mutations. Mutation rates are increased by adverse conditions in the laboratory, but there is no evidence that stressful environments that do not directly impact on DNA leave a mutational imprint on extant genomes. Mutational spectra in the laboratory are normally determined with unstressed cells but are unavailable with stressed bacteria. To by-pass problems with viability, selection effects, and growth rate differences due to stressful environments, in this study we used a set of genetically engineered strains to identify the mutational spectrum associated with nutritional stress. The strain set members each had a fixed level of the master regulator protein, RpoS, which controls the general stress response of Escherichia coli. By assessing mutations in cycA gene from 485 cycloserine resistant mutants collected from as many independent cultures with three distinct perceived stress (RpoS) levels, we were able establish a dose-dependent relationship between stress and mutational spectra. The altered mutational patterns included base pair substitutions, single base pair indels, longer indels, and transpositions of different insertion sequences. The mutational spectrum of low-RpoS cells closely matches the genome-wide spectrum previously generated in laboratory environments, while the spectra of high RpoS, high perceived stress cells more closely matches spectra found in comparisons of extant genomes. Our results offer an explanation of the uneven mutational profiles such as the transition-transversion biases observed in extant genomes and provide a framework for assessing the contribution of stress-induced mutagenesis to evolutionary transitions and the mutational emergence of antibiotic resistance and disease states.

Transcription-associated mutagenesis

Transcription requires unwinding complementary DNA strands, generating torsional stress, and sensitizing the exposed single strands to chemical reactions and endogenous damaging agents. In addition, transcription can occur concomitantly with the other major DNA metabolic processes (replication, repair, and recombination), creating opportunities for either cooperation or conflict. Genetic modifications associated with transcription are a global issue in the small genomes of microorganisms in which noncoding sequences are rare. Transcription likewise becomes significant when one considers that most of the human genome is transcriptionally active. In this review, we focus specifically on the mutagenic consequences of transcription. Mechanisms of transcription-associated mutagenesis in microorganisms are discussed, as is the role of transcription in somatic instability of the vertebrate immune system.

High incidence of multiple antibiotic resistant cells in cultures of in enterohemorrhagic Escherichia coli O157:H7

The spontaneous incidence of chloramphenicol (Cam) resistant mutant bacteria is at least ten-fold higher in cultures of enterohemorrhagic Escherichia coli O157:H7 strain EDL933 than in E. coli K-12. It is at least 100-fold higher in the dam (DNA adenine methyltransferase) derivative of EDL933, compared to the dam strain of E. coli K-12, thereby preventing the use of Cam resistance as a marker in gene replacement technology. Genome sequencing of Cam-resistant isolates of EDL933 and its dam derivatives showed that the marR (multiple antibiotic resistance) gene was mutated in every case but not in the Cam-sensitive parental strains. As expected from mutation in the marR gene, the Cam-resistant bacteria were also found to be resistant to tetracycline and nalidixic acid. The marR gene in strain EDL933 is annotated as a shorter open reading frame than that in E. coli K-12 but the longer marR(+) open reading frame was more efficient at complementing the marR antibiotic-resistance phenotype of strain EDL933. Beta-lactamase-tolerant derivatives were present at frequencies 10-100 times greater in cultures of marR derivatives of strain EDL933 than the parent strain. Spontaneous mutation frequency to rifampicin, spectinomycin and streptomycin resistance was the same in E. coli O157:H7 and E. coli K-12 strains.

On the mutational topology of the bacterial genome

By sequencing the genomes of 34 mutation accumulation lines of a mismatch-repair defective strain of Escherichia coli that had undergone a total of 12,750 generations, we identified 1625 spontaneous base-pair substitutions spread across the E. coli genome. These mutations are not distributed at random but, instead, fall into a wave-like spatial pattern that is repeated almost exactly in mirror image in the two separately replicated halves of the bacterial chromosome. The pattern is correlated to genomic features, with mutation densities greatest in regions predicted to have high superhelicity. Superimposed upon this pattern are regional hotspots, some of which are located where replication forks may collide or be blocked. These results suggest that, as they traverse the chromosome, the two replication forks encounter parallel structural features that change the fidelity of DNA replication.

Signal correlations in ecological niches can shape the organization and evolution of bacterial gene regulatory networks

Transcriptional regulation plays a significant role in the biological response of bacteria to changing environmental conditions. Therefore, mapping transcriptional regulatory networks is an important step not only in understanding how bacteria sense and interpret their environment but also to identify the functions involved in biological responses to specific conditions. Recent experimental and computational developments have facilitated the characterization of regulatory networks on a genome-wide scale in model organisms. In addition, the multiplication of complete genome sequences has encouraged comparative analyses to detect conserved regulatory elements and infer regulatory networks in other less well-studied organisms. However, transcription regulation appears to evolve rapidly, thus, creating challenges for the transfer of knowledge to nonmodel organisms. Nevertheless, the mechanisms and constraints driving the evolution of regulatory networks have been the subjects of numerous analyses, and several models have been proposed. Overall, the contributions of mutations, recombination, and horizontal gene transfer are complex. Finally, the rapid evolution of regulatory networks plays a significant role in the remarkable capacity of bacteria to adapt to new or changing environments. Conversely, the characteristics of environmental niches determine the selective pressures and can shape the structure of regulatory network accordingly.

High levels of transcription stimulate transversions at GC base pairs in yeast

High-levels of transcription through a gene stimulate spontaneous mutation rate, a phenomenon termed transcription-associated mutation (TAM). While transcriptional effects on specific mutation classes have been identified using forward mutation and frameshift-reversion assays, little is yet known about transcription-associated base substitutions in yeast. To address this issue, we developed a new base substitution reversion assay (the lys2-TAG allele). We report a 22-fold increase in overall reversion rate in the high- relative to the low-transcription strain (from 2.1- to 47- × 10(-9) ). While all detectable base substitution types increased in the high-transcription strain, G→T and G→C transversions increased disproportionately by 58- and 52-fold, respectively. To assess a potential role of DNA damage in the TAM events, we measured mutation rates and spectra in individual strains defective in the repair of specific DNA lesions or null for the error-prone translesion DNA polymerase zeta (Pol zeta). Results exclude a role of 8-oxoGuanine, general oxidative damage, or apurinic/apyrimidinic sites in the generation of TAM G→T and G→C transversions. In contrast, the TAM transversions at GC base pairs depend on Pol zeta for occurrence implicating DNA damage, other than oxidative lesions or AP sites, in the TAM mechanism. Results further indicate that transcription-dependent G→T transversions in yeast differ mechanistically from equivalent events in E. coli reported by others. Given their occurrences in repair-proficient cells, transcription-associated G→T and G→C events represent a novel type of transcription-associated mutagenesis in normal cells with potentially important implications for evolution and genetic disease.

Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing

Knowledge of the rate and nature of spontaneous mutation is fundamental to understanding evolutionary and molecular processes. In this report, we analyze spontaneous mutations accumulated over thousands of generations by wild-type Escherichia coli and a derivative defective in mismatch repair (MMR), the primary pathway for correcting replication errors. The major conclusions are (i) the mutation rate of a wild-type E. coli strain is ~1 × 10(-3) per genome per generation; (ii) mutations in the wild-type strain have the expected mutational bias for G:C > A:T mutations, but the bias changes to A:T > G:C mutations in the absence of MMR; (iii) during replication, A:T > G:C transitions preferentially occur with A templating the lagging strand and T templating the leading strand, whereas G:C > A:T transitions preferentially occur with C templating the lagging strand and G templating the leading strand; (iv) there is a strong bias for transition mutations to occur at 5'ApC3'/3'TpG5' sites (where bases 5'A and 3'T are mutated) and, to a lesser extent, at 5'GpC3'/3'CpG5' sites (where bases 5'G and 3'C are mutated); (v) although the rate of small (≤4 nt) insertions and deletions is high at repeat sequences, these events occur at only 1/10th the genomic rate of base-pair substitutions. MMR activity is genetically regulated, and bacteria isolated from nature often lack MMR capacity, suggesting that modulation of MMR can be adaptive. Thus, comparing results from the wild-type and MMR-defective strains may lead to a deeper understanding of factors that determine mutation rates and spectra, how these factors may differ among organisms, and how they may be shaped by environmental conditions.

Transcription as a source of genome instability

Alterations in genome sequence and structure contribute to somatic disease, affect the fitness of subsequent generations and drive evolutionary processes. The crucial roles of highly accurate replication and efficient repair in maintaining overall genome integrity are well-known, but the more localized stability costs that are associated with transcribing DNA into RNA molecules are less appreciated. Here we review the diverse ways in which the essential process of transcription alters the underlying DNA template and thereby modifies the genetic landscape.

Nucleotide and phylogenetic analyses of the Chlamydia trachomatis ompA gene indicates it is a hotspot for mutation

Background

Serovars of the human pathogen Chlamydia trachomatis occupy one of three specific tissue niches. Genomic analyses indicate that the serovars have a phylogeny congruent with their pathobiology and have an average substitution rate of less than one nucleotide per kilobase. In contrast, the gene that determines serovar specificity, ompA, has a phylogenetic association that is not congruent with tissue tropism and has a degree of nucleotide variability much higher than other genomic loci. The ompA gene encodes the major surface-exposed antigenic determinant, and the observed nucleotide diversity at the ompA locus is thought to be due to recombination and host immune selection pressure. The possible contribution of a localized increase in mutation rate, however, has not been investigated.

Results

Nucleotide diversity and phylogenetic relationships of the five constant and four variable domains of the ompA gene, as well as several loci surrounding ompA, were examined for each serovar. The loci flanking the ompA gene demonstrated that nucleotide diversity increased monotonically as ompA is approached and that their gene trees are not congruent with either ompA or tissue tropism. The variable domains of the ompA gene had a very high level of non-synonymous change, which is expected as these regions encode the surface-exposed epitopes and are under positive selection. However, the synonymous changes are clustered in the variable regions compared to the constant domains; if hitchhiking were to account for the increase in synonymous changes, these substitutions should be more evenly distributed across the gene. Recombination also cannot entirely account for this increase as the phylogenetic relationships of the constant and variable domains are congruent with each other.

Conclusions

The high number of synonymous substitutions observed within the variable domains of ompA appears to be due to an increased mutation rate within this region of the genome, whereas the increase in nucleotide substitution rate and the lack of phylogenetic congruence in the regions flanking ompA are characteristic motifs of gene conversion. Together, the increased mutation rate in the ompA gene, in conjunction with gene conversion and positive selection, results in a high degree of variability that promotes host immune evasion.

Transcriptional de-repression and Mfd are mutagenic in stressed Bacillus subtilis cells

In recent years, it has been proposed that conflicts between transcription and active chromosomal replication engender genome instability events. Furthermore, transcription elongation factors have been reported to prevent conflicts between transcription and replication and avoid genome instability. Here, we examined transcriptional de-repression as a genetic diversity-producing agent and showed, through the use of physiological and genetic means, that transcriptional de-represssion of a leuC defective allele leads to accumulation of Leu(+) mutations. We also showed, by using riboswitches that activate transcription in conditions of tyrosine or methionine starvation, that the effect of transcriptional de-repression of the leuC construct on the accumulation of Leu(+) mutations was independent of selection. We examined the role of Mfd, a transcription elongation factor involved in DNA repair, in this process and showed that proficiency of this factor promotes mutagenic events. These results are in stark contrast to previous reports in Escherichia coli, which showed that Mfd prevents replication fork collapses. Because our assays place cells under non-growing conditions, by starving them for two amino acids, we surmised that the Mfd mutagenic process associated with transcriptional de-repression does not result from conflicts with chromosomal replication. These results raise the interesting concept that transcription elongation factors may serve two functions in cells. In growing conditions, these factors prevent the generation of mutations, while in stress or non-growing conditions they mediate the production of genetic diversity.

Changes in transcriptional orientation are associated with increases in evolutionary rates of enterobacterial genes

Background

Changes in transcriptional orientation (“CTOs”) occur frequently in prokaryotic genomes. Such changes usually result from genomic inversions, which may cause a conflict between the directions of replication and transcription and an increase in mutation rate. However, CTOs do not always lead to the replication-transcription confrontation. Furthermore, CTOs may cause deleterious disruptions of operon structure and/or gene regulations. The currently existing CTOs may indicate relaxation of selection pressure. Therefore, it is of interest to investigate whether CTOs have an independent effect on the evolutionary rates of the affected genes, and whether these genes are subject to any type of selection pressure in prokaryotes.

Methods

Three closely related enterbacteria, Escherichia coli, Klebsiella pneumoniae and Salmonella enterica serovar Typhimurium, were selected for comparisons of synonymous (dS) and nonsynonymous (dN) substitution rate between the genes that have experienced changes in transcriptional orientation (changed-orientation genes, “COGs”) and those that do not (same-orientation genes, “SOGs”). The dN/dS ratio was also derived to evaluate the selection pressure on the analyzed genes. Confounding factors in the estimation of evolutionary rates, such as gene essentiality, gene expression level, replication-transcription confrontation, and decreased dS at gene terminals were controlled in the COG-SOG comparisons.

Results

We demonstrate that COGs have significantly higher dN and dS than SOGs when a series of confounding factors are controlled. However, the dN/dS ratios are similar between the two gene groups, suggesting that the increase in dS can sufficiently explain the increase in dN in COGs. Therefore, the increases in evolutionary rates in COGs may be mainly mutation-driven.

Conclusions

Here we show that CTOs can increase the evolutionary rates of the affected genes. This effect is independent of the replication-transcription confrontation, which is suggested to be the major cause of inversion-associated evolutionary rate increases. The real cause of such evolutionary rate increases remains unclear but is worth further explorations.

Evidence that mutation is universally biased towards AT in bacteria

Mutation is the engine that drives evolution and adaptation forward in that it generates the variation on which natural selection acts. Mutation is a random process that nevertheless occurs according to certain biases. Elucidating mutational biases and the way they vary across species and within genomes is crucial to understanding evolution and adaptation. Here we demonstrate that clonal pathogens that evolve under severely relaxed selection are uniquely suitable for studying mutational biases in bacteria. We estimate mutational patterns using sequence datasets from five such clonal pathogens belonging to four diverse bacterial clades that span most of the range of genomic nucleotide content. We demonstrate that across different types of sites and in all four clades mutation is consistently biased towards AT. This is true even in clades that have high genomic GC content. In all studied cases the mutational bias towards AT is primarily due to the high rate of C/G to T/A transitions. These results suggest that bacterial mutational biases are far less variable than previously thought. They further demonstrate that variation in nucleotide content cannot stem entirely from variation in mutational biases and that natural selection and/or a natural selection-like process such as biased gene conversion strongly affect nucleotide content.

Transcription-associated mutagenesis increases protein sequence diversity more effectively than does random mutagenesis in Escherichia coli

Background

During transcription, the nontranscribed DNA strand becomes single-stranded DNA (ssDNA), which can form secondary structures. Unpaired bases in the ssDNA are less protected from mutagens and hence experience more mutations than do paired bases. These mutations are called transcription-associated mutations. Transcription-associated mutagenesis is increased under stress and depends on the DNA sequence. Therefore, selection might significantly influence protein-coding sequences in terms of the transcription-associated mutability per transcription event under stress to improve the survival of Escherichia coli.

Methodology/Principal Findings

The mutability index (MI) was developed by Wright et al. to estimate the relative transcription-associated mutability of bases per transcription event. Using the most stable fold of each ssDNA that have an average length n, MI was defined as (the number of folds in which the base is unpaired)/n×(highest –ΔG of all n folds in which the base is unpaired), where ΔG is the free energy. The MI values show a significant correlation with mutation data under stress but not with spontaneous mutations in E. coli. Protein sequence diversity is preferred under stress but not under favorable conditions. Therefore, we evaluated the selection pressure on MI in terms of the protein sequence diversity for all the protein-coding sequences in E. coli. The distributions of the MI values were lower at bases that could be substituted with each of the other three bases without affecting the amino acid sequence than at bases that could not be so substituted. Start codons had lower distributions of MI values than did nonstart codons.

Conclusions/Significance

Our results suggest that the majority of protein-coding sequences have evolved to promote protein sequence diversity and to reduce gene knockout under stress. Consequently, transcription-associated mutagenesis increases protein sequence diversity more effectively than does random mutagenesis under stress. Nonrandom transcription-associated mutagenesis under stress should improve the survival of E. coli.

Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria

Background

Because bacteria do not have a robust fossil record, attempts to infer the timing of events in their evolutionary history requires comparisons of molecular sequences. This use of molecular clocks is based on the assumptions that substitution rates for homologous genes or sites are fairly constant through time and across taxa. Violation of these conditions can lead to erroneous inferences and result in estimates that are off by orders of magnitude. In this study, we examine the consistency of substitution rates among a set of conserved genes in diverse bacterial lineages, and address the questions regarding the validity of molecular dating.

Results

By examining the evolution of 16S rRNA gene in obligate endosymbionts, which can be calibrated by the fossil record of their hosts, we found that the rates are consistent within a clade but varied widely across different bacterial lineages. Genome-wide estimates of nonsynonymous and synonymous substitutions suggest that these two measures are highly variable in their rates across bacterial taxa. Genetic drift plays a fundamental role in determining the accumulation of substitutions in 16S rRNA genes and at nonsynonymous sites. Moreover, divergence estimates based on a set of universally conserved protein-coding genes also exhibit low correspondence to those based on 16S rRNA genes.

Conclusion

Our results document a wide range of substitution rates across genes and bacterial taxa. This high level of variation cautions against the assumption of a universal molecular clock for inferring divergence times in bacteria. However, by applying relative-rate tests to homologous genes, it is possible to derive reliable local clocks that can be used to calibrate bacterial evolution.

Reviewers

This article was reviewed by Adam Eyre-Walker, Simonetta Gribaldo and Tal Pupko (nominated by Dan Graur).

Whole-genome mutational biases in bacteria

A fundamental biological question is what forces shape the guanine plus cytosine (GC) content of genomes. We studied the specificity and rate of different mutational biases in real time in the bacterium Salmonella typhimurium under conditions of strongly reduced selection and in the absence of the major DNA repair systems involved in repairing common spontaneous mutations caused by oxidized and deaminated DNA bases. The mutational spectrum was determined by whole-genome sequencing of two S. typhimurium mutants that were serially passaged for 5,000 generations. Analysis of 943 identified base pair substitutions showed that 91% were GC-to-TA transversions and 7% were GC-to-AT transitions, commonly associated with 8-oxoG- and deamination-induced damages, respectively. Other types of base pair substitutions constituted the remaining 2% of the mutations. With regard to mutational biases, there was a significant increase in C-to-T transitions on the nontranscribed strand, and for highly expressed genes, C/G-to-T mutations were more common than expected; however, no significant mutational bias with regard to leading and lagging strands of replication or chromosome position were found. These results suggest that, based on the experimentally determined mutational rates and specificities, a bacterial genome lacking the relevant DNA repair systems could, as a consequence of these underlying mutational biases, very rapidly reduce its GC content.

Selection shapes malaria genomes and drives divergence between pathogens infecting hominids versus rodents

BACKGROUND

Malaria kills more people worldwide than all inherited human genetic disorders combined. To characterize how the parasites causing this disease adapt to different host environments, we compared the evolutionary genomics of two distinct groups of malaria pathogens in order to identify critical properties associated with infection of different hosts: those parasites infecting hominids (Plasmodium falciparum and P. reichenowi) versus parasites infecting rodent hosts (P. yoelii yoelii, P. berghei, and P. chabaudi). Adaptation by the parasite to its host is likely highly critical to the evolution of these species.

RESULTS

Our comparative analysis suggests that patterns of molecular evolution in the hominid parasite lineage are generally similar to those of the rodent lineage but distinct in several aspects. The most rapidly evolving genes in both lineages are those involved in host-parasite interactions as well as those that show the lowest expression levels. However, we found that, similar to their respective mammal host lineages, parasite genomes infecting hominids are generally less constrained, evolving at faster rates, and accumulating more deleterious mutations than those infecting murids, which may reflect an historical lower effective size of the hominid lineage and relaxed host-driven selective pressures.

CONCLUSION

Our study highlights for the first time the differences in trends and rates of evolution in Plasmodium lineages infecting different hosts and emphasizes the potential importance of the variation in effective size between lineages to explain variation in selective constraints among genomes.

Similar compositional biases are caused by very different mutational effects

Compositional replication strand bias, commonly referred to as GC skew, is present in many genomes of prokaryotes, eukaryotes, and viruses. Although cytosine deamination in ssDNA (resulting in C-->T changes on the leading strand) is often invoked as its major cause, the precise contributions of this and other substitution types are currently unknown. It is also unclear if the underlying mutational asymmetries are the same among taxa, are stable over time, or how closely the observed biases are to mutational equilibrium. We analyzed nearly neutral sites of seven taxa each with between three and six complete bacterial genomes, and inferred the substitution spectra of fourfold degenerate positions in nonhighly expressed genes. Using a bootstrap procedure, we extracted compositional biases associated with replication and identified the significant asymmetries. Although all taxa showed an overrepresentation of G relative to C on the leading strand (and imbalances between A and T), widely variable substitution asymmetries are noted. Surprisingly, all substitution types show significant asymmetry in at least one taxon, but none were universally biased in all taxa. Notably, in the two most biased genomes, A-->G, rather than C-->T, shapes the compositional bias. Given the variability in these biases, we propose that the process is multifactorial. Finally, we also find that most genomes are not at compositional equilibrium, and suggest that mutational-based heterotachy is deeply imprinted in the history of biological macromolecules. This shows that similar compositional biases associated with the same essential well-conserved process, replication, do not reflect similar mutational processes in different genomes, and that caution is required in inferring the roles of specific mutational biases on the basis of contemporary patterns of sequence composition.

The ompA gene in Chlamydia trachomatis differs in phylogeny and rate of evolution from other regions of the genome

Strains of Chlamydia trachomatis are classified into serovars based on nucleotide sequence differences in ompA, the gene that encodes the major outer membrane protein. Phylogenetic characterization of strains based on ompA, however, results in serovar groupings that are inconsistent with the distinguishing features of C. trachomatis pathobiology, e.g., tissue tropisms and disease presentation. We have compared nucleotide sequences at multiple sites distributed around the chlamydial genome from 18 strains representing 16 serovars; sampled regions included genes encoding housekeeping enzymes (totaling 2,073 bp), intergenic noncoding segments (1,612 bp), and a gene encoding a second outer membrane protein (porB; 1,023 bp), with the ompA sequence (1,194 bp) used for reference. These comparative analyses revealed substantial variation in nucleotide substitution patterns among the sampled regions, with average pairwise sequence differences ranging from 0.15% for the housekeeping genes to 12.1% for ompA. Phylogenetic characterization of the sampled genomic sequences yielded a strongly supported tree that divides the strains into groupings consistent with C. trachomatis biology and which has a topology quite distinct from the ompA tree. This phylogenetic incongruity can be accounted for by recombination of the ompA gene between different genomic backgrounds. We found, however, no evidence of recombination within or between any of the sampled regions around the C. trachomatis genome apart from ompA. Parallel analysis of published sequence data on four members of the pmp gene family are consistent with the phylogenetic analyses reported here.

Transcription promotes guanine to thymine mutations in the non-transcribed strand of an Escherichia coli gene

Transcription of DNA opens the chromatin, causes topological changes in DNA and transiently exposes the two strands to different biochemical environments. Consequently, it has long been argued that transcription may promote damage to DNA and there are data in Escherichia coli and yeast supporting a correlation between high transcription and mutations. We examined the transcription-dependence of the reversion of a nonsense codon (TGA) in E. coli and found that there was a strong dependence of mutations on transcription in strains defective in the repair of 8-oxoguanine in DNA. Under conditions of high transcription there was a three to five-fold increase in mutations that changed TGA in the non-transcribed strand to a sense codon. Furthermore, in both mutY and mutM mutY backgrounds the mutations were overwhelmingly G:C to T:A. In contrast, when the TGA was in the transcribed strand in relation with the inducible promoter, high transcription decreased the rate of reversion. Similar results were obtained in a strain defective in the transcription-repair coupling factor, Mfd, suggesting that transcription dependent increase in base substitutions does not require transcription-dependent DNA repair. However, Mfd does modulate the magnitude of the mutagenic effect of transcription. These data are consistent with a model in which the non-transcribed strand is more susceptible to oxidative damage during transcription than the transcribed strand. These results suggest that the magnitudes of individual base substitutions and their relative numbers in other studies of mutational spectra may also be affected by transcription.

Natural antisense LHCGR could make sense of hypogonadism, male-limited precocious puberty and pre-eclampsia

The pleiotropic effects of human chorionic gonadotrophin (hCG), the key regulator of human pregnancy, are dependent upon cell surface expression of its functional cognate receptor LHCGR in the placental trophoblasts, corpus luteum, uterus, vascular endothelial and smooth muscle cells. Additionally, lutenizing hormone-mediated signalling failure has often been linked to activating/inactivating mutations in LHCGR. One of the intriguing aspects of these studies is that the mutations are most frequently located within C-terminal 200-350 residues of the receptor protein. In an attempt to reconcile the mechanistic basis of LHCGR regulation and mutations, we have carried out bioinformatic analyses to identify the CpG-rich regions and the major potential scaffold/matrix attachment sites (S/MARs) in LHCGR and neighbouring gene (ALF) at human chromosome 2p21. Based on these analyses, we propose a chromatin-loop model, which may explain the temporal regulation and susceptibility to mutation of the human LHCGR. One of the characteristic features of the model, is that the major potential S/MAR sequences of the human LHCGR gene (68 kb) are located at the 3' end of the gene, and unlike mouse, the transmembrane and C-terminal protein coding sequences at exon 11 are embedded in this S/MAR site. Moreover, this region is subject to antisense transcription from the neighbouring gene ALF, which is gonad-specific and is only activated in meiotic spermatocytes and oocytes. Together, these analyses suggest that exon 11 of human LHCGR could be more susceptible to mutation than the other 10 exons together and that activation of LHCGR, contingent to the somatic silencing of neighbouring ALF, could be linked to male-limited precocious puberty and pre-eclampsia.

Adjusting for selection on synonymous sites in estimates of evolutionary distance

Evolution at silent sites is often used to estimate the pace of selectively neutral processes or to infer differences in divergence times of genes. However, silent sites are subject to selection in favor of preferred codons, and the strength of such selection varies dramatically across genes. Here, we use the relationship between codon bias and synonymous divergence observed in four species of the genus Saccharomyces to provide a simple correction for selection on silent sites.