Evolvability and hierarchy in rewired bacterial gene networks

Engineering Tunable, Low Latency Spatial Computation with Dual Input Quorum Sensing Promoters

Quorum sensing signals have evolved for population-level signaling in bacterial communities and are versatile tools for engineering cell-cell signaling in synthetic biology projects. Here, we characterize the spatial diffusion of a palette of quorum sensing signals and find that their diffusion in agar can be predicted from their molecular weight with a simple power law. We also engineer novel dual- and multi-input promoters that respond to quorum-sensing diffusive signals for use in engineered genetic systems. We engineer a promoter scaffold that can be adapted for activation and repression by multiple diffusers simultaneously. Lastly, we combine the knowledge on diffusion dynamics with the novel genetic components to build a new generation of spatial, stripe-forming systems with a simplified design, improved robustness, tuneability, and response time.

S-layer is a key element in metabolic response and entry into the stationary phase in Bacillus cereus AH187

Bacillus cereus is a food-borne Gram-positive pathogen. The emetic reference strain B. cereus AH187 is surrounded by a proteinaceous surface layer (S-layer) that contributes to its physico-chemical surface properties, and promotes its adhesion in response to starvation conditions. The S-layer produced by B. cereus AH187 is composed of two proteins, SL2 and EA1, which are incorporated at different growth stages. Here, we showed that deletion of the genes encoding SL2 and EA1 produced viable cells, but decreased the glucose uptake rate at the start of growth, and induced extensive reorganization of the cellular and exoproteomes upon entry into the stationary phase. As a consequence, stationary cells were less resistant to abiotic stress. Taken together, our data indicate that the S-layer is crucial but comes at a metabolic cost that modulates the stationary phase response. SIGNIFICANCE: The emetic strains of Bacillus cereus are known to cause severe food poisoning, making it crucial to understand the factors contributing to their selective enrichment in foods. Most emetic strains are surrounded by a crystalline S-layer, which is a costly protein structure to produce. In this study, we used high-throughput proteomics to investigate how S-layer synthesis affects the allocation of cellular resources in the emetic B. cereus strain AH187. Our results demonstrate that the synthesis of the S-layer plays a crucial role in the pathogen's ability to thrive under stationary growth phase conditions by modulating the stress response, thereby promoting its lifestyle as an emetic pathogen. We conclude that the synthesis of the S-layer is a critical adaptation for emetic B. cereus to successfully colonize specific niches.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Genome-wide promoter responses to CRISPR perturbations of regulators reveal regulatory networks in Escherichia coli

Elucidating genome-scale regulatory networks requires a comprehensive collection of gene expression profiles, yet measuring gene expression responses for every transcription factor (TF)-gene pair in living prokaryotic cells remains challenging. Here, we develop pooled promoter responses to TF perturbation sequencing (PPTP-seq) via CRISPR interference to address this challenge. Using PPTP-seq, we systematically measure the activity of 1372 Escherichia coli promoters under single knockdown of 183 TF genes, illustrating more than 200,000 possible TF-gene responses in one experiment. We perform PPTP-seq for E. coli growing in three different media. The PPTP-seq data reveal robust steady-state promoter activities under most single TF knockdown conditions. PPTP-seq also enables identifications of, to the best of our knowledge, previously unknown TF autoregulatory responses and complex transcriptional control on one-carbon metabolism. We further find context-dependent promoter regulation by multiple TFs whose relative binding strengths determined promoter activities. Additionally, PPTP-seq reveals different promoter responses in different growth media, suggesting condition-specific gene regulation. Overall, PPTP-seq provides a powerful method to examine genome-wide transcriptional regulatory networks and can be potentially expanded to reveal gene expression responses to other genetic elements.

Transcription factor expression levels and environmental signals constrain transcription factor innovation

Evolutionary innovation of transcription factors frequently drives phenotypic diversification and adaptation to environmental change. Transcription factors can gain or lose connections to target genes, resulting in novel regulatory responses and phenotypes. However the frequency of functional adaptation varies between different regulators, even when they are closely related. To identify factors influencing propensity for innovation, we utilise a Pseudomonas fluorescens SBW25 strain rendered incapable of flagellar mediated motility in soft-agar plates via deletion of the flagellar master regulator (fleQ ). This bacterium can evolve to rescue flagellar motility via gene regulatory network rewiring of an alternative transcription factor to rescue activity of FleQ. Previously, we have identified two members (out of 22) of the RpoN-dependent enhancer binding protein (RpoN-EBP) family of transcription factors (NtrC and PFLU1132) that are capable of innovating in this way. These two transcription factors rescue motility repeatably and reliably in a strict hierarchy – with NtrC the only route in a ∆fleQ background, and PFLU1132 the only route in a ∆fleQ ∆ntrC background. However, why other members in the same transcription factor family have not been observed to rescue flagellar activity is unclear. Previous work shows that protein homology cannot explain this pattern within the protein family (RpoN-EBPs), and mutations in strains that rescued motility suggested high levels of transcription factor expression and activation drive innovation. We predict that mutations that increase expression of the transcription factor are vital to unlock evolutionary potential for innovation. Here, we construct titratable expression mutant lines for 11 of the RpoN-EBPs in P. fluorescens . We show that in five additional RpoN-EBPs (FleR, HbcR, GcsR, DctD, AauR and PFLU2209), high expression levels result in different mutations conferring motility rescue, suggesting alternative rewiring pathways. Our results indicate that expression levels (and not protein homology) of RpoN-EBPs are a key constraining factor in determining evolutionary potential for innovation. This suggests that transcription factors that can achieve high expression through few mutational changes, or transcription factors that are active in the selective environment, are more likely to innovate and contribute to adaptive gene regulatory network evolution.

Robustness and innovation in synthetic genotype networks

