Next Article in Journal
Impact of Cellulose-Rich Organic Soil Amendments on Growth Dynamics and Pathogenicity of Rhizoctonia solani
Next Article in Special Issue
Functional Redundancy and Specialization of the Conserved Cold Shock Proteins in Bacillus subtilis
Previous Article in Journal
MALDI-TOF MS in a Medical Mycology Laboratory: On Stage and Backstage
Previous Article in Special Issue
Molecular Physiological Characterization of a High Heat Resistant Spore Forming Bacillus subtilis Food Isolate
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 6
10.3390/microorganisms9061284
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Tatiana Ermi
1,†,
Carmen Vallin
1,†,
Ana Gabriela Regalado García
2,
Moises Bravo
1,
Ismaray Fernandez Cordero
1,
Holly Anne Martin
1,
Mario Pedraza-Reyes
2 and
Eduardo Robleto
1,*
1
School of Life Sciences, University of Nevada, Las Vegas, 4505 S Maryland Pkwy, Las Vegas, NV 89154, USA
2
Department of Biology, Division of Natural and Exact Sciences, University of Guanajuato, P.O. Box 187, Guanajuato Gto. 36050, Mexico
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2021, 9(6), 1284; https://doi.org/10.3390/microorganisms9061284
Submission received: 12 May 2021 / Revised: 8 June 2021 / Accepted: 9 June 2021 / Published: 12 June 2021
(This article belongs to the Special Issue Bacillus subtilis as a Model Organism to Study Basic Cell Processes)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Experiments in the 1950s established that mutations are replication-dependent events that occur randomly during growth [1]. However, studies in the past few decades have uncovered transcription-dependent mutagenic mechanisms and their underlying factors. These types of mechanisms affect mutagenesis at specific genomic regions as a function of transcriptional activity [2,3,4,5,6]. One enzyme important in transcription-dependent mutagenesis is Mfd, the transcription-repair coupling factor [4,7,8]. Mfd gets recruited to stalled RNA Polymerase (RNAP) at damaged DNA sites and initiates high-fidelity transcription-coupled repair (TCR). Recent structural studies indicate that Mfd is a more dynamic protein than previously described and can function outside of TCR and in the absence of exogenous damage. Mfd was shown to bind [9] and translocate on double-stranded DNA [10] as well as to discriminate between different transcription elongation-stalling events [11]. Moreover, Mfd is proposed to be a pro-mutagenic/evolvability factor [12,13,14] that can associate with RNAP [9] and interact with different repair proteins to promote mutations [15,16]. More recently, it has been described as a factor affecting global transcription [17,18]. An exciting aspect raised by the transcription-induced mutagenesis mechanism is whether intrinsic sequences with the potential to form non-B DNA structures influence such a process in an Mfd-dependent manner.
Non-B DNA sequences can form structures that deviate from the canonical right-handed Watson and Crick structure. In vitro studies found that non-B DNA structures are promoted and stabilized by transcription [19,20] and can block the Mammalian and T7 phage RNA polymerases [21,22]. In eukaryotic cells, non-B DNA-forming sequences associate with genomic instability [23,24] as well as transcription and replication disruptions [25,26,27]. Many studies have shown that these types of events lead to recognition by and recruitment of several DNA repair factors such as CSB, the functional human homolog of Mfd, helicases, and error-prone polymerases [26,28,29]. However, we know less about the role of non-B DNA in bacteria, mainly whether it influences genome instability and the factors that promote or prevent mutagenesis at these sites.
Research on the influence of non-B DNA structures in mutagenesis in bacteria has focused on two types of sequence motifs predicted to form G-quadruplex (G4) DNA and hairpin structures [30,31]. A G4 DNA sequence has four equal guanine tracts separated by three loop regions: G3+NY1G3+NY2G3+NY3G3+. The four guanine tracts come together to form a tetrad plane via Hoogsteen base-pairing and stack vertically, creating a four-stranded structure [32]. Hairpin sequences are formed by inverted repeats, making a paired stem and a single-stranded DNA loop region. A recent study on E. coli growing cells found that the presence of a G4 DNA sequence in an actively transcribed reporter gene strongly increased mutation rates. Interestingly, they found that Hfq promoted this process by binding and stabilizing the G4 structure [31]. Hfq is a regulatory nucleic acid-binding protein whose function has broadened recently [33], and now extends to interactions with alternate DNA structures and a role in B. subtilis stationary-phase survival [34]. Similarly, hairpin sequences can promote genetic instability through slippage events that result in insertions or deletions (indels) during replication [35]. Notably, when reporter genes containing hairpin-forming sequences were transcriptionally active, mutation frequencies increased and correlated with promoter strength [36]. These studies in E. coli suggest that non-B DNA structures may have a role in mutagenesis, particularly in transcribed regions, but does this occur in other bacterial species? If so, what factors, besides Hfq, promote such a mechanism?
This study investigated if G4 DNA and hairpin-forming motifs influence stationary-phase mutagenesis (SPM) in Bacillus subtilis. B. subtilis is a Gram-positive sporulating bacterium adapted to feast and famine cycles in the soil’s upper layers [37]. Upon the onset of stress, such as starvation, B. subtilis cells activate gene expression programs, enter a non-growing state, and differentiate into distinct subpopulations [38]. Our previous observations have demonstrated that increased mutagenesis occurs in highly transcribed regions [39] and in distinct cell subpopulations [40,41]. Mfd also acts in coordination with nucleotide excision repair (NER), base excision repair (BER), and the transcription factor GreA to promote mutagenesis during stationary phase [15,42]. These observations make B. subtilis an attractive model to test our hypothesis that non-B DNA-forming motifs influence Mfd-dependent stationary-phase mutagenesis at transcribed regions.
We developed a system to measure non-B DNA’s influence 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 the amyE locus. To design these argF alleles, we used in silico tools to find endogenous non-B DNA sequences in the argF gene and introduced point mutations to the open reading frame to either increase or decrease their predicted structure stability. We then introduced a nonsense codon in a region corresponding to a loop in the structure, which resulted in a non-functional truncated ArgF (confers auxotrophy) and a system to measure reversions to Arg+ (Figure 1). Using this system, we found that sequences predicted to form non-B DNA structures promoted mutagenesis in stationary-phase cells; such a response did not occur in growing conditions. We also showed that Mfd promotes the formation of mutations within the predicted structures.

2. Materials and Methods

2.1. Strain Construction

The bacterial strains used in this study are shown in Table 1. For detailed strain construction, see supplemental material Non-B DNA strain construction. Briefly, Mfold [43] and QGRS mapper [44] were used to find sequences with the potential to form hairpin and G4 DNA structures, respectively, in the argF ORF. Once the sequences were found, synonymous point mutations were introduced into positions that were predicted to strengthen or destabilize the structure. A nonsense codon was then engineered into a region of the structure corresponding to a loop. For the +/−hairpin strains, we introduced an ochre nonsense codon at position 37Q (CAA→TAA) corresponding to the loop part of the predicted hairpin before disruption. For the +/−G4 strains, we introduced an amber nonsense codon at position 70Q (CAG→TAG) corresponding to the second loop of the predicted G4 structure. This nonsense codon resulted in a truncated ArgF protein and arginine auxotrophy. The new defective argF alleles were then cloned into pDR111 downstream of the IPTG-inducible Pspac promoter, transformed, and recombined into the amyE locus of the argF deletion strain CV4000 (see Figure 1). The Mfd− strains were constructed by transforming genomic DNA from YB9801 [8] into the +/−hairpin and +/−G4 strains and selecting on TBAB plates containing tetracycline.

