Nucleoid-associated proteins affect mutation dynamics in E. coli in a growth phase-specific manner

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

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

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

Mutation bias and adaptation in bacteria

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

Genomic inversions in Escherichia coli alter gene expression and are associated with nucleoid protein binding sites

Genomic reorganization, like rearrangements and inversions, influence how genetic information is organized within bacterial genomes. Inversions in particular, facilitate genome evolution through gene gain and loss, and can alter gene expression. Previous studies investigating the impact inversions have on gene expression induced inversions targeting specific genes or examine inversions between distantly related species. This fails to encompass a genome wide perspective on naturally occurring inversions and their post adaptation impact on gene expression. Here we use bioinformatic techniques and multiple RNA-seq datasets to investigate the short- and long-range impact inversions have on genomic gene expression within <i>Escherichia coli</i>. We observed differences in gene expression between homologous inverted and non-inverted genes, even after long term exposure to adaptive selection. In 4% of inversions representing 33 genes, differential gene expression between inverted and non-inverted homologs was detected, with nearly two thirds (71%) of differentially expressed inverted genes having 9.4-85.6 fold higher gene expression. The identified inversions had more overlap than expected with nucleoid associated protein binding sites, which assist in genomic gene expression regulation. Some inversions can drastically impact gene expression even between different strains of <i>E.coli</i>, and could provide a mechanism for the diversification of genetic content through controlled expression changes.

Special Issue: Role of Bacterial Chromatin in Environmental Sensing, Adaptation and Evolution

A typical bacterial cell is micron-sized and contains a genome several million base pairs in length [...].

Xenogeneic Silencing and Bacterial Genome Evolution: Mechanisms for DNA Recognition Imply Multifaceted Roles of Xenogeneic Silencers

Horizontal gene transfer (HGT) is a major driving force for bacterial evolution. To avoid the deleterious effects due to the unregulated expression of newly acquired foreign genes, bacteria have evolved specific proteins named xenogeneic silencers to recognize foreign DNA sequences and suppress their transcription. As there is considerable diversity in genomic base compositions among bacteria, how xenogeneic silencers distinguish self- from nonself DNA in different bacteria remains poorly understood. This review summarizes the progress in studying the DNA binding preferences and the underlying molecular mechanisms of known xenogeneic silencer families, represented by H-NS of Escherichia coli, Lsr2 of Mycobacterium, MvaT of Pseudomonas, and Rok of Bacillus. Comparative analyses of the published data indicate that the differences in DNA recognition mechanisms enable these xenogeneic silencers to have clear characteristics in DNA sequence preferences, which are further correlated with different host genomic features. These correlations provide insights into the mechanisms of how these xenogeneic silencers selectively target foreign DNA in different genomic backgrounds. Furthermore, it is revealed that the genomic AT contents of bacterial species with the same xenogeneic silencer family proteins are distributed in a limited range and are generally lower than those species without any known xenogeneic silencers in the same phylum/class/genus, indicating that xenogeneic silencers have multifaceted roles on bacterial genome evolution. In addition to regulating horizontal gene transfer, xenogeneic silencers also act as a selective force against the GC to AT mutational bias found in bacterial genomes and help the host genomic AT contents maintained at relatively low levels.

Mutation and Recombination Rates Vary Across Bacterial Chromosome

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

Integration Host Factor IHF facilitates homologous recombination and mutagenic processes in Pseudomonas putida

Nucleoid-associated proteins (NAPs) such as IHF, HU, Fis, and H-NS alter the topology of bound DNA and may thereby affect accessibility of DNA to repair and recombination processes. To examine this possibility, we investigated the effect of IHF on the frequency of homologous recombination (HR) and point mutations in soil bacterium Pseudomonas putida by using plasmidial and chromosomal assays. We observed positive effect of IHF on the frequency of HR, whereas this effect varied depending both on the chromosomal location of the HR target and the type of plasmid used in the assay. The occurrence of point mutations in plasmid was also facilitated by IHF, whereas in the chromosome the positive effect of IHF appeared only at certain DNA sequences and/or chromosomal positions. We did not observe any significant effects of IHF on the spectrum of mutations. However, despite of the presence or absence of IHF, different mutational hot spots appeared both in plasmid and in chromosome. Additionally, the frequency of frameshift mutations in the chromosome was also strongly affected by the location of the mutational target sequence. Taking together, our results indicate that IHF facilitates the occurrence of genetic changes in P. putida, whereas the location of the target sequence affects both the IHF-dependent and IHF-independent mechanisms.

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia. Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.

That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

Considerations on bacterial nucleoids

The classic genome organization of the bacterial chromosome is normally envisaged with all its genetic markers linked, thus forming a closed genetic circle of duplex stranded DNA (dsDNA) and several proteins in what it is called as "the bacterial nucleoid." This structure may be more or less corrugated depending on the physiological state of the bacterium (i.e., resting state or active growth) and is not surrounded by a double membrane as in eukayotic cells. The universality of the closed circle model in bacteria is however slowly changing, as new data emerge in different bacterial groups such as in Planctomycetes and related microorganisms, species of Borrelia, Streptomyces, Agrobacterium, or Phytoplasma. In these and possibly other microorganisms, the existence of complex formations of intracellular membranes or linear chromosomes is typical; all of these situations contributing to weakening the current cellular organization paradigm, i.e., prokaryotic vs eukaryotic cells.

