A Bacterial Cell-Cycle Regulatory Network Operating in Time and Space

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

DeepVelo: Single-cell transcriptomic deep velocity field learning with neural ordinary differential equations

Recent advances in single-cell sequencing technologies have provided unprecedented opportunities to measure the gene expression profile and RNA velocity of individual cells. However, modeling transcriptional dynamics is computationally challenging because of the high-dimensional, sparse nature of the single-cell gene expression measurements and the nonlinear regulatory relationships. Here, we present DeepVelo, a neural network-based ordinary differential equation that can model complex transcriptome dynamics by describing continuous-time gene expression changes within individual cells. We apply DeepVelo to public datasets from different sequencing platforms to (i) formulate transcriptome dynamics on different time scales, (ii) measure the instability of cell states, and (iii) identify developmental driver genes via perturbation analysis. Benchmarking against the state-of-the-art methods shows that DeepVelo can learn a more accurate representation of the velocity field. Furthermore, our perturbation studies reveal that single-cell dynamical systems could exhibit chaotic properties. In summary, DeepVelo allows data-driven discoveries of differential equations that delineate single-cell transcriptome dynamics.

Independent control of mean and noise by convolution of gene expression distributions

Gene expression noise can reduce cellular fitness or facilitate processes such as alternative metabolism, antibiotic resistance, and differentiation. Unfortunately, efforts to study the impacts of noise have been hampered by a scaling relationship between noise and expression level from individual promoters. Here, we use theory to demonstrate that mean and noise can be controlled independently by expressing two copies of a gene from separate inducible promoters in the same cell. We engineer low and high noise inducible promoters to validate this result in Escherichia coli, and develop a model that predicts the experimental distributions. Finally, we use our method to reveal that the response of a promoter to a repressor is less sensitive with higher repressor noise and explain this result using a law from probability theory. Our approach can be applied to investigate the effects of noise on diverse biological pathways or program cellular heterogeneity for synthetic biology applications.

Multi-omics Study of Planobispora rosea, Producer of the Thiopeptide Antibiotic GE2270A

Planobispora rosea is the natural producer of the potent thiopeptide antibiotic GE2270A. Here, we present the results of a metabolomics and transcriptomics analysis of P. rosea during production of GE2270A. The data generated provides useful insights into the biology of this genetically intractable bacterium. We characterize the details of the shutdown of protein biosynthesis and the respiratory chain associated with the end of the exponential growth phase. We also provide the first description of the phosphate regulon in P. rosea. Based on the transcriptomics data, we show that both phosphate and iron are limiting P. rosea growth in our experimental conditions. Additionally, we identified and validated a new biosynthetic gene cluster associated with the production of the siderophores benarthin and dibenarthin in P. rosea. Together, the metabolomics and transcriptomics data are used to inform and refine the very first genome-scale metabolic model for P. rosea, which will be a valuable framework for the interpretation of future studies of the biology of this interesting but poorly characterized species.

IMPORTANCEPlanobispora rosea is a genetically intractable bacterium used for the production of GE2270A on an industrial scale. GE2270A is a potent thiopeptide antibiotic currently used as a precursor for the synthesis of two compounds under clinical studies for the treatment of Clostridium difficile infection and acne. Here, we present the very first systematic multi-omics investigation of this important bacterium, which provides a much-needed detailed picture of the dynamics of metabolism of P. rosea while producing GE2270A.

Author Video: An author video summary of this article is available.

The Apparent Requirement for Protein Synthesis during G2 Phase Is due to Checkpoint Activation

Protein synthesis inhibitors (e.g., cycloheximide) block mitotic entry, suggesting that cell cycle progression requires protein synthesis until right before mitosis. However, cycloheximide is also known to activate p38 mitogen-activated protein kinase (MAPK), which can delay mitotic entry through a G2/M checkpoint. Here, we ask whether checkpoint activation or a requirement for protein synthesis is responsible for the cycloheximide effect. We find that p38 inhibitors prevent cycloheximide-treated cells from arresting in G2 phase and that G2 duration is normal in approximately half of these cells. The Wee1 inhibitor MK-1775 and Wee1/Myt1 inhibitor PD0166285 also prevent cycloheximide from blocking mitotic entry, raising the possibility that Wee1 and/or Myt1 mediate the cycloheximide-induced G2 arrest. Thus, protein synthesis during G2 phase is not required for mitotic entry, at least when the p38 checkpoint pathway is abrogated. However, M phase progression is delayed in cycloheximide-plus-kinase-inhibitor-treated cells, emphasizing the different requirements of protein synthesis for timely entry and completion of mitosis.

Protein synthesis inhibitors have long been known to prevent G2 phase cells from entering mitosis. Lockhead et al. demonstrate that this G2 arrest is due to the activation of p38 MAPK, not insufficient protein synthesis, arguing that protein synthesis in G2 phase is not absolutely required for mitotic entry.

Chromatin-remodeling links metabolic signaling to gene expression

BACKGROUND

ATP-dependent chromatin remodelers are evolutionarily conserved complexes that alter nucleosome positioning to influence many DNA-templated processes, such as replication, repair, and transcription. In particular, chromatin remodeling can dynamically regulate gene expression by altering accessibility of chromatin to transcription factors.

SCOPE OF REVIEW

This review provides an overview of the importance of chromatin remodelers in the regulation of metabolic gene expression. Particular emphasis is placed on the INO80 and SWI/SNF (BAF/PBAF) chromatin remodelers in both yeast and mammals. This review details discoveries from the initial identification of chromatin remodelers in Saccharomyces cerevisiae to recent discoveries in the metabolic requirements of developing embryonic tissues in mammals.

MAJOR CONCLUSIONS

INO80 and SWI/SNF (BAF/PBAF) chromatin remodelers regulate the expression of energy metabolism pathways in S. cerevisiae and mammals in response to diverse nutrient environments. In particular, the INO80 complex organizes the temporal expression of gene expression in the metabolically synchronized S. cerevisiae system. INO80-mediated chromatin remodeling is also needed to constrain cell division during metabolically favorable conditions. Conversely, the BAF/PBAF remodeler regulates tissue-specific glycolytic metabolism and is disrupted in cancers that are dependent on glycolysis for proliferation. The role of chromatin remodeling in metabolic gene expression is downstream of the metabolic signaling pathways, such as the TOR pathway, a critical regulator of metabolic homeostasis. Furthermore, the INO80 and BAF/PBAF chromatin remodelers have both been shown to regulate heart development, the tissues of which have unique requirements for energy metabolism during development. Collectively, these results demonstrate that chromatin remodelers communicate metabolic status to chromatin and are a central component of homeostasis pathways that optimize cell fitness, organismal development, and prevent disease.

Caulobacter crescentus β sliding clamp employs a noncanonical regulatory model of DNA replication

The eubacterial β sliding clamp (DnaN) plays a crucial role in DNA metabolism through direct interactions with DNA, polymerases, and a variety of protein factors. A canonical protein-DnaN interaction has been identified in Escherichia coli and some other species, during which protein partners are tethered into the conserved canonical hydrophobic crevice of DnaN via the consensus β-binding motif. Caulobacter crescentus is an excellent research model for use in the investigation of DNA replication and cell-cycle regulation due to its unique asymmetric cell division pattern with restricted replication initiation; however, little is known about the specific features of C. crescentus DnaN (CcDnaN). Here, we report a significant divergence in the association of CcDnaN with proteins based on docking analysis and crystal structures that show that the β-binding motifs of its protein partners bind a novel pocket instead of the canonical site. Pull-down and isothermal titration calorimetry results revealed that mutations within the novel pocket disrupt protein-CcDnaN interactions. It was also shown by replication and regulatory inactivation of DnaA assays that mediation of protein interaction by the novel pocket is closely related to the performance of CcDnaN during replication and the DnaN-mediated regulation process. Moreover, assessments of clamp competition showed that DNA does not compete with protein partners when binding to the novel pocket. Overall, our structural and biochemical analyses provide strong evidence that CcDnaN employs a noncanonical protein association pattern.