2.2. Stationary-Phase Mutagenesis Assay

We performed a stationary-phase mutagenesis assay on each strain [47]. Strains were streaked out on TBAB plates with the appropriate antibiotic and incubated overnight at 37 °C. The next day, a colony from the TBAB plate was used to inoculate 2 mL of Penassay Broth (PAB). The cells were grown overnight at 37 °C with aeration. An aliquot of the overnight was then added to a 250 mL Erlenmeyer culture flask containing 10 mL of PAB medium and 10 μL of 1000× Ho-Le trace elements. Cells were grown at 37 °C with aeration (250 rpm) to stationary phase. Growth was monitored with a spectrophotometer measuring O.D. at 600 nm (OD600). Cells were harvested 90 min after cessation of exponential growth and centrifuged at 4000 rpm for 5 min, the supernatant was decanted, and the cells were resuspended in Spizizen minimal salts (SMS); this process was repeated twice. To determine initial cell titers, we serial diluted and spread plated 100 μL of cells on minimal media supplemented with 100 μg/mL arginine. Colonies were counted after 48 h of incubation at 37 °C. Then, 100 μL aliquots were plated on minimal media containing 0.1 mM IPTG and no arginine (n = 5). The plates were incubated for nine days at 37 °C. Each day plates were scored for the appearance of Arg+ colonies. The initial titers were used to normalize the cumulative number of revertants per day to the number of the CFU plated. To assay the background survival of cells, the viability was tracked by taking agar plugs every other day, serial diluting, and plating on minima media supplemented with 100 μg/mL arginine.

2.3. Arg+ Mutant Analysis

To characterize Arg+ mutants, we performed a tRNA suppressor analysis and DNA sequencing [47]. Suppressor analysis was conducted by first patching between 20 and 25 Arg+ colonies collected from days 5–9 of each replicate trial onto the same media used in the SPM experiment (minimal media lacking arginine supplemented with 0.1 mM IPTG) to eliminate transient Arg+ prototrophs. Previous studies determined that analyzing mutants that arose after day 5 would avoid selecting slow-growing mutants from growth [47]. Arg+ mutants that grew were stocked and then patched on minimal media lacking either methionine or histidine. Our background B. subtilis YB955 strain is also auxotrophic for methionine and histidine due to ochre and amber nonsense codon in the metB and hisC genes, respectively [47]. Growth of the Arg+ mutants on either trace met and trace his plates was scored after 24 h and 48 h. If growth occurred, those revertants were considered tRNA suppressors. For those Arg+ mutants determined not to be tRNA suppressors, we performed a colony PCR amplifying and sequencing the complete argF gene.

2.4. Fluctuation Test

We conducted a fluctuation test to determine the mutation rates for arginine prototrophy during growth for each strain. Briefly, an overnight aerated culture grown overnight in 2 mL of PAB at 37 °C were diluted 1:10 in 1× Spizizen minimal salts (SMS), vortexed, and diluted 1:10 again in this solution. Aliquots of 0.5 mL of the second 1× SMS solution were added to 49.5 mL of PAB supplemented with 0.1 mM IPTG. The cells were vortexed, and 1 mL of the 10−4 mixture was dispensed into 40 different 18 mm test tubes. Cells were grown for 12 h at 37 °C with aeration (250 rpm). Each test tube culture was pelleted, washed, resuspended in 100 μL of 1× SMS, and individually plated on minimal media lacking arginine with 0.1 mM IPTG. In addition, three individual test tubes were serial diluted and plated on minimal media supplemented with 100 μg/mL arginine to determine cell count. Plates were incubated at 37 °C for 48 h, and then colonies were counted. The mutation rates were determined using the MSS maximum likelihood method [48,49].

2.5. Statistical Analysis

To determined what type of significance test to use, we first checked the data for normality and variance using a Shapiro–Wilk test and Brown–Forsythe test, respectively. Significance was determined by performing a Student’s t-test to compare two strains or a one-way ANOVA with a post-hoc Tukey’s test to compare more than two strains.

2.6. Protein Alignment Analysis

The B. subtilis 319 amino acid ArgF protein sequence was inputted into pBLAST as the query sequence and blasted into the UniProtKB/Swiss-Prot (SwissProt) database. The output identified 100 ArgF amino acid sequences. The sequences were exported and aligned. Then, we analyzed position 37 and 70 using Weblogo. Weblogo reported the conservation at that site. In addition, Weblogo reported the frequency of each amino acid found at that site within the 100 ArgF sequences.

3. Results

3.1. Non-B DNA-Forming Motifs Promote Mutagenesis

3.1.1. Hairpins and G4 DNA-Forming Motifs Accumulate Mutations in Stationary Phase B. subtilis Cells

Our stationary-phase mutagenesis assay revealed that the strains containing argF alleles with sequences predicted to form non-B DNA structures accumulated significantly more mutants over nine days than those argF alleles less likely to form such structures (Figure 2). Specifically, by day 6, the +G4 strain containing the sequence G3N(1–13)G3N(1–13)G3N(1–13)G3 accumulated ~3 times more mutants than the −G4, which was not predicted to form any G4 structure. Arg+ accumulation for the +G4 strain continued after day 6, but the number of revertants became too numerous to count compared to what was observed in the −G4. The +hairpin (120 b) strain showed a significant increase (~two fold) in the accumulation of Arg+ colonies compared to the −hairpin. Since window size is an important parameter for predicting hairpin structures [36], our approach included disrupting a hairpin structure using a smaller window size in the same region of argF. We found mutagenesis of the −hairpin (40 b) strain indistinguishable from the W.T. in our experiments (Supplemental Figure S6). In addition, mutagenic differences observed in all the strains were not the result of viability differences during the experiment (Supplementary Figure S7). Interestingly, the strains differing in G4 stability had more mutants over time compared to the strains differing in hairpin stability. These results are consistent with our previous report [47] and we attribute this difference to the type of nonsense codon used in the different constructs (see Supplemental TAG vs. TAA section).

3.1.2. Analysis of the Arg+ Population

