The three-dimensional architecture of a bacterial genome and its alteration by genetic perturbation.

We have determined the three-dimensional (3D) architecture of the Caulobacter crescentus genome by combining genome-wide chromatin interaction detection, live-cell imaging, and computational modeling. Using chromosome conformation capture carbon copy (5C), we derive ~13 kb resolution 3D models of the Caulobacter genome. The resulting models illustrate that the genome is ellipsoidal with periodically arranged arms. The parS sites, a pair of short contiguous sequence elements known to be involved in chromosome segregation, are positioned at one pole, where they anchor the chromosome to the cell and contribute to the formation of a compact chromatin conformation. Repositioning these elements resulted in rotations of the chromosome that changed the subcellular positions of most genes. Such rotations did not lead to large-scale changes in gene expression, indicating that genome folding does not strongly affect gene regulation. Collectively, our data suggest that genome folding is globally dictated by the parS sites and chromosome segregation.

Cell cycle progression in Caulobacter requires a nucleoid-associated protein with high AT sequence recognition.

Faithful cell cycle progression in the dimorphic bacterium Caulobacter crescentus requires spatiotemporal regulation of gene expression and cell pole differentiation. We discovered an essential DNA-associated protein, GapR, that is required for Caulobacter growth and asymmetric division. GapR interacts with adenine and thymine (AT)-rich chromosomal loci, associates with the promoter regions of cell cycle-regulated genes, and shares hundreds of recognition sites in common with known master regulators of cell cycle-dependent gene expression. GapR target loci are especially enriched in binding sites for the transcription factors GcrA and CtrA and overlap with nearly all of the binding sites for MucR1, a regulator that controls the establishment of swarmer cell fate. Despite constitutive synthesis, GapR accumulates preferentially in the swarmer compartment of the predivisional cell. Homologs of GapR, which are ubiquitous among the α-proteobacteria and are encoded on multiple bacteriophage genomes, also accumulate in the predivisional cell swarmer compartment when expressed in Caulobacter The Escherichia coli nucleoid-associated protein H-NS, like GapR, selectively associates with AT-rich DNA, yet it does not localize preferentially to the swarmer compartment when expressed exogenously in Caulobacter, suggesting that recognition of AT-rich DNA is not sufficient for the asymmetric accumulation of GapR. Further, GapR does not silence the expression of H-NS target genes when expressed in E. coli, suggesting that GapR and H-NS have distinct functions. We propose that Caulobacter has co-opted a nucleoid-associated protein with high AT recognition to serve as a mediator of cell cycle progression.

Dynamic translation regulation in Caulobacter cell cycle control.

Progression of the Caulobacter cell cycle requires temporal and spatial control of gene expression, culminating in an asymmetric cell division yielding distinct daughter cells. To explore the contribution of translational control, RNA-seq and ribosome profiling were used to assay global transcription and translation levels of individual genes at six times over the cell cycle. Translational efficiency (TE) was used as a metric for the relative rate of protein production from each mRNA. TE profiles with similar cell cycle patterns were found across multiple clusters of genes, including those in operons or in subsets of operons. Collections of genes associated with central cell cycle functional modules (e.g., biosynthesis of stalk, flagellum, or chemotaxis machinery) have consistent but different TE temporal patterns, independent of their operon organization. Differential translation of operon-encoded genes facilitates precise cell cycle-timing for the dynamic assembly of multiprotein complexes, such as the flagellum and the stalk and the correct positioning of regulatory proteins to specific cell poles. The cell cycle-regulatory pathways that produce specific temporal TE patterns are separate from-but highly coordinated with-the transcriptional cell cycle circuitry, suggesting that the scheduling of translational regulation is organized by the same cyclical regulatory circuit that directs the transcriptional control of the Caulobacter cell cycle.

CauloBrowser: A systems biology resource for Caulobacter crescentus.