INO80 Chromatin Remodeling Coordinates Metabolic Homeostasis with Cell Division

Adaptive survival requires the coordination of nutrient availability with expenditure of cellular resources. For example, in nutrient-limited environments, 50% of all S. cerevisiae genes synchronize and exhibit periodic bursts of expression in coordination with respiration and cell division in the yeast metabolic cycle (YMC). Despite the importance of metabolic and proliferative synchrony, the majority of YMC regulators are currently unknown. Here, we demonstrate that the INO80 chromatin-remodeling complex is required to coordinate respiration and cell division with periodic gene expression. Specifically, INO80 mutants have severe defects in oxygen consumption and promiscuous cell division that is no longer coupled with metabolic status. In mutant cells, chromatin accessibility of periodic genes, including TORC1-responsive genes, is relatively static, concomitant with severely attenuated gene expression. Collectively, these results reveal that the INO80 complex mediates metabolic signaling to chromatin to restrict proliferation to metabolically optimal states.

Chemical synthesis rewriting of a bacterial genome to achieve design flexibility and biological functionality

The fundamental biological functions of a living cell are stored within the DNA sequence of its genome. Classical genetic approaches dissect the functioning of biological systems by analyzing individual genes, yet uncovering the essential gene set of an organism has remained very challenging. It is argued that the rewriting of entire genomes through the process of chemical synthesis provides a powerful and complementary research concept to understand how essential functions are programed into genomes.

Understanding how to program biological functions into artificial DNA sequences remains a key challenge in synthetic genomics. Here, we report the chemical synthesis and testing of Caulobacter ethensis-2.0 (C. eth-2.0), a rewritten bacterial genome composed of the most fundamental functions of a bacterial cell. We rebuilt the essential genome of Caulobacter crescentus through the process of chemical synthesis rewriting and studied the genetic information content at the level of its essential genes. Within the 785,701-bp genome, we used sequence rewriting to reduce the number of encoded genetic features from 6,290 to 799. Overall, we introduced 133,313 base substitutions, resulting in the rewriting of 123,562 codons. We tested the biological functionality of the genome design in C. crescentus by transposon mutagenesis. Our analysis revealed that 432 essential genes of C. eth-2.0, corresponding to 81.5% of the design, are equal in functionality to natural genes. These findings suggest that neither changing mRNA structure nor changing the codon context have significant influence on biological functionality of synthetic genomes. Discovery of 98 genes that lost their function identified essential genes with incorrect annotation, including a limited set of 27 genes where we uncovered noncoding control features embedded within protein-coding sequences. In sum, our results highlight the promise of chemical synthesis rewriting to decode fundamental genome functions and its utility toward the design of improved organisms for industrial purposes and health benefits.

Programming Bacteria With Light-Sensors and Applications in Synthetic Biology

Photo-receptors are widely present in both prokaryotic and eukaryotic cells, which serves as the foundation of tuning cell behaviors with light. While practices in eukaryotic cells have been relatively established, trials in bacterial cells have only been emerging in the past few years. A number of light sensors have been engineered in bacteria cells and most of them fall into the categories of two-component and one-component systems. Such a sensor toolbox has enabled practices in controlling synthetic circuits at the level of transcription and protein activity which is a major topic in synthetic biology, according to the central dogma. Additionally, engineered light sensors and practices of tuning synthetic circuits have served as a foundation for achieving light based real-time feedback control. Here, we review programming bacteria cells with light, introducing engineered light sensors in bacteria and their applications, including tuning synthetic circuits and achieving feedback controls over microbial cell culture.

LexA Binds to Transcription Regulatory Site of Cell Division Gene ftsZ in Toxic Cyanobacterium Microcystis aeruginosa

Previously, we showed that DNA replication and cell division in toxic cyanobacterium Microcystis aeruginosa are coordinated by transcriptional regulation of cell division gene ftsZ and that an unknown protein specifically bound upstream of ftsZ (BpFz; DNA-binding protein to an upstream site of ftsZ) during successful DNA replication and cell division. Here, we purified BpFz from M. aeruginosa strain NIES-298 using DNA-affinity chromatography and gel-slicing combined with gel electrophoresis mobility shift assay (EMSA). The N-terminal amino acid sequence of BpFz was identified as TNLESLTQ, which was identical to that of transcription repressor LexA from NIES-843. EMSA analysis using mutant probes showed that the sequence GTACTAN₃GTGTTC was important in LexA binding. Comparison of the upstream regions of lexA in the genomes of closely related cyanobacteria suggested that the sequence TASTRNNNNTGTWC could be a putative LexA recognition sequence (LexA box). Searches for TASTRNNNNTGTWC as a transcriptional regulatory site (TRS) in the genome of M. aeruginosa NIES-843 showed that it was present in genes involved in cell division, photosynthesis, and extracellular polysaccharide biosynthesis. Considering that BpFz binds to the TRS of ftsZ during normal cell division, LexA may function as a transcriptional activator of genes related to cell reproduction in M. aeruginosa, including ftsZ. This may be an example of informality in the control of bacterial cell division.

Influences of adding negative couplings between cliques of Kuramoto-like oscillators

We study the dynamics in a clustered network of coupled oscillators by considering positive and negative coupling schemes. Second order oscillators can be interpreted as a model of consumers and generators working in a power network. Numerical results indicate that coupling strategies play an important role in the synchronizability of the clustered power network. It is found that the synchronizability can be enhanced as the positive intragroup connections increase. Meanwhile, when the intragroup interactions are positive and the probability p that two nodes belonging to different clusters are connected is increased, the synchronization has better performance. Besides, when the intragroup connections are negative, it is observed that the power network has poor synchronizability as the probability p increases. Our simulation results can help us understand the collective behavior of the power network with positive and negative couplings.

A Single-Domain Response Regulator Functions as an Integrating Hub To Coordinate General Stress Response and Development in Alphaproteobacteria

The alphaproteobacterial general stress response is governed by a conserved partner-switching mechanism that is triggered by phosphorylation of the response regulator PhyR. In the model organism Caulobacter crescentus, PhyR was proposed to be phosphorylated by the histidine kinase PhyK, but biochemical evidence in support of such a role of PhyK is missing. Here, we identify a single-domain response regulator, MrrA, that is essential for general stress response activation in C. crescentus We demonstrate that PhyK does not function as a kinase but accepts phosphoryl groups from MrrA and passes them on to PhyR, adopting the role of a histidine phosphotransferase. MrrA is phosphorylated by at least six histidine kinases that likely serve as stress sensors. MrrA also transfers phosphate to LovK, a histidine kinase involved in C. crescentus holdfast production and attachment, which also negatively regulates the general stress response. We show that LovK together with the response regulator LovR acts as a phosphate sink to redirect phosphate flux away from the PhyKR branch. In agreement with the biochemical data, an mrrA mutant is unable to activate the general stress response and shows a hyperattachment phenotype, which is linked to decreased expression of the major holdfast inhibitory protein HfiA. We propose that MrrA serves as a central phosphorylation hub that coordinates the general stress response with C. crescentus development and other adaptive behaviors. The characteristic bow-tie architecture of this phosphorylation network with MrrA as the central knot may expedite the evolvability and species-specific niche adaptation of this group of bacteria.IMPORTANCE Two-component systems (TCSs) consisting of a histidine kinase and a cognate response regulator are predominant signal transduction systems in bacteria. To avoid cross talk, TCSs are generally thought to be highly insulated from each other. However, this notion is based largely on studies of the HisKA subfamily of histidine kinases, while little information is available for the HWE and HisKA2 subfamilies. The latter have been implicated in the alphaproteobacterial general stress response. Here, we show that in the model organism Caulobacter crescentus an atypical FATGUY-type single-domain response regulator, MrrA, is highly promiscuous in accepting and transferring phosphoryl groups from and to multiple up- and downstream kinases, challenging the current view of strictly insulated TCSs. Instead, we propose that FATGUY response regulators have evolved in alphaproteobacteria as central phosphorylation hubs to broadly sample information and distribute phosphoryl groups between the general stress response pathway and other TCSs, thereby coordinating multiple cellular behaviors.

