Chromosomal Arrangement of Phosphorelay Genes Couples Sporulation and DNA Replication

Biophysical modeling reveals the transcriptional regulatory mechanism of Spo0A, the master regulator in starving Bacillus subtilis

In starving Bacillus subtilis bacteria, the initiation of two survival programs-biofilm formation and sporulation-is controlled by the same phosphorylated master regulator, Spo0A~P. Its gene, spo0A, is transcribed from two promoters, Pv and Ps, that are, respectively, regulated by RNA polymerase (RNAP) holoenzymes bearing σA and σH. Notably, transcription is directly autoregulated by Spo0A~P binding sites known as 0A1, 0A2, and 0A3 box, located in between the two promoters. It remains unclear whether, at the onset of starvation, these boxes activate or repress spo0A expression, and whether the Spo0A~P transcriptional feedback plays a role in the increase in spo0A expression. Based on the experimental data of the promoter activities under systematic perturbation of the promoter architecture, we developed a biophysical model of transcriptional regulation of spo0A by Spo0A~P binding to each of the 0A boxes. The model predicts that Spo0A~P binding to its boxes does not affect the RNAP recruitment to the promoters but instead affects the transcriptional initiation rate. Moreover, the effects of Spo0A~P binding to 0A boxes are mainly repressive and saturated early at the onset of starvation. Therefore, the increase in spo0A expression is mainly driven by the increase in RNAP holoenzyme levels. Additionally, we reveal that Spo0A~P affinity to 0A boxes is strongest at 0A3 and weakest at 0A2 and that there are attractive forces between the occupied 0A boxes. Our findings, in addition to clarifying how the sporulation master regulator is controlled, offer a framework to predict regulatory outcomes of complex gene-regulatory mechanisms.

IMPORTANCE

Cell differentiation is often critical for survival. In bacteria, differentiation decisions are controlled by transcriptional master regulators under transcriptional feedback control. Therefore, understanding how master regulators are transcriptionally regulated is required to understand differentiation. However, in many cases, the underlying regulation is complex, with multiple transcription factor binding sites and multiple promoters, making it challenging to dissect the exact mechanisms. Here, we address this problem for the Bacillus subtilis master regulator Spo0A. Using a biophysical model, we quantitatively characterize the effect of individual transcription factor binding sites on each spo0A promoter. Furthermore, the model allows us to identify the specific transcription step that is affected by transcription factor binding. Such a model is promising for the quantitative study of a wide range of master regulators involved in transcriptional feedback.

Most bacterial gene families are biased toward specific chromosomal positions

The arrangement of genes along bacterial chromosomes influences their expression through growth rate-dependent gene copy number changes during DNA replication. Although translation- and transcription-related genes often cluster near the origin of replication, the extent of positional biases across gene families remains unclear. We hypothesized that natural selection broadly favors specific chromosomal positions to optimize growth rate-dependent expression. Analyzing 910 bacterial species and proteomics data from Escherichia coli and Bacillus subtilis, we found that about two-thirds of bacterial gene families are positionally biased. Natural selection drives genes mainly toward the origin or terminus of replication, with the strongest selection in fast-growing species. Our findings reveal chromosomal positioning as a fundamental mechanism for coordinating gene expression with growth rate, highlighting evolutionary constraints on bacterial genome architecture.

Autoregulation of the Master Regulator Spo0A Controls Cell-Fate Decisions in Bacillus subtilis

Spo0A in Bacillus subtilis is activated by phosphorylation (Spo0A~P) upon starvation and differentially controls a set of genes involved in biofilm formation and sporulation. The spo0A gene is transcribed by two distinct promoters, a σA-recognized upstream promoter Pv during growth, and a σH-recognized downstream promoter Ps during starvation, and appears to be autoregulated by four Spo0A~P binding sites (0A1-4 boxes) localized between two promoters. However, the autoregulatory mechanisms and their impact on differentiation remain elusive. Here, we determined the relative affinity of Spo0A~P for each 0A box and dissected each promoter in combination with the systematic 0A box mutations. The data revealed that (1) the Pv and Ps promoters are on and off, respectively, under nutrient-rich conditions without Spo0A~P, (2) the Ps promoter is activated by first 0A3 and then 0A1 during early starvation with low Spo0A~P, (3) during later starvation with high Spo0A~P, the Pv promoter is repressed by first 0A1 and then 0A2 and 0A4, and (4) during prolonged starvation, both promoters are silenced by all 0A boxes with very high Spo0A~P. Our results indicate that the autoregulation of spo0A is one of the key determinants to achieve a developmental increase in Spo0A~P, leading to a temporal window for entry into biofilm formation or sporulation.

Regulation of biofilm gene expression by DNA replication in Bacillus subtilis

Bacillus subtilis relies on biofilms for survival in harsh environments. Extracellular polymeric substance (EPS) is a crucial component of biofilms, yet the dynamics of EPS production in single cells remain elusive. To unveil the modulation of EPS synthesis, we built a minimal network model comprising the SinI-SinR-SlrR module, Spo0A, and EPS. Stochastic simulations revealed that antagonistic interplay between SinI and SinR enables EPS production in bursts. SlrR widens these bursts and increases their frequency by stabilizing SinR-SlrR complexes and depleting free SinR. DNA replication and chromosomal positioning of key genes dictate pulsatile changes in the slrR:sinR gene dosage ratio (gr) and Spo0A-P levels, each promoting EPS production in distinct phases of the cell cycle. As the cell cycle lengthens with nutrient stress, the duty cycle of gr pulsing decreases, whereas the amplitude of Spo0A-P pulses elevates. This coordinated response facilitates keeping a constant proportion of EPS-secreting cells within colonies across diverse nutrient conditions. Our results suggest that bacteria may 'encode' eps expression through strategic chromosomal organization. This work illuminates how stochastic protein interactions, gene copy number imbalance, and cell-cycle dynamics orchestrate EPS synthesis, offering a deeper understanding of biofilm formation.

Artificial activation of both σH and Spo0A in Bacillus subtilis enforced initiation of spore development at the vegetatively growing phase

When Bacillus subtilis cells face environmental deterioration, such as exhaustion of nutrients and an increase in cell density, they form spores. It is known that phosphorylation of Spo0A and activation of σH are key events at the initiation of sporulation. However, the initiation of sporulation is an extremely complicated process, and the relationship between these two events remains to be elucidated. To determine the minimum requirements for triggering sporulation initiation, we attempted to induce cell sporulation at the log phase, regardless of nutrients and cell density. In rich media such as Luria-Bertani (LB) medium, the cells of B. subtilis do not sporulate efficiently, possibly because of excess nutrition. When the amount of xylose in the LB medium was limited, σH -dependent transcription of the strain, in which sigA was under the control of the xylose-inducible promoter, was induced, and the frequency of sporulation was elevated according to the decreased level of σA. We also employed a fusion of sad67, which codes for an active form of Spo0A, and the IPTG-inducible promoter. The combination of lowered σA expression and activated Spo0A allowed the cells in the log phase to stop growing and rush into spore development. This observation of enforced initiation of sporulation in the mutant strain was detected even in the presence of the wild-type strain, suggesting that only intracellular events initiate and fulfill spore development regardless of extracellular conditions. Under natural sporulation conditions, the amount of σA did not change drastically throughout growth. Mechanisms that sequester σA from the core RNA polymerase and help σH to become active exist, but this has not yet been elucidated.

