High-throughput identification of protein localization dependency networks

Large-scale identification of genes involved in septal pore plugging in multicellular fungi

Multicellular filamentous fungi have septal pores that allow cytoplasmic exchange, and thus connectivity, between neighboring cells in the filament. Hyphal wounding and other stress conditions induce septal pore closure to minimize cytoplasmic loss. However, the composition of the septal pore and the mechanisms underlying its function are not well understood. Here, we set out to identify new septal components by determining the subcellular localization of 776 uncharacterized proteins in a multicellular ascomycete, Aspergillus oryzae. The set of 776 uncharacterized proteins was selected on the basis that their genes were present in the genomes of multicellular, septal pore-bearing ascomycetes (three Aspergillus species, in subdivision Pezizomycotina) and absent/divergent in the genomes of septal pore-lacking ascomycetes (yeasts). Upon determining their subcellular localization, 62 proteins were found to localize to the septum or septal pore. Deletion of the encoding genes revealed that 23 proteins are involved in regulating septal pore plugging upon hyphal wounding. Thus, this study determines the subcellular localization of many uncharacterized proteins in A. oryzae and, in particular, identifies a set of proteins involved in septal pore function.

Experimental and Analytical Approaches for Improving the Resolution of Randomly Barcoded Transposon Insertion Sequencing (RB-TnSeq) Studies

Randomly barcoded transposon insertion sequencing (RB-TnSeq) is an efficient, multiplexed method to determine microbial gene function during growth under a selection condition of interest. This technique applies to growth, tolerance, and persistence studies in a variety of hosts, but the wealth of data generated can complicate the identification of the most critical gene targets. Experimental and analytical methods for improving the resolution of RB-TnSeq are proposed, using Pseudomonas putida KT2440 as an example organism. Several key parameters, such as baseline media selection, substantially influence the determination of gene fitness. We also present options to increase statistical confidence in gene fitness, including increasing the number of biological replicates and passaging the baseline culture in parallel with selection conditions. These considerations provide practitioners with several options to identify genes of importance in TnSeq data sets, thereby streamlining metabolic characterization.

Protein engineering strategies to stimulate the functions of bacterial pseudokinases

Enzymes orchestrate an array of concerted functions that often culminate in the chemical conversion of substrates into products. In the bacterial kingdom, histidine kinases autophosphorylate, then transfer that phosphate to a second protein called a response regulator. Bacterial genomes can encode large numbers of histidine kinases that provide surveillance of environmental and cytosolic stresses through signal stimulation of histidine kinase activity. Pseudokinases lack these hallmark catalytic functions but often retain binding interactions and allostery. Characterization of bacterial pseudokinases then takes a fundamentally different approach than their enzymatic counterparts. Here we discuss models for how bacterial pseudokinases can utilize protein-protein interactions and allostery to serve as crucial signaling pathway regulators. Then we describe a protein engineering strategy to interrogate these models, emphasizing how signals flow within bacterial pseudokinases. This description includes design considerations, cloning strategies, and the purification of leucine zippers fused to pseudokinases. We then describe two assays to interrogate this approach. First is a C. crescentus swarm plate assay to track motility phenotypes related to a bacterial pseudokinase. Second is an in vitro coupled-enzyme assay that can be applied to test if and how a pseudokinase regulates an active kinase. Together these approaches provide a blueprint for dissecting the mechanisms of cryptic bacterial pseudokinases.

Computational modeling of unphosphorylated CtrA:Cori binding in the Caulobacter cell cycle

In the alphaproteobacterium, Caulobacter crescentus, phosphorylated CtrA (CtrA∼P), a master regulatory protein, binds directly to the chromosome origin (Cori) to inhibit DNA replication. Using a mathematical model of CtrA binding at Cori site [d], we provide computational evidence that CtrA_U can displace CtrA∼P from Cori at the G1-S transition. Investigation of this interaction within a detailed model of the C. crescentus cell cycle suggests that CckA phosphatase may clear Cori of CtrA∼P by altering the [CtrA_U]/[CtrA∼P] ratio rather than by completely depleting CtrA∼P. Model analysis reveals that the mechanism allows for a speedier transition into S phase, stabilizes the timing of chromosome replication under fluctuating rates of CtrA proteolysis, and may contribute to the viability of numerous mutant strains. Overall, these results suggest that CtrA_U enhances the robustness of chromosome replication. More generally, our proposed regulation of CtrA:Cori dynamics may represent a novel motif for molecular signaling in cell physiology.

Dissecting Cellular Function and Distribution of β-Glucosidases in Trichoderma reesei

