Escherichia coli DNA distributions measured by flow cytometry and compared with theoretical computer simulations

Fifty-Five Years of Research on B, C and D in Escherichia coli

The basic properties of the Escherichia coli duplication process can be defined by two time periods: C, the time for a round of chromosome replication, and D, the time between the end of a round of replication and cell division. Given the durations of these periods, the pattern of chromosome replication during the cell cycle can be determined for cells growing with any doubling time. In the 55 years since these parameters were identified, there have been numerous investigations into their durations and into the elements that determine their initiations. In this review, I discuss the history of our involvement in these studies from the very beginning, some of what has been learned over the years by measuring the durations of C and D, and what might be learned with additional investigations.

Interdependent progression of bidirectional sister replisomes in E. coli

Bidirectional DNA replication complexes initiated from the same origin remain colocalized in a factory configuration for part or all their lifetimes. However, there is little evidence that sister replisomes are functionally interdependent, and the consequence of factory replication is unknown. Here, we investigated the functional relationship between sister replisomes in Escherichia coli, which naturally exhibits both factory and solitary configurations in the same replication cycle. Using an inducible transcription factor roadblocking system, we found that blocking one replisome caused a significant decrease in overall progression and velocity of the sister replisome. Remarkably, progression was impaired only if the block occurred while sister replisomes were still in a factory configuration - blocking one fork had no significant effect on the other replisome when sister replisomes were physically separate. Disruption of factory replication also led to increased fork stalling and requirement of fork restart mechanisms. These results suggest that physical association between sister replisomes is important for establishing an efficient and uninterrupted replication program. We discuss the implications of our findings on mechanisms of replication factory structure and function, and cellular strategies of replicating problematic DNA such as highly transcribed segments.

Involvement of multiple influx and efflux transporters in the accumulation of cationic fluorescent dyes by Escherichia coli

Background

It is widely believed that most xenobiotics cross biomembranes by diffusing through the phospholipid bilayer, and that the use of protein transporters is an occasional adjunct. According to an alternative view, phospholipid bilayer transport is negligible, and several different transporters may be involved in the uptake of an individual molecular type. We recognise here that the availability of gene knockout collections allows one to assess the contributions of all potential transporters, and flow cytometry based on fluorescence provides a convenient high-throughput assay for xenobiotic uptake in individual cells.

Results

We used high-throughput flow cytometry to assess the ability of individual gene knockout strains of E coli to take up two membrane-permeable, cationic fluorescent dyes, namely the carbocyanine diS-C3(5) and the DNA dye SYBR Green. Individual strains showed a large range of distributions of uptake. The range of modal steady-state uptakes for the carbocyanine between the different strains was 36-fold. Knockouts of the ATP synthase α- and β-subunits greatly inhibited uptake, implying that most uptake was ATP-driven rather than being driven by a membrane potential. Dozens of transporters changed the steady-state uptake of the dye by more than 50% with respect to that of the wild type, in either direction (increased or decreased); knockouts of known influx and efflux transporters behaved as expected, giving credence to the general strategy. Many of the knockouts with the most reduced uptake were transporter genes of unknown function (‘y-genes’). Similarly, several overexpression variants in the ‘ASKA’ collection had the anticipated, opposite effects. Similar results were obtained with SYBR Green (the range being approximately 69-fold). Although it too contains a benzothiazole motif there was negligible correlation between its uptake and that of the carbocyanine when compared across the various strains (although the membrane potential is presumably the same in each case).

Conclusions

Overall, we conclude that the uptake of these dyes may be catalysed by a great many transporters of putatively broad and presently unknown specificity, and that the very large range between the ‘lowest’ and the ‘highest’ levels of uptake, even in knockouts of just single genes, implies strongly that phospholipid bilayer transport is indeed negligible. This work also casts serious doubt upon the use of such dyes as quantitative stains for representing either bioenergetic parameters or the amount of cellular DNA in unfixed cells (in vivo). By contrast, it opens up their potential use as transporter assay substrates in high-throughput screening.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1561-0) contains supplementary material, which is available to authorized users.

Neutral mechanisms and niche differentiation in steady-state insular microbial communities revealed by single cell analysis

In completely insular microbial communities, evolution of community structure cannot be shaped by the immigration of new members. In addition, when those communities are run in steady state, the influence of environmental factors on their assembly is reduced. Therefore, one would expect similar community structures under steady‐state conditions. Yet, in parallel setups, variability does occur. To reveal ecological mechanisms behind this phenomenon, five parallel reactors were studied at the single‐cell level for about 100 generations and community structure variations were quantified by ecological measures. Whether community variability can be controlled was tested by implementing soft temperature stressors as potential synchronizers. The low slope of the lognormal rank‐order abundance curves indicated a predominance of neutral mechanisms, i.e., where species identity plays no role. Variations in abundance ranks of subcommunities and increase in inter‐community pairwise β‐diversity over time support this. Niche differentiation was also observed, as indicated by steeper geometric‐like rank‐order abundance curves and increased numbers of correlations between abiotic and biotic parameters during initial adaptation and after disturbances. Still, neutral forces dominated community assembly. Our findings suggest that complex microbial communities in insular steady‐state environments can be difficult to synchronize and maintained in their original or desired structure, as they are non‐equilibrium systems.

Fundamental principles in bacterial physiology-history, recent progress, and the future with focus on cell size control: a review

Bacterial physiology is a branch of biology that aims to understand overarching principles of cellular reproduction. Many important issues in bacterial physiology are inherently quantitative, and major contributors to the field have often brought together tools and ways of thinking from multiple disciplines. This article presents a comprehensive overview of major ideas and approaches developed since the early 20th century for anyone who is interested in the fundamental problems in bacterial physiology. This article is divided into two parts. In the first part (sections 1-3), we review the first 'golden era' of bacterial physiology from the 1940s to early 1970s and provide a complete list of major references from that period. In the second part (sections 4-7), we explain how the pioneering work from the first golden era has influenced various rediscoveries of general quantitative principles and significant further development in modern bacterial physiology. Specifically, section 4 presents the history and current progress of the 'adder' principle of cell size homeostasis. Section 5 discusses the implications of coarse-graining the cellular protein composition, and how the coarse-grained proteome 'sectors' re-balance under different growth conditions. Section 6 focuses on physiological invariants, and explains how they are the key to understanding the coordination between growth and the cell cycle underlying cell size control in steady-state growth. Section 7 overviews how the temporal organization of all the internal processes enables balanced growth. In the final section 8, we conclude by discussing the remaining challenges for the future in the field.

