Transcriptomic and phylogenetic analysis of a bacterial cell cycle reveals strong associations between gene co-expression and evolution

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.

CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

The noncoding RNA CcnA modulates the master cell cycle regulators CtrA and GcrA in Caulobacter crescentus

Bacteria are powerful models for understanding how cells divide and accomplish global regulatory programs. In Caulobacter crescentus, a cascade of essential master regulators supervises the correct and sequential activation of DNA replication, cell division, and development of different cell types. Among them, the response regulator CtrA plays a crucial role coordinating all those functions. Here, for the first time, we describe the role of a novel factor named CcnA (cell cycle noncoding RNA A), a cell cycle–regulated noncoding RNA (ncRNA) located at the origin of replication, presumably activated by CtrA, and responsible for the accumulation of CtrA itself. In addition, CcnA may be also involved in the inhibition of translation of the S-phase regulator, GcrA, by interacting with its 5′ untranslated region (5′ UTR). Performing in vitro experiments and mutagenesis, we propose a mechanism of action of CcnA based on liberation (ctrA) or sequestration (gcrA) of their ribosome-binding site (RBS). Finally, its role may be conserved in other alphaproteobacterial species, such as Sinorhizobium meliloti, representing indeed a potentially conserved process modulating cell cycle in Caulobacterales and Rhizobiales.

During cell cycle progression in the bacterium Caulobacter crescentus, the master cell cycle regulator CtrA is controlled by CcnA, a cell cycle-regulated non-coding RNA transcribed from a gene located at the origin of replication.

The Lon protease temporally restricts polar cell differentiation events during the Caulobacter cell cycle

The highly conserved protease Lon has important regulatory and protein quality control functions in cells from the three domains of life. Despite many years of research on Lon, only a few specific protein substrates are known in most organisms. Here, we used a quantitative proteomics approach to identify novel substrates of Lon in the dimorphic bacterium Caulobacter crescentus. We focused our study on proteins involved in polar cell differentiation and investigated the developmental regulator StaR and the flagella hook length regulator FliK as specific Lon substrates in detail. We show that Lon recognizes these proteins at their C-termini, and that Lon-dependent degradation ensures their temporally restricted accumulation in the cell cycle phase when their function is needed. Disruption of this precise temporal regulation of StaR and FliK levels in a Δlon mutant contributes to defects in stalk biogenesis and motility, respectively, revealing a critical role of Lon in coordinating developmental processes with cell cycle progression. Our work underscores the importance of Lon in the regulation of complex temporally controlled processes by adjusting the concentrations of critical regulatory proteins. Furthermore, this study includes the first characterization of FliK in C. crescentus and uncovers a dual role of the C-terminal amino acids of FliK in protein function and degradation.

Promoter Architecture Differences among Alphaproteobacteria and Other Bacterial Taxa

Much of our knowledge of bacterial transcription initiation has been derived from studying the promoters of Escherichia coli and Bacillus subtilis. Given the expansive diversity across the bacterial phylogeny, it is unclear how much of this knowledge can be applied to other organisms. Here, we report on bioinformatic analyses of promoter sequences of the primary σ factor (σ70) by leveraging publicly available transcription start site (TSS) sequencing data sets for nine bacterial species spanning five phyla. This analysis identifies previously unreported differences in the -35 and -10 elements of σ70-dependent promoters in several groups of bacteria. We found that Actinobacteria and Betaproteobacteria σ70-dependent promoters lack the TTG triad in their -35 element, which is predicted to be conserved across the bacterial phyla. In addition, the majority of the Alphaproteobacteria σ70-dependent promoters analyzed lacked the thymine at position -7 that is highly conserved in other phyla. Bioinformatic examination of the Alphaproteobacteria σ70-dependent promoters identifies a significant overrepresentation of essential genes and ones encoding proteins with common cellular functions downstream of promoters containing an A, C, or G at position -7. We propose that transcription of many σ70-dependent promoters in Alphaproteobacteria depends on the transcription factor CarD, which is an essential protein in several members of this phylum. Our analysis expands the knowledge of promoter architecture across the bacterial phylogeny and provides new information that can be used to engineer bacteria for use in medical, environmental, agricultural, and biotechnological processes. IMPORTANCE Transcription of DNA to RNA by RNA polymerase is essential for cells to grow, develop, and respond to stress. Understanding the process and control of transcription is important for health, disease, the environment, and biotechnology. Decades of research on a few bacteria have identified promoter DNA sequences that are recognized by the σ subunit of RNA polymerase. We used bioinformatic analyses to reveal previously unreported differences in promoter DNA sequences across the bacterial phylogeny. We found that many Actinobacteria and Betaproteobacteria promoters lack a sequence in their -35 DNA recognition element that was previously assumed to be conserved and that Alphaproteobacteria lack a thymine residue at position -7, also previously assumed to be conserved. Our work reports important new information about bacterial transcription, illustrates the benefits of studying bacteria across the phylogenetic tree, and proposes new lines of future investigation.