Genome Partitioner: A web tool for multi-level partitioning of large-scale DNA constructs for synthetic biology applications

Recent advances in lower-cost DNA synthesis techniques have enabled new innovations in the field of synthetic biology. Still, efficient design and higher-order assembly of genome-scale DNA constructs remains a labor-intensive process. Given the complexity, computer assisted design tools that fragment large DNA sequences into fabricable DNA blocks are needed to pave the way towards streamlined assembly of biological systems. Here, we present the Genome Partitioner software implemented as a web-based interface that permits multi-level partitioning of genome-scale DNA designs. Without the need for specialized computing skills, biologists can submit their DNA designs to a fully automated pipeline that generates the optimal retrosynthetic route for higher-order DNA assembly. To test the algorithm, we partitioned a 783 kb Caulobacter crescentus genome design. We validated the partitioning strategy by assembling a 20 kb test segment encompassing a difficult to synthesize DNA sequence. Successful assembly from 1 kb subblocks into the 20 kb segment highlights the effectiveness of the Genome Partitioner for reducing synthesis costs and timelines for higher-order DNA assembly. The Genome Partitioner is broadly applicable to translate DNA designs into ready to order sequences that can be assembled with standardized protocols, thus offering new opportunities to harness the diversity of microbial genomes for synthetic biology applications. The Genome Partitioner web tool can be accessed at https://christenlab.ethz.ch/GenomePartitioner.

A Spatial Control for Correct Timing of Gene Expression during the Escherichia coli Cell Cycle

Temporal transcriptions of genes are achieved by different mechanisms such as dynamic interaction of activator and repressor proteins with promoters, and accumulation and/or degradation of key regulators as a function of cell cycle. We find that the TorR protein localizes to the old poles of the Escherichia coli cells, forming a functional focus. The TorR focus co-localizes with the nucleoid in a cell-cycle-dependent manner, and consequently regulates transcription of a number of genes. Formation of one TorR focus at the old poles of cells requires interaction with the MreB and DnaK proteins, and ATP, suggesting that TorR delivery requires cytoskeleton organization and ATP. Further, absence of the protein–protein interactions and ATP leads to loss in function of TorR as a transcription factor. We propose a mechanism for timing of cell-cycle-dependent gene transcription, where a transcription factor interacts with its target genes during a specific period of the cell cycle by limiting its own spatial distribution.

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.

The Landscape of Evolution: Reconciling Structural and Dynamic Properties of Metabolic Networks in Adaptive Diversifications

The network of the interactions among genes, proteins, and metabolites delineates a range of potential phenotypic diversifications in a lineage, and realized phenotypic changes are the result of differences in the dynamics of the expression of the elements and interactions in this deterministic network. Regulatory mechanisms, such as hormones, mediate the relationship between the structural and dynamic properties of networks by determining how and when the elements are expressed and form a functional unit or state. Changes in regulatory mechanisms lead to variable expression of functional states of a network within and among generations. Functional properties of network elements, and the magnitude and direction of evolutionary change they determine, depend on their location within a network. Here, we examine the relationship between network structure and the dynamic mechanisms that regulate flux through a metabolic network. We review the mechanisms that control metabolic flux in enzymatic reactions and examine structural properties of the network locations that are targets of flux control. We aim to establish a predictive framework to test the contributions of structural and dynamic properties of deterministic networks to evolutionary diversifications.

OmpW of Caulobacter crescentus Functions as an Outer Membrane Channel for Cations

Caulobacter crescentus is an oligotrophic bacterium that lives in dilute organic environments such as soil and freshwater. This bacterium represents an interesting model for cellular differentiation and regulation because daughter cells after division have different forms: one is motile while the other is non-motile and can adhere to surfaces. Interestingly, the known genome of C. crescentus does not contain genes predicted to code for outer membrane porins of the OmpF/C general diffusion type present in enteric bacteria or those coding for specific porins selective for classes of substrates. Instead, genes coding for 67 TonB-dependent outer membrane receptors have been identified, suggesting that active transport of specific nutrients may be the norm. Here, we report that high channel-forming activity was observed with crude outer membrane extracts of C. crescentus in lipid bilayer experiments, indicating that the outer membrane of C. crescentus contained an ion-permeable channel with a single-channel conductance of about 120 pS in 1M KCl. The channel-forming protein with an apparent molecular mass of about 20 kDa was purified to homogeneity. Partial protein sequencing of the protein indicated it was a member of the OmpW family of outer membrane proteins from Gram-negative bacteria. This channel was not observed in reconstitution experiments with crude outer membrane extracts of an OmpW deficient C. crescentus mutant. Biophysical analysis of the C. crescentus OmpW suggested that it has features that are special for general diffusion porins of Gram-negative outer membranes because it was not a wide aqueous channel. Furthermore, OmpW of C. crescentus seems to be different to known OmpW porins and has a preference for ions, in particular cations. A putative model for OmpW of C. crescentus was built on the basis of the known 3D-structures of OmpW of Escherichia coli and OprG of Pseudomonas aeruginosa using homology modeling. A comparison of the two known structures with the model of OmpW of C. crescentus suggested that it has a more hydrophilic interior and possibly a larger diameter.

Genome Calligrapher: A Web Tool for Refactoring Bacterial Genome Sequences for de Novo DNA Synthesis

Recent advances in synthetic biology have resulted in an increasing demand for the de novo synthesis of large-scale DNA constructs. Any process improvement that enables fast and cost-effective streamlining of digitized genetic information into fabricable DNA sequences holds great promise to study, mine, and engineer genomes. Here, we present Genome Calligrapher, a computer-aided design web tool intended for whole genome refactoring of bacterial chromosomes for de novo DNA synthesis. By applying a neutral recoding algorithm, Genome Calligrapher optimizes GC content and removes obstructive DNA features known to interfere with the synthesis of double-stranded DNA and the higher order assembly into large DNA constructs. Subsequent bioinformatics analysis revealed that synthesis constraints are prevalent among bacterial genomes. However, a low level of codon replacement is sufficient for refactoring bacterial genomes into easy-to-synthesize DNA sequences. To test the algorithm, 168 kb of synthetic DNA comprising approximately 20 percent of the synthetic essential genome of the cell-cycle bacterium Caulobacter crescentus was streamlined and then ordered from a commercial supplier of low-cost de novo DNA synthesis. The successful assembly into eight 20 kb segments indicates that Genome Calligrapher algorithm can be efficiently used to refactor difficult-to-synthesize DNA. Genome Calligrapher is broadly applicable to recode biosynthetic pathways, DNA sequences, and whole bacterial genomes, thus offering new opportunities to use synthetic biology tools to explore the functionality of microbial diversity. The Genome Calligrapher web tool can be accessed at https://christenlab.ethz.ch/GenomeCalligrapher  .

Regulation of DNA Replication Initiation by Chromosome Structure