Predicting the Dynamics and Heterogeneity of Genomic DNA Content within Bacterial Populations across Variable Growth Regimes

For many applications in microbial synthetic biology, optimizing a desired function requires careful tuning of the degree to which various genes are expressed. One challenge for predicting such effects or interpreting typical characterization experiments is that in bacteria such as E. coli, genome copy number varies widely across different phases and rates of growth, which also impacts how and when genes are expressed from different loci. While such phenomena are relatively well-understood at a mechanistic level, our quantitative understanding of such processes is essentially limited to ideal exponential growth. In contrast, common experimental phenomena such as growth on heterogeneous media, metabolic adaptation, and oxygen restriction all cause substantial deviations from ideal exponential growth, particularly as cultures approach the higher densities at which industrial biomanufacturing and even routine screening experiments are conducted. To meet the need for predicting and explaining how gene dosage impacts cellular functions outside of exponential growth, we here report a novel modeling strategy that leverages agent-based simulation and high performance computing to robustly predict the dynamics and heterogeneity of genomic DNA content within bacterial populations across variable growth regimes. We show that by feeding routine experimental data, such as optical density time series, into our heterogeneous multiphasic growth simulator, we can predict genomic DNA distributions over a range of nonexponential growth conditions. This modeling strategy provides an important advance in the ability of synthetic biologists to evaluate the role of genomic DNA content and heterogeneity in affecting the performance of existing or engineered microbial functions.

Lagrangian Trajectories to Predict the Formation of Population Heterogeneity in Large-Scale Bioreactors

Successful scale-up of bioprocesses requires that laboratory-scale performance is equally achieved during large-scale production to meet economic constraints. In industry, heuristic approaches are often applied, making use of physical scale-up criteria that do not consider cellular needs or properties. As a consequence, large-scale productivities, conversion yields, or product purities are often deteriorated, which may prevent economic success. The occurrence of population heterogeneity in large-scale production may be the reason for underperformance. In this study, an in silico method to predict the formation of population heterogeneity by combining computational fluid dynamics (CFD) with a cell cycle model of Pseudomonas putida KT2440 was developed. The glucose gradient and flow field of a 54,000 L stirred tank reactor were generated with the Euler approach, and bacterial movement was simulated as Lagrange particles. The latter were statistically evaluated using a cell cycle model. Accordingly, 72% of all cells were found to switch between standard and multifork replication, and 10% were likely to undergo massive, transcriptional adaptations to respond to extracellular starving conditions. At the same time, 56% of all cells replicated very fast, with µ ≥ 0.3 h^-1 performing multifork replication. The population showed very strong heterogeneity, as indicated by the observation that 52.9% showed higher than average adenosine triphosphate (ATP) maintenance demands (12.2%, up to 1.5 fold). These results underline the potential of CFD linked to structured cell cycle models for predicting large-scale heterogeneity in silico and ab initio.

Threshold effect of growth rate on population variability of Escherichia coli cell lengths

A long-standing question in biology is the effect of growth on cell size. Here, we estimate the effect of Escherichia coli growth rate (r) on population cell size distributions by estimating the coefficient of variation of cell lengths (CV_L) from image analysis of fixed cells in DIC microscopy. We find that the CV_L is constant at growth rates less than one division per hour, whereas above this threshold, CV_L increases with an increase in the growth rate. We hypothesize that stochastic inhibition of cell division owing to replication stalling by a RecA-dependent mechanism, combined with the growth rate threshold of multi-fork replication (according to Cooper and Helmstetter), could form the basis of such a threshold effect. We proceed to test our hypothesis by increasing the frequency of stochastic stalling of replication forks with hydroxyurea (HU) treatment and find that cell length variability increases only when the growth rate exceeds this threshold. The population effect is also reproduced in single-cell studies using agar-pad cultures and 'mother machine'-based experiments to achieve synchrony. To test the role of RecA, critical for the repair of stalled replication forks, we examine the CV_L of E. coli ΔrecA cells. We find cell length variability in the mutant to be greater than wild-type, a phenotype that is rescued by plasmid-based RecA expression. Additionally, we find that RecA-GFP protein recruitment to nucleoids is more frequent at growth rates exceeding the growth rate threshold and is further enhanced on HU treatment. Thus, we find growth rates greater than a threshold result in increased E. coli cell lengths in the population, and this effect is, at least in part, mediated by RecA recruitment to the nucleoid and stochastic inhibition of division.

Threshold effect of growth rate on population variability of Escherichia coli cell lengths

A long-standing question in biology is the effect of growth on cell size. Here, we estimate the effect of Escherichia coli growth rate (r) on population cell size distributions by estimating the coefficient of variation of cell lengths (CV_L) from image analysis of fixed cells in DIC microscopy. We find that the CV_L is constant at growth rates less than one division per hour, whereas above this threshold, CV_L increases with an increase in the growth rate. We hypothesize that stochastic inhibition of cell division owing to replication stalling by a RecA-dependent mechanism, combined with the growth rate threshold of multi-fork replication (according to Cooper and Helmstetter), could form the basis of such a threshold effect. We proceed to test our hypothesis by increasing the frequency of stochastic stalling of replication forks with hydroxyurea (HU) treatment and find that cell length variability increases only when the growth rate exceeds this threshold. The population effect is also reproduced in single-cell studies using agar-pad cultures and 'mother machine'-based experiments to achieve synchrony. To test the role of RecA, critical for the repair of stalled replication forks, we examine the CV_L of E. coli ΔrecA cells. We find cell length variability in the mutant to be greater than wild-type, a phenotype that is rescued by plasmid-based RecA expression. Additionally, we find that RecA-GFP protein recruitment to nucleoids is more frequent at growth rates exceeding the growth rate threshold and is further enhanced on HU treatment. Thus, we find growth rates greater than a threshold result in increased E. coli cell lengths in the population, and this effect is, at least in part, mediated by RecA recruitment to the nucleoid and stochastic inhibition of division.

Automated Multi-Peak Tracking Kymography (AMTraK): A Tool to Quantify Sub-Cellular Dynamics with Sub-Pixel Accuracy