Environmental Conditions Modulate the Transcriptomic Response of Both Caulobacter crescentus Morphotypes to Cu Stress

Bacteria encounter elevated copper (Cu) concentrations in multiple environments, varying from mining wastes to antimicrobial applications of copper. As the role of the environment in the bacterial response to Cu ion exposure remains elusive, we used a tagRNA-seq approach to elucidate the disparate responses of two morphotypes of Caulobacter crescentus NA1000 to moderate Cu stress in a complex rich (PYE) medium and a defined poor (M2G) medium. The transcriptome was more responsive in M2G, where we observed an extensive oxidative stress response and reconfiguration of the proteome, as well as the induction of metal resistance clusters. In PYE, little evidence was found for an oxidative stress response, but several transport systems were differentially expressed, and an increased need for histidine was apparent. These results show that the Cu stress response is strongly dependent on the cellular environment. In addition, induction of the extracytoplasmic function sigma factor SigF and its regulon was shared by the Cu stress responses in both media, and its central role was confirmed by the phenotypic screening of a sigF::Tn5 mutant. In both media, stalked cells were more responsive to Cu stress than swarmer cells, and a stronger basal expression of several cell protection systems was noted, indicating that the swarmer cell is inherently more Cu resistant. Our approach also allowed for detecting several new transcription start sites, putatively indicating small regulatory RNAs, and additional levels of Cu-responsive regulation.

The Protein Quality Control Network in Caulobacter crescentus

The asymmetric life cycle of Caulobacter crescentus has provided a model in which to study how protein quality control (PQC) networks interface with cell cycle and developmental processes, and how the functions of these systems change during exposure to stress. As in most bacteria, the PQC network of Caulobacter contains highly conserved ATP-dependent chaperones and proteases as well as more specialized holdases. During growth in optimal conditions, these systems support a regulated circuit of protein synthesis and degradation that drives cell differentiation and cell cycle progression. When stress conditions threaten the proteome, most components of the Caulobacter proteostasis network are upregulated and switch to survival functions that prevent, revert, and remove protein damage, while simultaneously pausing the cell cycle in order to regain protein homeostasis. The specialized physiology of Caulobacter influences how it copes with proteotoxic stress, such as in the global management of damaged proteins during recovery as well as in cell type-specific stress responses. Our mini-review highlights the discoveries that have been made in how Caulobacter utilizes its PQC network for regulating its life cycle under optimal and proteotoxic stress conditions, and discusses open research questions in this model.

A first perturbome of Pseudomonas aeruginosa: Identification of core genes related to multiple perturbations by a machine learning approach

Tolerance to stress conditions is vital for organismal survival, including bacteria under specific environmental conditions, antibiotics, and other perturbations. Some studies have described common modulation and shared genes during stress response to different types of disturbances (termed as perturbome), leading to the idea of central control at the molecular level. We implemented a robust machine learning approach to identify and describe genes associated with multiple perturbations or perturbome in a Pseudomonas aeruginosa PAO1 model. Using microarray datasets from the Gene Expression Omnibus (GEO), we evaluated six approaches to rank and select genes: using two methodologies, data single partition (SP method) or multiple partitions (MP method) for training and testing datasets, we evaluated three classification algorithms (SVM Support Vector Machine, KNN K-Nearest neighbor and RF Random Forest). Gene expression patterns and topological features at the systems level were included to describe the perturbome elements. We were able to select and describe 46 core response genes associated with multiple perturbations in P. aeruginosa PAO1 and it can be considered a first report of the P. aeruginosa perturbome. Molecular annotations, patterns in expression levels, and topological features in molecular networks revealed biological functions of biosynthesis, binding, and metabolism, many of them related to DNA damage repair and aerobic respiration in the context of tolerance to stress. We also discuss different issues related to implemented and assessed algorithms, including data partitioning, classification approaches, and metrics. Altogether, this work offers a different and robust framework to select genes using a machine learning approach.

Transcriptomic determinants of the response of ST-111 Pseudomonas aeruginosa AG1 to ciprofloxacin identified by a top-down systems biology approach

