The transcriptional landscape of a rewritten bacterial genome reveals control elements and genome design principles.

Sequence rewriting enables low-cost genome synthesis and the design of biological systems with orthogonal genetic codes. The error-free, robust rewriting of nucleotide sequences can be achieved with a complete annotation of gene regulatory elements. Here, we compare transcription in Caulobacter crescentus to transcription from plasmid-borne segments of the synthesized genome of C. ethensis 2.0. This rewritten derivative contains an extensive amount of supposedly neutral mutations, including 123'562 synonymous codon changes. The transcriptional landscape refines 60 promoter annotations, exposes 18 termination elements and links extensive transcription throughout the synthesized genome to the unintentional introduction of sigma factor binding motifs. We reveal translational regulation for 20 CDS and uncover an essential translational regulatory element for the expression of ribosomal protein RplS. The annotation of gene regulatory elements allowed us to formulate design principles that improve design schemes for synthesized DNA, en route to a bright future of iteration-free programming of biological systems.

Co-catabolism of arginine and succinate drives symbiotic nitrogen fixation.

Biological nitrogen fixation emerging from the symbiosis between bacteria and crop plants holds promise to increase the sustainability of agriculture. One of the biggest hurdles for the engineering of nitrogen-fixing organisms is an incomplete knowledge of metabolic interactions between microbe and plant. In contrast to the previously assumed supply of only succinate, we describe here the CATCH-N cycle as a novel metabolic pathway that co-catabolizes plant-provided arginine and succinate to drive the energy-demanding process of symbiotic nitrogen fixation in endosymbiotic rhizobia. Using systems biology, isotope labeling studies and transposon sequencing in conjunction with biochemical characterization, we uncovered highly redundant network components of the CATCH-N cycle including transaminases that interlink the co-catabolism of arginine and succinate. The CATCH-N cycle uses N2 as an additional sink for reductant and therefore delivers up to 25% higher yields of nitrogen than classical arginine catabolism-two alanines and three ammonium ions are secreted for each input of arginine and succinate. We argue that the CATCH-N cycle has evolved as part of a synergistic interaction to sustain bacterial metabolism in the microoxic and highly acid environment of symbiosomes. Thus, the CATCH-N cycle entangles the metabolism of both partners to promote symbiosis. Our results provide a theoretical framework and metabolic blueprint for the rational design of plants and plant-associated organisms with new properties to improve nitrogen fixation.

Import of Aspartate and Malate by DcuABC Drives H2/Fumarate Respiration to Promote Initial Salmonella Gut-Lumen Colonization in Mice.

Initial enteropathogen growth in the microbiota-colonized gut is poorly understood. Salmonella Typhimurium is metabolically adaptable and can harvest energy by anaerobic respiration using microbiota-derived hydrogen (H2) as an electron donor and fumarate as an electron acceptor. As fumarate is scarce in the gut, the source of this electron acceptor is unclear. Here, transposon sequencing analysis along the colonization trajectory of S. Typhimurium implicates the C4-dicarboxylate antiporter DcuABC in early murine gut colonization. In competitive colonization assays, DcuABC and enzymes that convert the C4-dicarboxylates aspartate and malate into fumarate (AspA, FumABC), are required for fumarate/H2-dependent initial growth. Thus, S. Typhimurium obtains fumarate by DcuABC-mediated import and conversion of L-malate and L-aspartate. Fumarate reduction yields succinate, which is exported by DcuABC in exchange for L-aspartate and L-malate. This cycle allows S. Typhimurium to harvest energy by H2/fumarate respiration in the microbiota-colonized gut. This strategy may also be relevant for commensal E. coli diminishing the S. Typhimurium infection.

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.

YestroSens, a field-portable S. cerevisiae biosensor device for the detection of endocrine-disrupting chemicals: Reliability and stability.

Farming, industry and urbanization lead to increases in the concentrations of potentially harmful compounds in waste, surface and drinking waters. One example of such pollution are estrogens, the steroidal female reproductive hormones. Already at a few nanograms per litre, these hormones can trigger endocrine disruption and cause acute and chronic health problems in humans and wildlife. Here, we present a Saccharomyces cerevisiae estrogen biosensor capable of detecting estradiol, as well as ethinylestradiol, at concentrations of 1 nM. After an initial characterization of the sensor strain performance in an optimal laboratory setting, we focused on developing a biosensor device. We addressed current limitations of biosensors, such as the requirement of the cells for a liquid growth matrix, controlled storage conditions required to preserve cell viability, and the usually required bulky, as well as expensive, laboratory equipment. Our study provides significant new insights into the field of applied biosensors. The system presented in this work takes microorganism-based analytics one step closer to field application in decentralized locations.

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

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

Identification of Small-Molecule Modulators of Diguanylate Cyclase by FRET-Based High-Throughput Screening.