Kymographs or space-time plots are widely used in cell biology to reduce the dimensions of a time-series in microscopy for both qualitative and quantitative insight into spatio-temporal dynamics. While multiple tools for image kymography have been described before, quantification remains largely manual. Here, we describe a novel software tool for automated multi-peak tracking kymography (AMTraK), which uses peak information and distance minimization to track and automatically quantify kymographs, integrated in a GUI. The program takes fluorescence time-series data as an input and tracks contours in the kymographs based on intensity and gradient peaks. By integrating a branch-point detection method, it can be used to identify merging and splitting events of tracks, important in separation and coalescence events. In tests with synthetic images, we demonstrate sub-pixel positional accuracy of the program. We test the program by quantifying sub-cellular dynamics in rod-shaped bacteria, microtubule (MT) transport and vesicle dynamics. A time-series of E. coli cell division with labeled nucleoid DNA is used to identify the time-point and rate at which the nucleoid segregates. The mean velocity of microtubule (MT) gliding motility due to a recombinant kinesin motor is estimated as 0.5 μm/s, in agreement with published values, and comparable to estimates using software for nanometer precision filament-tracking. We proceed to employ AMTraK to analyze previously published time-series microscopy data where kymographs had been manually quantified: clathrin polymerization kinetics during vesicle formation and anterograde and retrograde transport in axons. AMTraK analysis not only reproduces the reported parameters, it also provides an objective and automated method for reproducible analysis of kymographs from in vitro and in vivo fluorescence microscopy time-series of sub-cellular dynamics.

A mechanistic stochastic framework for regulating bacterial cell division

How exponentially growing cells maintain size homeostasis is an important fundamental problem. Recent single-cell studies in prokaryotes have uncovered the adder principle, where cells add a fixed size (volume) from birth to division, irrespective of their size at birth. To mechanistically explain the adder principle, we consider a timekeeper protein that begins to get stochastically expressed after cell birth at a rate proportional to the volume. Cell-division time is formulated as the first-passage time for protein copy numbers to hit a fixed threshold. Consistent with data, the model predicts that the noise in division timing increases with size at birth. Intriguingly, our results show that the distribution of the volume added between successive cell-division events is independent of the newborn cell size. This was dramatically seen in experimental studies, where histograms of the added volume corresponding to different newborn sizes collapsed on top of each other. The model provides further insights consistent with experimental observations: the distribution of the added volume when scaled by its mean becomes invariant of the growth rate. In summary, our simple yet elegant model explains key experimental findings and suggests a mechanism for regulating both the mean and fluctuations in cell-division timing for controlling size.

Quantifying Impact of Chromosome Copy Number on Recombination in Escherichia coli

The ability to precisely and efficiently recombineer synthetic DNA into organisms of interest in a quantitative manner is a key requirement in genome engineering. Even though considerable effort has gone into the characterization of recombination in Escherichia coli, there is still substantial variability in reported recombination efficiencies. We hypothesized that this observed variability could, in part, be explained by the variability in chromosome copy number as well as the location of the replication forks relative to the recombination site. During rapid growth, E. coli cells may contain several pairs of open replication forks. While recombineered forks are resolving and segregating within the population, changes in apparent recombineering efficiency should be observed. In the case of dominant phenotypes, we predicted and then experimentally confirmed that the apparent recombination efficiency declined during recovery until complete segregation of recombineered and wild-type genomes had occurred. We observed the reverse trend for recessive phenotypes. The observed changes in apparent recombination efficiency were found to be in agreement with mathematical calculations based on our proposed mechanism. We also provide a model that can be used to estimate the total segregated recombination efficiency based on an initial efficiency and growth rate. These results emphasize the importance of employing quantitative strategies in the design of genome-scale engineering efforts.

Studies of Antibiotic Resistance of Beta-Lactamase Bacteria under Different Nutrition Limitations at the Single-Cell Level

Drug resistance involves many biological processes, including cell growth, cell communication, and cell cooperation. In the last few decades, bacterial drug resistance studies have made substantial progress. However, a major limitation of the traditional resistance study still exists: most of the studies have concentrated on the average behavior of enormous amounts of cells rather than surveying single cells with different phenotypes or genotypes. Here, we report our study of beta-lactamase bacterial drug resistance in a well-designed microfluidic device, which allows us to conduct more controllable experiments, such as controlling the nutrient concentration, switching the culture media, performing parallel experiments, observing single cells, and acquiring time-lapse images. By using GFP as a beta-lactamase indicator and acquiring time-lapse images at the single-cell level, we observed correlations between the bacterial heterogeneous phenotypes and their behavior in different culture media. The feedback loop between the growth rate and the beta-lactamase production suggests that the beta-lactamase bacteria are more resistant in a rich medium than in a relatively poor medium. In the poorest medium, the proportion of dormant cells may increase, which causes a lower death rate in the same generation. Our work may contribute to assaying the antibiotic resistance of pathogenic bacteria in heterogeneous complex media.

Topological Behavior of Plasmid DNA

The discovery of the B-form structure of DNA by Watson and Crick led to an explosion of research on nucleic acids in the fields of biochemistry, biophysics, and genetics. Powerful techniques were developed to reveal a myriad of different structural conformations that change B-DNA as it is transcribed, replicated, and recombined and as sister chromosomes are moved into new daughter cell compartments during cell division. This article links the original discoveries of superhelical structure and molecular topology to non-B form DNA structure and contemporary biochemical and biophysical techniques. The emphasis is on the power of plasmids for studying DNA structure and function. The conditions that trigger the formation of alternative DNA structures such as left-handed Z-DNA, inter- and intra-molecular triplexes, triple-stranded DNA, and linked catenanes and hemicatenanes are explained. The DNA dynamics and topological issues are detailed for stalled replication forks and for torsional and structural changes on DNA in front of and behind a transcription complex and a replisome. The complex and interconnected roles of topoisomerases and abundant small nucleoid association proteins are explained. And methods are described for comparing in vivo and in vitro reactions to probe and understand the temporal pathways of DNA and chromosome chemistry that occur inside living cells.

Chromosome segregation in Vibrio cholerae

The study of chromosome segregation is currently one of the most exciting research frontiers in cell biology. In this review, we discuss our current knowledge of the chromosome segregation process in Vibrio cholerae, based primarily on findings from fluorescence microscopy experiments. This bacterium is of special interest because of its eukaryotic feature of having a divided genome, a feature shared with 10% of known bacteria. We also discuss how the segregation mechanisms of V. cholerae compare with those in other bacteria, and highlight some of the remaining questions regarding the process of bacterial chromosome segregation.