Genotype networks are sets of genotypes connected by small mutational changes that share the same phenotype. They facilitate evolutionary innovation by enabling the exploration of different neighborhoods in genotype space. Genotype networks, first suggested by theoretical models, have been empirically confirmed for proteins and RNAs. Comparative studies also support their existence for gene regulatory networks (GRNs), but direct experimental evidence is lacking. Here, we report the construction of three interconnected genotype networks of synthetic GRNs producing three distinct phenotypes in Escherichia coli. Our synthetic GRNs contain three nodes regulating each other by CRISPR interference and governing the expression of fluorescent reporters. The genotype networks, composed of over twenty different synthetic GRNs, provide robustness in face of mutations while enabling transitions to innovative phenotypes. Through realistic mathematical modeling, we quantify robustness and evolvability for the complete genotype-phenotype map and link these features mechanistically to GRN motifs. Our work thereby exemplifies how GRN evolution along genotype networks might be driving evolutionary innovation.

A primer to directed evolution: current methodologies and future directions

This review summarises the methods available for directed evolution, including mutagenesis and variant selection techniques. The advantages and disadvantages of each technique are presented, and future challenges in the field are discussed.

Robustness in phenotypic plasticity and heterogeneity patterns enabled by EMT networks

Epithelial-mesenchymal plasticity (EMP) is a key arm of cancer metastasis and is observed across many contexts. Cells undergoing EMP can reversibly switch between three classes of phenotypes: epithelial (E), mesenchymal (M), and hybrid E/M. While a large number of multistable regulatory networks have been identified to be driving EMP in various contexts, the exact mechanisms and design principles that enable robustness in driving EMP across contexts are not yet fully understood. Here, we investigated dynamic and structural robustness in EMP networks with regard to phenotypic heterogeneity and plasticity. We use two different approaches to simulate these networks: a computationally inexpensive, parameter-independent continuous state space Boolean model, and an ODE-based parameter-agnostic framework (RACIPE), both of which yielded similar phenotypic distributions. While the latter approach is useful for measurements of plasticity, the former model enabled us to extensively investigate robustness in phenotypic heterogeneity. Using perturbations to network topology and by varying network parameters, we show that multistable EMP networks are structurally and dynamically more robust compared with their randomized counterparts, thereby highlighting their topological hallmarks. These features of robustness are governed by a balance of positive and negative feedback loops embedded in these networks. Using a combination of the number of negative and positive feedback loops weighted by their lengths, we identified a metric that can explain the structural and dynamical robustness of these networks. This metric enabled us to compare networks across multiple sizes, and the network principles thus obtained can be used to identify fragilities in large networks without simulating their dynamics. Our analysis highlights a network topology-based approach to quantify robustness in the phenotypic heterogeneity and plasticity emergent from EMP networks.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Mechanisms of antibiotic action shape the fitness landscapes of resistance mutations

Antibiotic-resistant pathogens are a major public health threat. A deeper understanding of how an antibiotic's mechanism of action influences the emergence of resistance would aid in the design of new drugs and help to preserve the effectiveness of existing ones. To this end, we developed a model that links bacterial population dynamics with antibiotic-target binding kinetics. Our approach allows us to derive mechanistic insights on drug activity from population-scale experimental data and to quantify the interplay between drug mechanism and resistance selection. We find that both bacteriostatic and bactericidal agents can be equally effective at suppressing the selection of resistant mutants, but that key determinants of resistance selection are the relationships between the number of drug-inactivated targets within a cell and the rates of cellular growth and death. We also show that heterogeneous drug-target binding within a population enables resistant bacteria to evolve fitness-improving secondary mutations even when drug doses remain above the resistant strain's minimum inhibitory concentration. Our work suggests that antibiotic doses beyond this "secondary mutation selection window" could safeguard against the emergence of high-fitness resistant strains during treatment.

Emergent expression of fitness-conferring genes by phenotypic selection

Genotypic and phenotypic adaptation is the consequence of ongoing natural selection in populations and is key to predicting and preventing drug resistance. Whereas classic antibiotic persistence is all-or-nothing, here we demonstrate that an antibiotic resistance gene displays linear dose-responsive selection for increased expression in proportion to rising antibiotic concentration in growing Escherichia coli populations. Furthermore, we report the potentially wide-spread nature of this form of emergent gene expression (EGE) by instantaneous phenotypic selection process under bactericidal and bacteriostatic antibiotic treatment, as well as an amino acid synthesis pathway enzyme under a range of auxotrophic conditions. We propose an analogy to Ohm's law in electricity (V = IR), where selection pressure acts similarly to voltage (V), gene expression to current (I), and resistance (R) to cellular machinery constraints and costs. Lastly, mathematical modeling using agent-based models of stochastic gene expression in growing populations and Bayesian model selection reveal that the EGE mechanism requires variability in gene expression within an isogenic population, and a cellular "memory" from positive feedbacks between growth and expression of any fitness-conferring gene. Finally, we discuss the connection of the observed phenomenon to a previously described general fluctuation-response relationship in biology.

Reducing metabolic burden in the PACEmid evolver system by remastering high-copy phagemid vectors

Orthogonal or non-cross-reacting transcription factors are used in synthetic biology as components of genetic circuits. Brödel et al. (2016) engineered 12 such cIλ transcription factor variants using a directed evolution 'PACEmid' system. The variants operate as dual activator/repressors and expand gene circuit construction possibilities. However, the high-copy phagemid vectors carrying the cIλ variants imposed high metabolic burden upon cells. Here, the authors 'remaster' the phagemid backbones to relieve their burden substantially, exhibited by a recovery in Escherichia coli growth. The remastered phagemids' ability to function within the PACEmid evolver system is maintained, as is the cIλ transcription factors' activity within these vectors. The low-burden phagemid versions are more suitable for use in PACEmid experiments and synthetic gene circuits; the authors have, therefore, replaced the original high-burden phagemids on the Addgene repository. The authors' work emphasises the importance of understanding metabolic burden and incorporating it into design steps in future synthetic biology ventures.