Trichoderma reesei has 11 putative β-glucosidases in its genome, playing key parts in the induction and production of cellulase. Nevertheless, the reason why the T. reesei genome encodes so many β-glucosidases and the distinct role each β-glucosidase plays in cellulase production remain unknown. In the present study, the cellular function and distribution of 10 known β-glucosidases (CEL3B, CEL3E, CEL3F, CEL3H, CEL3J, CEL1A, CEL3C, CEL1B, CEL3G, and CEL3D) were explored in T. reesei, leaving out BGL1 (CEL3A), which has been well investigated. We found that the overexpression of cel3b or cel3g significantly enhanced extracellular β-glucosidase production, whereas the overexpression of cel1b severely inhibited cellulase production by cellulose, resulting in nearly no growth of T. reesei. Four types of cellular distribution patterns were observed for β-glucosidases in T. reesei: (i) CEL3B, CEL3E, CEL3F, and CEL3G forming clearly separated protein secretion vesicles in the cytoplasm; (ii) CEL3H and CEL3J diffusing the whole endomembrane as well as the cell membrane with protein aggregation, like a reticular network; (iii) CEL1A and CEL3D in vacuoles; (iv) and CEL3C in the nucleus. β-glucosidases CEL1A, CEL3B, CEL3E, CEL3F, CEL3G, CEL3H, and CEL3J were identified as extracellular, CEL3C and CEL3D as intracellular, and CEL1B as unknown. The extracellular β-glucosidases CEL3B, CEL3E, CEL3F, CEL3H, and CEL3G were secreted through a tip-directed conventional secretion pathway, and CEL1A, via a vacuole-mediated pathway that was achieved without any signal peptide, while CEL3J was secreted via an unconventional protein pathway bypassing the endoplasmic reticulum (ER) and Golgi.

The type IV pilin PilA couples surface attachment and cell-cycle initiation in Caulobacter crescentus

Understanding how bacteria colonize surfaces and regulate cell-cycle progression in response to cellular adhesion is of fundamental importance. Here, we use transposon sequencing in conjunction with fluorescence resonance energy transfer (FRET) microscopy to uncover the molecular mechanism for how surface sensing drives cell-cycle initiation in Caulobacter crescentus We identify the type IV pilin protein PilA as the primary signaling input that couples surface contact to cell-cycle initiation via the second messenger cyclic di-GMP (c-di-GMP). Upon retraction of pili filaments, the monomeric pilin reservoir in the inner membrane is sensed by the 17-amino acid transmembrane helix of PilA to activate the PleC-PleD two-component signaling system, increase cellular c-di-GMP levels, and signal the onset of the cell cycle. We termed the PilA signaling sequence CIP for "cell-cycle initiating pilin" peptide. Addition of the chemically synthesized CIP peptide initiates cell-cycle progression and simultaneously inhibits surface attachment. The broad conservation of the type IV pili and their importance in pathogens for host colonization suggests that CIP peptide mimetics offer strategies to inhibit surface sensing, prevent biofilm formation and control persistent infections.

Discovery and Therapeutic Targeting of Differentiated Biofilm Subpopulations

The association of microorganisms into biofilms produces functionally organized microbial structures that promote community survival in a wide range of environments. Much like when individual cells within a multicellular organism express different genes from the same DNA blueprint, individual microbial cells located within different regions of a biofilm structure can exhibit distinct genetic programs. These spatially defined regions of physiologically differentiated cells are reminiscent of the role of tissues in multicellular organisms, with specific subpopulations in the microbial community serving defined roles to promote the overall health of the biofilm. The functions of these subpopulations are quite diverse and can range from dormant cells that can withstand antibiotic onslaughts to cells actively producing extracellular polymeric substances providing integrity to the entire community. The purpose of this review is to discuss the diverse roles of subpopulations in the stability and function of clonal biofilms, the methods for studying these subpopulations, and the ways these subpopulations can potentially be exploited for therapeutic intervention.

Combinatorial Sensor Design in Caulobacter crescentus for Selective Environmental Uranium Detection

The ability to detect uranium (U) through environmental monitoring is of critical importance for informing water resource protection and nonproliferation efforts. While technologies exist for environmental U detection, wide-area environmental monitoring, i.e. sampling coverage over large areas not known to possess U contamination, remains a challenging prospect that necessitates the development of novel detection approaches. Herein, we describe the development of a whole-cell U sensor by integrating two functionally independent, native U-responsive two-component signaling systems (TCS), UzcRS and UrpRS, within an AND gate circuit in the bacterium Caulobacter crescentus. Through leverage of the distinct but imperfect selectivity profiles of both TCS, this combinatorial approach enabled greater selectivity relative to a prior biosensor developed with UzcRS alone; no cross-reactivity was observed with most common environmental metals (e.g, Fe, As, Cu, Ca, Mg, Cd, Cr, Al) or the U decay-chain product Th, and the selectivity against Zn and Pb was significantly improved. In addition, integration of the UzcRS signal amplifier protein UzcY within the AND gate circuit further enhanced overall sensitivity and selectivity for U. The functionality of the sensor in an environmental context was confirmed by detection of U concentrations as low as 1 μM in groundwater samples. The results highlight the value of a combinatorial approach for constructing whole-cell sensors for the selective detection of analytes for which there are no known evolved regulators.