Subpopulation-proteomics reveal growth rate, but not cell cycling, as a major impact on protein composition in Pseudomonas putida KT2440

Population heterogeneity occurring in industrial microbial bioprocesses is regarded as a putative effector causing performance loss in large scale. While the existence of subpopulations is a commonly accepted fact, their appearance and impact on process performance still remains rather unclear. During cell cycling, distinct subpopulations differing in cell division state and DNA content appear which contribute individually to the efficiency of the bioprocess. To identify stressed or impaired subpopulations, we analyzed the interplay of growth rate, cell cycle and phenotypic profile of subpopulations by using flow cytometry and cell sorting in conjunction with mass spectrometry based global proteomics. Adjusting distinct growth rates in chemostats with the model strain Pseudomonas putida KT2440, cells were differentiated by DNA content reflecting different cell cycle stages. The proteome of separated subpopulations at given growth rates was found to be highly similar, while different growth rates caused major changes of the protein inventory with respect to e.g. carbon storage, motility, lipid metabolism and the translational machinery. In conclusion, cells in various cell cycle stages at the same growth rate were found to have similar to identical proteome profiles showing no significant population heterogeneity on the proteome level. In contrast, the growth rate clearly determines the protein composition and therefore the metabolic strategy of the cells.

AspC-mediated aspartate metabolism coordinates the Escherichia coli cell cycle

Background

The fast-growing bacterial cell cycle consists of at least two independent cycles of chromosome replication and cell division. To ensure proper cell cycles and viability, chromosome replication and cell division must be coordinated. It has been suggested that metabolism could affect the Escherichia coli cell cycle, but the idea is still lacking solid evidences.

Methodology/Principle Findings

We found that absence of AspC, an aminotransferase that catalyzes synthesis of aspartate, led to generation of small cells with less origins and slow growth. In contrast, excess AspC was found to exert the opposite effect. Further analysis showed that AspC-mediated aspartate metabolism had a specific effect in the cell cycle, as only extra aspartate of the 20 amino acids triggered production of bigger cells with more origins per cell and faster growth. The amount of DnaA protein per cell was found to be changed in response to the availability of AspC. Depletion of (p)ppGpp by ΔrelAΔspoT led to a slight delay in initiation of replication, but did not change the replication pattern found in the ΔaspC mutant.

Conclusion/Significances

The results suggest that AspC-mediated metabolism of aspartate coordinates the E. coli cell cycle through altering the amount of the initiator protein DnaA per cell and the division signal UDP-glucose. Furthermore, AspC sequence conservation suggests similar functions in other organisms.

A Drosophila model to image phagosome maturation

Phagocytosis involves the internalization of extracellular material by invagination of the plasma membrane to form intracellular vesicles called phagosomes, which have functions that include pathogen degradation. The degradative properties of phagosomes are thought to be conferred by sequential fusion with endosomes and lysosomes; however, this maturation process has not been studied in vivo. We employed Drosophila hemocytes, which are similar to mammalian professional macrophages, to establish a model of phagosome maturation. Adult Drosophila females, carrying transgenic Rab7-GFP endosome and Lamp1-GFP lysosome markers, were injected with E. coli DH5α and the hemocytes were collected at 15, 30, 45 and 60 minutes after infection. In wild-type females, E. coli were detected within enlarged Rab7-GFP positive phagosomes at 15 to 45 minutes after infection; and were also observed in enlarged Lamp1-GFP positive phagolysosomes at 45 minutes. Two-photon imaging of hemocytes in vivo confirmed this vesicle morphology, including enlargement of Rab7-GFP and Lamp1-GFP structures that often appeared to protrude from hemocytes. The interaction of endosomes and lysosomes with E. coli phagosomes observed in Drosophila hemocytes was consistent with that previously described for phagosome maturation in human ex vivo macrophages. We also tested our model as a tool for genetic analysis using 14-3-3ε mutants, and demonstrated altered phagosome maturation with delayed E. coli internalization, trafficking and/or degradation. These findings demonstrate that Drosophila hemocytes provide an appropriate, genetically amenable, model for analyzing phagosome maturation ex vivo and in vivo.

An easy-to-use simulation program demonstrates variations in bacterial cell cycle parameters depending on medium and temperature

Many studies are performed on chromosome replication and segregation in Escherichia coli and other bacteria capable of complex replication with C phases spanning several generations. For such investigations an understanding of the replication patterns, including copy numbers of origins and replication forks, is crucial for correct interpretation of the results.Flow cytometry is an important tool for generation of experimental DNA distributions of cell populations. Here, a Visual Basic based simulation program was written for the computation of theoretical DNA distributions for different choices of cell cycle parameters (C and D phase durations, doubling time etc). These cell cycle parameters can be iterated until the best fit between the experimental and theoretical DNA histograms is obtained. The Excel file containing the simulation software is attached as supporting information.Cultures of Escherichia coli were grown at twelve different media and temperature conditions, with following measurements by flow cytometry and simulation of the DNA distributions. A good fit was found for each growth condition by use of our simulation program. The resulting cell cycle parameters displayed clear inter-media differences in replication patterns, but indicated a high degree of temperature independence for each medium. The exception was the poorest medium (acetate), where the cells grew with overlapping replication cycles at 42 °C, but without at the lower temperatures.We have developed an easy-to-use tool for determination of bacteria's cell cycle parameters, and consequently the cells' chromosome configurations. The procedure only requires DNA distribution measurements by flow cytometry. Use of this simulation program for E. coli cultures shows that even cells growing quite slowly can have overlapping replication cycles. It is therefore always important not only to assume cells' replication patterns, but to actually determine the cell cycle parameters when changing growth conditions.

Genome engineering in Vibrio cholerae: a feasible approach to address biological issues

