Skip to main content
Advertisement
  • Loading metrics

Co-Orientation of Replication and Transcription Preserves Genome Integrity

  • Anjana Srivatsan,

    Affiliation Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America

  • Ashley Tehranchi,

    Affiliation Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America

  • David M. MacAlpine,

    Affiliation Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, North Carolina, United States of America

  • Jue D. Wang

    jdwang@bcm.tmc.edu

    Affiliation Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America

Abstract

In many bacteria, there is a genome-wide bias towards co-orientation of replication and transcription, with essential and/or highly-expressed genes further enriched co-directionally. We previously found that reversing this bias in the bacterium Bacillus subtilis slows replication elongation, and we proposed that this effect contributes to the evolutionary pressure selecting the transcription-replication co-orientation bias. This selection might have been based purely on selection for speedy replication; alternatively, the slowed replication might actually represent an average of individual replication-disruption events, each of which is counter-selected independently because genome integrity is selected. To differentiate these possibilities and define the precise forces driving this aspect of genome organization, we generated new strains with inversions either over ∼1/4 of the chromosome or at ribosomal RNA (rRNA) operons. Applying mathematical analysis to genomic microarray snapshots, we found that replication rates vary dramatically within the inverted genome. Replication is moderately impeded throughout the inverted region, which results in a small but significant competitive disadvantage in minimal medium. Importantly, replication is strongly obstructed at inverted rRNA loci in rich medium. This obstruction results in disruption of DNA replication, activation of DNA damage responses, loss of genome integrity, and cell death. Our results strongly suggest that preservation of genome integrity drives the evolution of co-orientation of replication and transcription, a conserved feature of genome organization.

Author Summary

An important feature of genome organization is that transcription and replication are selectively co-oriented. This feature helps to avoid conflicts between head-on replication and transcription. The precise consequences of the conflict and how it affects genome organization remain to be understood. We previously found that reversing the transcription bias slows replication in the Bacillus subtilis genome. Here we engineered new inversions to avoid changes in other aspects of genome organization. We found that the reversed transcription bias is sufficient to decrease replication speed, and it results in lowered fitness of the inversion strains and a competitive disadvantage relative to wild-type cells in minimal medium. Further, by analyzing genomic copy-number snapshots to obtain replication speed as a function of genome position, we found that inversion of the strongly-transcribed rRNA genes obstructs replication during growth in rich medium. This confers a strong growth disadvantage to cells in rich medium, turns on DNA damage responses, and leads to cell death in a subpopulation of cells, while the surviving cells are more sensitive to genotoxic agents. Our results strongly support the hypothesis that evolution has favored co-orientation of transcription with replication, mainly to avoid these effects.

Introduction

The fundamental processes of replication and transcription take place on the same template efficiently and accurately, requiring them to be coordinated with each other to avoid potential conflicts. In cells growing rapidly, both replication and transcription of ribosomal RNA (rRNA) genes, and many other genes, are initiated more frequently, further elevating this potential conflict [1][4]. Due to the asymmetry of the replisome and the transcription complex, the outcome of their encounter should depend strongly on their relative directionality. RNA polymerase (RNAP) is dislodged by replication in either direction [5],[6]. On the other hand, replication is affected mostly by head-on transcription [6][12].

Preventing or resolving this conflict not only requires numerous protein factors [13][16] but may also underlie several non-random aspects of genome organization [17],[18]. First, the highly-expressed rRNA and tRNA genes are transcribed almost exclusively co-directionally with replication across numerous species [19],[20]. Chromosomes of the bacteria Bacillus subtilis and Escherichia coli are replicated by bi-directional replication forks initiated from a single origin (oriC), and all rRNA operons are oriented away from oriC [21][25]. In yeast, replication fork barriers at the end of ribosomal DNA operons prevent replication from entering head-on into these strongly-transcribed regions [26]. Second, other highly-transcribed genes are also significantly enriched in the leading strand of replication in bacteria, ensuring that their transcription is co-oriented with replication [27]. This feature may be conserved in certain regions of the human genome [28]. Third, longer transcription units are enriched in the leading strand [27],[29]. Fourth, essential genes are enriched to a greater extent than non-essential genes in the leading strand [19]. Finally, there is a general bias for co-directionality of replication and transcription. In B. subtilis and E. coli, this bias is 75% and 55% of all genes, respectively [22][24].

Despite a general theme of avoiding head-on transcription and replication, the precise evolutionary forces shaping these inter-connected aspects of genome organization are not understood. The effect of head-on replication on transcription is proposed to impact fitness negatively by interrupting the expression of highly-transcribed genes [27], or in the case of essential genes, by leading to the formation of incomplete transcripts, which subsequently results in toxic truncated polypeptides [19]. However, the effects on replication are also deleterious. In E. coli, replication rate is largely unaffected by co-directional transcription, but is significantly slowed when it occurs head-on to a strong transcription unit [5],[30]. In addition, reversing transcription bias over an extended segment of the B. subtilis genome leads to a significant (30%) decrease of replication rate, extending the time required to replicate the chromosome and potentially impeding the cell cycle [31]. Head-on orientation of replication and transcription has been shown to result in genome instability, which can be due to obstructed replication or disrupted transcription [32][35]. It is proposed that the transcription of essential genes is preferentially co-oriented to lower their rate of mutagenesis [30]. Finally, apart from effects on replication and transcription, the transcription bias is also proposed to promote chromosome segregation [36],[37]. Is there a single evolutionary advantage associated with the co-directional bias? Alternatively, is the orientation of each gene selected in its own right? One challenge in understanding the evolutionary bases of orientation biases is dissecting how different aspects of genome organization are important in different circumstances and how they impact cellular fitness.

Here we report that the extent of the impact of head-on transcription on replication differs between genes within the same organism B. subtilis. This was dissected by creating new inversions of either an extensive region of the genome, or a localized region containing strongly-transcribed rRNA genes. Using quantitative genomic approaches, we observed differential rates of replication throughout the genome of the inversion strains—normal replication in intact genomic positions, impedance of replication elongation by ∼30% within the head-on region, and strong blockage of replication at inverted rRNA operons. We further characterized the fitness cost and found that inversion of the oriC-proximal half of a replichore results in a small decrease in growth rate in minimal medium, but is sufficient to confer a significant competitive disadvantage. On the other hand, the replication block at rRNA operons leads to major disruption of replication, induction of the DNA damage response and cell death. We also observed that the rate of mutation of the gene rpoB is increased when it is transcribed head-on to replication within an extended chromosomal inversion, specifically in rich medium. Our results strongly suggest that preservation of genome integrity has contributed to evolution of the genome-wide co-directional bias and its further enrichment in highly-expressed and essential genes.

Results

Inversion of ∼1/4 of the Bacillus subtilis chromosome to reverse its transcription bias

We previously moved the origin of replication (oriC) away from its endogenous position at 0° (Figure 1A) to 257° (not shown) or 94° (Figure 1B) to reverse the genomic transcription bias in an extended region of the chromosome, and observed that replication elongation was slowed moderately between 0° and the ectopic oriC position due to transcription [31]. This raised the intriguing question: what would be the potential impact of reversed transcription bias on cellular fitness and genome integrity? However, this question cannot be answered using these strains, because other aspects of their genome organization were also altered, including location of oriC and symmetry of the replichores (one spanning ¾ of the chromosome, the other ¼ of the chromosome). Such alterations have been shown to strongly impact cellular fitness in both E. coli and B. subtilis [38][40].

thumbnail
Figure 1. Alteration of different aspects of genome organization results in different growth defects.

(A) Schematic diagram of the wild type B. subtilis chromosome (black circle). oriC at 0° and terC at 172° represent the origin and terminus of DNA replication, respectively. Orange arrows: replication; grey arrows: predominant direction of transcription; blue arrowheads: rRNA operons. (B) Schematic diagram of a mutant chromosome with oriC relocated to 94° [31], resulting in head-on transcription and unequal replichores (HT+UR). Green arrow: replication head-on to transcription between 0° and 94°. (C) Schematic diagrams of homologous recombination events leading to inversions of ¼ of the chromosome. L and R: sequences flanking oriC that are repeated at 0° and 94°; spc: spectinomycin resistance gene, which was inserted to replace the endogenous oriC at 0° [31]. Recombination at L (dashed lines) and R (dotted lines) give rise to the chromosomes shown schematically in (D,E), respectively. (D) Schematic diagram of the chromosome with an inversion between 0° and 94° resulting in head-on transcription (HT). (E) Schematic diagram of the chromosome with an inversion between 0° and 94° resulting in unequal replichores (UR). (F) Doubling times at 37°C in liquid LB (red bars), minimal medium with (CAA, grey bars) or without (Min, black bars) casamino acids, calculated by measuring OD600. Data shown are for strains in the JH642 background.