Identification of a U/Zn/Cu responsive global regulatory two-component system in Caulobacter crescentus

Despite the well-known toxicity of uranium (U) to bacteria, little is known about how cells sense and respond to U. The recent finding of a U-specific stress response in Caulobacter crescentus has provided a foundation for studying the mechanisms of U- perception in bacteria. To gain insight into this process, we used a forward genetic screen to identify the regulatory components governing expression of the urcA promoter (P_urcA ) that is strongly induced by U. This approach unearthed a previously uncharacterized two-component system, named UzcRS, which is responsible for U-dependent activation of P_urcA . UzcRS is also highly responsive to zinc and copper, revealing a broader specificity than previously thought. Using ChIP-seq, we found that UzcR binds extensively throughout the genome in a metal-dependent manner and recognizes a noncanonical DNA-binding site. Coupling the genome-wide occupancy data with RNA-seq analysis revealed that UzcR is a global regulator of transcription, predominately activating genes encoding proteins that are localized to the cell envelope; these include metallopeptidases, multidrug-resistant efflux (MDR) pumps, TonB-dependent receptors and many proteins of unknown function. Collectively, our data suggest that UzcRS couples the perception of U, Zn and Cu with a novel extracytoplasmic stress response.

Spatiotemporal Models of the Asymmetric Division Cycle of Caulobacter crescentus

The spatial localization of proteins within the cytoplasm of bacteria is an underappreciated but critical aspect of cell cycle regulation for many prokaryotes. In Caulobacter crescentus-a model organism for the study of asymmetric cell reproduction in prokaryotes-heterogeneous localization of proteins has been identified as the underlying cause of asymmetry in cell morphology, DNA replication, and cell division. However, significant questions remain. Firstly, the mechanisms by which proteins localize in the organelle-free prokaryotic cytoplasm remain obscure. Furthermore, how variations in the spatial and temporal dynamics of cell fate determinants regulate signaling pathways and orchestrate the complex programs of asymmetric cell division and differentiation are subjects of ongoing research. In this chapter, we review current efforts in investigating these two questions. We describe how mathematical models of spatiotemporal protein dynamics are being used to generate and test competing hypotheses and provide complementary insight about the control mechanisms that regulate asymmetry in protein localization and cell division.

SuperSegger: robust image segmentation, analysis and lineage tracking of bacterial cells

Many quantitative cell biology questions require fast yet reliable automated image segmentation to identify and link cells from frame-to-frame, and characterize the cell morphology and fluorescence. We present SuperSegger, an automated MATLAB-based image processing package well-suited to quantitative analysis of high-throughput live-cell fluorescence microscopy of bacterial cells. SuperSegger incorporates machine-learning algorithms to optimize cellular boundaries and automated error resolution to reliably link cells from frame-to-frame. Unlike existing packages, it can reliably segment microcolonies with many cells, facilitating the analysis of cell-cycle dynamics in bacteria as well as cell-contact mediated phenomena. This package has a range of built-in capabilities for characterizing bacterial cells, including the identification of cell division events, mother, daughter and neighbouring cells, and computing statistics on cellular fluorescence, the location and intensity of fluorescent foci. SuperSegger provides a variety of postprocessing data visualization tools for single cell and population level analysis, such as histograms, kymographs, frame mosaics, movies and consensus images. Finally, we demonstrate the power of the package by analyzing lag phase growth with single cell resolution.

MicrobeJ, a tool for high throughput bacterial cell detection and quantitative analysis

Single cell analysis of bacteria and subcellular protein localization dynamics has shown that bacteria have elaborate life cycles, cytoskeletal protein networks, and complex signal transduction pathways driven by localized proteins. The volume of multi-dimensional images generated in such experiments and the computation time required to detect, associate, and track cells and subcellular features pose considerable challenges, especially for high-throughput experiments. Therefore, there is a need for a versatile, computationally efficient image analysis tool capable of extracting the desired relationships from images in a meaningful and unbiased way. Here we present MicrobeJ, a plug-in for the open-source platform ImageJ. MicrobeJ provides a comprehensive framework to process images derived from a wide variety of microscopy experiments with special emphasis on large image sets. It performs the most common intensity and morphology measurements as well as customized detection of poles, septa, fluorescent foci, and organelles, determines their sub-cellular localization with sub-pixel resolution, and tracks them over time. Because a dynamic link is maintained between the images, measurements, and all data representations derived from them, the editor and suite of advanced data presentation tools facilitates the image analysis process and provides a robust way to verify the accuracy and veracity of the data.

New Technologies for Studying Biofilms