A response regulator promotes Francisella tularensis intramacrophage growth by repressing an anti-virulence factor

The orphan response regulator PmrA is essential for the intramacrophage growth and survival of Francisella tularensis. PmrA was thought to promote intramacrophage growth by binding directly to promoters on the Francisella Pathogenicity Island (FPI) and positively regulating the expression of FPI genes, which encode a Type VI secretion system required for intramacrophage growth. Using both ChIP-Seq and RNA-Seq we identify those regions of the F. tularensis chromosome occupied by PmrA and those genes that are regulated by PmrA. We find that PmrA associates with 252 distinct regions of the F. tularensis chromosome, but exerts regulatory effects at only a few of these locations. Rather than by functioning directly as an activator of FPI gene expression we present evidence that PmrA promotes intramacrophage growth by repressing the expression of a single target gene we refer to as priM (PmrA-repressed inhibitor of intramacrophage growth). Our findings thus indicate that the role of PmrA in facilitating intracellular growth is to repress a previously unknown anti-virulence factor. PriM is the first bacterially encoded factor to be described that can interfere with the intramacrophage growth and survival of F. tularensis.

Bacterial chromatin: converging views at different scales

Bacterial genomes are functionally organized and compactly folded into a structure referred to as bacterial chromatin or the nucleoid. An important role in genome folding is attributed to Nucleoid-Associated Proteins, also referred to as bacterial chromatin proteins. Although a lot of molecular insight in the mechanisms of operation of these proteins has been generated in the test tube, knowledge on genome organization in the cellular context is still lagging behind severely. Here, we discuss important advances in the understanding of three-dimensional genome organization due to the application of Chromosome Conformation Capture and super-resolution microscopy techniques. We focus on bacterial chromatin proteins whose proposed role in genome organization is supported by these approaches. Moreover, we discuss recent insights into the interrelationship between genome organization and genome activity/stability in bacteria.

H-NS Facilitates Sequence Diversification of Horizontally Transferred DNAs during Their Integration in Host Chromosomes

Bacteria can acquire new traits through horizontal gene transfer. Inappropriate expression of transferred genes, however, can disrupt the physiology of the host bacteria. To reduce this risk, Escherichia coli expresses the nucleoid-associated protein, H-NS, which preferentially binds to horizontally transferred genes to control their expression. Once expression is optimized, the horizontally transferred genes may actually contribute to E. coli survival in new habitats. Therefore, we investigated whether and how H-NS contributes to this optimization process. A comparison of H-NS binding profiles on common chromosomal segments of three E. coli strains belonging to different phylogenetic groups indicated that the positions of H-NS-bound regions have been conserved in E. coli strains. The sequences of the H-NS-bound regions appear to have diverged more so than H-NS-unbound regions only when H-NS-bound regions are located upstream or in coding regions of genes. Because these regions generally contain regulatory elements for gene expression, sequence divergence in these regions may be associated with alteration of gene expression. Indeed, nucleotide substitutions in H-NS-bound regions of the ybdO promoter and coding regions have diversified the potential for H-NS-independent negative regulation among E. coli strains. The ybdO expression in these strains was still negatively regulated by H-NS, which reduced the effect of H-NS-independent regulation under normal growth conditions. Hence, we propose that, during E. coli evolution, the conservation of H-NS binding sites resulted in the diversification of the regulation of horizontally transferred genes, which may have facilitated E. coli adaptation to new ecological niches.

Horizontal gene transfer among bacteria is the major means of acquiring genetic diversity and has been a central factor in bacterial evolution. The expression of horizontally transferred genes could potentially be optimized to permit the host bacteria to expand their habitat. The results of our study suggest that DNA regions bound by the nucleoid-associated protein, H-NS, which preferentially binds to horizontally transferred genes, have been conserved during Escherichia coli evolution. Interestingly, H-NS-bound regions have evolved faster than H-NS-unbound regions, but only in gene regulatory and coding regions. We show that DNA sequence substitutions in H-NS-bound regions actually alter the regulation of gene expression in different E. coli strains. Thus, our results support the hypothesis that H-NS accelerates the diversification of the regulation of horizontally transferred genes such that their selective expression could potentially allow E. coli strains to adapt to new habitats.

Regulation of Bacterial DNA Packaging in Early Stationary Phase by Competitive DNA Binding of Dps and IHF