https://doi.org/10.1371/journal.pgen.1000810.g001

To examine exclusively the biological impact of reversed transcription bias, we constructed several new strains. We took advantage of the fact that the strain with repositioned oriC (Figure 1B) has oriC-flanking sequences present both at 0° and the ectopic location (Figure 1C). We reasoned that homologous recombination might occur at these repeats, and screened for such progenies (Figure 1D and 1E). Homologous recombination of repeats upstream of oriC (Figure 1C- repeats marked L, ∼400 bp) repositioned oriC to 0°, with concurrent inversion of the 0°–94° portion of the chromosome. The resulting strain had head-on transcription (HT) between 0° and 94° and equal replichore lengths (Figure 1D). Using the same strategy, we also obtained strains in which homologous recombination had taken place between repeats downstream of oriC (Figure 1C- repeats marked R, ∼300 bp). The resulting chromosomes have oriC positioned at 94° and unequal replichores (UR) but without extended regions of head-on transcription (Figure 1E). Finally, to minimize the possibility of reversion of the inversion, we removed a portion of the remnant homology region at 94°.

Using the HT strain, we evaluated the impact of head-on transcription on fitness by first examining its exponential growth (Figure 1F). The HT strain grew slowly in rich medium (LB) with a doubling time of 70 minutes, compared to 20 minutes for the control. In minimal medium however, the doubling time of the HT strain was similar to that of the control (53 and 52 minutes, respectively). In contrast, the UR strain was sicker than the control in both LB and minimal media (doubling times of 52 and 77 minutes, respectively), indicating a general growth defect due to differing replichore lengths and/or ectopic positioning of oriC. Therefore, unlike the general growth defect introduced by uneven replichores, the growth defect caused by inverting transcription bias over ¼ of the chromosome is nutrient-dependent.

Head-on transcription decreases replication fork speed on a genomic scale

We next examined whether head-on transcription has an effect on replication in the HT strain. We monitored synchronized fork progression in this strain using genomic microarrays (Figure 2A), in minimal medium where no significant loss of growth rate was observed (Figure 1F). Cells were synchronized for their replication cycles using a temperature-sensitive allele of the replication protein DnaB [41],[42]. The gene dosage profile obtained 30 minutes after replication initiation indicates that ∼50% of cells initiated replication. The average position of replication forks can be estimated as the midpoint of the transition between replicated and unreplicated genomic positions [43]. This position is ∼0.68 Mbp from oriC on the left replichore and ∼0.56 Mbp on the inverted right replichore (Figure 2A, blue arrows), indicating that replication forks move slower within the inverted region. We inhibited transcription initiation by adding the drug rifampicin 4 minutes after synchronized replication began, and found that in rifampicin-treated cells replication forks progressed further in the inverted region (∼0.72 Mbp) (Figure 2A, red arrows) compared to untreated cells (∼0.56 Mbp), demonstrating that the reduction in replication fork speed is due to transcription. This reduction in fork movement does not lead to proportionally slower growth likely because B. subtilis has flexible cell division cycles and can compensate for slower fork progression via multifork replication [4].

thumbnail
Figure 2. Genomic microarray profiles of the HT inversion strain.

(A) Synchronized replication profile of the HT strain. Cells containing the dnaB134ts allele were grown in Min at 30°C and their replication cycle was synchronized by shifting to 45°C for 60 minutes, and then back to 30°C to allow replication initiation. Genomic profiles were obtained 30 minutes after initiation, relative to pre-initiation reference DNA. Profiles were obtained without (blue) or with (red) addition of rifampicin (rif) 4 minutes after initiation to inhibit transcription. Elevation of the baseline of the gene dosage profile on the right replichore is likely due to incomplete replication of the pre-initiation reference DNA of the HT strain, despite the 60 minutes incubation at 45°C. Grey shaded region: inversion; blue and red arrows: average positions of replication forks in the absence (blue) or presence (red) of rif. Insets: schematic diagram of chromosomes of the HT strain undergoing synchronized replication with (top), and without (bottom) transcription (grey arrows). Purple ovals: replisomes. (B) Overlay of the asynchronous genomic profiles of the HT strain grown in Min (blue), and CAA (orange). Profiles were obtained from asynchronous cultures grown at 37°C to OD600∼0.5. Average gene dosage ratios (log2) are plotted relative to gene positions adjusted according to deletion of the phage SPβ and the integrative and conjugative element ICEBs1 clusters. The prophage-like skin element is not removed in this strain, which likely explains the small peak at −1.5 Mbp (black arrow). Green arrows: sharp changes in slope at rRNA loci in the CAA profile. Inset: schematic diagram of chromosomes of the HT strain undergoing asynchronous replication. (C) Expanded view of the ratio of gene dosage (log2) in CAA versus Min within the region outlined in B (−1.2 Mbp to 1.2 Mbp). Green arrows: positions of inverted rRNA operons; orange arrow: position of co-directional rrnB; grey shaded region: inversion.

https://doi.org/10.1371/journal.pgen.1000810.g002

Nutrient-dependent changes in the asynchronous gene dosage profile

We next examined whether the decreased replication rate in the inverted region varied depending on nutrient status and genomic position. To this end, we obtained the genomic microarray profile of the HT strain during exponential growth (Figure 2B). Cells are not synchronized and hence the positions of the replication forks would vary from cell to cell (Figure 2B inset). Importantly, by assuming that these cells are in a steady state and their genomic profile is time-invariant on a population basis, we can use this profile to calculate replication speed at every position on the chromosome. This speed is inversely correlated with the local slope of log values of the gene dosage with respect to gene positions (see Materials and Methods).

We first observed that in minimal medium, the genomic profile was smooth but asymmetric (Figure 2B, blue). The rates of replication were similar in the unaltered regions of the chromosome, indicated by similar slopes on the left replichore (172° to 360°, 0.408±0.002/Mbp) and the non-inverted region on the right replichore (94° to 153°, 0.462±0.007/Mbp) (Table 1). In contrast, within the inverted segment on the right replichore (0° to 94°), the slope was 0.646±0.004/Mbp, indicating a ∼30% decrease in fork speed within the head-on region, in agreement with our previous results using an ectopic oriC [31].

thumbnail
Table 1. Slopes of asynchronous gene dosage profiles in the HT inversion strain.

https://doi.org/10.1371/journal.pgen.1000810.t001

Interestingly, when cells were grown in a relatively rich medium (minimal medium supplemented with casamino acids, hereafter referred to as CAA), the gene dosage profile changed sharply at specific locations on the chromosome (Figure 2B, orange). These transition points could be differentiated more clearly when we derived the relative gene dosages between the profiles in CAA vs. minimal medium (Figure 2C), and clearly corresponded to the positions of the rRNA operons within the inverted segment (Figure 2B and 2C, green arrows). The steep slopes at these transitions suggest that replication progression is strongly impeded and even stalled at these locations. Other than at rRNA loci, the genomic profile in CAA was very similar to that of minimal medium including at the intact rRNA operon on the left replichore (Figure 2C, orange arrow). Our results indicate that the genome-wide impedance of replication by reversed transcription bias other than at rRNA loci is largely uniform and unaltered by nutrient conditions. Importantly, we identify rRNA loci as positions of strong, nutrient-dependent obstruction of fork movement when they are transcribed head-on.

Head-on transcription of rRNA disrupts replication and induces the DNA damage response

We next examined whether the forks obstructed by head-on transcription during growth in CAA were also disrupted. The recombination protein RecA localizes to stalled replication forks as foci only when the forks are disrupted [44]. Hence, we examined the sub-cellular localization of RecA using a recA-gfp fusion construct [45]. We visualized microscopically cells carrying this allele grown in CAA, and found that 97% of the cells of the HT strain had RecA foci compared to 27% in the control (Figure 3A and 3B), indicating that replication forks are disrupted by head-on transcription.