Bacteria have traditionally been studied as single-cell organisms. In laboratory settings, aerobic bacteria are usually cultured in aerated flasks, where the cells are considered essentially homogenous. However, in many natural environments, bacteria and other microorganisms grow in mixed communities, often associated with surfaces. Biofilms are comprised of surface-associated microorganisms, their extracellular matrix material, and environmental chemicals that have adsorbed to the bacteria or their matrix material. While this definition of a biofilm is fairly simple, biofilms are complex and dynamic. Our understanding of the activities of individual biofilm cells and whole biofilm systems has developed rapidly, due in part to advances in molecular, analytical, and imaging tools and the miniaturization of tools designed to characterize biofilms at the enzyme level, cellular level, and systems level.

Quantitative Selection Analysis of Bacteriophage φCbK Susceptibility in Caulobacter crescentus

Classical molecular genetics uses stringent selective conditions to identify mutants with distinct phenotypic responses. Mutations giving rise to less pronounced phenotypes are often missed. However, to gain systems-level insights into complex genetic interaction networks requires genome-wide assignment of quantitative phenotypic traits. In this paper, we present a quantitative selection approach coupled with transposon sequencing (QS-TnSeq) to globally identify the cellular components that orchestrate susceptibility of the cell cycle model bacterium Caulobacter crescentus toward bacteriophage φCbK infection. We found that 135 genes representing 3.30% of the Caulobacter genome exhibit significant accumulation of transposon insertions upon φCbK selection. More than 85% thereof consist of new factors not previously associated with phage φCbK susceptibility. Using hierarchical clustering of dose-dependent TnSeq datasets, we grouped these genes into functional modules that correlate with different stages of the φCbK infection process. We assign φCbK susceptibility to eight new genes that represent novel components of the pilus secretion machinery. Further, we demonstrate that, from 86 motility genes, only seven genes encoding structural and regulatory components of the flagellar hook increase phage resistance when disrupted by transposons, suggesting a link between flagellar hook assembly and pili biogenesis. In addition, we observe high recovery of Tn5 insertions within regulatory sequences of the genes encoding the essential NADH:ubiquinone oxidoreductase complex indicating that intact proton motive force is crucial for effective phage propagation. In sum, QS-TnSeq is broadly applicable to perform quantitative and genome-wide systems-genetics analysis of complex phenotypic traits.

Dynamical Localization of DivL and PleC in the Asymmetric Division Cycle of Caulobacter crescentus: A Theoretical Investigation of Alternative Models

Cell-fate asymmetry in the predivisional cell of Caulobacter crescentus requires that the regulatory protein DivL localizes to the new pole of the cell where it up-regulates CckA kinase, resulting in a gradient of CtrA~P across the cell. In the preceding stage of the cell cycle (the "stalked" cell), DivL is localized uniformly along the cell membrane and maintained in an inactive form by DivK~P. It is unclear how DivL overcomes inhibition by DivK~P in the predivisional cell simply by changing its location to the new pole. It has been suggested that co-localization of DivL with PleC phosphatase at the new pole is essential to DivL's activity there. However, there are contrasting views on whether the bifunctional enzyme, PleC, acts as a kinase or phosphatase at the new pole. To explore these ambiguities, we formulated a mathematical model of the spatiotemporal distributions of DivL, PleC and associated proteins (DivJ, DivK, CckA, and CtrA) during the asymmetric division cycle of a Caulobacter cell. By varying localization profiles of DivL and PleC in our model, we show how the physiologically observed spatial distributions of these proteins are essential for the transition from a stalked cell to a predivisional cell. Our simulations suggest that PleC is a kinase in predivisional cells, and that, by sequestering DivK~P, the kinase form of PleC enables DivL to be reactivated at the new pole. Hence, co-localization of PleC kinase and DivL is essential to establishing cellular asymmetry. Our simulations reproduce the experimentally observed spatial distribution and phosphorylation status of CtrA in wild-type and mutant cells. Based on the model, we explore novel combinations of mutant alleles, making predictions that can be tested experimentally.

Effects of (p)ppGpp on the progression of the cell cycle of Caulobacter crescentus

