An underlying mechanism for the increased mutagenesis of lagging-strand genes in Bacillus subtilis

Pathogenic bacteria experience pervasive RNA polymerase backtracking during infection

IMPORTANCE

Eukaryotic hosts have defense mechanisms that may disrupt molecular transactions along the pathogen's chromosome through excessive DNA damage. Given that DNA damage stalls RNA polymerase (RNAP) thereby increasing mutagenesis, investigating how host defense mechanisms impact the movement of the transcription machinery on the pathogen chromosome is crucial. Using a new methodology we developed, we elucidated the dynamics of RNAP movement and association with the chromosome in the pathogenic bacterium Salmonella enterica during infection. We found that dynamics of RNAP movement on the chromosome change significantly during infection genome-wide, including at regions that encode for key virulence genes. In particular, we found that there is pervasive RNAP backtracking on the bacterial chromosome during infections and that anti-backtracking factors are critical for pathogenesis. Altogether, our results suggest that, interestingly, the host environment can promote the development of antimicrobial resistance and hypervirulence as stalled RNAPs can accelerate evolution through increased mutagenesis.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

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.

Transcription-coupled global genomic repair in E. coli

The nucleotide excision repair (NER) pathway removes helix-distorting lesions from DNA in all organisms. Escherichia coli has long been a model for understanding NER, which is traditionally divided into major and minor subpathways known as global genome repair (GGR) and transcription-coupled repair (TCR), respectively. TCR has been assumed to be mediated exclusively by Mfd, a DNA translocase of minimal NER phenotype. This review summarizes the evidence that shaped the traditional view of NER in bacteria, and reviews data supporting a new model in which GGR and TCR are inseparable. In this new model, RNA polymerase serves both as the essential primary sensor of bulky DNA lesions genome-wide and as the delivery platform for the assembly of functional NER complexes in living cells.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome wide

RNA:DNA hybrids compromise replication fork progression and genome integrity in all cells. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate strongly in noncoding RNAs and 5'-UTRs of coding sequences. For ΔrnhB, hybrids accumulate preferentially in untranslated regions and early in coding sequences. We show that hybrid accumulation is particularly sensitive to gene expression in ΔrnhC cells. DNA replication in ΔrnhC cells is disrupted, leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts cause mutagenesis and shape genome organization.

Pathogenic bacteria experience pervasive RNA polymerase backtracking during infection

Pathogenic bacteria and their eukaryotic hosts are in a constant arms race. Hosts have numerous defense mechanisms at their disposal that not only challenge the bacterial invaders, but have the potential to disrupt molecular transactions along the bacterial chromosome. However, it is unclear how the host impacts association of proteins with the bacterial chromosome at the molecular level during infection. This is partially due to the lack of a method that could detect these events in pathogens while they are within host cells. We developed and optimized a system capable of mapping and measuring levels of bacterial proteins associated with the chromosome while they are actively infecting the host (referred to as PIC-seq). Here, we focused on the dynamics of RNA polymerase (RNAP) movement and association with the chromosome in the pathogenic bacterium Salmonella enterica as a model system during infection. Using PIC-seq, we found that RNAP association patterns with the chromosome change during infection genome-wide, including at regions that encode for key virulence genes. Importantly, we found that infection of a host significantly increases RNAP backtracking on the bacterial chromosome. RNAP backtracking is the most common form of disruption to RNAP progress on the chromosome. Interestingly, we found that the resolution of backtracked RNAPs via the anti-backtracking factors GreA and GreB is critical for pathogenesis, revealing a new class of virulence genes. Altogether, our results strongly suggest that infection of a host significantly impacts transcription by disrupting RNAP movement on the chromosome within the bacterial pathogen. The increased backtracking events have important implications not only for efficient transcription, but also for mutation rates as stalled RNAPs increase the levels of mutagenesis.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome-wide

RNA:DNA hybrids such as R-loops affect genome integrity and DNA replication fork progression. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome-wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation, but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate most strongly in non-coding RNAs and 5'-UTRs of coding sequences. For Δ rnhB , hybrids accumulate preferentially in untranslated regions and early in coding sequences. Hybrid accumulation is particularly sensitive to gene expression in Δ rnhC ; in cells lacking RnhC, DNA replication is disrupted leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts interact with replication to cause mutagenesis and shape genome organization.

The in vivo measurement of replication fork velocity and pausing by lag-time analysis

An important step towards understanding the mechanistic basis of the central dogma is the quantitative characterization of the dynamics of nucleic-acid-bound molecular motors in the context of the living cell. To capture these dynamics, we develop lag-time analysis, a method for measuring in vivo dynamics. Using this approach, we provide quantitative locus-specific measurements of fork velocity, in units of kilobases per second, as well as replisome pause durations, some with the precision of seconds. The measured fork velocity is observed to be both locus and time dependent, even in wild-type cells. In this work, we quantitatively characterize known phenomena, detect brief, locus-specific pauses at ribosomal DNA loci in wild-type cells, and observe temporal fork velocity oscillations in three highly-divergent bacterial species.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

The Impact of RNA-DNA Hybrids on Genome Integrity in Bacteria

During the essential processes of DNA replication and transcription, RNA-DNA hybrid intermediates are formed that pose significant risks to genome integrity when left unresolved. To manage RNA-DNA hybrids, all cells rely on RNase H family enzymes that specifically cleave the RNA portion of the many different types of hybrids that form in vivo. Recent experimental advances have provided new insight into how RNA-DNA hybrids form and the consequences to genome integrity that ensue when persistent hybrids remain unresolved. Here we review the types of RNA-DNA hybrids, including R-loops, RNA primers, and ribonucleotide misincorporations, that form during DNA replication and transcription and discuss how each type of hybrid can contribute to genome instability in bacteria. Further, we discuss how bacterial RNase HI, HII, and HIII and bacterial FEN enzymes contribute to genome maintenance through the resolution of hybrids.

Dynamic switching of crotonylation to ubiquitination of H2A at lysine 119 attenuates transcription-replication conflicts caused by replication stress

The reversible post-translational modification (PTM) of proteins plays an important role in many cellular processes. Lysine crotonylation (Kcr) is a newly identified PTM, but its functional significance remains unclear. Here, we found that Kcr is involved in the replication stress response. We show that crotonylation of histone H2A at lysine 119 (H2AK119) and ubiquitination of H2AK119 are reversibly regulated by replication stress. Decrotonylation of H2AK119 by SIRT1 is a prerequisite for subsequent ubiquitination of H2AK119 by BMI1. Accumulation of ubiquitinated H2AK119 at reversed replication forks leads to the release of RNA Polymerase II and transcription repression in the vicinity of stalled replication forks. These effects attenuate transcription–replication conflicts (TRCs) and TRC-associated R-loop formation and DNA double-strand breaks. These findings suggest that decrotonylation and ubiquitination of H2A at lysine 119 act together to resolve replication stress-induced TRCs and protect genome stability.