Caulobacter crescentus is a premier model organism for studying the molecular basis of cellular asymmetry. The Caulobacter community has generated a wealth of high-throughput spatiotemporal databases including data from gene expression profiling experiments (microarrays, RNA-seq, ChIP-seq, ribosome profiling, LC-ms proteomics), gene essentiality studies (Tn-seq), genome wide protein localization studies, and global chromosome methylation analyses (SMRT sequencing). A major challenge involves the integration of these diverse data sets into one comprehensive community resource. To address this need, we have generated CauloBrowser (www.caulobrowser.org), an online resource for Caulobacter studies. This site provides a user-friendly interface for quickly searching genes of interest and downloading genome-wide results. Search results about individual genes are displayed as tables, graphs of time resolved expression profiles, and schematics of protein localization throughout the cell cycle. In addition, the site provides a genome viewer that enables customizable visualization of all published high-throughput genomic data. The depth and diversity of data sets collected by the Caulobacter community makes CauloBrowser a unique and valuable systems biology resource.

The global regulatory architecture of transcription during the Caulobacter cell cycle.

Each Caulobacter cell cycle involves differentiation and an asymmetric cell division driven by a cyclical regulatory circuit comprised of four transcription factors (TFs) and a DNA methyltransferase. Using a modified global 5' RACE protocol, we globally mapped transcription start sites (TSSs) at base-pair resolution, measured their transcription levels at multiple times in the cell cycle, and identified their transcription factor binding sites. Out of 2726 TSSs, 586 were shown to be cell cycle-regulated and we identified 529 binding sites for the cell cycle master regulators. Twenty-three percent of the cell cycle-regulated promoters were found to be under the combinatorial control of two or more of the global regulators. Previously unknown features of the core cell cycle circuit were identified, including 107 antisense TSSs which exhibit cell cycle-control, and 241 genes with multiple TSSs whose transcription levels often exhibited different cell cycle timing. Cumulatively, this study uncovered novel new layers of transcriptional regulation mediating the bacterial cell cycle.

The coding and noncoding architecture of the Caulobacter crescentus genome.

Caulobacter crescentus undergoes an asymmetric cell division controlled by a genetic circuit that cycles in space and time. We provide a universal strategy for defining the coding potential of bacterial genomes by applying ribosome profiling, RNA-seq, global 5'-RACE, and liquid chromatography coupled with tandem mass spectrometry (LC-MS) data to the 4-megabase C. crescentus genome. We mapped transcript units at single base-pair resolution using RNA-seq together with global 5'-RACE. Additionally, using ribosome profiling and LC-MS, we mapped translation start sites and coding regions with near complete coverage. We found most start codons lacked corresponding Shine-Dalgarno sites although ribosomes were observed to pause at internal Shine-Dalgarno sites within the coding DNA sequence (CDS). These data suggest a more prevalent use of the Shine-Dalgarno sequence for ribosome pausing rather than translation initiation in C. crescentus. Overall 19% of the transcribed and translated genomic elements were newly identified or significantly improved by this approach, providing a valuable genomic resource to elucidate the complete C. crescentus genetic circuitry that controls asymmetric cell division.

The functions of DNA methylation by CcrM in Caulobacter crescentus: a global approach.

DNA methylation is involved in a diversity of processes in bacteria, including maintenance of genome integrity and regulation of gene expression. Here, using Caulobacter crescentus as a model, we exploit genome-wide experimental methods to uncover the functions of CcrM, a DNA methyltransferase conserved in most Alphaproteobacteria. Using single molecule sequencing, we provide evidence that most CcrM target motifs (GANTC) switch from a fully methylated to a hemi-methylated state when they are replicated, and back to a fully methylated state at the onset of cell division. We show that DNA methylation by CcrM is not required for the control of the initiation of chromosome replication or for DNA mismatch repair. By contrast, our transcriptome analysis shows that >10% of the genes are misexpressed in cells lacking or constitutively over-expressing CcrM. Strikingly, GANTC methylation is needed for the efficient transcription of dozens of genes that are essential for cell cycle progression, in particular for DNA metabolism and cell division. Many of them are controlled by promoters methylated by CcrM and co-regulated by other global cell cycle regulators, demonstrating an extensive cross talk between DNA methylation and the complex regulatory network that controls the cell cycle of C. crescentus and, presumably, of many other Alphaproteobacteria.