Bacteria must control the progression of their cell cycle in response to nutrient availability. This regulation can be mediated by guanosine tetra- or pentaphosphate [(p)ppGpp], which are synthesized by enzymes of the RelA/SpoT homologue (Rsh) family, particularly under starvation conditions. Here, we study the effects of (p)ppGpp on the cell cycle of Caulobacter crescentus, an oligotrophic bacterium with a dimorphic life cycle. C. crescentus divides asymmetrically, producing a motile swarmer cell that cannot replicate its chromosome and a sessile stalked cell that is replication competent. The swarmer cell rapidly differentiates into a stalked cell in appropriate conditions. An artificial increase in the levels of (p)ppGpp in nonstarved C. crescentus cells was achieved by expressing a truncated relA gene from Escherichia coli, encoding a constitutively active (p)ppGpp synthetase. By combining single-cell microscopy, flow cytometry approaches, and swarming assays, we show that an increase in the intracellular concentration of (p)ppGpp is sufficient to slow down the swarmer-to-stalked cell differentiation process and to delay the initiation of chromosome replication. We also present evidence that the intracellular levels of two master regulators of the cell cycle of C. crescentus, DnaA and CtrA, are modulated in response to (p)ppGpp accumulation, even in the absence of actual starvation. CtrA proteolysis and DnaA synthesis seem indirectly inhibited by (p)ppGpp accumulation. By extending the life span of the motile nonreproductive swarmer cell and thus promoting dispersal and foraging functions over multiplication under starvation conditions, (p)ppGpp may play a central role in the ecological adaptation of C. crescentus to nutritional stresses.

Increasing the Content of High-Content Screening: An Overview

Target-based high-throughput screening (HTS) has recently been critiqued for its relatively poor yield compared to phenotypic screening approaches. One type of phenotypic screening, image-based high-content screening (HCS), has been seen as particularly promising.

In this article, we assess whether HCS is as high content as it can be. We analyze HCS publications and find that although the number of HCS experiments published each year continues to grow steadily, the information content lags behind. We find that a majority of high-content screens published so far (60−80%) made use of only one or two image-based features measured from each sample and disregarded the distribution of those features among each cell population. We discuss several potential explanations, focusing on the hypothesis that data analysis traditions are to blame. This includes practical problems related to managing large and multidimensional HCS data sets as well as the adoption of assay quality statistics from HTS to HCS. Both may have led to the simplification or systematic rejection of assays carrying complex and valuable phenotypic information.

We predict that advanced data analysis methods that enable full multiparametric data to be harvested for entire cell populations will enable HCS to finally reach its potential.

Proteomic profiling of the surface-exposed cell envelope proteins of Caulobacter crescentus

Biotinylation of intact cells, avidin enrichment of derivatized peptides, and shotgun proteomics were employed to reveal the composition of the surface-exposed proteome of the aquatic bacterium, Caulobacter crescentus. Ninety-one unique proteins were identified with the majority originating from the outer membrane, periplasm, and inner membrane, subcellular regions that comprise the Gram-negative bacterium cell envelope. Many of these proteins were described as 'conserved hypothetical protein' or 'hypothetical protein'; and so, the actual expression of these gene products was confirmed. Others did not have any known function or lacked annotation. However, this investigation of the Caulobacter surfaceome did reveal the unanticipated presence of a number of enzymes involved in protein degradation.

BIOLOGICAL SIGNIFICANCE

The results presented here can provide a starting point for hypothesis-driven research projects focused on this bacterium in particular and centered on understanding Gram-negative cell architecture and outer membrane biogenesis broadly. The detected protein degradation enzymes anchored on or located within the outer membrane suggest that Caulobacter has nutrient sources larger than small molecules and/or further processes surface proteins once secreted to this location. Additionally, confirmation of outer membrane residency of those proteins predicted to be periplasmic or whose location prediction was not definitive could potentially elucidate the identities of Gram-negative specific anchorless surface proteins. This article is part of a Special Issue entitled: Trends in Microbial Proteomics.

Cell sorting enriches Escherichia coli mutants that rely on peptidoglycan endopeptidases to suppress highly aberrant morphologies

Bacterial morphology imparts physiological advantages to cells in different environments and, judging by the fidelity with which shape is passed to daughter cells, is a tightly regulated characteristic. Surprisingly, only in the past 10 to 15 years has significant headway been made in identifying the mechanisms by which cells create and maintain particular shapes. One reason for this is that the relevant discoveries have relied heavily on the arduous, somewhat subjective process of manual microscopy. Here, we show that flow cytometry, coupled with the sorting capability of fluorescence-activated cell sorting (FACS), can detect, quantify, and enrich bacteria with morphological alterations. The light scattering properties of several highly aberrant morphological mutants of Escherichia coli were characterized by flow cytometry. Cells from a region that overlapped the distribution of normal rod-shaped cells were collected by FACS and reincubated. After 4 to 15 iterations of this enrichment process, suppressor mutants were isolated that returned almost all the population to a near-normal shape. Suppressors were successfully isolated from strains lacking three or four penicillin binding proteins (PBPs) but not from a mutant lacking a total of seven PBPs. The peptidoglycan endopeptidase, AmpH, was identified as being important for the suppression process, as was a related endopeptidase, MepA. The results validate the use of cell sorting as a means for studying bacterial morphology and identify at least one new class of enzymes required for the suppression of cell shape defects.

Analysis of gene expression levels in individual bacterial cells without image segmentation