Stress-Induced Mutagenesis, Gambler Cells, and Stealth Targeting Antibiotic-Induced Evolution

Mechanisms of evolution and evolution of antibiotic resistance are both fundamental and world health problems. Stress-induced mutagenesis defines mechanisms of mutagenesis upregulated by stress responses, which drive adaptation when cells are maladapted to their environments—when stressed. Work in mutagenesis induced by antibiotics had produced tantalizing clues but not coherent mechanisms. We review recent advances in antibiotic-induced mutagenesis that integrate how reactive oxygen species (ROS), the SOS and general stress responses, and multichromosome cells orchestrate a stress response-induced switch from high-fidelity to mutagenic repair of DNA breaks. Moreover, while sibling cells stay stable, a mutable “gambler” cell subpopulation is induced by differentially generated ROS, which signal the general stress response. We discuss other evolvable subpopulations and consider diverse evolution-promoting molecules as potential targets for drugs to slow evolution of antibiotic resistance, cross-resistance, and immune evasion. An FDA-approved drug exemplifies “stealth” evolution-slowing drugs that avoid selecting resistance to themselves or antibiotics.

2-Dimensional Genetic Algorithm Exhibited an Essentiality of Gene Interaction for Evolution

Organisms consist of several genetic factors differing between species. However, the evolutionary effects of gene interactions on the evolutionary rate, adaptation, and divergence of organisms remain unknown. In a previous study, the 2-dimensional genetic algorithm (2DGA) program, including a gene interaction parameter, could simulate punctuated equilibrium under the disparity mode. Following this, we verified the effect of the number of gene interactions (gene cluster size) on evolution speed, adaptation, and divergence using the advanced 2DGA program. In this program, the population was replicated, mutated, and selected for 200,000 generations, and the fitness score, divergence, number of population, and genotype were output and plotted. The genotype data were used for evaluating the phylogenetic relations among the population. The gene cluster size 1) affected the disparity and parity mutagenesis modes differently, 2) determined the growth/exclusion rate and error threshold, and 3) accelerated or decelerated the population's speed of evolutionary advancement. In particular, when the gene cluster size expanded, the rate of increase in fitness scores decreased independently of the mutation rate and mode of mutation (disparity mode/parity mode). The mutation rate at the error threshold was also decreased by expanding the gene cluster size. Dendrograms traced the genotypes of the simulated population, indicating that the disparity mode caused the evolutionary process to enter 1) a stun mode, 2) an evolution mode, or 3) a divergence mode based on the mutation rate and gene cluster size, while the parity mode did not cause the population to enter a stun mode. Based on the above findings, we compared the predictions of the present study with evolution observed in the laboratory or the natural world and the processes of ongoing virus evolution, suggesting that our findings possibly explained the real evolution.

OUP accepted manuscript

Cellular DNA is continuously transcribed into RNA by multisubunit RNA polymerases (RNAPs). The continuity of transcription can be disrupted by DNA lesions that arise from the activities of cellular enzymes, reactions with endogenous and exogenous chemicals or irradiation. Here, we review available data on translesion RNA synthesis by multisubunit RNAPs from various domains of life, define common principles and variations in DNA damage sensing by RNAP, and consider existing controversies in the field of translesion transcription. Depending on the type of DNA lesion, it may be correctly bypassed by RNAP, or lead to transcriptional mutagenesis, or result in transcription stalling. Various lesions can affect the loading of the templating base into the active site of RNAP, or interfere with nucleotide binding and incorporation into RNA, or impair RNAP translocation. Stalled RNAP acts as a sensor of DNA damage during transcription-coupled repair. The outcome of DNA lesion recognition by RNAP depends on the interplay between multiple transcription and repair factors, which can stimulate RNAP bypass or increase RNAP stalling, and plays the central role in maintaining the DNA integrity. Unveiling the mechanisms of translesion transcription in various systems is thus instrumental for understanding molecular pathways underlying gene regulation and genome stability.

DNA Replication-Transcription Conflicts Do Not Significantly Contribute to Spontaneous Mutations Due to Replication Errors in Escherichia coli

Encounters between DNA replication and transcription can cause genomic disruption, particularly when the two meet head-on. Whether these conflicts produce point mutations is debated. This paper presents detailed analyses of a large collection of mutations generated during mutation accumulation experiments with mismatch repair (MMR)-defective Escherichia coli. With MMR absent, mutations are primarily due to DNA replication errors. Overall, there were no differences in the frequencies of base pair substitutions or small indels (i.e., insertion and deletions of ≤4 bp) in the coding sequences or promoters of genes oriented codirectionally versus head-on to replication. Among a subset of highly expressed genes, there was a 2- to 3-fold bias for indels in genes oriented head-on to replication, but this difference was almost entirely due to the asymmetrical genomic locations of tRNA genes containing mononucleotide runs, which are hot spots for indels. No additional orientation bias in mutation frequencies occurred when MMR⁻ strains were also defective for transcription-coupled repair (TCR). However, in contrast to other reports, loss of TCR slightly increased the overall mutation rate, meaning that TCR is antimutagenic. There was no orientation bias in mutation frequencies among the stress response genes that are regulated by RpoS or induced by DNA damage. Thus, biases in the locations of mutational targets can account for most, if not all, apparent biases in mutation frequencies between genes oriented head-on versus codirectional to replication. In addition, the data revealed a strong correlation of the frequency of base pair substitutions with gene length but no correlation with gene expression levels.

Mfd - at the crossroads of bacterial DNA repair, transcriptional regulation and molecular evolvability

For survival, bacteria need to continuously evolve and adapt to complex environments, including those that may impact the integrity of the DNA, the repository of genetic information to be passed on to future generations. The multiple factors of DNA repair share the substrate on which they operate with other key cellular machineries, principally those of replication and transcription, implying a high degree of coordination of DNA-based activities. In this review, I focus on progress made in the understanding of the protein factors operating at the crossroads of these three fundamental processes, with emphasis on the mutation frequency decline protein (Mfd, aka TRCF). Although Mfd research has a rich history that goes back in time for more than half a century, recent reports hint that much remains to be uncovered. I argue that besides being a transcription-repair coupling factor (TRCF), Mfd is also a global regulator of transcription and a pro-mutagenic factor, and that the way it interfaces with transcription, replication and nucleotide excision repair makes it an attractive candidate for the development of strategies to curb molecular evolution, hence, antibiotic resistance.

Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis

Transcription-induced mutagenic mechanisms limit genetic changes to times when expression happens and to coding DNA. It has been hypothesized that intrinsic sequences that have the potential to form alternate DNA structures, such as non-B DNA structures, influence these mechanisms. Non-B DNA structures are promoted by transcription and induce genome instability in eukaryotic cells, but their impact in bacterial genomes is less known. Here, we investigated if G4 DNA- and hairpin-forming motifs influence stationary-phase mutagenesis in Bacillus subtilis. We developed a system to measure the influence of non-B DNA on B. subtilis stationary-phase mutagenesis by deleting the wild-type argF at its chromosomal position and introducing IPTG-inducible argF alleles differing in their ability to form hairpin and G4 DNA structures into an ectopic locus. Using this system, we found that sequences predicted to form non-B DNA structures promoted mutagenesis in B. subtilis stationary-phase cells; such a response did not occur in growing conditions. We also found that the transcription-coupled repair factor Mfd promoted mutagenesis at these predicted structures. In summary, we showed that non-B DNA-forming motifs promote genetic instability, particularly in coding regions in stressed cells; therefore, non-B DNA structures may have a spatial and temporal mutagenic effect in bacteria. This study provides insights into mechanisms that prevent or promote mutagenesis and advances our understanding of processes underlying bacterial evolution.

Mfd regulates RNA polymerase association with hard-to-transcribe regions in vivo, especially those with structured RNAs

RNA polymerase (RNAP) encounters various roadblocks during transcription. These obstacles can impede RNAP movement and influence transcription, ultimately necessitating the activity of RNAP-associated factors. One such factor is the bacterial protein Mfd, a highly conserved DNA translocase and evolvability factor that interacts with RNAP. Although Mfd is thought to function primarily in the repair of DNA lesions that stall RNAP, increasing evidence suggests that it may also be important for transcription regulation. However, this is yet to be fully characterized. To shed light on Mfd's in vivo functions, we identified the chromosomal regions where it associates. We analyzed Mfd's impact on RNAP association and transcription regulation genome-wide. We found that Mfd represses RNAP association at many chromosomal regions. We found that these regions show increased RNAP pausing, suggesting that they are hard to transcribe. Interestingly, we noticed that the majority of the regions where Mfd regulates transcription contain highly structured regulatory RNAs. The RNAs identified regulate a myriad of biological processes, ranging from metabolism to transfer RNA regulation to toxin-antitoxin (TA) functions. We found that cells lacking Mfd are highly sensitive to toxin overexpression. Finally, we found that Mfd promotes mutagenesis in at least one toxin gene, suggesting that its function in regulating transcription may promote evolution of certain TA systems and other regions containing strong RNA secondary structures. We conclude that Mfd is an RNAP cofactor that is important, and at times critical, for transcription regulation at hard-to-transcribe regions, especially those that express structured regulatory RNAs.

Cotranscriptional R-loop formation by Mfd involves topological partitioning of DNA

R-loops are nucleic acid hybrids which form when an RNA invades duplex DNA to pair with its template sequence. Although they are implicated in a growing number of gene regulatory processes, their mechanistic origins remain unclear. We here report real-time observations of cotranscriptional R-loop formation at single-molecule resolution and propose a mechanism for their formation. We show that the bacterial Mfd protein can simultaneously interact with both elongating RNA polymerase and upstream DNA, tethering the two together and partitioning the DNA into distinct supercoiled domains. A highly negatively supercoiled domain forms in between Mfd and RNA polymerase, and compensatory positive supercoiling appears in front of the RNA polymerase and behind Mfd. The nascent RNA invades the negatively supercoiled domain and forms a stable R-loop that can drive mutagenesis. This mechanism theoretically enables any protein that simultaneously binds an actively translocating RNA polymerase and upstream DNA to stimulate R-loop formation.

Topological stress is responsible for the detrimental outcomes of head-on replication-transcription conflicts

Conflicts between the replication and transcription machineries have profound effects on chromosome duplication, genome organization, and evolution across species. Head-on conflicts (lagging-strand genes) are significantly more detrimental than codirectional conflicts (leading-strand genes). The fundamental reason for this difference is unknown. Here, we report that topological stress significantly contributes to this difference. We find that head-on, but not codirectional, conflict resolution requires the relaxation of positive supercoils by the type II topoisomerases DNA gyrase and Topo IV, at least in the Gram-positive model bacterium Bacillus subtilis. Interestingly, our data suggest that after positive supercoil resolution, gyrase introduces excessive negative supercoils at head-on conflict regions, driving pervasive R-loop formation. Altogether, our results reveal a fundamental mechanistic difference between the two types of encounters, addressing a long-standing question in the field of replication-transcription conflicts.

Lang and Merrikh show that resolution of head-on, but not codirectional, conflicts between replication and transcription machineries requires type II topoisomerases, suggesting that a fundamental difference between the two types of conflicts is supercoil buildup in DNA. Furthermore, they show that supercoil resolution at head-on conflict regions drives R-loop formation.

Asymmetron: a toolkit for the identification of strand asymmetry patterns in biological sequences

DNA strand asymmetries can have a major effect on several biological functions, including replication, transcription and transcription factor binding. As such, DNA strand asymmetries and mutational strand bias can provide information about biological function. However, a versatile tool to explore this does not exist. Here, we present Asymmetron, a user-friendly computational tool that performs statistical analysis and visualizations for the evaluation of strand asymmetries. Asymmetron takes as input DNA features provided with strand annotation and outputs strand asymmetries for consecutive occurrences of a single DNA feature or between pairs of features. We illustrate the use of Asymmetron by identifying transcriptional and replicative strand asymmetries of germline structural variant breakpoints. We also show that the orientation of the binding sites of 45% of human transcription factors analyzed have a significant DNA strand bias in transcribed regions, that is also corroborated in ChIP-seq analyses, and is likely associated with transcription. In summary, we provide a novel tool to assess DNA strand asymmetries and show how it can be used to derive new insights across a variety of biological disciplines.

Visualizing mutagenic repair: novel insights into bacterial translesion synthesis

DNA repair is essential for cell survival. In all domains of life, error-prone and error-free repair pathways ensure maintenance of genome integrity under stress. Mutagenic, low-fidelity repair mechanisms help avoid potential lethality associated with unrepaired damage, thus making them important for genome maintenance and, in some cases, the preferred mode of repair. However, cells carefully regulate pathway choice to restrict activity of these pathways to only certain conditions. One such repair mechanism is translesion synthesis (TLS), where a low-fidelity DNA polymerase is employed to synthesize across a lesion. In bacteria, TLS is a potent source of stress-induced mutagenesis, with potential implications in cellular adaptation as well as antibiotic resistance. Extensive genetic and biochemical studies, predominantly in Escherichia coli, have established a central role for TLS in bypassing bulky DNA lesions associated with ongoing replication, either at or behind the replication fork. More recently, imaging-based approaches have been applied to understand the molecular mechanisms of TLS and how its function is regulated. Together, these studies have highlighted replication-independent roles for TLS as well. In this review, we discuss the current status of research on bacterial TLS, with emphasis on recent insights gained mostly through microscopy at the single-cell and single-molecule level.