A synthetic synthesis to explore animal evolution and development

Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Evolution enhances mutational robustness and suppresses the emergence of a new phenotype: A new computational approach for studying evolution

The aim of this paper is two-fold. First, we propose a new computational method to investigate the particularities of evolution. Second, we apply this method to a model of gene regulatory networks (GRNs) and explore the evolution of mutational robustness and bistability. Living systems have developed their functions through evolutionary processes. To understand the particularities of this process theoretically, evolutionary simulation (ES) alone is insufficient because the outcomes of ES depend on evolutionary pathways. We need a reference system for comparison. An appropriate reference system for this purpose is an ensemble of the randomly sampled genotypes. However, generating high-fitness genotypes by simple random sampling is difficult because such genotypes are rare. In this study, we used the multicanonical Monte Carlo method developed in statistical physics to construct a reference ensemble of GRNs and compared it with the outcomes of ES. We obtained the following results. First, mutational robustness was significantly higher in ES than in the reference ensemble at the same fitness level. Second, the emergence of a new phenotype, bistability, was delayed in evolution. Third, the bistable group of GRNs contains many mutationally fragile GRNs compared with those in the non-bistable group. This suggests that the delayed emergence of bistability is a consequence of the mutation-selection mechanism.

Living systems are products of evolution, and their present forms reflect their evolutionary history. Thus, to investigate the particularity of the evolutionary process by computer simulations, an appropriate reference system is needed for comparison with the outcomes of evolutionary simulations. In this study, we considered a model of gene regulatory networks (GRNs). Our idea was to construct a reference ensemble comprising randomly generated GRNs. To produce GRNs with high fitness values, which are rare, we employed a “rare event sampling” method developed in statistical physics. In particular, we focused on the evolution of mutational robustness. Living systems do not lose viability readily, even when some genes are mutated. This trait, called mutational robustness, has developed throughout evolution, along with functionality. Using the abovementioned method, we found that mutational robustness resulting from evolution exceeded that of the reference set. Therefore, mutational robustness is enhanced by evolution. We also found that the emergence of a new phenotype was significantly delayed in evolution. Our results suggest that this delay is a consequence of the fact that mutationally robust GRNs are favored by evolution.

Metacode: One code to rule them all

The code of codes or metacode is a microcosm where biological layers, as well as their codes, interact together allowing the continuity of information flow in organisms by increasing biological entities' complexity. Through this novel organic code, biological systems scale towards niches with higher informatic freedom building structures that increase the entropy in the universe. Code biology has developed a novel informational framework where biological entities strive themselves through the information flow carried out through organic codes consisting of two molecular or functional landscapes intertwined through arbitrary linkages via an adaptor whose nature is autonomous from molecular determinism. Here we will integrate genomic and epigenomic codes according to the evidence released in ENCODE (phase 3), psychENCODE and GTEx project, outlining the principles of the metacode, to address the continuous nature of biological systems and their inter-layered information flow. This novel complex metacode maps from very constrained sets of elements (i.e., regulation sites modulating gene expression) to new ones with greater freedom of decoding (i.e., a continuous cell phenotypic space). This leads to a new domain in code biology where biological systems are informatic attractors that navigate an energy metaspace through a complexity-noise balance, stalling in emergent niches where organic codes take meaning.

Engineering Pseudomonas putida for efficient aromatic conversion to bioproduct using high throughput screening in a bioreactor

Pseudomonas putida KT2440 is an emerging biomanufacturing host amenable for use with renewable carbon streams including aromatics such as para-coumarate. We used a pooled transposon library disrupting nearly all (4,778) non-essential genes to characterize this microbe under common stirred-tank bioreactor parameters with quantitative fitness assays. Assessing differential fitness values by monitoring changes in mutant strain abundance identified 33 gene mutants with improved fitness across multiple stirred-tank bioreactor formats. Twenty-one deletion strains from this subset were reconstructed, including GacA, a regulator, TtgB, an ABC transporter, and PP_0063, a lipid A acyltransferase. Thirteen deletion strains with roles in varying cellular functions were evaluated for conversion of para-coumarate, to a heterologous bioproduct, indigoidine. Several mutants, such as the ΔgacA strain improved fitness in a bioreactor by 35 fold and showed an 8-fold improvement in indigoidine production (4.5 g/L, 0.29 g/g, 23% of maximum theoretical yield) from para-coumarate as the carbon source.

Elucidation of the Pathogenicity-Associated Regulatory Network in Xanthomonas oryzae pv. oryzae

Xanthomonas is a notorious plant pathogen causing serious diseases in hundreds of plant hosts. Xanthomonas species are equipped with an array of signal transduction systems that regulate gene expression to survive in various harsh environments and successfully infect hosts. Although certain pathogenicity-associated regulators have been functionally characterized, signal transduction systems always function as a regulatory network which remains to be elucidated in Xanthomonas. This study used a systematic approach to characterize all identified pathogenicity-associated regulators in Xanthomonas oryzae pv. oryzae (Xoo), including a transcriptional regulator with unknown function, and their interactive regulatory network. RNA sequencing was used in elucidating the patterns of the 10 pathogenicity-associated regulators identified. Results revealed that each pathogenicity-associated regulator has cross talk with others and all these regulators function as a regulatory network, with VemR and PXO_RS20790 being the master pathogenicity-associated regulators and HrpX being the final executant. Moreover, regulome analysis showed that numerous genes other than genes in pathogenicity islands are finely regulated within the regulatory network. Given that most of the pathogenicity-associated regulators are conserved in Xanthomonadales, our findings suggest a global network of gene regulation in this evolutionarily conserved pathogen. In conclusion, our study provides essential basic information about the regulatory network in Xoo, suggesting that this complicated regulatory network is one of the reasons for the robustness and fitness of Xanthomonas spp.