Studies of stochasticity in gene expression typically make use of fluorescent protein reporters, which permit the measurement of expression levels within individual cells by fluorescence microscopy. Analysis of such microscopy images is almost invariably based on a segmentation algorithm, where the image of a cell or cluster is analyzed mathematically to delineate individual cell boundaries. However segmentation can be ineffective for studying bacterial cells or clusters, especially at lower magnification, where outlines of individual cells are poorly resolved. Here we demonstrate an alternative method for analyzing such images without segmentation. The method employs a comparison between the pixel brightness in phase contrast vs fluorescence microscopy images. By fitting the correlation between phase contrast and fluorescence intensity to a physical model, we obtain well-defined estimates for the different levels of gene expression that are present in the cell or cluster. The method reveals the boundaries of the individual cells, even if the source images lack the resolution to show these boundaries clearly.

Image analysis in fluorescence microscopy: bacterial dynamics as a case study

Fluorescence microscopy is the primary tool for studying complex processes inside individual living cells. Technical advances in both molecular biology and microscopy have made it possible to image cells from many genetic and environmental backgrounds. These images contain a vast amount of information, which is often hidden behind various sources of noise, convoluted with other information and stochastic in nature. Accessing the desired biological information therefore requires new tools of computational image analysis and modeling. Here, we review some of the recent advances in computational analysis of images obtained from fluorescence microscopy, focusing on bacterial systems. We emphasize techniques that are readily available to molecular and cell biologists but also point out examples where problem-specific image analyses are necessary. Thus, image analysis is not only a toolkit to be applied to new images but also an integral part of the design and implementation of a microscopy experiment.

Structure of the pilus assembly protein TadZ from Eubacterium rectale: implications for polar localization

The tad (tight adherence) locus encodes a protein translocation system that produces a novel variant of type IV pili. The pilus assembly protein TadZ (called CpaE in Caulobacter crescentus) is ubiquitous in tad loci, but is absent in other type IV pilus biogenesis systems. The crystal structure of TadZ from Eubacterium rectale (ErTadZ), in complex with ATP and Mg(2+) , was determined to 2.1 Å resolution. ErTadZ contains an atypical ATPase domain with a variant of a deviant Walker-A motif that retains ATP binding capacity while displaying only low intrinsic ATPase activity. The bound ATP plays an important role in dimerization of ErTadZ. The N-terminal atypical receiver domain resembles the canonical receiver domain of response regulators, but has a degenerate, stripped-down 'active site'. Homology modelling of the N-terminal atypical receiver domain of CpaE indicates that it has a conserved protein-protein binding surface similar to that of the polar localization module of the social mobility protein FrzS, suggesting a similar function. Our structural results also suggest that TadZ localizes to the pole through the atypical receiver domain during an early stage of pili biogenesis, and functions as a hub for recruiting other pili components, thus providing insights into the Tad pilus assembly process.

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

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

Evidence-based annotation of gene function in Shewanella oneidensis MR-1 using genome-wide fitness profiling across 121 conditions

Most genes in bacteria are experimentally uncharacterized and cannot be annotated with a specific function. Given the great diversity of bacteria and the ease of genome sequencing, high-throughput approaches to identify gene function experimentally are needed. Here, we use pools of tagged transposon mutants in the metal-reducing bacterium Shewanella oneidensis MR-1 to probe the mutant fitness of 3,355 genes in 121 diverse conditions including different growth substrates, alternative electron acceptors, stresses, and motility. We find that 2,350 genes have a pattern of fitness that is significantly different from random and 1,230 of these genes (37% of our total assayed genes) have enough signal to show strong biological correlations. We find that genes in all functional categories have phenotypes, including hundreds of hypotheticals, and that potentially redundant genes (over 50% amino acid identity to another gene in the genome) are also likely to have distinct phenotypes. Using fitness patterns, we were able to propose specific molecular functions for 40 genes or operons that lacked specific annotations or had incomplete annotations. In one example, we demonstrate that the previously hypothetical gene SO_3749 encodes a functional acetylornithine deacetylase, thus filling a missing step in S. oneidensis metabolism. Additionally, we demonstrate that the orphan histidine kinase SO_2742 and orphan response regulator SO_2648 form a signal transduction pathway that activates expression of acetyl-CoA synthase and is required for S. oneidensis to grow on acetate as a carbon source. Lastly, we demonstrate that gene expression and mutant fitness are poorly correlated and that mutant fitness generates more confident predictions of gene function than does gene expression. The approach described here can be applied generally to create large-scale gene-phenotype maps for evidence-based annotation of gene function in prokaryotes.