Pseudomonas aeruginosa is an opportunistic pathogen that thrives in diverse environments and causes a variety of human infections. Pseudomonas aeruginosa AG1 (PaeAG1) is a high-risk sequence type 111 (ST-111) strain isolated from a Costa Rican hospital in 2010. PaeAG1 has both blaVIM-2 and blaIMP-18 genes encoding for metallo-β-lactamases, and it is resistant to β-lactams (including carbapenems), aminoglycosides, and fluoroquinolones. Ciprofloxacin (CIP) is an antibiotic commonly used to treat P. aeruginosa infections, and it is known to produce DNA damage, triggering a complex molecular response. In order to evaluate the effects of a sub-inhibitory CIP concentration on PaeAG1, growth curves using increasing CIP concentrations were compared. We then measured gene expression using RNA-Seq at three time points (0, 2.5 and 5 h) after CIP exposure to identify the transcriptomic determinants of the response (i.e. hub genes, gene clusters and enriched pathways). Changes in expression were determined using differential expression analysis and network analysis using a top–down systems biology approach. A hybrid model using database-based and co-expression analysis approaches was implemented to predict gene–gene interactions. We observed a reduction of the growth curve rate as the sub-inhibitory CIP concentrations were increased. In the transcriptomic analysis, we detected that over time CIP treatment resulted in the differential expression of 518 genes, showing a complex impact at the molecular level. The transcriptomic determinants were 14 hub genes, multiple gene clusters at different levels (associated to hub genes or as co-expression modules) and 15 enriched pathways. Down-regulation of genes implicated in several metabolism pathways, virulence elements and ribosomal activity was observed. In contrast, amino acid catabolism, RpoS factor, proteases, and phenazines genes were up-regulated. Remarkably, > 80 resident-phage genes were up-regulated after CIP treatment, which was validated at phenomic level using a phage plaque assay. Thus, reduction of the growth curve rate and increasing phage induction was evidenced as the CIP concentrations were increased. In summary, transcriptomic and network analyses, as well as the growth curves and phage plaque assays provide evidence that PaeAG1 presents a complex, concentration-dependent response to sub-inhibitory CIP exposure, showing pleiotropic effects at the systems level. Manipulation of these determinants, such as phage genes, could be used to gain more insights about the regulation of responses in PaeAG1 as well as the identification of possible therapeutic targets. To our knowledge, this is the first report of the transcriptomic analysis of CIP response in a ST-111 high-risk P. aeruginosa strain, in particular using a top-down systems biology approach.

Untargeted metabolomics links glutathione to bacterial cell cycle progression

Cell cycle progression requires the coordination of cell growth, chromosome replication, and division. Consequently, a functional cell cycle must be coupled with metabolism. However, direct measurements of metabolome dynamics remained scarce, in particular in bacteria. Here, we describe an untargeted metabolomics approach with synchronized Caulobacter crescentus cells to monitor the relative abundance changes of ~400 putative metabolites as a function of the cell cycle. While the majority of metabolite pools remains homeostatic, ~14% respond to cell cycle progression. In particular, sulfur metabolism is redirected during the G1-S transition, and glutathione levels periodically change over the cell cycle with a peak in late S phase. A lack of glutathione perturbs cell size by uncoupling cell growth and division through dysregulation of KefB, a K+/H+ antiporter. Overall, we here describe the impact of the C. crescentus cell cycle progression on metabolism, and in turn relate glutathione and potassium homeostasis to timely cell division.

The Two Chemotaxis Clusters in Caulobacter crescentus Play Different Roles in Chemotaxis and Biofilm Regulation

The holdfast polysaccharide adhesin is crucial for irreversible cell adhesion and biofilm formation in Caulobacter crescentus Holdfast production is tightly controlled via developmental regulators, as well as via environmental and physical signals. Here, we identify a novel mode of regulation of holdfast synthesis that involves chemotaxis proteins. We characterized the two identified chemotaxis clusters of C. crescentus and showed that only the previously characterized major cluster is involved in the chemotactic response toward different carbon sources. However, both chemotaxis clusters encoded in the C. crescentus genome play a role in biofilm formation and holdfast production by regulating the expression of hfiA, the gene encoding the holdfast inhibitor HfiA. We show that CheA and CheB proteins act in an antagonistic manner, as follows: while the two CheA proteins negatively regulate hfiA expression, the CheB proteins are positive regulators, thus providing a modulation of holdfast synthesis and surface attachment.IMPORTANCE Chemosensory systems constitute major signal transduction pathways in bacteria. These systems are involved in chemotaxis and other cell responses to environment conditions, such as the production of adhesins to enable irreversible adhesion to a surface and surface colonization. The C. crescentus genome encodes two complete chemotaxis clusters. Here, we characterized the second novel chemotaxis-like cluster. While only the major chemotaxis cluster is involved in chemotaxis, both chemotaxis systems modulate C. crescentus adhesion by controlling expression of the holdfast synthesis inhibitor HfiA. Here, we identify a new level in holdfast regulation, providing new insights into the control of adhesin production that leads to the formation of biofilms in response to the environment.