IMPORTANCE The host plant infection process of pathogenic bacteria is a coordinating cellular behavior, which requires dynamic regulation at several levels in response to variations in host plants or fluctuations in the external environment. As one of the most important genera of plant-pathogenic bacteria, Xanthomonas has been studied as a model. Although certain pathogenicity-associated regulators have been functionally characterized, interactions among them remain to be elucidated. This study systematically characterized pathogenicity-associated regulators in Xoo and revealed that cross talk exists among pathogenicity-associated regulators and function as a regulatory network in which a hierarchy exists among the regulators. Our study elucidated the landscape of the pathogenicity-associated regulatory network in Xanthomonas, promoting understanding of the infection process of pathogenic bacteria.

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.

A quantitative method for proteome reallocation using minimal regulatory interventions

Engineering resource allocation in biological systems is an ongoing challenge. Organisms allocate resources for ensuring survival, reducing the productivity of synthetic biology functions. Here we present a new approach for engineering the resource allocation of Escherichia coli by rationally modifying its transcriptional regulatory network. Our method (ReProMin) identifies the minimal set of genetic interventions that maximizes the savings in cell resources. To this end, we categorized transcription factors according to the essentiality of its targets and we used proteomic data to rank them. We designed the combinatorial removal of transcription factors that maximize the release of resources. Our resulting strain containing only three mutations, theoretically releasing 0.5% of its proteome, had higher proteome budget, increased production of an engineered metabolic pathway and showed that the regulatory interventions are highly specific. This approach shows that combining proteomic and regulatory data is an effective way of optimizing strains using conventional molecular methods.

Unraveling the Complex Interplay of Fis and IHF Through Synthetic Promoter Engineering

Bacterial promoters are usually formed by multiple cis-regulatory elements recognized by a plethora of transcriptional factors (TFs). From those, global regulators are key elements since these TFs are responsible for the regulation of hundreds of genes in the bacterial genome. For instance, Fis and IHF are global regulators that play a major role in gene expression control in Escherichia coli, and usually, multiple cis-regulatory elements for these proteins are present at target promoters. Here, we investigated the relationship between the architecture of the cis-regulatory elements for Fis and IHF in E. coli. For this, we analyze 42 synthetic promoter variants harboring consensus cis-elements for Fis and IHF at different distances from the core -35/-10 region and in various numbers and combinations. We first demonstrated that although Fis preferentially recognizes its consensus cis-element, it can also recognize, to some extent, the consensus-binding site for IHF, and the same was true for IHF, which was also able to recognize Fis binding sites. However, changing the arrangement of the cis-elements (i.e., the position or number of sites) can completely abolish the non-specific binding of both TFs. More remarkably, we demonstrated that combining cis-elements for both TFs could result in Fis and IHF repressed or activated promoters depending on the final architecture of the promoters in an unpredictable way. Taken together, the data presented here demonstrate how small changes in the architecture of bacterial promoters could result in drastic changes in the final regulatory logic of the system, with important implications for the understanding of natural complex promoters in bacteria and their engineering for novel applications.

Emergence of cooperative bistability and robustness of gene regulatory networks

Gene regulatory networks (GRNs) are complex systems in which many genes regulate mutually to adapt the cell state to environmental conditions. In addition to function, the GRNs possess several kinds of robustness. This robustness means that systems do not lose their functionality when exposed to disturbances such as mutations or noise, and is widely observed at many levels in living systems. Both function and robustness have been acquired through evolution. In this respect, GRNs utilized in living systems are rare among all possible GRNs. In this study, we explored the fitness landscape of GRNs and investigated how robustness emerged in highly-fit GRNs. We considered a toy model of GRNs with one input gene and one output gene. The difference in the expression level of the output gene between two input states, "on" and "off", was considered as fitness. Thus, the determination of the fitness of a GRN was based on how sensitively it responded to the input. We employed the multicanonical Monte Carlo method, which can sample GRNs randomly in a wide range of fitness levels, and classified the GRNs according to their fitness. As a result, the following properties were found: (1) Highly-fit GRNs exhibited bistability for intermediate input between "on" and "off". This means that such GRNs responded to two input states by using different fixed points of dynamics. This bistability emerges necessarily as fitness increases. (2) These highly-fit GRNs were robust against noise because of their bistability. In other words, noise robustness is a byproduct of high fitness. (3) GRNs that were robust against mutations were not extremely rare among the highly-fit GRNs. This implies that mutational robustness is readily acquired through the evolutionary process. These properties are universal irrespective of the evolutionary pathway, because the results do not rely on evolutionary simulation.

Programmable CRISPR-Cas transcriptional activation in bacteria

Programmable gene activation enables fine-tuned regulation of endogenous and synthetic gene circuits to control cellular behavior. While CRISPR-Cas-mediated gene activation has been extensively developed for eukaryotic systems, similar strategies have been difficult to implement in bacteria. Here, we present a generalizable platform for screening and selection of functional bacterial CRISPR-Cas transcription activators. Using this platform, we identified a novel CRISPR activator, dCas9-AsiA, that could activate gene expression by more than 200-fold across genomic and plasmid targets with diverse promoters after directed evolution. The evolved dCas9-AsiA can simultaneously mediate activation and repression of bacterial regulons in E. coli. We further identified hundreds of promoters with varying basal expression that could be induced by dCas9-AsiA, which provides a rich resource of genetic parts for inducible gene activation. Finally, we show that dCas9-AsiA can be ported to other bacteria of clinical and bioindustrial relevance, thus enabling bacterial CRISPRa in more application areas. This work expands the toolbox for programmable gene regulation in bacteria and provides a useful resource for future engineering of other bacterial CRISPR-based gene regulators.

Accelerated evolution of a minimal 63-amino acid dual transcription factor