Recent advancements in fluorescence imaging have shown that the bacterial nucleoid is surprisingly dynamic in terms of both behavior (movement and organization) and structure (density and supercoiling). Links between chromosome structure and replication initiation have been made in a number of species, and it is universally accepted that favorable chromosome structure is required for initiation in all cells. However, almost nothing is known about whether cells use changes in chromosome structure as a regulatory mechanism for initiation. Such changes could occur during natural cell cycle or growth phase transitions, or they could be manufactured through genetic switches of topoisomerase and nucleoid structure genes. In this review, we explore the relationship between chromosome structure and replication initiation and highlight recent work implicating structure as a regulatory mechanism. A three-component origin activation model is proposed in which thermal and topological structural elements are balanced with trans-acting control elements (DnaA) to allow efficient initiation control under a variety of nutritional and environmental conditions. Selective imbalances in these components allow cells to block replication in response to cell cycle impasse, override once-per-cell-cycle programming during growth phase transitions, and promote reinitiation when replication forks fail to complete.

Synchronization of Caulobacter crescentus for investigation of the bacterial cell cycle

The cell cycle is important for growth, genome replication, and development in all cells. In bacteria, studies of the cell cycle have focused largely on unsynchronized cells making it difficult to order the temporal events required for cell cycle progression, genome replication, and division. Caulobacter crescentus provides an excellent model system for the bacterial cell cycle whereby cells can be rapidly synchronized in a G0 state by density centrifugation. Cell cycle synchronization experiments have been used to establish the molecular events governing chromosome replication and segregation, to map a genetic regulatory network controlling cell cycle progression, and to identify the establishment of polar signaling complexes required for asymmetric cell division. Here we provide a detailed protocol for the rapid synchronization of Caulobacter NA1000 cells. Synchronization can be performed in a large-scale format for gene expression profiling and western blot assays, as well as a small-scale format for microscopy or FACS assays. The rapid synchronizability and high cell yields of Caulobacter make this organism a powerful model system for studies of the bacterial cell cycle.

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 generation of diverse cell types occurs through two fundamental processes; asymmetric cell division and cell differentiation. Cells progress through these developmental changes guided by complex and layered genetic programs that lead to differential expression of the genome. To explore how a genetic program directs cell cycle progression, we examined the global activity of promoters at distinct stages of the cell cycle of the bacterium Caulobacter crescentus, that undergoes cellular differentiation and divides asymmetrically at each cell division. We found that approximately 21% of transcription start sites are cell cycle-regulated, driving the transcription of both mRNAs and non-coding and antisense RNAs. In addition, 102 cell cycle-regulated genes are transcribed from multiple promoters, allowing multiple regulatory inputs to control the logic of gene activation. We found combinatorial control by the five master transcription regulators that provide the core regulation for the genetic circuitry controlling the cell cycle. Much of this combinatorial control appears to be directed at refinement of temporal expression of various genes over the cell cycle, and at tighter control of asymmetric gene expression between the swarmer and stalked daughter cells.

Maintenance of chromosome structure in Pseudomonas aeruginosa

Replication and segregation of genetic information are the activities central to the well-being of all living cells. Concerted mechanisms have evolved that ensure that each cellular chromosome is replicated once and only once per cell cycle and then faithfully segregated into daughter cells. Despite remarkable taxonomic diversity, these mechanisms are largely conserved across eubacteria, although species-specific distinctions can often be noted. Here, we provide an overview of the current state of knowledge about maintenance of the chromosome structure in Pseudomonas aeruginosa. We focus on global chromosome organization and its dynamics during DNA replication and cell division. Special emphasis is made on contrasting these activities in P. aeruginosa and other bacteria. Among unique P. aeruginosa, features are the presence of two distinct autonomously replicating sequences and multiple condensins, which suggests existence of novel regulatory mechanisms.

Haem-based sensors: a still growing old superfamily

The haem-based sensors are chimeric multi-domain proteins responsible for the cellular adaptive responses to environmental changes. The signal transduction is mediated by the sensing capability of the haem-binding domain, which transmits a usable signal to the cognate transmitter domain, responsible for providing the adequate answer. Four major families of haem-based sensors can be recognized, depending on the nature of the haem-binding domain: (i) the haem-binding PAS domain, (ii) the CO-sensitive carbon monoxide oxidation activator, (iii) the haem NO-binding domain, and (iv) the globin-coupled sensors. The functional classification of the haem-binding sensors is based on the activity of the transmitter domain and, traditionally, comprises: (i) sensors with aerotactic function; (ii) sensors with gene-regulating function; and (iii) sensors with unknown function. We have implemented this classification with newly identified proteins, that is, the Streptomyces avermitilis and Frankia sp. that present a C-terminal-truncated globin fused to an N-terminal cofactor-free monooxygenase, the structural-related class of non-haem globins in Bacillus subtilis, Moorella thermoacetica, and Bacillus anthracis, and a haemerythrin-coupled diguanylate cyclase in Vibrio cholerae. This review summarizes the structures, the functions, and the structure-function relationships known to date on this broad protein family. We also propose unresolved questions and new possible research approaches.

The CtrA phosphorelay integrates differentiation and communication in the marine alphaproteobacterium Dinoroseobacter shibae

Background

Dinoroseobacter shibae, a member of the Roseobacter clade abundant in marine environments, maintains morphological heterogeneity throughout growth, with small cells dividing by binary fission and large cells dividing by budding from one or both cell poles. This morphological heterogeneity is lost if the quorum sensing (QS) system is silenced, concurrent with a decreased expression of the CtrA phosphorelay, a regulatory system conserved in Alphaproteobacteria and the master regulator of the Caulobacter crescentus cell cycle. It consists of the sensor histidine kinase CckA, the phosphotransferase ChpT and the transcriptional regulator CtrA. Here we tested if the QS induced differentiation of D. shibae is mediated by the CtrA phosphorelay.

Results

Mutants for ctrA, chpT and cckA showed almost homogeneous cell morphology and divided by binary fission. For ctrA and chpT, expression in trans on a plasmid caused the fraction of cells containing more than two chromosome equivalents to increase above wild-type level, indicating that gene copy number directly controls chromosome number. Transcriptome analysis revealed that CtrA is a master regulator for flagellar biosynthesis and has a great influence on the transition to stationary phase. Interestingly, the expression of the autoinducer synthase genes luxI₂ and luxI₃ was strongly reduced in all three mutants, resulting in loss of biosynthesis of acylated homoserine-lactones with C14 side-chain, but could be restored by expressing these genes in trans. Several phylogenetic clusters of Alphaproteobacteria revealed a CtrA binding site in the promoters of QS genes, including Roseobacters and Rhizobia.

Conclusions

The CtrA phosphorelay induces differentiation of a marine Roseobacter strain that is strikingly different from that of C. crescentus. Instead of a tightly regulated cell cycle and a switch between two morphotypes, the morphology and cell division of Dinoroseobacter shibae are highly heterogeneous. We discovered for the first time that the CtrA phosphorelay controls the biosynthesis of signaling molecules. Thus cell-cell communication and differentiation are interlinked in this organism. This may be a common strategy, since we found a similar genetic set-up in other species in the ecologically relevant group of Alphaproteobacteria. D. shibae will be a valuable model organism to study bacterial differentiation into pleomorphic cells.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-130) contains supplementary material, which is available to authorized users.

Immunization with a DNA adenine methylase over-producing Yersinia pseudotuberculosis vaccine confers robust cross-protection against heterologous pathogenic serotypes