High-resolution 3D models of Caulobacter crescentus chromosome reveal genome structural variability and organization

High-resolution three-dimensional models of Caulobacter crescentus nucleoid structures were generated via a multi-scale modeling protocol. Models were built as a plectonemically supercoiled circular DNA and by incorporating chromosome conformation capture based data to generate an ensemble of base pair resolution models consistent with the experimental data. Significant structural variability was found with different degrees of bending and twisting but with overall similar topologies and shapes that are consistent with C. crescentus cell dimensions. The models allowed a direct mapping of the genomic sequence onto the three-dimensional nucleoid structures. Distinct spatial distributions were found for several genomic elements such as AT-rich sequence elements where nucleoid associated proteins (NAPs) are likely to bind, promoter sites, and some genes with common cellular functions. These findings shed light on the correlation between the spatial organization of the genome and biological functions.

Periplasmic protein EipA determines envelope stress resistance and virulence in Brucella abortus

Molecular components of the Brucella abortus cell envelope play a major role in its ability to infect, colonize and survive inside mammalian host cells. In this study, we have defined a role for a conserved gene of unknown function in B. abortus envelope stress resistance and infection. Expression of this gene, which we name eipA, is directly activated by the essential cell cycle regulator, CtrA. eipA encodes a soluble periplasmic protein that adopts an unusual eight-stranded β-barrel fold. Deletion of eipA attenuates replication and survival in macrophage and mouse infection models, and results in sensitivity to treatments that compromise the cell envelope integrity. Transposon disruption of genes required for LPS O-polysaccharide biosynthesis is synthetically lethal with eipA deletion. This genetic connection between O-polysaccharide and eipA is corroborated by our discovery that eipA is essential in Brucella ovis, a naturally rough species that harbors mutations in several genes required for O-polysaccharide production. Conditional depletion of eipA expression in B. ovis results in a cell chaining phenotype, providing evidence that eipA directly or indirectly influences cell division in Brucella. We conclude that EipA is a molecular determinant of Brucella virulence that functions to maintain cell envelope integrity and influences cell division.

Cell division control in Caulobacter crescentus

Caulobacter crescentus is a free-living Alphaproteobacterium that thrives in oligotrophic environments. This review focuses on the regulatory network used by this bacterium to control the levels of cell division proteins, their organization inside the cell and their activity as a function of the cell cycle. Strikingly, C. crescentus makes frequent use of master transcriptional regulators and epigenetic signals, most likely to synchronize cell division with other events of the cell cycle. In addition, cellular metabolism and DNA damage sensors emerge as central players regulating cell division in response to changing environmental conditions.

Gene network analysis identifies a central post-transcriptional regulator of cellular stress survival

Cells adapt to shifts in their environment by remodeling transcription. Measuring changes in transcription at the genome scale is now routine, but defining the functional significance of individual genes within large gene expression datasets remains a major challenge. We applied a network-based algorithm to interrogate publicly available gene expression data to predict genes that serve major functional roles in Caulobacter crescentus stress survival. This approach identified GsrN, a conserved small RNA that is directly activated by the general stress sigma factor, σ^T, and functions as a potent post-transcriptional regulator of survival across distinct conditions including osmotic and oxidative stress. Under hydrogen peroxide stress, GsrN protects cells by base pairing with the leader of katG mRNA and activating expression of KatG catalase/peroxidase protein. We conclude that GsrN convenes a post-transcriptional layer of gene expression that serves a central functional role in Caulobacter stress physiology.

More than a Tad: spatiotemporal control of Caulobacter pili

The Type IV pilus (T4P) is a powerful and sophisticated bacterial nanomachine involved in numerous cellular processes, including adhesion, DNA uptake and motility. Aside from the well-described subtype T4aP of the Gram-negative genera, including Myxococcus, Pseudomonas and Neisseria, the Tad (tight adherence) pilus secretion system re-shuffles homologous parts from other secretion systems along with uncharacterized components into a new type of protein translocation apparatus. A representative of the Tad apparatus, the Caulobacter crescentus pilus assembly (Cpa) machine is built exclusively at the newborn cell pole once per cell cycle. Recent comprehensive genetic analyses unearthed a myriad of spatiotemporal determinants acting on the Tad/Cpa system, many of which are conserved in other α-proteobacteria, including obligate intracellular pathogens and symbionts.

Systems level analysis of the Chlamydomonas reinhardtii metabolic network reveals variability in evolutionary co-conservation

