Chromosome structuring limits genome plasticity in Escherichia coli

A bipartite interaction with the processivity clamp potentiates Pol IV-mediated TLS

Processivity clamps mediate polymerase switching for translesion synthesis (TLS). All three E. coli TLS polymerases interact with the β2 processivity clamp through a conserved clamp-binding motif (CBM), which is indispensable for TLS. Notably, Pol IV also makes a unique secondary contact with the clamp through non-CBM residues. However, the role of this "rim contact" in Pol IV-mediated TLS remains poorly understood. Here we show that the rim contact is critical for TLS past strong replication blocks. In in vitro reconstituted Pol IV-mediated TLS, ablating the rim contact compromises TLS past 3-methyl dA, a strong block, while barely affecting TLS past N2-furfuryl dG, a weak block. Similar observations are also made in E. coli cells bearing a single copy of these lesions in the genome. Within lesion-stalled replication forks, the rim interaction and ssDNA binding protein cooperatively poise Pol IV to better compete with Pol III for binding to a cleft through its CBM. We propose that this bipartite clamp interaction enables Pol IV to rapidly resolve lesion-stalled replication through TLS at the fork, which reduces damage induced mutagenesis.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Lethal and mutagenic effects of different LET radiations on Bacillus subtilis spores

New, useful microorganism resources have been generated by ionizing radiation breeding technology. However, the mutagenic effects of ionizing radiation on microorganisms have not been systematically clarified. For a deeper understanding and characterization of ionizing radiation-induced mutations in microorganisms, we investigated the lethal effects of seven different linear energy transfer (LET) radiations based on the survival fraction (SF) and whole-genome sequencing analysis of the mutagenic effects of a dose resulting in an SF of around 1% in Bacillus subtilis spores. Consequently, the lower LET radiations (gamma [surface LET: 0.2 keV/µm] and 4He2+ [24 keV/µm]) showed low lethality and high mutation frequency (MF), resulting in the major induction of single-base substitutions. Whereas higher LET radiations (12C5+ [156 keV/µm] and 12C6+ [179 keV/µm]) showed high lethality and low MF, resulting in the preferential induction of deletion mutations. In addition, 12C6+ (111) ion beams likely possess characteristics of both low- and high-LET radiations simultaneously. A decrease in the relative biological effectiveness and an evaluation of the inactivation cross section indicated that 20Ne8+ (468 keV/µm) and 40Ar13+ (2214 keV/µm) ion beams had overkill effects. In conclusion, in the mutation breeding of microorganisms, it should be possible to regulate the proportions, types, and frequencies of induced mutations by selecting an ionizing radiation of an appropriate LET in accordance with the intended purpose.

Chromosomal Position of Ribosomal Protein Genes Affects Long-Term Evolution of Vibrio cholerae

It is unclear how gene order within the chromosome influences genome evolution. Bacteria cluster transcription and translation genes close to the replication origin (oriC). In Vibrio cholerae, relocation of s10-spc-α locus (S10), the major locus of ribosomal protein genes, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction in growth rate, fitness, and infectivity. To test the long-term impact of this trait, we evolved 12 populations of V. cholerae strains bearing S10 at an oriC-proximal or an oriC-distal location for 1,000 generations. During the first 250 generations, positive selection was the main force driving mutation. After 1,000 generations, we observed more nonadaptative mutations and hypermutator genotypes. Populations fixed inactivating mutations at many genes linked to virulence: flagellum, chemotaxis, biofilm, and quorum sensing. Throughout the experiment, all populations increased their growth rates. However, those bearing S10 close to oriC remained the fittest, indicating that suppressor mutations cannot compensate for the genomic position of the main ribosomal protein locus. Selection and sequencing of the fastest-growing clones allowed us to characterize mutations inactivating, among other sites, flagellum master regulators. Reintroduction of these mutations into the wild-type context led to a ≈10% growth improvement. In conclusion, the genomic location of ribosomal protein genes conditions the evolutionary trajectory of V. cholerae. While genomic content is highly plastic in prokaryotes, gene order is an underestimated factor that conditions cellular physiology and evolution. A lack of suppression enables artificial gene relocation as a tool for genetic circuit reprogramming. IMPORTANCE The bacterial chromosome harbors several entangled processes such as replication, transcription, DNA repair, and segregation. Replication begins bidirectionally at the replication origin (oriC) until the terminal region (ter) organizing the genome along the ori-ter axis gene order along this axis could link genome structure to cell physiology. Fast-growing bacteria cluster translation genes near oriC. In Vibrio cholerae, moving them away was feasible but at the cost of losing fitness and infectivity. Here, we evolved strains harboring ribosomal genes close or far from oriC. Growth rate differences persisted after 1,000 generations. No mutation was able to compensate for the growth defect, showing that ribosomal gene location conditions their evolutionary trajectory. Despite the high plasticity of bacterial genomes, evolution has sculpted gene order to optimize the ecological strategy of the microorganism. We observed growth rate improvement throughout the evolution experiment that occurred at expense of energetically costly processes such the flagellum biosynthesis and virulence-related functions. From the biotechnological point of view, manipulation of gene order enables altering bacterial growth with no escape events.

Enzymatic Supercoiling of Bacterial Chromosomes Facilitates Genome Manipulation

The physical stability of bacterial chromosomes is important for their in vitro manipulation, while genetic stability is important in vivo. However, extracted naked chromosomes in the open circular form are fragile due to nicks and gaps. Using a nick/gap repair and negative supercoiling reaction (named SCR), we first achieved the negative supercoiling of the whole genomes extracted from Escherichia coli and Vibrio natriegens cells. Supercoiled chromosomes of 0.2-4.6 megabase (Mb) were separated by size using a conventional agarose gel electrophoresis and served as DNA size markers. We also achieved the enzymatic replication of 1-2 Mb chromosomes using the reconstituted E. coli replication-cycle reaction (RCR). Electroporation-ready 1 Mb chromosomes were prepared by a modified SCR performed at a low salt concentration (L-SCR) and directly introduced into commercial electrocompetent E. coli cells. Since successful electroporation relies on the genetic stability of a chromosome in cells, genetically stable 1 Mb chromosomes were developed according to a portable chromosome format (PCF). Using physically and genetically stabilized chromosomes, the democratization of genome synthetic biology will be greatly accelerated.

Compartmentalization of the replication fork by single-stranded DNA-binding protein regulates translesion synthesis

Processivity clamps tether DNA polymerases to DNA, allowing their access to the primer-template junction. In addition to DNA replication, DNA polymerases also participate in various genome maintenance activities, including translesion synthesis (TLS). However, owing to the error-prone nature of TLS polymerases, their association with clamps must be tightly regulated. Here we show that fork-associated ssDNA-binding protein (SSB) selectively enriches the bacterial TLS polymerase Pol IV at stalled replication forks. This enrichment enables Pol IV to associate with the processivity clamp and is required for TLS on both the leading and lagging strands. In contrast, clamp-interacting proteins (CLIPs) lacking SSB binding are spatially segregated from the replication fork, minimally interfering with Pol IV-mediated TLS. We propose that stalling-dependent structural changes within clusters of fork-associated SSB establish hierarchical access to the processivity clamp. This mechanism prioritizes a subset of CLIPs with SSB-binding activity and facilitates their exchange at the replication fork.