Transcription factors control gene expression in all life. This raises the question of what is the smallest protein that can support such activity. In nature, Cro from bacteriophage λ is one of the smallest known repressors (66 amino acids), and activators are typically much larger (e.g., λ cI, 237 amino acids). Previous efforts to engineer a minimal activator from λ Cro resulted in no activity in vivo in cells. In this study, we show that directed evolution results in a new Cro activator-repressor that functions as efficiently as λ cI in vivo. To achieve this, we develop phagemid-assisted continuous evolution (PACEmid). We find that a peptide as small as 63 amino acids functions efficiently as an activator and/or repressor. To our knowledge, this is the smallest protein activator that enables polymerase recruitment, highlighting the capacity of transcription factors to evolve from very short peptide sequences.

The Phenotype Paradox: Lessons From Natural Transcriptome Evolution on How to Engineer Plants

Plants have evolved genome complexity through iterative rounds of single gene and whole genome duplication. This has led to substantial expansion in transcription factor numbers following preferential retention and subsequent functional divergence of these regulatory genes. Here we review how this simple evolutionary network rewiring process, regulatory gene duplication followed by functional divergence, can be used to inspire synthetic biology approaches that seek to develop novel phenotypic variation for future trait based breeding programs in plants.

Evolutionary regain of lost gene circuit function

Natural or synthetic genetic modules can lose their function over long-term evolution if the function is costly. How populations can evolve to restore such broken function is poorly understood. To test the reversibility of evolutionary breakdown, we use yeast cell populations with a chromosomally integrated synthetic gene circuit. In previous evolution experiments the gene circuit lost its costly function through various mutations. By exposing such mutant populations to conditions where regaining gene circuit function would be beneficial, we find adaptation scenarios with or without repairing lost gene circuit function. These results are important for drug resistance or future synthetic biology applications where evolutionary loss and regain of function play a significant role.

Evolutionary reversibility—the ability to regain a lost function—is an important problem both in evolutionary and synthetic biology, where repairing natural or synthetic systems broken by evolutionary processes may be valuable. Here, we use a synthetic positive-feedback (PF) gene circuit integrated into haploid Saccharomyces cerevisiae cells to test if the population can restore lost PF function. In previous evolution experiments, mutations in a gene eliminated the fitness costs of PF activation. Since PF activation also provides drug resistance, exposing such compromised or broken mutants to both drug and inducer should create selection pressure to regain drug resistance and possibly PF function. Indeed, evolving 7 PF mutant strains in the presence of drug revealed 3 adaptation scenarios through genomic, PF-external mutations that elevate PF basal expression, possibly by affecting transcription, translation, degradation, and other fundamental cellular processes. Nonfunctional mutants gained drug resistance without ever developing high expression, while quasifunctional and dysfunctional PF mutants developed high expression nongenetically, which then diminished, although more slowly for dysfunctional mutants where revertant clones arose. These results highlight how intracellular context, such as the growth rate, can affect regulatory network dynamics and evolutionary dynamics, which has important consequences for understanding the evolution of drug resistance and developing future synthetic biology applications.

Adaptive walks on high-dimensional fitness landscapes and seascapes with distance-dependent statistics

The dynamics of evolution is intimately shaped by epistasis - interactions between genetic elements which cause the fitness-effect of combinations of mutations to be non-additive. Analyzing evolutionary dynamics that involves large numbers of epistatic mutations is intrinsically difficult. A crucial feature is that the fitness landscape in the vicinity of the current genome depends on the evolutionary history. A key step is thus developing models that enable study of the effects of past evolution on future evolution. In this work, we introduce a broad class of high-dimensional random fitness landscapes for which the correlations between fitnesses of genomes are a general function of genetic distance. Their Gaussian character allows for tractable computational as well as analytic understanding. We study the properties of these landscapes focusing on the simplest evolutionary process: random adaptive (uphill) walks. Conventional measures of "ruggedness" are shown to not much affect such adaptive walks. Instead, the long-distance statistics of epistasis cause all properties to be highly conditional on past evolution, determining the statistics of the local landscape (the distribution of fitness-effects of available mutations and combinations of these), as well as the global geometry of evolutionary trajectories. In order to further explore the effects of conditioning on past evolution, we model the effects of slowly changing environments. At long times, such fitness "seascapes" cause a statistical steady state with highly intermittent evolutionary dynamics: populations undergo bursts of rapid adaptation, interspersed with periods in which adaptive mutations are rare and the population waits for more new directions to be opened up by changes in the environment. Finally, we discuss prospects for studying more complex evolutionary dynamics and on broader classes of high-dimensional landscapes and seascapes.

Determination of the Gene Regulatory Network of a Genome-Reduced Bacterium Highlights Alternative Regulation Independent of Transcription Factors

Here, we determined the relative importance of different transcriptional mechanisms in the genome-reduced bacterium Mycoplasma pneumoniae, by employing an array of experimental techniques under multiple genetic and environmental perturbations. Of the 143 genes tested (21% of the bacterium’s annotated proteins), only 55% showed an altered phenotype, highlighting the robustness of biological systems. We identified nine transcription factors (TFs) and their targets, representing 43% of the genome, and 16 regulators that indirectly affect transcription. Only 20% of transcriptional regulation is mediated by canonical TFs when responding to perturbations. Using a Random Forest, we quantified the non-redundant contribution of different mechanisms such as supercoiling, metabolic control, RNA degradation, and chromosome topology to transcriptional changes. Model-predicted gene changes correlate well with experimental data in 95% of the tested perturbations, explaining up to 70% of the total variance when also considering noise. This analysis highlights the importance of considering non-TF-mediated regulation when engineering bacteria.

•

Full comprehensive reconstruction of a bacterial gene regulatory network achieved

•

Genome-reduced bacterium Mycoplasma pneumoniae is robust to genetic perturbations

•

Large part of transcription regulation in bacteria is transcription-factor independent

•

Transcription-factor-independent regulation has a smaller dynamic range

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

Inferring Gene Regulatory Networks from Multiple Datasets