Transcription-replication interactions reveal bacterial genome regulation

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

The motility-matrix production switch in Bacillus subtilis-a modeling perspective

Phenotype switching can be triggered by external stimuli and by intrinsic stochasticity. Here, we focus on the motility-matrix production switch in Bacillus subtilis. We use modeling to describe the SinR-SlrR bistable switch and its regulation by SinI and to distinguish different sources of stochasticity. Our simulations indicate that intrinsic fluctuations in the synthesis of SinI are insufficient to drive spontaneous switching and suggest that switching is triggered by upstream noise from the Spo0A phosphorelay. IMPORTANCE The switch from motility to matrix production is the first step toward biofilm formation and, thus, to multicellular behavior in Bacillus subtilis. The transition is governed by a bistable switch based on the interplay of the regulators SinR and SlrR, while SinI transmits upstream signals to that switch. Quantitative modeling can be used to study the switching dynamics. Here, we build such a model step by step to describe the dynamics of the switch and its regulation and to study how spontaneous switching is triggered by upstream noise from the Spo0A phosphorelay.

Ploidy in Vibrio natriegens: Very Dynamic and Rapidly Changing Copy Numbers of Both Chromosomes

Vibrio natriegens is the fastest-growing bacterium, with a doubling time of approximately 12-14 min. It has a high potential for basic research and biotechnological applications, e.g., it can be used for the cell-free production of (labeled) heterologous proteins, for synthetic biological applications, and for the production of various compounds. However, the ploidy level in V. natriegens remains unknown. At nine time points throughout the growth curve, we analyzed the numbers of origins and termini of both chromosomes with qPCR and the relative abundances of all genomic sites with marker frequency analyses. During the lag phase until early exponential growth, the origin copy number and origin/terminus ratio of chromosome 1 increased severalfold, but the increase was lower for chromosome 2. This increase was paralleled by an increase in cell volume. During the exponential phase, the origin/terminus ratio and cell volume decreased again. This highly dynamic and fast regulation has not yet been described for any other species. In this study, the gene dosage increase in origin-adjacent genes during the lag phase is discussed together with the nonrandom distribution of genes on the chromosomes of V. natriegens. Taken together, the results of this study provide the first comprehensive overview of the chromosome dynamics in V. natriegens and will guide the optimization of molecular biological characterization and biotechnological applications.

Transcription-replication interactions reveal principles of bacterial genome regulation

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. These modes interact with a changing cellular environment to yield highly dynamic expression patterns2. In bacteria, the relationship between a gene's regulatory architecture and its expression is well understood for individual model gene circuits3,4. However, a broader perspective of these dynamics at the genome-scale is lacking, in part because bacterial transcriptomics have hitherto captured only a static snapshot of expression averaged across millions of cells5. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on each gene's transcriptional response to its own replication, which we term the Transcription-Replication Interaction Profile (TRIP). We found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal a gene's local regulatory context. While the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, including altered timing or amplitude of expression, and this is shaped by factors such as intra-operon position, repression state, or presence on mobile genetic elements. Our transcriptome analysis also simultaneously captures global properties, such as the rates of replication and transcription, as well as the nestedness of replication patterns. This work challenges previous notions of the drivers of expression heterogeneity within a population of cells, and unearths a previously unseen world of gene transcription dynamics.

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

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

The Slowdown of Growth Rate Controls the Single-Cell Distribution of Biofilm Matrix Production via an SinI-SinR-SlrR Network

In Bacillus subtilis, master regulator Spo0A controls several cell-differentiation pathways. Under moderate starvation, phosphorylated Spo0A (Spo0A~P) induces biofilm formation by indirectly activating genes controlling matrix production in a subpopulation of cells via an SinI-SinR-SlrR network. Under severe starvation, Spo0A~P induces sporulation by directly and indirectly regulating sporulation gene expression. However, what determines the heterogeneity of individual cell fates is not fully understood. In particular, it is still unclear why, despite being controlled by a single master regulator, biofilm matrix production and sporulation seem mutually exclusive on a single-cell level. In this work, with mathematical modeling, we showed that the fluctuations in the growth rate and the intrinsic noise amplified by the bistability in the SinI-SinR-SlrR network could explain the single-cell distribution of matrix production. Moreover, we predicted an incoherent feed-forward loop; the decrease in the cellular growth rate first activates matrix production by increasing in Spo0A phosphorylation level but then represses it via changing the relative concentrations of SinR and SlrR. Experimental data provide evidence to support model predictions. In particular, we demonstrate how the degree to which matrix production and sporulation appear mutually exclusive is affected by genetic perturbations. IMPORTANCE The mechanisms of cell-fate decisions are fundamental to our understanding of multicellular organisms and bacterial communities. However, even for the best-studied model systems we still lack a complete picture of how phenotypic heterogeneity of genetically identical cells is controlled. Here, using B. subtilis as a model system, we employ a combination of mathematical modeling and experiments to explain the population-level dynamics and single-cell level heterogeneity of matrix gene expression. The results demonstrate how the two cell fates, biofilm matrix production and sporulation, can appear mutually exclusive without explicitly inhibiting one another. Such a mechanism could be used in a wide range of other biological systems.

Transcriptome Dynamics of Pseudomonas aeruginosa during Transition from Overlapping To Non-Overlapping Cell Cycles

Bacteria either duplicate their chromosome once per cell division or a new round of replication is initiated before the cells divide, thus cell cycles overlap. Here, we show that the opportunistic pathogen Pseudomonas aeruginosa switches from fast growth with overlapping cell cycles to sustained slow growth with only one replication round per cell division when cultivated under standard laboratory conditions. The transition was characterized by fast-paced, sequential changes in transcriptional activity along the ori-ter axis of the chromosome reflecting adaptation to the metabolic needs during both growth phases. Quorum sensing (QS) activity was highest at the onset of the slow growth phase with non-overlapping cell cycles. RNA sequencing of subpopulations of these cultures sorted based on their DNA content, revealed a strong gene dosage effect as well as specific expression patterns for replicating and nonreplicating cells. Expression of flagella and mexE, involved in multidrug efflux was restricted to cells that did not replicate, while those that did showed a high activity of the cell division locus and recombination genes. A possible role of QS in the formation of these subpopulations upon switching to non-overlapping cell cycles could be a subject of further research. IMPORTANCE The coordination of gene expression with the cell cycle has so far been studied only in a few bacteria, the bottleneck being the need for synchronized cultures. Here, we determined replication-associated effects on transcription by comparing Pseudomonas aeruginosa cultures that differ in their growth mode and number of replicating chromosomes. We further show that cell cycle-specific gene regulation can be principally identified by RNA sequencing of subpopulations from cultures that replicate only once per cell division and that are sorted according to their DNA content. Our approach opens the possibility to study asynchronously growing bacteria from a wide phylogenetic range and thereby enhance our understanding of the evolution of cell cycle control on the transcriptional level.