The bacterial nucleoid, a bacterial genome packed by nucleoid binding proteins, forms the physical basis for cellular processes such as gene transcription and DNA replication. Bacteria need to dynamically modulate their nucleoid structures at different growth phases and in response to environmental changes. At the nutrients deficient stationary phase, DNA-binding proteins from starved cells (Dps) and Integration host factors (IHF) are the two most abundant nucleoid associated proteins in E. coli. Yet, it remains unclear how the nucleoid architecture is controlled by the interplay between these two proteins, as well as the nucleoid's response to environmental changes. This question is addressed here using single DNA manipulation approach. Our results reveal that the two proteins are differentially selected for DNA binding, which can be tuned by changing environmental factors over physiological ranges including KCl (50-300 mM), MgCl2 (0-10 mM), pH (6.5-8.5) and temperature (23-37 °C). Increasing pH and MgCl2 concentrations switch from Dps-binding to IHF-binding. Stable Dps-DNA and IHF-DNA complexes are insensitive to temperature changes for the range tested. The environment dependent selection between IHF and Dps results in different physical organizations of DNA. Overall, our findings provide important insights into E. coli nucleoid architecture.

Roles of Nucleoid-Associated Proteins in Stress-Induced Mutagenic Break Repair in Starving Escherichia coli

The mutagenicity of DNA double-strand break repair in Escherichia coli is controlled by DNA-damage (SOS) and general (RpoS) stress responses, which let error-prone DNA polymerases participate, potentially accelerating evolution during stress. Either base substitutions and indels or genome rearrangements result. Here we discovered that most small basic proteins that compact the genome, nucleoid-associated proteins (NAPs), promote or inhibit mutagenic break repair (MBR) via different routes. Of 15 NAPs, H-NS, Fis, CspE, and CbpA were required for MBR; Dps inhibited MBR; StpA and Hha did neither; and five others were characterized previously. Three essential genes were not tested. Using multiple tests, we found the following: First, Dps, which reduces reactive oxygen species (ROS), inhibited MBR, implicating ROS in MBR. Second, CbpA promoted F' plasmid maintenance, allowing MBR to be measured in an F'-based assay. Third, Fis was required for activation of the SOS DNA-damage response and could be substituted in MBR by SOS-induced levels of DinB error-prone DNA polymerase. Thus, Fis promoted MBR by allowing SOS activation. Fourth, H-NS represses ROS detoxifier sodB and was substituted in MBR by deletion of sodB, which was not otherwise mutagenic. We conclude that normal ROS levels promote MBR and that H-NS promotes MBR by maintaining ROS. CspE positively regulates RpoS, which is required for MBR. Four of five previously characterized NAPs promoted stress responses that enhance MBR. Hence, most NAPs affect MBR, the majority via regulatory functions. The data show that a total of six NAPs promote MBR by regulating stress responses, indicating the importance of nucleoid structure and function to the regulation of MBR and of coupling mutagenesis to stress, creating genetic diversity responsively.

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.

Mutation bias, rather than binding preference, underlies the nucleosome-associated G+C% variation in eukaryotes

The effects of genetic content on epigenetic status have been extensively studied, but how epigenetic status affects genetic content is not well understood. As a key epigenetic factor the nucleosome structure is highly correlated with local G+C% in eukaryotic genomes. The prevailing explanation to the pattern is that nucleosome occupancy favors higher G+C% sequences more than lower G+C% sequences. However, recent observation of a biased mutation spectrum caused by nucleosome occupancy suggests that the higher G+C% of nucleosomal DNA might be the evolutionary consequence of nucleosome occupancy. To distinguish the two explanations, we examined data from an in vitro nucleosome reconstitution experiment in which histones are incubated with yeast Saccharomyces cerevisiae and Escherichia coli genomic DNA, the former has been shaped by nucleosome structure while the latter has not. There is a strong positive correlation between nucleosome density and G+C% for the yeast DNA, an observation consistent with in vivo data, and such a pattern nearly vanishes for E. coli genomic DNA, suggesting that biased mutation, rather than biased occupancy, explains the most nucleosome-associated G+C% variation in eukaryotic genomes.

Strong purifying selection at synonymous sites in D. melanogaster

Synonymous sites are generally assumed to be subject to weak selective constraint. For this reason, they are often neglected as a possible source of important functional variation. We use site frequency spectra from deep population sequencing data to show that, contrary to this expectation, 22% of four-fold synonymous (4D) sites in Drosophila melanogaster evolve under very strong selective constraint while few, if any, appear to be under weak constraint. Linking polymorphism with divergence data, we further find that the fraction of synonymous sites exposed to strong purifying selection is higher for those positions that show slower evolution on the Drosophila phylogeny. The function underlying the inferred strong constraint appears to be separate from splicing enhancers, nucleosome positioning, and the translational optimization generating canonical codon bias. The fraction of synonymous sites under strong constraint within a gene correlates well with gene expression, particularly in the mid-late embryo, pupae, and adult developmental stages. Genes enriched in strongly constrained synonymous sites tend to be particularly functionally important and are often involved in key developmental pathways. Given that the observed widespread constraint acting on synonymous sites is likely not limited to Drosophila, the role of synonymous sites in genetic disease and adaptation should be reevaluated.