Analysis of a sample of the Arg+ colonies that arose after day 5 indicated that two categories of mutants restored arginine prototrophy in the tested strains: (1) nonsense suppressor tRNAs and (2) mutations in argF (Figure 3A,D). Nonsense tRNA suppressors represented <50% of the Arg+ colonies sampled. Sequencing of non-tRNA suppressor mutants showed that mutations in argF were limited to base substitutions that changed the stop codon’s first or third position (Figure 3B,E). We did not observe deletion events in our system.
In the strains differing in hairpin stability with the ochre nonsense codon TAA at position 37, transitions in the first position from T⟶C restored the W.T. polar glutamine-Q codon, and T⟶G transversions resulted in the acidic glutamate-E codon (Figure 3B). A protein analysis that looked at 100 ArgF sequences across different phyla confirmed that this codon position had some flexibility and that glutamate-E was indeed found at this position in some bacterial species (Figure 3C). Only transitions in the first positions were observed in the +hairpin strain, and no mutations occurred in the second or third position in either of the hairpin strains. Interestingly, the proportion of tRNA suppressor mutants appeared higher in our sample of the +hairpin than in the −hairpin one (Figure 3A).
In the strains differing in G4 stability carrying the amber nonsense codon TAG at position 70, revertants displayed T→C transitions in the first codon position, which restored the W.T. polar glutamine-Q triplet, and G→T transversions in the third codon position, which resulted in the aromatic amino acid tyrosine-Y triplet (Figure 3E). Protein analysis of this codon position found that it exhibited some flexibility (Figure 3F). Additionally, glutamate-E was the second most frequent amino acid found at this position in other ArgF proteins, while tyrosine was not observed in our sample of 100 ArgF proteins. (Figure 3F). Moreover, the proportions of the revertants with these genetic changes were similar in both strains (Figure 3E). The proportion of tRNA suppressors was similar in our sample of the Arg+ population between the +/−G4 strains (Figure 3D).

3.2. Hairpin- and G4 DNA-Forming Motifs Do Not Influence Mutation Rates in B. subtilis

To test whether the results observed were exclusive to stationary-phase conditions, we performed a fluctuation test on the strains employed in this study. We found that the +/−hairpin and +/−G4 DNA’s predicted stability did not correlate with mutation rates in B. subtilis growing cells. Additionally, we did not observe differences between the sets of strains as we did in stationary phase (hairpin vs. G4) (Figure 4). This result suggests that cell physiology influences the role non-B DNA structures have in mutagenesis.

3.3. Mfd Promotes Mutations at Non-B DNA Sequences

3.3.1. Accumulation of Arg+ Mutations at Non-B DNA Sequences Was Decreased in the Absence of Mfd

We tested if the transcription-repair coupling factor Mfd influenced SPM at non-B DNA sequences in B. subtilis. We transformed the mfd mutation into the strains carrying argF alleles differing in their predicted ability to form either hairpin or G4 DNA structures and performed a stationary-phase mutagenesis assay. Disrupting mfd significantly decreased the accumulation of mutations by ~two-fold in the +hairpin and +G4 strains (Figure 5). This response was not observed in the −hairpin, and there was a small but significant increase in the accumulation of Arg+ mutants in the −G4 strain deficient in Mfd (Figure 5B). These results were not due to differences in viability (Supplementary Figure S8).

3.3.2. Analysis of the Arg+ Population in Mfd− Cells

Subsequent analysis of the Arg+ mutants found that Mfd deficiency abolished mutations in argF as almost 100% of the mutants sampled were nonsense tRNA suppressors (Figure 6A,D). Thus, the spectrum of mutations changed drastically in the Mfd deficient strains. We estimated the total number of tRNA suppressors in this Arg+ population based on the sampling of mutants that arose at days 5–9 during the SPM trials. This analysis suggested that nonsense tRNA suppressors significantly increase in the +/−G4 lacking but not in the +/−hairpin strains in the Mfd− background (Figure 6C,F). The genetic changes that confer arginine prototrophy in argF (Figure 3E) suggest that mutagenesis of tyrosine (first position in the G/ATA anticodon) and glutamine (third position in the TTG anticodon) tRNA genes are increased in the absence of Mfd.

4. Discussion

4.1. An In Vivo System to Measure the Effects of Non-B DNA-Forming Motifs on Bacterial Mutagenesis

We investigated if non-B DNA-forming motifs, such as hairpins and G4 DNA sequences, contribute to stationary-phase mutagenesis in B. subtilis. First, we exploited the fact that this bacterium does not have a strong codon bias [50,51,52] to develop a system that altered endogenous non-B DNA sequences in the argF gene by introducing point mutations that increase or decrease their predicted stability. We then tested the effect of these changes on mutagenesis during growth and stationary phase.
Our approach to construct the non-B DNA sequence variants included the use of in silico tools that predict their stability based on extensive thermodynamic studies. We used Mfold [43] and QGRS mapper [44] to measure hairpins and G4 DNA stability, respectively. Previous in vitro [21] and in vivo studies [31] used optimized sequence motifs in tandem arrangements. These arrangements are not commonly found in coding regions, likely due to functional constraints. In contrast, we worked with an endogenous system by altering sequence motifs within a coding region. The G4 motif used in our experiments had longer asymmetrical loops (G3N(1–13)G3N(1–13)G3N(1–13)G3) than those commonly used (G3N(7)G3N(7)G3N(7)G3), and thermodynamic studies have shown that sequences having asymmetrical loops can still form stable G-quadruplex structures [24,53]. Additionally, G4 DNA containing a single-nucleotide loop, like the one used in this study, was shown to be thermally stable and to incite genomic instability in eukaryotic genomes [24]. In this regard, we performed a G4 motif search across the B. subtilis genome and found that the canonical G4 motif occurred four times in coding regions, while the G4 motif used here appeared in 39 genes (Supplemental Table S2). Further characterization of these genes showed they were distributed throughout the B. subtilis chromosome and many coded for membrane-bound proteins (Supplemental Figures S9 and S10).
To design the constructs that tested the effect of hairpins within coding regions on SPM, we used two window sizes because this parameter was reported to affect their stability, particularly in the context of their transcription [36]. High levels of transcription generate more supercoiling, which increases the likelihood of forming single-stranded DNA and intrastrand base pairing between distant DNA residues [54,55]. We found that disrupting a hairpin using a 40 b window size did not affect mutation levels (Supplemental Figure S6); in contrast, using a 120 b window size did. The IPTG-inducible Pspac promoter used in our study is robust, suggesting that the design using a larger window size provided a better parameter for forming single-stranded hairpins. Studies have used window sizes ranging from 30–300 bases [25,36]. In summary, the features we used to design the non-B DNA forming sequences are consistent with what is observed within gene coding regions and account for conditions in which a region is robustly transcribed.

4.2. Non-B DNA-Forming Motifs in Transcribed Coding Regions Promote Mutations in Nutritionally Stressed Cells

Early in silico analysis in bacteria found that non-B DNA forming sequences were enriched in noncoding regions, including promoters and untranslated regions (UTR). This finding led researchers to focus on elucidating the role of these structures in gene regulation [56,57]. Indeed, a study in a polyextremophile, Deinococcus radiodurans, which contains G4-motifs in promoter regions, found that adding a G4 stabilizing ligand led to repression of DNA repair genes and decreased radioresistance [58].
The effect of DNA structures found in bacterial regulatory regions continues to be well studied [59]. However, very little is known about what happens to alternate DNA structures that are within an ORF. Wright and colleagues proposed that hairpin structures formed during transcription in coding regions were hotspots for mutagenesis in E. coli [30]. More recently, a study examined the role of G4 structures in genome instability in E. coli-transcribed regions and found that these sequences increased the rate of deletions [31]. Here, we investigated if G4 DNA- and hairpin-forming motifs as intrinsic elements that influence stationary-phase mutagenesis in B. subtilis. We found that non-B DNA forming sequences promoted mutagenesis in stationary-phase cells. Specifically, sequences increased or decreased in the likelihood to form a hairpin, and G4 DNA were increased and decreased in the accumulation of Arg+ mutants, respectively. Furthermore, unlike results in E. coli, this response was only observed in nutritionally stressed cells and did not occur in growing cells. This result underlines the contribution of mutagenic processes taking place in non-replicating conditions to the evolutionary process and the importance of using different model systems. Moreover, this work provides support to the concept that non-B DNA plays a role in bacterial transcription-dependent mutagenic processes.