Role for DNA double strand end-resection activity of RecBCD in control of aberrant chromosomal replication initiation in Escherichia coli

Replication of the circular bacterial chromosome is initiated from a locus oriC with the aid of an essential protein DnaA. One approach to identify factors acting to prevent aberrant oriC-independent replication initiation in Escherichia coli has been that to obtain mutants which survive loss of DnaA. Here, we show that a ΔrecD mutation, associated with attenuation of RecBCD’s DNA double strand end-resection activity, provokes abnormal replication and rescues ΔdnaA lethality in two situations: (i) in absence of 5′-3′ single-strand DNA exonuclease RecJ, or (ii) when multiple two-ended DNA double strand breaks (DSBs) are generated either by I-SceI endonucleolytic cleavages or by radiomimetic agents phleomycin or bleomycin. One-ended DSBs in the ΔrecD mutant did not rescue ΔdnaA lethality. With two-ended DSBs in the ΔrecD strain, ΔdnaA viability was retained even after linearization of the chromosome. Data from genome-wide DNA copy number determinations in ΔdnaA-rescued cells lead us to propose a model that nuclease-mediated DNA resection activity of RecBCD is critical for prevention of a σ-mode of rolling-circle over-replication when convergent replication forks merge and fuse, as may be expected to occur during normal replication at the chromosomal terminus region or during repair of two-ended DSBs following ‘ends-in’ replication.

Genomic structural plasticity of rodent-associated Bartonella in nature

Rodent‐associated Bartonella species have shown a remarkable genetic diversity and pathogenic potential. To further explore the extent of the natural intraspecific genomic variation and its potential role as an evolutionary driver, we focused on a single genetically diverse Bartonella species, Bartonella krasnovii, which circulates among gerbils and their associated fleas. Twenty genomes from 16 different B. krasnovii genotypes were fully characterized through a genome sequencing assay (using short and long read sequencing), pulse field gel electrophoresis (PFGE), and PCR validation. Genomic analyses were performed in comparison to the B. krasnovii strain OE 1–1. While, single nucleotide polymorphism represented only a 0.3% of the genome variation, structural diversity was identified in these genomes, with an average of 51 ± 24 structural variation (SV) events per genome. Interestingly, a large proportion of the SVs (>40%) was associated with prophages. Further analyses revealed that most of the SVs, and prophage insertions were found at the chromosome replication termination site (ter), suggesting this site as a plastic zone of the B. krasnovii chromosome. Accordingly, six genomes were found to be unbalanced, and essential genes near the ter showed a shift between the leading and lagging strands, revealing the SV effect on these genomes. In summary, our findings demonstrate the extensive genomic diversity harbored by wild B. krasnovii strains and suggests that its diversification is initially promoted by structural changes, probably driven by phages. These events may constantly feed the system with novel genotypes that ultimately lead to inter‐ and intraspecies competition and adaptation.

Dynamic interactions between mega symbiosis ICEs and bacterial chromosomes maintain genome architecture

Acquisition of mobile genetic elements can confer novel traits to bacteria. Some integrative and conjugative elements confer upon members of Bradyrhizobium the capacity to fix nitrogen in symbiosis with legumes. These so-called symbiosis integrative conjugative elements (symICEs) can be extremely large and vary as monopartite and polypartite configurations within chromosomes of related strains. These features are predicted to impose fitness costs and have defied explanation. Here, we show that chromosome architecture is largely conserved despite diversity in genome composition, variations in locations of attachment sites recognized by integrases of symICEs, and differences in large-scale chromosomal changes that occur upon integration. Conversely, many simulated nonnative chromosome–symICE combinations are predicted to result in lethal deletions or disruptions to architecture. Findings suggest that there is compatibility between chromosomes and symICEs. We hypothesize that the size and structural flexibility of symICEs are important for generating combinations that maintain chromosome architecture across a genus of nitrogen-fixing bacteria with diverse and dynamic genomes.

Extended sister-chromosome catenation leads to massive reorganization of the E. coli genome

In bacteria, chromosome segregation occurs progressively from the origin to terminus within minutes of replication of each locus. Between replication and segregation, sister loci are held in an apparent cohesive state by topological links. The decatenation activity of topoisomerase IV (Topo IV) is required for segregation of replicated loci, yet little is known about the structuring of the chromosome maintained in a cohesive state. In this work, we investigated chromosome folding in cells with altered decatenation activities. Within minutes after Topo IV inactivation, massive chromosome reorganization occurs, associated with increased in contacts between nearby loci, likely trans-contacts between sister chromatids, and in long-range contacts between the terminus and distant loci. We deciphered the respective roles of Topo III, MatP and MukB when TopoIV activity becomes limiting. Topo III reduces short-range inter-sister contacts suggesting its activity near replication forks. MatP, the terminus macrodomain organizing system, and MukB, the Escherichia coli SMC, promote long-range contacts with the terminus. We propose that the large-scale conformational changes observed under these conditions reveal defective decatenation attempts involving the terminus area. Our results support a model of spatial and temporal partitioning of the tasks required for sister chromosome segregation.

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.

Grand scale genome manipulation via chromosome swapping in Escherichia coli programmed by three one megabase chromosomes

In bacterial synthetic biology, whole genome transplantation has been achieved only in mycoplasmas that contain a small genome and are competent for foreign genome uptake. In this study, we developed Escherichia coli strains programmed by three 1-megabase (Mb) chromosomes by splitting the 3-Mb chromosome of a genome-reduced strain. The first split-chromosome retains the original replication origin (oriC) and partitioning (par) system. The second one has an oriC and the par locus from the F plasmid, while the third one has the ori and par locus of the Vibrio tubiashii secondary chromosome. The tripartite-genome cells maintained the rod-shaped form and grew only twice as slowly as their parent, allowing their further genetic engineering. A proportion of these 1-Mb chromosomes were purified as covalently closed supercoiled molecules with a conventional alkaline lysis method and anion exchange columns. Furthermore, the second and third chromosomes could be individually electroporated into competent cells. In contrast, the first split-chromosome was not able to coexist with another chromosome carrying the same origin region. However, it was exchangeable via conjugation between tripartite-genome strains by using different selection markers. We believe that this E. coli-based technology has the potential to greatly accelerate synthetic biology and synthetic genomics.

Comparative in silico genome analysis of Clostridium perfringens unravels stable phylogroups with different genome characteristics and pathogenic potential