Caulobacter chromosome in vivo configuration matches model predictions for a supercoiled polymer in a cell-like confinement.

We measured the distance between fluorescent-labeled DNA loci of various interloci contour lengths in Caulobacter crescentus swarmer cells to determine the in vivo configuration of the chromosome. For DNA segments less than about 300 kb, the mean interloci distances, , scale as n(0.22), where n is the contour length, and cell-to-cell distribution of the interloci distance r is a universal function of r/n(0.22) with broad cell-to-cell variability. For DNA segments greater than about 300 kb, the mean interloci distances scale as n, in agreement with previous observations. The 0.22 value of the scaling exponent for short DNA segments is consistent with theoretical predictions for a branched DNA polymer structure. Predictions from Brownian dynamics simulations of the packing of supercoiled DNA polymers in an elongated cell-like confinement are also consistent with a branched DNA structure, and simulated interloci distance distributions predict that confinement leads to "freezing" of the supercoiled configuration. Lateral positions of labeled loci at comparable positions along the length of the cell are strongly correlated when the longitudinal locus positions differ by <0.16 µm. We conclude that the chromosome structure is supercoiled locally and elongated at large length scales and that substantial cell-to-cell variability in the interloci distances indicates that in vivo crowding prevents the chromosome from reaching an equilibrium arrangement. We suggest that the force causing rapid transport of loci remote from the parS centromere to the distal cell pole may arise from the release at the polar region of potential energy within the supercoiled DNA.

Global methylation state at base-pair resolution of the Caulobacter genome throughout the cell cycle.

The Caulobacter DNA methyltransferase CcrM is one of five master cell-cycle regulators. CcrM is transiently present near the end of DNA replication when it rapidly methylates the adenine in hemimethylated GANTC sequences. The timing of transcription of two master regulator genes and two cell division genes is controlled by the methylation state of GANTC sites in their promoters. To explore the global extent of this regulatory mechanism, we determined the methylation state of the entire chromosome at every base pair at five time points in the cell cycle using single-molecule, real-time sequencing. The methylation state of 4,515 GANTC sites, preferentially positioned in intergenic regions, changed progressively from full to hemimethylation as the replication forks advanced. However, 27 GANTC sites remained unmethylated throughout the cell cycle, suggesting that these protected sites could participate in epigenetic regulatory functions. An analysis of the time of activation of every cell-cycle regulatory transcription start site, coupled to both the position of a GANTC site in their promoter regions and the time in the cell cycle when the GANTC site transitions from full to hemimethylation, allowed the identification of 59 genes as candidates for epigenetic regulation. In addition, we identified two previously unidentified N(6)-methyladenine motifs and showed that they maintained a constant methylation state throughout the cell cycle. The cognate methyltransferase was identified for one of these motifs as well as for one of two 5-methylcytosine motifs.

Direct inference of protein-DNA interactions using compressed sensing methods.

Compressed sensing has revolutionized signal acquisition, by enabling complex signals to be measured with remarkable fidelity using a small number of so-called incoherent sensors. We show that molecular interactions, e.g., protein-DNA interactions, can be analyzed in a directly analogous manner and with similarly remarkable results. Specifically, mesoscopic molecular interactions act as incoherent sensors that measure the energies of microscopic interactions between atoms. We combine concepts from compressed sensing and statistical mechanics to determine the interatomic interaction energies of a molecular system exclusively from experimental measurements, resulting in a "de novo" energy potential. In contrast, conventional methods for estimating energy potentials are based on theoretical models premised on a priori assumptions and extensive domain knowledge. We determine the de novo energy potential for pairwise interactions between protein and DNA atoms from (i) experimental measurements of the binding affinity of protein-DNA complexes and (ii) crystal structures of the complexes. We show that the de novo energy potential can be used to predict the binding specificity of proteins to DNA with approximately 90% accuracy, compared to approximately 60% for the best performing alternative computational methods applied to this fundamental problem. This de novo potential method is directly extendable to other biomolecule interaction domains (enzymes and signaling molecule interactions) and to other classes of molecular interactions.