4.3. Mfd Promotes Mutations at Non-B DNA Sequences

To identify factors governing the process by which non-B DNA-forming motifs promote the accumulation of mutations, we investigated if the transcription-repair coupling factor Mfd promoted mutations at non-B DNA sites. Mfd is a factor shown to promote transcription-dependent mutagenesis and is recruited to genomic regions containing non-B-DNA-forming motifs [17,28]. Disrupting mfd led to significant decreases in mutagenesis levels in the +hairpin and +G4 DNA strains. Subsequent analysis of the Arg+ population that arose in the absence of Mfd showed that the vast majority of revertants were tRNA suppressors in all strains; argF mutants were abolished in the absence of Mfd. This result is consistent with our previous studies in stationary-phase mutagenesis in B. subtilis [8,60]. This result further supports the promutagenic role of Mfd at transcribed regions and suggests that non-B DNA sites in bacteria influence Mfd-dependent mutagenesis.
Mfd associates with RNAP when transcription elongation is blocked by a bulky DNA lesion. Upon contacting the blocked RNAP, Mfd undergoes a series of conformational changes that recruit nucleotide excision repair proteins [11]. However, recent single-molecule resolution and structural studies suggest that Mfd and RNAP interact, even in the absence of DNA damage, and that Mfd works to modulate transcription [9,10]. In this context, in vivo and in vitro reports showed that non-B DNA structures formed in the DNA or RNA could interrupt transcription, especially when they are found in the coding strand [21,61,62,63]. In B. subtilis, loss of Mfd led to increases in RNAP association at genes enriched with hairpins [17]. Similarly, deficiencies in CSB, the functional homolog of Mfd, caused RNAPII pausing at G4 DNA motifs in a neuroblastoma cell line [28]. Considering these observations, one can propose a model in which transcription promotes the formation of alternate structures, through changes in supercoiling, in the DNA or mRNA with the potential to impede transcription. We speculate that the stability of such structures defines whether the impediment is resolved by Mfd through transcription rescue (low structure stability) or disassembly of the elongation complex (high structure stability) and recruitment of error-prone repair of DNA. In E. coli, Mfd promotes mutagenesis in stressed cells via the formation of R-loops and recombination intermediates that generate point mutations or amplifications in a Lac+ mutagenesis system [4]. This system requires recombination functions and activation of the sigma S regulon (amplification) and the DinB (point mutations) [64]. Unlike in E. coli, Mfd-dependent mutagenesis in B. subtilis operates independently of recombination functions [47] and can operate by interacting with components of the NER and BER systems [15,16,65]. It is also possible that the loss of Mfd causes depression in mutagenesis indirectly. A recent RNAseq study showed profound changes in the transcriptome of Mfd− cells compared to wild-type stationary-phase cells. Almost half the genes in the B. subtilis genome were dysregulated, including factors previously identified to promote SPM, such as GreA [42], PolY [66], PolX [67], and others [15]. Therefore, it is possible that the Mfd effects observed in mutagenesis are the indirect product of dysregulation of a factor(s), such as a helicase or error-prone polymerase involved in mutagenesis at non-B DNA forming sequences. Interestingly, a previous study identified Hfq as a factor that promotes mutation at G4-DNA motifs in E. coli growing cells [31]. Like Mfd, the role of Hfq has been broadening [33], and like Mfd, cells deficient in Hfq have transcriptome changes and decreased survival during stress in B. subtilis [34]. Future work will determine if Hfq, along with Mfd, contribute to stationary-phase mutagenesis at non-B DNA motifs.

5. Conclusions

Non-B DNA sequences are present in all domains of life. These structures have a role in genome instability and disease in humans, but we know less about their function in bacteria. Here, we showed that non-B DNA-forming motifs promote genetic variation in B. subtilis, particularly in transcribed coding regions in stressed cells in an Mfd-dependent manner. Future work studying SPM in B. subtilis will discern between the competing but not necessarily exclusive hypotheses that Mfd works directly or indirectly to promote the formation of mutations at non-B DNA sequences in actively transcribed genes. This type of study provides insight into transcription-induced mutagenic mechanisms limiting mutations to times of stress and defined genomic regions, and it improves our view of evolution.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9061284/s1, Figure S1: Prediction of hairpin structures across the argF gene using a sliding window analysis and in silico tools. Figure S2: Hairpin-forming motifs and predicted structures of the WT and +/−Hairpin (120) strains. Figure S3: Hairpin-forming motif and predicted structure of WT and −Hairpin (40) strains. Figure S4: G4-forming motif and predicted structure of +G4 and −G4 strains. Figure S5: Codon analysis using the Kazusa Codon Usage Database of regions in argF that were changed to disrupt or stabilize an endogenous non-B DNA structure. Figure S6: The −Hairpin (40bp strain) is indistinguishable from the W.T. Figure S7: Viability of the non-revertant population for the non-B DNA strains during the SPM assay. Figure S8: Viability of the non-revertant population for the non-B DNA strains and their Mfd deficient counterpart during the SPM assay. Figure S9: A. Schematic showing the chromosomal distribution of the 39 genes that contained the G4 motif used in this study. Figure S10: Modified figure from Bsucyc.org showing the cellular function of the 39 genes that contained the G4 motif used in this study. Table S1: QGRS Pattern search in argF ORF. Table S2: QGRS Pattern search in B. subtilis coding regions. Table S3: QGRS Results. Table S4: Primers used in this study. Table S5: Summary of Phenotype Test.

Author Contributions

Conceptualization, C.V., T.E., M.P.-R. and E.R.; methodology, H.A.M., A.G.R.G., and M.B.; formal analysis, C.V., H.A.M., and T.E.; investigation, C.V., T.E., M.B., and I.F.C.; writing—original draft preparation, T.E. and C.V.; writing—review and editing, T.E., C.V., H.A.M., M.P.-R., and E.R.; visualization, C.V.; supervision, E.R.; project administration, E.R.; funding acquisition, C.V., E.R., and M.P.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation (NSF-GRFP), the National Institutes of Health (GM131410), and the Consejo Nacional de Ciencia y Tecnología (A-1S-27116).

Data Availability Statement