The roles of replication-transcription conflict in mutagenesis and evolution of genome organization

Replication-transcription conflicts promote mutagenesis and give rise to evolutionary signatures, with fundamental importance to genome stability ranging from bacteria to metastatic cancer cells. This review focuses on the interplay between replication-transcription conflicts and the evolution of gene directionality. In most bacteria, the majority of genes are encoded on the leading strand of replication such that their transcription is co-directional with the direction of DNA replication fork movement. This gene strand bias arises primarily due to negative selection against deleterious consequences of head-on replication-transcription conflict. However, many genes remain head-on. Can head-on orientation provide some benefit? We combine insights from both mechanistic and evolutionary studies, review published work, and analyze gene expression data to evaluate an emerging model that head-on genes are temporal targets for adaptive mutagenesis during stress. We highlight the alternative explanation that genes in the head-on orientation may simply be the result of genomic inversions and relaxed selection acting on nonessential genes. We seek to clarify how the mechanisms of replication-transcription conflict, in concert with other mutagenic mechanisms, balanced by natural selection, have shaped bacterial genome evolution.

Molecular determinants for dsDNA translocation by the transcription-repair coupling and evolvability factor Mfd

Mfd couples transcription to nucleotide excision repair, and acts on RNA polymerases when elongation is impeded. Depending on impediment severity, this action results in either transcription termination or elongation rescue, which rely on ATP-dependent Mfd translocation on DNA. Due to its role in antibiotic resistance, Mfd is also emerging as a prime target for developing anti-evolution drugs. Here we report the structure of DNA-bound Mfd, which reveals large DNA-induced structural changes that are linked to the active site via ATPase motif VI. These changes relieve autoinhibitory contacts between the N- and C-termini and unmask UvrA recognition determinants. We also demonstrate that translocation relies on a threonine in motif Ic, widely conserved in translocases, and a family-specific histidine near motif IVa, reminiscent of the “arginine clamp” of RNA helicases. Thus, Mfd employs a mode of DNA recognition that at its core is common to ss/ds translocases that act on DNA or RNA.

Transcription-repair coupling factors (TRCFs) are large ATPases that mediate the preferential repair of the transcribed DNA strand. Here the authors reveal the cryo-EM structure of DNA-bound Mfd, the bacterial TRCF, and provide molecular insights into its mode of action.

Targeting evolution to inhibit antibiotic resistance

Drug-resistant bacterial infections have led to a global health crisis. Although much effort is placed on the development of new antibiotics or variants that are less subject to existing resistance mechanisms, history shows that this strategy by itself is unlikely to solve the problem of drug resistance. Here, we discuss inhibiting evolution as a strategy that, in combination with antibiotics, may resolve the problem. Although mutagenesis is the main driver of drug resistance development, attacking the drivers of genetic diversification in pathogens has not been well explored. Bacteria possess active mechanisms that increase the rate of mutagenesis, especially at times of stress, such as during replication within eukaryotic host cells, or exposure to antibiotics. We highlight how the existence of these promutagenic proteins (evolvability factors) presents an opportunity that can be capitalized upon for the effective inhibition of drug resistance development. To help move this idea from concept to execution, we first describe a set of criteria that an 'optimal' evolvability factor would likely have to meet to be a viable therapeutic target. We then discuss the intricacies of some of the known mutagenic mechanisms and evaluate their potential as drug targets to inhibit evolution. In principle, and as suggested by recent studies, we argue that the inhibition of these and other evolvability factors should reduce resistance development. Finally, we discuss the challenges of transitioning anti-evolution drugs from the laboratory to the clinic.

Role of Mfd and GreA in Bacillus subtilis Base Excision Repair-Dependent Stationary-Phase Mutagenesis

We report that the absence of an oxidized guanine (GO) system or the apurinic/apyrimidinic (AP) endonucleases Nfo, ExoA, and Nth promoted stress-associated mutagenesis (SAM) in Bacillus subtilis YB955 (hisC952 metB5 leuC427). Moreover, MutY-promoted SAM was Mfd dependent, suggesting that transcriptional transactions over nonbulky DNA lesions promoted error-prone repair. Here, we inquired whether Mfd and GreA, which control transcription-coupled repair and transcription fidelity, influence the mutagenic events occurring in nutritionally stressed B. subtilis YB955 cells deficient in the GO or AP endonuclease repair proteins. To this end, mfd and greA were disabled in genetic backgrounds defective in the GO and AP endonuclease repair proteins, and the strains were tested for growth-associated and stress-associated mutagenesis. The results revealed that disruption of mfd or greA abrogated the production of stress-associated amino acid revertants in the GO and nfo exoA nth strains, respectively. These results suggest that in nutritionally stressed B. subtilis cells, spontaneous nonbulky DNA lesions are processed in an error-prone manner with the participation of Mfd and GreA. In support of this notion, stationary-phase ΔytkD ΔmutM ΔmutY (referred to here as ΔGO) and Δnfo ΔexoA Δnth (referred to here as ΔAP) cells accumulated 8-oxoguanine (8-OxoG) lesions, which increased significantly following Mfd disruption. In contrast, during exponential growth, disruption of mfd or greA increased the production of His⁺, Met⁺, or Leu⁺ prototrophs in both DNA repair-deficient strains. Thus, in addition to unveiling a role for GreA in mutagenesis, our results suggest that Mfd and GreA promote or prevent mutagenic events driven by spontaneous genetic lesions during the life cycle of B. subtilis IMPORTANCE In this paper, we report that spontaneous genetic lesions of an oxidative nature in growing and nutritionally stressed B. subtilis strain YB955 (hisC952 metB5 leuC427) cells drive Mfd- and GreA-dependent repair transactions. However, whereas Mfd and GreA elicit faithful repair events during growth to maintain genome fidelity, under starving conditions, both factors promote error-prone repair to produce genetic diversity, allowing B. subtilis to escape from growth-limiting conditions.

Tumor treating fields cause replication stress and interfere with DNA replication fork maintenance: Implications for cancer therapy