Three enhancements to the inference of statistical protein-DNA potentials.

The energetics of protein-DNA interactions are often modeled using so-called statistical potentials, that is, energy models derived from the atomic structures of protein-DNA complexes. Many statistical protein-DNA potentials based on differing theoretical assumptions have been investigated, but little attention has been paid to the types of data and the parameter estimation process used in deriving the statistical potentials. We describe three enhancements to statistical potential inference that significantly improve the accuracy of predicted protein-DNA interactions: (i) incorporation of binding energy data of protein-DNA complexes, in conjunction with their X-ray crystal structures, (ii) use of spatially-aware parameter fitting, and (iii) use of ensemble-based parameter fitting. We apply these enhancements to three widely-used statistical potentials and use the resulting enhanced potentials in a structure-based prediction of the DNA binding sites of proteins. These enhancements are directly applicable to all statistical potentials used in protein-DNA modeling, and we show that they can improve the accuracy of predicted DNA binding sites by up to 21%.

An essential transcription factor, SciP, enhances robustness of Caulobacter cell cycle regulation.

A cyclical control circuit composed of four master regulators drives the Caulobacter cell cycle. We report that SciP, a helix-turn-helix transcription factor, is an essential component of this circuit. SciP is cell cycle-controlled and co-conserved with the global cell cycle regulator CtrA in the α-proteobacteria. SciP is expressed late in the cell cycle and accumulates preferentially in the daughter swarmer cell. At least 58 genes, including many flagellar and chemotaxis genes, are regulated by a type 1 incoherent feedforward motif in which CtrA activates sciP, followed by SciP repression of ctrA and CtrA target genes. We demonstrate that SciP binds to DNA at a motif distinct from the CtrA binding motif that is present in the promoters of genes co-regulated by SciP and CtrA. SciP overexpression disrupts the balance between activation and repression of the CtrA-SciP coregulated genes yielding filamentous cells and loss of viability. The type 1 incoherent feedforward circuit motif enhances the pulse-like expression of the downstream genes, and the negative feedback to ctrA expression reduces peak CtrA accumulation. The presence of SciP in the control network enhances the robustness of the cell cycle to varying growth rates.

High-throughput identification of protein localization dependency networks.

Bacterial cells are highly organized with many protein complexes and DNA loci dynamically positioned to distinct subcellular sites over the course of a cell cycle. Such dynamic protein localization is essential for polar organelle development, establishment of asymmetry, and chromosome replication during the Caulobacter crescentus cell cycle. We used a fluorescence microscopy screen optimized for high-throughput to find strains with anomalous temporal or spatial protein localization patterns in transposon-generated mutant libraries. Automated image acquisition and analysis allowed us to identify genes that affect the localization of two polar cell cycle histidine kinases, PleC and DivJ, and the pole-specific pili protein CpaE, each tagged with a different fluorescent marker in a single strain. Four metrics characterizing the observed localization patterns of each of the three labeled proteins were extracted for hundreds of cell images from each of 854 mapped mutant strains. Using cluster analysis of the resulting set of 12-element vectors for each of these strains, we identified 52 strains with mutations that affected the localization pattern of the three tagged proteins. This information, combined with quantitative localization data from epitasis experiments, also identified all previously known proteins affecting such localization. These studies provide insights into factors affecting the PleC/DivJ localization network and into regulatory links between the localization of the pili assembly protein CpaE and the kinase localization pathway. Our high-throughput screening methodology can be adapted readily to any sequenced bacterial species, opening the potential for databases of localization regulatory networks across species, and investigation of localization network phylogenies.