Metabolic networks, which are mathematical representations of organismal metabolism, are reconstructed to provide computational platforms to guide metabolic engineering experiments and explore fundamental questions on metabolism. Systems level analyses, such as interrogation of phylogenetic relationships within the network, can provide further guidance on the modification of metabolic circuitries. Chlamydomonas reinhardtii, a biofuel relevant green alga that has retained key genes with plant, animal, and protist affinities, serves as an ideal model organism to investigate the interplay between gene function and phylogenetic affinities at multiple organizational levels. Here, using detailed topological and functional analyses, coupled with transcriptomics studies on a metabolic network that we have reconstructed for C. reinhardtii, we show that network connectivity has a significant concordance with the co-conservation of genes; however, a distinction between topological and functional relationships is observable within the network. Dynamic and static modes of co-conservation were defined and observed in a subset of gene-pairs across the network topologically. In contrast, genes with predicted synthetic interactions, or genes involved in coupled reactions, show significant enrichment for both shorter and longer phylogenetic distances. Based on our results, we propose that the metabolic network of C. reinhardtii is assembled with an architecture to minimize phylogenetic profile distances topologically, while it includes an expansion of such distances for functionally interacting genes. This arrangement may increase the robustness of C. reinhardtii's network in dealing with varied environmental challenges that the species may face. The defined evolutionary constraints within the network, which identify important pairings of genes in metabolism, may offer guidance on synthetic biology approaches to optimize the production of desirable metabolites.

Phylogenomic analysis of gene co-expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large-scale genomic and phylogenetic data with gene co-expression networks to extensively study the evolutionary make-up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co-expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de.

A New Essential Cell Division Protein in Caulobacter crescentus

Bacterial cell division is a complex process that relies on a multiprotein complex composed of a core of widely conserved and generally essential proteins and on accessory proteins that vary in number and identity in different bacteria. The assembly of this complex and, particularly, the initiation of constriction are regulated processes that have come under intensive study. In this work, we characterize the function of DipI, a protein conserved in Alphaproteobacteria and Betaproteobacteria that is essential in Caulobacter crescentus Our results show that DipI is a periplasmic protein that is recruited late to the division site and that it is required for the initiation of constriction. The recruitment of the conserved cell division proteins is not affected by the absence of DipI, but localization of DipI to the division site occurs only after a mature divisome has formed. Yeast two-hybrid analysis showed that DipI strongly interacts with the FtsQLB complex, which has been recently implicated in regulating constriction initiation. A possible role of DipI in this process is discussed.IMPORTANCE Bacterial cell division is a complex process for which most bacterial cells assemble a multiprotein complex that consists of conserved proteins and of accessory proteins that differ among bacterial groups. In this work, we describe a new cell division protein (DipI) present only in a group of bacteria but essential in Caulobacter crescentus Cells devoid of DipI cannot constrict. Although a mature divisome is required for DipI recruitment, DipI is not needed for recruiting other division proteins. These results, together with the interaction of DipI with a protein complex that has been suggested to regulate cell wall synthesis during division, suggest that DipI may be part of the regulatory mechanism that controls constriction initiation.

Cohesive Properties of the Caulobacter crescentus Holdfast Adhesin Are Regulated by a Novel c-di-GMP Effector Protein

When encountering surfaces, many bacteria produce adhesins to facilitate their initial attachment and to irreversibly glue themselves to the solid substrate. A central molecule regulating the processes of this motile-sessile transition is the second messenger c-di-GMP, which stimulates the production of a variety of exopolysaccharide adhesins in different bacterial model organisms. In Caulobacter crescentus, c-di-GMP regulates the synthesis of the polar holdfast adhesin during the cell cycle, yet the molecular and cellular details of this control are currently unknown. Here we identify HfsK, a member of a versatile N-acetyltransferase family, as a novel c-di-GMP effector involved in holdfast biogenesis. Cells lacking HfsK form highly malleable holdfast structures with reduced adhesive strength that cannot support surface colonization. We present indirect evidence that HfsK modifies the polysaccharide component of holdfast to buttress its cohesive properties. HfsK is a soluble protein but associates with the cell membrane during most of the cell cycle. Coincident with peak c-di-GMP levels during the C. crescentus cell cycle, HfsK relocalizes to the cytosol in a c-di-GMP-dependent manner. Our results indicate that this c-di-GMP-mediated dynamic positioning controls HfsK activity, leading to its inactivation at high c-di-GMP levels. A short C-terminal extension is essential for the membrane association, c-di-GMP binding, and activity of HfsK. We propose a model in which c-di-GMP binding leads to the dispersal and inactivation of HfsK as part of holdfast biogenesis progression.