Gaussian process dynamical systems (GPDS) represent Bayesian nonparametric approaches to inference of nonlinear dynamical systems, and provide a principled framework for the learning of biological networks from multiple perturbed time series measurements of gene or protein expression. Such approaches are able to capture the full richness of complex ODE models, and can be scaled for inference in moderately large systems containing hundreds of genes. Related hierarchical approaches allow for inference from multiple datasets in which the underlying generative networks are assumed to have been rewired, either by context-dependent changes in network structure, evolutionary processes, or synthetic manipulation. These approaches can also be used to leverage experimentally determined network structures from one species into another where the network structure is unknown. Collectively, these methods provide a comprehensive and flexible platform for inference from a diverse range of data, with applications in systems and synthetic biology, as well as spatiotemporal modelling of embryo development. In this chapter we provide an overview of GPDS approaches and highlight their applications in the biological sciences, with accompanying tutorials available as a Jupyter notebook from https://github.com/cap76/GPDS .

Evolutionary potential of transcription factors for gene regulatory rewiring

Gene regulatory networks evolve through rewiring of individual components-that is, through changes in regulatory connections. However, the mechanistic basis of regulatory rewiring is poorly understood. Using a canonical gene regulatory system, we quantify the properties of transcription factors that determine the evolutionary potential for rewiring of regulatory connections: robustness, tunability and evolvability. In vivo repression measurements of two repressors at mutated operator sites reveal their contrasting evolutionary potential: while robustness and evolvability were positively correlated, both were in trade-off with tunability. Epistatic interactions between adjacent operators alleviated this trade-off. A thermodynamic model explains how the differences in robustness, tunability and evolvability arise from biophysical characteristics of repressor-DNA binding. The model also uncovers that the energy matrix, which describes how mutations affect repressor-DNA binding, encodes crucial information about the evolutionary potential of a repressor. The biophysical determinants of evolutionary potential for regulatory rewiring constitute a mechanistic framework for understanding network evolution.

Understanding Biological Regulation Through Synthetic Biology

Engineering synthetic gene regulatory circuits proceeds through iterative cycles of design, building, and testing. Initial circuit designs must rely on often-incomplete models of regulation established by fields of reductive inquiry—biochemistry and molecular and systems biology. As differences in designed and experimentally observed circuit behavior are inevitably encountered, investigated, and resolved, each turn of the engineering cycle can force a resynthesis in understanding of natural network function. Here, we outline research that uses the process of gene circuit engineering to advance biological discovery. Synthetic gene circuit engineering research has not only refined our understanding of cellular regulation but furnished biologists with a toolkit that can be directed at natural systems to exact precision manipulation of network structure. As we discuss, using circuit engineering to predictively reorganize, rewire, and reconstruct cellular regulation serves as the ultimate means of testing and understanding how cellular phenotype emerges from systems-level network function.

A Postgenomic Perspective on Molecular Cytogenetics

BACKGROUND

The postgenomic era is featured by massive data collection and analyses from various large scale-omics studies. Despite the promising capability of systems biology and bioinformatics to handle large data sets, data interpretation, especially the translation of -omics data into clinical implications, has been challenging.

DISCUSSION

In this perspective, some important conceptual and technological limitations of current systems biology are discussed in the context of the ultimate importance of the genome beyond the collection of all genes. Following a brief summary of the contributions of molecular cytogenetics/cytogenomics in the pre- and post-genomic eras, new challenges for postgenomic research are discussed. Such discussion leads to a call to search for a new conceptual framework and holistic methodologies.

CONCLUSION

Throughout this synthesis, the genome theory of somatic cell evolution is highlighted in contrast to gene theory, which ignores the karyotype-mediated higher level of genetic information. Since "system inheritance" is defined by the genome context (gene content and genomic topology) while "parts inheritance" is defined by genes/epigenes, molecular cytogenetics and cytogenomics (which directly study genome structure, function, alteration and evolution) will play important roles in this postgenomic era.

Dynamics of Bacterial Gene Regulatory Networks

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

Re-thinking our understanding of immunity: Robustness in the tissue reconstruction system

Robustness, understood as the maintenance of specific functionalities of a given system against internal and external perturbations, is pervasive in today's biology. Yet precise applications of this notion to the immune system have been scarce. Here we show that the concept of robustness sheds light on tissue repair, and particularly on the crucial role the immune system plays in this process. We describe the specific mechanisms, including plasticity and redundancy, by which robustness is achieved in the tissue reconstruction system (TRS). In turn, tissue repair offers a very important test case for assessing the usefulness of the concept of robustness, and identifying different varieties of robustness.

Emergent Properties in Complex Synthetic Bacterial Promoters

Regulation of gene expression in bacteria results from the interplay between hundreds of transcriptional factors (TFs) at target promoters. However, how the arrangement of binding sites for TFs generates the regulatory logic of promoters is not well-known. Here, we generated and fully characterized a library of synthetic complex promoters for the global regulators, CRP and IHF, in Escherichia coli, which are formed by a weak -35/-10 consensus sequence preceded by four combinatorial binding sites for these two TFs. Using this approach, we found that while cis-elements for CRP preferentially activate promoters when located immediately upstream of the promoter consensus, binding sites for IHF mainly function as "UP" elements and stimulate transcription in several different architectures in the absence of this protein. However, the combination of CRP- and IHF-binding sites resulted in emergent properties in these complex promoters, where the activity of combinatorial promoters cannot be predicted from the individual behavior of its components. Taken together, the results presented here add to the information on architecture-logic of complex promoters in bacteria.

Systematic and synthetic approaches to rewire regulatory networks

Microbial gene regulatory networks are composed of cis- and trans-components that in concert act to control essential and adaptive cellular functions. Regulatory components and interactions evolve to adopt new configurations through mutations and network rewiring events, resulting in novel phenotypes that may benefit the cell. Advances in high-throughput DNA synthesis and sequencing have enabled the development of new tools and approaches to better characterize and perturb various elements of regulatory networks. Here, we highlight key recent approaches to systematically dissect the sequence space of cis-regulatory elements and trans-regulators as well as their inter-connections. These efforts yield fundamental insights into the architecture, robustness, and dynamics of gene regulation and provide models and design principles for building synthetic regulatory networks for a variety of practical applications.