thumbnail
Figure 3. The HT strain exhibits disruption of DNA replication and loss of genome integrity.

(A) RecA localization and nucleoid morphology in the HT and isogenic control strains in CAA. Upper and lower panels: phase contrast images (blue) overlaid with RecA-GFP (green) or DAPI (white) fluorescence images, respectively. (B) Percentages of cells with RecA-GFP foci in the indicated media. Cells were grown in CAA or Min and stained with DAPI and the membrane dye FM4-64 (not shown) to count individual cells. Light grey bars: control; dark grey bars: HT strain. The HT strain with recA-gfp was extremely sick in LB. (C) Induction of the SOS DNA damage response in the HT and isogenic control strains in the indicated media, as monitored using a TagC-GFP reporter. The total length of SOS-positive cells was divided by the total length of cells in each field. Light grey bars: control; dark grey bars: HT strain. (D) Live/dead staining of the HT and control strains in LB and Min. The fluorescent dyes SYTO9 (cyan) and propidium iodide (red) were added to exponential phase cultures (OD600∼0.2–0.6) to label live and dead cells, respectively. (E) Average percentages of dead cells in the control (light grey bars) and HT (dark grey bars) strains in different growth media (LB, CAA and Min). The cell length that was stained with propidium iodide was divided by the total length of cells in each field. (F) Mitomycin C (MMC) sensitivity of the HT and isogenic control strains in LB and Min. Cells were grown to OD600 = 0.3 and serial dilutions ranging from 10−2 to 10−6 were spotted on LB and Min plates with or without MMC (0.0625 µg/ml). Plates were scanned after overnight incubation at 37°C. Scale bar in A and D: 5 µm.

https://doi.org/10.1371/journal.pgen.1000810.g003

We also noticed that the HT strain had abnormal nucleoid morphology, which appeared filamented and even fragmented in some cases (Figure 3A, lower panels). Cell lengths also significantly increased (not shown) and chain lengths doubled (Table 2), which might explain the ∼2-fold decrease in the number of colony-forming units in strain HT compared to its control (Table 3).

thumbnail
Table 2. Increased length of chains of cells in rich medium.

https://doi.org/10.1371/journal.pgen.1000810.t002

thumbnail
Table 3. Colony-formation by inversion strains in different media.

https://doi.org/10.1371/journal.pgen.1000810.t003

In addition to RecA foci formation, a subpopulation of cells exhibits the SOS DNA damage response (Figure 3C). Using a GFP-fusion reporter of tagC, a member of the SOS regulon [46], we found that the SOS response is induced in greater than 6% of single cells of the HT strain in LB medium, but in less than 1% in minimal medium (Figure 3C). The increase in SOS response in rich medium is accompanied by increased cell death. We performed live/dead staining in which the nucleoids of dead cells with permeable membranes stain with propidium iodide (red), while live cells stain with SYTO9 (cyan) (Figure 3D). There was a marked increase in the fraction of dead cells in the HT strain relative to the isogenic control, again specifically in rich medium (12% versus 1% of cells in LB) (Figure 3E). This suggests that failure to repair replication forks disrupted by strong head-on transcription might lead to failure to complete replication and cell death. In agreement with this hypothesis, we found that HT cells show higher sensitivity to the genotoxic agent mitomycin C in rich medium (Figure 3F), suggesting that exogenous DNA damage adds further demand on their already overwhelmed DNA repair capacity, and dramatically elevates cell death in the population.

Increased mutation rate of a gene transcribed head-on to replication

Having obtained evidence of disruption of replication by head-on transcription, we next examined whether it also has a consequence on genome stability by measuring the rates of mutations conferring resistance to the drug rifampicin (rifR). In B. subtilis, rifR mutations map to rpoB [47],[48], which is transcribed co-directionally to replication in the wild-type strain but head-on in the HT strain (Figure 4A). Since rpoB encodes a sub-unit of RNA polymerase, its mutation might confer a growth advantage in this scenario. Hence we analyzed the results of the fluctuation test with the P0 method, which measures the rate of mutation independently of its effect on growth rate [49]. We observed that the mutation rate increased ∼3-fold in rich medium in the HT strain compared to an isogenic control with no inversion (Figure 4C). This could be due to a global cellular response to disruption of replication caused by the chromosomal inversion. To examine this possibility, we also monitored the rifR mutation rate in a strain with inversion of ½ of the left replichore, leaving the rpoB region unaltered (Figure 4B). This strain has the same extent of reversed transcription bias as the HT strain, exhibits a similar nutrient-dependent growth defect (Table 4, Figure S1A), and has one inverted rRNA operon located near oriC after chromosome inversion, which causes replication blockage in rich medium, as monitored by microarrays (Figure S1B and S1C). However, there was no increase in the rifR mutation rate in this strain (Figure 4C). Therefore the presence of an inversion alone is not sufficient to cause an increase in rifR mutation rate, rather it is specific to the strain in which rpoB is within the inverted region. In minimal medium all strains had similar rifR mutation rates (Figure 4D).

thumbnail
Figure 4. Increased mutation rate of a gene within the inversion in the HT strain.

(A,B) Schematic diagrams of chromosomes (black circles) with inversion of ½ of either the right (A) or the left (B) replichore, resulting in head-on transcription (HT) within the inversion. Orange and green arrows: replication co-directional and head-on to transcription, respectively; grey arrows: predominant direction of transcription. Mutations conferring resistance to rifampicin (rifR) map to the rpoB gene (red arrow) which is transcribed co-directional to replication in the wild type and left replichore HT strains, but transcribed head-on in the right replichore HT strain. (C) Spontaneous rifR mutation rates as measured by fluctuation tests in LB (dark grey bars) in the control, right and left replichore HT strains. *: rifR rate in the right replichore inversion strain is significantly different from control (P<0.05). (D) RifR mutation rates in Min (light grey bars).

https://doi.org/10.1371/journal.pgen.1000810.g004

Inversion of ribosomal RNA operons is sufficient to disrupt replication

The most dramatic reduction of replication speed due to head-on transcription occurs at the rRNA operons within the inversion (Figure 2C). This suggests that the nutrient-dependent effect on replication fork progression is mostly due to inversion of the strongly-transcribed rRNA operons. We tested this hypothesis by examining the consequences of specifically inverting rRNA operons. We inverted the rrnIHG cluster (∼17 kbp) that contains 3 rRNA operons and 6 tRNA genes, by inserting two overlapping halves of the neomycin resistance gene (neo) flanking the cluster, and selecting for recombination events that created a complete neo gene, similar to [50] (Figure 5A and 5B). The rrn inversion strain was inviable in LB and had a strong growth defect in CAA compared to the pre-inversion control (doubling times of 44 and 28 minutes, respectively), while their doubling times in minimal medium were similar (44 and 42 minutes, respectively) (Figure 5C). These results indicate that the growth defect of the HT strain was mostly due to inverted rRNA genes.

thumbnail
Figure 5. Inversion of rrnIHG impedes replication, triggers RecA recruitment, and elevates cell death.

(A,B) Schematic diagrams of the chromosomes (black circles) of the pre-inversion control strain (A) and the rrnIHG (rrn) inversion strain (B). Blue arrowheads: rRNA operons; large green arrowhead: rrnIHG cluster. (C) Doubling times of the pre-inversion and rrn inversion strains at 37°C in CAA (grey bars) and Min (black bars). (D,E) Asynchronous gene dosage profiles of the rrn inversion strain in CAA (D) and Min (E). Profiles were obtained from asynchronous cultures grown at 37°C in the indicated media to OD600∼0.5. Average ratios of copy number (log2) in the rrn inversion strain relative to fully-replicated control are plotted against the genomic position. Arrows indicate the position of the rrnIHG inversion. (F) RecA recruitment and nucleoid morphology in the pre-inversion and rrn inversion strains in CAA, obtained as described in Figure 3A. Blue: phase contrast images; green: RecA foci; white: nucleoids stained with DAPI. (G) Average percentage of cells with RecA recruitment in the indicated media. Light grey bars: pre-inversion control; dark grey bars: rrn inversion. (H) Live/dead staining of the pre-inversion and rrn inversion strains grown in CAA and Min. Images were obtained as described in Figure 3D. Cyan: live cells; red: dead cells. (I) Average percentages of dead cells in the pre-inversion (light grey bars) and rrn inversion (dark grey bars) strains. The cell length that was stained with propidium iodide was divided by the total length of cells in each field. Scale bar in F and H: 5 µm.