Exopolysaccharide (EPS) adhesins are important determinants of bacterial surface colonization and biofilm formation. Biofilms are a major cause of chronic infections and are responsible for biofouling on water-exposed surfaces. To tackle these problems, it is essential to dissect the processes leading to surface colonization at the molecular and cellular levels. Here we describe a novel c-di-GMP effector, HfsK, that contributes to the cohesive properties and stability of the holdfast adhesin in C. crescentus. We demonstrate for the first time that c-di-GMP, in addition to its role in the regulation of the rate of EPS production, also modulates the physicochemical properties of bacterial adhesins. By demonstrating how c-di-GMP coordinates the activity and subcellular localization of HfsK, we provide a novel understanding of the cellular processes involved in adhesin biogenesis control. Homologs of HfsK are found in representatives of different bacterial phyla, suggesting that they play important roles in various EPS synthesis systems.

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.

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

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

Caulobacter crescentus intrinsic dimorphism provides a prompt bimodal response to copper stress

Stress response to fluctuating environments often implies a time-consuming reprogramming of gene expression. In bacteria, the so-called bet hedging strategy, which promotes phenotypic stochasticity within a cell population, is the only fast stress response described so far(1). Here, we show that Caulobacter crescentus asymmetrical cell division allows an immediate bimodal response to a toxic metals-rich environment by allocating specific defence strategies to morphologically and functionally distinct siblings. In this context, a motile swarmer cell favours negative chemotaxis to flee from a copper source, whereas a sessile stalked sibling engages a ready-to-use PcoAB copper homeostasis system, providing evidence of a prompt stress response through intrinsic bacterial dimorphism.

The bacterial cell cycle regulator GcrA is a σ70 cofactor that drives gene expression from a subset of methylated promoters

Haakonsen et al. find that the essential cell cycle regulator GcrA in Caulobacter crescentus forms a stable complex with RNA polymerase and localizes to almost all active σ⁷⁰-dependent promoters in vivo but activates transcription primarily at promoters harboring certain DNA methylation sites. GcrA could stabilize RNA polymerase binding and directly stimulate open complex formation to activate transcription.

Cell cycle progression in most organisms requires tightly regulated programs of gene expression. The transcription factors involved typically stimulate gene expression by binding specific DNA sequences in promoters and recruiting RNA polymerase. Here, we found that the essential cell cycle regulator GcrA in Caulobacter crescentus activates the transcription of target genes in a fundamentally different manner. GcrA forms a stable complex with RNA polymerase and localizes to almost all active σ⁷⁰-dependent promoters in vivo but activates transcription primarily at promoters harboring certain DNA methylation sites. Whereas most transcription factors that contact σ⁷⁰ interact with domain 4, GcrA interfaces with domain 2, the region that binds the −10 element during strand separation. Using kinetic analyses and a reconstituted in vitro transcription assay, we demonstrated that GcrA can stabilize RNA polymerase binding and directly stimulate open complex formation to activate transcription. Guided by these studies, we identified a regulon of ∼200 genes, providing new insight into the essential functions of GcrA. Collectively, our work reveals a new mechanism for transcriptional regulation, and we discuss the potential benefits of activating transcription by promoting RNA polymerase isomerization rather than recruitment exclusively.

Cell cycle control in Alphaproteobacteria

Alphaproteobacteria include many medically and environmentally important organisms. Despite the diversity of their niches and lifestyles, from free-living to host-associated, they usually rely on very similar mechanisms to control their cell cycles. Studies on Caulobacter crescentus still lay the foundation for understanding the molecular details of pathways regulating DNA replication and cell division and coordinating these two processes with other events of the cell cycle. This review highlights recent discoveries on the regulation and the mode of action of conserved global regulators and small molecules like c-di-GMP and (p)ppGpp, which play key roles in cell cycle control. It also describes several newly identified mechanisms that modulate cell cycle progression in response to stresses or environmental conditions.

Fundamental constraints on the abundances of chemotaxis proteins

Flagellated bacteria, such as Escherichia coli, perform directed motion in gradients of concentration of attractants and repellents in a process called chemotaxis. The E. coli chemotaxis signaling pathway is a model for signal transduction, but it has unique features. We demonstrate that the need for fast signaling necessitates high abundances of the proteins involved in this pathway. We show that further constraints on the abundances of chemotaxis proteins arise from the requirements of self-assembly both of flagellar motors and of chemoreceptor arrays. All these constraints are specific to chemotaxis, and published data confirm that chemotaxis proteins tend to be more highly expressed than their homologs in other pathways. Employing a chemotaxis pathway model, we show that the gain of the pathway at the level of the response regulator CheY increases with overall chemotaxis protein abundances. This may explain why, at least in one E. coli strain, the abundance of all chemotaxis proteins is higher in media with lower nutrient content. We also demonstrate that the E. coli chemotaxis pathway is particularly robust to abundance variations of the motor protein FliM.