Patterns of abundance, chromosomal localization, and domain organization among c-di-GMP-metabolizing genes revealed by comparative genomics of five alphaproteobacterial orders

BACKGROUND

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is a bacterial second messenger that affects diverse processes in different bacteria, including the cell cycle, motility, and biofilm formation. Its cellular levels are controlled by the opposing activities of two types of enzymes, with synthesis by diguanylate cyclases containing a GGDEF domain and degradation by phosphodiesterases containing either an HD-GYP or an EAL domain. These enzymes are ubiquitous in bacteria with up to 50 encoded in some genomes, the specific functions of which are mostly unknown.

RESULTS

We used comparative analyses to identify genomic patterns among genes encoding proteins with GGDEF, EAL, and HD-GYP domains in five orders of the class Alphaproteobacteria. GGDEF-containing sequences and GGDEF-EAL hybrids were the most abundant and had the highest diversity of co-occurring auxiliary domains while EAL and HD-GYP containing sequences were less abundant and less diverse with respect to auxiliary domains. There were striking patterns in the chromosomal localizations of the genes found in two of the orders. The Rhodobacterales' EAL-encoding genes and Rhizobiales' GGDEF-EAL-encoding genes showed opposing patterns of distribution compared to the GGDEF-encoding genes. In the Rhodobacterales, the GGDEF-encoding genes showed a tri-modal distribution with peaks mid-way between the origin (ori) and terminus (ter) of replication and at ter while the EAL-encoding genes peaked near ori. The patterns were more complex in the Rhizobiales, but the GGDEF-encoding genes were biased for localization near ter.

CONCLUSIONS

The observed patterns in the chromosomal localizations of these genes suggest a coupling of synthesis and hydrolysis of c-di-GMP with the cell cycle. Moreover, the higher proportions and diversities of auxiliary domains associated with GGDEF domains and GGDEF-EAL hybrids compared to EAL or HD-GYP domains could indicate that more stimuli affect synthesis compared to hydrolysis of c-di-GMP.

Real-time detection of response regulator phosphorylation dynamics in live bacteria

Bacteria utilize two-component system (TCS) signal transduction pathways to sense and adapt to changing environments. In a typical TCS, a stimulus induces a sensor histidine kinase (SHK) to phosphorylate a response regulator (RR), which then dimerizes and activates a transcriptional response. Here, we demonstrate that oligomerization-dependent depolarization of excitation light by fused mNeonGreen fluorescent protein probes enables real-time monitoring of RR dimerization dynamics in live bacteria. Using inducible promoters to independently express SHKs and RRs, we detect RR dimerization within seconds of stimulus addition in several model pathways. We go on to combine experiments with mathematical modeling to reveal that TCS phosphosignaling accelerates with SHK expression but decelerates with RR expression and SHK phosphatase activity. We further observe pulsatile activation of the SHK NarX in response to addition and depletion of the extracellular electron acceptor nitrate when the corresponding TCS is expressed from both inducible systems and the native chromosomal operon. Finally, we combine our method with polarized light microscopy to enable single-cell measurements of RR dimerization under changing stimulus conditions. Direct in vivo characterization of RR oligomerization dynamics should enable insights into the regulation of bacterial physiology.

Bacillus subtilis Histidine Kinase KinC Activates Biofilm Formation by Controlling Heterogeneity of Single-Cell Responses

In Bacillus subtilis, biofilm and sporulation pathways are both controlled by a master regulator, Spo0A, which is activated by phosphorylation via a phosphorelay-a cascade of phosphotransfer reactions commencing with autophosphorylation of histidine kinases KinA, KinB, KinC, KinD, and KinE. However, it is unclear how the kinases, despite acting via the same regulator, Spo0A, differentially regulate downstream pathways, i.e., how KinA mainly activates sporulation genes and KinC mainly activates biofilm genes. In this work, we found that KinC also downregulates sporulation genes, suggesting that KinC has a negative effect on Spo0A activity. To explain this effect, with a mathematical model of the phosphorelay, we revealed that unlike KinA, which always activates Spo0A, KinC has distinct effects on Spo0A at different growth stages: during fast growth, KinC acts as a phosphate source and activates Spo0A, whereas during slow growth, KinC becomes a phosphate sink and contributes to decreasing Spo0A activity. However, under these conditions, KinC can still increase the population-mean biofilm matrix production activity. In a population, individual cells grow at different rates, and KinC would increase the Spo0A activity in the fast-growing cells but reduce the Spo0A activity in the slow-growing cells. This mechanism reduces single-cell heterogeneity of Spo0A activity, thereby increasing the fraction of cells that activate biofilm matrix production. Thus, KinC activates biofilm formation by controlling the fraction of cells activating biofilm gene expression. IMPORTANCE In many bacterial and eukaryotic systems, multiple cell fate decisions are activated by a single master regulator. Typically, the activities of the regulators are controlled posttranslationally in response to different environmental stimuli. The mechanisms underlying the ability of these regulators to control multiple outcomes are not understood in many systems. By investigating the regulation of Bacillus subtilis master regulator Spo0A, we show that sensor kinases can use a novel mechanism to control cell fate decisions. By acting as a phosphate source or sink, kinases can interact with one another and provide accurate regulation of the phosphorylation level. Moreover, this mechanism affects the cell-to-cell heterogeneity of the transcription factor activity and eventually determines the fraction of different cell types in the population. These results demonstrate the importance of intercellular heterogeneity for understanding the effects of genetic perturbations on cell fate decisions. Such effects can be applicable to a wide range of cellular systems.

Measurement of circular intensity differential scattering (CIDS) from single airborne aerosol particles for bioaerosol detection and identification

The circular intensity differential scattering (CIDS), i.e. the normalized Mueller matrix element -S₁₄/S₁₁, can be used to detect the helical structures of DNA molecules in biological systems, however, no CIDS measurement from single particles has been reported to date. We report an innovative method for measuring CIDS phase functions from single particles individually flowing through a scattering laser beam. CIDS signals were obtained from polystyrene latex (PSL) microspheres with or without coating of DNA molecules, tryptophan particles, and aggregates of B. subtilis spores, at the size of 3 μm in diameter. Preliminary results show that this method is able to measure CIDS phase function in tens of microseconds from single particles, and has the ability to identify particles containing biological molecules.