Data is contained within the article or supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luria, S.E.; Delbrück, M. Mutations of Bacteria from Virus Sensitivity to Virus Resistance. Genetics 1943, 28, 491–511. [Google Scholar] [CrossRef]
  2. Reimers, J.M.; Schmidt, K.H.; Longacre, A.; Reschke, D.K.; Wright, B.E. Increased transcription rates correlate with increased reversion rates in leuB and argH Escherichia coli auxotrophs. Microbiology 2004, 150, 1457–1466. [Google Scholar] [CrossRef] [Green Version]
  3. Lang, K.S.; Hall, A.N.; Merrikh, C.N.; Ragheb, M.; Tabakh, H.; Pollock, A.J.; Woodward, J.J.; Dreifus, J.E.; Merrikh, H. Replication-Transcription Conflicts Generate R-Loops that Orchestrate Bacterial Stress Survival and Pathogenesis. Cell 2017, 170, 787–799. [Google Scholar] [CrossRef] [Green Version]
  4. Wimberly, H.; Shee, C.; Thornton, P.C.; Sivaramakrishnan, P.; Rosenberg, S.M.; Hastings, P.J. R-loops and nicks initiate DNA breakage and genome instability in non-growing Escherichia coli. Nat. Commun. 2013, 4, 2115. [Google Scholar] [CrossRef] [Green Version]
  5. Merrikh, H. Spatial and Temporal Control of Evolution through Replication–Transcription Conflicts. Trends Microbiol. 2017, 25, 515–521. [Google Scholar] [CrossRef] [PubMed]
  6. Million-Weaver, S.; Samadpour, A.N.; Moreno-Habel, D.A.; Nugent, P.; Brittnacher, M.J.; Weiss, E.; Hayden, H.S.; Miller, S.I.; Liachko, I.; Merrikh, H. An underlying mechanism for the increased mutagenesis of lagging-strand genes in Bacillus subtilis. Proc. Natl. Acad. Sci. USA 2015, 112, E1096–E1105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Paul, S.; Million-Weaver, S.; Chattopadhyay, S.; Sokurenko, E.V.; Merrikh, H. Accelerated gene evolution through replication–transcription conflicts. Nature 2013, 495, 512–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ross, C.; Pybus, C.; Pedraza-Reyes, M.; Sung, H.-M.; Yasbin, R.E.; Robleto, E. Novel Role of mfd: Effects on Stationary-Phase Mutagenesis in Bacillus subtilis. J. Bacteriol. 2006, 188, 7512–7520. [Google Scholar] [CrossRef] [Green Version]
  9. Ho, H.N.; Van Oijen, A.M.; Ghodke, H. The transcription-repair coupling factor Mfd associates with RNA polymerase in the absence of exogenous damage. Nat. Commun. 2018, 9, 1570. [Google Scholar] [CrossRef] [Green Version]
  10. Le, T.T.; Yang, Y.; Tan, C.; Suhanovsky, M.M.; Fulbright, R.M.; Inman, J.T.; Li, M.; Lee, J.; Perelman, S.; Roberts, J.W.; et al. Mfd Dynamically Regulates Transcription via a Release and Catch-Up Mechanism. Cell 2018, 172, 344–357. [Google Scholar] [CrossRef] [Green Version]
  11. Kang, J.Y.; Llewellyn, E.; Chen, J.; Olinares, P.D.B.; Brewer, J.; Chait, B.T.; Campbell, E.A.; Darst, S.A. Structural basis for transcription complex disruption by the Mfd translocase. Elife 2021, 10, 1–28. [Google Scholar] [CrossRef]
  12. Ragheb, M.; Thomason, M.K.; Hsu, C.; Nugent, P.; Gage, J.; Samadpour, A.N.; Kariisa, A.; Merrikh, C.N.; Miller, S.I.; Sherman, D.R.; et al. Inhibiting the Evolution of Antibiotic Resistance. Mol. Cell 2019, 73, 157–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Han, S.; Gong, Z.; Liang, T.; Chen, Y.; Xie, J. The role of Mfd in Mycobacterium tuberculosis physiology and underlying regulatory network. Microbiol. Res. 2021, 246, 126718. [Google Scholar] [CrossRef]
  14. Ukkivi, K.; Kivisaar, M. Involvement of transcription-coupled repair factor Mfd and DNA helicase UvrD in mutational processes in Pseudomonas putida. DNA Repair 2018, 72, 18–27. [Google Scholar] [CrossRef]
  15. Gómez-Marroquín, M.; Martin, H.A.; Pepper, A.; Girard, M.E.; Kidman, A.A.; Vallin, C.; Yasbin, R.E.; Pedraza-Reyes, M.; Robleto, E.A. Stationary-Phase Mutagenesis in Stressed Bacillus subtilis Cells Operates by Mfd-Dependent Mutagenic Pathways. Genes 2016, 7, 33. [Google Scholar] [CrossRef]
  16. Suárez, V.P.; Martínez, L.E.; Leyva-Sánchez, H.C.; Valenzuela-García, L.I.; Lara-Martínez, R.; Jiménez-García, L.F.; Ramírez-Ramírez, N.; Obregon-Herrera, A.; Cuéllar-Cruz, M.; Robleto, E.A.; et al. Transcriptional coupling and repair of 8-OxoG activate a RecA—Dependent checkpoint that controls the onset of sporulation in Bacillus subtilis. Sci. Rep. 2021, 11, 1–12. [Google Scholar] [CrossRef] [PubMed]
  17. Ragheb, M.N.; Merrikh, C.; Browning, K.; Merrikh, H. Mfd regulates RNA polymerase association with hard-to-transcribe regions in vivo, especially those with structured RNAs. Proc. Natl. Acad. Sci. USA 2021, 118, 1–10. [Google Scholar] [CrossRef] [PubMed]
  18. Martin, H.A.; Sundararajan, A.; Ermi, T.S.; Heron, R.; Gonzales, J.; Lee, K.; Anguiano-Mendez, D.; Schilkey, F.; Pedraza-Reyes, M.; Robleto, E.A. Mfd Affects Global Transcription and the Physiology of Stressed Bacillus subtilis Cells. Front. Microbiol. 2021, 12, 1–14. [Google Scholar] [CrossRef]
  19. Zhang, C.; Liu, H.-H.; Zheng, K.-W.; Hao, Y.-H.; Tan, Z. DNA G-quadruplex formation in response to remote downstream transcription activity: Long-range sensing and signal transducing in DNA double helix. Nucleic Acids Res. 2013, 41, 7144–7152. [Google Scholar] [CrossRef]
  20. Duquette, M.L.; Handa, P.; Vincent, J.A.; Taylor, A.F.; Maizels, N. Intracellular transcription of G-rich DNAs induces formation of G-loops, novel structures containing G4 DNA. Genes Dev. 2004, 18, 1618–1629. [Google Scholar] [CrossRef] [Green Version]
  21. Tornaletti, S.; Park-Snyder, S.; Hanawalt, P.C. G4-forming Sequences in the Non-transcribed DNA Strand Pose Blocks to T7 RNA Polymerase and Mammalian RNA Polymerase II. J. Biol. Chem. 2008, 283, 12756–12762. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pandey, S.; Ogloblina, A.M.; Belotserkovskii, B.P.; Dolinnaya, N.; Yakubovskaya, M.G.; Mirkin, S.M.; Hanawalt, P.C. Transcription blockage by stable H-DNA analogs in vitro. Nucleic Acids Res. 2015, 43, 6994–7004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Du, X.; Gertz, E.M.; Wójtowicz, D.; Zhabinskaya, D.; Levens, D.; Benham, C.J.; Schäffer, A.A.; Przytycka, T.M. Potential non-B DNA regions in the human genome are associated with higher rates of nucleotide mutation and expression variation. Nucleic Acids Res. 2014, 42, 12367–12379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Piazza, A.; Adrian, M.; Samazan, F.; Heddi, B.; Hamon, F.; Serero, A.; Lopes, J.; Teulade-Fichou, M.; Phan, A.T.; Nicolas, A. Short loop length and high thermal stability determine genomic instability induced by G-quadruplex-forming minisatellites. EMBO J. 2015, 34, 1718–1734. [Google Scholar] [CrossRef] [Green Version]
  25. Szlachta, K.; Thys, R.G.; Atkin, N.D.; Pierce, L.C.T.; Bekiranov, S.; Wang, Y.-H. Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human. Genome Biol. 2018, 19, 89. [Google Scholar] [CrossRef]
  26. Northam, M.R.; Moore, E.A.; Mertz, T.M.; Binz, S.K.; Stith, C.M.; Stepchenkova, E.I.; Wendt, K.L.; Burgers, P.M.J.; Shcherbakova, P.V. DNA polymerases ζ and Rev1 mediate error-prone bypass of non-B DNA structures. Nucleic Acids Res. 2014, 42, 290–306. [Google Scholar] [CrossRef] [Green Version]
  27. Teng, F.-Y.; Hou, X.-M.; Fan, S.-H.; Rety, S.; Dou, S.-X.; Xi, X.-G. Escherichia coli DNA polymerase I can disrupt G-quadruplex structures during DNA replication. FEBS J. 2017, 284, 4051–4065. [Google Scholar] [CrossRef] [Green Version]
  28. Scheibye-Knudsen, M.; Tseng, A.; Jensen, M.B.; Scheibye-Alsing, K.; Fang, E.F.; Iyama, T.; Bharti, S.K.; Marosi, K.; Froetscher, L.; Kassahun, H.; et al. Cockayne syndrome group A and B proteins converge on transcription-linked resolution of non-B DNA. Proc. Natl. Acad. Sci. USA 2016, 113, 12502–12507. [Google Scholar] [CrossRef] [Green Version]
  29. Huber, M.D.; Lee, D.C.; Maizels, N. G4 DNA unwinding by BLM and Sgs1p: Substrate specificity and substrate-specific inhibition. Nucleic Acids Res. 2002, 30, 3954–3961. [Google Scholar] [CrossRef] [Green Version]
  30. Wright, B.E.; Schmidt, K.H.; Davis, N.; Hunt, A.T.; Minnick, M.F., II. Correlations between secondary structure stability and mutation frequency during somatic hypermutation. Mol. Immunol. 2008, 45, 3600–3608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Parekh, V.J.; Niccum, B.A.; Shah, R.; Rivera, M.A.; Novak, M.J.; Geinguenaud, F.; Wien, F.; Arluison, V.; Sinden, R.R. Role of Hfq in Genome Evolution: Instability of G-Quadruplex Sequences in E. coli. Microorganisms 2019, 8, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Maizels, N.; Gray, L.T. The G4 Genome. PLoS Genet. 2013, 9, e1003468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ricardo, F.; Cecília, S.; José, M.A. New molecular interactions broaden the functions of the RNA chaperone Hfq. Curr. Genet. 2019, 65, 1313–1319. [Google Scholar]
  34. Rochat, T.; Delumeau, O.; Figueroa-Bossi, N.; Noirot, P.; Bossi, L.; Dervyn, E.; Bouloc, P. Tracking the Elusive Function of Bacillus subtilis Hfq. PLoS ONE 2015, 10, e0124977. [Google Scholar] [CrossRef] [Green Version]
  35. Sinden, R.R.; Zheng, G.X.; Brankamp, R.G.; Allen, K.N. On the deletion of inverted repeated DNA in Escherichia coli: Effects of length, thermal stability, and cruciform formation in vivo. Genetics 1991, 129, 991–1005. [Google Scholar] [CrossRef]
  36. Schmidt, K.H.; Reimers, J.M.; Wright, B.E. The effect of promoter strength, supercoiling and secondary structure on mutation rates in Escherichia coli. Mol. Microbiol. 2006, 60, 1251–1261. [Google Scholar] [CrossRef]
  37. Bernhardt, J.; Weibezahn, J.; Scharf, C.; Hecker, M. Bacillus subtilis During Feast and Famine: Visualization of the Overall Regulation of Protein Synthesis during Glucose Starvation by Proteome Analysis. Genome Res. 2003, 13, 224–237. [Google Scholar] [CrossRef] [Green Version]
  38. López, D.; Kolter, R. Extracellular signals that define distinct and coexisting cell fates in Bacillus subtilis. FEMS Microbiol. Rev. 2010, 34, 134–149. [Google Scholar] [CrossRef] [Green Version]
  39. Pybus, C.; Pedraza-Reyes, M.; Ross, C.A.; Martin, H.; Ona, K.; Yasbin, R.E.; Robleto, E. Transcription-Associated Mutation in Bacillus subtilis Cells under Stress. J. Bacteriol. 2010, 192, 3321–3328. [Google Scholar] [CrossRef] [Green Version]
  40. Martin, H.A.; Kidman, A.A.; Socea, J.; Vallin, C.; Pedraza-Reyes, M.; Robleto, E.A. The Bacillus subtilis K-State Promotes Stationary-Phase Mutagenesis via Oxidative Damage. Genes 2020, 11, 190. [Google Scholar] [CrossRef] [Green Version]
  41. Ambriz-Aviña, V.; Yasbin, R.E.; Robleto, E.A.; Pedraza-Reyes, M. Role of Base Excision Repair (BER) in Transcription-associated Mutagenesis of Nutritionally Stressed Nongrowing Bacillus subtilis Cell Subpopulations. Curr. Microbiol. 2016, 73, 721–726. [Google Scholar] [CrossRef] [Green Version]
  42. Leyva-Sánchez, H.C.; Villegas-Negrete, N.; Abundiz-Yañez, K.; Yasbin, R.E.; Robleto, E.A.; Pedraza-Reyes, M. Role of Mfd and GreA in Bacillus subtilis Base Excision Repair-Dependent Stationary-Phase Mutagenesis. J. Bacteriol. 2020, 202, 1–17. [Google Scholar] [CrossRef] [PubMed]
  43. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef] [PubMed]
  44. Kikin, O.; D’Antonio, L.; Bagga, P.S. QGRS Mapper: A web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res. 2006, 34, W676–W682. [Google Scholar] [CrossRef] [PubMed]
  45. Martin, H.A.; Porter, K.E.; Vallin, C.; Ermi, T.; Contreras, N.; Pedraza-Reyes, M.; Robleto, E.A. Mfd protects against oxidative stress in Bacillus subtilis independently of its canonical function in DNA repair. BMC Microbiol. 2019, 19, 26. [Google Scholar] [CrossRef] [PubMed]
  46. Itaya, M.; Kondo, K.; Tanaka, T. A neomycin resistance gene cassette selectable in a single copy state in the Bacillus subtilis chromosome. Nucleic Acids Res. 1989, 17, 4410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Sung, H.-M.; Yasbin, R.E. Adaptive, or Stationary-Phase, Mutagenesis, a Component of Bacterial Differentiation in Bacillus subtilis. J. Bacteriol. 2002, 184, 5641–5653. [Google Scholar] [CrossRef] [Green Version]
  48. Foster, P.L. Methods for Determining Spontaneous Mutation Rates. Methods Enzymol. 2006, 409, 195–213. [Google Scholar]
  49. Rosche, W.A.; Foster, P.L. Determining Mutation Rates in Bacterial Populations. Methods 2000, 20, 4–17. [Google Scholar] [CrossRef] [Green Version]
  50. Ogasawara, N. Markedly unbiased codon usage in Bacillus subtilis. Gene 1985, 40, 145–150. [Google Scholar]
  51. Li, G.-W.; Oh, E.; Weissman, J.S. The anti-Shine–Dalgarno sequence drives translational pausing and codon choice in bacteria. Nature 2012, 484, 538–541. [Google Scholar] [CrossRef]
  52. Shields, D.C.; Sharp, P.M. Synonymous codon usage in Bacillus subtilis reflects both translational selection and mutational biases. Nucleic Acids Res. 1987, 15, 8023–8040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Guédin, A.; Gros, J.; Alberti, P.; Mergny, J.-L. How long is too long? Effects of loop size on G-quadruplex stability. Nucleic Acids Res. 2010, 38, 7858–7868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kouzine, F.; Gupta, A.; Baranello, L.; Wojtowicz, D.; Ben-Aissa, K.; Liu, J.; Przytycka, M.T.; Levens, D. Transcription-dependent dynamic supercoiling is a short-range genomic force. Nat. Struct. Mol. Biol. 2013, 20, 396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zhabinskaya, D.; Benham, C.J. Competitive superhelical transitions involving cruciform extrusion. Nucleic Acids Res. 2013, 41, 9610–9621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Du, X.; Wójtowicz, D.; Bowers, A.A.; Levens, D.; Benham, C.J.; Przytycka, T.M. The genome-wide distribution of non-B DNA motifs is shaped by operon structure and suggests the transcriptional importance of non-B DNA structures in Escherichia coli. Nucleic Acids Res. 2013, 41, 5965–5977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Rawal, P.; Kummarasetti, V.B.R.; Ravindran, J.; Kumar, N.; Halder, K.; Sharma, R.; Mukerji, M.; Das, S.K.; Chowdhury, S. Genome-wide prediction of G4 DNA as regulatory motifs: Role in Escherichia coli global regulation. Genome Res. 2006, 16, 644–655. [Google Scholar] [CrossRef] [Green Version]
  58. Beaume, N.; Pathak, R.; Yadav, V.; Kota, S.; Misra, H.S.; Gautam, H.K.; Chowdhury, S. Genome-wide study predicts promoter-G4 DNA motifs regulate selective functions in bacteria: Radioresistance of D. radiodurans involves G4 DNA-mediated regulation. Nucleic Acids Res. 2012, 41, 76–89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Dey, U.; Sarkar, S.; Teronpi, V.; Yella, V.R.; Kumar, A. G-quadruplex motifs are functionally conserved in cis-regulatory regions of pathogenic bacteria: An in-silico evaluation. Biochimie 2021, 184, 40–51. [Google Scholar] [CrossRef]
  60. Vidales, L.E.; Cárdenas, L.C.; Robleto, E.; Yasbin, R.E.; Pedraza-Reyes, M. Defects in the error prevention oxidized guanine system potentiate stationary-phase mutagenesis in Bacillus subtilis. J. Bacteriol. 2009, 191, 506–513. [Google Scholar] [CrossRef] [Green Version]
  61. Agarwal, T.; Roy, S.; Kumar, S.; Chakraborty, T.K.; Maiti, S. In the Sense of Transcription Regulation by G-Quadruplexes: Asymmetric Effects in Sense and Antisense Strands. Biochemistry 2014, 53, 3711–3718. [Google Scholar] [CrossRef]
  62. Eddy, J.; Vallur, A.C.; Varma, S.; Liu, H.; Reinhold, W.C.; Pommier, Y.; Maizels, N. G4 motifs correlate with promoter-proximal transcriptional pausing in human genes. Nucleic Acids Res. 2011, 39, 4975–4983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Endoh, T.; Kawasaki, Y.; Sugimoto, N. Suppression of Gene Expression by G-Quadruplexes in Open Reading Frames Depends on G-Quadruplex Stability. Angew. Chem. 2013, 125, 5522–5526. [Google Scholar] [CrossRef] [PubMed]
  64. Rosenberg, S.M.; Shee, C.; Frisch, R.L.; Hastings, P.J. Stress-induced mutation via DNA breaks in Escherichia coli: A molecular mechanism with implications for evolution and medicine. BioEssays 2012, 34, 885–892. [Google Scholar] [CrossRef] [Green Version]
  65. Gómez-Marroquín, M.; Vidales, L.E.; Debora, B.N.; Santos-Escobar, F.; Obregón-Herrera, A.; Robleto, E.A.; Pedraza-Reyes, M. Role of Bacillus subtilis DNA Glycosylase MutM in Counteracting Oxidatively Induced DNA Damage and in Stationary-Phase-Associated Mutagenesis. J. Bacteriol. 2015, 197, 1963–1971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Sung, H.M.; Yeamans, G.; Ross, C.a.; Yasbin, R.E. Roles of YqjH and YqjW, homologs of the Escherichia coli UmuC/DinB or Y superfamily of DNA polymerases, in stationary-phase mutagenesis and UV-induced mutagenesis of Bacillus subtilis. J. Bacteriol. 2003, 185, 2153–2160. [Google Scholar] [CrossRef] [Green Version]
  67. Barajas-Ornelas, R.D.C.; Ramírez-Guadiana, F.H.; Juárez-Godínez, R.; Ayala-García, V.M.; Robleto, E.A.; Yasbin, R.E.; Pedraza-Reyes, M. Error-Prone Processing of Apurinic/Apyrimidinic (AP) Sites by PolX Underlies a Novel Mechanism That Promotes Adaptive Mutagenesis in Bacillus subtilis. J. Bacteriol. 2014, 196, 3012–3022. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 09 01284 g001 550