Mapping the CgrA regulon of Rhodospirillum centenum reveals a hierarchal network controlling Gram-negative cyst development

Background

Several Gram-negative species undergo development leading to the formation of metabolically dormant desiccation resistant cysts. Recent analysis of cyst development has revealed that ~20 % of the Rhodospirillum centenum transcriptome undergo temporal changes in expression as cells transition from vegetative to cyst forms. It has also been established that one trigger for cyst formation is the synthesis of the signaling nucleotide 3‘, 5‘- cyclic guanosine monophosphate (cGMP) that is sensed by a homolog of the catabolite repressor protein called CgrA. CgrA in the presence of cGMP initiate a cascade of gene expression leading to the development of cysts.

Results

In this study, we have used RNA-seq and chromatin immunoprecipitation (ChIP-Seq) techniques to define the CgrA-cGMP regulon. Our results indicate that disruption of CgrA leads to altered expression of 258 genes, 131 of which have been previously reported to be involved in cyst development. ChIP-seq analysis combined with transcriptome data also demonstrates that CgrA directly regulates the expression of numerous sigma factors and transcription factors several of which are known to be involved in cyst cell development.

Conclusions

This analysis reveals the presence of CgrA binding sites upstream of many developmentally regulated genes including many transcription factors and signal transduction components. CgrA thus functions as master controller of the cyst development by initiating a hierarchal cascade of downstream transcription factors that induces temporal expression of encystment genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2248-z) contains supplementary material, which is available to authorized users.

CauloBrowser: A systems biology resource for Caulobacter crescentus

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

Transcriptome analysis of cyst formation in Rhodospirillum centenum reveals large global changes in expression during cyst development

Background

Rhodospirillum centenum is a photosynthetic member of the Gram-negative Azospirillum clade members of which exhibit a complex developmental life-cycle featuring morphologically distinct cell types. Under periods of nutrient deprivation, replicative vegetative cells differentiate into metabolically dormant cysts that survive harsh environmental stresses such as desiccation. Encystment involves a multi-stage developmental process that includes the rounding of cells, production of large intracellular storage granules of poly-hydroxybutyrate (PHB) and the excretion of a protective exopolysaccharide coating that envelops dormant cysts.

Results

To study the process of cyst development, we performed RNA-seq studies on cells that were induced to undergo cyst development. To assay for temporal changes in gene expression, RNA was extracted at 4, 24, 48, 72, 96 hours during development and subjected to deep sequence analysis. These results show that 812 genes exhibit log₂ ≥ 1.5-fold changes in expression over a 96 hour cyst induction period demonstrating large global changes in gene expression during cyst development.

Conclusions

Notable changes in expression occurred in numerous genes involved in cell wall and lipid biosynthesis, metabolic enzymes, and numerous regulatory genes such as histidine kinases and transcription factors. Many genes involved in protein synthesis and DNA replication were also significantly reduced during late stages of cyst development. Genes previously identified by genetic screens as being critical for cyst development also exhibited changes of expression during cyst induction. This study provides the first transcriptome profile of global changes in gene expression that occur during development of cysts in a Gram-negative species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1250-9) contains supplementary material, which is available to authorized users.

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.

Versatility of global transcriptional regulators in alpha-Proteobacteria: from essential cell cycle control to ancillary functions

Recent data indicate that cell cycle transcription in many alpha-Proteobacteria is executed by at least three conserved functional modules in which pairs of antagonistic regulators act jointly, rather than in isolation, to control transcription in S-, G2- or G1-phase. Inactivation of module components often results in pleiotropic defects, ranging from cell death and impaired cell division to fairly benign deficiencies in motility. Expression of module components can follow systemic (cell cycle) or external (nutritional/cell density) cues and may be implemented by auto-regulation, ancillary regulators or other (unknown) mechanisms. Here, we highlight the recent progress in understanding the molecular events and the genetic relationships of the module components in environmental, pathogenic and/or symbiotic alpha-proteobacterial genera. Additionally, we take advantage of the recent genome-wide transcriptional analyses performed in the model alpha-Proteobacterium Caulobacter crescentus to illustrate the complexity of the interactions of the global regulators at selected cell cycle-regulated promoters and we detail the consequences of (mis-)expression when the regulators are absent. This review thus provides the first detailed mechanistic framework for understanding orthologous operational principles acting on cell cycle-regulated promoters in other alpha-Proteobacteria.

Release of nonstop ribosomes is essential