Bacteriophage self-counting in the presence of viral replication

When host cells are in low abundance, temperate bacteriophages opt for dormant (lysogenic) infection. Phage lambda implements this strategy by increasing the frequency of lysogeny at higher multiplicity of infection (MOI). However, it remains unclear how the phage reliably counts infecting viral genomes even as their intracellular number increases because of replication. By combining theoretical modeling with single-cell measurements of viral copy number and gene expression, we find that instead of hindering lambda's decision, replication facilitates it. In a nonreplicating mutant, viral gene expression simply scales with MOI rather than diverging into lytic (virulent) and lysogenic trajectories. A similar pattern is followed during early infection by wild-type phage. However, later in the infection, the modulation of viral replication by the decision genes amplifies the initially modest gene expression differences into divergent trajectories. Replication thus ensures the optimal decision-lysis upon single-phage infection and lysogeny at higher MOI.

Maximization of information transmission influences selection of native phosphorelay architectures

Phosphorelays are signal transduction circuits that sense environmental changes and adjust cellular metabolism. Five different circuit architectures account for 99% of all phosphorelay operons annotated in over 9,000 fully sequenced genomes. Here we asked what biological design principles, if any, could explain selection among those architectures in nature. We began by studying kinetically well characterized phosphorelays (Spo0 of Bacillus subtilis and Sln1 of Saccharomyces cerevisiae). We find that natural circuit architecture maximizes information transmission in both cases. We use mathematical models to compare information transmission among the architectures for a realistic range of concentration and parameter values. Mapping experimentally determined phosphorelay protein concentrations onto that range reveals that the native architecture maximizes information transmission in sixteen out of seventeen analyzed phosphorelays. These results suggest that maximization of information transmission is important in the selection of native phosphorelay architectures, parameter values and protein concentrations.

Connection Between Chromosomal Location and Function of CtrA Phosphorelay Genes in Alphaproteobacteria

Most bacterial chromosomes are circular, with replication starting at one origin (ori) and proceeding on both replichores toward the terminus (ter). Several studies have shown that the location of genes relative to ori and ter can have profound effects on regulatory networks and physiological processes. The CtrA phosphorelay is a gene regulatory system conserved in most alphaproteobacteria. It was first discovered in Caulobacter crescentus where it controls replication and division into a stalked and a motile cell in coordination with other factors. The locations of the ctrA gene and targets of this response regulator on the chromosome affect their expression through replication-induced DNA hemi-methylation and specific positioning along a CtrA activity gradient in the dividing cell, respectively. Here we asked to what extent the location of CtrA regulatory network genes might be conserved in the alphaproteobacteria. We determined the locations of the CtrA phosphorelay and associated genes in closed genomes with unambiguously identifiable ori from members of five alphaproteobacterial orders. The location of the phosphorelay genes was the least conserved in the Rhodospirillales followed by the Sphingomonadales. In the Rhizobiales a trend toward certain chromosomal positions could be observed. Compared to the other orders, the CtrA phosphorelay genes were conserved closer to ori in the Caulobacterales. In contrast, the genes were highly conserved closer to ter in the Rhodobacterales. Our data suggest selection pressure results in differential positioning of CtrA phosphorelay and associated genes in alphaproteobacteria, particularly in the orders Rhodobacterales, Caulobacterales and Rhizobiales that is worth deeper investigation.

RNA localization in prokaryotes: Where, when, how, and why

Only recently has it been recognized that the transcriptome of bacteria and archaea can be spatiotemporally regulated. All types of prokaryotic transcripts-rRNAs, tRNAs, mRNAs, and regulatory RNAs-may acquire specific localization and these patterns can be temporally regulated. In some cases bacterial RNAs reside in the vicinity of the transcription site, but in many others, transcripts show distinct localizations to the cytoplasm, the inner membrane, or the pole of rod-shaped species. This localization, which often overlaps with that of the encoded proteins, can be achieved either in a translation-dependent or translation-independent fashion. The latter implies that RNAs carry sequence-level features that determine their final localization with the aid of RNA-targeting factors. Localization of transcripts regulates their posttranscriptional fate by affecting their degradation and processing, translation efficiency, sRNA-mediated regulation, and/or propensity to undergo RNA modifications. By facilitating complex assembly and liquid-liquid phase separation, RNA localization is not only a consequence but also a driver of subcellular spatiotemporal complexity. We foresee that in the coming years the study of RNA localization in prokaryotes will produce important novel insights regarding the fundamental understanding of membrane-less subcellular organization and lead to practical outputs with biotechnological and therapeutic implications. This article is categorized under: RNA Export and Localization > RNA Localization Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Abhyankar W, N. Swarge B, van der Wel N, Kramer G, Brul S, J. de Koning L. Artificial Sporulation Induction (ASI) by kinA Overexpression Affects the Proteomes and Properties of Bacillus subtilis Spores

To facilitate more accurate spore proteomic analysis, the current study focuses on inducing homogeneous sporulation by overexpressing kinA and assesses the effect of synchronized sporulation initiation on spore resistance, structures, the germination behavior at single-spore level and the proteome. The results indicate that, in our set up, the sporulation by overexpressing kinA can generate a spore yield of 70% within 8 h. The procedure increases spore wet heat resistance and thickness of the spore coat and cortex layers, whilst delaying the time to spore phase-darkening and burst after addition of germinant. The proteome analysis reveals that the upregulated proteins in the kinA induced spores, compared to spores without kinA induction, as well as the 'wildtype' spores, are mostly involved in spore formation. The downregulated proteins mostly belong to the categories of coping with stress, carbon and nitrogen metabolism, as well as the regulation of sporulation. Thus, while kinA overexpression enhances synchronicity in sporulation initiation, it also has profound effects on the central equilibrium of spore formation and spore germination, through modulation of the spore molecular composition and stress resistance physiology.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Stochastic pulsing of gene expression enables the generation of spatial patterns in Bacillus subtilis biofilms

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells, but it is unclear whether it has a role in the development of multicellular systems. Here, we show how stochastic pulsing of gene expression enables spatial patterns to form in a model multicellular system, Bacillus subtilis bacterial biofilms. We use quantitative microscopy and time-lapse imaging to observe pulses in the activity of the general stress response sigma factor σ^B in individual cells during biofilm development. Both σ^B and sporulation activity increase in a gradient, peaking at the top of the biofilm, even though σ^B represses sporulation. As predicted by a simple mathematical model, increasing σ^B expression shifts the peak of sporulation to the middle of the biofilm. Our results demonstrate how stochastic pulsing of gene expression can play a key role in pattern formation during biofilm development.

Stochastic pulsing of gene expression can generate phenotypic diversity in a genetically identical population of cells. Here, the authors show that stochastic pulsing in the expression of a sigma factor enables the formation of spatial patterns in a multicellular system, Bacillus subtilis bacterial biofilms.

Models of protein production along the cell cycle: An investigation of possible sources of noise