Although bacteria with multipartite genomes are prevalent, our knowledge of the mechanisms maintaining their genome is very limited, and much remains to be learned about the structural and functional interrelationships of multiple chromosomes. Owing to its bi-chromosomal genome architecture and its importance in public health, Vibrio cholerae, the causative agent of cholera, has become a preferred model to study bacteria with multipartite genomes. However, most in vivo studies in V. cholerae have been hampered by its genome architecture, as it is difficult to give phenotypes to a specific chromosome. This difficulty was surmounted using a unique and powerful strategy based on massive rearrangement of prokaryotic genomes. We developed a site-specific recombination-based engineering tool, which allows targeted, oriented, and reciprocal DNA exchanges. Using this genetic tool, we obtained a panel of V. cholerae mutants with various genome configurations: one with a single chromosome, one with two chromosomes of equal size, and one with both chromosomes controlled by identical origins. We used these synthetic strains to address several biological questions--the specific case of the essentiality of Dam methylation in V. cholerae and the general question concerning bacteria carrying circular chromosomes--by looking at the effect of chromosome size on topological issues. In this article, we show that Dam, RctB, and ParA2/ParB2 are strictly essential for chrII origin maintenance, and we formally demonstrate that the formation of chromosome dimers increases exponentially with chromosome size.

Replication and segregation of an Escherichia coli chromosome with two replication origins

Characterized bacteria, unlike eukaryotes and some archaea, initiate replication bidirectionally from a single replication origin contained within a circular or linear chromosome. We constructed Escherichia coli cells with two WT origins separated by 1 Mb in their 4.64-Mb chromosome. Productive bidirectional replication initiated synchronously at both spatially separate origins. Newly replicated DNA from both origins was segregated sequentially as replication progressed, with two temporally and spatially separate replication termination events. Replication initiation occurred at a cell volume identical to that of cells with a single WT origin, showing that initiation control is independent of cellular and chromosomal oriC concentration. Cells containing just the ectopic origin initiated bidirectional replication at the expected cell mass and at the normal cellular location of that region. In all strains, spatial separation of sister loci adjacent to active origins occurred shortly after their replication, independently of whether replication initiated at the normal origin, the ectopic origin, or both origins.

RecA-dependent replication in the nrdA101(Ts) mutant of Escherichia coli under restrictive conditions

Cells carrying the thermosensitive nrdA101 allele are able to replicate entire chromosomes at 42°C when new DNA initiation events are inhibited. We investigated the role of the recombination enzymes on the progression of the DNA replication forks in the nrdA101 mutant at 42°C in the presence of rifampin. Using pulsed-field gel electrophoresis (PFGE), we demonstrated that the replication forks stalled and reversed during the replication progression under this restrictive condition. DNA labeling and flow cytometry experiments supported this finding as the deleterious effects found in the RecB-deficient background were suppressed specifically by the absence of RuvABC; however, this did not occur in a RecG-deficient background. Furthermore, we show that the RecA protein is absolutely required for DNA replication in the nrdA101 mutant at restrictive temperature when the replication forks are reversed. The detrimental effect of the recA deletion is not related to the chromosomal degradation caused by the absence of RecA. The inhibition of DNA replication observed in the nrdA101 recA mutant at 42°C in the presence of rifampin was reverted by the presence of the wild-type RecA protein expressed ectopically but only partially suppressed by the RecA protein with an S25P mutation [RecA(S25P)], deficient in the rescue of the stalled replication forks. We propose that RecA is required to maintain the integrity of the reversed forks in the nrdA101 mutant under certain restrictive conditions, supporting the relationship between DNA replication and recombination enzymes through the stabilization and repair of the stalled replication forks.

Overlap of replication rounds disturbs the progression of replicating forks in a ribonucleotide reductase mutant of Escherichia coli

Ribonucleotide reductase (RNR) is the only enzyme specifically required for the synthesis of deoxyribonucleotides (dNTPs). Surprisingly, Escherichia coli cells carrying the nrdA101 allele, which codes for a thermosensitive RNR101, are able to replicate entire chromosomes at 42 °C under RNA or protein synthesis inhibition. Here we show that the RNR101 protein is unstable at 42 °C and that its degradation under restrictive conditions is prevented by the presence of rifampicin. Nevertheless, the mere stability of the RNR protein at 42 °C cannot explain the completion of chromosomal DNA replication in the nrdA101 mutant. We found that inactivation of the DnaA protein by using several dnaAts alleles allows complete chromosome replication in the absence of rifampicin and suppresses the nucleoid segregation and cell division defects observed in the nrdA101 mutant at 42 °C. As both inactivation of the DnaA protein and inhibition of RNA synthesis block the occurrence of new DNA initiations, the consequent decrease in the number of forks per chromosome could be related to those effects. In support of this notion, we found that avoiding multifork replication rounds by the presence of moderate extra copies of datA sequence increases the relative amount of DNA synthesis of the nrdA101 mutant at 42 °C. We propose that a lower replication fork density results in an improvement of the progression of DNA replication, allowing replication of the entire chromosome at the restrictive temperature. The mechanism related to this effect is also discussed.

DNA replication defects in a mutant deficient in the thioredoxin homolog YbbN

