Co-Orientation of Replication and Transcription Preserves Genome Integrity
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.
Published in the journal:
Co-Orientation of Replication and Transcription Preserves Genome Integrity. PLoS Genet 6(1): e32767. doi:10.1371/journal.pgen.1000810
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1000810
Summary
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.
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].
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].
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].
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.
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).
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).
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.
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
Zdroje
1. MiuraA
KruegerJH
ItohS
de BoerHA
NomuraM
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. GourseRL
GaalT
BartlettMS
ApplemanJA
RossW
1996 rRNA transcription and growth rate-dependent regulation of ribosome synthesis in Escherichia coli. Annu Rev Microbiol 50 645 677
3. BremerH
DennisP
1996 Modulation of Chemical Composition and Other Parameters of the Cell by Growth Rate.
NeidhardtFC
CurtissRI
IngrahamJL
LinECC
LowKB
Escherichia coli and Salmonella: Cellular and Molecular Biology Washington, DC ASM Press 1553 1569
4. CooperS
HelmstetterCE
1968 Chromosome replication and the division cycle of Escherichia coli B/r. J Mol Biol 31 519 540
5. FrenchS
1992 Consequences of replication fork movement through transcription units in vivo. Science 258 1362 1365
6. PomerantzRT
O'DonnellM
2008 The replisome uses mRNA as a primer after colliding with RNA polymerase. Nature 456 762 766
7. YaoNY
O'DonnellM
2009 Replisome structure and conformational dynamics underlie fork progression past obstacles. Curr Opin Cell Biol 21 336 343
8. MirkinEV
MirkinSM
2007 Replication fork stalling at natural impediments. Microbiol Mol Biol Rev 71 13 35
9. RochaEP
2008 The organization of the bacterial genome. Annu Rev Genet 42 211 233
10. LiuB
AlbertsBM
1995 Head-on collision between a DNA replication apparatus and RNA polymerase transcription complex. Science 267 1131 1137
11. DeshpandeAM
NewlonCS
1996 DNA replication fork pause sites dependent on transcription. Science 272 1030 1033
12. Elias-ArnanzM
SalasM
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. AzvolinskyA
GiresiPG
LiebJD
ZakianVA
2009 Highly transcribed RNA polymerase II genes are impediments to replication fork progression in Saccharomyces cerevisiae. Mol Cell 34 722 734
14. WellingerRE
PradoF
AguileraA
2006 Replication fork progression is impaired by transcription in hyperrecombinant yeast cells lacking a functional THO complex. Mol Cell Biol 26 3327 3334
15. TrautingerBW
JaktajiRP
RusakovaE
LloydRG
2005 RNA polymerase modulators and DNA repair activities resolve conflicts between DNA replication and transcription. Mol Cell 19 247 258
16. BoubakriH
de SeptenvilleAL
VigueraE
MichelB
2009 The helicases DinG, Rep and UvrD cooperate to promote replication across transcription units in vivo. Embo J
17. BrewerBJ
1988 When polymerases collide: replication and the transcriptional organization of the E. coli chromosome. Cell 53 679 686
18. RochaEP
2004 The replication-related organization of bacterial genomes. Microbiology 150 1609 1627
19. RochaEP
DanchinA
2003 Essentiality, not expressiveness, drives gene-strand bias in bacteria. Nat Genet 34 377 378
20. GuyL
RotenCA
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. EllwoodM
NomuraM
1982 Chromosomal locations of the genes for rRNA in Escherichia coli K-12. J Bacteriol 149 458 468
22. KunstF
OgasawaraN
MoszerI
AlbertiniAM
AlloniG
1997 The complete genome sequence of the gram-positive bacterium Bacillus subtilis. Nature 390 249 256
23. BlattnerFR
PlunkettG3rd
BlochCA
PernaNT
BurlandV
1997 The complete genome sequence of Escherichia coli K-12. Science 277 1453 1474
24. RochaE
2002 Is there a role for replication fork asymmetry in the distribution of genes in bacterial genomes? Trends Microbiol 10 393 395
25. ZeiglerDR
DeanDH
1990 Orientation of genes in the Bacillus subtilis chromosome. Genetics 125 703 708
26. BrewerBJ
FangmanWL
1988 A replication fork barrier at the 3′ end of yeast ribosomal RNA genes. Cell 55 637 643
27. PriceMN
AlmEJ
ArkinAP
2005 Interruptions in gene expression drive highly expressed operons to the leading strand of DNA replication. Nucleic Acids Res 33 3224 3234
28. HuvetM
NicolayS
TouchonM
AuditB
d'Aubenton-CarafaY
2007 Human gene organization driven by the coordination of replication and transcription. Genome Res 17 1278 1285
29. OmontN
KepesF
2004 Transcription/replication collisions cause bacterial transcription units to be longer on the leading strand of replication. Bioinformatics 20 2719 2725
30. MirkinEV
MirkinSM
2005 Mechanisms of transcription-replication collisions in bacteria. Mol Cell Biol 25 888 895
31. WangJD
BerkmenMB
GrossmanAD
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. PradoF
AguileraA
2005 Impairment of replication fork progression mediates RNA polII transcription-associated recombination. Embo J 24 1267 1276
33. TorresJZ
BesslerJB
ZakianVA
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. ViletteD
EhrlichSD
MichelB
1996 Transcription-induced deletions in plasmid vectors: M13 DNA replication as a source of instability. Mol Gen Genet 252 398 403
35. KimN
AbdulovicAL
GealyR
LippertMJ
Jinks-RobertsonS
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. DworkinJ
LosickR
2002 Does RNA polymerase help drive chromosome segregation in bacteria? Proc Natl Acad Sci U S A 99 14089 14094
37. RochaEP
FralickJ
VediyappanG
DanchinA
NorrisV
2003 A strand-specific model for chromosome segregation in bacteria. Mol Microbiol 49 895 903
38. EsnaultE
ValensM
EspeliO
BoccardF
2007 Chromosome structuring limits genome plasticity in Escherichia coli. PLoS Genet 3 e226
39. LesterlinC
PagesC
DubarryN
DasguptaS
CornetF
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. BerkmenMB
GrossmanAD
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. MendelsonNH
GrossJD
1967 Characterization of a temperature-sensitive mutant of Bacillus subtilis defective in deoxyribonucleic acid replication. J Bacteriol 94 1603 1608
42. KaramataD
GrossJD
1970 Isolation and genetic analysis of temperature-sensitive mutants of B. subtilis defective in DNA synthesis. Mol Gen Genet 108 277 287
43. BreierAM
WeierHU
CozzarelliNR
2005 Independence of replisomes in Escherichia coli chromosomal replication. Proc Natl Acad Sci U S A 102 3942 3947
44. WangJD
SandersGM
GrossmanAD
2007 Nutritional control of elongation of DNA replication by (p)ppGpp. Cell 128 865 875
45. SimmonsLA
GrossmanAD
WalkerGC
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. BrittonRA
Kuster-SchockE
AuchtungTA
GrossmanAD
2007 SOS induction in a subpopulation of structural maintenance of chromosome (Smc) mutant cells in Bacillus subtilis. J Bacteriol 189 4359 4366
47. MaughanH
GaleanoB
NicholsonWL
2004 Novel rpoB mutations conferring rifampin resistance on Bacillus subtilis: global effects on growth, competence, sporulation, and germination. J Bacteriol 186 2481 2486
48. InghamCJ
FurneauxPA
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. RoscheWA
FosterPL
2000 Determining mutation rates in bacterial populations. Methods 20 4 17
50. TodaT
TanakaT
ItayaM
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. SegallA
MahanMJ
RothJR
1988 Rearrangement of the bacterial chromosome: forbidden inversions. Science 241 1314 1318
52. RebolloJE
FrancoisV
LouarnJM
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. KonradEB
1977 Method for the isolation of Escherichia coli mutants with enhanced recombination between chromosomal duplications. J Bacteriol 130 167 172
54. HillCW
GrayJA
1988 Effects of chromosomal inversion on cell fitness in Escherichia coli K-12. Genetics 119 771 778
55. CampoN
DiasMJ
Daveran-MingotML
RitzenthalerP
Le BourgeoisP
2004 Chromosomal constraints in Gram-positive bacteria revealed by artificial inversions. Mol Microbiol 51 511 522
56. LesterlinC
MercierR
BoccardF
BarreFX
CornetF
2005 Roles for replichores and macrodomains in segregation of the Escherichia coli chromosome. EMBO Rep 6 557 562
57. EisenJA
HeidelbergJF
WhiteO
SalzbergSL
2000 Evidence for symmetric chromosomal inversions around the replication origin in bacteria. Genome Biol 1 RESEARCH0011
58. SegreD
DelunaA
ChurchGM
KishonyR
2005 Modular epistasis in yeast metabolism. Nat Genet 37 77 83
59. ElenaSF
LenskiRE
1997 Test of synergistic interactions among deleterious mutations in bacteria. Nature 390 395 398
60. ChandlerMG
PritchardRH
1975 The effect of gene concentration and relative gene dosage on gene output in Escherichia coli. Mol Gen Genet 138 127 141
61. SullivanNL
MarquisKA
RudnerDZ
2009 Recruitment of SMC by ParB-parS organizes the origin region and promotes efficient chromosome segregation. Cell 137 697 707
62. GruberS
ErringtonJ
2009 Recruitment of condensin to replication origin regions by ParB/SpoOJ promotes chromosome segregation in B. subtilis. Cell 137 685 696
63. LinDC
GrossmanAD
1998 Identification and characterization of a bacterial chromosome partitioning site. Cell 92 675 685
64. MurrayHD
SchneiderDA
GourseRL
2003 Control of rRNA expression by small molecules is dynamic and nonredundant. Mol Cell 12 125 134
65. KrasilnikovaMM
SamadashwilyGM
KrasilnikovAS
MirkinSM
1998 Transcription through a simple DNA repeat blocks replication elongation. Embo J 17 5095 5102
66. Gomez-GonzalezB
Felipe-AbrioI
AguileraA
2009 The S-Phase checkpoint is required to respond to R-loops accumulated in THO mutants. Mol Cell Biol
67. OlavarrietaL
HernandezP
KrimerDB
SchvartzmanJB
2002 DNA knotting caused by head-on collision of transcription and replication. J Mol Biol 322 1 6
68. LiuLF
WangJC
1987 Supercoiling of the DNA template during transcription. Proc Natl Acad Sci U S A 84 7024 7027
69. PossozC
FilipeSR
GraingeI
SherrattDJ
2006 Tracking of controlled Escherichia coli replication fork stalling and restart at repressor-bound DNA in vivo. Embo J 25 2596 2604
70. YasbinRE
1977 DNA repair in Bacillus subtilis. I. The presence of an inducible system. Mol Gen Genet 153 211 218
71. SimmonsLA
GoranovAI
KobayashiH
DaviesBW
YuanDS
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. KuzminovA
1995 Collapse and repair of replication forks in Escherichia coli. Mol Microbiol 16 373 384
73. BidnenkoV
EhrlichSD
MichelB
2002 Replication fork collapse at replication terminator sequences. Embo J 21 3898 3907
74. CouturierE
RochaEP
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. SavesonCJ
LovettST
1997 Enhanced deletion formation by aberrant DNA replication in Escherichia coli. Genetics 146 457 470
76. PonderRG
FonvilleNC
RosenbergSM
2005 A switch from high-fidelity to error-prone DNA double-strand break repair underlies stress-induced mutation. Mol Cell 19 791 804
77. JiangQ
KarataK
WoodgateR
CoxMM
GoodmanMF
2009 The active form of DNA polymerase V is UmuD'2C-RecA-ATP. Nature 460 359 363
78. ReimersJM
SchmidtKH
LongacreA
ReschkeDK
WrightBE
2004 Increased transcription rates correlate with increased reversion rates in leuB and argH Escherichia coli auxotrophs. Microbiology 150 1457 1466
79. WrightBE
LongacreA
ReimersJM
1999 Hypermutation in derepressed operons of Escherichia coli K12. Proc Natl Acad Sci U S A 96 5089 5094
80. RossC
PybusC
Pedraza-ReyesM
SungHM
YasbinRE
2006 Novel role of mfd: effects on stationary-phase mutagenesis in Bacillus subtilis. J Bacteriol 188 7512 7520
81. AguileraA
Gomez-GonzalezB
2008 Genome instability: a mechanistic view of its causes and consequences. Nat Rev Genet 9 204 217
82. CohenSE
GodoyVG
WalkerGC
2009 Transcriptional modulator NusA interacts with translesion DNA polymerases in Escherichia coli. J Bacteriol 191 665 672
83. VasanthaN
FreeseE
1980 Enzyme changes during Bacillus subtilis sporulation caused by deprivation of guanine nucleotides. J Bacteriol 144 1119 1125
84. HarwoodCR
CuttingSM
1990 Molecular Biological Methods for Bacillus subtilis Chichester, U.K. Wiley
85. DeanDR
HochJA
AronsonAI
1977 Alteration of the Bacillus subtilis glutamine synthetase results in overproduction of the enzyme. J Bacteriol 131 981 987
86. YasbinRE
FieldsPI
AndersenBJ
1980 Properties of Bacillus subtilis 168 derivatives freed of their natural prophages. Gene 12 155 159
87. AuchtungJM
LeeCA
MonsonRE
LehmanAP
GrossmanAD
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. WuLJ
ErringtonJ
2002 A large dispersed chromosomal region required for chromosome segregation in sporulating cells of Bacillus subtilis. Embo J 21 4001 4011
89. YoungmanP
PothH
GreenB
YorkK
OlmedoG
1989 Methods for genetic manipulation, cloning, and functional analysis of sporulation genes in Bacillus subtilis.
SmithI
SlepeckyR
SetlowP
Regulation of procaryotic development Washington, D.C. American Society for Microbiology 65 87
90. ItayaM
KondoK
TanakaT
1989 A neomycin resistance gene cassette selectable in a single copy state in the Bacillus subtilis chromosome. Nucleic Acids Res 17 4410
Štítky
Genetika Reprodukčná medicínaČlánok vyšiel v časopise
PLOS Genetics
2010 Číslo 1
- Je „freeze-all“ pro všechny? Odborníci na fertilitu diskutovali na virtuálním summitu
- Gynekologové a odborníci na reprodukční medicínu se sejdou na prvním virtuálním summitu
Najčítanejšie v tomto čísle
- A Major Role of the RecFOR Pathway in DNA Double-Strand-Break Repair through ESDSA in
- Kidney Development in the Absence of and Requires
- The Werner Syndrome Protein Functions Upstream of ATR and ATM in Response to DNA Replication Inhibition and Double-Strand DNA Breaks
- Alternative Epigenetic Chromatin States of Polycomb Target Genes