Many computationally predicted gene annotations in bacteria are incomplete or wrong. Consequently, experimental methods to systematically determine gene function in bacteria are required. Here, we describe a genetic approach to meet this challenge. We constructed a large transposon mutant library in the metal-reducing bacterium Shewanella oneidensis MR-1 and profiled the fitness of this collection in more than 100 diverse experimental conditions. In addition to identifying a phenotype for more than 2,000 genes, we demonstrate that mutant fitness profiles can be used to assign “evidence-based” gene annotations for enzymes, signaling proteins, transporters, and transcription factors, a subset of which we verify experimentally.

An SMC ATPase mutant disrupts chromosome segregation in Caulobacter

Accurate replication and segregation of the bacterial genome are essential for cell cycle progression. We have identified a single amino acid substitution in the Caulobacter structural maintenance of chromosomes (SMC) protein that disrupts chromosome segregation and cell division. The E1076Q point mutation in the SMC ATPase domain caused a dominant-negative phenotype in which DNA replication was able to proceed, but duplicated parS centromeres, normally found at opposite cell poles, remained at one pole. The cellular positions of other chromosomal loci were in the wild-type order relative to the parS centromere, but chromosomes remained unsegregated and appeared to be stacked upon one another. Purified SMC-E1076Q was deficient in ATP hydrolysis and exhibited abnormally stable binding to DNA. We propose that SMC spuriously links the duplicated chromosome immediately after passage of the replication fork. In wild-type cells, ATP hydrolysis opens the SMC dimer, freeing one chromosome to segregate to the opposite pole. The loss of ATP hydrolysis causes the SMC-E1076Q dimer to remain bound to both chromosomes, inhibiting segregation.

The essential genome of a bacterium

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Mutations in the nucleotide binding pocket of MreB can alter cell curvature and polar morphology in Caulobacter

The maintenance of cell shape in Caulobacter crescentus requires the essential gene mreB, which encodes a member of the actin superfamily and the target of the antibiotic, A22. We isolated 35 unique A22-resistant Caulobacter strains with single amino acid substitutions near the nucleotide binding site of MreB. Mutations that alter cell curvature and mislocalize the intermediate filament crescentin cluster on the back surface of MreB's structure. Another subset have variable cell widths, with wide cell bodies and actively growing thin extensions of the cell poles that concentrate fluorescent MreB. We found that the extent to which MreB localization is perturbed is linearly correlated with the development of pointed cell poles and variable cell widths. Further, we find that a mutation to glycine of two conserved aspartic acid residues that are important for nucleotide hydrolysis in other members of the actin superfamily abolishes robust midcell recruitment of MreB but supports a normal rate of growth. These mutant strains provide novel insight into how MreB's protein structure, subcellular localization, and activity contribute to its function in bacterial cell shape.

Assembly of the Caulobacter cell division machine

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

Regulatory response to carbon starvation in Caulobacter crescentus

Bacteria adapt to shifts from rapid to slow growth, and have developed strategies for long-term survival during prolonged starvation and stress conditions. We report the regulatory response of C. crescentus to carbon starvation, based on combined high-throughput proteome and transcriptome analyses. Our results identify cell cycle changes in gene expression in response to carbon starvation that involve the prominent role of the FixK FNR/CAP family transcription factor and the CtrA cell cycle regulator. Notably, the SigT ECF sigma factor mediates the carbon starvation-induced degradation of CtrA, while activating a core set of general starvation-stress genes that respond to carbon starvation, osmotic stress, and exposure to heavy metals. Comparison of the response of swarmer cells and stalked cells to carbon starvation revealed four groups of genes that exhibit different expression profiles. Also, cell pole morphogenesis and initiation of chromosome replication normally occurring at the swarmer-to-stalked cell transition are uncoupled in carbon-starved cells.

High-throughput, subpixel precision analysis of bacterial morphogenesis and intracellular spatio-temporal dynamics

Bacteria display various shapes and rely on complex spatial organization of their intracellular components for many cellular processes. This organization changes in response to internal and external cues. Quantitative, unbiased study of these spatio-temporal dynamics requires automated image analysis of large microscopy datasets. We have therefore developed MicrobeTracker, a versatile and high-throughput image analysis program that outlines and segments cells with subpixel precision, even in crowded images and mini-colonies, enabling cell lineage tracking. MicrobeTracker comes with an integrated accessory tool, SpotFinder, which precisely tracks foci of fluorescently labelled molecules inside cells. Using MicrobeTracker, we discover that the dynamics of the extensively studied Escherichia coli Min oscillator depends on Min protein concentration, unveiling critical limitations in robustness within the oscillator. We also find that the fraction of MinD proteins oscillating increases with cell length, indicating that the oscillator has evolved to be most effective when cells attain an appropriate length. MicrobeTracker was also used to uncover novel aspects of morphogenesis and cell cycle regulation in Caulobacter crescentus. By tracking filamentous cells, we show that the chromosomal origin at the old-pole is responsible for most replication/separation events while the others remain largely silent despite contiguous cytoplasm. This surprising position-dependent silencing is regulated by division.