Yersinia pseudotuberculosis is a foodborne pathogen that can cause serious human illness. Although the source and route of transmission often remain obscure, livestock have been implicated in some cases. The diversity of yersiniae present on farms and their widespread distribution in animal and environmental reservoirs necessitates the use of broad prophylactic strategies that are efficacious against many serotypes simultaneously. Herein, immunization of mice with a modified, live attenuated Y. pseudotuberculosis vaccine that overproduces the DNA adenine methylase (Dam(OP)) conferred robust protection against virulent challenge (150-fold LD50) with homologous and heterologous serotypes that have been associated with human disease (O:1, O:1a, O:3). Further, the dam gene was shown to be essential for cell viability in all (7 of 7) Y. pseudotuberculosis strains tested. Direct selection for the inheritance of dam mutant alleles in Y. pseudotuberculosis resulted in dam strain variants that contained compensatory (second-site suppressor) mutations in genes encoding methyl-directed mismatch repair proteins (mutHLS) that are involved in suppression of the non-viable cell phenotype in all (19/19) strains tested. Such dam mutH variants exhibited a significant increase in virulence and spontaneous mutation frequency relative to that of a Dam(OP) vaccine strain. These studies indicate that Y. pseudotuberculosis Dam(OP) strains conferred potent cross-protective efficacy as well as decreased virulence and spontaneous mutation frequency relative to those that lack Dam, which have compensatory mutations in mutHLS loci. These data suggest that development of yersiniae livestock vaccines based on Dam overproduction is a viable mitigation strategy to reduce these potential foodborne contaminants.

Signal-specific temporal response by the Salmonella PhoP/PhoQ regulatory system

The two-component system PhoP/PhoQ controls a large number of genes responsible for a variety of physiological and virulence functions in Salmonella enterica serovar Typhimurium. Here we describe a mechanism whereby the transcriptional activator PhoP elicits expression of dissimilar gene sets when its cognate sensor PhoQ is activated by different signals in the periplasm. We determine that full transcription of over half of the genes directly activated by PhoP requires the Mg(2+) transporter MgtA when the PhoQ inducing signal is low Mg(2+) , but not when PhoQ is activated by mildly acidic pH or the antimicrobial peptide C18G. MgtA promotes the active (i.e. phosphorylated) form of PhoP by removing Mg(2+) from the periplasm, where it functions as a repressing signal for PhoQ. MgtA-dependent expression enhances resistance to the cationic antibiotic polymyxin B. Production of the MgtA protein requires cytoplasmic Mg(2+) levels to drop below a certain threshold, thereby creating a two-tiered temporal response among PhoP-dependent genes.

Induction of olfaction and cancer-related genes in mice fed a high-fat diet as assessed through the mode-of-action by network identification analysis

The pathophysiological mechanisms underlying the development of obesity and metabolic diseases are not well understood. To gain more insight into the genetic mediators associated with the onset and progression of diet-induced obesity and metabolic diseases, we studied the molecular changes in response to a high-fat diet (HFD) by using a mode-of-action by network identification (MNI) analysis. Oligo DNA microarray analysis was performed on visceral and subcutaneous adipose tissues and muscles of male C57BL/6N mice fed a normal diet or HFD for 2, 4, 8, and 12 weeks. Each of these data was queried against the MNI algorithm, and the lists of top 5 highly ranked genes and gene ontology (GO)-annotated pathways that were significantly overrepresented among the 100 highest ranked genes at each time point in the 3 different tissues of mice fed the HFD were considered in the present study. The 40 highest ranked genes identified by MNI analysis at each time point in the different tissues of mice with diet-induced obesity were subjected to clustering based on their temporal patterns. On the basis of the above-mentioned results, we investigated the sequential induction of distinct olfactory receptors and the stimulation of cancer-related genes during the development of obesity in both adipose tissues and muscles. The top 5 genes recognized using the MNI analysis at each time point and gene cluster identified based on their temporal patterns in the peripheral tissues of mice provided novel and often surprising insights into the potential genetic mediators for obesity progression.

Bioimage informatics for understanding spatiotemporal dynamics of cellular processes

The inner environment of the cell is highly dynamic and heterogeneous yet exquisitely organized. Successful completion of cellular processes within this environment depends on the right molecules or molecular complexes to function at the right place at the right time. Understanding spatiotemporal behaviors of cellular processes is therefore essential to understanding their molecular mechanisms at the systems level. These behaviors are usually visualized and recorded using imaging techniques. However, to infer from them systems-level molecular mechanisms, computational analysis and understanding of recorded image data is crucial, not only for acquiring quantitative behavior measurements but also for comprehending complex interactions among the molecules or molecular complexes involved. The technology of computational analysis and understanding of biological images is often referred to simply as bioimage informatics. This article introduces fundamentals of bioimage informatics for understanding spatiotemporal dynamics of cellular processes and reviews recent advances on this topic. Basic bioimage informatics concepts and techniques for characterizing spatiotemporal cell dynamics are introduced first. Studies on specific cellular processes such as cell migration and signal transduction are then used as examples to analyze and summarize recent advances, with the focus on transforming quantitative measurements of spatiotemporal cellular behaviors into knowledge of underlying molecular mechanisms. Despite the advances made, substantial technological challenges remain, especially in representation of spatiotemporal cellular behaviors and inference of systems-level molecular mechanisms. These challenges are briefly discussed. Overall, understanding spatiotemporal cell dynamics will provide critical insights into how specific cellular processes as well as the entire inner cellular environment are dynamically organized and regulated.

Phase resetting reveals network dynamics underlying a bacterial cell cycle

Genomic and proteomic methods yield networks of biological regulatory interactions but do not provide direct insight into how those interactions are organized into functional modules, or how information flows from one module to another. In this work we introduce an approach that provides this complementary information and apply it to the bacterium Caulobacter crescentus, a paradigm for cell-cycle control. Operationally, we use an inducible promoter to express the essential transcriptional regulatory gene ctrA in a periodic, pulsed fashion. This chemical perturbation causes the population of cells to divide synchronously, and we use the resulting advance or delay of the division times of single cells to construct a phase resetting curve. We find that delay is strongly favored over advance. This finding is surprising since it does not follow from the temporal expression profile of CtrA and, in turn, simulations of existing network models. We propose a phenomenological model that suggests that the cell-cycle network comprises two distinct functional modules that oscillate autonomously and couple in a highly asymmetric fashion. These features collectively provide a new mechanism for tight temporal control of the cell cycle in C. crescentus. We discuss how the procedure can serve as the basis for a general approach for probing network dynamics, which we term chemical perturbation spectroscopy (CPS).

During the cell cycle, the cell progresses through a series of stages that are associated with various cell cycle events such as replication of genetic materials. Genetic and molecular dissections have revealed that the cell cycle is regulated by a network of interacting molecules that produces oscillatory dynamics. The major cell cycle regulators have been identified previously in different species and the activity of these regulators oscillates. However, the question of how cell cycle regulators coordinate different cell cycle events during the cell cycle remains controversial. Here, we investigate this question in a model bacterial system for cell cycle, Caulobacter crescentus. We perturb the expression of the master cell cycle regulator ctrA in a pulsatile fashion and quantify the response of the cell cycle to such perturbations. The measured response is contradictory to the existing mechanism of Caulobacter cell cycle control, which views the cell cycle progression as a sequential activation/inhibition process. We propose a new model that involves coupling of multiple oscillators and show the quantitative agreement between this new model and our measurements. We expect this procedure to be generalized and applied to a broad range of systems to obtain information that complements that obtained from other methods.

The Caulobacter crescentus phage phiCbK: genomics of a canonical phage

Background

The bacterium Caulobacter crescentus is a popular model for the study of cell cycle regulation and senescence. The large prolate siphophage phiCbK has been an important tool in C. crescentus biology, and has been studied in its own right as a model for viral morphogenesis. Although a system of some interest, to date little genomic information is available on phiCbK or its relatives.

Results

