Validation of bacterial replication termination models using simulation of genomic mutations

Interplay between chromosomal architecture and termination of DNA replication in bacteria

Faithful transmission of the genome from one generation to the next is key to life in all cellular organisms. In the majority of bacteria, the genome is comprised of a single circular chromosome that is normally replicated from a single origin, though additional genetic information may be encoded within much smaller extrachromosomal elements called plasmids. By contrast, the genome of a eukaryote is distributed across multiple linear chromosomes, each of which is replicated from multiple origins. The genomes of archaeal species are circular, but are predominantly replicated from multiple origins. In all three cases, replication is bidirectional and terminates when converging replication fork complexes merge and 'fuse' as replication of the chromosomal DNA is completed. While the mechanics of replication initiation are quite well understood, exactly what happens during termination is far from clear, although studies in bacterial and eukaryotic models over recent years have started to provide some insight. Bacterial models with a circular chromosome and a single bidirectional origin offer the distinct advantage that there is normally just one fusion event between two replication fork complexes as synthesis terminates. Moreover, whereas termination of replication appears to happen in many bacteria wherever forks happen to meet, termination in some bacterial species, including the well-studied bacteria Escherichia coli and Bacillus subtilis, is more restrictive and confined to a 'replication fork trap' region, making termination even more tractable. This region is defined by multiple genomic terminator (ter) sites, which, if bound by specific terminator proteins, form unidirectional fork barriers. In this review we discuss a range of experimental results highlighting how the fork fusion process can trigger significant pathologies that interfere with the successful conclusion of DNA replication, how these pathologies might be resolved in bacteria without a fork trap system and how the acquisition of a fork trap might have provided an alternative and cleaner solution, thus explaining why in bacterial species that have acquired a fork trap system, this system is remarkably well maintained. Finally, we consider how eukaryotic cells can cope with a much-increased number of termination events.

A soft Tus-Ter interaction is hiding a fail-safe lock in the replication fork trap of Dickeya paradisiaca

A variety of replication fork traps have recently been characterised in Enterobacterales, unveiling two different types of architecture. Of these, the degenerate type II fork traps are commonly found in Enterobacteriaceae such as Escherichia coli. The newly characterised type I fork traps are found almost exclusively outside Enterobacteriaceae within Enterobacterales and include several archetypes of possible ancestral architectures. Dickeya paradisiaca harbours a somewhat degenerate type I fork trap with a unique Ter1 adjacent to tus gene on one side of the circular chromosome and three putative Ter2-4 sites on the other side of the fork trap. The two innermost Ter1 and Ter2 sites are only separated by 18 kb, which is the shortest distance between two innermost Ter sites of any chromosomal fork trap identified so far. Of note, the dif site is located between these two sites, coinciding with a sharp GC-skew flip. Here we examined and compared the binding modalities of E. coli and D. paradisiaca Tus proteins for these Ter sites. Surprisingly, while Ter1-3 were functional, no significant Tus binding was observed for Ter4 even in low salt conditions, which is in stark contrast with the significant non-specific protein-DNA interactions that occur with E. coli Tus. Even more surprising was the finding that D. paradisiaca Tus has a relatively moderate binding affinity to double-stranded Ter while retaining an extremely high affinity to Ter-lock sequences. Our data revealed major differences in the salt resistance and stability between the D. paradisiaca and E. coli Tus protein complexes, suggesting that while Tus protein evolution can be quite flexible regarding the initial Ter binding step, it requires a highly stringent purifying selection for its final locked complex formation.

Delineation of the Ancestral Tus-Dependent Replication Fork Trap

In Escherichia coli, DNA replication termination is orchestrated by two clusters of Ter sites forming a DNA replication fork trap when bound by Tus proteins. The formation of a ‘locked’ Tus–Ter complex is essential for halting incoming DNA replication forks. However, the absence of replication fork arrest at some Ter sites raised questions about their significance. In this study, we examined the genome-wide distribution of Tus and found that only the six innermost Ter sites (TerA–E and G) were significantly bound by Tus. We also found that a single ectopic insertion of TerB in its non-permissive orientation could not be achieved, advocating against a need for ‘back-up’ Ter sites. Finally, examination of the genomes of a variety of Enterobacterales revealed a new replication fork trap architecture mostly found outside the Enterobacteriaceae family. Taken together, our data enabled the delineation of a narrow ancestral Tus-dependent DNA replication fork trap consisting of only two Ter sites.

A Fork Trap in the Chromosomal Termination Area Is Highly Conserved across All Escherichia coli Phylogenetic Groups

Termination of DNA replication, the final stage of genome duplication, is surprisingly complex, and failures to bring DNA synthesis to an accurate conclusion can impact genome stability and cell viability. In Escherichia coli, termination takes place in a specialised termination area opposite the origin. A 'replication fork trap' is formed by unidirectional fork barriers via the binding of Tus protein to genomic ter sites. Such a fork trap system is found in some bacterial species, but it appears not to be a general feature of bacterial chromosomes. The biochemical properties of fork trap systems have been extensively characterised, but little is known about their precise physiological roles. In this study, we compare locations and distributions of ter terminator sites in E. coli genomes across all phylogenetic groups, including Shigella. Our analysis shows that all ter sites are highly conserved in E. coli, with slightly more variability in the Shigella genomes. Our sequence analysis of ter sites and Tus proteins shows that the fork trap is likely to be active in all strains investigated. In addition, our analysis shows that the dif chromosome dimer resolution site is consistently located between the innermost ter sites, even if rearrangements have changed the location of the innermost termination area. Our data further support the idea that the replication fork trap has an important physiological role that provides an evolutionary advantage.

Genome-wide analysis of the synonymous codon usage pattern of Streptococcus suis

Streptococcus suis (S. suis) is a gram-positive coccus that causes disease in humans and animals. The codon usage pattern of bacteria reveals a range of evolutionary changes that assist them to enhance tolerance to environments. To better understand the genetic features during the evolution of S. suis, we performed codon usage analysis. Nine pathogenic strains of different serotypes and different geographical distribution were analyzed to better understand the differences in their evolutionary process. Nucleotide compositions and relative synonymous codon usage (RSCU) analysis revealed that A/T-ending codons are dominant in S. suis. Neutrality analysis, correspondence analysis and ENC-plot results revealed that natural selection is the predominant element prompting codon usage. Cluster analysis based on RSCU was roughly consistent with the dendrogram rooted genomic BLAST analysis. Comparison of synonymous codon usage pattern between S. suis and susceptible hosts (H. sapiens and S. scrofa) revealed that the codon usage of S. suis is separated from the synonymous codon usage of susceptible hosts. The CAI values implied that S. suis includes a series of predicted highly expressed coding sequences contained in metabolism and transcriptional regulation, revealing the necessity of this pathogen to deal with various environmental conditions. The study of codon usage in S. suis may provide evidence involving the molecular evolution of bacteria and a better understanding of evolutionary relationships between S. suis and its corresponding hosts.

Accelerated Laboratory Evolution Reveals the Influence of Replication on the GC Skew in Escherichia coli

Most bacterial genomes display contrasting strand asymmetry in a variety of features, such as nucleotide composition and gene orientation, of the two replichores separated by the replication origin and terminus. The cause for the polarization is often attributed to mutations arising from the asymmetric replication machinery. Notably, a base compositional bias known as a GC skew is focused on as a footprint of the bacterial genome evolution driven by DNA replication. Previously, although a replication driven mutation pattern responsible for the GC skew formation or the related mathematical models have been well reported, an exact impact of the replication-related elements on the genomic structure is yet actively debated, and not confirmed experimentally. However, the GC skew formation is very time consuming and challenging in the laboratory. We, therefore, used cytosine deaminase as a DNA mutator, and by monitoring the mutations during an accelerated laboratory evolution procedure with Illumina sequencing, we enabled the trial and error of the GC skew formation in high resolution. Using this technology, we succeeded in reconfirming the influence of bacterial replication machinery on the genomic structure at high resolution.

