Conserved patterns in bacterial genomes: a conundrum physically tailored by evolutionary tinkering

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Computational Tools for the Multiscale Analysis of Hi-C Data in Bacterial Chromosomes

Just as in eukaryotes, high-throughput chromosome conformation capture (Hi-C) data have revealed nested organizations of bacterial chromosomes into overlapping interaction domains. In this chapter, we present a multiscale analysis framework aiming at capturing and quantifying these properties. These include both standard tools (e.g., contact laws) and novel ones such as an index that allows identifying loci involved in domain formation independently of the structuring scale at play. Our objective is twofold. On the one hand, we aim at providing a full, understandable Python/Jupyter-based code which can be used by both computer scientists and biologists with no advanced computational background. On the other hand, we discuss statistical issues inherent to Hi-C data analysis, focusing more particularly on how to properly assess the statistical significance of results. As a pedagogical example, we analyze data produced in Pseudomonas aeruginosa, a model pathogenetic bacterium. All files (codes and input data) can be found on a GitHub repository. We have also embedded the files into a Binder package so that the full analysis can be run on any machine through Internet.

Supercoiled DNA and non-equilibrium formation of protein complexes: A quantitative model of the nucleoprotein ParBS partition complex

ParABS, the most widespread bacterial DNA segregation system, is composed of a centromeric sequence, parS, and two proteins, the ParA ATPase and the ParB DNA binding proteins. Hundreds of ParB proteins assemble dynamically to form nucleoprotein parS-anchored complexes that serve as substrates for ParA molecules to catalyze positioning and segregation events. The exact nature of this ParBS complex has remained elusive, what we address here by revisiting the Stochastic Binding model (SBM) introduced to explain the non-specific binding profile of ParB in the vicinity of parS. In the SBM, DNA loops stochastically bring loci inside a sharp cluster of ParB. However, previous SBM versions did not include the negative supercoiling of bacterial DNA, leading to use unphysically small DNA persistences to explain the ParB binding profiles. In addition, recent super-resolution microscopy experiments have revealed a ParB cluster that is significantly smaller than previous estimations and suggest that it results from a liquid-liquid like phase separation. Here, by simulating the folding of long (≥ 30 kb) supercoiled DNA molecules calibrated with realistic DNA parameters and by considering different possibilities for the physics of the ParB cluster assembly, we show that the SBM can quantitatively explain the ChIP-seq ParB binding profiles without any fitting parameter, aside from the supercoiling density of DNA, which, remarkably, is in accord with independent measurements. We also predict that ParB assembly results from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., ParB clusters behave like liquid-like protein condensates with unconventional “leaky” boundaries.

In bacteria, faithful genome inheritance requires the two replicated DNA molecules to be segregated at the opposite halves of the cell. ParABS, the most widespread bacterial DNA segregation system, is composed of a centromere sequence, parS, and two proteins, the ParA ATPase and the ParB DNA binding protein. Hundreds of ParB assemble dynamically to form clusters around parS, which then serve as substrates for ParA molecules to catalyze the positioning and segregation events. The nature of these clusters and their interaction with DNA have remained elusive. Here, we propose a realistic minimal model that captures quantitatively the peculiar DNA binding profile of ParB in the vicinity of parS in Escherichia coli. From the viewpoint of DNA, the only fitting parameter is the in vivo supercoiling density resulting from the removal of DNA helices by toposiomerases, which is in accord with previous independent estimations. From the viewpoint of ParB clusters, we predict that they behave like liquid-like protein condensates with unconventional boundaries. Namely, we predict boundaries to be leaky (i.e. not sharp) as a result of the non-equilibrium protein production, diffusion and dilution. Altogether, our work provides novel insights into bacterial DNA organization and intracellular liquid-liquid phase separation.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Transcriptional Regulation Buffers Gene Dosage Effects on a Highly Expressed Operon in Salmonella

A feature of bacterial chromosomes is that highly expressed essential genes are usually located close to the origin of replication. Because bacteria have overlapping cycles of replication, genes located close to the origin will often be present in multiple copies, and this is thought to be of selective benefit where high levels of expression support high growth rate. However, the magnitude of this selective effect and whether other forces could be at play are poorly understood. To study this, we translocated a highly expressed essential operon, tufB, to different locations and measured growth fitness. We found that transcriptional regulation buffered the effects of translocation and that even under conditions where growth rate was reduced, genetic changes that increased the expression of tufB were easily and rapidly selected. We conclude, at least for tufB, that forces other than gene dosage may be significant in selecting for chromosomal location.

Highly expressed genes are commonly located close to the origin of replication of bacterial chromosomes (OriC). This location skew is thought to reflect selective advantages associated with gene dosage effects during the replication cycle. The expression of constitutively expressed genes can vary up to fivefold based on chromosomal location, but it is not clear what level of variation would occur in naturally regulated operons. We tested the magnitude of the chromosome location effect using EF-Tu (tufA, tufB), an abundant protein whose cellular level correlates with, and limits, the maximum growth rate. We translocated the Salmonella tufB operon to four locations across the chromosome. The distance from OriC had only a small effect on growth rate, consistent with this operon having the natural ability to upregulate expression and compensate for reduced gene dosage. In contrast, when the total EF-Tu concentration was limiting for the growth rate (tufA deleted), we observed a strong gene dosage effect when tufB was located further from OriC. However, only a short period of experimental evolution was required before the bacteria adapted to this EF-Tu starvation situation by acquiring genetic changes that increased expression levels from the translocated tufB gene, restoring growth rates. Our findings demonstrate that, at least for the tufB operon, gene dosage is probably not the dominant force selecting for a chromosomal location close to OriC. We suggest that the colocation of highly expressed genes close to OriC might instead be selected because it enhances their coregulation during various growth states, with gene dosage being a secondary benefit.

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.

Insights into the Mechanisms of Basal Coordination of Transcription Using a Genome-Reduced Bacterium

Coordination of transcription in bacteria occurs at supra-operonic scales, but the extent, specificity, and mechanisms of such regulation are poorly understood. Here, we tackle this problem by profiling the transcriptome of the model organism Mycoplasma pneumoniae across 115 growth conditions. We identify three qualitatively different levels of co-expression corresponding to distinct relative orientations and intergenic properties of adjacent genes. We reveal that the degree of co-expression between co-directional adjacent operons, and more generally between genes, is tightly related to their capacity to be transcribed en bloc into the same mRNA. We further show that this genome-wide pervasive transcription of adjacent genes and operons is specifically repressed by DNA regions preferentially bound by RNA polymerases, by intrinsic terminators, and by large intergenic distances. Taken together, our findings suggest that the basal coordination of transcription is mediated by the physical entities and mechanical properties of the transcription process itself, and that operon-like behaviors may strongly vary from condition to condition.

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

From structure to function of bacterial chromosomes: Evolutionary perspectives and ideas for new experiments

The link between chromosome structure and function is a challenging open question because chromosomes in vivo are highly dynamic and arduous to manipulate. Here, we examine several promising approaches to tackle this question specifically in bacteria, by integrating knowledge from different sources. Toward this end, we first provide a brief overview of experimental tools that have provided insights into the description of the bacterial chromosome, including genetic, biochemical and fluorescence microscopy techniques. We then explore the possibility of using comparative genomics to isolate functionally important features of chromosome organization, exploiting the fact that features shared between phylogenetically distant bacterial species reflect functional significance. Finally, we discuss possible future perspectives from the field of experimental evolution. Specifically, we propose novel experiments in which bacteria could be screened and selected on the basis of the structural properties of their chromosomes.

Principles for the organization of gene-sets

A gene-set, an important concept in microarray expression analysis and systems biology, is a collection of genes and/or their products (i.e. proteins) that have some features in common. There are many different ways to construct gene-sets, but a systematic organization of these ways is lacking. Gene-sets are mainly organized ad hoc in current public-domain databases, with group header names often determined by practical reasons (such as the types of technology in obtaining the gene-sets or a balanced number of gene-sets under a header). Here we aim at providing a gene-set organization principle according to the level at which genes are connected: homology, physical map proximity, chemical interaction, biological, and phenotypic-medical levels. We also distinguish two types of connections between genes: actual connection versus sharing of a label. Actual connections denote direct biological interactions, whereas shared label connection denotes shared membership in a group. Some extensions of the framework are also addressed such as overlapping of gene-sets, modules, and the incorporation of other non-protein-coding entities such as microRNAs.