Five novel phiCbK-like C. crescentus bacteriophages, CcrMagneto, CcrSwift, CcrKarma, CcrRogue and CcrColossus, were isolated from the environment. The genomes of phage phiCbK and these five environmental phage isolates were obtained by 454 pyrosequencing. The phiCbK-like phage genomes range in size from 205 kb encoding 318 proteins (phiCbK) to 280 kb encoding 448 proteins (CcrColossus), and were found to contain nonpermuted terminal redundancies of 10 to 17 kb. A novel method of terminal ligation was developed to map genomic termini, which confirmed termini predicted by coverage analysis. This suggests that sequence coverage discontinuities may be useable as predictors of genomic termini in phage genomes. Genomic modules encoding virion morphogenesis, lysis and DNA replication proteins were identified. The phiCbK-like phages were also found to encode a number of intriguing proteins; all contain a clearly T7-like DNA polymerase, and five of the six encode a possible homolog of the C. crescentus cell cycle regulator GcrA, which may allow the phage to alter the host cell’s replicative state. The structural proteome of phage phiCbK was determined, identifying the portal, major and minor capsid proteins, the tail tape measure and possible tail fiber proteins. All six phage genomes are clearly related; phiCbK, CcrMagneto, CcrSwift, CcrKarma and CcrRogue form a group related at the DNA level, while CcrColossus is more diverged but retains significant similarity at the protein level.

Conclusions

Due to their lack of any apparent relationship to other described phages, this group is proposed as the founding cohort of a new phage type, the phiCbK-like phages. This work will serve as a foundation for future studies on morphogenesis, infection and phage-host interactions in C. crescentus.

Regulatory cohesion of cell cycle and cell differentiation through interlinked phosphorylation and second messenger networks

In Caulobacter crescentus, phosphorylation of key regulators is coordinated with the second messenger cyclic di-GMP to drive cell-cycle progression and differentiation. The diguanylate cyclase PleD directs pole morphogenesis, while the c-di-GMP effector PopA initiates degradation of the replication inhibitor CtrA by the AAA+ protease ClpXP to license S phase entry. Here, we establish a direct link between PleD and PopA reliant on the phosphodiesterase PdeA and the diguanylate cyclase DgcB. PdeA antagonizes DgcB activity until the G1-S transition, when PdeA is degraded by the ClpXP protease. The unopposed DgcB activity, together with PleD activation, upshifts c-di-GMP to drive PopA-dependent CtrA degradation and S phase entry. PdeA degradation requires CpdR, a response regulator that delivers PdeA to the ClpXP protease in a phosphorylation-dependent manner. Thus, CpdR serves as a crucial link between phosphorylation pathways and c-di-GMP metabolism to mediate protein degradation events that irreversibly and coordinately drive bacterial cell-cycle progression and development.

The essential genome of a bacterium

This study reports the essential Caulobacter genome at 8 bp resolution determined by saturated transposon mutagenesis and high-throughput sequencing. This strategy is applicable to full genome essentiality studies in a broad class of bacterial species.

The essential Caulobacter genome was determined at 8 bp resolution using hyper-saturated transposon mutagenesis coupled with high-throughput sequencing.

Essential protein-coding sequences comprise 90% of the essential genome; the remaining 10% comprising essential non-coding RNA sequences, gene regulatory elements and essential genome replication features.

Of the 3876 annotated open reading frames (ORFs), 480 (12.4%) were essential ORFs, 3240 (83.6%) were non-essential ORFs and 156 (4.0%) were ORFs that severely impacted fitness when mutated.

The essential elements are preferentially positioned near the origin and terminus of the Caulobacter chromosome.

This high-resolution strategy is applicable to high-throughput, full genome essentiality studies and large-scale genetic perturbation experiments in a broad class of bacterial species.

The regulatory events that control polar differentiation and cell-cycle progression in the bacterium Caulobacter crescentus are highly integrated, and they have to occur in the proper order (McAdams and Shapiro, 2011). Components of the core regulatory circuit are largely known. Full discovery of its essential genome, including non-coding, regulatory and coding elements, is a prerequisite for understanding the complete regulatory network of this bacterial cell. We have identified all the essential coding and non-coding elements of the Caulobacter chromosome using a hyper-saturated transposon mutagenesis strategy that is scalable and can be readily extended to obtain rapid and accurate identification of the essential genome elements of any sequenced bacterial species at a resolution of a few base pairs.

We engineered a Tn5 derivative transposon (Tn5Pxyl) that carries at one end an inducible outward pointing Pxyl promoter (Christen et al, 2010). We showed that this transposon construct inserts into the genome randomly where it can activate or disrupt transcription at the site of integration, depending on the insertion orientation. DNA from hundred of thousands of transposon insertion sites reading outward into flanking genomic regions was parallel PCR amplified and sequenced by Illumina paired-end sequencing to locate the insertion site in each mutant strain (Figure 1). A single sequencing run on DNA from a mutagenized cell population yielded 118 million raw sequencing reads. Of these, >90 million (>80%) read outward from the transposon element into adjacent genomic DNA regions and the insertion site could be mapped with single nucleotide resolution. This yielded the location and orientation of 428 735 independent transposon insertions in the 4-Mbp Caulobacter genome.

Within non-coding sequences of the Caulobacter genome, we detected 130 non-disruptable DNA segments between 90 and 393 bp long in addition to all essential promoter elements. Among 27 previously identified and validated sRNAs (Landt et al, 2008), three were contained within non-disruptable DNA segments and another three were partially disruptable, that is, insertions caused a notable growth defect. Two additional small RNAs found to be essential are the transfer-messenger RNA (tmRNA) and the ribozyme RNAseP (Landt et al, 2008). In addition to the 8 non-disruptable sRNAs, 29 out of the 130 intergenic essential non-coding sequences contained non-redundant tRNA genes; duplicated tRNA genes were non-essential. We also identified two non-disruptable DNA segments within the chromosomal origin of replication. Thus, we resolved essential non-coding RNAs, tRNAs and essential replication elements within the origin region of the chromosome. An additional 90 non-disruptable small genome elements of currently unknown function were identified. Eighteen of these are conserved in at least one closely related species. Only 2 could encode a protein of over 50 amino acids.

For each of the 3876 annotated open reading frames (ORFs), we analyzed the distribution, orientation, and genetic context of transposon insertions. There are 480 essential ORFs and 3240 non-essential ORFs. In addition, there were 156 ORFs that severely impacted fitness when mutated. The 8-bp resolution allowed a dissection of the essential and non-essential regions of the coding sequences. Sixty ORFs had transposon insertions within a significant portion of their 3′ region but lacked insertions in the essential 5′ coding region, allowing the identification of non-essential protein segments. For example, transposon insertions in the essential cell-cycle regulatory gene divL, a tyrosine kinase, showed that the last 204 C-terminal amino acids did not impact viability, confirming previous reports that the C-terminal ATPase domain of DivL is dispensable for viability (Reisinger et al, 2007; Iniesta et al, 2010). In addition, we found that 30 out of 480 (6.3%) of the essential ORFs appear to be shorter than the annotated ORF, suggesting that these are probably mis-annotated.

Among the 480 ORFs essential for growth on rich media, there were 10 essential transcriptional regulatory proteins, including 5 previously identified cell-cycle regulators (McAdams and Shapiro, 2003; Holtzendorff et al, 2004; Collier and Shapiro, 2007; Gora et al, 2010; Tan et al, 2010) and 5 uncharacterized predicted transcription factors. In addition, two RNA polymerase sigma factors RpoH and RpoD, as well as the anti-sigma factor ChrR, which mitigates rpoE-dependent stress response under physiological growth conditions (Lourenco and Gomes, 2009), were also found to be essential. Thus, a set of 10 transcription factors, 2 RNA polymerase sigma factors and 1 anti-sigma factor are the core essential transcriptional regulators for growth on rich media. To further characterize the core components of the Caulobacter cell-cycle control network, we identified all essential regulatory sequences and operon transcripts. Altogether, the 480 essential protein-coding and 37 essential RNA-coding Caulobacter genes are organized into operons such that 402 individual promoter regions are sufficient to regulate their expression. Of these 402 essential promoters, the transcription start sites (TSSs) of 105 were previously identified (McGrath et al, 2007).