In this article, we quantitatively study, through stochastic models, the effects of several intracellular phenomena, such as cell volume growth, cell division, gene replication as well as fluctuations of available RNA polymerases and ribosomes. These phenomena are indeed rarely considered in classic models of protein production and no relative quantitative comparison among them has been performed. The parameters for a large and representative class of proteins are determined using experimental measures. The main important and surprising conclusion of our study is to show that despite the significant fluctuations of free RNA polymerases and free ribosomes, they bring little variability to protein production contrary to what has been previously proposed in the literature. After verifying the robustness of this quite counter-intuitive result, we discuss its possible origin from a theoretical view, and interpret it as the result of a mean-field effect.

Overcoming the Cost of Positive Autoregulation by Accelerating the Response with a Coupled Negative Feedback

A fundamental trade-off between rapid response and optimal expression of genes below cytotoxic levels exists for many signaling circuits, particularly for positively autoregulated systems with an inherent response delay. Here, we describe a regulatory scheme in the E. coli PhoB-PhoR two-component system, which overcomes the cost of positive feedback and achieves both fast and optimal steadystate response for maximal fitness across different environments. Quantitation of the cellular activities enables accurate modeling of the response dynamics to describe how requirements for optimal protein concentrations place limits on response speed. An observed fast response that exceeds the limit led to the prediction and discovery of a coupled negative autoregulation, which allows fast gene expression without increasing steady-state levels. We demonstrate the fitness advantages for the coupled feedbacks in both dynamic and stable environments. Such regulatory schemes offer great flexibility for accurate control of gene expression levels and dynamics upon environmental changes.

Positive autoregulation of transcription produces a delayed response. Gao and Stock describe the limit of response delay caused by requirements of optimal protein levels in the PhoBR twocomponent system. Coupled negative autoregulation is discovered to allow a strong promoter for fast response without incurring cost of increasing protein expression levels.

Enhancement of Bacillus subtilis Growth and Sporulation by Two-Stage Solid-State Fermentation Strategy

Two-stage solid-state fermentation strategy was exploited and systematically optimized to enhance Bacillus subtilis growth and sporulation for increasing effective cell number in B. subtilis microbial ecological agents. The first stage focused on improving cell growth followed by the second stage aiming to enhance both cell growth and sporulation. The optimal fermentation condition was that temperature changed from 37 °C to 47 °C at a fermentation time of 48 h and Mn2+ content in medium was 4.9 mg MnSO4/g dry medium. Solid medium properties were improved by the optimal two-stage fermentation. HPLC results demonstrated that glucose utilization was facilitated and low-field nuclear magnetic resonance (LF-NMR) results showed that more active sites in medium for microbial cells were generated during the optimal two-stage fermentation. Moreover, microbial growth and sporulation were enhanced simultaneously during the second stage of fermentation through delaying microbial decline phase and increasing sporulation rate. As a result, effective cell number of B. subtilis reached 1.79 × 1010/g dry medium after fermentation for 72 h, which was 29.7% and 8.48% higher than that of conventional fermentation for 72 h and 48 h, respectively. Therefore, the optimal two-stage fermentation could increase the effective cell number of B. subtilis microbial ecological agents efficiently.

An Engineered B. subtilis Inducible Promoter System with over 10 000-Fold Dynamic Range

Bacillus subtilis is the leading model Gram-positive bacterium, and a widely used chassis for industrial protein production. However, B. subtilis research is limited by a lack of inducible promoter systems with low leakiness and high dynamic range. Here, we engineer an inducible promoter system based on the T7 RNA Polymerase (T7 RNAP), the lactose repressor LacI, and the chimeric promoter P_T7lac, integrated as a single copy in the B. subtilis genome. In the absence of IPTG, LacI strongly represses T7 RNAP and P_T7lac and minimizes leakiness. Addition of IPTG derepresses P_T7lac and simultaneously induces expression of T7RNAP, which results in very high output expression. Using green fluorescent and β-galactosidase reporter proteins, we estimate that this LacI-T7 system can regulate expression with a dynamic range of over 10 000, by far the largest reported for an inducible B. subtilis promoter system. Furthermore, LacI-T7 responds to similar IPTG concentrations and with similar kinetics as the widely used P_hy-spank IPTG-inducible system, which we show has a dynamic range of at most 300 in a similar genetic context. Due to its superior performance, our LacI-T7 system should have broad applications in fundamental B. subtilis biology studies and biotechnology.

Optogenetic control of Bacillus subtilis gene expression

The Gram-positive bacterium Bacillus subtilis exhibits complex spatial and temporal gene expression signals. Although optogenetic tools are ideal for studying such processes, none has been engineered for this organism. Here, we port a cyanobacterial light sensor pathway comprising the green/red photoreversible two-component system CcaSR, two metabolic enzymes for production of the chromophore phycocyanobilin (PCB), and an output promoter to control transcription of a gene of interest into B. subtilis. Following an initial non-functional design, we optimize expression of pathway genes, enhance PCB production via a translational fusion of the biosynthetic enzymes, engineer a strong chimeric output promoter, and increase dynamic range with a miniaturized photosensor kinase. Our final design exhibits over 70-fold activation and rapid response dynamics, making it well-suited to studying a wide range of gene regulatory processes. In addition, the synthetic biology methods we develop to port this pathway should make B. subtilis easier to engineer in the future.

Electrical Polarization Enables Integrative Quality Control during Bacterial Differentiation into Spores

Quality control of offspring is important for the survival of cells. However, the mechanisms by which quality of offspring cells may be checked while running genetic programs of cellular differentiation remain unclear. Here we investigated quality control during sporulating in Bacillus subtilis by combining single-cell time-lapse microscopy, molecular biology, and mathematical modeling. Our results revealed that the quality control via premature germination is coupled with the electrical polarization of outer membranes of developing forespores. The forespores that accumulate fewer cations on their surface are more likely to be aborted. This charge accumulation enables the projection of multi-dimensional information about the external environment and morphological development of the forespore into one-dimensional information of cation accumulation. We thus present a paradigm of cellular regulation by bacterial electrical signaling. Moreover, based on the insight we gain, we propose an electrophysiology-based approach of reducing the yield and quality of Bacillus endospores.

•

Quality control during bacterial sporulation is coupled with cation accumulation

•

Cation accumulation prevents premature germination

•

Cation accumulation integrates information on morphological defects and environments

•

Spores are less fit when sporulated with Thioflavin T

Are There Circadian Clocks in Non-Photosynthetic Bacteria?