The bacterial second messenger cyclic diguanosine monophosphate (c-di-GMP) is a key regulator of cellular motility, the cell cycle, and biofilm formation with its resultant antibiotic tolerance, which can make chronic infections difficult to treat. Therefore, diguanylate cyclases, which regulate the spatiotemporal production of c-di-GMP, might be attractive drug targets for control of biofilm formation that is part of chronic infections. We present a FRET-based biochemical high-throughput screening approach coupled with detailed structure-activity studies to identify synthetic small-molecule modulators of the diguanylate cyclase DgcA from Caulobacter crescentus. We identified a set of seven small molecules that regulate DgcA enzymatic activity in the low-micromolar range. Subsequent structure-activity studies on selected scaffolds revealed a remarkable diversity of modulatory behavior, including slight chemical substitutions that reverse the effects from allosteric enzyme inhibition to activation. The compounds identified represent new chemotypes and are potentially developable into chemical genetic tools for the dissection of c-di-GMP signaling networks and alteration of c-di-GMP-associated phenotypes. In sum, our studies underline the importance of detailed mechanism-of-action studies for inhibitors of c-di-GMP signaling and demonstrate the complex interplay between synthetic small molecules and the regulatory mechanisms that control the activity of diguanylate cyclases.

Transposon Sequencing of Brucella abortus Uncovers Essential Genes for Growth In Vitro and Inside Macrophages.

Brucella abortus is a class III zoonotic bacterial pathogen able to survive and replicate inside host cells, including macrophages. Here we report a multidimensional transposon sequencing analysis to identify genes essential for Brucella abortus growth in rich medium and replication in RAW 264.7 macrophages. The construction of a dense transposon mutant library and mapping of 929,769 unique mini-Tn5 insertion sites in the genome allowed identification of 491 essential coding sequences and essential segments in the B. abortus genome. Chromosome II carries a lower proportion (5%) of essential genes than chromosome I (19%), supporting the hypothesis of a recent acquisition of a megaplasmid as the origin of chromosome II. Temporally resolved transposon sequencing analysis as a function of macrophage infection stages identified 79 genes with a specific attenuation phenotype in macrophages, at either 2, 5, or 24 h postinfection, and 86 genes for which the attenuated mutant phenotype correlated with a growth defect on plates. We identified 48 genes required for intracellular growth, including the virB operon, encoding the type IV secretion system, which supports the validity of the screen. The remaining genes encode amino acid and pyrimidine biosynthesis, electron transfer systems, transcriptional regulators, and transporters. In particular, we report the need of an intact pyrimidine nucleotide biosynthesis pathway in order for B. abortus to proliferate inside RAW 264.7 macrophages.

Transposon Sequencing Uncovers an Essential Regulatory Function of Phosphoribulokinase for Methylotrophy.

Methylotrophy is the ability of organisms to grow at the expense of reduced one-carbon compounds, such as methanol or methane. Here, we used transposon sequencing combining hyper-saturated transposon mutagenesis with high-throughput sequencing to define the essential methylotrophy genome of Methylobacterium extorquens PA1, a model methylotroph. To distinguish genomic regions required for growth only on methanol from general required genes, we contrasted growth on methanol with growth on succinate, a non-methylotrophic reference substrate. About 500,000 insertions were mapped for each condition, resulting in a median insertion distance of five base pairs. We identified 147 genes and 76 genes as specific for growth on methanol and succinate, respectively, and a set of 590 genes as required under both growth conditions. For the integration of metabolic functions, we reconstructed a genome-scale metabolic model and performed in silico essentiality analysis. In total, the approach uncovered 95 genes not previously described as crucial for methylotrophy, including genes involved in respiration, carbon metabolism, transport, and regulation. Strikingly, regardless of the absence of the Calvin cycle in the methylotroph, the screen led to the identification of the gene for phosphoribulokinase as essential during growth on methanol, but not during growth on succinate. Genetic experiments in addition to metabolomics and proteomics revealed that phosphoribulokinase serves a key regulatory function. Our data support a model according to which ribulose-1,5-bisphosphate is an essential metabolite that induces a transcriptional regulator driving one-carbon assimilation.

Gene Transfer Agent Promotes Evolvability within the Fittest Subpopulation of a Bacterial Pathogen.

The Bartonella gene transfer agent (BaGTA) is an archetypical example for domestication of a phage-derived element to permit high-frequency genetic exchange in bacterial populations. Here we used multiplexed transposon sequencing (TnSeq) and single-cell reporters to globally define the core components and transfer dynamics of BaGTA. Our systems-level analysis has identified inner- and outer-circle components of the BaGTA system, including 55 regulatory components, as well as an additional 74 and 107 components mediating donor transfer and recipient uptake functions. We show that the stringent response signal guanosine-tetraphosphate (ppGpp) restricts BaGTA induction to a subset of fast-growing cells, whereas BaGTA particle uptake depends on a functional Tol-Pal trans-envelope complex that mediates outer-membrane invagination upon cell division. Our findings suggest that Bartonella evolved an efficient strategy to promote genetic exchange within the fittest subpopulation while disfavoring exchange of deleterious genetic information, thereby facilitating genome integrity and rapid host adaptation.

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

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

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.