The essential genome features are non-uniformly distributed on the Caulobacter genome and enriched near the origin and the terminus regions. In contrast, the chromosomal positions of the published E. coli essential coding sequences (Rocha, 2004) are preferentially located at either side of the origin (Figure 4A). This indicates that there are selective pressures on chromosomal positioning of some essential elements (Figure 4A).

The strategy described in this report could be readily extended to quickly determine the essential genome for a large class of bacterial species.

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.

RepC protein of the octopine-type Ti plasmid binds to the probable origin of replication within repC and functions only in cis

Vegetative replication and partitioning of many plasmids and some chromosomes of alphaproteobacteria are directed by their repABC operons. RepA and RepB proteins direct the partitioning of replicons to daughter cells, while RepC proteins are replication initiators, although they do not resemble any characterized replication initiation protein. Here we show that the replication origin of an Agrobacterium tumefaciens Ti plasmid resides fully within its repC gene. Purified RepC bound to a site within repC with moderate affinity, high specificity and with twofold cooperativity. The binding site was localized to an AT-rich region that contains a large number of GANTC sites, which have been implicated in replication regulation in related organisms. A fragment of RepC containing residues 26-158 was sufficient to bind DNA, although with limited sequence specificity. This portion of RepC is predicted to have structural homology to members of the MarR family of transcription factors. Overexpression of RepC in A. tumefaciens caused large increases in copy number in cis but did not change the copy number of plasmids containing the same oriV sequence in trans, confirming other observations that RepC functions only in cis.

The architecture and conservation pattern of whole-cell control circuitry

The control circuitry that directs and paces Caulobacter cell cycle progression involves the entire cell operating as an integrated system. This control circuitry monitors the environment and the internal state of the cell, including the cell topology, as it orchestrates orderly activation of cell cycle subsystems and Caulobacter's asymmetric cell division. The proteins of the Caulobacter cell cycle control system and its internal organization are co-conserved across many alphaproteobacteria species, but there are great differences in the regulatory apparatus' functionality and peripheral connectivity to other cellular subsystems from species to species. This pattern is similar to that observed for the "kernels" of the regulatory networks that regulate development of metazoan body plans. The Caulobacter cell cycle control system has been exquisitely optimized as a total system for robust operation in the face of internal stochastic noise and environmental uncertainty. When sufficient details accumulate, as for Caulobacter cell cycle regulation, the system design has been found to be eminently rational and indeed consistent with good design practices for human-designed asynchronous control systems.

Single-gene tuning of Caulobacter cell cycle period and noise, swarming motility, and surface adhesion

We established that the sensor histidine kinase DivJ has an important role in the regulation of C. crescentus cell cycle period and noise. This was accomplished by designing and conducting single-cell experiments to probe the dependence of cell cycle noise on divJ expression and constructing a simplified cell cycle model that captures the dependence of cell cycle noise on DivJ with molecular details.

In addition to its role in regulating the cell cycle, DivJ also affects polar cell development in C. crescentus, regulating swarming motility and surface adhesion. We propose that pleiotropic control of polar cell development by the DivJ–DivK–PleC signaling pathway underlies divJ-dependent tuning of cell swarming and adhesion behaviors.

We have integrated the study of single-cell fluorescence dynamics with a kinetic model simulation to provide direct quantitative evidence that the DivJ histidine kinase is localized to the cell pole through a dynamic diffusion-and-capture mechanism during the C. crescentus cell cycle.

Temporally-coordinated localization of various structural and signaling proteins is critical for proper cell cycle regulation and polar cell development in the bacterium, Caulobacter crescentus. Included among these dynamically-localized regulatory proteins is the sensor histidine kinase, DivJ (Wheeler and Shapiro, 1999). Co-localized with DivJ in the early stalked phase is the phosphorylated response regulator DivK∼P (Jacobs et al, 2001), and the protease ClpXP (McGrath et al, 2006), which degrades the master cell cycle regulator, CtrA (Jenal and Fuchs, 1998). Recent single-cell measurements of surface attached C. crescentus cells have revealed an intriguing role for DivJ in the control of noise in cell division period (Siegal-Gaskins and Crosson, 2008). The noise of the cell cycle increases significantly upon disruption of the divJ gene, with a relatively small accompanying increase in the mean cell cycle time. The deterministic nature of the existing cell cycle models (Li et al, 2008, 2009; Shen et al, 2008) cannot explain the measured increase in cell cycle period and noise in a divJ null strain. Moreover, mechanistic descriptions of how DivJ and its signaling partners are localized and how these proteins underlie the control of polar cell development and cell adhesion in C. crescentus remain immature.

The single-cell experiments and analysis presented herein reveal that C. crescentus cell cycle period and noise can be tuned by DivJ (Figure 2). Specifically, in the case of low (or no) divJ expression the cell cycle is perturbed, and this is quantified by way of the (measured) noise in the cell cycle period. The level of noise is readily controlled through regulated expression of the divJ gene (Figure 2B). A simplified protein interaction network of stalked C. crescentus cell cycle regulation involving minimal components (CtrA, CtrA∼P, DivK, DivK∼P, and DivJ) was constructed to explore such tunability at the molecular level. The agreement of our model with our (and other) experiments suggests this simplified protein regulatory network is sufficient to explain the major features of the C. crescentus cell cycle. Indeed, stochastic simulations of this model using the Gillespie method (Gillespie, 1976) establish the importance of robust DivJ-mediated phosphorylation of its cognate receiver protein, DivK, in regulating the variance of cell cycle oscillations. Increased variability in the concentration of DivK∼P at the single cell level under divJ depletion subsequently leads to increased noise in the regulation of CtrA phosphorylation and degradation. Our experiments and simulations provide evidence that the steady state level of DivK∼P at the single-cell level (as maintained by DivJ) is essential in maintaining regular timing of the cell division period in C. crescentus.

In addition to its role in regulating cell cycle, divJ expression also affects polar cell development in C. crescentus. Specifically, the capacity of swarmer cells to adhere to a glass surface is suppressed at high levels of divJ expression. The effect of elevated divJ expression on the adhesive capacity of the cell is reflected in a reduced rate of two-dimensional biofilm formation. This effect is quantitatively captured by our mathematical model that relates single-cell surface adhesion physiology and biofilm formation dynamics. This result, and our observation that divJ expression tunes swarming motility in semi-solid growth medium, suggests a model in which increased DivJ concentration in the swarmer compartment (due to constitutive overexpression) ultimately results in improper development of polar organelles that are required for adhesion of swarming motility.

Despite the appreciated significance of protein localization for bacterial physiological functions, the molecular mechanism of how polar protein localization is achieved has only been tested in a few cases (Shapiro et al, 2002; Thanbichler and Shapiro, 2008). Mechanisms such as the polar insertion model and diffusion-and-capture have been proposed but the community's knowledge is limited to very few examples (Charles et al, 2001; Rudner et al, 2002). We provide direct evidence from experiments and simulations that the DivJ histidine kinase becomes localized to the cell pole through a dynamic diffusion-and-capture mechanism during the C. crescentus cell cycle (Figure 7). We show that a kinetic model based on a Langmuir adsorption/desorption relationship (Figure 7D) is sufficient to explain the time evolution of the single cell fluorescence time traces (Figure 7C and E) and allows establishing quantitative correspondences between the simulated dynamics and experimentally determined DivJ–EGFP dynamics. This localization mechanism is consistent with a diffusion-and-capture model. In short, the model posits that proteins are randomly distributed and are freely diffusing until they are captured at the site where they ultimately reside (Rudner et al, 2002; Shapiro et al, 2002; Bardy and Maddock, 2007). With a diffusion-and-capture pathway, it has been argued that proteins can be adsorbed either dynamically or statically (Shapiro et al, 2009). Our analysis of DivJ–EGFP in single cells supports a dynamic diffuse-and-capture mechanism for DivJ localization.