Compaction and transport properties of newly replicated Caulobacter crescentus DNA.

Upon initiating replication of the Caulobacter chromosome, one copy of the parS centromere remains at the stalked pole; the other moves to the distal pole. We identified the segregation dynamics and compaction characteristics of newly replicated Caulobacter DNA during transport (highly variable from cell to cell) using time-lapse fluorescence microscopy. The parS centromere and a length (also highly variable) of parS proximal DNA on each arm of the chromosome are segregated with the same relatively slow transport pattern as the parS locus. Newly replicated DNA further than about 100 kb from parS segregates with a different and faster pattern, while loci at 48 kb from parS segregate with the slow pattern in some cells and the fast pattern in others. The observed parS-proximal DNA compaction characteristics have scaling properties that suggest the DNA is branched. HU2-deletion strains exhibited a reduced compaction phenotype except near the parS site where only the ΔHU1ΔHU2 double mutant had a compaction phenotype. The chromosome shows speed-dependent extension during translocation suggesting the DNA polymer is under tension. While DNA segregation is highly reliable and succeeds in virtually all wild-type cells, the high degree of cell to cell variation in the segregation process is noteworthy.

Assembly of the Caulobacter cell division machine.

Cytokinesis in Gram-negative bacteria is mediated by a multiprotein machine (the divisome) that invaginates and remodels the inner membrane, peptidoglycan and outer membrane. Understanding the order of divisome assembly would inform models of the interactions among its components and their respective functions. We leveraged the ability to isolate synchronous populations of Caulobacter crescentus cells to investigate assembly of the divisome and place the arrival of each component into functional context. Additionally, we investigated the genetic dependence of localization among divisome proteins and the cell cycle regulation of their transcript and protein levels to gain insight into the control mechanisms underlying their assembly. Our results revealed a picture of divisome assembly with unprecedented temporal resolution. Specifically, we observed (i) initial establishment of the division site, (ii) recruitment of early FtsZ-binding proteins, (iii) arrival of proteins involved in peptidoglycan remodelling, (iv) arrival of FtsA, (v) assembly of core divisome components, (vi) initiation of envelope invagination, (vii) recruitment of polar markers and cytoplasmic compartmentalization and (viii) cell separation. Our analysis revealed differences in divisome assembly among Caulobacter and other bacteria that establish a framework for identifying aspects of bacterial cytokinesis that are widely conserved from those that are more variable.

The essential genome of a bacterium.

Caulobacter crescentus is a model organism for the integrated circuitry that runs a bacterial cell cycle. Full discovery of its essential genome, including non-coding, regulatory and coding elements, is a prerequisite for understanding the complete regulatory network of a bacterial cell. Using hyper-saturated transposon mutagenesis coupled with high-throughput sequencing, we determined the essential Caulobacter genome at 8 bp resolution, including 1012 essential genome features: 480 ORFs, 402 regulatory sequences and 130 non-coding elements, including 90 intergenic segments of unknown function. The essential transcriptional circuitry for growth on rich media includes 10 transcription factors, 2 RNA polymerase sigma factors and 1 anti-sigma factor. We identified all essential promoter elements for the cell cycle-regulated genes. The essential elements are preferentially positioned near the origin and terminus of the chromosome. The high-resolution strategy used here is applicable to high-throughput, full genome essentiality studies and large-scale genetic perturbation experiments in a broad class of bacterial species.

Caulobacter PopZ forms a polar subdomain dictating sequential changes in pole composition and function.