Circadian clocks in plants, animals, fungi, and in photosynthetic bacteria have been well-described. Observations of circadian rhythms in non-photosynthetic Eubacteria have been sporadic, and the molecular basis for these potential rhythms remains unclear. Here, we present the published experimental and bioinformatical evidence for circadian rhythms in these non-photosynthetic Eubacteria. From this, we suggest that the timekeeping functions of these organisms will be best observed and studied in their appropriate complex environments. Given the rich temporal changes that exist in these environments, it is proposed that microorganisms both adapt to and contribute to these daily dynamics through the process of temporal mutualism. Understanding the timekeeping and temporal interactions within these systems will enable a deeper understanding of circadian clocks and temporal programs and provide valuable insights for medicine and agriculture.

Differentiation of Vegetative Cells into Spores: a Kinetic Model Applied to Bacillus subtilis

Spore-forming bacteria are natural contaminants of food raw materials, and sporulation can occur in many environments from farm to fork. In order to characterize and to predict spore formation over time, we developed a model that describes both the kinetics of growth and the differentiation of vegetative cells into spores. The model is based on a classical growth model and enables description of the kinetics of sporulation with the addition of three parameters specific to sporulation. Two parameters are related to the probability of each vegetative cell to commit to sporulation and to form a spore, and the last one is related to the time needed to form a spore once the cell is committed to sporulation. The goodness of fit of this growth-sporulation model was assessed using growth-sporulation kinetics at various temperatures in laboratory medium or in whey for Bacillus subtilis, Bacillus cereus, and Bacillus licheniformis The model accurately describes the kinetics in these different conditions, with a mean error lower than 0.78 log₁₀ CFU/ml for the growth and 1.08 log₁₀ CFU/ml for the sporulation. The biological meaning of the parameters was validated with a derivative strain of Bacillus subtilis 168 which produces green fluorescent protein at the initiation of sporulation. This model provides physiological information on the spore formation and on the temporal abilities of vegetative cells to differentiate into spores and reveals the heterogeneity of spore formation during and after growth.IMPORTANCE The growth-sporulation model describes the progressive transition from vegetative cells to spores with sporulation parameters describing the sporulation potential of each vegetative cell. Consequently, the model constitutes an interesting tool to assess the sporulation potential of a bacterial population over time with accurate parameters such as the time needed to obtain one resistant spore and the probability of sporulation. Further, this model can be used to assess these data under various environmental conditions in order to better identify the conditions favorable for sporulation regarding the time to obtain the first spore and/or the concentrations of spores which could be reached during a food process.

Integrated Transcriptional Regulatory Network of Quorum Sensing, Replication Control, and SOS Response in Dinoroseobacter shibae

Quorum sensing (QS) coordinates population wide gene expression of bacterial species. Highly adaptive traits like gene transfer agents (GTA), morphological heterogeneity, type 4 secretion systems (T4SS), and flagella are QS controlled in Dinoroseobacter shibae, a Roseobacter model organism. Its QS regulatory network is integrated with the CtrA phosphorelay that controls cell division in alphaproteobacteria. To elucidate the network topology, we analyzed the transcriptional response of the QS-negative D. shibae strain ΔluxI₁ toward externally added autoinducer (AI) over a time period of 3 h. The signaling cascade is initiated by the CtrA phosphorelay, followed by the QS genes and other target genes, including the second messenger c-di-GMP, competence, flagella and pili. Identification of transcription factor binding sites in promoters of QS induced genes revealed the integration of QS, CtrA phosphorelay and the SOS stress response mediated by LexA. The concentration of regulatory genes located close to the origin or terminus of replication suggests that gene regulation and replication are tightly coupled. Indeed, addition of AI first stimulates and then represses replication. The restart of replication comes along with increased c-di-GMP levels. We propose a model in which QS induces replication followed by differentiation into GTA producing and non-producing cells. CtrA-activity is controlled by the c-di-GMP level, allowing some of the daughter cells to replicate again. The size of the GTA producing subpopulation is tightly controlled by QS via the AI Synthase LuxI₂. Finally, induction of the SOS response allows for integration of GTA DNA into the host chromosome.

Synthetic biology approaches for chromosomal integration of genes and pathways in industrial microbial systems

Industrial biotechnology is reliant on native pathway engineering or foreign pathway introduction for efficient biosynthesis of target products. Chromosomal integration, with intrinsic genetic stability, is an indispensable step for reliable expression of homologous or heterologous genes and pathways in large-scale and long-term fermentation. With advances in synthetic biology and CRISPR-based genome editing approaches, a wide variety of novel enabling technologies have been developed for single-step, markerless, multi-locus genomic integration of large biochemical pathways, which significantly facilitate microbial overproduction of chemicals, pharmaceuticals and other value-added biomolecules. Notably, the newly discovered homology-mediated end joining strategy could be widely applicable for high-efficiency genomic integration in a number of homologous recombination-deficient microbes. In this review, we explore the fundamental principles and characteristics of genomic integration, and highlight the development and applications of targeted integration approaches in the three representative industrial microbial systems, including Escherichia coli, actinomycetes and yeasts.

Cell size control and gene expression homeostasis in single-cells

Growth of a cell and its subsequent division into daughters is a fundamental aspect of all cellular living systems. During these processes, how do individual cells correct size aberrations so that they do not grow abnormally large or small? How do cells ensure that the concentration of essential gene products are maintained at desired levels, in spite of dynamic/stochastic changes in cell size during growth and division? Both these questions have fascinated researchers for over a century. We review how advances in singe-cell technologies and measurements are providing unique insights into these questions across organisms from prokaryotes to human cells. More specifically, diverse strategies based on timing of cell-cycle events, regulating growth, and number of daughters are employed to maintain cell size homeostasis. Interestingly, size homeostasis often results in size optimality - proliferation of individual cells in a population is maximized at an optimal cell size. We further discuss how size-dependent expression or gene-replication timing can buffer concentration of a gene product from cell-to-cell size variations within a population. Finally, we speculate on an intriguing hypothesis that specific size control strategies may have evolved as a consequence of gene-product concentration homeostasis.

Dynamics of Bacterial Gene Regulatory Networks

The ability of bacterial cells to adjust their gene expression program in response to environmental perturbation is often critical for their survival. Recent experimental advances allowing us to quantitatively record gene expression dynamics in single cells and in populations coupled with mathematical modeling enable mechanistic understanding on how these responses are shaped by the underlying regulatory networks. Here, we review how the combination of local and global factors affect dynamical responses of gene regulatory networks. Our goal is to discuss the general principles that allow extrapolation from a few model bacteria to less understood microbes. We emphasize that, in addition to well-studied effects of network architecture, network dynamics are shaped by global pleiotropic effects and cell physiology.

Reconfigurable Analog Signal Processing by Living Cells

Living cells are known for their capacity for versatile signal processing, particularly the ability to respond differently to the same stimuli using biochemical networks that integrate environmental signals and reconfigure their dynamic responses. However, the complexity of natural biological networks confounds the discovery of fundamental mechanisms behind versatile signaling. Here, we study one specific aspect of reconfigurable signal processing in which a minimal biological network integrates two signals, using one to reconfigure the network's transfer function with respect to the other, producing an emergent switch between induction and repression. In contrast to known mechanisms, the new mechanism reconfigures transfer functions through genetic networks without extensive protein-protein interactions. These results provide a novel explanation for the versatility of genetic programs, and suggest a new mechanism of signal integration that may govern flexibility and plasticity of gene expression.