https://doi.org/10.1371/journal.pgen.1000810.g005

We next examined asynchronous replication in the rrn inversion strain using genomic microarrays, and observed impedance specifically at the inverted loci, where the gene dosage profile showed sharp discontinuity (Figure 5D and 5E). This indicates that a significant number of replication forks were stalled in this short segment. This effect was much stronger in cells replicating asynchronously in CAA than in minimal medium. Further, the majority of cells with the rrn inversion had RecA-GFP foci/filaments when grown in CAA (Figure 5F and 5G), indicating that strong head-on transcription at this cluster is sufficient to disrupt replication forks. Notably, the nucleoid morphology of the rrn inversion strain was largely normal (Figure 5F, lower panels), indicating that the gross nucleoid defects of the HT strain were not due to disruption of replication by rrn inversion.

Finally, the rrn inversion lowered cell viability. The rrn inversion strain had a much higher fraction of dead cells relative to the pre-inversion control in CAA (10% versus 0.7%, respectively) (Figure 5H and 5I). Thus strong head-on transcription of just rRNA and tRNA operons drastically impacts cell viability through disruption of replication.

Discussion

Genome organization has evolved to enhance fitness, as evidenced by observations that certain genome rearrangements are not tolerated, or cause growth defects [17], [38], [39], [51][56], while others do not and might even be prevalent [57]. One feature of genome organization is that it precludes extensive head-on transcription [5],[17],[18],[27],[31]. In this study we engineered B. subtilis strains with the co-directional bias of replication and transcription reversed over either an extended segment of the genome or at a localized rRNA gene cluster. We employed microarray-based copy number profiling that enabled us to visualize directly the replication status of the entire genome, to identify positions at which there are significant perturbations and to quantify the extent of such perturbations. We found that replication is affected by head-on transcription in at least two ways: one is the apparently uniform deceleration throughout the extended region of reversed transcription bias, and the second and stronger is the disruption of replication at highly-transcribed rRNA loci. Disruption of replication at rRNA genes activates DNA repair pathways and results in sensitivity to genotoxic stress and loss of viability. Together these observations support the hypothesis that head-on collisions between transcription and replication result in loss of genome integrity, and that avoidance of this consequence contributes to the evolution of co-orientation bias in genomes.

Other aspects of genome organization

Previously we had engineered B. subtilis strains in which the replichores were unequally distributed, and a significant portion of the genome was replicated by forks traveling in the opposite direction to transcription. We found that DNA replication elongation was impeded within the region of reversed transcription bias and these strains had strong growth defects [31]. However it was not possible to attribute the growth defects to reversed transcription bias alone since uneven replichores themselves have been shown in E. coli to lead to strong growth defects and dependence on recombination and/or DNA translocation machineries for viability [38],[39]. Therefore, we have constructed several new strains that separate the alteration of replichore symmetry from that of the transcription-replication bias (Figure 1).

We found that the B. subtilis HT strain with head-on transcription and replication exhibits a strong growth defect only in rich medium, and the UR strain with asymmetric replichores exhibits growth defects regardless of growth medium. The effect of both aspects of genome alteration on fitness is approximately multiplicative (Figure 1F and Table S1) [58],[59]. There are small deviations but this is significant only in minimal medium (p = 0.01). Thus we were able to largely separate two aspects of genome organization–the impact of unequal replichore size, and the impact of colliding transcription and replication forks, on genome integrity and cellular fitness.

We found that the HT strain with the extended inversion also exhibits strong disruption of replication and induction of DNA damage response in rich medium. In addition, it has altered nucleoid morphology and a long and twisted cell shape especially in rich medium (Figure 3A and 3D). We further inverted only an rrn cluster, and demonstrated that it is sufficient to cause disruption of replication, but not the nucleoid morphology defect. The change in morphology in the HT strain can be either associated with head-on transcription at locations other than rRNA operons, or with additional effects of inverting ¼ of the chromosome. The latter could include the alteration of gene positions relative to the origin, which affects their dosage [4],[60], ectopic localization of the parS sites which affect chromosome organization [61][63], or defects in chromosome segregation which is proposed to be facilitated by transcription [36],[37].

Relationship between genome-wide replication rate and transcription

Using synchronized microarrays, we demonstrate that inverting the transcription bias over ¼ of the chromosome decreases fork progression rate in a transcription-dependent manner (Figure 2A), confirming our previous results obtained using an engineered strain with an ectopic oriC [31]. However, transcription throughout the genome is not uniform. To examine whether different inverted transcription units have different impacts on replication, we obtained asynchronous microarray profiles of exponentially growing cells (Figure 2B). Using mathematical analysis to obtain the rate of replication as a function of genomic position, we confirmed a modest genome-wide impedance of replication throughout the inverted region and showed that it is mostly independent of gene position. An important exception is at rRNA and tRNA operons, evidenced by the punctuated pattern of replication stalling at these clusters. Further, this stalling is strongly potentiated by growth in rich medium.

In general, growth in rich medium results in higher initiation frequencies of both replication and rRNA transcription [1][4], either of which could elevate the conflict between transcription and replication, thereby accounting for the observed increase in replication stalling. However, we observed that in the right replichore HT inversion strain, replication is not initiated more frequently in rich medium than in minimal medium. In contrast to wild type cells where gene dosage at oriC is much higher in rich medium than in poor medium, indicating higher rate of replication initiation [31], gene dosage at oriC in the HT strain is similar in CAA and minimal media (Figure 2B). The genomic profiles of the HT strain in the two media are almost identical except at the inverted rrn loci (Figure 2C). It is possible that failure of replication elongation prevents subsequent replication initiation; alternatively, replication initiation frequency could be lower because it is coupled to growth [3],[4], which is slower for the HT strain in rich medium (Figure 1F). Regardless of the reason, it is clear that in the inversion strains, the effect on replication observed at rRNA operons in rich medium is not due to increased replication but is exclusively due to stronger rRNA transcription, which is initiated more frequently because of higher iNTP and lower (p)ppGpp levels [64].

Several models exist to explain why replication is stalled by strong head-on transcription of rRNA operons. The replisome might be capable of bypassing a single head-on RNAP, but the presence of multiple RNAPs on the long and highly-transcribed rrn region could make it harder for the replisome to proceed. In addition, since rRNA operons are highly structured regions, their transcription might obstruct replication forks, as proposed for other unusually structured regions [8],[65]. RNAPs might also stall upon head-on replication to form backed-up RNAP arrays. Backed-up RNAP arrays can create a barrier to replication [15]. Finally, head-on transcription might create RNA-DNA duplexes or supercoiling of DNA that poses a barrier to replication [66][68], and this barrier might strengthen to the extent of blocking replication when transcription is sufficiently strong.

Impact of head-on transcription on fitness and genome integrity

Inverting the transcriptional bias of 1/4 of the B. subtilis chromosome slows replication rate within this region by ∼30%. The growth rate of the HT inversion strain in minimal medium is not significantly affected (Table 4). However we discovered that the HT strain indeed has a significant growth disadvantage even in minimal medium when competing with wild type cells, with its relative fitness being 0.92 (+/−0.07) (after factoring in the marker effect) (Table S2, Text S1). This selective effect enables the wild type strain to take over after multiple generations and clearly is sufficient to shape genome evolution.

The impact of inversions on replication and cellular fitness is much stronger when cells are grown in rich media. Inverted rRNA genes in both the HT and rrn strains result in replication blocks (Figure 2 and Figure 5) and likely lead to extensively delayed cell-cycle progression, which explains the dramatic increase of doubling time in rich medium (Figure 1F and Figure 5C). Indeed, blockage of replication elongation has been shown to prevent cell proliferation in E. coli, which can only be reversed upon removal of the barrier [69]. More importantly, we obtained strong evidence that the obstruction created by inverted rRNA transcription also leads to disruption of replication. First, RecA forms foci/filaments in the majority of these cells in rich medium, indicating generation of single-stranded DNA or double-stranded ends (DSEs). Second, there is a significant increase in induction of the SOS DNA damage response [70] in the inversion strains in rich medium (Figure 3C). In B. subtilis the SOS response is not robustly turned on by DSEs due to efficient repair by RecN, and our observation of ∼6% of cells of the HT strain showing SOS induction agrees with the reported value [71]. Third, an increased number of cell deaths occur in the inversion strains especially in rich medium, likely due to failure to repair damaged replication forks. Finally, the inversion renders cells more sensitive to the genotoxic agent mitomycin C especially in rich medium, suggesting that the DNA repair capacity in these cells is highly compromised due to overwhelming demand, leading to detrimental consequences upon challenge by external DNA damage.