Clostridium perfringens causes a plethora of devastating infections, with toxin production being the underlying mechanism of pathogenicity in various hosts. Genomic analyses of 206 public-available C. perfringens strains´ sequence data identified a substantial degree of genomic variability in respect to episome content, chromosome size and mobile elements. However, the position and order of the local collinear blocks on the chromosome showed a considerable degree of preservation. The strains were divided into five stable phylogroups (I–V). Phylogroup I contained human food poisoning strains with chromosomal enterotoxin (cpe) and a Darmbrand strain characterized by a high frequency of mobile elements, a relatively small genome size and a marked loss of chromosomal genes, including loss of genes encoding virulence traits. These features might correspond to the adaptation of these strains to a particular habitat, causing human foodborne illnesses. This contrasts strains that belong to phylogroup II where the genome size points to the acquisition of genetic material. Most strains of phylogroup II have been isolated from enteric lesions in horses and dogs. Phylogroups III, IV and V are heterogeneous groups containing a variety of different strains, with phylogroup III being the most abundant (65.5%). In conclusion, C. perfringens displays five stable phylogroups reflecting different disease involvements, prompting further studies on the evolution of this highly important pathogen.

Characterization of ori and parS-like functions in secondary genome replicons in Deinococcus radiodurans

The mechanisms underlying multipartite genome maintenance and its functional significance in extraordinary radioresistance of Deinococcus radiodurans are not well understood. The sequences upstream to parAB operons in chrII (cisII) and MP (cisMP) could stabilize an otherwise, non-replicative colE1 plasmid, in D. radiodurans DnaA and cognate ParB proteins bound specifically with cisII and cisMP elements. The ΔcisII and ΔcisMP cells showed the reduced copy number of cognate replicons and radioresistance as compared with wild type. Fluorescent reporter-operator system inserted in chrI, chrII, and MP in wild type and cisII mutants showed the presence of all three replicons in wild-type cells. Although chrI was present in all the ΔcisII and ΔcisMP cells, nearly half of these cells had chrII and MP, respectively, and the other half had the reduced number of foci representing these replications. These results suggested that cisII and cisMP elements contain both origin of replication and parS-like functions and the secondary genome replicons (chrII and MP) are maintained independent of chrI and have roles in radioresistance of D. radiodurans.

Overcoming the Challenges of Megabase-Sized Plasmid Construction in Escherichia coli

Although Escherichia coli has been a popular tool for plasmid construction, this bacterium was believed to be "unsuitable" for constructing a large plasmid whose size exceeds 500 kilobases. We assumed that traditional plasmid vectors may lack some regulatory DNA elements required for the stable replication and segregation of such a large plasmid. In addition, the use of a few site-specific recombination systems may facilitate cloning of large DNA segments. Here we show two strategies for constructing 1-megabase (1-Mb) secondary chromosomes by using new bacterial artificial chromosome (BAC) vectors. First, the 3-Mb genome of a genome-reduced E. coli strain was split into two chromosomes (2-Mb and 1-Mb), of which the smaller one has the origin of replication and the partitioning locus of the Vibrio tubiashii secondary chromosome. This chromosome fission method (Flp-POP cloning) works via flippase-mediated excision, which coincides with the reassembly of a split chloramphenicol resistance gene, allowing chloramphenicol selection. Next, we developed a new cloning method (oriT-POP cloning) and a fully equipped BAC vector (pMegaBAC1H) for developing a 1-Mb plasmid. Two 0.5-Mb genomic regions were sequentially transferred from two donor strains to a recipient strain via conjugation and captured by pMegaBAC1H in the recipient strain to produce a 1-Mb plasmid. This 1-Mb plasmid was transmissible to another E. coli strain via conjugation. Furthermore, these 1-Mb secondary chromosomes were amplifiable in vitro by using the reconstituted E. coli chromosome replication cycle reaction (RCR). These strategies and technologies would make popular E. coli cells a productive factory for designer chromosome engineering.

Macromolecular crowding links ribosomal protein gene dosage to growth rate in Vibrio cholerae

BACKGROUND

In fast-growing bacteria, the genomic location of ribosomal protein (RP) genes is biased towards the replication origin (oriC). This trait allows optimizing their expression during exponential phase since oriC neighboring regions are in higher dose due to multifork replication. Relocation of s10-spc-α locus (S10), which codes for most of the RP, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction on its dosage, its expression, and bacterial growth rate. However, a mechanism linking S10 dosage to cell physiology has still not been determined.

RESULTS

We hypothesized that S10 dosage perturbations impact protein synthesis capacity. Strikingly, we observed that in Vibrio cholerae, protein production capacity was independent of S10 position. Deep sequencing revealed that S10 relocation altered chromosomal replication dynamics and genome-wide transcription. Such changes increased as a function of oriC-S10 distance. Since RP constitutes a large proportion of cell mass, lower S10 dosage could lead to changes in macromolecular crowding, impacting cell physiology. Accordingly, cytoplasm fluidity was higher in mutants where S10 is most distant from oriC. In hyperosmotic conditions, when crowding differences are minimized, the growth rate and replication dynamics were highly alleviated in these strains.

CONCLUSIONS

The genomic location of RP genes ensures its optimal dosage. However, besides of its essential function in translation, their genomic position sustains an optimal macromolecular crowding essential for maximizing growth. Hence, this could be another mechanism coordinating DNA replication to bacterial growth.

Exception to the exception rule: synthetic and naturally occurring single chromosome Vibrio cholerae

The genome of Vibrio cholerae, the etiological agent of cholera, is an exception to the single chromosome rule found in the vast majority of bacteria and has its genome partitioned between two unequally sized chromosomes. This unusual two-chromosome arrangement in V. cholerae has sparked considerable research interest since its discovery. It was demonstrated that the two chromosomes could be fused by deliberate genome engineering or forced to fuse spontaneously by blocking the replication of Chr2, the secondary chromosome. Recently, natural isolates of V. cholerae with chromosomal fusion have been found. Here, we summarize the pertinent findings on this exception to the exception rule and discuss the potential utility of single-chromosome V. cholerae to address fundamental questions on chromosome biology in general and DNA replication in particular.

Early fate of exogenous promoters in E. coli

Gene gain by horizontal gene transfer is a major pathway of genome innovation in bacteria. The current view posits that acquired genes initially need to be silenced and that a bacterial chromatin protein, H-NS, plays a role in this silencing. However, we lack direct observation of the early fate of a horizontally transferred gene to prove this theory. We combine sequencing, flow cytometry and sorting, followed by microscopy to monitor gene expression and its variability after large-scale random insertions of a reporter gene in a population of Escherichia coli bacteria. We find that inserted promoters have a wide range of gene-expression variability related to their location. We find that high-expression clones carry insertions that are not correlated with H-NS binding. Conversely, binding of H-NS correlates with silencing. Finally, while most promoters show a common level of extrinsic noise, some insertions show higher noise levels. Analysis of these high-noise clones supports a scenario of switching due to transcriptional interference from divergent ribosomal promoters. Altogether, our findings point to evolutionary pathways where newly-acquired genes are not necessarily silenced, but may immediately explore a wide range of expression levels to probe the optimal ones.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Laboratory Evolution Experiments Help Identify a Predominant Region of Constitutive Stable DNA Replication Initiation

