High-Resolution Mapping of the Escherichia coli Chromosome Reveals Positions of High and Low Transcription

Spatio-temporal organization of the E. coli chromosome from base to cellular length scales

Escherichia coli has been a vital model organism for studying chromosomal structure, thanks, in part, to its small and circular genome (4.6 million base pairs) and well-characterized biochemical pathways. Over the last several decades, we have made considerable progress in understanding the intricacies of the structure and subsequent function of the E. coli nucleoid. At the smallest scale, DNA, with no physical constraints, takes on a shape reminiscent of a randomly twisted cable, forming mostly random coils but partly affected by its stiffness. This ball-of-spaghetti-like shape forms a structure several times too large to fit into the cell. Once the physiological constraints of the cell are added, the DNA takes on overtwisted (negatively supercoiled) structures, which are shaped by an intricate interplay of many proteins carrying out essential biological processes. At shorter length scales (up to about 1 kb), nucleoid-associated proteins organize and condense the chromosome by inducing loops, bends, and forming bridges. Zooming out further and including cellular processes, topological domains are formed, which are flanked by supercoiling barriers. At the megabase-scale both large, highly self-interacting regions (macrodomains) and strong contacts between distant but co-regulated genes have been observed. At the largest scale, the nucleoid forms a helical ellipsoid. In this review, we will explore the history and recent advances that pave the way for a better understanding of E. coli chromosome organization and structure, discussing the cellular processes that drive changes in DNA shape, and what contributes to compaction and formation of dynamic structures, and in turn how bacterial chromatin affects key processes such as transcription and replication.

Transcription activation in Escherichia coli and Salmonella

Promoter-specific activation of transcript initiation provides an important regulatory device in Escherichia coli and Salmonella. Here, we describe the different mechanisms that operate, focusing on how they have evolved to manage the "housekeeping" bacterial transcription machinery. Some mechanisms involve assisting the bacterial DNA-dependent RNA polymerase or replacing or remodeling one of its subunits. Others are directed to chromosomal DNA, improving promoter function, or relieving repression. We discuss how different activators work together at promoters and how the present complex network of transcription factors evolved.

Dual-Plasmid Mini-Tn5 System to Stably Integrate Multicopy of Target Genes in Escherichia coli

The efficiency of valuable metabolite production by engineered microorganisms underscores the importance of stable and controllable gene expression. While plasmid-based methods offer flexibility, integrating genes into host chromosomes can establish stability without selection pressure. However, achieving site-directed multicopy integration presents challenges, including site selection and stability. We introduced a stable multicopy integration method by using a novel dual-plasmid mini-Tn5 system to insert genes into Escherichia coli's genome. The gene of interest was combined with a removable antibiotic resistance gene. After the selection of bacteria with inserted genes, the antibiotic resistance gene was removed. Optimizations yielded an integration efficiency of approximately 5.5 × 10-3 per recipient cell in a single round. Six rounds of integration resulted in 19 and 5 copies of the egfp gene in the RecA+ strain MG1655 and the RecA- strain XL1-Blue MRF', respectively. Additionally, we integrated a polyhydroxybutyrate (PHB) synthesis gene cluster into E. coli MG1655, yielding an 8-copy integration strain producing more PHB than strains with the cluster on a high-copy plasmid. The method was efficient in generating gene insertions in various E. coli strains, and the inserted genes were stable after extended culture. This stable, high-copy integration tool offers potential for diverse applications in synthetic biology.

Integration Site Library for Efficient Construction of Plasmid-Free Microbial Cell Factories in Escherichia coli

Enhanced production stability and efficiency along with a decrease in production costs are required to build efficient microbial cell factories. Target genes can be integrated into the genome to enhance genetic stability, reduce reliance on antibiotics, and alleviate the metabolic burden. However, selecting the optimal insertion site for the desired gene expression levels remains challenging. Therefore, 18 commonly usedEscherichia coliintegration sites were systematically characterized in this study. Promoters of different strengths were combined with integration sites, yielding a differential intensity range of up to 93-fold. This indicated the versatility and precision of this approach for controlling gene expression levels. Referring to the library, pathway genes were strategically integrated into theE. coligenome based on their respective expression levels. Genetically stable and highly efficient engineered strains that could biosynthesize arbutin and p-aminobenzoic acid were constructed.

2024: A "nucleoid space" odyssey featuring H-NS

The three-dimensional architecture of the bacterial chromosome is intertwined with genome processes such as transcription and replication. Conspicuously so, that the structure of the chromosome permits accurate prediction of active genome processes. Although appreciation of this interplay has developed rapidly in the past two decades, our understanding of this subject is still in its infancy, with research primarily focusing on how the process of transcription regulates and is regulated by chromosome structure. Here, we summarize the latest developments in the field with a focus on the interplay between chromosome structure and transcription in Escherichia coli (E. coli) as mediated by H-NS-a model nucleoid structuring protein. We describe how the organization of chromosomes at the global and local scales is dependent on transcription, and how transcription is regulated by chromosome structure. Finally, we take note of studies that highlight our limited knowledge of structure-function relationships in the chromosome, and we point out research tracks that will improve our insight in the topic.

Plasmid-free production of the plant lignan pinoresinol in growing Escherichia coli cells

BACKGROUND

The high-value aryl tetralin lignan (+)-pinoresinol is the main precursor of many plant lignans including (-)-podophyllotoxin, which is used for the synthesis of chemotherapeutics. As (-)-podophyllotoxin is traditionally isolated from endangered and therefore limited natural sources, there is a particular need for biotechnological production. Recently, we developed a reconstituted biosynthetic pathway from (+)-pinoresinol to (-)-deoxypodophyllotoxin, the direct precursor of (-)-podophyllotoxin, in the recombinant host Escherichia coli. However, the use of the expensive substrate (+)-pinoresinol limits its application from the economic viewpoint. In addition, the simultaneous expression of multiple heterologous genes from different plasmids for a multi-enzyme cascade can be challenging and limits large-scale use.

RESULTS

In this study, recombinant plasmid-free E. coli strains for the multi-step synthesis of pinoresinol from ferulic acid were constructed. To this end, a simple and versatile plasmid toolbox for CRISPR/Cas9-assisted chromosomal integration has been developed, which allows the easy transfer of genes from the pET vector series into the E. coli chromosome. Two versions of the developed toolbox enable the efficient integration of either one or two genes into intergenic high expression loci in both E. coli K-12 and B strains. After evaluation of this toolbox using the fluorescent reporter mCherry, genes from Petroselinum crispum and Zea mays for the synthesis of the monolignol coniferyl alcohol were integrated into different E. coli strains. The product titers achieved with plasmid-free E. coli W3110(T7) were comparable to those of the plasmid-based expression system. For the subsequent oxidative coupling of coniferyl alcohol to pinoresinol, a laccase from Corynebacterium glutamicum was selected. Testing of different culture media as well as optimization of gene copy number and copper availability for laccase activity resulted in the synthesis of 100 mg/L pinoresinol using growing E. coli cells.