Sensor histidine kinases underlie the regulation of a range of physiological processes in bacterial cells, from chemotaxis to cell division. In the gram-negative bacterium Caulobacter crescentus, the membrane-bound histidine kinase, DivJ, is a polar-localized regulator of cell cycle progression and development. We show that DivJ localizes to the cell pole through a dynamic diffusion and capture mechanism rather than by active localization. Analysis of single C. crescentus cells in microfluidic culture demonstrates that controlled expression of divJ permits facile tuning of both the mean and noise of the cell division period. Simulations of the cell cycle that use a simplified protein interaction network capture previously measured oscillatory protein profiles, and recapitulate the experimental observation that deletion of divJ increases the cell cycle period and noise. We further demonstrate that surface adhesion and swarming motility of C. crescentus in semi-solid media can also be tuned by divJ expression. We propose a model in which pleiotropic control of polar cell development by the DivJ–DivK–PleC signaling pathway underlies divJ-dependent tuning of cell swarming and adhesion behaviors.

Spatial aspects in biological system simulations

Mathematical models of the dynamical properties of biological systems aim to improve our understanding of the studied system with the ultimate goal of being able to predict system responses in the absence of experimentation. Despite the enormous advances that have been made in biological modeling and simulation, the inherently multiscale character of biological systems and the stochasticity of biological processes continue to present significant computational and conceptual challenges. Biological systems often consist of well-organized structural hierarchies, which inevitably lead to multiscale problems. This chapter introduces and discusses the advantages and shortcomings of several simulation methods that are being used by the scientific community to investigate the spatiotemporal properties of model biological systems. We first describe the foundations of the methods and then describe their relevance and possible application areas with illustrative examples from our own research. Possible ways to address the encountered computational difficulties are also discussed.

Network biology: a direct approach to study biological function

In this paper we discuss the dualism of gene networks and their role in systems biology. We argue that gene networks (1) can serve as a conceptual framework, forming a fundamental level of a phenomenological description, and (2) are a means to represent and analyze data. The latter point does not only allow a systems analysis but is even amenable for a direct approach to study biological function. Here we focus on the clarity of our main arguments and conceptual meaning of gene networks, rather than the causal inference of gene networks from data. WIREs Syst Biol Med 2011 3 379-391 DOI: 10.1002/wsbm.134 For further resources related to this article, please visit the WIREs website.

Architecture-dependent robustness and bistability in a class of genetic circuits

Understanding the relationship between genotype and phenotype is a challenge in systems biology. An interesting yet related issue is why a particular circuit topology is present in a cell when the same function can supposedly be obtained from an alternative architecture. Here we analyzed two topologically equivalent genetic circuits of coupled positive and negative feedback loops, named NAT and ALT circuits, respectively. The computational search for the oscillation volume of the entire biologically reasonable parameter region through large-scale random samplings shows that the NAT circuit exhibits a distinctly larger fraction of the oscillatory region than the ALT circuit. Such a global robustness difference between two circuits is supplemented by analyzing local robustness, including robustness to parameter perturbations and to molecular noise. In addition, detailed dynamical analysis shows that the molecular noise of both circuits can induce transient switching of the different mechanism between a stable steady state and a stable limit cycle. Our investigation on robustness and dynamics through examples provides insights into the relationship between network architecture and its function.

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.

Origin and analysis of microbial population heterogeneity in bioprocesses

Heterogeneity of industrial production cultures is accepted to a certain degree; however, the underlying mechanisms are seldom perceived or included in the development of new bioprocess control strategies. Population heterogeneity and its basics, perceptible in the diverse proficiency of cells, begins with asymmetric birth and is found to recess during the life cycle. Since inefficient subpopulations have significant impact on the productivity of industrial cultures, cellular heterogeneity needs to be detected and quantified by using high speed detection tools like flow cytometry. Possible origins of population heterogeneity, sophisticated fluorescent techniques for detection of individual cell states, and cutting-edge Omics-technologies for extended information beyond the resolution of fluorescent labelling are highlighted.

Network synchronization landscape reveals compensatory structures, quantization, and the positive effect of negative interactions

Synchronization, in which individual dynamical units keep in pace with each other in a decentralized fashion, depends both on the dynamical units and on the properties of the interaction network. Yet, the role played by the network has resisted comprehensive characterization within the prevailing paradigm that interactions facilitating pairwise synchronization also facilitate collective synchronization. Here we challenge this paradigm and show that networks with best complete synchronization, least coupling cost, and maximum dynamical robustness, have arbitrary complexity but quantized total interaction strength, which constrains the allowed number of connections. It stems from this characterization that negative interactions as well as link removals can be used to systematically improve and optimize synchronization properties in both directed and undirected networks. These results extend the recently discovered compensatory perturbations in metabolic networks to the realm of oscillator networks and demonstrate why "less can be more" in network synchronization.

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.

System-level design of bacterial cell cycle control

Understanding of the cell cycle control logic in Caulobacter has progressed to the point where we now have an integrated view of the operation of an entire bacterial cell cycle system functioning as a state machine. Oscillating levels of a few temporally-controlled master regulator proteins in a cyclical circuit drive cell cycle progression. To a striking degree, the cell cycle regulation is a whole cell phenomenon. Phospho-signaling proteins and proteases dynamically deployed to specific locations on the cell wall are vital. An essential phospho-signaling system integral to the cell cycle circuitry is central to accomplishing asymmetric cell division.

Three (and more) component regulatory systems - auxiliary regulators of bacterial histidine kinases

Two-component signal transduction (TCST) is the most prevalent mechanism employed by microbes to sense and respond to environmental changes. It is characterized by the signal-induced transfer of phosphate from a sensor histidine kinase (HK) to a response regulator (RR), resulting in a cellular response. An emerging theme in the field of TCST signalling is the discovery of auxiliary factors, distinct from the HK and RR, which are capable of influencing phosphotransfer. One group of TCST auxiliary proteins accomplishes this task by acting on HKs. Auxiliary regulators of HKs are widespread and have been identified in all cellular compartments, where they can influence HK activity through interactions with the sensing, transmembrane or enzymatic domains of the HK. The effects of an auxiliary regulator are controlled by its regulated expression, modification and/or through ligand binding. Ultimately, auxiliary regulators can connect a given TCST system to other regulatory networks in the cell or result in regulation of the TCST system in response to an expanded range of stimuli. The studies highlighted in this review draw attention to an emerging view of bacterial TCST systems as core signalling units upon which auxiliary factors act.

Temporal controls of the asymmetric cell division cycle in Caulobacter crescentus

The asymmetric cell division cycle of Caulobacter crescentus is orchestrated by an elaborate gene-protein regulatory network, centered on three major control proteins, DnaA, GcrA and CtrA. The regulatory network is cast into a quantitative computational model to investigate in a systematic fashion how these three proteins control the relevant genetic, biochemical and physiological properties of proliferating bacteria. Different controls for both swarmer and stalked cell cycles are represented in the mathematical scheme. The model is validated against observed phenotypes of wild-type cells and relevant mutants, and it predicts the phenotypes of novel mutants and of known mutants under novel experimental conditions. Because the cell cycle control proteins of Caulobacter are conserved across many species of alpha-proteobacteria, the model we are proposing here may be applicable to other genera of importance to agriculture and medicine (e.g., Rhizobium, Brucella).