The bacterium Escherichia coli can initiate replication in the absence of the replication initiator protein DnaA and/or the canonical origin of replication oriC in a ΔrnhA background. This phenomenon, which can be primed by R-loops, is called constitutive stable DNA replication (cSDR). Whether DNA replication during cSDR initiates in a stochastic manner through the length of the chromosome or at specific sites and how E. coli can find adaptations to loss of fitness caused by cSDR remain inadequately answered. We use laboratory evolution experiments of ΔrnhA-ΔdnaA strains followed by deep sequencing to show that DNA replication preferentially initiates within a broad region located ∼0.4 to 0.7 Mb clockwise of oriC. This region includes many bisulfite-sensitive sites, which have been previously defined as R-loop-forming regions, and includes a site containing sequence motifs that favor R-loop formation. Initiation from this region would result in head-on replication-transcription conflicts at rRNA loci. Inversions of these rRNA loci, which can partly resolve these conflicts, help the bacterium suppress the fitness defects of cSDR. These inversions partially restore the gene expression changes brought about by cSDR. The inversion, however, increases the possibility of conflicts at essential mRNA genes, which would utilize only a minuscule fraction of RNA polymerase molecules, most of which transcribe rRNA genes. Whether subsequent adaptive strategies would attempt to resolve these conflicts remains an open question.IMPORTANCE The bacterium E. coli can replicate its DNA even in the absence of the molecules that are required for canonical replication initiation. This often requires the formation of RNA-DNA hybrid structures and is referred to as constitutive stable DNA replication (cSDR). Where on the chromosome does cSDR initiate? We answer this question using laboratory evolution experiments and genomics and show that selection favors cSDR initiation predominantly at a region ∼0.6 Mb clockwise of oriC. Initiation from this site will result in more head-on collisions of DNA polymerase with RNA polymerase operating on rRNA loci. The bacterium adapts to this problem by inverting a region of the genome including several rRNA loci such that head-on collisions between the two polymerases are minimized. Understanding such evolutionary strategies in the context of cSDR can provide insights into the potential causes of resistance to antibiotics that target initiation of DNA replication.

Cell morphology and nucleoid dynamics in dividing Deinococcus radiodurans

Our knowledge of bacterial nucleoids originates mostly from studies of rod- or crescent-shaped bacteria. Here we reveal that Deinococcus radiodurans, a relatively large spherical bacterium with a multipartite genome, constitutes a valuable system for the study of the nucleoid in cocci. Using advanced microscopy, we show that D. radiodurans undergoes coordinated morphological changes at both the cellular and nucleoid level as it progresses through its cell cycle. The nucleoid is highly condensed, but also surprisingly dynamic, adopting multiple configurations and presenting an unusual arrangement in which oriC loci are radially distributed around clustered ter sites maintained at the cell centre. Single-particle tracking and fluorescence recovery after photobleaching studies of the histone-like HU protein suggest that its loose binding to DNA may contribute to this remarkable plasticity. These findings demonstrate that nucleoid organization is complex and tightly coupled to cell cycle progression in this organism.

DNA lesions proximity modulates damage tolerance pathways in Escherichia coli

The genome of all organisms is constantly threatened by numerous agents that cause DNA damage. When the replication fork encounters an unrepaired DNA lesion, two DNA damage tolerance pathways are possible: error-prone translesion synthesis (TLS) that requires specialized DNA polymerases, and error-free damage avoidance that relies on homologous recombination (HR). The balance between these two mechanisms is essential since it defines the level of mutagenesis during lesion bypass, allowing genetic variability and adaptation to the environment, but also introduces the risk of generating genome instability. Here we report that the mere proximity of replication-blocking lesions that arise in Escherichia coli's genome during a genotoxic stress leads to a strong increase in the use of the error-prone TLS. We show that this increase is caused by the local inhibition of HR due to the overlapping of single-stranded DNA regions generated downstream of the lesions. This increase in TLS is independent of SOS activation, but its mutagenic effect is additive with the one of SOS. Hence, the combination of SOS induction and lesions proximity leads to a strong increase in TLS that becomes the main lesion tolerance pathway used by the cell during a genotoxic stress.

Replication termination without a replication fork trap

Bacterial chromosomes harbour a unique origin of bidirectional replication, oriC. They are almost always circular, with replication terminating in a region diametrically opposite to oriC, the terminus. The oriC-terminus organisation is reflected by the orientation of the genes and by the disposition of DNA-binding protein motifs implicated in the coordination of chromosome replication and segregation with cell division. Correspondingly, the E. coli and B. subtilis model bacteria possess a replication fork trap system, Tus/ter and RTP/ter, respectively, which enforces replication termination in the terminus region. Here, we show that tus and rtp are restricted to four clades of bacteria, suggesting that tus was recently domesticated from a plasmid gene. We further demonstrate that there is no replication fork system in Vibrio cholerae, a bacterium closely related to E. coli. Marker frequency analysis showed that replication forks originating from ectopic origins were not blocked in the terminus region of either of the two V. cholerae chromosomes, but progressed normally until they encountered an opposite fork. As expected, termination synchrony of the two chromosomes is disrupted by these ectopic origins. Finally, we show that premature completion of the primary chromosome replication did not modify the choreography of segregation of its terminus region.

Transcriptome and Comparative Genomics Analyses Reveal New Functional Insights on Key Determinants of Pathogenesis and Interbacterial Competition in Pectobacterium and Dickeya spp

Soft-rot Enterobacteriaceae (SRE), typified by Pectobacterium and Dickeya genera, are phytopathogenic bacteria inflicting soft-rot disease in crops worldwide. By combining genomic information from 100 SRE with whole-transcriptome data sets, we identified novel genomic and transcriptional associations among key pathogenicity themes in this group. Comparative genomics revealed solid linkage between the type I secretion system (T1SS) and the carotovoricin bacteriophage (Ctv) conserved in 96.7% of Pectobacterium genomes. Moreover, their coactivation during infection indicates a novel functional association involving T1SS and Ctv. Another bacteriophage-borne genomic region, mostly confined to less than 10% of Pectobacterium strains, was found, presumably comprising a novel lineage-specific prophage in the genus. We also detected the transcriptional coregulation of a previously predicted toxin/immunity pair (WHH and SMI1_KNR4 families), along with the type VI secretion system (T6SS), which includes hcp and/or vgrG genes, suggesting a role in disease development as T6SS-dependent effectors. Further, we showed that another predicted T6SS-dependent endonuclease (AHH family) exhibited toxicity in ectopic expression assays, indicating antibacterial activity. Additionally, we report the striking conservation of the group 4 capsule (GFC) cluster in 100 SRE strains which consistently features adjacently conserved serotype-specific gene arrays comprising a previously unknown organization in GFC clusters. Also, extensive sequence variations found in gfcA orthologs suggest a serotype-specific role in the GfcABCD machinery.IMPORTANCE Despite the considerable loss inflicted on important crops yearly by Pectobacterium and Dickeya diseases, investigations on key virulence and interbacterial competition assets relying on extensive comparative genomics are still surprisingly lacking for these genera. Such approaches become more powerful over time, underpinned by the growing amount of genomic information in public databases. In particular, our findings point to new functional associations among well-known genomic themes enabling alternative means of neutralizing SRE diseases through disruption of pivotal virulence programs. By elucidating novel transcriptional and genomic associations, this study adds valuable information on virulence candidates that could be decisive in molecular applications in the near future. The utilization of 100 genomes of Pectobacterium and Dickeya strains in this study is unprecedented for comparative analyses in these taxa, and it provides novel insights on the biology of economically important plant pathogens.