Tumor treating fields (TTFields) is a noninvasive physical modality of cancer therapy that applies low-intensity, intermediate frequency, and alternating electric fields to a tumor. Interference with mitosis was the first mechanism describing the effects of TTFields on cancer cells; however, TTFields was shown to not only reduce the rejoining of radiation-induced DNA double-strand breaks (DSBs), but to also induce DNA DSBs. The mechanism(s) by which TTFields generates DNA DSBs is related to the generation of replication stress including reduced expression of the DNA replication complex genes MCM6 and MCM10 and the Fanconi's Anemia pathway genes. When markers of DNA replication stress as a result of TTFields exposure were examined, newly replicated DNA length was reduced with TTFields exposure time and there was increased R-loop formation. Furthermore, as cells were exposed to TTFields a conditional vulnerability environment developed which rendered cells more susceptible to DNA damaging agents or agents that interfere with DNA repair or replication fork maintenance. The effect of TTFields exposure with concomitant exposure to cisplatin or PARP inhibition, the combination of TTFields plus concomitant PARP inhibition followed by radiation, or radiation alone at the end of a TTFields exposure were all synergistic. Finally, gene expression analysis of 47 key mitosis regulator genes suggested that TTFields-induced mitotic aberrations and DNA damage/replication stress events, although intimately linked to one another, are likely initiated independently of one another. This suggests that enhanced replication stress and reduced DNA repair capacity are also major mechanisms of TTFields effects, effects for which there are therapeutic implications.

H3K4 methylation at active genes mitigates transcription-replication conflicts during replication stress

Transcription-replication conflicts (TRCs) occur when intensive transcriptional activity compromises replication fork stability, potentially leading to gene mutations. Transcription-deposited H3K4 methylation (H3K4me) is associated with regions that are susceptible to TRCs; however, the interplay between H3K4me and TRCs is unknown. Here we show that H3K4me aggravates TRC-induced replication failure in checkpoint-defective cells, and the presence of methylated H3K4 slows down ongoing replication. Both S-phase checkpoint activity and H3K4me are crucial for faithful DNA synthesis under replication stress, especially in highly transcribed regions where the presence of H3K4me is highest and TRCs most often occur. H3K4me mitigates TRCs by decelerating ongoing replication, analogous to how speed bumps slow down cars. These findings establish the concept that H3K4me defines the transcriptional status of a genomic region and defends the genome from TRC-mediated replication stress and instability.

Transcription-replication conflicts (TRC) can contribute to genome instability. Here the authors reveal that under replication stress H3K4 methylation can play a role in TRC prevention.

Role of RNase H enzymes in maintaining genome stability in Escherichia coli expressing a steric-gate mutant of pol VICE391

pol V_ICE391 (RumA'₂B) is a low-fidelity polymerase that promotes considerably higher levels of spontaneous "SOS-induced" mutagenesis than the related E. coli pol V (UmuD'₂C). The molecular basis for the enhanced mutagenesis was previously unknown. Using single molecule fluorescence microscopy to visualize pol V enzymes, we discovered that the elevated levels of mutagenesis are likely due, in part, to prolonged binding of RumB to genomic DNA leading to increased levels of DNA synthesis compared to UmuC. We have generated a steric gate pol V_ICE391 variant (pol V_ICE391_Y13A) that readily misincorporates ribonucleotides into the E. coli genome and have used the enzyme to investigate the molecular mechanisms of Ribonucleotide Excision Repair (RER) under conditions of increased ribonucleotide-induced stress. To do so, we compared the extent of spontaneous mutagenesis promoted by pol V and pol V_ICE391 to that of their respective steric gate variants. Levels of mutagenesis promoted by the steric gate variants that are lower than that of the wild-type enzyme are indicative of active RER that removes misincorporated ribonucleotides, but also misincorporated deoxyribonucleotides from the genome. Using such an approach, we confirmed that RNase HII plays a pivotal role in RER. In the absence of RNase HII, Nucleotide Excision Repair (NER) proteins help remove misincorporated ribonucleotides. However, significant RER occurs in the absence of RNase HII and NER. Most of the RNase HII and NER-independent RER occurs on the lagging strand during genome duplication. We suggest that this is most likely due to efficient RNase HI-dependent RER which recognizes the polyribonucleotide tracts generated by pol V_ICE391_Y13A. These activities are critical for the maintenance of genomic integrity when RNase HII is overwhelmed, or inactivated, as ΔrnhB or ΔrnhB ΔuvrA strains expressing pol V_ICE391_Y13A exhibit genome and plasmid instability in the absence of RNase HI.

The enigmatic role of Mfd in replication-transcription conflicts in bacteria

Conflicts between replication and transcription can have life-threatening consequences. RNA polymerase (RNAP) is the major impediment to replication progression, and its efficient removal from DNA should mitigate the consequences of collisions with replication. Cells have various proteins that can resolve conflicts by removing stalled (or actively translocating) RNAP from DNA. It would therefore seem logical that RNAP-associated factors, such as the bacterial DNA translocase Mfd, would minimize the effects of conflicts. Despite seemingly conclusive statements in most textbooks, the role of Mfd in conflicts remains an enigma. In this review, we will discuss the different physical states of RNAP during transcription, and how each distinct state can influence conflict severity and potentially trigger the involvement of Mfd. We propose models to explain the contradictory conclusions from published studies on the potential role of Mfd in resolving conflicts.

Transcription-Coupled Repair: From Cells to Single Molecules and Back Again

Transcription-coupled repair is mediated by the Mfd protein. TCR is defined as the preferential repair of DNA lesions in the transcribed strand of actively transcribed genes, and is opposed to the strand-aspecific global genome repair. The Mfd protein mediates TCR by binding to and displacing RNA polymerase, which is stalled at a DNA lesion on the transcribed strand of DNA, then recruiting UvrA and UvrB. The repair cascade results in the recruitment of, and DNA excision by, UvrC; removal of the damage-bearing oligonucleotide by UvrD; "filling-in" of the DNA by DNA polymerase; and sealing of the strands by DNA ligase. The gene required for Mfd was originally identified as a gene needed for the "mutation frequency decline" phenotype in which the repair of certain UV-induced lesions in the transcribed strand of tRNA genes is increased when cells are forced to delay replication immediately following UV exposure. This review will focus on the genetics that led to the discovery of the Mfd gene; summarize the subsequent biochemical, structural and single-molecule interrogations of the Mfd protein; and explore the more recent findings of Mfd in mutagenesis.

New insights into the structures and interactions of bacterial Y-family DNA polymerases

Bacterial Y-family DNA polymerases are usually classified into DinB (Pol IV), UmuC (the catalytic subunit of Pol V) and ImuB, a catalytically dead essential component of the ImuA-ImuB-DnaE2 mutasome. However, the true diversity of Y-family polymerases is unknown. Furthermore, for most of them the structures are unavailable and interactions are poorly characterized. To gain a better understanding of bacterial Y-family DNA polymerases, we performed a detailed computational study. It revealed substantial diversity, far exceeding traditional classification. We found that a large number of subfamilies feature a C-terminal extension next to the common Y-family region. Unexpectedly, in most C-terminal extensions we identified a region homologous to the N-terminal oligomerization motif of RecA. This finding implies a universal mode of interaction between Y-family members and RecA (or ImuA), in the case of Pol V strongly supported by experimental data. In gram-positive bacteria, we identified a putative Pol V counterpart composed of a Y-family polymerase, a YolD homolog and RecA. We also found ImuA-ImuB-DnaE2 variants lacking ImuA, but retaining active or inactive Y-family polymerase, a standalone ImuB C-terminal domain and/or DnaE2. In summary, our analyses revealed that, despite considerable diversity, bacterial Y-family polymerases share previously unanticipated similarities in their structural domains/motifs and interactions.