Late-Arriving Signals Contribute Less to Cell-Fate Decisions

Gene regulatory networks are largely responsible for cellular decision-making. These networks sense diverse external signals and respond by adjusting gene expression, enabling cells to reach environment-dependent decisions crucial for their survival or reproduction. However, information-carrying signals may arrive at variable times. Besides the intrinsic strength of these signals, their arrival time (timing) may also carry information about the environment and can influence cellular decision-making in ways that are poorly understood. For example, it is unclear how the timing of individual phage infections affects the lysis-lysogeny decision of bacteriophage λ despite variable infection times being likely in the wild and even in laboratory conditions. In this work, we combine mathematical modeling with experimentation to address this question. We develop an experimentally testable theory, which reveals that late-infecting phages contribute less to cellular decision-making. This implies that infection delays lower the probability of lysogeny compared to simultaneous infections. Furthermore, we show that infection delays reduce lysogenization by providing insufficient CII for threshold crossing during the critical decision-making period. We find evidence for a cutoff time after which subsequent infections cannot influence the cellular decision. We derive an intuitive formula that approximates the probability of lysogeny for variable infection times by a time-weighted average of probabilities for simultaneous infections. We validate these theoretical predictions experimentally. Similar concepts and simplifying modeling approaches may help elucidate the mechanisms underlying other cellular decisions.

A model of cell-wall dynamics during sporulation in Bacillus subtilis

To survive starvation, Bacillus subtilis forms durable spores. After asymmetric cell division, the septum grows around the forespore in a process called engulfment, but the mechanism of force generation is unknown. Here, we derived a novel biophysical model for the dynamics of cell-wall remodeling during engulfment based on a balancing of dissipative, active, and mechanical forces. By plotting phase diagrams, we predict that sporulation is promoted by a line tension from the attachment of the septum to the outer cell wall, as well as by an imbalance in turgor pressures in the mother-cell and forespore compartments. We also predict that significant mother-cell growth hinders engulfment. Hence, relatively simple physical principles may guide this complex biological process.

Genetic networks controlled by the bacterial replication initiator and transcription factor DnaA in Bacillus subtilis

DnaA is the widely conserved bacterial AAA+ ATPase that functions as both the replication initiator and a transcription factor. In many organisms, DnaA controls expression of its own gene and likely several others during growth and in response to replication stress. To evaluate the effects of DnaA on gene expression, separate from its role in replication initiation, we analyzed changes in mRNA levels in Bacillus subtilis cells with and without dnaA, using engineered strains in which dnaA is not essential. We found that dnaA was required for many of the changes in gene expression in response to replication stress. We also found that dnaA indirectly affected expression of several regulons during growth, including those controlled by the transcription factors Spo0A, AbrB, PhoP, SinR, RemA, Rok and YvrH. These effects were largely mediated by the effects of DnaA on expression of sda. DnaA activates transcription of sda, and Sda inhibits histidine protein kinases required for activation of the transcription factor Spo0A. We also found that loss of dnaA caused a decrease in the development of genetic competence. Together, our results indicate that DnaA plays an important role in modulating cell physiology, separate from its role in replication initiation.

A role for Rtt109 in buffering gene-dosage imbalance during DNA replication

Chromatin can function as an integrator of DNA-related processes, allowing communication, for example, between DNA replication and gene transcription. Such communication is needed to overcome the gene-dosage imbalance introduced during DNA replication, when certain genes are replicated prior to others. Increased transcription of early replicating genes could alter regulatory balances. This does not occur, suggesting a mechanism that suppresses expression from newly replicated DNA. Critical to this buffering is Rtt109, which acetylates the internal K56 residue of newly synthesized histone H3 prior to incorporation onto DNA. H3K56ac distinguishes replicated from non-replicated DNA, communicating this information to the transcription machinery to ensure expression homeostasis during S phase.

The PAS domains of the major sporulation kinase in Bacillus subtilis play a role in tetramer formation that is essential for the autokinase activity

Sporulation in Bacillus subtilis is induced upon starvation. In a widely accepted model, an N‐terminal “sensor” domain of the major sporulation kinase KinA recognizes a hypothetical starvation signal(s) and autophosphorylates a histidine residue to activate the master regulator Spo0A via a multicomponent phosphorelay. However, to date no confirmed signal has been found. Here, we demonstrated that PAS‐A, the most N‐terminal of the three PAS domains (PAS‐ABC), is dispensable for the activity, contrary to a previous report. Our data indicated that the autokinase activity is dependent on the formation of a functional tetramer, which is mediated by, at least, PAS‐B and PAS‐C. Additionally, we ruled out the previously proposed notion that NAD⁺/NADH ratio controls KinA activity through the PAS‐A domain by demonstrating that the cofactors show no effects on the kinase activity in vitro. In support of these data, we found that the cofactors exist in approximately 1000‐fold excess of KinA in the cell and the cofactors’ ratio does not change significantly during growth and sporulation, suggesting that changes in the cofactor ratio might not play a role in controlling KinA activity. These data may refute the widely‐held belief that the activity of KinA is regulated in response to an unknown starvation signal(s).

Image analysis driven single-cell analytics for systems microbiology

BACKGROUND

Time-lapse microscopy is an essential tool for capturing and correlating bacterial morphology and gene expression dynamics at single-cell resolution. However state-of-the-art computational methods are limited in terms of the complexity of cell movies that they can analyze and lack of automation. The proposed Bacterial image analysis driven Single Cell Analytics (BaSCA) computational pipeline addresses these limitations thus enabling high throughput systems microbiology.

RESULTS

BaSCA can segment and track multiple bacterial colonies and single-cells, as they grow and divide over time (cell segmentation and lineage tree construction) to give rise to dense communities with thousands of interacting cells in the field of view. It combines advanced image processing and machine learning methods to deliver very accurate bacterial cell segmentation and tracking (F-measure over 95%) even when processing images of imperfect quality with several overcrowded colonies in the field of view. In addition, BaSCA extracts on the fly a plethora of single-cell properties, which get organized into a database summarizing the analysis of the cell movie. We present alternative ways to analyze and visually explore the spatiotemporal evolution of single-cell properties in order to understand trends and epigenetic effects across cell generations. The robustness of BaSCA is demonstrated across different imaging modalities and microscopy types.

CONCLUSIONS

BaSCA can be used to analyze accurately and efficiently cell movies both at a high resolution (single-cell level) and at a large scale (communities with many dense colonies) as needed to shed light on e.g. how bacterial community effects and epigenetic information transfer play a role on important phenomena for human health, such as biofilm formation, persisters' emergence etc. Moreover, it enables studying the role of single-cell stochasticity without losing sight of community effects that may drive it.