Engineering microbial phenotypes through rewiring of genetic networks

The ability to program cellular behaviour is a major goal of synthetic biology, with applications in health, agriculture and chemicals production. Despite efforts to build 'orthogonal' systems, interactions between engineered genetic circuits and the endogenous regulatory network of a host cell can have a significant impact on desired functionality. We have developed a strategy to rewire the endogenous cellular regulatory network of yeast to enhance compatibility with synthetic protein and metabolite production. We found that introducing novel connections in the cellular regulatory network enabled us to increase the production of heterologous proteins and metabolites. This strategy is demonstrated in yeast strains that show significantly enhanced heterologous protein expression and higher titers of terpenoid production. Specifically, we found that the addition of transcriptional regulation between free radical induced signalling and nitrogen regulation provided robust improvement of protein production. Assessment of rewired networks revealed the importance of key topological features such as high betweenness centrality. The generation of rewired transcriptional networks, selection for specific phenotypes, and analysis of resulting library members is a powerful tool for engineering cellular behavior and may enable improved integration of heterologous protein and metabolite pathways.

The adaptive landscape of wildtype and glycosylation-deficient populations of the industrial yeast Pichia pastoris

Background

The effects of long-term environmental adaptation and the implications of major cellular malfunctions are still poorly understood for non-model but biotechnologically relevant species. In this study we performed a large-scale laboratory evolution experiment with 48 populations of the yeast Pichia pastoris in order to establish a general adaptive landscape upon long-term selection in several glucose-based growth environments. As a model for a cellular malfunction the implications of OCH1 mannosyltransferase knockout-mediated glycosylation-deficiency were analyzed.

Results

In-depth growth profiling of evolved populations revealed several instances of genotype-dependent growth trade-off/cross-benefit correlations in non-evolutionary growth conditions. On the genome level a high degree of mutational convergence was observed among independent populations. Environment-dependent mutational hotspots were related to osmotic stress-, Rim - and cAMP signaling pathways. In agreement with the observed growth phenotypes, our data also suggest diverging compensatory mutations in glycosylation-deficient populations. High osmolarity glycerol (HOG) pathway loss-of-functions mutations, including genes such as SSK2 and SSK4, represented a major adaptive strategy during environmental adaptation. However, genotype-specific HOG-related mutations were predominantly observed in opposing environmental conditions. Surprisingly, such mutations emerged during salt stress adaptation in OCH1 knockout populations and led to growth trade-offs in non-adaptive conditions that were distinct from wildtype HOG-mutants. Further environment-dependent mutations were identified for a hitherto uncharacterized species-specific Gal4-like transcriptional regulator involved in environmental sensing.

Conclusion

We show that metabolic constraints such as glycosylation-deficiency can contribute to evolution on the molecular level, even in non-diverging growth environments. Our dataset suggests universal adaptive mechanisms involving cellular stress response and cAMP/PKA signaling but also the existence of highly species-specific strategies involving unique transcriptional regulators, improving our biological understanding of distinct Ascomycetes species.

Electronic supplementary material

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

Synthetic Biology: Engineering Living Systems from Biophysical Principles

Synthetic biology was founded as a biophysical discipline that sought explanations for the origins of life from chemical and physical first principles. Modern synthetic biology has been reinvented as an engineering discipline to design new organisms as well as to better understand fundamental biological mechanisms. However, success is still largely limited to the laboratory and transformative applications of synthetic biology are still in their infancy. Here, we review six principles of living systems and how they compare and contrast with engineered systems. We cite specific examples from the synthetic biology literature that illustrate these principles and speculate on their implications for further study. To fully realize the promise of synthetic biology, we must be aware of life's unique properties.

CRISPR Perturbation of Gene Expression Alters Bacterial Fitness under Stress and Reveals Underlying Epistatic Constraints

The evolution of antibiotic resistance has engendered an impending global health crisis that necessitates a greater understanding of how resistance emerges. The impact of nongenetic factors and how they influence the evolution of resistance is a largely unexplored area of research. Here we present a novel application of CRISPR-Cas9 technology for investigating how gene expression governs the adaptive pathways available to bacteria during the evolution of resistance. We examine the impact of gene expression changes on bacterial adaptation by constructing a library of deactivated CRISPR-Cas9 synthetic devices to tune the expression of a set of stress-response genes in Escherichia coli. We show that artificially inducing perturbations in gene expression imparts significant synthetic control over fitness and growth during stress exposure. We present evidence that these impacts are reversible; strains with synthetically perturbed gene expression regained wild-type growth phenotypes upon stress removal, while maintaining divergent growth characteristics under stress. Furthermore, we demonstrate a prevailing trend toward negative epistatic interactions when multiple gene perturbations are combined simultaneously, thereby posing an intrinsic constraint on gene expression underlying adaptive trajectories. Together, these results emphasize how CRISPR-Cas9 can be employed to engineer gene expression changes that shape bacterial adaptation, and present a novel approach to synthetically control the evolution of antimicrobial resistance.

Engineering orthogonal dual transcription factors for multi-input synthetic promoters

Synthetic biology has seen an explosive growth in the capability of engineering artificial gene circuits from transcription factors (TFs), particularly in bacteria. However, most artificial networks still employ the same core set of TFs (for example LacI, TetR and cI). The TFs mostly function via repression and it is difficult to integrate multiple inputs in promoter logic. Here we present to our knowledge the first set of dual activator-repressor switches for orthogonal logic gates, based on bacteriophage λ cI variants and multi-input promoter architectures. Our toolkit contains 12 TFs, flexibly operating as activators, repressors, dual activator–repressors or dual repressor–repressors, on up to 270 synthetic promoters. To engineer non cross-reacting cI variants, we design a new M13 phagemid-based system for the directed evolution of biomolecules. Because cI is used in so many synthetic biology projects, the new set of variants will easily slot into the existing projects of other groups, greatly expanding current engineering capacities.