Mfd protects against oxidative stress in Bacillus subtilis independently of its canonical function in DNA repair

BACKGROUND

Previous reports showed that mutagenesis in nutrient-limiting conditions is dependent on Mfd in Bacillus subtilis. Mfd initiates one type of transcription-coupled repair (TCR); this type of repair is known to target bulky lesions, like those associated with UV exposure. Interestingly, the roles of Mfd in repair of oxidative-promoted DNA damage and regulation of transcription differ. Here, we used a genetic approach to test whether Mfd protected B. subtilis from exposure to two different oxidants.

RESULTS

Wild-type cells survived tert-butyl hydroperoxide (t-BHP) exposure significantly better than Mfd-deficient cells. This protective effect was independent of UvrA, a component of the canonical TCR/nucleotide excision repair (NER) pathway. Further, our results suggest that Mfd and MutY, a DNA glycosylase that processes 8-oxoG DNA mismatches, work together to protect cells from lesions generated by oxidative damage. We also tested the role of Mfd in mutagenesis in starved cells exposed to t-BHP. In conditions of oxidative stress, Mfd and MutY may work together in the formation of mutations. Unexpectedly, Mfd increased survival when cells were exposed to the protein oxidant diamide. Under this type of oxidative stress, cells survival was not affected by MutY or UvrA.

CONCLUSIONS

These results are significant because they show that Mfd mediates error-prone repair of DNA and protects cells against oxidation of proteins by affecting gene expression; Mfd deficiency resulted in increased gene expression of the OhrR repressor which controls the cellular response to organic peroxide exposure. These observations point to Mfd functioning beyond a DNA repair factor in cells experiencing oxidative stress.

Inhibiting the Evolution of Antibiotic Resistance

Efforts to battle antimicrobial resistance (AMR) are generally focused on developing novel antibiotics. However, history shows that resistance arises regardless of the nature or potency of new drugs. Here, we propose and provide evidence for an alternate strategy to resolve this problem: inhibiting evolution. We determined that the DNA translocase Mfd is an “evolvability factor” that promotes mutagenesis and is required for rapid resistance development to all antibiotics tested across highly divergent bacterial species. Importantly, hypermutator alleles that accelerate AMR development did not arise without Mfd, at least during evolution of trimethoprim resistance. We also show that Mfd’s role in AMR development depends on its interactions with the RNA polymerase subunit RpoB and the nucleotide excision repair protein UvrA. Our findings suggest that AMR development can be inhibited through inactivation of evolvability factors (potentially with “anti-evolution” drugs)—in particular, Mfd—providing an unexplored route toward battling the AMR crisis.

•

The bacterial transcription-coupled repair (TCR) factor Mfd promotes mutagenesis

•

Mfd-driven mutagenesis accelerates the evolution of antimicrobial resistance (AMR)

•

The rapid evolution of AMR requires Mfd’s interaction with RpoB and UvrA

•

Mfd may be an ideal target for “anti-evolution” drugs that inhibit AMR development

Gene inversion potentiates bacterial evolvability and virulence

Most bacterial genes are encoded on the leading strand, co-orienting the movement of the replication machinery with RNA polymerases. This bias reduces the frequency of detrimental head-on collisions between the two machineries. The negative outcomes of these collisions should lead to selection against head-on alleles, maximizing genome co-orientation. Our findings challenge this model. Using the GC skew calculation, we reveal the evolutionary inversion record of all chromosomally encoded genes in multiple divergent bacterial pathogens. Against expectations, we find that a large number of co-oriented genes have inverted to the head-on orientation, presumably increasing the frequency of head-on replication-transcription conflicts. Furthermore, we find that head-on genes, (including key antibiotic resistance and virulence genes) have higher rates of non-synonymous mutations and are more frequently under positive selection (dN/dS > 1). Based on these results, we propose that spontaneous gene inversions can increase the evolvability and pathogenic capacity of bacteria through head-on replication-transcription collisions.

Head-on replication-transcription collisions occur within genes encoded on the lagging DNA strand. Here, the authors show that a large number of originally co-oriented (leading strand) genes have inverted to the head-on orientation, increasing both gene-specific mutation rates, and the overall evolvability of several bacterial pathogens.

Transcriptional coupling (Mfd) and DNA damage scanning (DisA) coordinate excision repair events for efficient Bacillus subtilis spore outgrowth

The absence of base excision repair (BER) proteins involved in processing ROS-promoted genetic insults activates a DNA damage scanning (DisA)-dependent checkpoint event in outgrowing Bacillus subtilis spores. Here, we report that genetic disabling of transcription-coupled repair (TCR) or nucleotide excision repair (NER) pathways severely affected outgrowth of ΔdisA spores, and much more so than the effects of these mutations on log phase growth. This defect delayed the first division of spore's nucleoid suggesting that unrepaired lesions affected transcription and/or replication during outgrowth. Accordingly, return to life of spores deficient in DisA/Mfd or DisA/UvrA was severely affected by a ROS-inducer or a replication blocking agent, hydrogen peroxide and 4-nitroquinoline-oxide, respectively. Mutation frequencies to rifampin resistance (Rifr ) revealed that DisA allowed faithful NER-dependent DNA repair but activated error-prone repair in TCR-deficient outgrowing spores. Sequencing analysis of rpoB from spontaneous Rifr colonies revealed that mutations resulting from base deamination predominated in outgrowing wild-type spores. Interestingly, a wide range of base substitutions promoted by oxidized DNA bases were detected in ΔdisA and Δmfd outgrown spores. Overall, our results suggest that Mfd and DisA coordinate excision repair events in spore outgrowth to eliminate DNA lesions that interfere with replication and transcription during this developmental period.

Involvement of transcription-coupled repair factor Mfd and DNA helicase UvrD in mutational processes in Pseudomonas putida

Stalled RNA polymerases (RNAPs) pose an obstacle for the replicating complexes, which could lead to transcription-replication conflicts and result in genetic instability. Stalled RNAPs and DNA lesions blocking RNAP elongation are removed by transcription-coupled repair (TCR), the process which in bacteria is mediated by TCR factor Mfd and helicase UvrD. Although the mechanism of TCR has been extensively studied, its role in mutagenesis is still obscure. In the current study we have investigated the role of Mfd and UvrD in mutational processes in soil bacterium Pseudomonas putida. Our results revealed that UvrD helicase is essential to prevent the emergence of mutations, as the loss of uvrD resulted in elevated mutant frequency both in exponential- and stationary-phase bacterial cultures. UvrD was also found to be necessary to survive DNA damage, but NER or MMR pathways are not completely abolished in UvrD-deficient P. putida. Mfd-deficiency had a moderate impact on surviving DNA damage and did not influence the frequency of mutations occurred in exponentially growing bacteria. However, the absence of Mfd caused approximately a two-fold decline in stationary-phase mutant frequency compared to the P. putida wild-type strain and suppressed the elevated mutant frequency observed in the ΔuvrD strain. Remarkably, the Mfd-deficient strain also formed less UV-induced mutants. These results suggest that in P. putida the Mfd-mediated TCR could be associated with UV- and stationary-phase mutagenesis.