CONCLUSIONS

For efficient and simple transfer of genes from pET vectors into the E. coli chromosome, an easy-to-handle molecular toolbox was developed and successfully tested on several E. coli strains. By combining heterologous and endogenous enzymes of the host, a plasmid-free recombinant E. coli growing cell system has been established that enables the synthesis of the key lignan pinoresinol.

CRISPR-COPIES: an in silico platform for discovery of neutral integration sites for CRISPR/Cas-facilitated gene integration

The CRISPR/Cas system has emerged as a powerful tool for genome editing in metabolic engineering and human gene therapy. However, locating the optimal site on the chromosome to integrate heterologous genes using the CRISPR/Cas system remains an open question. Selecting a suitable site for gene integration involves considering multiple complex criteria, including factors related to CRISPR/Cas-mediated integration, genetic stability, and gene expression. Consequently, identifying such sites on specific or different chromosomal locations typically requires extensive characterization efforts. To address these challenges, we have developed CRISPR-COPIES, a COmputational Pipeline for the Identification of CRISPR/Cas-facilitated intEgration Sites. This tool leverages ScaNN, a state-of-the-art model on the embedding-based nearest neighbor search for fast and accurate off-target search, and can identify genome-wide intergenic sites for most bacterial and fungal genomes within minutes. As a proof of concept, we utilized CRISPR-COPIES to characterize neutral integration sites in three diverse species: Saccharomyces cerevisiae, Cupriavidus necator, and HEK293T cells. In addition, we developed a user-friendly web interface for CRISPR-COPIES (https://biofoundry.web.illinois.edu/copies/). We anticipate that CRISPR-COPIES will serve as a valuable tool for targeted DNA integration and aid in the characterization of synthetic biology toolkits, enable rapid strain construction to produce valuable biochemicals, and support human gene and cell therapy applications.

Elucidation of Spatial Positioning of Ribosomes around Chromosome in Escherichia coli Cytoplasm via a Data-Informed Polymer-Based Model

The spatial arrangement of ribosomes and chromosome in Escherichia coli's cytoplasm challenges conventional wisdom. Contrary to the notion of ribosomes acting as inert crowders to the chromosome in the cytoplasm, here we propose a nuanced view by integrating a wide array of experimental data sets into a polymer-based computer model. A set of data-informed computer simulations determines that a delicate balance of attractive and repulsive interactions between ribosomes and the chromosome is required in order to reproduce experimentally obtained linear densities and brings forth the view that ribosomes are not mere inert crowders in the cytoplasm. The model finds that the ribosomes represent themselves as a poor solvent for the chromosome with a 50 nm mesh size, consistent with previous experimental analysis. Our multidimensional analysis of ribosome distribution, both free (30S and 50S) and bound (70S polysome), uncovers a relatively less pronounced segregation pattern than previously thought. Notably, we identify a ribosome-rich central region within the innermost core of the nucleoid. Moreover, our exploration of the chromosome mesh size and the conformation of bound ribosomes suggests that these ribosomes maintain elongated shapes, enabling them to navigate through the chromosome mesh and access the central core. This dynamic localization challenges the static segregation model and underscores the pivotal role of ribosome-chromosome interactions in cellular media.

Development of a Data-Driven Integrative Model of a Bacterial Chromosome

The chromosome of archetypal bacteria E. coli is known for a complex topology with a 4.6 × 106 base pairs (bp) long sequence of nucleotides packed within a micrometer-sized cellular confinement. The inherent organization underlying this chromosome eludes general consensus due to the lack of a high-resolution picture of its conformation. Here we present our development of an integrative model of E. coli at a 500 bp resolution (https://github.com/JMLab-tifrh/ecoli_finer), which optimally combines a set of multiresolution genome-wide experimentally measured data within a framework of polymer based architecture. In particular the model is informed with an intragenome contact probability map at 5000 bp resolution derived via the Hi-C experiment and RNA-sequencing data at 500 bp resolution. Via dynamical simulations, this data-driven polymer based model generates an appropriate conformational ensemble commensurate with chromosome architectures that E. coli adopts. As a key hallmark of the E. coli chromosome the model spontaneously self-organizes into a set of nonoverlapping macrodomains and suitably locates plectonemic loops near the cell membrane. As novel extensions, it predicts a contact probability map simulated at a higher resolution than precedent experiments and can demonstrate segregation of chromosomes in a partially replicating cell. Finally, the modular nature of the model helps us devise control simulations to quantify the individual role of key features in hierarchical organization of the bacterial chromosome.

Standardized Iterative Genome Editing Method for Escherichia coli Based on CRISPR-Cas9

The introduction of complex biosynthetic pathways into the hosts' chromosomes is gaining attention with the development of synthetic biology. While CRISPR-Cas9 has been widely employed for gene knock-in, the process of multigene insertion remains cumbersome due to laborious and empirical gene cloning procedures. To address this, we devised a standardized iterative genome editing system for Escherichia coli, harnessing the power of CRISPR-Cas9 and MetClo assembly. This comprehensive toolkit comprises two fundamental elements based on the Golden Gate standard for modular assembly of sgRNA or CRISPR arrays and donor DNAs. We achieved a gene insertion efficiency of up to 100%, targeting a single locus. Expression of tracrRNA using a strong promoter enhances multiplex genomic insertion efficiency to 7.3%, compared with 0.76% when a native promoter is used. To demonstrate the robust capabilities of this genome editing toolbox, we successfully integrated 5-10 genes from the coenzyme B12 biosynthetic pathway ranging from 5.3 to 8 Kb in length into the chromosome of E. coli chassis cells, resulting in 14 antibiotic-free, plasmid-free producers. Following an extensive screening process involving genes from diverse sources, cistronic design modifications, and chromosome repositioning, we obtained a recombinant strain yielding 1.49 mg L-1 coenzyme B12, the highest known titer achieved by using E. coli as the producer. Illuminating its user-friendliness, this genome editing system is an exceedingly versatile tool for expediently integrating complex biosynthetic pathway genes into hosts' genomes, thus facilitating pathway optimization for chemical production.

Enabling pathway design by multiplex experimentation and machine learning

The remarkable metabolic diversity observed in nature has provided a foundation for sustainable production of a wide array of valuable molecules. However, transferring the biosynthetic pathway to the desired host often runs into inherent failures that arise from intermediate accumulation and reduced flux resulting from competing pathways within the host cell. Moreover, the conventional trial and error methods utilized in pathway optimization struggle to fully grasp the intricacies of installed pathways, leading to time-consuming and labor-intensive experiments, ultimately resulting in suboptimal yields. Considering these obstacles, there is a pressing need to explore the enzyme expression landscape and identify the optimal pathway configuration for enhanced production of molecules. This review delves into recent advancements in pathway engineering, with a focus on multiplex experimentation and machine learning techniques. These approaches play a pivotal role in overcoming the limitations of traditional methods, enabling exploration of a broader design space and increasing the likelihood of discovering optimal pathway configurations for enhanced production of molecules. We discuss several tools and strategies for pathway design, construction, and optimization for sustainable and cost-effective microbial production of molecules ranging from bulk to fine chemicals. We also highlight major successes in academia and industry through compelling case studies.

Massive integration of large gene libraries in the chromosome of Escherichia coli

Large gene libraries are frequently created in Escherichia coli plasmids, which can induce cell toxicity and expression instability due to the high gene dosage. To address these limitations, gene libraries can be integrated in a single copy into the bacterial chromosome. Here, we describe an efficient system for the massive integration (MAIN) of large gene libraries in the E. coli chromosome that generates in-frame gene fusions that are expressed stably. MAIN uses a thermosensitive integrative plasmid that is linearized in vivo to promote extensive integration of the gene library via homologous recombination. Positive and negative selections efficiently remove bacteria lacking gene integration in the target site. We tested MAIN with a library of 107 VHH genes that encode nanobodies (Nbs). The integration of VHH genes into a custom target locus of the E. coli chromosome enabled stable expression and surface display of the Nbs. Next-generation DNA sequencing confirmed that MAIN preserved the diversity of the gene library after integration. Finally, we screened the integrated library to select Nbs that bind a specific antigen using magnetic and fluorescence-activated cell sorting. This allowed us to identify Nbs binding the epidermal growth factor receptor that were not previously isolated in a similar screening of a multicopy plasmid library. Our results demonstrate that MAIN enables large gene library integration into the E. coli chromosome, creating stably expressed in-frame fusions for functional screening.

The roles of nucleoid-associated proteins and topoisomerases in chromosome structure, strand segregation, and the generation of phenotypic heterogeneity in bacteria

How to adapt to a changing environment is a fundamental, recurrent problem confronting cells. One solution is for cells to organize their constituents into a limited number of spatially extended, functionally relevant, macromolecular assemblies or hyperstructures, and then to segregate these hyperstructures asymmetrically into daughter cells. This asymmetric segregation becomes a particularly powerful way of generating a coherent phenotypic diversity when the segregation of certain hyperstructures is with only one of the parental DNA strands and when this pattern of segregation continues over successive generations. Candidate hyperstructures for such asymmetric segregation in prokaryotes include those containing the nucleoid-associated proteins (NAPs) and the topoisomerases. Another solution to the problem of creating a coherent phenotypic diversity is by creating a growth-environment-dependent gradient of supercoiling generated along the replication origin-to-terminus axis of the bacterial chromosome. This gradient is modulated by transcription, NAPs, and topoisomerases. Here, we focus primarily on two topoisomerases, TopoIV and DNA gyrase in Escherichia coli, on three of its NAPs (H-NS, HU, and IHF), and on the single-stranded binding protein, SSB. We propose that the combination of supercoiling-gradient-dependent and strand-segregation-dependent topoisomerase activities result in significant differences in the supercoiling of daughter chromosomes, and hence in the phenotypes of daughter cells.

In Vivo Sampling of Intracellular Heterogeneity of Pseudomonas putida Enables Multiobjective Optimization of Genetic Devices

The inner physicochemical heterogeneity of bacterial cells generates three-dimensional (3D)-dependent variations of resources for effective expression of given chromosomally located genes. This fact has been exploited for adjusting the most favorable parameters for implanting a complex device for optogenetic control of biofilm formation in the soil bacterium Pseudomonas putida. To this end, a DNA segment encoding a superactive variant of the Caulobacter crescendus diguanylate cyclase PleD expressed under the control of the cyanobacterial light-responsive CcaSR system was placed in a mini-Tn5 transposon vector and randomly inserted through the chromosome of wild-type and biofilm-deficient variants of P. putida lacking the wsp gene cluster. This operation delivered a collection of clones covering a whole range of biofilm-building capacities and dynamic ranges in response to green light. Since the phenotypic output of the device depends on a large number of parameters (multiple promoters, RNA stability, translational efficacy, metabolic precursors, protein folding, etc.), we argue that random chromosomal insertions enable sampling the intracellular milieu for an optimal set of resources that deliver a preset phenotypic specification. Results thus support the notion that the context dependency can be exploited as a tool for multiobjective optimization, rather than a foe to be suppressed in Synthetic Biology constructs.

Bacterial Subcellular Architecture, Structural Epistasis, and Antibiotic Resistance

Epistasis refers to the way in which genetic interactions between some genetic loci affect phenotypes and fitness. In this study, we propose the concept of "structural epistasis" to emphasize the role of the variable physical interactions between molecules located in particular spaces inside the bacterial cell in the emergence of novel phenotypes. The architecture of the bacterial cell (typically Gram-negative), which consists of concentrical layers of membranes, particles, and molecules with differing configurations and densities (from the outer membrane to the nucleoid) determines and is in turn determined by the cell shape and size, depending on the growth phases, exposure to toxic conditions, stress responses, and the bacterial environment. Antibiotics change the bacterial cell's internal molecular topology, producing unexpected interactions among molecules. In contrast, changes in shape and size may alter antibiotic action. The mechanisms of antibiotic resistance (and their vectors, as mobile genetic elements) also influence molecular connectivity in the bacterial cell and can produce unexpected phenotypes, influencing the action of other antimicrobial agents.

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

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

Systematic metabolic engineering of Escherichia coli for the enhanced production of cinnamaldehyde

Cinnamaldehyde (CAD) derived from cinnamon bark has received much attention for its potential as a nematicide and food additive. Previously, we have succeeded in developing an Escherichia coli strain (YHP05) capable of synthesizing cinnamaldehyde; however, the production titer (75 mg/L) was not sufficient for commercialization. Herein, to develop an economical and sustainable production bioprocess, we further engineered the YHP05 strain for non-auxotrophic, antibiotic-free, inducer-free hyperproduction of CAD using systematic metabolic engineering. First, the conversion of trans-cinnamic acid (t-CA) to CAD was improved by the co-expression of carboxylic acid reductase and phosphopantetheinyl transferase (PPTase) genes. Second, to prevent the spontaneous conversion of CAD to cinnamyl alcohol, 10 endogenous reductase and dehydrogenase genes were deleted. Third, all expression cassettes were integrated into the chromosomal DNA using an auto-inducible system for antibiotic- and inducer-free production. Subsequently, to facilitate CAD production, available pools of cofactors (NADPH, CoA, and ATP) were increased, and acetate pathways were deleted. With the final antibiotic-, plasmid-, and inducer-free strain (H-11MPmR), fed-batch cultivations combined with in situ product recovery (ISPR) were performed, and the production titer of CAD as high as 3.8 g/L could be achieved with 49.1 mg/L/h productivity, which is the highest CAD titer ever reported.

Strand specificity of ribonucleotide excision repair in Escherichia coli

In Escherichia coli, replication of both strands of genomic DNA is carried out by a single replicase-DNA polymerase III holoenzyme (pol III HE). However, in certain genetic backgrounds, the low-fidelity TLS polymerase, DNA polymerase V (pol V) gains access to undamaged genomic DNA where it promotes elevated levels of spontaneous mutagenesis preferentially on the lagging strand. We employed active site mutants of pol III (pol IIIα_S759N) and pol V (pol V_Y11A) to analyze ribonucleotide incorporation and removal from the E. coli chromosome on a genome-wide scale under conditions of normal replication, as well as SOS induction. Using a variety of methods tuned to the specific properties of these polymerases (analysis of lacI mutational spectra, lacZ reversion assay, HydEn-seq, alkaline gel electrophoresis), we present evidence that repair of ribonucleotides from both DNA strands in E. coli is unequal. While RNase HII plays a primary role in leading-strand Ribonucleotide Excision Repair (RER), the lagging strand is subject to other repair systems (RNase HI and under conditions of SOS activation also Nucleotide Excision Repair). Importantly, we suggest that RNase HI activity can also influence the repair of single ribonucleotides incorporated by the replicase pol III HE into the lagging strand.

Effective 5-aminolevulinic acid production via T7 RNA polymerase and RuBisCO equipped Escherichia coli W3110

Chromosome-based engineering is a superior approach for gene integration generating a stable and robust chassis. Therefore, an effective amplifier, T7 RNA polymerase (T7RNAP) from bacteriophage, has been incorporated into Escherichia coli W3110 by site-specific integration. Herein, we performed the 5-aminolevulinic acid (5-ALA) production in four T7RNAP-equipped W3110 strains using recombinant 5-aminolevulinic synthase and further explored the metabolic difference in best strain. The fastest glucose consumption resulted in the highest biomass and the 5-ALA production reached to 5.5 g/L; thus, the least by-product of acetate was shown in RH strain in which T7RNAP was inserted at HK022 phage attack site. Overexpression of phosphoenolpyruvate (PEP) carboxylase would pull PEP to oxaloacetic acid in tricarboxylic acid cycle, leading to energy conservation and even no acetate production, thus, 6.53 g/L of 5-ALA was achieved. Amino acid utilization in RH deciphered the major metabolic flux in α-ketoglutaric acid dominating 5-ALA production. Finally, the ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and phosphoribulokinase were expressed for carbon dioxide recycling; a robust and efficient chassis toward low-carbon assimilation and high-level of 5-ALA production up to 11.2 g/L in fed-batch fermentation was established.

The role of replication-induced chromosomal copy numbers in spatio-temporal gene regulation and evolutionary chromosome plasticity

For a coherent response to environmental changes, bacterial evolution has formed a complex transcriptional regulatory system comprising classical DNA binding proteins sigma factors and modulation of DNA topology. In this study, we investigate replication-induced gene copy numbers - a regulatory concept that is unlike the others not based on modulation of promoter activity but on replication dynamics. We show that a large fraction of genes are predominantly affected by transient copy numbers and identify cellular functions and central pathways governed by this mechanism in Escherichia coli. Furthermore, we show quantitatively that the previously observed spatio-temporal expression pattern between different growth phases mainly emerges from transient chromosomal copy numbers. We extend the analysis to the plant pathogen Dickeya dadantii and the biotechnologically relevant organism Vibrio natriegens. The analysis reveals a connection between growth phase dependent gene expression and evolutionary gene migration in these species. A further extension to the bacterial kingdom indicates that chromosome evolution is governed by growth rate related transient copy numbers.

Ribosomal RNA operons define a central functional compartment in the Streptomyces chromosome

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. These bacteria possess a large linear chromosome genetically compartmentalized: core genes are grouped in the central part, while terminal regions are populated by poorly conserved genes. In exponentially growing cells, chromosome conformation capture unveiled sharp boundaries formed by ribosomal RNA (rrn) operons that segment the chromosome into multiple domains. Here we further explore the link between the genetic distribution of rrn operons and Streptomyces genetic compartmentalization. A large panel of genomes of species representative of the genus diversity revealed that rrn operons and core genes form a central skeleton, the former being identifiable from their core gene environment. We implemented a new nomenclature for Streptomyces genomes and trace their rrn-based evolutionary history. Remarkably, rrn operons are close to pericentric inversions. Moreover, the central compartment delimited by rrn operons has a very dense, nearly invariant core gene content. Finally, this compartment harbors genes with the highest expression levels, regardless of gene persistence and distance to the origin of replication. Our results highlight that rrn operons are structural boundaries of a central functional compartment prone to transcription in Streptomyces.

A conundrum in 3-D genome organization and expression?

Recent advances in our understanding of how the genome is folded within the nucleus have included cases in which this positioning correlates with gene expression, either positively or negatively. But is the 3-D location of a gene a cause or an effect of its expression? In this Perspective I articulate the problem and then cite as guideposts recent cases where causation has indeed been arguably established. The hope is to critically illuminate this issue for continued consideration in this important, evolving field.

Efficient Bacterial Genome Engineering throughout the Central Dogma Using the Dual-Selection Marker tetAOPT

Engineering of bacterial genomes is a fundamental craft in contemporary biotechnology. The ability to precisely edit chromosomes allows for the development of cells with specific phenotypes for metabolic engineering and for the creation of minimized genomes. Genetic tools are needed to select for cells that underwent editing, and dual-selection markers that enable both positive and negative selection are highly useful. Here, we present an optimized and easy-to-use version of the tetA dual-selection marker and demonstrate how this tetAOPT can be used efficiently to engineer at different stages of the central dogma of molecular biology. On the DNA level, tetAOPT can be used to create scarless knockouts across the Escherichia coli genome with efficiency above 90%, whereas recombinant gene integrations can be achieved with approximately 50% efficiency. On the RNA and protein level, we show that tetAOPT enables advanced genome engineering of both gene translation and transcription by introducing sequence variation in the translation initiation region or by exchanging promoters. Finally, we demonstrate the use of tetAOPT for genome engineering in the industrially relevant probiotic strain E. coli Nissle.

Genetic context effects can override canonical cis regulatory elements in Escherichia coli

Recent experiments have shown that in addition to control by cis regulatory elements, the local chromosomal context of a gene also has a profound impact on its transcription. Although this chromosome-position dependent expression variation has been empirically mapped at high-resolution, the underlying causes of the variation have not been elucidated. Here, we demonstrate that 1 kb of flanking, non-coding synthetic sequences with a low frequency of guanosine and cytosine (GC) can dramatically reduce reporter expression compared to neutral and high GC-content flanks in Escherichia coli. Natural and artificial genetic context can have a similarly strong effect on reporter expression, regardless of cell growth phase or medium. Despite the strong reduction in the maximal expression level from the fully-induced reporter, low GC synthetic flanks do not affect the time required to reach the maximal expression level after induction. Overall, we demonstrate key determinants of transcriptional propensity that appear to act as tunable modulators of transcription, independent of regulatory sequences such as the promoter. These findings provide insight into the regulation of naturally occurring genes and an independent control for optimizing expression of synthetic biology constructs.

A spatially resolved stochastic model reveals the role of supercoiling in transcription regulation

In Escherichia coli, translocation of RNA polymerase (RNAP) during transcription introduces supercoiling to DNA, which influences the initiation and elongation behaviors of RNAP. To quantify the role of supercoiling in transcription regulation, we developed a spatially resolved supercoiling model of transcription. The integrated model describes how RNAP activity feeds back with the local DNA supercoiling and how this mechanochemical feedback controls transcription, subject to topoisomerase activities and stochastic topological domain formation. This model establishes that transcription-induced supercoiling mediates the cooperation of co-transcribing RNAP molecules in highly expressed genes, and this cooperation is achieved under moderate supercoiling diffusion and high topoisomerase unbinding rates. It predicts that a topological domain could serve as a transcription regulator, generating substantial transcriptional noise. It also shows the relative orientation of two closely arranged genes plays an important role in regulating their transcription. The model provides a quantitative platform for investigating how genome organization impacts transcription.

The context matrix: Navigating biological complexity for advanced biodesign

Synthetic biology offers many solutions in healthcare, production, sensing and agriculture. However, the ability to rationally engineer synthetic biosystems with predictable and robust functionality remains a challenge. A major reason is the complex interplay between the synthetic genetic construct, its host, and the environment. Each of these contexts contains a number of input factors which together can create unpredictable behaviours in the engineered biosystem. It has become apparent that for the accurate assessment of these contextual effects a more holistic approach to design and characterisation is required. In this perspective article, we present the context matrix, a conceptual framework to categorise and explore these contexts and their net effect on the designed synthetic biosystem. We propose the use and community-development of the context matrix as an aid for experimental design that simplifies navigation through the complex design space in synthetic biology.

Escherichia coli transcription factors of unknown function: sequence features and possible evolutionary relationships

Organisms need mechanisms to perceive the environment and respond accordingly to environmental changes or the presence of hazards. Transcription factors (TFs) are required for cells to respond to the environment by controlling the expression of genes needed. Escherichia coli has been the model bacterium for many decades, and still, there are features embedded in its genome that remain unstudied. To date, 58 TFs remain poorly characterized, although their binding sites have been experimentally determined. This study showed that these TFs have sequence variation at the third codon position G+C content but maintain the same Codon Adaptation Index (CAI) trend as annotated functional transcription factors. Most of these transcription factors are in areas of the genome where abundant repetitive and mobile elements are present. Sequence divergence points to groups with distinctive sequence signatures but maintaining the same type of DNA binding domain. Finally, the analysis of the promoter sequences of the 58 TFs showed A+T rich regions that agree with the features of horizontally transferred genes. The findings reported here pave the way for future research of these TFs that may uncover their role as spare factors in case of lose-of-function mutations in core TFs and trace back their evolutionary history.

Relationship between the Chromosome Structural Dynamics and Gene Expression—A Chicken and Egg Dilemma?

Prokaryotic transcription was extensively studied over the last half-century. A great deal of data has been accumulated regarding the control of gene expression by transcription factors regulating their target genes by binding at specific DNA sites. However, there is a significant gap between the mechanistic description of transcriptional control obtained from in vitro biochemical studies and the complexity of transcriptional regulation in the context of the living cell. Indeed, recent studies provide ample evidence for additional levels of complexity pertaining to the regulation of transcription in vivo, such as, for example, the role of the subcellular localization and spatial organization of different molecular components involved in the transcriptional control and, especially, the role of chromosome configurational dynamics. The question as to how the chromosome is dynamically reorganized under the changing environmental conditions and how this reorganization is related to gene expression is still far from being clear. In this article, we focus on the relationships between the chromosome structural dynamics and modulation of gene expression during bacterial adaptation. We argue that spatial organization of the bacterial chromosome is of central importance in the adaptation of gene expression to changing environmental conditions and vice versa, that gene expression affects chromosome dynamics.

The effect of droplet size on syntrophic dynamics in droplet-enabled microbial co-cultivation

Co-cultivation in microfluidic droplets has emerged as a versatile tool for the study of natural and synthetic microbial communities. In particular, the identification and characterization of syntrophic interactions in these communities is attracting increasing interest due to their critical importance for the functioning of environmental and host-associated communities as well as new biotechnological applications. However, one critical parameter in droplet-enabled co-cultivation that has evaded appropriate evaluation is the droplet size. Given the same number of initial cells, a larger droplet size can increase the length scale secreted metabolites must diffuse as well as dilute the initial concentration of cells and exchanged metabolites, impacting the community dynamics. To evaluate the effect of droplet size on a spectrum of syntrophic interactions, we cultivated a synthetic model system consisting of two E. coli auxotrophs, whose interactions could be modulated through supplementation of related amino acids in the medium. Our results demonstrate that the droplet size impacts substantially numerous aspects of the growth of a cross-feeding bi-culture, particularly the growth capacity, maximum specific growth rate, and lag time, depending on the degree of the interaction. This work heavily suggests that one droplet size does not fit all types of interactions; this parameter should be carefully evaluated and chosen in experimental studies that aim to utilize droplet-enabled co-cultivation to characterize or elucidate microbial interactions.

The zinc-finger bearing xenogeneic silencer MucR in α-proteobacteria balances adaptation and regulatory integrity

Foreign AT-rich genes drive bacterial adaptation to new niches while challenging the existing regulation network. Here we report that MucR, a conserved regulator in α-proteobacteria, balances adaptation and regulatory integrity in Sinorhizobium fredii, a facultative microsymbiont of legumes. Chromatin immunoprecipitation sequencing coupled with transcriptomic data reveal that average transcription levels of both target and non-target genes, under free-living and symbiotic conditions, increase with their conservation levels. Targets involved in environmental adaptation and symbiosis belong to genus or species core and can be repressed or activated by MucR in a condition-dependent manner, implying regulatory integrations. However, most targets are enriched in strain-specific genes of lower expression levels and higher AT%. Within each conservation levels, targets have higher AT% and average transcription levels than non-target genes and can be further up-regulated in the mucR mutant. This is consistent with higher AT% of spacers between -35 and -10 elements of promoters for target genes, which enhances transcription. The MucR recruitment level linearly increases with AT% and the number of a flexible pattern (with periodic repeats of Ts) of target sequences. Collectively, MucR directly represses AT-rich foreign genes with predisposed high transcription potential while progressive erosions of its target sites facilitate regulatory integrations of foreign genes.

CRISPRi enables fast growth followed by stable aerobic pyruvate formation in Escherichia coli without auxotrophy

CRISPR interference (CRISPRi) was applied to enable the aerobic production of pyruvate in Escherichia coli MG1655 under glucose excess conditions by targeting the promoter regions of aceE or pdhR. Knockdown strains were cultivated in aerobic shaking flasks and the influence of inducer concentration and different sgRNA binding sites on the production of pyruvate was measured. Targeting the promoter regions of aceE or pdhR triggered pyruvate production during the exponential phase and reduced expression of aceE. In lab-scale bioreactor fermentations, an aceE silenced strain successfully produced pyruvate under fully aerobic conditions during the exponential phase, but loss of productivity occurred during a subsequent nitrogen-limited phase. Targeting the promoter region of pdhR enabled pyruvate production during the growth phase of cultivations, and a continued low-level accumulation during the nitrogen-limited production phase. Combinatorial targeting of the promoter regions of both aceE and pdhR in E. coli MG1655 pdCas9 psgRNA_aceE_234_pdhR_329 resulted in the stable aerobic production of pyruvate with non-growing cells at YP/S  =  0.36 ± 0.029 gPyruvate/gGlucose in lab-scale bioreactors throughout an extended nitrogen-limited production phase.

Distinct heterochromatin-like domains promote transcriptional memory and silence parasitic genetic elements in bacteria

There is increasing evidence that prokaryotes maintain chromosome structure, which in turn impacts gene expression. We recently characterized densely occupied, multi-kilobase regions in the E. coli genome that are transcriptionally silent, similar to eukaryotic heterochromatin. These extended protein occupancy domains (EPODs) span genomic regions containing genes encoding metabolic pathways as well as parasitic elements such as prophages. Here, we investigate the contributions of nucleoid-associated proteins (NAPs) to the structuring of these domains, by examining the impacts of deleting NAPs on EPODs genome-wide in E. coli and B. subtilis. We identify key NAPs contributing to the silencing of specific EPODs, whose deletion opens a chromosomal region for RNA polymerase binding at genes contained within that region. We show that changes in E. coli EPODs facilitate an extra layer of transcriptional regulation, which prepares cells for exposure to exotic carbon sources. Furthermore, we distinguish novel xenogeneic silencing roles for the NAPs Fis and Hfq, with the presence of at least one being essential for cell viability in the presence of domesticated prophages. Our findings reveal previously unrecognized mechanisms through which genomic architecture primes bacteria for changing metabolic environments and silences harmful genomic elements.

Genomic inversions in Escherichia coli alter gene expression and are associated with nucleoid protein binding sites

Genomic reorganization, like rearrangements and inversions, influence how genetic information is organized within bacterial genomes. Inversions in particular, facilitate genome evolution through gene gain and loss, and can alter gene expression. Previous studies investigating the impact inversions have on gene expression induced inversions targeting specific genes or examine inversions between distantly related species. This fails to encompass a genome wide perspective on naturally occurring inversions and their post adaptation impact on gene expression. Here we use bioinformatic techniques and multiple RNA-seq datasets to investigate the short- and long-range impact inversions have on genomic gene expression within <i>Escherichia coli</i>. We observed differences in gene expression between homologous inverted and non-inverted genes, even after long term exposure to adaptive selection. In 4% of inversions representing 33 genes, differential gene expression between inverted and non-inverted homologs was detected, with nearly two thirds (71%) of differentially expressed inverted genes having 9.4-85.6 fold higher gene expression. The identified inversions had more overlap than expected with nucleoid associated protein binding sites, which assist in genomic gene expression regulation. Some inversions can drastically impact gene expression even between different strains of <i>E.coli</i>, and could provide a mechanism for the diversification of genetic content through controlled expression changes.

The Dynamics of Mycoplasma gallisepticum Nucleoid Structure at the Exponential and Stationary Growth Phases

The structure and dynamics of bacterial nucleoids play important roles in regulating gene expression. Bacteria of class Mollicutes and, in particular, mycoplasmas feature extremely reduced genomes. They lack multiple structural proteins of the nucleoid, as well as regulators of gene expression. We studied the organization of Mycoplasma gallisepticum nucleoids in the stationary and exponential growth phases at the structural and protein levels. The growth phase transition results in the structural reorganization of M. gallisepticum nucleoid. In particular, it undergoes condensation and changes in the protein content. The observed changes corroborate with the previously identified global rearrangement of the transcriptional landscape in this bacterium during the growth phase transition. In addition, we identified that the glycolytic enzyme enolase functions as a nucleoid structural protein in this bacterium. It is capable of non-specific DNA binding and can form fibril-like complexes with DNA.

Intelligent host engineering for metabolic flux optimisation in biotechnology

Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.

A standardized genome architecture for bacterial synthetic biology (SEGA)

Chromosomal recombinant gene expression offers a number of advantages over plasmid-based synthetic biology. However, the methods applied for bacterial genome engineering are still challenging and far from being standardized. Here, in an attempt to realize the simplest recombinant genome technology imaginable and facilitate the transition from recombinant plasmids to genomes, we create a simplistic methodology and a comprehensive strain collection called the Standardized Genome Architecture (SEGA). In its simplest form, SEGA enables genome engineering by combining only two reagents: a DNA fragment that can be ordered from a commercial vendor and a stock solution of bacterial cells followed by incubation on agar plates. Recombinant genomes are identified by visual inspection using green-white colony screening akin to classical blue-white screening for recombinant plasmids. The modular nature of SEGA allows precise multi-level control of transcriptional, translational, and post-translational regulation. The SEGA architecture simultaneously supports increased standardization of genetic designs and a broad application range by utilizing well-characterized parts optimized for robust performance in the context of the bacterial genome. Ultimately, its adaption and expansion by the scientific community should improve predictability and comparability of experimental outcomes across different laboratories.

The spatial position effect: synthetic biology enters the era of 3D genomics

Microbial cell factories are critical to achieving green biomanufacturing. A position effect occurs when a synthetic gene circuit is expressed from different positions in the chassis strain genome. Here, we propose the concept of the 'spatial position effect,' which uses technologies in 3D genomics to reveal the spatial structure characteristics of the 3D genome of the chassis. On this basis, we propose to rationally design the integration sites of synthetic gene circuits, use reporter genes for preliminary screening, and integrate synthetic gene circuits into promising sites for further experiments. This approach can produce stable and efficient chassis strains for green biomanufacturing. The proposed spatial position effect brings synthetic biology into the era of 3D genomics.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Nucleoid-associated proteins shape chromatin structure and transcriptional regulation across the bacterial kingdom

Genome architecture has proven to be critical in determining gene regulation across almost all domains of life. While many of the key components and mechanisms of eukaryotic genome organization have been described, the interplay between bacterial DNA organization and gene regulation is only now being fully appreciated. An increasing pool of evidence has demonstrated that the bacterial chromosome can reasonably be thought of as chromatin, and that bacterial chromosomes contain transcriptionally silent and transcriptionally active regions analogous to heterochromatin and euchromatin, respectively. The roles played by histones in eukaryotic systems appear to be shared across a range of nucleoid-associated proteins (NAPs) in bacteria, which function to compact, structure, and regulate large portions of bacterial chromosomes. The broad range of extant NAPs, and the extent to which they differ from species to species, has raised additional challenges in identifying and characterizing their roles in all but a handful of model bacteria. Here we review the regulatory roles played by NAPs in several well-studied bacteria and use the resulting state of knowledge to provide a working definition for NAPs, based on their function, binding pattern, and expression levels. We present a screening procedure which can be applied to any species for which transcriptomic data are available. Finally, we note that NAPs tend to play two major regulatory roles - xenogeneic silencers and developmental regulators - and that many unrecognized potential NAPs exist in each bacterial species examined.

Modular optimization in metabolic engineering

There is an increasing demand for bioproducts produced by metabolically engineered microbes, such as pharmaceuticals, biofuels, biochemicals and other high value compounds. In order to meet this demand, modular optimization, the optimizing of subsections instead of the whole system, has been adopted to engineer cells to overproduce products. Research into modularity has focused on traditional approaches such as DNA, RNA, and protein-level modularity of intercellular machinery, by optimizing metabolic pathways for enhanced production. While research into these traditional approaches continues, limitations such as scale-up and time cost hold them back from wider use, while at the same time there is a shift to more novel methods, such as moving from episomal expression to chromosomal integration. Recently, nontraditional approaches such as co-culture systems and cell-free metabolic engineering (CFME) are being investigated for modular optimization. Co-culture modularity looks to optimally divide the metabolic burden between different hosts. CFME seeks to modularly optimize metabolic pathways in vitro, both speeding up the design of such systems and eliminating the issues associated with live hosts. In this review we will examine both traditional and nontraditional approaches for modular optimization, examining recent developments and discussing issues and emerging solutions for future research in metabolic engineering.

Dynamic landscape of protein occupancy across the Escherichia coli chromosome

Free-living bacteria adapt to environmental change by reprogramming gene expression through precise interactions of hundreds of DNA-binding proteins. A predictive understanding of bacterial physiology requires us to globally monitor all such protein-DNA interactions across a range of environmental and genetic perturbations. Here, we show that such global observations are possible using an optimized version of in vivo protein occupancy display technology (in vivo protein occupancy display-high resolution, IPOD-HR) and present a pilot application to Escherichia coli. We observe that the E. coli protein-DNA interactome organizes into 2 distinct prototypic features: (1) highly dynamic condition-dependent transcription factor (TF) occupancy; and (2) robust kilobase scale occupancy by nucleoid factors, forming silencing domains analogous to eukaryotic heterochromatin. We show that occupancy dynamics across a range of conditions can rapidly reveal the global transcriptional regulatory organization of a bacterium. Beyond discovery of previously hidden regulatory logic, we show that these observations can be utilized to computationally determine sequence specificity models for the majority of active TFs. Our study demonstrates that global observations of protein occupancy combined with statistical inference can rapidly and systematically reveal the transcriptional regulatory and structural features of a bacterial genome. This capacity is particularly crucial for non-model bacteria that are not amenable to routine genetic manipulation.

A Hi-C data-integrated model elucidates E. coli chromosome's multiscale organization at various replication stages

The chromosome of Escherichia coli is riddled with multi-faceted complexity. The emergence of chromosome conformation capture techniques are providing newer ways to explore chromosome organization. Here we combine a beads-on-a-spring polymer-based framework with recently reported Hi-C data for E. coli chromosome, in rich growth condition, to develop a comprehensive model of its chromosome at 5 kb resolution. The investigation focuses on a range of diverse chromosome architectures of E. coli at various replication states corresponding to a collection of cells, individually present in different stages of cell cycle. The Hi-C data-integrated model captures the self-organization of E. coli chromosome into multiple macrodomains within a ring-like architecture. The model demonstrates that the position of oriC is dependent on architecture and replication state of chromosomes. The distance profiles extracted from the model reconcile fluorescence microscopy and DNA-recombination assay experiments. Investigations into writhe of the chromosome model reveal that it adopts helix-like conformation with no net chirality, earlier hypothesized in experiments. A genome-wide radius of gyration map captures multiple chromosomal interaction domains and identifies the precise locations of rrn operons in the chromosome. We show that a model devoid of Hi-C encoded information would fail to recapitulate most genomic features unique to E. coli.

Engineering Transcriptional Interference through RNA Polymerase Processivity Control

Antisense transcription is widespread in all kingdoms of life and has been shown to influence gene expression through transcriptional interference (TI), a phenomenon in which one transcriptional process negatively influences another in cis. The processivity, or uninterrupted transcription, of an RNA polymerase (RNAP) is closely tied to levels of antisense transcription in bacterial genomes, but its influence on TI, while likely important, is not well-characterized. Here, we show that TI can be tuned through processivity control via three distinct antitermination strategies: the antibiotic bicyclomycin, phage protein Psu, and ribosome-RNAP coupling. We apply these methods toward TI and tune ribosome-RNAP coupling to produce 38-fold transcription-level gene repression due to both RNAP collisions and antisense RNA interference. We then couple protein roadblock and TI to design minimal genetic NAND and NOR logic gates. Together, these results show the importance of processivity control for strong TI and demonstrate TI's potential for synthetic biology.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Consequences of producing DNA gyrase from a synthetic gyrBA operon in Salmonella enterica serovar Typhimurium

DNA gyrase is an essential type II topoisomerase that is composed of two subunits, GyrA and GyrB, and has an A₂B₂ structure. Although the A and B subunits are required in equal proportions to form DNA gyrase, the gyrA and gyrB genes that encode them in Salmonella (and in many other bacteria) are at separate locations on the chromosome, are under separate transcriptional control, and are present in different copy numbers in rapidly growing bacteria. In wild‐type Salmonella, gyrA is near the chromosome's replication terminus, while gyrB is near the origin. We generated a synthetic gyrBA operon at the oriC‐proximal location of gyrB to test the significance of the gyrase gene position for Salmonella physiology. Although the strain producing gyrase from an operon had a modest alteration to its DNA supercoiling set points, most housekeeping functions were unaffected. However, its SPI‐2 virulence genes were expressed at a reduced level and its survival was reduced in macrophage. Our data reveal that the horizontally acquired SPI‐2 genes have a greater sensitivity to disturbance of DNA topology than the core genome and we discuss its significance in the context of Salmonella genome evolution and the gyrA and gyrB gene arrangements found in other bacteria.

Optimized gene expression from bacterial chromosome by high-throughput integration and screening

Chromosomal integration of recombinant genes is desirable compared with expression from plasmids due to increased stability, reduced cell-to-cell variability, and elimination of the need for antibiotics for plasmid maintenance. Here, we present a new approach for tuning pathway gene expression levels via random integration and high-throughput screening. We demonstrate multiplexed gene integration and expression-level optimization for isobutanol production in Escherichia coli The integrated strains could, with far lower expression levels than plasmid-based expression, produce high titers (10.0 ± 0.9 g/liter isobutanol in 48 hours) and yields (69% of the theoretical maximum). Close examination of pathway expression in the top-performing, as well as other isolates, reveals the complexity of cellular metabolism and regulation, underscoring the need for precise optimization while integrating pathway genes into the chromosome. We expect this method for pathway integration and optimization can be readily extended to a wide range of pathways and chassis to create robust and efficient production strains.

Divergently Transcribed ncRNAs in Escherichia coli: Refinement of the Transcription Starts Assumes Functional Diversification

Non-coding regulatory RNAs (ncRNAs) comprise specialized group of essential genetically encoded biological molecules involved in the wide variety of cellular metabolic processes. The progressive increase in the number of newly identified ncRNAs and the defining of their genome location indicate their predominant nesting in intergenic regions and expression under the control of their own regulatory elements. At the same time, the regulation of ncRNA's transcription cannot be considered in isolation from the processes occurring in the immediate genetic environment. A number of experimental data indicate the notable impact of positional regulation of gene expression mediated by dynamic temporal DNA rearrangements accompanying transcription events in the vicinity of neighboring genes. This issue can be perceived as particularly significant for divergently transcribed ncRNAs being actually subjected to double regulatory pressure. Based on available results of RNAseq experiments for Escherichia coli, we screened out divergent ncRNAs and the adjacent genes for the exact positions of transcription start sites (TSSs) and relative efficiency of RNA production. This analysis revealed extension or shortening of some previously annotated ncRNAs resulting in modified secondary structure, confirmed stable expression of four ncRNAs annotated earlier as putative, and approved the possibility of expression of divergently transcribed ncRNAs containing repetitive extragenic palindromic (REP) elements. The biogenesis of secreted ncRNAs from divergently transcribed ffs, chiX, ralA, and ryhB is discussed taking into account positions of TSSs. Refinement of TSSs for the neighboring genes renders some ncRNAs as true antisense overlapping with 5'UTR of divergently transcribed mRNAs.

Random Chromosomal Integration and Screening Yields E. coli K-12 Derivatives Capable of Efficient Sucrose Utilization

Chromosomal expression of heterologous genes offers stability and maintenance advantages over episomal expression, yet remains difficult to optimize through site-specific integration. The challenge has in large part been due to the variability of chromosomal gene expression, which has only recently been shown to be affected by multiple factors, including the local genomic context. In this work we utilize Tn5 transposase to randomly integrate a three-gene csc operon encoding nonphosphotransferase sucrose catabolism into the E. coli K-12 chromosome. Isolates from the transposon library yielded a range of growth rates on sucrose as the sole carbon source, including some that were comparable to that of E. coli K-12 on glucose (μ_max = 0.70 ± 0.03 h^-1). Narrowness of the growth rate distributions and faster growth compared to plasmids indicate that efficient csc expression is attainable. Furthermore, enhanced growth rate upon transduction into strains that underwent adaptive laboratory evolution indicate that sucrose catabolism is not limiting to cellular growth. We also show that transduction of a csc fast-growth locus into an isobutanol production strain yields high titer (7.56 ± 0.25 g/L) on sucrose as the sole carbon source. Our results demonstrate that random integration is an effective strategy for optimizing heterologous expression within the context of cellular metabolism for both fast growth and biochemical production phenotypes.

Genome Complexity Browser: Visualization and quantification of genome variability

Comparative genomics studies may be used to acquire new knowledge regarding genome architecture, which defines the rules for combining sets of genes in the genome of living organisms. Hundreds of thousands of prokaryotic genomes have been sequenced and assembled. However, computational tools capable of simultaneously comparing large numbers of genomes are lacking. We developed the Genome Complexity Browser, a tool that allows the visualization of gene contexts, in a graph-based format, and the quantification of variability for different segments of a genome. The graph-based visualization allows the inspection of changes in gene contents and neighborhoods across hundreds of genomes, simultaneously, which may facilitate the identification of conserved and variable segments of operons or the estimation of the overall variability associated with a particular genome locus. We introduced a measure called complexity, to quantify genome variability. Intraspecies and interspecies comparisons revealed that regions with high complexity values tended to be located in areas that are conserved across different strains and species.

The comparison of genomes among different bacteria and archaea species has revealed that many species frequently exchange genes. Occasionally, such horizontal gene transfer events result in the acquisition of pathogenic properties or antibiotic resistance in the recipient organism. Previously, the probabilities of gene insertions were found to vary, with unequal distributions along a chromosome. At some loci, referred to as hotspots, changes occur with much higher frequencies compared with other regions of the chromosome. We developed a computational method and a software tool, called Genome Complexity Browser, that allows the identification of genome variability hotspots and the visualization of changes. We compared the localization of various hotspots and revealed that some demonstrate conserved localizations, even across species, whereas others are transient. Our tool allows users to visually inspect the patterns of gene changes in graph-based format, which presents the visualization in a format that is both compact and informative.