Poles apart: prokaryotic polar organelles and their spatial regulation

While polar organelles hold the key to understanding the fundamentals of cell polarity and cell biological principles in general, they have served in the past merely for taxonomical purposes. Here, we highlight recent efforts in unraveling the molecular basis of polar organelle positioning in bacterial cells. Specifically, we detail the role of members of the Ras-like GTPase superfamily and coiled-coil-rich scaffolding proteins in modulating bacterial cell polarity and in recruiting effector proteins to polar sites. Such roles are well established for eukaryotic cells, but not for bacterial cells that are generally considered diffusion-limited. Studies on spatial regulation of protein positioning in bacterial cells, though still in their infancy, will undoubtedly experience a surge of interest, as comprehensive localization screens have yielded an extensive list of (polarly) localized proteins, potentially reflecting subcellular sites of functional specialization predicted for organelles.

Imaging-based identification of a critical regulator of FtsZ protofilament curvature in Caulobacter

FtsZ is an essential bacterial GTPase that polymerizes at midcell, recruits the division machinery, and may generate constrictive forces necessary for cytokinesis. However, many of the mechanistic details underlying these functions are unknown. We sought to identify FtsZ-binding proteins that influence FtsZ function in Caulobacter crescentus. Here, we present a microscopy-based screen through which we discovered two FtsZ-binding proteins, FzlA and FzlC. FzlA is conserved in α-proteobacteria and was found to be functionally critical for cell division in Caulobacter. FzlA altered FtsZ structure both in vivo and in vitro, forming stable higher-order structures that were resistant to depolymerization by MipZ, a spatial determinant of FtsZ assembly. Electron microscopy revealed that FzlA organizes FtsZ protofilaments into striking helical bundles. The degree of curvature induced by FzlA depended on the nucleotide bound to FtsZ. Induction of FtsZ curvature by FzlA carries implications for regulating FtsZ function by modulating its superstructure.

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

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

Automated quantitative live cell fluorescence microscopy

Advances in microscopy automation and image analysis have given biologists the tools to attempt large scale systems-level experiments on biological systems using microscope image readout. Fluorescence microscopy has become a standard tool for assaying gene function in RNAi knockdown screens and protein localization studies in eukaryotic systems. Similar high throughput studies can be attempted in prokaryotes, though the difficulties surrounding work at the diffraction limit pose challenges, and targeting essential genes in a high throughput way can be difficult. Here we will discuss efforts to make live-cell fluorescent microscopy based experiments using genetically encoded fluorescent reporters an automated, high throughput, and quantitative endeavor amenable to systems-level experiments in bacteria. We emphasize a quantitative data reduction approach, using simulation to help develop biologically relevant cell measurements that completely characterize the cell image. We give an example of how this type of data can be directly exploited by statistical learning algorithms to discover functional pathways.

Polar remodeling and histidine kinase activation, which is essential for Caulobacter cell cycle progression, are dependent on DNA replication initiation

Caulobacter crescentus initiates a single round of DNA replication during each cell cycle. Following the initiation of DNA replication, the essential CckA histidine kinase is activated by phosphorylation, which (via the ChpT phosphotransferase) enables the phosphorylation and activation of the CtrA global regulator. CtrA approximately P then blocks the reinitiation of replication while regulating the transcription of a large number of cell cycle-controlled genes. It has been shown that DNA replication serves as a checkpoint for flagellar biosynthesis and cell division and that this checkpoint is mediated by the availability of active CtrA. Because CckA approximately P promotes the activation of CtrA, we addressed the question of what controls the temporal activation of CckA. We found that the initiation of DNA replication is a prerequisite for remodeling the new cell pole, which includes the localization of the DivL protein kinase to that pole and, consequently, the localization, autophosphorylation, and activation of CckA at that pole. Thus, CckA activation is dependent on polar remodeling and a DNA replication initiation checkpoint that is tightly integrated with the polar phospho-signaling cascade governing cell cycle progression.

Cell pole-specific activation of a critical bacterial cell cycle kinase

Caulobacter crescentus integrates phospho-signaling pathways and transcription factor regulatory cascades to drive the cell cycle. Despite the essential role of the CckA histidine kinase in the control of cell cycle events, the factors that signal its activation at a specific time in the cell cycle have remained elusive. A conditional genetic screen for CckA mislocalization mutants, using automated fluorescence microscopy and an image processing platform, revealed that the essential DivL protein kinase promotes CckA localization, autophosphorylation, and activity at the new cell pole. The transient accumulation of DivL at the new cell pole, but not its kinase activity, is required for the localization and activation of CckA. Because DivL and CckA accumulate at the same cell pole after the initiation of DNA replication and were found to interact in vivo, we propose that DivL recruits CckA to the pole, thereby promoting its autophosphorylation and activity.