Escherichia coli contains two thioredoxins, Trx1 and Trx2, and a thioredoxin-like protein, YbbN, that displays both redox and chaperone properties. Since three out of the six proteins of the YbbN interactome (Butland et al., 2005) are components of DNA polymerase 3 holoenzyme (i.e. the β-clamp DnaN, the θ subunit HolE and the δ' subunit HolB), we investigated whether the ybbN mutant presents DNA replication defects. We found that this mutant incorporates (3)H-thymidine at higher rates than the parental strain and displays overinitiation, hypermutator and filamentation phenotypes with the occurrence of anucleated cells. Moreover, YbbN functions as a bona fide chaperone in the refolding of the urea-unfolded β-clamp. These results suggest that the DNA replication and cell division defects of the ybbN mutant might best be explained by chaperone functions of YbbN in the biogenesis of DNA polymerase 3 holoenzyme.

Replication patterns and organization of replication forks in Vibrio cholerae

We have investigated the replication patterns of the two chromosomes of the bacterium Vibrio cholerae grown in four different media. By combining flow cytometry and quantitative real-time PCR with computer simulations, we show that in rich media, V. cholerae cells grow with overlapping replication cycles of both the large chromosome (ChrI) and the small chromosome (ChrII). In Luria-Bertani (LB) medium, initiation occurs at four copies of the ChrI origin and two copies of the ChrII origin. Replication of ChrII was found to occur at the end of the ChrI replication period in all four growth conditions. Novel cell-sorting experiments with marker frequency analysis support these conclusions. Incubation with protein synthesis inhibitors indicated that the potential for initiation of replication of ChrII was present at the same time as that of ChrI, but was actively delayed until much of ChrI was replicated. Investigations of the localization of SeqA bound to new DNA at replication forks indicated that the forks were co-localized in pairs when cells grew without overlapping replication cycles and in higher-order structures during more rapid growth. The increased degree of fork organization during rapid growth may be a means by which correct segregation of daughter molecules is facilitated.

Interplay of DNA repair, homologous recombination, and DNA polymerases in resistance to the DNA damaging agent 4-nitroquinoline-1-oxide in Escherichia coli

Escherichia coli has three DNA damage-inducible DNA polymerases: DNA polymerase II (Pol II), DNA polymerase IV (Pol IV), and DNA polymerase V (Pol V). While the in vivo function of Pol V is well understood, the precise roles of Pol IV and Pol II in DNA replication and repair are not as clear. Study of these polymerases has largely focused on their participation in the recovery of failed replication forks, translesion DNA synthesis, and origin-independent DNA replication. However, their roles in other repair and recombination pathways in E. coli have not been extensively examined. This study investigated how E. coli's inducible DNA polymerases and various DNA repair and recombination pathways function together to convey resistance to 4-nitroquinoline-1-oxide (NQO), a DNA damaging agent that produces replication blocking DNA base adducts. The data suggest that full resistance to this compound depends upon an intricate interplay among the activities of the inducible DNA polymerases and recombination. The data also suggest new relationships between the different pathways that process recombination intermediates.

DnaA-ATP acts as a molecular switch to control levels of ribonucleotide reductase expression in Escherichia coli

Ribonucleotide reductase (RNR) is the bottleneck enzyme in the synthesis of dNTPs required for DNA replication. In order to avoid the mutagenic effects of imbalances in dNTPs the amount and activity of RNR enzyme in the cell is tightly regulated. RNR expression from the nrdAB operon is thus coupled to coincide with the initiation of DNA replication. However, the mechanism for the co-ordination of gene transcription and DNA replication remains to be elucidated. The timing and synchrony of DNA replication initiation in Escherichia coli is controlled in part by the binding of the DnaA protein to the origin of replication. DnaA is also a transcription factor of the nrdAB operon and could thus be the link between these two processes. Here we show that RNA polymerase can form a stable transcription initiation complex at the nrdAB promoter by direct interaction with the far upstream sites required for the timing of expression as a function of DNA replication. In addition, we show that the binding of DnaA on the promoter can either activate or repress transcription as a function of its concentration and its nucleotide-bound state. However, transcription regulation by DnaA does not significantly affect the timing of expression of RNR from the nrdAB operon.

A reduction in ribonucleotide reductase activity slows down the chromosome replication fork but does not change its localization

Background

It has been proposed that the enzymes of nucleotide biosynthesis may be compartmentalized or concentrated in a structure affecting the organization of newly replicated DNA. Here we have investigated the effect of changes in ribonucleotide reductase (RNR) activity on chromosome replication and organization of replication forks in Escherichia coli.

Methodology/Principal Findings

Reduced concentrations of deoxyribonucleotides (dNTPs) obtained by reducing the activity of wild type RNR by treatment with hydroxyurea or by mutation, resulted in a lengthening of the replication period. The replication fork speed was found to be gradually reduced proportionately to moderate reductions in nucleotide availability. Cells with highly extended C periods showed a “delay” in cell division i.e. had a higher cell mass. Visualization of SeqA structures by immunofluorescence indicated no change in organization of the new DNA upon moderate limitation of RNR activity. Severe nucleotide limitation led to replication fork stalling and reversal. Well defined SeqA structures were not found in situations of extensive replication fork repair. In cells with stalled forks obtained by UV irradiation, considerable DNA compaction was observed, possibly indicating a reorganization of the DNA into a “repair structure” during the initial phase of the SOS response.

Conclusion/Significance

The results indicate that the replication fork is slowed down in a controlled manner during moderate nucleotide depletion and that a change in the activity of RNR does not lead to a change in the organization of newly replicated DNA. Control of cell division but not control of initiation was affected by the changes in replication elongation.

Growth rate dependent numbers of SeqA structures organize the multiple replication forks in rapidly growing Escherichia coli

When the bacterium Escherichia coli is grown in rich medium, the replication and segregation periods may span two, three or four generations and cells may contain up to 24 replication forks. The newly synthesized, hemimethylated DNA at each fork is bound by SeqA protein. The SeqA-DNA structures form distinct foci that can be observed by immunofluorescence microscopy. The numbers of foci were lower than the numbers of replication forks indicating fork co-localization. The extent of co-localization correlated with the extent of replication cycle overlap in wild-type cells. No abrupt increase in the numbers of foci occurred at the time of initiation of replication, suggesting that new replication forks bind to existing SeqA structures. Manipulations with replication control mechanisms that led to extension or reduction of the replication period and number of forks, did not lead to changes in the numbers of SeqA foci per cell. The results indicate that the number of SeqA foci is not directly governed by the number of replication forks, and supports the idea that new DNA may be 'captured' by existing SeqA structures.

Plasmid DNA vaccine vector design: impact on efficacy, safety and upstream production

Critical molecular and cellular biological factors impacting design of licensable DNA vaccine vectors that combine high yield and integrity during bacterial production with increased expression in mammalian cells are reviewed. Food and Drug Administration (FDA), World Health Organization (WHO) and European Medical Agencies (EMEA) regulatory guidance's are discussed, as they relate to vector design and plasmid fermentation. While all new vectors will require extensive preclinical testing to validate safety and performance prior to clinical use, regulatory testing burden for follow-on products can be reduced by combining carefully designed synthetic genes with existing validated vector backbones. A flowchart for creation of new synthetic genes, combining rationale design with bioinformatics, is presented. The biology of plasmid replication is reviewed, and process engineering strategies that reduce metabolic burden discussed. Utilizing recently developed low metabolic burden seed stock and fermentation strategies, optimized vectors can now be manufactured in high yields exceeding 2 g/L, with specific plasmid yields of 5% total dry cell weight.

The initiator protein DnaA contributes to keeping new origins inactivated by promoting the presence of hemimethylated DNA

The Escherichia coli replication origin oriC and other regions with high numbers of GATC sites remain hemimethylated after replication much longer than regions with average numbers of GATC sites. The prolonged period of hemimethylation has been attributed to the presence of bound SeqA protein. Here, it was found that a GATC cluster inserted at the datA site, which binds large amounts of DnaA in vivo, did not become remethylated at all, unless the availability of the DnaA protein was severely reduced. Sequestration of oriC was also found to be affected by the availability of DnaA. The period of origin hemimethylation was reduced by approximately 30% upon a reduction in the availability of DnaA. The result shows that not only SeqA binding but also DnaA binding to newly replicated origins contributes to keeping them hemimethylated. It was also found that the number of SeqA foci increased in cells with a combination of DnaA-mediated protection and sequestration at the GATC::datA cluster.

New model for assembly dynamics of bacterial tubulin in relation to the stages of DNA replication

How living cells receive their genome through cell division has been one of the important questions of biology. In prokaryotes, cell division starts with formation of a ring-shaped microtubule-like structure, FtsZ-ring, at the potential division site. All the previous models suggested that FtsZ-ring is formed coupling to termination or far after initiation of DNA replication. In contrast, we demonstrated that a close communication with DNA replication is maintained throughout the cell cycle. FtsZ starts to assemble to the cell center coupling to initiation of DNA replication, and stabilizes as FtsZ-ring at its termination, but does not constrict before separation of nucleoids. This combination of a positive and a negative control would guarantee that a successful replication event would inevitably induce one cell division such that each of the daughter cells would receive one and only one daughter nucleoid.

Chromosomal model for analysis of a long CTG/CAG tract stability in wild-type Escherichia coli and its nucleotide excision repair mutants

Many human hereditary neurological diseases, including fragile X syndrome, myotonic dystrophy, and Friedreich's ataxia, are associated with expansions of the triplet repeat sequences (TRS) (CGG/CCG, CTG/CAG, and GAA/TTC) within or near specific genes. Mechanisms that mediate mutations of TRS include DNA replication, repair, and gene conversion and (or) recombination. The involvement of the repair systems in TRS instability was investigated in Escherichia coli on plasmid models, and the results showed that the deficiency of some nucleotide excision repair (NER) functions dramatically affects the stability of long CTG inserts. In such models in which there are tens or hundreds of plasmid molecules in each bacterial cell, repetitive sequences may interact between themselves and according to a recombination hypothesis, which may lead to expansions and deletions within such repeated tracts. Since one cannot control interaction between plasmids, it is also sometimes difficult to give precise interpretation of the results. Therefore, using modified lambda phage (lambdaInCh), we have constructed a chromosomal model to study the instability of trinucleotide repeat sequences in E. coli. We have shown that the stability of (CTG/CAG)68 tracts in the bacterial chromosome is influenced by mutations in NER genes in E. coli. The absence of the uvrC or uvrD gene products greatly enhances the instability of the TRS in the chromosome, whereas the lack of the functional UvrA or UvrB proteins causes substantial stabilization of (CTG/CAG) tracts.

Organization of sister origins and replisomes during multifork DNA replication in Escherichia coli

The replication period of Escherichia coli cells grown in rich medium lasts longer than one generation. Initiation thus occurs in the 'mother-' or 'grandmother generation'. Sister origins in such cells were found to be colocalized for an entire generation or more, whereas sister origins in slow-growing cells were colocalized for about 0.1-0.2 generations. The role of origin inactivation (sequestration) by the SeqA protein in origin colocalization was studied by comparing sequestration-deficient mutants with wild-type cells. Cells with mutant, non-sequesterable origins showed wild-type colocalization of sister origins. In contrast, cells unable to sequester new origins due to loss of SeqA, showed aberrant localization of origins indicating a lack of organization of new origins. In these cells, aberrant replisome organization was also found. These results suggest that correct organization of sister origins and sister replisomes is dependent on the binding of SeqA protein to newly formed DNA at the replication forks, but independent of origin sequestration. In agreement, in vitro experiments indicate that SeqA is capable of pairing newly replicated DNA molecules.

The two chromosomes of Vibrio cholerae are initiated at different time points in the cell cycle

The bacterium Vibrio cholerae, the cause of the diarrhoeal disease cholera, has its genome divided between two chromosomes, a feature uncommon for bacteria. The two chromosomes are of different sizes and different initiator molecules control their replication independently. Using novel methods for analysing flow cytometry data and marker frequency analysis, we show that the small chromosome II is replicated late in the C period of the cell cycle, where most of chromosome I has been replicated. Owing to the delay in initiation of chromosome II, the two chromosomes terminate replication at approximately the same time and the average number of replication origins per cell is higher for chromosome I than for chromosome II. Analysis of cell-cycle parameters shows that chromosome replication and segregation is exceptionally fast in V. cholerae. The divided genome and delayed replication of chromosome II may reduce the metabolic burden and complexity of chromosome replication by postponing DNA synthesis to the last part of the cell cycle and reducing the need for overlapping replication cycles during rapid proliferation.

The size and form of chromosomes are constant in the nucleus, but highly variable in bacteria, mitochondria and chloroplasts

From cytological examination, the size and form of the chromosomes in the eukaryotic nucleus are invariant across generations, leading to the expectation that constancy of inheritance likely depends on constancy of the chromosomal DNA molecule conveying the constant phenotype. Indeed, except for rare mutations, major phenotypic traits appear largely without change from generation to generation. Thus, when it was discovered that the inheritance of traits for bacteria, mitochondria and chloroplasts was also constant, it was assumed that chromosomes in those locations were also constant. Such has not turned out to be the case, however; those chromosomes are highly variable in structure. I propose, therefore, that only for the nucleus is there a requirement that a chromosome be "finished" (contain only fully replicated genomes) before it may segregate to daughter cells. This requirement does not apply to the variable chromosomes among chloroplasts, mitochondria and bacteria.

Defective ribonucleoside diphosphate reductase impairs replication fork progression in Escherichia coli

The observed lengthening of the C period in the presence of a defective ribonucleoside diphosphate reductase has been assumed to be due solely to the low deoxyribonucleotide supply in the nrdA101 mutant strain. We show here that the nrdA101 mutation induces DNA double-strand breaks at the permissive temperature in a recB-deficient background, suggesting an increase in the number of stalled replication forks that could account for the slowing of replication fork progression observed in the nrdA101 strain in a Rec(+) context. These DNA double-strand breaks require the presence of the Holliday junction resolvase RuvABC, indicating that they have been generated from stalled replication forks that were processed by the specific reaction named "replication fork reversal." Viability results supported the occurrence of this process, as specific lethality was observed in the nrdA101 recB double mutant and was suppressed by the additional inactivation of ruvABC. None of these effects seem to be due to the limitation of the deoxyribonucleotide supply in the nrdA101 strain even at the permissive temperature, as we found the same level of DNA double-strand breaks in the nrdA(+) strain growing under limited (2-microg/ml) or under optimal (5-microg/ml) thymidine concentrations. We propose that the presence of an altered NDP reductase, as a component of the replication machinery, impairs the progression of the replication fork, contributing to the lengthening of the C period in the nrdA101 mutant at the permissive temperature.

Sequence characteristics required for cooperative binding and efficient in vivo titration of the replication initiator protein DnaA in E. coli

Plasmids carrying the mioC promoter region, which contains two DnaA boxes, R5 and R6 with one misfit to the consensus TT(A)/(T)TNCACA, are as efficient in in vivo titration of the DnaA protein as plasmids carrying a replication-inactivated oriC region with its eight DnaA boxes. Three additional DnaA boxes around the promoter proximal R5 DnaA box were identified and shown by mutational analysis to be necessary for the cooperative binding of DnaA required for titration. These four DnaA boxes are located in the same orientation and with a spacing of two or three base-pairs. The cooperative binding was eliminated by insertion of half a helical turn between any of the DnaA boxes. Titration strongly depends on the presence and orientation of the promoter distal R6 DnaA box located 104 bp upstream of the R5 box as well as neighbouring sequences downstream of R6. Titration depends on the integrity of a 43 bp region containing the R5 DnaA box, while repression of mioC transcription by DnaA, which is dependent on the R5 DnaA box, was independent of the two DnaA boxes downstream of R5. Repression was also independent of the spacing between the two upstream DnaA boxes and the promoter as long as a DnaA box was located less than 20 bp upstream of the -35 sequence. Thus, the architectural requirements for titration and for repression of transcription are different. A new set of rules for identifying efficiently titrating DnaA box regions was formulated and used to analyse sequences for which good titration data are available.

Differences in the degree of inhibition of NDP reductase by chemical inactivation and by the thermosensitive mutation nrdA101 in Escherichia coli suggest an effect on chromosome segregation

NDP reductase activity can be inhibited either by treatment with hydroxyurea or by incubation of an nrdA_ts mutant strain at the non-permissive temperature. Both methods inhibit replication, but experiments on these two types of inhibition yielded very different results. The chemical treatment immediately inhibited DNA synthesis but did not affect the cell and nucleoid appearance, while the incubation of an nrdA101 mutant strain at the non-permissive temperature inhibited DNA synthesis after more than 50 min, and resulted in aberrant chromosome segregation, long filaments, and a high frequency of anucleate cells. These phenotypes are not induced by SOS. In view of these results, we suggest there is an indirect relationship between NDP reductase and the chromosome segregation machinery through the maintenance of the proposed replication hyperstructure.

Tracking of controlled Escherichia coli replication fork stalling and restart at repressor-bound DNA in vivo

We report an efficient, controllable, site-specific replication roadblock that blocks cell proliferation, but which can be rapidly and efficiently reversed, leading to recovery of viability. Escherichia coli replication forks of both polarities stalled in vivo within the first 500 bp of a 10 kb repressor-bound array of operator DNA-binding sites. Controlled release of repressor binding led to rapid restart of the blocked replication fork without the participation of homologous recombination. Cytological tracking of fork stalling and restart showed that the replisome-associated SSB protein remains associated with the blocked fork for extended periods and that duplication of the fluorescent foci associated with the blocked operator array occurs immediately after restart, thereby demonstrating a lack of sister cohesion in the region of the array. Roadblocks positioned near oriC or the dif site did not prevent replication and segregation of the rest of the chromosome.

Generation of DNA-free Escherichia coli cells by 2-aminopurine requires mismatch repair and nonmethylated DNA

Undirected mismatch repair initiated by the incorporation of the base analog 2-aminopurine kills DNA-methylation-deficient Escherichia coli dam cells by DNA double-strand breakage. Subsequently, the chromosomal DNA is totally degraded, resulting in DNA-free cells.

Dancing around the divisome: asymmetric chromosome segregation in Escherichia coli

By simultaneously tracking pairs of specific genetic regions and divisome proteins in live Escherichia coli, we develop a new scheme for the relationship between DNA replication-segregation, chromosome organization, and cell division. A remarkable asymmetric pattern of segregation of different loci in the replication termination region (ter) suggests that individual replichores segregate to distinct nucleoid positions, consistent with an asymmetric segregation of leading and lagging strand templates after replication. Cells growing with a generation time of 100 min are born with a nonreplicating chromosome and have their origin region close to mid-cell and their ter polar. After replication initiation, the two newly replicated origin regions move away from mid-cell to opposite cell halves. By mid-S phase, FtsZ forms a ring at mid-cell at the time of initiation of nucleoid separation; ter remains polar. In the latter half of S phase, ter moves quickly toward mid-cell. FtsK, which coordinates the late stages of chromosome segregation with cell division, forms a ring coincident with the FtsZ ring as S phase completes, approximately 50 min after its initiation. As ter duplicates at mid-cell, sister nucleoid separation appears complete. After initiation of invagination, the FtsZ ring disassembles, leaving FtsK to complete chromosome segregation and cytokinesis.

Localization of replication forks in wild-type and mukB mutant cells of Escherichia coli

To examine the subcellular localization of the replication machinery in Escherichia coli, we have developed an immunofluorescence method that allows us to determine the subcellular location of newly synthesized DNA pulse-labeled with 5-bromo-2'-deoxyuridine (BrdU). Using this technique, we have analyzed growing cells. In wild-type cells that showed a single BrdU fluorescence signal, the focus was located in the middle of the cell; in cells with two signals, the foci were localized at positions equivalent to 1/4 and 3/4 of the cell length. The formation of BrdU foci was dependent upon ongoing chromosomal replication. A mutant lacking MukB, which is required for proper partitioning of sister chromosomes, failed to maintain the ordered localization of BrdU foci: (1) a single BrdU focus tended to be localized at a pole-proximal region of the nucleoid, and (2) a focus was often found to consist of two replicating chromosomes. Thus, the positioning of replication forks is affected by the disruption of the mukB gene.