Bacterial ribosomes frequently translate to the 3′ end of an mRNA without terminating at a stop codon. Almost all bacteria use the transfer-messenger RNA (tmRNA)-based trans-translation pathway to release these “nonstop” ribosomes and maintain protein synthesis capacity. trans-translation is essential in some species, but in others, such as Caulobacter crescentus, trans-translation can be inactivated. To determine why trans-translation is dispensable in C. crescentus, a Tn-seq screen was used to identify genes that specifically alter growth in cells lacking ssrA, the gene encoding tmRNA. One of these genes, CC1214, was essential in ΔssrA cells. Purified CC1214 protein could release nonstop ribosomes in vitro. CC1214 is a homolog of the Escherichia coli ArfB protein, and using the CC1214 sequence, ArfB homologs were identified in the majority of bacterial phyla. Most species in which ssrA has been deleted contain an ArfB homolog, suggesting that release of nonstop ribosomes may be essential in most or all bacteria.

Genes that are conserved across large phylogenetic distances are expected to confer a selective advantage. The genes required for trans-translation, ssrA and smpB, have been found in >99% of sequenced bacterial genomes, suggesting that they are broadly important. However, these genes can be deleted in some species without loss of viability. The identification and characterization of C. crescentus ArfB reveals why trans-translation is not essential in C. crescentus and suggests that many other bacteria are likely to use ArfB to survive when trans-translation is compromised.

Localization of the outer membrane protein OmpA2 in Caulobacter crescentus depends on the position of the gene in the chromosome

The outer membrane of Gram-negative bacteria is an essential structure involved in nutrient uptake, protection against harmful substances, and cell growth. Different proteins keep the outer membrane from blebbing out by simultaneously interacting with it and with the cell wall. These proteins have been mainly studied in enterobacteria, where OmpA and the Braun and Pal lipoproteins stabilize the outer membrane. Some degree of functional redundancy exists between these proteins, since none of them is essential but the absence of two of them results in a severe phenotype. Caulobacter crescentus has a different strategy to maintain its outer membrane, since it lacks the Braun lipoprotein and Pal is essential. In this work, we characterized OmpA2, an OmpA-like protein, in this bacterium. Our results showed that this protein is required for normal stalk growth and that it plays a minor role in the stability of the outer membrane. An OmpA2 fluorescent fusion protein showed that the concentration of this protein decreases from the stalk to the new pole. This localization pattern is important for its function, and it depends on the position of the gene locus in the chromosome and, as a consequence, in the cell. This result suggests that little diffusion occurs from the moment that the gene is transcribed until the mature protein attaches to the cell wall in the periplasm. This mechanism reveals the integration of different levels of information from protein function down to genome arrangement that allows the cell to self-organize.

Correction of the Caulobacter crescentus NA1000 genome annotation

Bacterial genome annotations are accumulating rapidly in the GenBank database and the use of automated annotation technologies to create these annotations has become the norm. However, these automated methods commonly result in a small, but significant percentage of genome annotation errors. To improve accuracy and reliability, we analyzed the Caulobacter crescentus NA1000 genome utilizing computer programs Artemis and MICheck to manually examine the third codon position GC content, alignment to a third codon position GC frame plot peak, and matches in the GenBank database. We identified 11 new genes, modified the start site of 113 genes, and changed the reading frame of 38 genes that had been incorrectly annotated. Furthermore, our manual method of identifying protein-coding genes allowed us to remove 112 non-coding regions that had been designated as coding regions. The improved NA1000 genome annotation resulted in a reduction in the use of rare codons since noncoding regions with atypical codon usage were removed from the annotation and 49 new coding regions were added to the annotation. Thus, a more accurate codon usage table was generated as well. These results demonstrate that a comparison of the location of peaks third codon position GC content to the location of protein coding regions could be used to verify the annotation of any genome that has a GC content that is greater than 60%.

A cell cycle and nutritional checkpoint controlling bacterial surface adhesion

In natural environments, bacteria often adhere to surfaces where they form complex multicellular communities. Surface adherence is determined by the biochemical composition of the cell envelope. We describe a novel regulatory mechanism by which the bacterium, Caulobacter crescentus, integrates cell cycle and nutritional signals to control development of an adhesive envelope structure known as the holdfast. Specifically, we have discovered a 68-residue protein inhibitor of holdfast development (HfiA) that directly targets a conserved glycolipid glycosyltransferase required for holdfast production (HfsJ). Multiple cell cycle regulators associate with the hfiA and hfsJ promoters and control their expression, temporally constraining holdfast development to the late stages of G1. HfiA further functions as part of a 'nutritional override' system that decouples holdfast development from the cell cycle in response to nutritional cues. This control mechanism can limit surface adhesion in nutritionally sub-optimal environments without affecting cell cycle progression. We conclude that post-translational regulation of cell envelope enzymes by small proteins like HfiA may provide a general means to modulate the surface properties of bacterial cells.