Genetic circuits usually employ the same set of transcription factors which can act via repression or activation of the target promoter. Here the authors present dual activator-repressor switches, designed via directed evolution, for orthogonal logic gates and multi-input circuit architectures.

A genome wide dosage suppressor network reveals genomic robustness

Genomic robustness is the extent to which an organism has evolved to withstand the effects of deleterious mutations. We explored the extent of genomic robustness in budding yeast by genome wide dosage suppressor analysis of 53 conditional lethal mutations in cell division cycle and RNA synthesis related genes, revealing 660 suppressor interactions of which 642 are novel. This collection has several distinctive features, including high co-occurrence of mutant-suppressor pairs within protein modules, highly correlated functions between the pairs and higher diversity of functions among the co-suppressors than previously observed. Dosage suppression of essential genes encoding RNA polymerase subunits and chromosome cohesion complex suggests a surprising degree of functional plasticity of macromolecular complexes, and the existence of numerous degenerate pathways for circumventing the effects of potentially lethal mutations. These results imply that organisms and cancer are likely able to exploit the genomic robustness properties, due the persistence of cryptic gene and pathway functions, to generate variation and adapt to selective pressures.

Production of selenium nanoparticles in Pseudomonas putida KT2440

Selenium (Se) is an essential element for the cell that has multiple applications in medicine and technology; microorganisms play an important role in Se transformations in the environment. Here we report the previously unidentified ability of the soil bacterium Pseudomonas putida KT2440 to synthesize nanoparticles of elemental selenium (nano-Se) from selenite. Our results show that P. putida is able to reduce selenite aerobically, but not selenate, to nano-Se. Kinetic analysis indicates that, in LB medium supplemented with selenite (1 mM), reduction to nano-Se occurs at a rate of 0.444 mmol L⁻¹ h⁻¹ beginning in the middle-exponential phase and with a final conversion yield of 89%. Measurements with a transmission electron microscope (TEM) show that nano-Se particles synthesized by P. putida have a size range of 100 to 500 nm and that they are located in the surrounding medium or bound to the cell membrane. Experiments involving dynamic light scattering (DLS) show that, in aqueous solution, recovered nano-Se particles have a size range of 70 to 360 nm. The rapid kinetics of conversion, easy retrieval of nano-Se and the metabolic versatility of P. putida offer the opportunity to use this model organism as a microbial factory for production of selenium nanoparticles.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

From noise to synthetic nucleoli: can synthetic biology achieve new insights?

Synthetic biology aims to re-organise and control biological components to make functional devices. Along the way, the iterative process of designing and testing gene circuits has the potential to yield many insights into the functioning of the underlying chassis of cells. Thus, synthetic biology is converging with disciplines such as systems biology and even classical cell biology, to give a new level of predictability to gene expression, cell metabolism and cellular signalling networks. This review gives an overview of the contributions that synthetic biology has made in understanding gene expression, in terms of cell heterogeneity (noise), the coupling of growth and energy usage to expression, and spatiotemporal considerations. We mainly compare progress in bacterial and mammalian systems, which have some of the most-developed engineering frameworks. Overall, one view of synthetic biology can be neatly summarised as "creating in order to understand."

High-throughput mathematical analysis identifies Turing networks for patterning with equally diffusing signals

The Turing reaction-diffusion model explains how identical cells can self-organize to form spatial patterns. It has been suggested that extracellular signaling molecules with different diffusion coefficients underlie this model, but the contribution of cell-autonomous signaling components is largely unknown. We developed an automated mathematical analysis to derive a catalog of realistic Turing networks. This analysis reveals that in the presence of cell-autonomous factors, networks can form a pattern with equally diffusing signals and even for any combination of diffusion coefficients. We provide a software (available at http://www.RDNets.com) to explore these networks and to constrain topologies with qualitative and quantitative experimental data. We use the software to examine the self-organizing networks that control embryonic axis specification and digit patterning. Finally, we demonstrate how existing synthetic circuits can be extended with additional feedbacks to form Turing reaction-diffusion systems. Our study offers a new theoretical framework to understand multicellular pattern formation and enables the wide-spread use of mathematical biology to engineer synthetic patterning systems.

The propagation of perturbations in rewired bacterial gene networks

What happens to gene expression when you add new links to a gene regulatory network? To answer this question, we profile 85 network rewirings in E. coli. Here we report that concerted patterns of differential expression propagate from reconnected hub genes. The rewirings link promoter regions to different transcription factor and σ-factor genes, resulting in perturbations that span four orders of magnitude, changing up to ∼ 70% of the transcriptome. Importantly, factor connectivity and promoter activity both associate with perturbation size. Perturbations from related rewirings have more similar transcription profiles and a statistical analysis reveals ∼ 20 underlying states of the system, associating particular gene groups with rewiring constructs. We examine two large clusters (ribosomal and flagellar genes) in detail. These represent alternative global outcomes from different rewirings because of antagonism between these major cell states. This data set of systematically related perturbations enables reverse engineering and discovery of underlying network interactions.

Re-engineering cellular physiology by rewiring high-level global regulatory genes

Knowledge of global regulatory networks has been exploited to rewire the gene control programmes of the model bacterium Salmonella enterica serovar Typhimurium. The product is an organism with competitive fitness that is superior to that of the wild type but tuneable under specific growth conditions. The paralogous hns and stpA global regulatory genes are located in distinct regions of the chromosome and control hundreds of target genes, many of which contribute to stress resistance. The locations of the hns and stpA open reading frames were exchanged reciprocally, each acquiring the transcription control signals of the other. The new strain had none of the compensatory mutations normally associated with alterations to hns expression in Salmonella; instead it displayed rescheduled expression of the stress and stationary phase sigma factor RpoS and its regulon. Thus the expression patterns of global regulators can be adjusted artificially to manipulate microbial physiology, creating a new and resilient organism.