It remains unclear whether replication fork collapse [72] takes place soon after forks are stalled by head-on transcription, or only when a second round of replication forks collides with a prior round of stalled replication forks, as was demonstrated previously at the replication terminator sequence in E. coli [73]. It is possible that both types of collisions would lead to disruption of replication, with collision of subsequent replication forks being more costly. This might explain why inverting rrnIHG near the origin of replication results in a particularly strong growth defect, as the second fork encounters the stalled first fork soon after it initiates from the origin of replication. Indeed the oriC-proximal region is the chromosomal location of highly-expressed genes, especially those involved in macromolecular synthesis [74], and thus inversions here might be particularly detrimental.

Disruption of replication has been shown to lead to higher levels of genome instability especially in the vicinity of the disruption [73],[75],[76]. We observed an increase in rpoB mutation rate when the genomic region encoding rpoB is inverted (Figure 4). This increase is only observed when cells are grown in rich medium, implicating a dependence on the level of rpoB transcription or disruption of replication. It is unlikely to be solely due to a global cellular response to disruption of replication such as the SOS response, since inversion of a symmetric half of the other replichore without affecting rpoB, did not elevate its mutation rate despite the fact that replication is clearly inhibited by inverted rrnB transcription in this strain (Figure 4 and Figure S1). One possibility is that disruption of replication at rpoB results in recruitment of an error-prone DNA polymerase via RecA. RecA has been shown to activate directly an error-prone DNA polymerase in E. coli [77]. Further work will be required to differentiate this possibility from other remaining possibilities, e.g., that the increased mutation rate is due to altering rpoB location rather than orientation, or due to disruption of replication at the neighboring rRNA loci, or that the rpoB mutagenesis is due to a threshold level of the SOS response that is only met upon inversion of the right, but not the left replichore.

The impairment of replication by transcription has been shown to result in transcription- associated recombination [32] or deletion [34]. Replication orientation significantly influences the spectrum of point mutations in yeast [35]. This suggests that impairment of genome replication might also contribute to transcription-associated mutagenesis [78][81]. Another major mechanism which might be at play is the direct activation of the error-prone translesion polymerase via its interactions with transcription factors such as NusA [82].

Implications for the evolution of genome organization

Our observations offer strong support for the hypothesis that the effect of transcription on replication is an important driving force in the evolution of genome organization [18]. In addition, our work suggests that the precise cost might vary depending on both the gene type and the growth environment. First, inversion of the strand bias of transcription within an extended segment of the chromosome results in a small growth defect in nutrient-poor medium yet is sufficient to confer a strong competitive disadvantage. Second, the inversion of rRNA operons leads to disruption of DNA replication, which is especially costly if cells are grown in rich medium (Figure 2, Figure 3, and Figure 5). This explains why rRNA operons are all oriented co-directionally. Although the most prominent disruptive effects we observed were at inverted rRNA operons, it is possible that these effects can also be extended to other highly-expressed or long and structured genes. Third, we also observed that mutation rate is higher for the rpoB gene transcribed head-on, supporting a model that co-orientation of transcription and replication of essential genes might have evolved to avoid their mutagenesis [30] (Figure 4). In addition, highly-expressed non-essential genes are also intolerant of mutagenesis since it could result in reduction of their expressivity. Thus minimizing mutagenesis may also underlie the orientation bias of highly-expressed genes.

There are considerable differences in transcription orientation biases among organisms. While low G+C Firmicutes (such as B. subtilis) and the Mycoplasmas have strong transcription orientation biases, other bacteria do not [18]. The widely-studied bacterium E. coli has only 55% co-orientation bias. Interestingly, inversion of long segments of the E. coli replichores including rRNA operons results in no growth defect [38]. What causes such a drastic difference in the penalty of head-on transcription? There are at least three possible explanations: differences in the composition of their replication machineries [24], non-replicative helicases [16] or transcription factors (AT and JDW unpublished). These considerable differences might underlie the differential abilities of organisms to cope with conflict between transcription and replication and thus influence the evolution of their genome organization.

Materials and Methods

Media and growth conditions

Cells were grown in LB or defined minimal medium (50mM MOPS) [83] with 1% glucose, supplemented with 40 µg/ml tryptophan, 40 µg/ml methionine, 40 µg/ml phenylalanine, 100 µg/ml arginine, or casamino acids (CAA) (Difco) (0.5%). Cells containing inversion of rRNA operons were maintained on 5 µg/ml kanamycin (kan). Chloramphenicol (cm), erythromycin (erm), and spectinomycin (spc) were used as described at 5, 0.5 and 40 µg/ml, respectively.

Strains and plasmids

Standard techniques were used for genetic and molecular biological manipulations [84]. Strains used are listed in Table S3. Primer sequences are listed in Table S4. Strains were constructed in the JH642 background [85], and in the YB886 phage-defective background [86] because phage excision and duplication during stressful conditions often create localized gene copy number alterations. YB886 is cured of phage SPβ, defective for phage PBSX induction [86], and also lacks the transposon-like element ICEBs1 [87].

Strain JDW704 was created by introducing an ectopic oriC at aprE (94°) in YB886 and then deleting the endogenous oriC (0°), using genomic DNA from JDW258 [31] and MMB703 [40], respectively. Progenies of JDW704 were screened for recombination events in sequences flanking oriC by PCR using oJW114/oJW135, oJW115/oJW157, oJW112/oJW75 and oJW113/oJW146. The stabilized inversion strain JDW713 was obtained by transforming JDW704 with linearized plasmid pJW247. The isogenic control was generated by transforming YB886 with linearized plasmid pJW207. The left replichore inversion strain JDW605 was created by screening progenies of MMB703 [40] by PCR using primers oJW112/oJW388, and oJW113/oJW146.

Inversion of rRNA operons was performed similarly to the method described in [50] and [88]. Briefly, two halves of the neomycin resistance gene (neo) overlapping by 583 bp, were inserted flanking the rrnIHG region in strain YB886 using plasmids pJW260 and pJW261, respectively. The strain was plated on kanamycin to select for cells in which a complete neo gene was created by recombination between the two halves. The inversion junctions were tested by PCR using oJW450/oJW442, and oJW452/oJW436. Genomic coordinates [22] of the rrn inversion are 159778 to 176408 in JDW860 and 154793 to 176408 in JDW861.

Plasmid pJW247 was constructed by cloning sequences flanking the inversion junction into the vector pUC18, on either side of the cat gene. The sequences were amplified using oJW206/oJW360 and oJW210/oJW211, and cat was amplified from pGEMcat [89] using oJW208/oJW209. Plasmid pJW207 was constructed similarly, except that the PCR product of oJW204/oJW205 was used instead of oJW206/oJW360. Plasmid pJW260 was constructed by cloning sequences flanking the downstream rrn inversion junction into the vector pUC18erm, on either side of the ‘neo fragment. The sequences were amplified using oJW434/oJW435 and oJW438/439, and ‘neo was amplified from pBEST502 [90] using oJW436/oJW437. Plasmid pJW261 was constructed by cloning sequences flanking the upstream rrn inversion junction into the vector pBEST501 [90], on either side of the neo’ fragment. The sequences were amplified using oJW428/oJW429 and oJW432/433, and cat was amplified from pGEMcat using oJW485/oJW486.

Microarrays and data analysis

Strains with the dnaB134ts allele [41],[42] were grown in minimal medium at 30°C to OD600 = 0.2. Cells were shifted to 45°C for 60 minutes to prevent new initiation and allow ongoing replication to complete. The temperature was rapidly shifted down to 30°C to allow synchronized initiation of replication, and the culture was split into 2 flasks. 4 minutes after the down-shift, rifampicin was added to one flask to 0.25mg/ml. Cells were collected 30 minutes after the down-shift, mixed with an equal volume of 100% ice-cold methanol and processed for microarray analysis as described [44].