Figure 1. Schematic of the stationary-phase assay. (A) The argF gene was replaced with a selectable neomycin marker at its chromosomal locus. Then, IPTG-inducible argF alleles containing nonsense codons and differing in their potential to form non-B DNA structures were introduced into an ectopic chromosomal region. (B) Stationary-phase mutagenesis assay. Cells are grown to stationary phase and plated on minimal media lacking arginine. Then, we scored Arg+ mutants over time.
Figure 1. Schematic of the stationary-phase assay. (A) The argF gene was replaced with a selectable neomycin marker at its chromosomal locus. Then, IPTG-inducible argF alleles containing nonsense codons and differing in their potential to form non-B DNA structures were introduced into an ectopic chromosomal region. (B) Stationary-phase mutagenesis assay. Cells are grown to stationary phase and plated on minimal media lacking arginine. Then, we scored Arg+ mutants over time.
Microorganisms 09 01284 g001
Microorganisms 09 01284 g002 550
Figure 2. Hairpins and G4 DNA-forming motifs promote mutations in stationary phase B. subtilis cells. (A) Stationary-phase mutagenesis assay results for +/−hairpin strains. (B) Mfold predicted hairpin structures for 120 bp region in argF (positions 82–202). (C) Stationary-phase mutagenesis assay results for +/−G4 strains. (D) Schematic of predicted G4 DNA structures and G-score in argF (position 182–219). Please note the differences in scale between the two graphs. Bars represent the mean, and error bars are SEM (n = 3). Statistical significance was determined using a Student’s t-test * p < 0.05, *** p < 0.001.
Figure 2. Hairpins and G4 DNA-forming motifs promote mutations in stationary phase B. subtilis cells. (A) Stationary-phase mutagenesis assay results for +/−hairpin strains. (B) Mfold predicted hairpin structures for 120 bp region in argF (positions 82–202). (C) Stationary-phase mutagenesis assay results for +/−G4 strains. (D) Schematic of predicted G4 DNA structures and G-score in argF (position 182–219). Please note the differences in scale between the two graphs. Bars represent the mean, and error bars are SEM (n = 3). Statistical significance was determined using a Student’s t-test * p < 0.05, *** p < 0.001.
Microorganisms 09 01284 g002
Microorganisms 09 01284 g003 550
Figure 3. Arg+ revertant analysis. (A–C) +/−Hairpin results; (D–F) +/−G4 results. (A,D) The proportion of tRNA suppressors from Arg+ population sampled that arose from day 5–9. (B,E) argF sequencing results from Arg+ mutants that were determined not to be tRNA suppressors. (C,F) Protein alignment analysis of positions containing the nonsense codons 37 and 70.
Figure 3. Arg+ revertant analysis. (A–C) +/−Hairpin results; (D–F) +/−G4 results. (A,D) The proportion of tRNA suppressors from Arg+ population sampled that arose from day 5–9. (B,E) argF sequencing results from Arg+ mutants that were determined not to be tRNA suppressors. (C,F) Protein alignment analysis of positions containing the nonsense codons 37 and 70.
Microorganisms 09 01284 g003
Microorganisms 09 01284 g004 550
Figure 4. Hairpin (A) and G4 DNA (B) forming motifs did not affect mutation rates in growing B. subtilis cells. Bars represent mutation rates using the MSS maximum likelihood method, and error bars are confidence limits (CL95).
Figure 4. Hairpin (A) and G4 DNA (B) forming motifs did not affect mutation rates in growing B. subtilis cells. Bars represent mutation rates using the MSS maximum likelihood method, and error bars are confidence limits (CL95).
Microorganisms 09 01284 g004
Microorganisms 09 01284 g005 550
Figure 5. Mfd promotes mutations at hairpin (A) and G4 DNA (B) forming sequences. Bars represent the mean of at least three SPM trials, and error bars are SEM. Please note the differences in scales on the Y axis between graphs. Significance was determined using a Student’s t-test * p < 0.05.
Figure 5. Mfd promotes mutations at hairpin (A) and G4 DNA (B) forming sequences. Bars represent the mean of at least three SPM trials, and error bars are SEM. Please note the differences in scales on the Y axis between graphs. Significance was determined using a Student’s t-test * p < 0.05.
Microorganisms 09 01284 g005
Microorganisms 09 01284 g006 550
Figure 6. Mfd promotes mutagenesis at argF. (A–C) +/−Hairpin results; (D–F) +/−G4 results. (A,D) The proportion of tRNA suppressors from Arg+ population sampled that arose from day 5–9. (B,E) Arg+ mutants that accumulated from day 5–9. (C,F) tRNA suppressors that arose in the Arg+ population after day 5. Bars represent the mean with SEM. Significance was determined using ANOVA; * denotes p < 0.05, ** denotes p < 0.01, *** denotes p < 0.001, and **** denotes p < 0.0001.
Figure 6. Mfd promotes mutagenesis at argF. (A–C) +/−Hairpin results; (D–F) +/−G4 results. (A,D) The proportion of tRNA suppressors from Arg+ population sampled that arose from day 5–9. (B,E) Arg+ mutants that accumulated from day 5–9. (C,F) tRNA suppressors that arose in the Arg+ population after day 5. Bars represent the mean with SEM. Significance was determined using ANOVA; * denotes p < 0.05, ** denotes p < 0.01, *** denotes p < 0.001, and **** denotes p < 0.0001.
Microorganisms 09 01284 g006
Table
Table 1. Strains and plasmids used in this study.
Table 1. Strains and plasmids used in this study.
Strain or PlasmidRelevant GenotypeSource or Reference
CV1000 (+Hairpin)metB5, hisC952, leuC427, argF::neo,
amyE::pHS− argF+SLS		This study
CV2000 (−Hairpin)metB5, hisC952, leuC427, argF::neo,
amyE::pHS− argF−SLS		This study
CV1001 (+Hairpin Mfd−)metB5, hisC952, leuC427, argF::neo, mfd::tc,
amyE::pHS− argF+SLS		This study
CV2009 (−Hairpin Mfd−)metB5, hisC952, leuC427, argF::neo, mfd::tc,
amyE::pHS− argF−SLS		This study
TE300 (+G4)metB5, hisC952, leuC427, argF::neo,
amyE::pHS− argF+G4		This study
TE302 (−G4)metB5, hisC952, leuC427, argF::neo,
amyE::pHS− argF−G4		This study
TE400 (+G4 Mfd−)metB5, hisC952, leuC427, argF::neo, mfd::tc
amyE::pHS− argF+G4		This study
YB9801 (Mfd−)metB5, hisC952, leuC427, mfd::tet[8]
CV4000 (argF-)metB5, hisC952, leuC427, argF::neo[45]
TE402 (−G4 Mfd−)metB5, hisC952, leuC427, argF::neo, mfd::tc
amyE::pHS−argF−G4		This study
pBEST502neomycin resistance (Nmr) gene[46]
pDR111pdr111 amyE-hyper-SPANK (spec)Rudner lab
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ermi, T.; Vallin, C.; García, A.G.R.; Bravo, M.; Cordero, I.F.; Martin, H.A.; Pedraza-Reyes, M.; Robleto, E. Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis. Microorganisms 2021, 9, 1284. https://doi.org/10.3390/microorganisms9061284

AMA Style

Ermi T, Vallin C, García AGR, Bravo M, Cordero IF, Martin HA, Pedraza-Reyes M, Robleto E. Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis. Microorganisms. 2021; 9(6):1284. https://doi.org/10.3390/microorganisms9061284

Chicago/Turabian Style

Ermi, Tatiana, Carmen Vallin, Ana Gabriela Regalado García, Moises Bravo, Ismaray Fernandez Cordero, Holly Anne Martin, Mario Pedraza-Reyes, and Eduardo Robleto. 2021. "Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis" Microorganisms 9, no. 6: 1284. https://doi.org/10.3390/microorganisms9061284

APA Style

Ermi, T., Vallin, C., García, A. G. R., Bravo, M., Cordero, I. F., Martin, H. A., Pedraza-Reyes, M., & Robleto, E. (2021). Non-B DNA-Forming Motifs Promote Mfd-Dependent Stationary-Phase Mutagenesis in Bacillus subtilis. Microorganisms, 9(6), 1284. https://doi.org/10.3390/microorganisms9061284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop