The organization of bacterial genomes: Towards understanding the interplay between structure and function

Machine learning classification of trajectories from molecular dynamics simulations of chromosome segregation

In contrast to the well characterized mitotic machinery in eukaryotes it seems as if there is no universal mechanism organizing chromosome segregation in all bacteria. Apparently, some bacteria even use combinations of different segregation mechanisms such as protein machines or rely on physical forces. The identification of the relevant mechanisms is a difficult task. Here, we introduce a new machine learning approach to this problem. It is based on the analysis of trajectories of individual loci in the course of chromosomal segregation obtained by fluorescence microscopy. While machine learning approaches have already been applied successfully to trajectory classification in other areas, so far it has not been possible to use them to discriminate segregation mechanisms in bacteria. A main obstacle for this is the large number of trajectories required to train machine learning algorithms that we overcome here by using trajectories obtained from molecular dynamics simulations. We used these trajectories to train four different machine learning algorithms, two linear models and two tree-based classifiers, to discriminate segregation mechanisms and possible combinations of them. The classification was performed once using the complete trajectories as high-dimensional input vectors as well as on a set of features which were used to transform the trajectories into low-dimensional input vectors for the classifiers. Finally, we tested our classifiers on shorter trajectories with duration times comparable (or even shorter) than typical experimental trajectories and on trajectories measured with varying temporal resolutions. Our results demonstrate that machine learning algorithms are indeed capable of discriminating different segregation mechanisms in bacteria and to even resolve combinations of the mechanisms on rather short time scales.

Special Issue: Role of Bacterial Chromatin in Environmental Sensing, Adaptation and Evolution

A typical bacterial cell is micron-sized and contains a genome several million base pairs in length [...].

Vibrionaceae core, shell and cloud genes are non-randomly distributed on Chr 1: An hypothesis that links the genomic location of genes with their intracellular placement

Background

The genome of Vibrionaceae bacteria, which consists of two circular chromosomes, is replicated in a highly ordered fashion. In fast-growing bacteria, multifork replication results in higher gene copy numbers and increased expression of genes located close to the origin of replication of Chr 1 (ori1). This is believed to be a growth optimization strategy to satisfy the high demand of essential growth factors during fast growth. The relationship between ori1-proximate growth-related genes and gene expression during fast growth has been investigated by many researchers. However, it remains unclear which other gene categories that are present close to ori1 and if expression of all ori1-proximate genes is increased during fast growth, or if expression is selectively elevated for certain gene categories.

Results

We calculated the pangenome of all complete genomes from the Vibrionaceae family and mapped the four pangene categories, core, softcore, shell and cloud, to their chromosomal positions. This revealed that core and softcore genes were found heavily biased towards ori1, while shell genes were overrepresented at the opposite part of Chr 1 (i.e., close to ter1). RNA-seq of Aliivibrio salmonicida and Vibrio natriegens showed global gene expression patterns that consistently correlated with chromosomal distance to ori1. Despite a biased gene distribution pattern, all pangene categories contributed to a skewed expression pattern at fast-growing conditions, whereas at slow-growing conditions, softcore, shell and cloud genes were responsible for elevated expression.

Conclusion

The pangene categories were non-randomly organized on Chr 1, with an overrepresentation of core and softcore genes around ori1, and overrepresentation of shell and cloud genes around ter1. Furthermore, we mapped our gene distribution data on to the intracellular positioning of chromatin described for V. cholerae, and found that core/softcore and shell/cloud genes appear enriched at two spatially separated intracellular regions. Based on these observations, we hypothesize that there is a link between the genomic location of genes and their cellular placement.

The functional determinants in the organization of bacterial genomes

Bacterial genomes are now recognized as interacting intimately with cellular processes. Uncovering organizational mechanisms of bacterial genomes has been a primary focus of researchers to reveal the potential cellular activities. The advances in both experimental techniques and computational models provide a tremendous opportunity for understanding these mechanisms, and various studies have been proposed to explore the organization rules of bacterial genomes associated with functions recently. This review focuses mainly on the principles that shape the organization of bacterial genomes, both locally and globally. We first illustrate local structures as operons/transcription units for facilitating co-transcription and horizontal transfer of genes. We then clarify the constraints that globally shape bacterial genomes, such as metabolism, transcription and replication. Finally, we highlight challenges and opportunities to advance bacterial genomic studies and provide application perspectives of genome organization, including pathway hole assignment and genome assembly and understanding disease mechanisms.

Different Proteins Mediate Step-Wise Chromosome Architectures in Thermoplasma acidophilum and Pyrobaculum calidifontis