Hybridization was performed according to the Agilent Oligo aCGH protocol using custom 44K oligonucleotide Agilent microarrays. Microarrays were scanned using a GenePix 4000B scanner (Axon Instruments). Cy3 and Cy5 levels were quantified using Agilent's Feature Extraction software. Relative DNA content (log2 ratio of Cy3 to Cy5 levels) was plotted against the gene position on the chromosome, with the origin in the center and the terminus at each end of the x-axis. For the inversion strains, the genomic positions are rearranged to reflect genome reorganization. For synchronized microarrays, the rolling average of gene dosage ratios (log2) for every 200 consecutive positions was calculated from the raw data, and plotted against the mid-point of these positions.

Calculating replication rate as a function of gene-position

Genomic microarray profiles of cells grown in mid-exponential phase were obtained by hybridizing against a synchronized reference, such that the ratios were proportional to the actual gene dosage. The data for the right replichore were analyzed based on the following (for the left replichore, the equations are identical except with negative signs):

1. During exponential growth, the total number of cells increases exponentially with cell mass doubling. Define T as the mass doubling time, t as the time of measurement, and N(t) as the total number of cells at time t, then(E1)

2. Define x as the position of the gene, and the total gene dosage of the population of cells at position x as f(x, t). There are 4.2Mbp of nucleotides and over 108 cells, so we can approximate f as a continuous function with continuous variables x and t, despite that fact that each cell is undergoing discrete events including replication initiation and cell division.

Assuming that cells are in steady state and their genomic profile is time-invariant on a population basis, we have(E2)

Combining E1 and E2, we have(E3)

3. The rate of replication fork progression is a function of genome position x but not of time t during steady state growth. Therefore we can define the replication rate at position x as v(x). By this definition, for small Δx, we have:(E4)Combining E3 and E4, we haveDefine , thenwhich can be rearranged as(E5)which can also be written in differential form as(E6)

Once we obtain g(x, t), equations E5 and E6 give a precise definition of how to obtain v(x). Our microarray data show that, for a fixed t, g(x, t) can be approximated by a piecewise linear function over x. Let a discrete series (x1, x2….xm) be the connecting points of the piecewise linear function. If x falls between xi and xi+1, and , where ai is independent of t, then we can obtain v(x) as:(E7)

Measurement of mutation rates

Measurement of rifampicin-resistance (rifR) mutation rates was performed using the fluctuation test as described [49]. 50 parallel cultures of 1–2 ml each were set up for each strain, grown at 37°C to OD600∼0.5 and plated on minimal medium containing rifampicin (5 µg/ml). Serial dilutions were also plated on non-selective medium to count the number of colony-forming units. After incubation at 37°C for 36 hours, the number of plates with no rifR colonies was counted. RifR mutation rate was calculated using the P0 method [49]. Results from 2 or more independent experiments were averaged. Error bars represent the range of data for n = 2, and the standard error for n>2.

Live cell microscopy