Segregation but Not Replication of the Pseudomonas aeruginosa Chromosome Terminates at Dif

Coordination between chromosome replication and segregation is essential for equal partitioning of genetic material between daughter cells. In bacteria, this is achieved through the proximity of the origin of replication, oriC, and the chromosome partitioning site, parS We report here that in Pseudomonas aeruginosa, segregation but not replication is also controlled at the terminus region of the chromosome. Using the fluorescent repressor operator system (FROS), we investigated chromosome segregation in P. aeruginosa strain PAO1-UW, wherein the chromosome dimer resolution site, dif, is asymmetrically positioned relative to oriC In these cells, segregation proceeded sequentially along the two chromosomal arms and terminated at dif In contrast, chromosome replication terminated elsewhere, opposite from oriC We further found two large domains on the longer arm of the chromosome, wherein DNA segregated simultaneously. Notably, GC-skew, which reflects a bias in nucleotide usage between the leading and lagging strands of the chromosome, switches polarity at the dif locus but not necessarily at the terminus of replication. These data demonstrate that termination of chromosome replication and segregation can be physically separated without adverse effects on bacterial fitness. They also reveal the critical role of the dif region in defining the global layout of the chromosome and the progression of chromosome segregation and suggest that chromosome packing adapts to its subcellular layout.IMPORTANCE Segregation of genetic information is a central event in cellular life. In bacteria, chromosome segregation occurs concurrently with replication, sequentially along the two arms from oriC to dif How the two processes are coordinated is unknown. We explored here chromosome segregation in an opportunistic human pathogen, Pseudomonas aeruginosa, using its strain with markedly unequal chromosomal arms. We found that replication and segregation diverge in this strain and terminate at very different locations, whereas the longer chromosomal arm folds into large domains to align itself with the shorter arm. The significance of this research is in establishing that segregation and replication of bacterial chromosomes are largely uncoupled from each other and that the large-scale structure of the chromosome adapts to its subcellular layout.

Comparison of Intracellular "Ca. Endomicrobium Trichonymphae" Genomovars Illuminates the Requirement and Decay of Defense Systems against Foreign DNA

“Candidatus Endomicrobium trichonymphae” (Bacteria; Elusimicrobia) is an obligate intracellular symbiont of the cellulolytic protist genus Trichonympha in the termite gut. A previous genome analysis of “Ca. Endomicrobium trichonymphae” phylotype Rs-D17 (genomovar Ri2008), obtained from a Trichonympha agilis cell in the gut of the termite Reticulitermes speratus, revealed that its genome is small (1.1 Mb) and contains many pseudogenes; it is in the course of reductive genome evolution. Here we report the complete genome sequence of another Rs-D17 genomovar, Ti2015, obtained from a different T. agilis cell present in an R. speratus gut. These two genomovars share most intact protein-coding genes and pseudogenes, showing 98.6% chromosome sequence similarity. However, characteristic differences were found in their defense systems, which comprised restriction-modification and CRISPR/Cas systems. The repertoire of intact restriction-modification systems differed between the genomovars, and two of the three CRISPR/Cas loci in genomovar Ri2008 are pseudogenized or missing in genomovar Ti2015. These results suggest relaxed selection pressure for maintaining these defense systems. Nevertheless, the remaining CRISPR/Cas system in each genomovar appears to be active; none of the “spacer” sequences (112 in Ri2008 and 128 in Ti2015) were shared whereas the “repeat” sequences were identical. Furthermore, we obtained draft genomes of three additional endosymbiotic Endomicrobium phylotypes from different host protist species, and discovered multiple, intact CRISPR/Cas systems in each genome. Collectively, unlike bacteriome endosymbionts in insects, the Endomicrobium endosymbionts of termite-gut protists appear to require defense against foreign DNA, although the required level of defense has likely been reduced during their intracellular lives.

Replication Termination: Containing Fork Fusion-Mediated Pathologies in Escherichia coli