Archaeal species encode a variety of distinct lineage-specific chromosomal proteins. We have previously shown that in Thermococcus kodakarensis, histone, Alba, and TrmBL2 play distinct roles in chromosome organization. Although our understanding of individual archaeal chromosomal proteins has been advancing, how archaeal chromosomes are folded into higher-order structures and how they are regulated are largely unknown. Here, we investigated the primary and higher-order structures of archaeal chromosomes from different archaeal lineages. Atomic force microscopy of chromosome spreads out of Thermoplasma acidophilum and Pyrobaculum calidifontis cells revealed 10-nm fibers and 30–40-nm globular structures, suggesting the occurrence of higher-order chromosomal folding. Our results also indicated that chromosome compaction occurs toward the stationary phase. Micrococcal nuclease digestion indicated that fundamental structural units of the chromosome exist in T. acidophilum and T. kodakarensis but not in P. calidifontis or Sulfolobus solfataricus. In vitro reconstitution showed that, in T. acidophilum, the bacterial HU protein homolog HTa formed a 6-nm fiber by wrapping DNA, and that Alba was responsible for the formation of the 10-nm fiber by binding along the DNA without wrapping. Remarkably, Alba could form different higher-order complexes with histone or HTa on DNA in vitro. Mass spectrometry detected HTa and Rad50 in the T. acidophilum chromosome but not in other species. A putative transcriptional regulator of the AsnC/Lrp family (Pcal_1183) was detected on the P. calidifontis chromosome, but not on that of other species studied. Putative membrane-associated proteins were detected in the chromosomes of the three archaeal species studied, including T. acidophilum, P. calidifontis, and T. kodakarensis. Collectively, our data show that Archaea use different combinations of proteins to achieve chromosomal architecture and functional regulation.

Cross-linker mediated compaction and local morphologies in a model chromosome

Chromatin and associated proteins constitute the highly folded structure of chromosomes. We consider a self-avoiding polymer model of the chromatin, segments of which may get cross-linked via protein binders that repel each other. The binders cluster together via the polymer mediated attraction, in turn, folding the polymer. Using molecular dynamics simulations, and a mean field description, we explicitly demonstrate the continuous nature of the folding transition, characterized by unimodal distributions of the polymer size across the transition. At the transition point the chromatin size and cross-linker clusters display large fluctuations, and a maximum in their negative cross-correlation, apart from a critical slowing down. Along the transition, we distinguish the local chain morphologies in terms of topological loops, inter-loop gaps, and zippering. The topologies are dominated by simply connected loops at the criticality, and by zippering in the folded phase.

Bacterial chromosome organization. II. Few special cross-links, cell confinement, and molecular crowders play the pivotal roles

Using a coarse-grained bead-spring model of bacterial chromosomes of Caulobacter crescentus and Escherichia coli, we show that just 33 and 38 effective cross-links in 4017 and 4642 monomer chains at special positions along the chain contour can lead to the large-scale organization of the DNA polymer, where confinement effects of the cell walls play a key role in the organization. The positions of the 33/38 cross-links along the chain contour are chosen from the Hi-C contact map of bacteria C. crescentus and E. coli. We represent 1000 base pairs as a coarse-grained monomer in our bead-spring flexible ring polymer model of the DNA polymer. Thus, 4017/4642 beads on a flexible ring polymer represent the C. crescentus/E. coli DNA polymer with 4017/4642 kilo-base pairs. Choosing suitable parameters from Paper I, we also incorporate the role of compaction of the polymer coil due to the presence of molecular crowders and the ability of the chain to release topological constraints. We validate our prediction of the organization of the bacterial chromosomes with available experimental data and also give a prediction of the approximate positions of different segments within the cell. In the absence of confinement, the minimal number of effective cross-links required to organize the DNA chains of 4017/4642 monomers was 60/82 [Agarwal et al., Europhys. Lett. 121, 18004 (2018) and Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018)].

Bacterial chromosome organization. I. Crucial role of release of topological constraints and molecular crowders

We showed in our previous studies that just 3% cross-links (CLs), at special points along the contour of the bacterial DNA, help the DNA-polymer to get organized at micron length scales [T. Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018) and T. Agarwal et al., EPL (Europhys. Lett.) 121, 18004 (2018)]. In this work, we investigate how does the release of topological constraints help in the "organization" of the DNA-polymer. Furthermore, we show that the chain compaction induced by the crowded environment in the bacterial cytoplasm contributes to the organization of the DNA-polymer. We model the DNA chain as a flexible bead-spring ring polymer, where each bead represents 1000 base pairs. The specific positions of the CLs have been taken from the experimental contact maps of the bacteria Caulobacter crescentus and Escherichia coli. We introduce different extents of ease of release of topological constraints in our model by systematically changing the diameter of the monomer bead. It varies from the value where chain crossing can occur freely to the value where chain crossing is disallowed. We also study the role of compaction of the chain due to molecular crowders by introducing an "effective" weak Lennard-Jones attraction between the monomers. Using Monte Carlo simulations, we show that the release of topological constraints and the crowding environment play a crucial role to obtain a unique organization of the polymer.

Confinement and crowding control the morphology and dynamics of a model bacterial chromosome

Motivated by recent experiments probing the shape, size and dynamics of bacterial chromosomes in growing cells, we consider a polymer model consisting of a circular backbone to which side-loops are attached, confined to a cylindrical cell. Such a model chromosome spontaneously adopts a helical shape, which is further compacted by molecular crowders to occupy a nucleoid-like sub-volume of the cell. With increasing cell length, the longitudinal size of the chromosome increases in a non-linear fashion until finally saturating, its morphology gradually opening up while displaying a changing number of helical turns. For shorter cells, the chromosome extension varies non-monotonically with cell size, which we show is associated with a radial to longitudinal spatial reordering of the crowders. Confinement and crowders constrain chain dynamics leading to anomalous diffusion. While the scaling exponent for the mean squared displacement of center of mass grows and saturates with cell length, that of individual loci displays a broad distribution with a sharp maximum.