The recA::(recA-gfp spc) allele [45] was used to replace the endogenous recA gene in the inversion and control strains. The recA-gfp strains were grown at 30°C in CAA and minimal medium. At OD600 = 0.2–0.6, 250 µl aliquots were labeled with the membrane dye FM4-64 (0.05 µg/ml) and/or DAPI (0.1 µg/ml). Cells were spotted onto thin agarose pads (1% agarose in 1× Spizizen's salts) on multi-well slides, covered with a cover slip and imaged in a Zeiss Axiovert microscope using a 100× oil immersion objective. Images were captured using a Hamamatsu Digital CCD camera, and analyzed using the AxioVision software. The number of cells with RecA-GFP foci or filaments relative to the total number of cells in each image was counted.

SOS induction was monitored similarly in cells containing the tagC-gfp reporter grown at 37°C in LB, CAA and minimal medium. Data were analyzed by calculating SOS induction per cell length.

For estimating the number of dead cells microscopically, the Live/Dead BacLight Bacterial Viability Kit (Molecular Probes) was used, in which live cells are labeled with SYTO9 (green fluorescence, colored in cyan in Figure 3D and Figure 5H) and dead cells with propidium iodide (red fluorescence). Data were analyzed by calculating cell death per cell length.

For all microscopy experiments approximately 1000 cells were counted for each strain and each growth condition. Results from 2 or more independent experiments were averaged. Error bars represent the range of data for n = 2, and the standard error for n>2.

Measurement of number of colony-forming units (cfu)

Cells were grown at 37°C to OD600∼0.2–0.6 and serial dilutions were plated on minimal medium. Colonies were counted after 36 hours of incubation at 37°C. Data from 3 independent experiments were averaged.

Mitomycin C sensitivity assay

Cells were grown in either LB or minimal medium to OD600 = 0.3. 5 µl of 1∶10 serial dilutions ranging from 10−2 to 10−6 were spotted correspondingly on LB, Min and plates supplemented with mitomycin C to a concentration of 0.0625 µg/ml. Plates were incubated overnight at 37°C and photographed the next day. The experiment was repeated twice and representative images are shown.

Supporting Information

Figure S1.

Inversion of the oriC-proximal half of the left replichore (HT-left) also leads to impedance of replication fork progression. (A) HT-Left exhibits a strong growth defect especially in rich media. Doubling times at 37°C in liquid LB (red bars), minimal medium with (CAA, grey bars), or without (Min, black bars) casamino acids were calculated by measuring OD600. (B) Overlay of the asynchronous genomic profiles of the HT strain grown in Min (blue), and CAA (orange). Profiles were obtained from asynchronous cultures grown at 37°C to OD600∼0.5. Average gene dosage ratios (log2) are plotted relative to gene positions adjusted according to known deletions of the background strain JH642 [6]. (C) Ratios of gene dosage (log2) in CAA versus Min of the HT-Left strain, calculated similarly to Figure 2C. Green arrow: position of the inverted rRNA operon rrnB; orange arrows: positions of co-directional rRNA operons; grey shaded region: inversion.

https://doi.org/10.1371/journal.pgen.1000810.s001

(0.95 MB TIF)

Table S1.

Comparison between the observed and expected fitness of the HT and UR mutants in the indicated growth media, under the multiplicative null model.

https://doi.org/10.1371/journal.pgen.1000810.s002

(0.02 MB DOC)

Table S2.

Relative fitness of the HT strain in the indicated growth media (WHT).

https://doi.org/10.1371/journal.pgen.1000810.s003

(0.03 MB DOC)

Text S1.

Supplemental Materials and Methods.

https://doi.org/10.1371/journal.pgen.1000810.s006

(0.03 MB DOC)

Acknowledgments

We thank Allison Kriel for technical assistance; Alycia Bittner, Elicia Grace, Phil Hastings, Xiangwe He, Elizabeth Ostrowski, Susan Rosenberg, and Lizhao Zhang for helpful discussions; and the anonymous reviewers for critical comments.

Author Contributions

Conceived and designed the experiments: AS JDW. Performed the experiments: AS AT. Analyzed the data: AS JDW. Contributed reagents/materials/analysis tools: DMM JDW. Wrote the paper: AS JDW.

References

  1. 1. Miura A, Krueger JH, Itoh S, de Boer HA, Nomura M (1981) Growth-rate-dependent regulation of ribosome synthesis in E. coli: expression of the lacZ and galK genes fused to ribosomal promoters. Cell 25: 773–782.
  2. 2. Gourse RL, Gaal T, Bartlett MS, Appleman JA, Ross W (1996) rRNA transcription and growth rate-dependent regulation of ribosome synthesis in Escherichia coli. Annu Rev Microbiol 50: 645–677.
  3. 3. Bremer H, Dennis P (1996) Modulation of Chemical Composition and Other Parameters of the Cell by Growth Rate. In: Neidhardt FC, Curtiss RI, Ingraham JL, Lin ECC, Low KB, editors. Escherichia coli and Salmonella: Cellular and Molecular Biology. Washington, DC: ASM Press. pp. 1553–1569.
  4. 4. Cooper S, Helmstetter CE (1968) Chromosome replication and the division cycle of Escherichia coli B/r. J Mol Biol 31: 519–540.
  5. 5. French S (1992) Consequences of replication fork movement through transcription units in vivo. Science 258: 1362–1365.
  6. 6. Pomerantz RT, O'Donnell M (2008) The replisome uses mRNA as a primer after colliding with RNA polymerase. Nature 456: 762–766.
  7. 7. Yao NY, O'Donnell M (2009) Replisome structure and conformational dynamics underlie fork progression past obstacles. Curr Opin Cell Biol 21: 336–343.
  8. 8. Mirkin EV, Mirkin SM (2007) Replication fork stalling at natural impediments. Microbiol Mol Biol Rev 71: 13–35.
  9. 9. Rocha EP (2008) The organization of the bacterial genome. Annu Rev Genet 42: 211–233.
  10. 10. Liu B, Alberts BM (1995) Head-on collision between a DNA replication apparatus and RNA polymerase transcription complex. Science 267: 1131–1137.
  11. 11. Deshpande AM, Newlon CS (1996) DNA replication fork pause sites dependent on transcription. Science 272: 1030–1033.
  12. 12. Elias-Arnanz M, Salas M (1999) Resolution of head-on collisions between the transcription machinery and bacteriophage phi29 DNA polymerase is dependent on RNA polymerase translocation. Embo J 18: 5675–5682.
  13. 13. Azvolinsky A, Giresi PG, Lieb JD, Zakian VA (2009) Highly transcribed RNA polymerase II genes are impediments to replication fork progression in Saccharomyces cerevisiae. Mol Cell 34: 722–734.
  14. 14. Wellinger RE, Prado F, Aguilera A (2006) Replication fork progression is impaired by transcription in hyperrecombinant yeast cells lacking a functional THO complex. Mol Cell Biol 26: 3327–3334.
  15. 15. Trautinger BW, Jaktaji RP, Rusakova E, Lloyd RG (2005) RNA polymerase modulators and DNA repair activities resolve conflicts between DNA replication and transcription. Mol Cell 19: 247–258.
  16. 16. Boubakri H, de Septenville AL, Viguera E, Michel B (2009) The helicases DinG, Rep and UvrD cooperate to promote replication across transcription units in vivo. Embo J.
  17. 17. Brewer BJ (1988) When polymerases collide: replication and the transcriptional organization of the E. coli chromosome. Cell 53: 679–686.
  18. 18. Rocha EP (2004) The replication-related organization of bacterial genomes. Microbiology 150: 1609–1627.
  19. 19. Rocha EP, Danchin A (2003) Essentiality, not expressiveness, drives gene-strand bias in bacteria. Nat Genet 34: 377–378.
  20. 20. Guy L, Roten CA (2004) Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication. Gene 340: 45–52.
  21. 21. Ellwood M, Nomura M (1982) Chromosomal locations of the genes for rRNA in Escherichia coli K-12. J Bacteriol 149: 458–468.
  22. 22. Kunst F, Ogasawara N, Moszer I, Albertini AM, Alloni G, et al. (1997) The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature 390: 249–256.
  23. 23. Blattner FR, Plunkett G 3rd, Bloch CA, Perna NT, Burland V, et al. (1997) The complete genome sequence of Escherichia coli K-12. Science 277: 1453–1474.
  24. 24. Rocha E (2002) Is there a role for replication fork asymmetry in the distribution of genes in bacterial genomes? Trends Microbiol 10: 393–395.
  25. 25. Zeigler DR, Dean DH (1990) Orientation of genes in the Bacillus subtilis chromosome. Genetics 125: 703–708.
  26. 26. Brewer BJ, Fangman WL (1988) A replication fork barrier at the 3′ end of yeast ribosomal RNA genes. Cell 55: 637–643.
  27. 27. Price MN, Alm EJ, Arkin AP (2005) Interruptions in gene expression drive highly expressed operons to the leading strand of DNA replication. Nucleic Acids Res 33: 3224–3234.
  28. 28. Huvet M, Nicolay S, Touchon M, Audit B, d'Aubenton-Carafa Y, et al. (2007) Human gene organization driven by the coordination of replication and transcription. Genome Res 17: 1278–1285.
  29. 29. Omont N, Kepes F (2004) Transcription/replication collisions cause bacterial transcription units to be longer on the leading strand of replication. Bioinformatics 20: 2719–2725.
  30. 30. Mirkin EV, Mirkin SM (2005) Mechanisms of transcription-replication collisions in bacteria. Mol Cell Biol 25: 888–895.
  31. 31. Wang JD, Berkmen MB, Grossman AD (2007) Genome-wide coorientation of replication and transcription reduces adverse effects on replication in Bacillus subtilis. Proc Natl Acad Sci U S A 104: 5608–5613.
  32. 32. Prado F, Aguilera A (2005) Impairment of replication fork progression mediates RNA polII transcription-associated recombination. Embo J 24: 1267–1276.
  33. 33. Torres JZ, Bessler JB, Zakian VA (2004) Local chromatin structure at the ribosomal DNA causes replication fork pausing and genome instability in the absence of the S. cerevisiae DNA helicase Rrm3p. Genes Dev 18: 498–503.
  34. 34. Vilette D, Ehrlich SD, Michel B (1996) Transcription-induced deletions in plasmid vectors: M13 DNA replication as a source of instability. Mol Gen Genet 252: 398–403.
  35. 35. Kim N, Abdulovic AL, Gealy R, Lippert MJ, Jinks-Robertson S (2007) Transcription-associated mutagenesis in yeast is directly proportional to the level of gene expression and influenced by the direction of DNA replication. DNA Repair (Amst) 6: 1285–1296.
  36. 36. Dworkin J, Losick R (2002) Does RNA polymerase help drive chromosome segregation in bacteria? Proc Natl Acad Sci U S A 99: 14089–14094.
  37. 37. Rocha EP, Fralick J, Vediyappan G, Danchin A, Norris V (2003) A strand-specific model for chromosome segregation in bacteria. Mol Microbiol 49: 895–903.
  38. 38. Esnault E, Valens M, Espeli O, Boccard F (2007) Chromosome structuring limits genome plasticity in Escherichia coli. PLoS Genet 3: e226.
  39. 39. Lesterlin C, Pages C, Dubarry N, Dasgupta S, Cornet F (2008) Asymmetry of chromosome Replichores renders the DNA translocase activity of FtsK essential for cell division and cell shape maintenance in Escherichia coli. PLoS Genet 4: e1000288.
  40. 40. Berkmen MB, Grossman AD (2007) Subcellular positioning of the origin region of the Bacillus subtilis chromosome is independent of sequences within oriC, the site of replication initiation, and the replication initiator DnaA. Mol Microbiol 63: 150–165.
  41. 41. Mendelson NH, Gross JD (1967) Characterization of a temperature-sensitive mutant of Bacillus subtilis defective in deoxyribonucleic acid replication. J Bacteriol 94: 1603–1608.
  42. 42. Karamata D, Gross JD (1970) Isolation and genetic analysis of temperature-sensitive mutants of B. subtilis defective in DNA synthesis. Mol Gen Genet 108: 277–287.
  43. 43. Breier AM, Weier HU, Cozzarelli NR (2005) Independence of replisomes in Escherichia coli chromosomal replication. Proc Natl Acad Sci U S A 102: 3942–3947.
  44. 44. Wang JD, Sanders GM, Grossman AD (2007) Nutritional control of elongation of DNA replication by (p)ppGpp. Cell 128: 865–875.
  45. 45. Simmons LA, Grossman AD, Walker GC (2007) Replication is required for the RecA localization response to DNA damage in Bacillus subtilis. Proc Natl Acad Sci U S A 104: 1360–1365.
  46. 46. Britton RA, Kuster-Schock E, Auchtung TA, Grossman AD (2007) SOS induction in a subpopulation of structural maintenance of chromosome (Smc) mutant cells in Bacillus subtilis. J Bacteriol 189: 4359–4366.
  47. 47. Maughan H, Galeano B, Nicholson WL (2004) Novel rpoB mutations conferring rifampin resistance on Bacillus subtilis: global effects on growth, competence, sporulation, and germination. J Bacteriol 186: 2481–2486.
  48. 48. Ingham CJ, Furneaux PA (2000) Mutations in the ss subunit of the Bacillus subtilis RNA polymerase that confer both rifampicin resistance and hypersensitivity to NusG. Microbiology 146 Pt 12: 3041–3049.
  49. 49. Rosche WA, Foster PL (2000) Determining mutation rates in bacterial populations. Methods 20: 4–17.
  50. 50. Toda T, Tanaka T, Itaya M (1996) A method to invert DNA segments of the Bacillus subtilis 168 genome by recombination between two homologous sequences. Biosci Biotechnol Biochem 60: 773–778.
  51. 51. Segall A, Mahan MJ, Roth JR (1988) Rearrangement of the bacterial chromosome: forbidden inversions. Science 241: 1314–1318.
  52. 52. Rebollo JE, Francois V, Louarn JM (1988) Detection and possible role of two large nondivisible zones on the Escherichia coli chromosome. Proc Natl Acad Sci U S A 85: 9391–9395.
  53. 53. Konrad EB (1977) Method for the isolation of Escherichia coli mutants with enhanced recombination between chromosomal duplications. J Bacteriol 130: 167–172.
  54. 54. Hill CW, Gray JA (1988) Effects of chromosomal inversion on cell fitness in Escherichia coli K-12. Genetics 119: 771–778.
  55. 55. Campo N, Dias MJ, Daveran-Mingot ML, Ritzenthaler P, Le Bourgeois P (2004) Chromosomal constraints in Gram-positive bacteria revealed by artificial inversions. Mol Microbiol 51: 511–522.
  56. 56. Lesterlin C, Mercier R, Boccard F, Barre FX, Cornet F (2005) Roles for replichores and macrodomains in segregation of the Escherichia coli chromosome. EMBO Rep 6: 557–562.
  57. 57. Eisen JA, Heidelberg JF, White O, Salzberg SL (2000) Evidence for symmetric chromosomal inversions around the replication origin in bacteria. Genome Biol 1: RESEARCH0011.
  58. 58. Segre D, Deluna A, Church GM, Kishony R (2005) Modular epistasis in yeast metabolism. Nat Genet 37: 77–83.
  59. 59. Elena SF, Lenski RE (1997) Test of synergistic interactions among deleterious mutations in bacteria. Nature 390: 395–398.
  60. 60. Chandler MG, Pritchard RH (1975) The effect of gene concentration and relative gene dosage on gene output in Escherichia coli. Mol Gen Genet 138: 127–141.
  61. 61. Sullivan NL, Marquis KA, Rudner DZ (2009) Recruitment of SMC by ParB-parS organizes the origin region and promotes efficient chromosome segregation. Cell 137: 697–707.
  62. 62. Gruber S, Errington J (2009) Recruitment of condensin to replication origin regions by ParB/SpoOJ promotes chromosome segregation in B. subtilis. Cell 137: 685–696.
  63. 63. Lin DC, Grossman AD (1998) Identification and characterization of a bacterial chromosome partitioning site. Cell 92: 675–685.
  64. 64. Murray HD, Schneider DA, Gourse RL (2003) Control of rRNA expression by small molecules is dynamic and nonredundant. Mol Cell 12: 125–134.
  65. 65. Krasilnikova MM, Samadashwily GM, Krasilnikov AS, Mirkin SM (1998) Transcription through a simple DNA repeat blocks replication elongation. Embo J 17: 5095–5102.
  66. 66. Gomez-Gonzalez B, Felipe-Abrio I, Aguilera A (2009) The S-Phase checkpoint is required to respond to R-loops accumulated in THO mutants. Mol Cell Biol.
  67. 67. Olavarrieta L, Hernandez P, Krimer DB, Schvartzman JB (2002) DNA knotting caused by head-on collision of transcription and replication. J Mol Biol 322: 1–6.
  68. 68. Liu LF, Wang JC (1987) Supercoiling of the DNA template during transcription. Proc Natl Acad Sci U S A 84: 7024–7027.
  69. 69. Possoz C, Filipe SR, Grainge I, Sherratt DJ (2006) Tracking of controlled Escherichia coli replication fork stalling and restart at repressor-bound DNA in vivo. Embo J 25: 2596–2604.
  70. 70. Yasbin RE (1977) DNA repair in Bacillus subtilis. I. The presence of an inducible system. Mol Gen Genet 153: 211–218.
  71. 71. Simmons LA, Goranov AI, Kobayashi H, Davies BW, Yuan DS, et al. (2009) Comparison of responses to double-strand breaks between Escherichia coli and Bacillus subtilis reveals different requirements for SOS induction. J Bacteriol 191: 1152–1161.
  72. 72. Kuzminov A (1995) Collapse and repair of replication forks in Escherichia coli. Mol Microbiol 16: 373–384.
  73. 73. Bidnenko V, Ehrlich SD, Michel B (2002) Replication fork collapse at replication terminator sequences. Embo J 21: 3898–3907.
  74. 74. Couturier E, Rocha EP (2006) Replication-associated gene dosage effects shape the genomes of fast-growing bacteria but only for transcription and translation genes. Mol Microbiol 59: 1506–1518.
  75. 75. Saveson CJ, Lovett ST (1997) Enhanced deletion formation by aberrant DNA replication in Escherichia coli. Genetics 146: 457–470.
  76. 76. Ponder RG, Fonville NC, Rosenberg SM (2005) A switch from high-fidelity to error-prone DNA double-strand break repair underlies stress-induced mutation. Mol Cell 19: 791–804.
  77. 77. Jiang Q, Karata K, Woodgate R, Cox MM, Goodman MF (2009) The active form of DNA polymerase V is UmuD'2C-RecA-ATP. Nature 460: 359–363.
  78. 78. Reimers JM, Schmidt KH, Longacre A, Reschke DK, Wright BE (2004) Increased transcription rates correlate with increased reversion rates in leuB and argH Escherichia coli auxotrophs. Microbiology 150: 1457–1466.
  79. 79. Wright BE, Longacre A, Reimers JM (1999) Hypermutation in derepressed operons of Escherichia coli K12. Proc Natl Acad Sci U S A 96: 5089–5094.
  80. 80. Ross C, Pybus C, Pedraza-Reyes M, Sung HM, Yasbin RE, et al. (2006) Novel role of mfd: effects on stationary-phase mutagenesis in Bacillus subtilis. J Bacteriol 188: 7512–7520.
  81. 81. Aguilera A, Gomez-Gonzalez B (2008) Genome instability: a mechanistic view of its causes and consequences. Nat Rev Genet 9: 204–217.
  82. 82. Cohen SE, Godoy VG, Walker GC (2009) Transcriptional modulator NusA interacts with translesion DNA polymerases in Escherichia coli. J Bacteriol 191: 665–672.
  83. 83. Vasantha N, Freese E (1980) Enzyme changes during Bacillus subtilis sporulation caused by deprivation of guanine nucleotides. J Bacteriol 144: 1119–1125.
  84. 84. Harwood CR, Cutting SM (1990) Molecular Biological Methods for Bacillus subtilis. Chichester, U.K.: Wiley.
  85. 85. Dean DR, Hoch JA, Aronson AI (1977) Alteration of the Bacillus subtilis glutamine synthetase results in overproduction of the enzyme. J Bacteriol 131: 981–987.
  86. 86. Yasbin RE, Fields PI, Andersen BJ (1980) Properties of Bacillus subtilis 168 derivatives freed of their natural prophages. Gene 12: 155–159.
  87. 87. Auchtung JM, Lee CA, Monson RE, Lehman AP, Grossman AD (2005) Regulation of a Bacillus subtilis mobile genetic element by intercellular signaling and the global DNA damage response. Proc Natl Acad Sci U S A 102: 12554–12559.
  88. 88. Wu LJ, Errington J (2002) A large dispersed chromosomal region required for chromosome segregation in sporulating cells of Bacillus subtilis. Embo J 21: 4001–4011.
  89. 89. Youngman P, Poth H, Green B, York K, Olmedo G, et al. (1989) Methods for genetic manipulation, cloning, and functional analysis of sporulation genes in Bacillus subtilis. In: Smith I, Slepecky R, Setlow P, editors. Regulation of procaryotic development. Washington, D.C.: American Society for Microbiology. pp. 65–87.
  90. 90. Itaya M, Kondo K, Tanaka T (1989) A neomycin resistance gene cassette selectable in a single copy state in the Bacillus subtilis chromosome. Nucleic Acids Res 17: 4410.