Evolution of Genome Architecture in Archaea: Spontaneous Generation of a New Chromosome in Haloferax volcanii

The common ancestry of archaea and eukaryotes is evident in their genome architecture. All eukaryotic and several archaeal genomes consist of multiple chromosomes, each replicated from multiple origins. Three scenarios have been proposed for the evolution of this genome architecture: 1) mutational diversification of a multi-copy chromosome; 2) capture of a new chromosome by horizontal transfer; 3) acquisition of new origins and splitting into two replication-competent chromosomes. We report an example of the third scenario: the multi-origin chromosome of the archaeon Haloferax volcanii has split into two elements via homologous recombination. The newly generated elements are bona fide chromosomes, because each bears "chromosomal" replication origins, rRNA loci, and essential genes. The new chromosomes were stable during routine growth but additional genetic manipulation, which involves selective bottlenecks, provoked further rearrangements. To the best of our knowledge, rearrangement of a naturally evolved prokaryotic genome to generate two new chromosomes has not been described previously.

Inserting Site-Specific DNA Lesions into Whole Genomes

Here, we describe a methodology that allows the insertion of site-specific DNA lesions into genomes in living cells. The technique involves the integration of a plasmid containing a site-specific lesion engineered in vitro into a precise location in the genome via the site-specific recombination reaction from phage lambda. The notion of DNA lesion is not restricted to chemically modified nucleotides but also refers to unusual DNA structures. This method will be instrumental to study qualitatively and quantitatively the genetic consequences of site-specific lesions in vivo; moreover, it does also allow analyzing the molecular structure of stalled replication forks at well-defined locations.

Non-Random Inversion Landscapes in Prokaryotic Genomes Are Shaped by Heterogeneous Selection Pressures

Inversions are a major contributor to structural genome evolution in prokaryotes. Here, using a novel alignment-based method, we systematically compare 1,651 bacterial and 98 archaeal genomes to show that inversion landscapes are frequently biased toward (symmetric) inversions around the origin–terminus axis. However, symmetric inversion bias is not a universal feature of prokaryotic genome evolution but varies considerably across clades. At the extremes, inversion landscapes in Bacillus–Clostridium and Actinobacteria are dominated by symmetric inversions, while there is little or no systematic bias favoring symmetric rearrangements in archaea with a single origin of replication. Within clades, we find strong but clade-specific relationships between symmetric inversion bias and different features of adaptive genome architecture, including the distance of essential genes to the origin of replication and the preferential localization of genes on the leading strand. We suggest that heterogeneous selection pressures have converged to produce similar patterns of structural genome evolution across prokaryotes.

Broad-scale redistribution of mRNA abundance and transcriptional machinery in response to growth rate in Salmonella enterica serovar Typhimurium

We have investigated the connection between the four-dimensional architecture of the bacterial nucleoid and the organism's global gene expression programme. By localizing the transcription machinery and the transcriptional outputs across the genome of the model bacterium Salmonella enterica serovar Typhimurium at different stages of the growth cycle, a surprising disconnection between gene dosage and transcriptional output was revealed. During exponential growth, gene output occurred chiefly in the Ori (origin), Ter (terminus) and NSL (non-structured left) domains, whereas the Left macrodomain remained transcriptionally quiescent at all stages of growth. The apparently high transcriptional output in Ter was correlated with an enhanced stability of the RNA expressed there during exponential growth, suggesting that longer mRNA half-lives compensate for low gene dosage. During exponential growth, RNA polymerase (RNAP) was detected everywhere, whereas in stationary phase cells, RNAP was concentrated in the Ter macrodomain. The alternative sigma factors RpoE, RpoH and RpoN were not required to drive transcription in these growth conditions, consistent with their observed binding to regions away from RNAP and regions of active transcription. Specifically, these alternative sigma factors were found in the Ter macrodomain during exponential growth, whereas they were localized at the Ori macrodomain in stationary phase.

The Divided Bacterial Genome: Structure, Function, and Evolution

Approximately 10% of bacterial genomes are split between two or more large DNA fragments, a genome architecture referred to as a multipartite genome. This multipartite organization is found in many important organisms, including plant symbionts, such as the nitrogen-fixing rhizobia, and plant, animal, and human pathogens, including the genera Brucella, Vibrio, and Burkholderia. The availability of many complete bacterial genome sequences means that we can now examine on a broad scale the characteristics of the different types of DNA molecules in a genome. Recent work has begun to shed light on the unique properties of each class of replicon, the unique functional role of chromosomal and nonchromosomal DNA molecules, and how the exploitation of novel niches may have driven the evolution of the multipartite genome. The aims of this review are to (i) outline the literature regarding bacterial genomes that are divided into multiple fragments, (ii) provide a meta-analysis of completed bacterial genomes from 1,708 species as a way of reviewing the abundant information present in these genome sequences, and (iii) provide an encompassing model to explain the evolution and function of the multipartite genome structure. This review covers, among other topics, salient genome terminology; mechanisms of multipartite genome formation; the phylogenetic distribution of multipartite genomes; how each part of a genome differs with respect to genomic signatures, genetic variability, and gene functional annotation; how each DNA molecule may interact; as well as the costs and benefits of this genome structure.

A non-catalytic role of RecBCD in homology directed gap repair and translesion synthesis

The RecBCD complex is a key factor in DNA metabolism. This protein complex harbors a processive nuclease and two helicases activities that give it the ability to process duplex DNA ends. These enzymatic activities make RecBCD a major player in double strand break repair, conjugational recombination and degradation of linear DNA. In this work, we unravel a new role of the RecBCD complex in the processing of DNA single-strand gaps that are generated at DNA replication-blocking lesions. We show that independently of its nuclease or helicase activities, the entire RecBCD complex is required for recombinational repair of the gap and efficient translesion synthesis. Since none of the catalytic functions of RecBCD are required for those processes, we surmise that the complex acts as a structural element that stabilizes the blocked replication fork, allowing efficient DNA damage tolerance.

Long range chromosome organization in Escherichia coli: The position of the replication origin defines the non-structured regions and the Right and Left macrodomains