Non-canonical processing of DNA photodimers with Bacillus subtilis UV-endonuclease YwjD, 5'→3' exonuclease YpcP and low-fidelity DNA polymerases YqjH and YqjW

It has been shown that mutation frequency decline (Mfd) and nucleotide excision repair (NER) deficiencies promote UV-induced mutagenesis in B. subtilis sporangia. As replication is halted in sporulating B. subtilis cells, in this report, we investigated if this response may result from an error-prone repair event involving the UV-endonuclease YwjD and low fidelity (LF) DNA synthesis. Accordingly, disruption of YwjD generated B. subtilis sporangia that were more susceptible to UV-C radiation than sporangia of the WT strain and such susceptibility increased even more after the single or simultaneous inactivation of the LF DNA polymerases YqjH and YqjW. To further explore this concept, functional His₆-tagged YwjD and Y-DNA polymerases YqjH and YqjW were produced and purified to homogeneity. In vitro repair assays showed that YwjD hydrolyzed the phosphodiester bond immediately located 5´-end of a ds-DNA substrate bearing either, cyclobutane pyrimidine dimers (CPD), 6-4 photoproducts (6-4 PD) or Dewar isomers (DWI). Notably, the 6-4 PD and DWI but not the CPDs repair intermediaries of YwjD were efficiently processed by the LF polymerase YqjH suggesting that an additional 5'→3' exonuclease event was necessary to process PD. Accordingly, the LF polymerase YqjW efficiently processed the incision-excision repair products derived from YwjD and exonuclease YpcP attack over CPD-containing DNA. In summary, our results unveiled a novel non-canonical repair pathway that employs YwjD to incise PD-containing DNA and low fidelity synthesis contributing thus to mutagenesis, survival and spore morphogenesis in B. subtilis.

The Clash of Macromolecular Titans: Replication-Transcription Conflicts in Bacteria

Within the last decade, it has become clear that DNA replication and transcription are routinely in conflict with each other in growing cells. Much of the seminal work on this topic has been carried out in bacteria, specifically, Escherichia coli and Bacillus subtilis; therefore, studies of conflicts in these species deserve special attention. Collectively, the recent findings on conflicts have fundamentally changed the way we think about DNA replication in vivo. Furthermore, new insights on this topic have revealed that the conflicts between replication and transcription significantly influence many key parameters of cellular function, including genome organization, mutagenesis, and evolution of stress response and virulence genes. In this review, we discuss the consequences of replication-transcription conflicts on the life of bacteria and describe some key strategies cells use to resolve them. We put special emphasis on two critical aspects of these encounters: ( a) the consequences of conflicts on replisome stability and dynamics, and ( b) the resulting increase in spontaneous mutagenesis.

Sources of spontaneous mutagenesis in bacteria

Mutations in an organism's genome can arise spontaneously, that is, in the absence of exogenous stress and prior to selection. Mutations are often neutral or deleterious to individual fitness but can also provide genetic diversity driving evolution. Mutagenesis in bacteria contributes to the already serious and growing problem of antibiotic resistance. However, the negative impacts of spontaneous mutagenesis on human health are not limited to bacterial antibiotic resistance. Spontaneous mutations also underlie tumorigenesis and evolution of drug resistance. To better understand the causes of genetic change and how they may be manipulated in order to curb antibiotic resistance or the development of cancer, we must acquire a mechanistic understanding of the major sources of mutagenesis. Bacterial systems are particularly well-suited to studying mutagenesis because of their fast growth rate and the panoply of available experimental tools, but efforts to understand mutagenic mechanisms can be complicated by the experimental system employed. Here, we review our current understanding of mutagenic mechanisms in bacteria and describe the methods used to study mutagenesis in bacterial systems.

Bacillus subtilis DisA helps to circumvent replicative stress during spore revival

The mechanisms that allow to circumvent replicative stress, and to resume DNA synthesis are poorly understood in Bacillus subtilis. To study the role of the diadenylate cyclase DisA and branch migration translocase (BMT) RadA/Sms in restarting a stalled replication fork, we nicked and broke the circular chromosome of an inert mature haploid spore, damaged the bases, and measured survival of reviving spores. During undisturbed ripening, nicks and breaks should be repaired by pathways that do not invoke long-range end resection or genetic exchange by homologous recombination, after which DNA replication might be initiated. We found that DNA damage reduced the viability of spores that lacked DisA, BMT (RadA/Sms, RuvAB or RecG), the Holliday junction resolvase RecU, or the translesion synthesis DNA polymerases (PolY1 or PolY2). DisA and RadA/Sms, in concert with RuvAB, RecG, RecU, PolY1 or PolY2, are needed to bypass replication-blocking lesions. DisA, which binds to stalled or reversed forks, did not apparently affect initiation of PriA-dependent DNA replication in vitro. We propose that DisA is necessary to coordinate responses to replicative stress; it could help to circumvent damaged template bases that otherwise impede fork progression.

Selection for energy efficiency drives strand-biased gene distribution in prokaryotes

Lagging-strand genes accumulate more deleterious mutations. Genes are thus preferably located on the leading strand, an observation known as strand-biased gene distribution (SGD). Despite of this mechanistic understanding, a satisfactory quantitative model is still lacking. Replication-transcription-collisions induce stalling of the replication machinery, expose DNA to various attacks, and are followed by error-prone repairs. We found that mutational biases in non-transcribed regions can explain ~71% of the variations in SGDs in 1,552 genomes, supporting the mutagenesis origin of SGD. Mutational biases introduce energetically cheaper nucleotides on the lagging strand, and result in more expensive protein products; consistently, the cost difference between the two strands explains ~50% of the variance in SGDs. Protein costs decrease with increasing gene expression. At similar expression levels, protein products of leading-strand genes are generally cheaper than lagging-strand genes; however, highly-expressed lagging genes are still cheaper than lowly-expressed leading genes. Selection for energy efficiency thus drives some genes to the leading strand, especially those highly expressed and essential, but certainly not all genes. Stronger mutational biases are often associated with low-GC genomes; as low-GC genes encode expensive proteins, low-GC genomes thus tend to have stronger SGDs to alleviate the stronger pressure on efficient energy usage.

Replication-Transcription Conflicts Generate R-Loops that Orchestrate Bacterial Stress Survival and Pathogenesis