Duplication of bacterial chromosomes is initiated via the assembly of two replication forks at a single defined origin. Forks proceed bi-directionally until they fuse in a specialised termination area opposite the origin. This area is flanked by polar replication fork pause sites that allow forks to enter but not to leave. The precise function of this replication fork trap has remained enigmatic, as no obvious phenotypes have been associated with its inactivation. However, the fork trap becomes a serious problem to cells if the second fork is stalled at an impediment, as replication cannot be completed, suggesting that a significant evolutionary advantage for maintaining this chromosomal arrangement must exist. Recently, we demonstrated that head-on fusion of replication forks can trigger over-replication of the chromosome. This over-replication is normally prevented by a number of proteins including RecG helicase and 3’ exonucleases. However, even in the absence of these proteins it can be safely contained within the replication fork trap, highlighting that multiple systems might be involved in coordinating replication fork fusions. Here, we discuss whether considering the problems associated with head-on replication fork fusion events helps us to better understand the important role of the replication fork trap in cellular metabolism.

Quasispecies in population of compositional assemblies

BACKGROUND

The quasispecies model refers to information carriers that undergo self-replication with errors. A quasispecies is a steady-state population of biopolymer sequence variants generated by mutations from a master sequence. A quasispecies error threshold is a minimal replication accuracy below which the population structure breaks down. Theory and experimentation of this model often refer to biopolymers, e.g. RNA molecules or viral genomes, while its prebiotic context is often associated with an RNA world scenario. Here, we study the possibility that compositional entities which code for compositional information, intrinsically different from biopolymers coding for sequential information, could show quasispecies dynamics.

RESULTS

We employed a chemistry-based model, graded autocatalysis replication domain (GARD), which simulates the network dynamics within compositional molecular assemblies. In GARD, a compotype represents a population of similar assemblies that constitute a quasi-stationary state in compositional space. A compotype's center-of-mass is found to be analogous to a master sequence for a sequential quasispecies. Using single-cycle GARD dynamics, we measured the quasispecies transition matrix (Q) for the probabilities of transition from one center-of-mass Euclidean distance to another. Similarly, the quasispecies' growth rate vector (A) was obtained. This allowed computing a steady state distribution of distances to the center of mass, as derived from the quasispecies equation. In parallel, a steady state distribution was obtained via the GARD equation kinetics. Rewardingly, a significant correlation was observed between the distributions obtained by these two methods. This was only seen for distances to the compotype center-of-mass, and not to randomly selected compositions. A similar correspondence was found when comparing the quasispecies time dependent dynamics towards steady state. Further, changing the error rate by modifying basal assembly joining rate of GARD kinetics was found to display an error catastrophe, similar to the standard quasispecies model. Additional augmentation of compositional mutations leads to the complete disappearance of the master-like composition.

CONCLUSIONS

Our results show that compositional assemblies, as simulated by the GARD formalism, portray significant attributes of quasispecies dynamics. This expands the applicability of the quasispecies model beyond sequence-based entities, and potentially enhances validity of GARD as a model for prebiotic evolution.

Undesigned selection for replication termination of bacterial chromosomes

The oriC DNA replication origin in bacterial chromosomes, the location of which appears to be physically identified, is genetically regulated by relevant molecular machinery. In contrast, the location of the terminus remains obscure for many bacterial replicons, except for terC, the proposed and well-studied chromosome termination site in certain bacteria. The terC locus, which is composed of specific sequences for its binding protein, is located at a site opposite from oriC, exhibiting a symmetric structure around the oriC-terC axis. Here, we investigated Bacillus subtilis 168 strains whose axes were hindered and found that the native terC function was robust. However, eradication of terminus region specific binding protein resulted in the natural terC sites not being used for termination; instead, new termini were selected at a location exactly opposite to oriC. We concluded that replication generally terminates at the loci where the two approaching replisomes meet. This site was automatically selected, and two replisomes moving at the same rate supported symmetrical chromosome structures relative to oriC. The rule, which was even validated by artificial chromosomes irrespective of oriC, should be general for replicons administered by two replisomes.