The Escherichia coli chromosome is organized into four macrodomains (Ori, Ter, Right and Left) and two non-structured regions. This organization influences the segregation of sister chromatids, the mobility of chromosomal DNA, and the cellular localization of the chromosome. The organization of the Ter and Ori macrodomains relies on two specific systems, MatP/matS for the Ter domain and MaoP/maoS for the Ori domain, respectively. Here by constructing strains with chromosome rearrangements to reshuffle the distribution of chromosomal segments, we reveal that the difference between the non-structured regions and the Right and Left lateral macrodomains relies on their position on the chromosome. A change in the genetic location of oriC generated either by an inversion within the Ori macrodomain or by the insertion of a second oriC modifies the position of Right and Left macrodomains, as the chromosome region the closest to oriC are always non-structured while the regions further away behave as macrodomain regardless of their DNA sequence. Using fluorescent microscopy we estimated that loci belonging to a non-structured region are significantly closer to the Ori MD than loci belonging to a lateral MD. Altogether, our results suggest that the origin of replication plays a prominent role in chromosome organization in E. coli, as it determines structuring and localization of macrodomains in growing cell.

Chromosomes allow the genetic information to be stored, transmitted and organized inside the cell. In bacteria, chromosomes are generally circular and they are shaped and organized by several mechanisms allowing simultaneous transcription, replication and segregation. The way such fundamental processes are managed is still unclear, but in the Gram negative bacteria Escherichia coli, one level of chromosome organization relies on a large scale structuring of the chromosome in macrodomains and non-structured regions. Macrodomains have been defined as chromosomal domains of megabase (Mb) sized genetically isolated from each other. E. coli chromosome is divided in 4 macrodomains (Ter, Ori, Right and Left) and two right/left non-structured regions (NSR, and NSL). Factors that organize Ter and Ori MD have been identified and characterized previously: MatP structures the Ter MD by binding to 23 matS sequences disseminated in the Ter MD, and MaoP together with the unique sequence maoS organises the Ori MD by an unknown mechanism. In constrast, we show here by reshuffling the chromosome that the Right and Left MD as well as the NSR and NSL regions are defined by their chromosomal location and that the chromosomal position of oriC defines the position and the extent of the NS regions and Right/Left macrodomains.

The History of Bordetella pertussis Genome Evolution Includes Structural Rearrangement

Despite high pertussis vaccine coverage, reported cases of whooping cough (pertussis) have increased over the last decade in the United States and other developed countries. Although Bordetella pertussis is well known for its limited gene sequence variation, recent advances in long-read sequencing technology have begun to reveal genomic structural heterogeneity among otherwise indistinguishable isolates, even within geographically or temporally defined epidemics. We have compared rearrangements among complete genome assemblies from 257 B. pertussis isolates to examine the potential evolution of the chromosomal structure in a pathogen with minimal gene nucleotide sequence diversity. Discrete changes in gene order were identified that differentiated genomes from vaccine reference strains and clinical isolates of various genotypes, frequently along phylogenetic boundaries defined by single nucleotide polymorphisms. The observed rearrangements were primarily large inversions centered on the replication origin or terminus and flanked by IS481, a mobile genetic element with >240 copies per genome and previously suspected to mediate rearrangements and deletions by homologous recombination. These data illustrate that structural genome evolution in B. pertussis is not limited to reduction but also includes rearrangement. Therefore, although genomes of clinical isolates are structurally diverse, specific changes in gene order are conserved, perhaps due to positive selection, providing novel information for investigating disease resurgence and molecular epidemiology.IMPORTANCE Whooping cough, primarily caused by Bordetella pertussis, has resurged in the United States even though the coverage with pertussis-containing vaccines remains high. The rise in reported cases has included increased disease rates among all vaccinated age groups, provoking questions about the pathogen's evolution. The chromosome of B. pertussis includes a large number of repetitive mobile genetic elements that obstruct genome analysis. However, these mobile elements facilitate large rearrangements that alter the order and orientation of essential protein-encoding genes, which otherwise exhibit little nucleotide sequence diversity. By comparing the complete genome assemblies from 257 isolates, we show that specific rearrangements have been conserved throughout recent evolutionary history, perhaps by eliciting changes in gene expression, which may also provide useful information for molecular epidemiology.

Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra

Oxidative stress is a major cause of mutation but little is known about how growth in the absence of oxygen impacts the rate and spectrum of mutations. We employed long-term mutation accumulation experiments to directly measure the rates and spectra of spontaneous mutation events in Escherichia coli populations propagated under aerobic and anaerobic conditions. To detect mutations, whole genome sequencing was coupled with methods of analysis sufficient to identify a broad range of mutational classes, including structural variants (SVs) generated by movement of repetitive elements. The anaerobically grown populations displayed a mutation rate nearly twice that of the aerobic populations, showed distinct asymmetric mutational strand biases, and greater insertion element activity. Consistent with mutation rate and spectra observations, genes for transposition and recombination repair associated with SVs were up-regulated during anaerobic growth. Together, these results define differences in mutational spectra affecting the evolution of facultative anaerobes.

Novel Technologies for Optimal Strain Breeding

The implementation of a knowledge-based bioeconomy requires the rapid development of highly efficient microbial production strains that are able to convert renewable carbon sources to value-added products, such as bulk and fine chemicals, pharmaceuticals, or proteins at industrial scale. Starting from classical strain breeding by random mutagenesis and screening in the 1950s via rational design by metabolic engineering initiated in the 1970s, a range of powerful new technologies have been developed in the past two decades that can revolutionize future strain engineering. In particular, next-generation sequencing technologies combined with new methods of genome engineering and high-throughput screening based on genetically encoded biosensors have allowed for new concepts. In this chapter, selected new technologies relevant for breeding microbial production strains with a special emphasis on amino acid producers will be summarized.

The MaoP/maoS Site-Specific System Organizes the Ori Region of the E. coli Chromosome into a Macrodomain

The Ori region of bacterial genomes is segregated early in the replication cycle of bacterial chromosomes. Consequently, Ori region positioning plays a pivotal role in chromosome dynamics. The Ori region of the E. coli chromosome is organized as a macrodomain with specific properties concerning DNA mobility, segregation of loci and long distance DNA interactions. Here, by using strains with chromosome rearrangements and DNA mobility as a read-out, we have identified the MaoP/maoS system responsible for constraining DNA mobility in the Ori region and limiting long distance DNA interactions with other regions of the chromosome. MaoP belongs to a group of proteins conserved in the Enterobacteria that coevolved with Dam methylase including SeqA, MukBEF and MatP that are all involved in the control of chromosome conformation and segregation. Analysis of DNA rings excised from the chromosome demonstrated that the single maoS site is required in cis on the chromosome to exert its effect while MaoP can act both in cis and in trans. The position of markers in the Ori region was affected by inactivating maoP. However, the MaoP/maoS system was not sufficient for positioning the Ori region at the ¼–¾ regions of the cell. We also demonstrate that the replication and the resulting expansion of bulk DNA are localized centrally in the cell. Implications of these results for chromosome positioning and segregation in E. coli are discussed.