Replication-transcription collisions shape genomes, influence evolution, and promote genetic diseases. Although unclear why, head-on transcription (lagging strand genes) is especially disruptive to replication and promotes genomic instability. Here, we find that head-on collisions promote R-loop formation in Bacillus subtilis. We show that pervasive R-loop formation at head-on collision regions completely blocks replication, elevates mutagenesis, and inhibits gene expression. Accordingly, the activity of the R-loop processing enzyme RNase HIII at collision regions is crucial for stress survival in B. subtilis, as many stress response genes are head-on to replication. Remarkably, without RNase HIII, the ability of the intracellular pathogen Listeria monocytogenes to infect and replicate in hosts is weakened significantly, most likely because many virulence genes are head-on to replication. We conclude that the detrimental effects of head-on collisions stem primarily from excessive R-loop formation and that the resolution of these structures is critical for bacterial stress survival and pathogenesis.

Implementation of a loss-of-function system to determine growth and stress-associated mutagenesis in Bacillus subtilis

A forward mutagenesis system based on the acquisition of mutations that inactivate the thymidylate synthase gene (TMS) and confer a trimethoprim resistant (Tmpr) phenotype was developed and utilized to study transcription-mediated mutagenesis (TMM). In addition to thyA, Bacillus subtilis possesses thyB, whose expression occurs under conditions of cell stress; therefore, we generated a thyB- thyA+ mutant strain. Tmpr colonies of this strain were produced with a spontaneous mutation frequency of ~1.4 × 10-9. Genetic disruption of the canonical mismatch (MMR) and guanine oxidized (GO) repair pathways increased the Tmpr frequency of mutation by ~2-3 orders of magnitude. A wide spectrum of base substitutions as well as insertion and deletions in the ORF of thyA were found to confer a Tmpr phenotype. Stationary-phase-associated mutagenesis (SPM) assays revealed that colonies with a Tmpr phenotype, accumulated over a period of ten days with a frequency of ~ 60 ×10-7. The Tmpr system was further modified to study TMM by constructing a ΔthyA ΔthyB strain carrying an IPTG-inducible Pspac-thyA cassette. In conditions of transcriptional induction of thyA, the generation of Tmpr colonies increased ~3-fold compared to conditions of transcriptional repression. Further, the Mfd and GreA factors were necessary for the generation of Tmpr colonies in the presence of IPTG in B. subtilis. Because GreA and Mfd facilitate transcription-coupled repair, our results suggest that TMM is a mechanim to produce genetic diversity in highly transcribed regions in growth-limited B. subtilis cells.

The Accelerated Evolution of Lagging Strand Genes Is Independent of Sequence Context

We previously discovered that lagging strand genes evolve faster in Bacillus subtilis (and potentially other bacteria). Lagging strand genes are transcribed in the head-on orientation with respect to DNA replication, leading to collisions between the two machineries that stall replication and can destabilize genomes. Our previous work indicated that the increased mutagenesis of head-on genes depends on transcription-coupled repair and the activity of an error prone polymerase which is likely activated in response to these collisions. Recently, it was proposed that sequence context is a major contributor to the increased mutagenesis and evolution of head-on genes. These models are based on laboratory-based evolution experiments performed in B. subtilis. However, critical evolutionary analyses of naturally occurring single nucleotide polymorphisms (SNPs) in wild strains were not performed. Using the genomic sequences from nine closely related wild B. subtilis strains, we analyzed over 200,000 naturally occurring SNPs as a proxy for natural mutation patterns for all genes and in particular, head-on genes. Our analysis suggests that (frame-independent) triplet sequence context can impact mutation rates: certain triplet sequences (TAG, CCC, CTA, and ACC) accumulate SNPs at a higher rate and are depleted from the genome. However, the triplet sequences previously identified as mutagenic in laboratory experiments (CCG, GCG, and CAC) do not have an elevated rate of SNP accumulation and are not depleted from the genome. Importantly, dN/dS analyses indicate that the accelerated evolution of head-on genes is not dependent on any particular triplet sequence. Thus, in agreement with our previous results, mutagenic transcription-coupled repair, rather than sequence context, is sufficient to explain the accelerated evolution of head-on genes.

Spatial and Temporal Control of Evolution through Replication-Transcription Conflicts

Evolution could potentially be accelerated if an organism could selectively increase the mutation rate of specific genes that are actively under positive selection. Recently, a mechanism that cells can use to target rapid evolution to specific genes was discovered. This mechanism is driven by gene orientation-dependent encounters between DNA replication and transcription machineries. These encounters increase mutagenesis in lagging-strand genes, where replication-transcription conflicts are severe. Due to the orientation and transcription-dependent nature of this process, conflict-driven mutagenesis can be used by cells to spatially (gene-specifically) and temporally (only upon transcription induction) regulate the rate of gene evolution. Here, I summarize recent findings on this topic, and discuss the implications of increasing mutagenesis rates and accelerating evolution through active mechanisms.

Transcription leads to pervasive replisome instability in bacteria

The canonical model of DNA replication describes a highly-processive and largely continuous process by which the genome is duplicated. This continuous model is based upon in vitro reconstitution and in vivo ensemble experiments. Here, we characterize the replisome-complex stoichiometry and dynamics with single-molecule resolution in bacterial cells. Strikingly, the stoichiometries of the replicative helicase, DNA polymerase, and clamp loader complexes are consistent with the presence of only one active replisome in a significant fraction of cells (>40%). Furthermore, many of the observed complexes have short lifetimes (<8 min), suggesting that replisome disassembly is quite prevalent, possibly occurring several times per cell cycle. The instability of the replisome complex is conflict-induced: transcription inhibition stabilizes these complexes, restoring the second replisome in many of the cells. Our results suggest that, in contrast to the canonical model, DNA replication is a largely discontinuous process in vivo due to pervasive replication-transcription conflicts.

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

Transcription leads to pervasive replisome instability in bacteria

The canonical model of DNA replication describes a highly-processive and largely continuous process by which the genome is duplicated. This continuous model is based upon in vitro reconstitution and in vivo ensemble experiments. Here, we characterize the replisome-complex stoichiometry and dynamics with single-molecule resolution in bacterial cells. Strikingly, the stoichiometries of the replicative helicase, DNA polymerase, and clamp loader complexes are consistent with the presence of only one active replisome in a significant fraction of cells (>40%). Furthermore, many of the observed complexes have short lifetimes (<8 min), suggesting that replisome disassembly is quite prevalent, possibly occurring several times per cell cycle. The instability of the replisome complex is conflict-induced: transcription inhibition stabilizes these complexes, restoring the second replisome in many of the cells. Our results suggest that, in contrast to the canonical model, DNA replication is a largely discontinuous process in vivo due to pervasive replication-transcription conflicts.