The bacterium Caulobacter crescentus has morphologically and functionally distinct cell poles that undergo sequential changes during the cell cycle. We show that the PopZ oligomeric network forms polar ribosome exclusion zones that change function during cell cycle progression. The parS/ParB chromosomal centromere is tethered to PopZ at one pole prior to the initiation of DNA replication. During polar maturation, the PopZ-centromere tether is broken, and the PopZ zone at that pole then switches function to act as a recruitment factor for the ordered addition of multiple proteins that promote the transformation of the flagellated pole into a stalked pole. Stalked pole assembly, in turn, triggers the initiation of chromosome replication, which signals the formation of a new PopZ zone at the opposite cell pole, where it functions to anchor the newly duplicated centromere that has traversed the long axis of the cell. We propose that pole-specific control of PopZ function co-ordinates polar development and cell cycle progression by enabling independent assembly and tethering activities at the two cell poles.

Caulobacter requires a dedicated mechanism to initiate chromosome segregation.

Chromosome segregation in bacteria is rapid and directed, but the mechanisms responsible for this movement are still unclear. We show that Caulobacter crescentus makes use of and requires a dedicated mechanism to initiate chromosome segregation. Caulobacter has a single circular chromosome whose origin of replication is positioned at one cell pole. Upon initiation of replication, an 8-kb region of the chromosome containing both the origin and parS moves rapidly to the opposite pole. This movement requires the highly conserved ParABS locus that is essential in Caulobacter. We use chromosomal inversions and in vivo time-lapse imaging to show that parS is the Caulobacter site of force exertion, independent of its position in the chromosome. When parS is moved farther from the origin, the cell waits for parS to be replicated before segregation can begin. Also, a mutation in the ATPase domain of ParA halts segregation without affecting replication initiation. Chromosome segregation in Caulobacter cannot occur unless a dedicated parS guiding mechanism initiates movement.

Architecture and inherent robustness of a bacterial cell-cycle control system.

A closed-loop control system drives progression of the coupled stalked and swarmer cell cycles of the bacterium Caulobacter crescentus in a near-mechanical step-like fashion. The cell-cycle control has a cyclical genetic circuit composed of four regulatory proteins with tight coupling to processive chromosome replication and cell division subsystems. We report a hybrid simulation of the coupled cell-cycle control system, including asymmetric cell division and responses to external starvation signals, that replicates mRNA and protein concentration patterns and is consistent with observed mutant phenotypes. An asynchronous sequential digital circuit model equivalent to the validated simulation model was created. Formal model-checking analysis of the digital circuit showed that the cell-cycle control is robust to intrinsic stochastic variations in reaction rates and nutrient supply, and that it reliably stops and restarts to accommodate nutrient starvation. Model checking also showed that mechanisms involving methylation-state changes in regulatory promoter regions during DNA replication increase the robustness of the cell-cycle control. The hybrid cell-cycle simulation implementation is inherently extensible and provides a promising approach for development of whole-cell behavioral models that can replicate the observed functionality of the cell and its responses to changing environmental conditions.

The caulobacter Tol-Pal complex is essential for outer membrane integrity and the positioning of a polar localization factor.

Cell division in Caulobacter crescentus involves constriction and fission of the inner membrane (IM) followed about 20 min later by fission of the outer membrane (OM) and daughter cell separation. In contrast to Escherichia coli, the Caulobacter Tol-Pal complex is essential. Cryo-electron microscopy images of the Caulobacter cell envelope exhibited outer membrane disruption, and cells failed to complete cell division in TolA, TolB, or Pal mutant strains. In wild-type cells, components of the Tol-Pal complex localize to the division plane in early predivisional cells and remain predominantly at the new pole of swarmer and stalked progeny upon completion of division. The Tol-Pal complex is required to maintain the position of the transmembrane TipN polar marker, and indirectly the PleC histidine kinase, at the cell pole, but it is not required for the polar maintenance of other transmembrane and membrane-associated polar proteins tested. Coimmunoprecipitation experiments show that both TolA and Pal interact directly or indirectly with TipN. We propose that disruption of the trans-envelope Tol-Pal complex releases TipN from its subcellular position. The Caulobacter Tol-Pal complex is thus a key component of cell envelope structure and function, mediating OM constriction at the final step of cell division as well as the positioning of a protein localization factor.