Interrogating the topological robustness of gene regulatory circuits by randomization

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Cells are able to robustly carry out their essential biological functions, possibly because of multiple layers of tight regulation via complex, yet well-designed, gene regulatory networks involving a substantial number of genes. State-of-the-art genomics technology has enabled the mapping of these large gene networks, yet it remains a tremendous challenge to elucidate their design principles and the regulatory mechanisms underlying their biological functions such as signal processing and decision-making. One of the key barriers is the absence of accurate kinetics for the regulatory interactions, especially from in vivo experiments. To this end, we have developed a new computational modeling method, Random Circuit Perturbation (RACIPE), to explore the dynamic behaviors of gene regulatory circuits without the requirement of detailed kinetic parameters. RACIPE takes a network topology as the input, and generates an unbiased ensemble of models with varying kinetic parameters. Each model is subjected to simulation, followed by statistical analysis for the ensemble. We tested RACIPE on several gene circuits, and found that the predicted gene expression patterns from all of the models converge to experimentally observed gene state clusters. We expect RACIPE to be a powerful method to identify the role of network topology in determining network operating principles.

The Proximity of Ribosomal Protein Genes to oriC Enhances Vibrio cholerae Fitness in the Absence of Multifork Replication

Recent works suggest that bacterial gene order links chromosome structure to cell homeostasis. Comparative genomics showed that, in fast-growing bacteria, ribosomal protein genes (RP) locate near the replication origin (oriC). We recently showed that Vibrio cholerae employs this positional bias as a growth optimization strategy: under fast-growth conditions, multifork replication increases RP dosage and expression. However, RP location may provide advantages in a dosage-independent manner: for example, the physical proximity of the many ribosomal components, in the context of a crowded cytoplasm, may favor ribosome biogenesis. To uncover putative dosage-independent effects, we studied isogenic V. cholerae derivatives in which the major RP locus, S10-spc-α (S10), was relocated to alternative genomic positions. When bacteria grew fast, bacterial fitness was reduced according to the S10 relative distance to oriC The growth of wild-type V. cholerae could not be improved by additional copies of the locus, suggesting a physiologically optimized genomic location. Slow growth is expected to uncouple RP position from dosage, since multifork replication does not occur. Under these conditions, we detected a fitness impairment when S10 was far from oriC Deep sequencing followed by marker frequency analysis in the absence of multifork replication revealed an up to 30% S10 dosage reduction associated with its relocation that closely correlated with fitness alterations. Hence, the impact of S10 location goes beyond a growth optimization strategy during feast periods. RP location may be important during the whole life cycle of this pathogen.IMPORTANCE The role of gene order within the bacterial chromosome is poorly understood. In fast growers, the location of genes linked with the expression of genetic information (i.e., transcription and translation) is biased toward oriC It was proposed that the location of these genes helps to maximize their expression by recruiting multifork replication during fast growth. Our results show that such genomic positioning impacts cell fitness beyond fast-growth conditions, probably across the whole life cycle of fast growers. Thus, the genomic position of key highly expressed genes, such as RP, was finely tuned during the evolution of fast-growing bacteria and may also be important in slow growers. In the near future, many more genes whose genomic position impacts bacterial phenotype will be described. These studies will contribute to discovery the rules of genome organization and application of them for the design of synthetic chromosomes and the creation of artificial life forms.

Control of Initiation of DNA Replication in Bacillus subtilis and Escherichia coli

Initiation of DNA Replication is tightly regulated in all cells since imbalances in chromosomal copy number are deleterious and often lethal. In bacteria such as Bacillus subtilis and Escherichia coli, at the point of cytokinesis, there must be two complete copies of the chromosome to partition into the daughter cells following division at mid-cell during vegetative growth. Under conditions of rapid growth, when the time taken to replicate the chromosome exceeds the doubling time of the cells, there will be multiple initiations per cell cycle and daughter cells will inherit chromosomes that are already undergoing replication. In contrast, cells entering the sporulation pathway in B. subtilis can do so only during a short interval in the cell cycle when there are two, and only two, chromosomes per cell, one destined for the spore and one for the mother cell. Here, we briefly describe the overall process of DNA replication in bacteria before reviewing initiation of DNA replication in detail. The review covers DnaA-directed assembly of the replisome at oriC and the multitude of mechanisms of regulation of initiation, with a focus on the similarities and differences between E. coli and B. subtilis.

Role of Autoregulation and Relative Synthesis of Operon Partners in Alternative Sigma Factor Networks

Despite the central role of alternative sigma factors in bacterial stress response and virulence their regulation remains incompletely understood. Here we investigate one of the best-studied examples of alternative sigma factors: the σ^B network that controls the general stress response of Bacillus subtilis to uncover widely relevant general design principles that describe the structure-function relationship of alternative sigma factor regulatory networks. We show that the relative stoichiometry of the synthesis rates of σ^B, its anti-sigma factor RsbW and the anti-anti-sigma factor RsbV plays a critical role in shaping the network behavior by forcing the σ^B network to function as an ultrasensitive negative feedback loop. We further demonstrate how this negative feedback regulation insulates alternative sigma factor activity from competition with the housekeeping sigma factor for RNA polymerase and allows multiple stress sigma factors to function simultaneously with little competitive interference.

Understanding the regulation of bacterial stress response holds the key to tackling the problems of emerging resistance to anti-bacteria’s and antibiotics. To this end, here we study one of the longest serving model systems of bacterial stress response: the σ^B pathway of Bacillus subtilis. The sigma factor σ^B controls the general stress response of Bacillus subtilis to a variety of stress conditions including starvation, antibiotics and harmful environmental perturbations. Recent studies have demonstrated that an increase in stress triggers pulsatile activation of σ^B. Using mathematical modeling we identify the core structural design feature of the network that are responsible for its pulsatile response. We further demonstrate how the same core design features are common to a variety of stress response pathways. As a result of these features, cells can respond to multiple simultaneous stresses without interference or competition between the different pathways.

Effect of Genome Position on Heterologous Gene Expression in Bacillus subtilis: An Unbiased Analysis

A fixed gene copy number is important for the in silico construction of engineered synthetic networks. However, the copy number of integrated genes depends on their genomic location. This gene dosage effect is rarely addressed in synthetic biology. Two studies in Escherichia coli presented conflicting data on the impact of gene dosage. Here, we investigate how genome location and gene orientation influences expression in Bacillus subtilis. An important difference with the E. coli studies is that we used an unbiased genome integration approach mediated by random transposon insertion. We found that there is a strong gene dosage effect in fast growing B. subtilis cells, which can amount to a 5-fold difference in gene expression. In contrast, gene orientation with respect to DNA replication direction does not influence gene expression. Our study shows that gene dosage should be taken into account when designing synthetic circuits in B. subtilis and presumably other bacteria.