The Ori region from bacterial chromosomes plays a pivotal role in chromosome organization and segregation as it is replicated and segregated early in cell division cycle and its positioning impacts the cellular organization of the chromosome in the cell. The E. coli chromosome is divided into four macrodomains (MD) defined as large regions in which DNA interactions occurred preferentially. Here we have identified a new system responsible for specifying properties to the Ori MD. This system is composed of two elements: a cis-acting target sequence called maoS and a gene of unknown function acting in trans called maoP. Remarkably, MaoP belongs to a group of proteins conserved only in Enterobacteria that coevolved with the Dam DNA methylase and that includes the MatP protein structuring the Ter macrodomain and the SeqA and MukBEF proteins involved in the control of chromosome conformation and segregation. These results reveal the presence of a dedicated set of factors required in chromosome management in enterobacteria that might compensate, at least partially, for the absence of the ParABS system involved in the condensation and/or segregation of the Ori region in most bacteria.

DNA Replication Control Is Linked to Genomic Positioning of Control Regions in Escherichia coli

Chromosome replication in Escherichia coli is in part controlled by three non-coding genomic sequences, DARS1, DARS2, and datA that modulate the activity of the initiator protein DnaA. The relative distance from oriC to the non-coding regions are conserved among E. coli species, despite large variations in genome size. Here we use a combination of i) site directed translocation of each region to new positions on the bacterial chromosome and ii) random transposon mediated translocation followed by culture evolution, to show genetic evidence for the importance of position. Here we provide evidence that the genomic locations of these regulatory sequences are important for cell cycle control and bacterial fitness. In addition, our work shows that the functionally redundant DARS1 and DARS2 regions play different roles in replication control. DARS1 is mainly involved in maintaining the origin concentration, whether DARS2 is also involved in maintaining single cell synchrony.

Comparative Genomics of Interreplichore Translocations in Bacteria: A Measure of Chromosome Topology?

Genomes evolve not only in base sequence but also in terms of their architecture, defined by gene organization and chromosome topology. Whereas genome sequence data inform us about the changes in base sequences for a large variety of organisms, the study of chromosome topology is restricted to a few model organisms studied using microscopy and chromosome conformation capture techniques. Here, we exploit whole genome sequence data to study the link between gene organization and chromosome topology in bacteria. Using comparative genomics across ∼250 pairs of closely related bacteria we show that: (a) many organisms show a high degree of interreplichore translocations throughout the chromosome and not limited to the inversion-prone terminus (ter) or the origin of replication (oriC); (b) translocation maps may reflect chromosome topologies; and (c) symmetric interreplichore translocations do not disrupt the distance of a gene from oriC or affect gene expression states or strand biases in gene densities. In summary, we suggest that translocation maps might be a first line in defining a gross chromosome topology given a pair of closely related genome sequences.

A defect in homologous recombination leads to increased translesion synthesis in E. coli

DNA damage tolerance pathways allow cells to duplicate their genomes despite the presence of replication blocking lesions. Cells possess two major tolerance strategies, namely translesion synthesis (TLS) and homology directed gap repair (HDGR). TLS pathways involve specialized DNA polymerases that are able to synthesize past DNA lesions with an intrinsic risk of causing point mutations. In contrast, HDGR pathways are essentially error-free as they rely on the recovery of missing information from the sister chromatid by RecA-mediated homologous recombination. We have investigated the genetic control of pathway choice between TLS and HDGR in vivo in Escherichia coli In a strain with wild type RecA activity, the extent of TLS across replication blocking lesions is generally low while HDGR is used extensively. Interestingly, recA alleles that are partially impaired in D-loop formation confer a decrease in HDGR and a concomitant increase in TLS. Thus, partial defect of RecA's capacity to invade the homologous sister chromatid increases the lifetime of the ssDNA.RecA filament, i.e. the 'SOS signal'. This increase favors TLS by increasing both the TLS polymerase concentration and the lifetime of the TLS substrate, before it becomes sequestered by homologous recombination. In conclusion, the pathway choice between error-prone TLS and error-free HDGR is controlled by the efficiency of homologous recombination.

Fork restart protein, PriA, binds around oriC after depletion of nucleotide precursors: Replication fork arrest near the replication origin

Arrest of replication fork progression is one of the most common causes for increasing the genomic instability. In bacteria, PriA, a conserved DEXH-type helicase, plays a major role in recognition of the stalled forks and restart of DNA replication. We took advantage of PriA's ability to specifically recognize stalled replication forks to determine the genomic loci where replication forks are prone to stall on the Escherichia coli genome. We found that PriA binds around oriC upon thymine starvation which reduces the nucleotide supply and causes replication fork stalling. PriA binding quickly disappeared upon readdition of thymine. Furthermore, BrdU was incorporated at around oriC upon release from thymine starvation. Our results indicate that reduced supply of DNA replication precursors causes replication fork stalling preferentially in the 600 kb segment centered at oriC. This suggests that replication of the vicinity of oriC requires higher level of nucleotide precursors. The results also point to a possibility of slow fork movement and/or the presence of multiple fork arrest signals within this segment. Indeed, we have identified rather strong fork stall/pausing signals symmetrically located at ∼50 kb away from oriC. We speculate that replication pausing and fork-slow-down shortly after initiation may represent a novel checkpoint that ensures the presence of sufficient nucleotide supply prior to commitment to duplication of the entire genome.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

Bacterial Proliferation: Keep Dividing and Don't Mind the Gap

DNA Damage Tolerance (DDT) mechanisms help dealing with unrepaired DNA lesions that block replication and challenge genome integrity. Previous in vitro studies showed that the bacterial replicase is able to re-prime downstream of a DNA lesion, leaving behind a single-stranded DNA gap. The question remains of what happens to this gap in vivo. Following the insertion of a single lesion in the chromosome of a living cell, we showed that this gap is mostly filled in by Homology Directed Gap Repair in a RecA dependent manner. When cells fail to repair this gap, or when homologous recombination is impaired, cells are still able to divide, leading to the loss of the damaged chromatid, suggesting that bacteria lack a stringent cell division checkpoint mechanism. Hence, at the expense of losing one chromatid, cell survival and proliferation are ensured.

DNA Damage Tolerance (DDT) mechanisms help dealing with unrepaired DNA lesions that block replication, thus challenging genome integrity. Two DDT mechanisms have previously been described: error prone Translesion Synthesis operated by specialized DNA polymerases and error free bypass that uses the information of the sister chromatid to bypass the lesion. In this work, we set up a novel genetic system that allows to insert a single DNA blocking lesion in the chromosome of a living cell and to visualize the exchange of genetic information between the undamaged and the damaged strand. Using this system, we showed in vivo that the replication fork is able to re-prime downstream of the lesion, leaving a gap. This gap is mostly filled in by the error free pathway through the RecA homologous recombination mechanism. We show that when the gap is left unrepaired, cells are still able to divide by losing the damaged chromatid, which evidences the lack of a stringent cell division checkpoint system.

Control regions for chromosome replication are conserved with respect to sequence and location among Escherichia coli strains