Transposon Mutagenesis Paired with Deep Sequencing of Caulobacter crescentus under Uranium Stress Reveals Genes Essential for Detoxification and Stress Tolerance.

UNLABELLED: The ubiquitous aquatic bacterium Caulobacter crescentus is highly resistant to uranium (U) and facilitates U biomineralization and thus holds promise as an agent of U bioremediation. To gain an understanding of how C. crescentus tolerates U, we employed transposon (Tn) mutagenesis paired with deep sequencing (Tn-seq) in a global screen for genomic elements required for U resistance. Of the 3,879 annotated genes in the C. crescentus genome, 37 were found to be specifically associated with fitness under U stress, 15 of which were subsequently tested through mutational analysis. Systematic deletion analysis revealed that mutants lacking outer membrane transporters (rsaFa and rsaFb), a stress-responsive transcription factor (cztR), or a ppGpp synthetase/hydrolase (spoT) exhibited a significantly lower survival rate under U stress. RsaFa and RsaFb, which are homologues of TolC in Escherichia coli, have previously been shown to mediate S-layer export. Transcriptional analysis revealed upregulation of rsaFa and rsaFb by 4- and 10-fold, respectively, in the presence of U. We additionally show that rsaFa mutants accumulated higher levels of U than the wild type, with no significant increase in oxidative stress levels. Our results suggest a function for RsaFa and RsaFb in U efflux and/or maintenance of membrane integrity during U stress. In addition, we present data implicating CztR and SpoT in resistance to U stress. Together, our findings reveal novel gene targets that are key to understanding the molecular mechanisms of U resistance in C. crescentus. IMPORTANCE: Caulobacter crescentus is an aerobic bacterium that is highly resistant to uranium (U) and has great potential to be used in U bioremediation, but its mechanisms of U resistance are poorly understood. We conducted a Tn-seq screen to identify genes specifically required for U resistance in C. crescentus. The genes that we identified have previously remained elusive using other omics approaches and thus provide significant insight into the mechanisms of U resistance by C. crescentus. In particular, we show that outer membrane transporters RsaFa and RsaFb, previously known as part of the S-layer export machinery, may confer U resistance by U efflux and/or by maintaining membrane integrity during U stress.

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

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

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 essential genome of a bacterium.

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

Asymmetrical distribution of the second messenger c-di-GMP upon bacterial cell division.

The bacterial second messenger cyclic diguanosine monophosphate (c-di-GMP) regulates cellular motility and the synthesis of organelles and molecules that promote adhesion to a variety of biological and nonbiological surfaces. These properties likely require tight spatial and temporal regulation of c-di-GMP concentration. We have developed genetically encoded fluorescence resonance energy transfer (FRET)-based biosensors to monitor c-di-GMP concentrations within single bacterial cells by microscopy. Fluctuations of c-di-GMP were visualized in diverse Gram-negative bacterial species and observed to be cell cycle dependent. Asymmetrical distribution of c-di-GMP in the progeny correlated with the time of cell division and polarization for Caulobacter crescentus and Pseudomonas aeruginosa. Thus, asymmetrical distribution of c-di-GMP was observed as part of cell division, which may indicate an important regulatory step in extracellular organelle biosynthesis or function.

High-throughput identification of protein localization dependency networks.

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

Regulation of D-xylose metabolism in Caulobacter crescentus by a LacI-type repressor.

In the oligotrophic freshwater bacterium Caulobacter crescentus, D-xylose induces expression of over 50 genes, including the xyl operon, which encodes key enzymes for xylose metabolism. The promoter (P(xylX)) controlling expression of the xyl operon is widely used as a tool for inducible heterologous gene expression in C. crescentus. We show here that P(xylX) and at least one other promoter in the xylose regulon (P(xylE)) are controlled by the CC3065 (xylR) gene product, a LacI-type repressor. Electrophoretic gel mobility shift assays showed that operator binding by XylR is greatly reduced in the presence of D-xylose. The data support the hypothesis that there is a simple regulatory mechanism in which XylR obstructs xylose-inducible promoters in the absence of the sugar; the repressor is induced to release DNA upon binding D-xylose, thereby freeing the promoter for productive interaction with RNA polymerase. XylR also has an effect on glucose metabolism, as xylR mutants exhibit reduced expression of the Entner-Doudoroff operon and their ability to utilize glucose as a sole carbon and energy source is compromised.