In Escherichia coli, chromosome replication is initiated from oriC by the DnaA initiator protein associated with ATP. Three non-coding regions contribute to the activity of DnaA. The datA locus is instrumental in conversion of DnaA^ATP to DnaA^ADP (datA dependent DnaA^ATP hydrolysis) whereas DnaA rejuvenation sequences 1 and 2 (DARS1 and DARS2) reactivate DnaA^ADP to DnaA^ATP. The structural organization of oriC, datA, DARS1, and DARS2 were found conserved among 59 fully sequenced E. coli genomes, with differences primarily in the non-functional spacer regions between key protein binding sites. The relative distances from oriC to datA, DARS1, and DARS2, respectively, was also conserved despite of large variations in genome size, suggesting that the gene dosage of either region is important for bacterial growth. Yet all three regions could be deleted alone or in combination without loss of viability. Competition experiments during balanced growth in rich medium and during mouse colonization indicated roles of datA, DARS1, and DARS2 for bacterial fitness although the relative contribution of each region differed between growth conditions. We suggest that this fitness advantage has contributed to conservation of both sequence and chromosomal location for datA, DARS1, and DARS2.

Comparative genomics reveals conserved positioning of essential genomic clusters in highly rearranged Thermococcales chromosomes

The genomes of the 21 completely sequenced Thermococcales display a characteristic high level of rearrangements. As a result, the prediction of their origin and termination of replication on the sole basis of chromosomal DNA composition or skew is inoperative. Using a different approach based on biologically relevant sequences, we were able to determine oriC position in all 21 genomes. The position of dif, the site where chromosome dimers are resolved before DNA segregation could be predicted in 19 genomes. Computation of the core genome uncovered a number of essential gene clusters with a remarkably stable chromosomal position across species, in sharp contrast with the scrambled nature of their genomes. The active chromosomal reorganization of numerous genes acquired by horizontal transfer, mainly from mobile elements, could explain this phenomenon.

•

Thermococcales chromosomal landmarks were uncovered using biologically relevant sequences.

•

Core genomes procedures predict integration of mobile elements on Thermococcales chromosomes.

•

Thermococcales genomes are highly rearranged but core clusters positions remain invariable.

•

Thermococcales core genes are more expressed and predominantly encoded on the leading strand.

MEGA (Multiple Essential Genes Assembling) deletion and replacement method for genome reduction in Escherichia coli

Top-down reduction of the bacterial genome to construct desired chassis cells is important for synthetic biology. However, the current progress in the field of genome reduction is greatly hindered by indispensable life-essential genes that are interspersed throughout the chromosomal loci. Here, we described a new method designated as "MEGA (Multiple Essential Genes Assembling) deletion and replacement" that functions by assembling multiple essential genes in an E. coli-S. cerevisiae shuttle vector, removing targeted chromosomal regions containing essential and nonessential genes using a one-round deletion, and then integrating the cloned essential genes into the in situ chromosomal loci via I-SceI endonuclease cleavage. As a proof of concept, we separately generated three large deletions (80-205 kbp) in the E. coli MDS42 chromosome. We believe that the MEGA deletion and replacement method has potential to become widely used in large-scale genome reductions in other sequenced organisms in addition to E. coli.

The B. subtilis Accessory Helicase PcrA Facilitates DNA Replication through Transcription Units

In bacteria the concurrence of DNA replication and transcription leads to potentially deleterious encounters between the two machineries, which can occur in either the head-on (lagging strand genes) or co-directional (leading strand genes) orientations. These conflicts lead to replication fork stalling and can destabilize the genome. Both eukaryotic and prokaryotic cells possess resolution factors that reduce the severity of these encounters. Though Escherichia coli accessory helicases have been implicated in the mitigation of head-on conflicts, direct evidence of these proteins mitigating co-directional conflicts is lacking. Furthermore, the endogenous chromosomal regions where these helicases act, and the mechanism of recruitment, have not been identified. We show that the essential Bacillus subtilis accessory helicase PcrA aids replication progression through protein coding genes of both head-on and co-directional orientations, as well as rRNA and tRNA genes. ChIP-Seq experiments show that co-directional conflicts at highly transcribed rRNA, tRNA, and head-on protein coding genes are major targets of PcrA activity on the chromosome. Partial depletion of PcrA renders cells extremely sensitive to head-on conflicts, linking the essential function of PcrA to conflict resolution. Furthermore, ablating PcrA’s ATPase/helicase activity simultaneously increases its association with conflict regions, while incapacitating its ability to mitigate conflicts, and leads to cell death. In contrast, disruption of PcrA’s C-terminal RNA polymerase interaction domain does not impact its ability to mitigate conflicts between replication and transcription, its association with conflict regions, or cell survival. Altogether, this work establishes PcrA as an essential factor involved in mitigating transcription-replication conflicts and identifies chromosomal regions where it routinely acts. As both conflicts and accessory helicases are found in all domains of life, these results are broadly relevant.

In bacteria the concurrence of DNA replication and transcription leads to potentially deleterious encounters between the two machineries. These encounters can destabilize the genome and lead to mutations. Both eukaryotic and prokaryotic cells possess conflict resolution factors that reduce the detrimental effects of these collisions. In this study we show that without the essential Bacillus subtilis accessory DNA helicase, PcrA, the replication machinery slows down at certain regions of the chromosome in a transcription-dependent manner. PcrA is essential to life but incomplete depletion of PcrA only partially inhibits cell survival. We find that, under these conditions, partial survival defects are significantly exacerbated in the presence of a single severe conflict. In summary our work identifies a high degree of conservation for accessory helicase function in conflict resolution, directly establishes PcrA’s role in co-directional conflict resolution, and maps the natural chromosomal regions where such activities are routinely needed. Because both conflicts and accessory helicases are found in all domains of life, the results of this work are broadly relevant.

Francisella tularensis Subtype A.II Genomic Plasticity in Comparison with Subtype A.I

Although Francisella tularensis is considered a monomorphic intracellular pathogen, molecular genotyping and virulence studies have demonstrated important differences within the tularensis subspecies (type A). To evaluate genetic variation within type A strains, sequencing and assembly of a new subtype A.II genome was achieved for comparison to other completed F. tularensis type A genomes. In contrast with the F. tularensis A.I strains (SCHU S4, FSC198, NE061598, and TI0902), substantial genomic variation was observed between the newly sequenced F. tularensis A.II strain (WY-00W4114) and the only other publically available A.II strain (WY96-3418). Genome differences between WY-00W4114 and WY96-3418 included three major chromosomal translocations, 1580 indels, and 286 nucleotide substitutions of which 159 were observed in predicted open reading frames and 127 were located in intergenic regions. The majority of WY-00W4114 nucleotide deletions occurred in intergenic regions, whereas most of the insertions and substitutions occurred in predicted genes. Of the nucleotide substitutions, 48 (30%) were synonymous and 111 (70%) were nonsynonymous. WY-00W4114 and WY96-3418 nucleotide polymorphisms were predominantly G/C to A/T allelic mutations, with WY-00W4114 having more A+T enrichment. In addition, the A.II genomes contained a considerably higher number of intact genes and longer repetitive sequences, including transposon remnants than the A.I genomes. Together these findings support the premise that F. tularensis A.II may have a fitness advantage compared to the A.I subtype due to the higher abundance of functional genes and repeated chromosomal sequences. A better understanding of the selective forces driving F. tularensis genetic diversity and plasticity is needed.