A unified resource for transcriptional regulation in Escherichia coli K-12 incorporating high-throughput-generated binding data into RegulonDB version 10.0

Irreversibility in bacterial regulatory networks

Irreversibility, in which a transient perturbation leaves a system in a new state, is an emergent property in systems of interacting entities. This property has well-established implications in statistical physics but remains underexplored in biological networks, especially for bacteria and other prokaryotes whose regulation of gene expression occurs predominantly at the transcriptional level. Focusing on the reconstructed regulatory network of Escherichia coli, we examine network responses to transient single-gene perturbations. We predict irreversibility in numerous cases and find that the incidence of irreversibility increases with the proximity of the perturbed gene to positive circuits in the network. Comparison with experimental data suggests a connection between the predicted irreversibility to transient perturbations and the evolutionary response to permanent perturbations.

Pervasive gene deregulation underlies adaptation and maladaptation in trimethoprim-resistant E. coli

IMPORTANCE

Bacteria employ a number of mechanisms to adapt to antibiotics. Mutations in transcriptional regulators alter the expression levels of genes that can change the susceptibility of bacteria to antibiotics. Two-component signaling proteins are a major class of signaling molecule used by bacteria to regulate transcription. In previous work, we found that mutations in MgrB, a feedback regulator of the PhoQP two-component system, conferred trimethoprim tolerance to Escherichia coli. Here, we elucidate how mutations in MgrB have a domino-like effect on the gene regulatory network of E. coli. As a result, pervasive perturbation of gene regulation ensues. Depending on the environmental context, this pervasive deregulation is either adaptive or maladaptive. Our study sheds light on how deregulation of gene expression can be beneficial for bacteria when challenged with antibiotics, and why regulators like MgrB may have evolved in the first place.

Systems biology of competency in Vibrio natriegens is revealed by applying novel data analytics to the transcriptome

Vibrio natriegens regulates natural competence through the TfoX and QstR transcription factors, which are involved in external DNA capture and transport. However, the extensive genetic and transcriptional regulatory basis for competency remains unknown. We used a machine-learning approach to decompose Vibrio natriegens's transcriptome into 45 groups of independently modulated sets of genes (iModulons). Our findings show that competency is associated with the repression of two housekeeping iModulons (iron metabolism and translation) and the activation of six iModulons; including TfoX and QstR, a novel iModulon of unknown function, and three housekeeping iModulons (representing motility, polycations, and reactive oxygen species [ROS] responses). Phenotypic screening of 83 gene deletion strains demonstrates that loss of iModulon function reduces or eliminates competency. This database-iModulon-discovery cycle unveils the transcriptomic basis for competency and its relationship to housekeeping functions. These results provide the genetic basis for systems biology of competency in this organism.

Strengthening microbial cell factories for efficient production of bioactive molecules

High demand of bioactive molecules (food additives, antibiotics, plant growth enhancers, cosmetics, pigments and other commercial products) is the prime need for the betterment of human life where the applicability of the synthetic chemical product is on the saturation due to associated toxicity and ornamentations. It has been noticed that the discovery and productivity of such molecules in natural scenarios are limited due to low cellular yields as well as less optimized conventional methods. In this respect, microbial cell factories timely fulfilling the requirement of synthesizing bioactive molecules by improving production yield and screening more promising structural homologues of the native molecule. Where the robustness of the microbial host can be potentially achieved by taking advantage of cell engineering approaches such as tuning functional and adjustable factors, metabolic balancing, adapting cellular transcription machinery, applying high throughput OMICs tools, stability of genotype/phenotype, organelle optimizations, genome editing (CRISPER/Cas mediated system) and also by developing accurate model systems via machine-learning tools. In this article, we provide an overview from traditional to recent trends and the application of newly developed technologies, for strengthening the systemic approaches and providing future directions for enhancing the robustness of microbial cell factories to speed up the production of biomolecules for commercial purposes.

Exploring the mono-/bistability range of positively autoregulated signaling systems in the presence of competing transcription factor binding sites

Binding of transcription factor (TF) proteins to regulatory DNA sites is key to accurate control of gene expression in response to environmental stimuli. Theoretical modeling of transcription regulation is often focused on a limited set of genes of interest, while binding of the TF to other genomic sites is seldom considered. The total number of TF binding sites (TFBSs) affects the availability of TF protein molecules and sequestration of a TF by TFBSs can promote bistability. For many signaling systems where a graded response is desirable for continuous control over the input range, biochemical parameters of the regulatory proteins need be tuned to avoid bistability. Here we analyze the mono-/bistable parameter range for positively autoregulated two-component systems (TCSs) in the presence of different numbers of competing TFBSs. TCS signaling, one of the major bacterial signaling strategies, couples signal perception with output responses via protein phosphorylation. For bistability, competition for TF proteins by TFBSs lowers the requirement for high fold change of the autoregulated transcription but demands high phosphorylation activities of TCS proteins. We show that bistability can be avoided with a low phosphorylation capacity of TCSs, a high TF affinity for the autoregulated promoter or a low fold change in signaling protein levels upon induction. These may represent general design rules for TCSs to ensure uniform graded responses. Examining the mono-/bistability parameter range allows qualitative prediction of steady-state responses, which are experimentally validated in the E. coli CusRS system.

RegulonDB 11.0: Comprehensive high-throughput datasets on transcriptional regulation in Escherichia coli K-12

Genomics has set the basis for a variety of methodologies that produce high-throughput datasets identifying the different players that define gene regulation, particularly regulation of transcription initiation and operon organization. These datasets are available in public repositories, such as the Gene Expression Omnibus, or ArrayExpress. However, accessing and navigating such a wealth of data is not straightforward. No resource currently exists that offers all available high and low-throughput data on transcriptional regulation in Escherichia coli K-12 to easily use both as whole datasets, or as individual interactions and regulatory elements. RegulonDB (https://regulondb.ccg.unam.mx) began gathering high-throughput dataset collections in 2009, starting with transcription start sites, then adding ChIP-seq and gSELEX in 2012, with up to 99 different experimental high-throughput datasets available in 2019. In this paper we present a radical upgrade to more than 2000 high-throughput datasets, processed to facilitate their comparison, introducing up-to-date collections of transcription termination sites, transcription units, as well as transcription factor binding interactions derived from ChIP-seq, ChIP-exo, gSELEX and DAP-seq experiments, besides expression profiles derived from RNA-seq experiments. For ChIP-seq experiments we offer both the data as presented by the authors, as well as data uniformly processed in-house, enhancing their comparability, as well as the traceability of the methods and reproducibility of the results. Furthermore, we have expanded the tools available for browsing and visualization across and within datasets. We include comparisons against previously existing knowledge in RegulonDB from classic experiments, a nucleotide-resolution genome viewer, and an interface that enables users to browse datasets by querying their metadata. A particular effort was made to automatically extract detailed experimental growth conditions by implementing an assisted curation strategy applying Natural language processing and machine learning. We provide summaries with the total number of interactions found in each experiment, as well as tools to identify common results among different experiments. This is a long-awaited resource to make use of such wealth of knowledge and advance our understanding of the biology of the model bacterium E. coli K-12.

Sensory Systems and Transcriptional Regulation in Escherichia coli

In free-living bacteria, the ability to regulate gene expression is at the core of adapting and interacting with the environment. For these systems to have a logic, a signal must trigger a genetic change that helps the cell to deal with what implies its presence in the environment; briefly, the response is expected to include a feedback to the signal. Thus, it makes sense to think of genetic sensory mechanisms of gene regulation. Escherichia coli K-12 is the bacterium model for which the largest number of regulatory systems and its sensing capabilities have been studied in detail at the molecular level. In this special issue focused on biomolecular sensing systems, we offer an overview of the transcriptional regulatory corpus of knowledge for E. coli that has been gathered in our database, RegulonDB, from the perspective of sensing regulatory systems. Thus, we start with the beginning of the information flux, which is the signal's chemical or physical elements detected by the cell as changes in the environment; these signals are internally transduced to transcription factors and alter their conformation. Signals transduced to effectors bind allosterically to transcription factors, and this defines the dominant sensing mechanism in E. coli. We offer an updated list of the repertoire of known allosteric effectors, as well as a list of the currently known different mechanisms of this sensing capability. Our previous definition of elementary genetic sensory-response units, GENSOR units for short, that integrate signals, transport, gene regulation, and the biochemical response of the regulated gene products of a given transcriptional factor fit perfectly with the purpose of this overview. We summarize the functional heterogeneity of their response, based on our updated collection of GENSORs, and we use them to identify the expected feedback as part of their response. Finally, we address the question of multiple sensing in the regulatory network of E. coli. This overview introduces the architecture of sensing and regulation of native components in E.coli K-12, which might be a source of inspiration to bioengineering applications.

proChIPdb: a chromatin immunoprecipitation database for prokaryotic organisms

The transcriptional regulatory network in prokaryotes controls global gene expression mostly through transcription factors (TFs), which are DNA-binding proteins. Chromatin immunoprecipitation (ChIP) with DNA sequencing methods can identify TF binding sites across the genome, providing a bottom-up, mechanistic understanding of how gene expression is regulated. ChIP provides indispensable evidence toward the goal of acquiring a comprehensive understanding of cellular adaptation and regulation, including condition-specificity. ChIP-derived data's importance and labor-intensiveness motivate its broad dissemination and reuse, which is currently an unmet need in the prokaryotic domain. To fill this gap, we present proChIPdb (prochipdb.org), an information-rich, interactive web database. This website collects public ChIP-seq/-exo data across several prokaryotes and presents them in dashboards that include curated binding sites, nucleotide-resolution genome viewers, and summary plots such as motif enrichment sequence logos. Users can search for TFs of interest or their target genes, download all data, dashboards, and visuals, and follow external links to understand regulons through biological databases and the literature. This initial release of proChIPdb covers diverse organisms, including most major TFs of Escherichia coli, and can be expanded to support regulon discovery across the prokaryotic domain.

A balancing act in transcription regulation by response regulators: titration of transcription factor activity by decoy DNA binding sites

Studies of transcription regulation are often focused on binding of transcription factors (TFs) to a small number of promoters of interest. It is often assumed that TFs are in great excess to their binding sites (TFBSs) and competition for TFs between DNA sites is seldom considered. With increasing evidence that TFBSs are exceedingly abundant for many TFs and significant variations in TF and TFBS numbers occur during growth, the interplay between a TF and all TFBSs should not be ignored. Here, we use additional decoy DNA sites to quantitatively analyze how the relative abundance of a TF to its TFBSs impacts the steady-state level and onset time of gene expression for the auto-activated Escherichia coli PhoB response regulator. We show that increasing numbers of decoy sites progressively delayed transcription activation and lowered promoter activities. Perturbation of transcription regulation by additional TFBSs did not require extreme numbers of decoys, suggesting that PhoB is approximately at capacity for its DNA sites. Addition of decoys also converted a graded response to a bi-modal response. We developed a binding competition model that captures the major features of experimental observations, providing a quantitative framework to assess how variations in TFs and TFBSs influence transcriptional responses.

Responses of Escherichia coli and Listeria monocytogenes to ozone treatment on non-host tomato: Efficacy of intervention and evidence of induced acclimation

Because of the continuous rise of foodborne illnesses caused by the consumption of raw fruits and vegetables, effective post-harvest anti-microbial strategies are necessary. The aim of this study was to evaluate the anti-microbial efficacy of ozone (O₃) against two common causes of fresh produce contamination, the Gram-negative Escherichia coli O157:H7 and Gram-positive Listeria monocytogenes, and to relate its effects to potential mechanisms of xenobiosis by transcriptional network modeling. The study on non-host tomato environment correlated the dose × time aspects of xenobiosis by examining the correlation between bacterial survival in terms of log-reduction and defense responses at the level of gene expression. In E. coli, low (1 μg O₃/g of fruit) and moderate (2 μg O₃/g of fruit) doses caused insignificant reduction in survival, while high dose (3 μg/g of fruit) caused significant reduction in survival in a time-dependent manner. In L. monocytogenes, moderate dose caused significant reduction even with short-duration exposure. Distinct responses to O₃ xenobiosis between E. coli and L. monocytogenes are likely related to differences in membrane and cytoplasmic structure and components. Transcriptome profiling by RNA-Seq showed that primary defenses in E. coli were attenuated after exposure to a low dose, while the responses at moderate dose were characterized by massive upregulation of pathogenesis and stress-related genes, which implied the activation of defense responses. More genes were downregulated during the first hour at high dose, with a large number of such genes getting significantly upregulated after 2 hr and 3 hr. This trend suggests that prolonged exposure led to potential adaptation. In contrast, massive downregulation of genes was observed in L. monocytogenes regardless of dose and exposure duration, implying a mechanism of defense distinct from that of E. coli. The nature of bacterial responses revealed by this study should guide the selection of xenobiotic agents for eliminating bacterial contamination on fresh produce without overlooking the potential risks of adaptation.

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.

Direct Regulons of AtxA, the Master Virulence Regulator of Bacillus anthracis

AtxA, the master virulence regulator of Bacillus anthracis, regulates the expression of three toxins and genes for capsule formation that are required for the pathogenicity of B. anthracis. Recent transcriptome analyses showed that AtxA affects a large number of genes on the chromosome and plasmids, suggesting a role as a global regulator. However, information on genes directly regulated by AtxA is scarce. In this work, we conducted genome-wide analyses and cataloged the binding sites of AtxA in vivo and transcription start sites on the B. anthracis genome. By integrating these results, we detected eight genes as direct regulons of AtxA. These consisted of five protein-coding genes, including two of the three toxin genes, and three genes encoding the small RNAs XrrA and XrrB and a newly discovered 95-nucleotide small RNA, XrrC. Transcriptomes from single-knockout mutants of these small RNAs revealed changes in the transcription levels of genes related to the aerobic electron transport chain, heme biosynthesis, and amino acid metabolism, suggesting their function for the control of cell physiology. These results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide a data set for the further study of the genomics and genetics of B. anthracis. IMPORTANCE Bacillus anthracis is the Gram-positive bacterial species that causes anthrax. Anthrax is still prevalent in countries mainly in Asia and Africa, where it causes economic damage and remains a public health issue. The mechanism of pathogenicity is mainly explained by the three toxin proteins expressed from the pXO1 plasmid and by proteins involved in capsule formation expressed from the pXO2 plasmid. AtxA is a protein expressed from the pXO1 plasmid that is known to upregulate genes involved in toxin production and capsule formation and is thus considered the master virulence regulator of B. anthracis. Therefore, understanding the detailed mechanism of gene regulation is important for the control of anthrax. The significance of this work lies in the identification of genes that are directly regulated by AtxA via genome-wide analyses. The results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide useful resources for a further understanding of B. anthracis.

LuxAB-Based Microbial Cell Factories for the Sensing, Manufacturing and Transformation of Industrial Aldehydes

The application of genetically encoded biosensors enables the detection of small molecules in living cells and has facilitated the characterization of enzymes, their directed evolution and the engineering of (natural) metabolic pathways. In this work, the LuxAB biosensor system from Photorhabdus luminescens was implemented in Escherichia coli to monitor the enzymatic production of aldehydes from primary alcohols and carboxylic acid substrates. A simple high-throughput assay utilized the bacterial luciferase—previously reported to only accept aliphatic long-chain aldehydes—to detect structurally diverse aldehydes, including aromatic and monoterpene aldehydes. LuxAB was used to screen the substrate scopes of three prokaryotic oxidoreductases: an alcohol dehydrogenase (Pseudomonas putida), a choline oxidase variant (Arthrobacter chlorophenolicus) and a carboxylic acid reductase (Mycobacterium marinum). Consequently, high-value aldehydes such as cinnamaldehyde, citral and citronellal could be produced in vivo in up to 80% yield. Furthermore, the dual role of LuxAB as sensor and monooxygenase, emitting bioluminescence through the oxidation of aldehydes to the corresponding carboxylates, promises implementation in artificial enzyme cascades for the synthesis of carboxylic acids. These findings advance the bio-based detection, preparation and transformation of industrially important aldehydes in living cells.

The EcoCyc Database in 2021

The EcoCyc model-organism database collects and summarizes experimental data for Escherichia coli K-12. EcoCyc is regularly updated by the manual curation of individual database entries, such as genes, proteins, and metabolic pathways, and by the programmatic addition of results from select high-throughput analyses. Updates to the Pathway Tools software that supports EcoCyc and to the web interface that enables user access have continuously improved its usability and expanded its functionality. This article highlights recent improvements to the curated data in the areas of metabolism, transport, DNA repair, and regulation of gene expression. New and revised data analysis and visualization tools include an interactive metabolic network explorer, a circular genome viewer, and various improvements to the speed and usability of existing tools.

Recent trends in biocatalysis

Biocatalysis has undergone revolutionary progress in the past century. Benefited by the integration of multidisciplinary technologies, natural enzymatic reactions are constantly being explored. Protein engineering gives birth to robust biocatalysts that are widely used in industrial production. These research achievements have gradually constructed a network containing natural enzymatic synthesis pathways and artificially designed enzymatic cascades. Nowadays, the development of artificial intelligence, automation, and ultra-high-throughput technology provides infinite possibilities for the discovery of novel enzymes, enzymatic mechanisms and enzymatic cascades, and gradually complements the lack of remaining key steps in the pathway design of enzymatic total synthesis. Therefore, the research of biocatalysis is gradually moving towards the era of novel technology integration, intelligent manufacturing and enzymatic total synthesis.

Widespread divergent transcription from bacterial and archaeal promoters is a consequence of DNA-sequence symmetry

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Independent component analysis recovers consistent regulatory signals from disparate datasets

The availability of bacterial transcriptomes has dramatically increased in recent years. This data deluge could result in detailed inference of underlying regulatory networks, but the diversity of experimental platforms and protocols introduces critical biases that could hinder scalable analysis of existing data. Here, we show that the underlying structure of the E. coli transcriptome, as determined by Independent Component Analysis (ICA), is conserved across multiple independent datasets, including both RNA-seq and microarray datasets. We subsequently combined five transcriptomics datasets into a large compendium containing over 800 expression profiles and discovered that its underlying ICA-based structure was still comparable to that of the individual datasets. With this understanding, we expanded our analysis to over 3,000 E. coli expression profiles and predicted three high-impact regulons that respond to oxidative stress, anaerobiosis, and antibiotic treatment. ICA thus enables deep analysis of disparate data to uncover new insights that were not visible in the individual datasets.

Cells adapt to diverse environments by regulating gene expression. Genome-wide measurements of gene expression levels have exponentially increased in recent years, but successful integration and analysis of these datasets are limited. Recently, we showed that independent component analysis (ICA), a signal deconvolution algorithm, can separate a large bacterial gene expression dataset into groups of co-regulated genes. This previous study focused on data generated by a standardized pipeline and did not address whether ICA extracts the same quantitative co-expression signals across expression profiling platforms. In this study, we show that ICA finds similar co-regulation patterns underlying multiple gene expression datasets and can be used as a tool to integrate and interpret diverse datasets. Using a dataset containing over 3,000 expression profiles, we predicted three new regulons and characterized their activities. Since large, standardized expression datasets only exist for a few bacterial strains, these results broaden the possible applications of this tool to better understand transcriptional regulation across a wide range of microbes.

How Retroactivity Affects the Behavior of Incoherent Feedforward Loops

An incoherent feedforward loop (IFFL) is a network motif known for its ability to accelerate responses and generate pulses. It remains an open question to understand the behavior of IFFLs in contexts with high levels of retroactivity, where an upstream transcription factor binds to numerous downstream binding sites. Here we study the behavior of IFFLs by simulating and comparing ODE models with different levels of retroactivity. We find that increasing retroactivity in an IFFL can increase, decrease, or keep the network's response time and pulse amplitude constant. This suggests that increasing retroactivity, traditionally considered an impediment to designing robust synthetic systems, could be exploited to improve the performance of IFFLs. In contrast, we find that increasing retroactivity in a negative autoregulated circuit can only slow the response. The ability of an IFFL to flexibly handle retroactivity may have contributed to its significant abundance in both bacterial and eukaryotic regulatory networks.

Elucidation of Regulatory Modes for Five Two-Component Systems in Escherichia coli Reveals Novel Relationships

Escherichia coli uses two-component systems (TCSs) to respond to environmental signals. TCSs affect gene expression and are parts of E. coli's global transcriptional regulatory network (TRN). Here, we identified the regulons of five TCSs in E. coli MG1655: BaeSR and CpxAR, which were stimulated by ethanol stress; KdpDE and PhoRB, induced by limiting potassium and phosphate, respectively; and ZraSR, stimulated by zinc. We analyzed RNA-seq data using independent component analysis (ICA). ChIP-exo data were used to validate condition-specific target gene binding sites. Based on these data, we do the following: (i) identify the target genes for each TCS; (ii) show how the target genes are transcribed in response to stimulus; and (iii) reveal novel relationships between TCSs, which indicate noncognate inducers for various response regulators, such as BaeR to iron starvation, CpxR to phosphate limitation, and PhoB and ZraR to cell envelope stress. Our understanding of the TRN in E. coli is thus notably expanded.IMPORTANCE E. coli is a common commensal microbe found in the human gut microenvironment; however, some strains cause diseases like diarrhea, urinary tract infections, and meningitis. E. coli's two-component systems (TCSs) modulate target gene expression, especially related to virulence, pathogenesis, and antimicrobial peptides, in response to environmental stimuli. Thus, it is of utmost importance to understand the transcriptional regulation of TCSs to infer bacterial environmental adaptation and disease pathogenicity. Utilizing a combinatorial approach integrating RNA sequencing (RNA-seq), independent component analysis, chromatin immunoprecipitation coupled with exonuclease treatment (ChIP-exo), and data mining, we suggest five different modes of TCS transcriptional regulation. Our data further highlight noncognate inducers of TCSs, which emphasizes the cross-regulatory nature of TCSs in E. coli and suggests that TCSs may have a role beyond their cognate functionalities. In summary, these results can lead to an understanding of the metabolic capabilities of bacteria and correctly predict complex phenotype under diverse conditions, especially when further incorporated with genome-scale metabolic models.

Redefining fundamental concepts of transcription initiation in bacteria

Despite enormous progress in understanding the fundamentals of bacterial gene regulation, our knowledge remains limited when compared with the number of bacterial genomes and regulatory systems to be discovered. Derived from a small number of initial studies, classic definitions for concepts of gene regulation have evolved as the number of characterized promoters has increased. Together with discoveries made using new technologies, this knowledge has led to revised generalizations and principles. In this Expert Recommendation, we suggest precise, updated definitions that support a logical, consistent conceptual framework of bacterial gene regulation, focusing on transcription initiation. The resulting concepts can be formalized by ontologies for computational modelling, laying the foundation for improved bioinformatics tools, knowledge-based resources and scientific communication. Thus, this work will help researchers construct better predictive models, with different formalisms, that will be useful in engineering, synthetic biology, microbiology and genetics.

A review of methods for the reconstruction and analysis of integrated genome-scale models of metabolism and regulation

The current survey aims to describe the main methodologies for extending the reconstruction and analysis of genome-scale metabolic models and phenotype simulation with Flux Balance Analysis mathematical frameworks, via the integration of Transcriptional Regulatory Networks and/or gene expression data. Although the surveyed methods are aimed at improving phenotype simulations obtained from these models, the perspective of reconstructing integrated genome-scale models of metabolism and gene expression for diverse prokaryotes is still an open challenge.

Xenobiotic Effects of Chlorine Dioxide to Escherichia coli O157:H7 on Non-host Tomato Environment Revealed by Transcriptional Network Modeling: Implications to Adaptation and Selection

Escherichia coli serotype O157:H7 is one of the major agents of pathogen outbreaks associated with fresh fruits and vegetables. Gaseous chlorine dioxide (ClO2) has been reported to be an effective intervention to eliminate bacterial contamination on fresh produce. Although remarkable positive effects of low doses of ClO2 have been reported, the genetic regulatory machinery coordinating the mechanisms of xenobiotic effects and the potential bacterial adaptation remained unclear. This study examined the temporal transcriptome profiles of E. coli O157:H7 during exposure to different doses of ClO2 in order to elucidate the genetic mechanisms underlying bacterial survival under such harsh conditions. Dosages of 1 μg, 5 μg, and 10 μg ClO2 per gram of tomato fruits cause different effects with dose-by-time dynamics. The first hour of exposure to 1 μg and 5 μg ClO2 caused only partial killing with significant growth reduction starting at the second hour, and without further significant reduction at the third hour. However, 10 μg ClO2 exposure led to massive bacterial cell death at 1 h with further increase in cell death at 2 and 3 h. The first hour exposure to 1 μg ClO2 caused activation of primary defense and survival mechanisms. However, the defense response was attenuated during the second and third hours. Upon treatment with 5 μg ClO2, the transcriptional networks showed massive downregulation of pathogenesis and stress response genes at the first hour of exposure, with decreasing number of differentially expressed genes at the second and third hours. In contrast, more genes were further downregulated with exposure to 10 μg ClO2 at the first hour, with the number of both upregulated and downregulated genes significantly decreasing at the second hour. A total of 810 genes were uniquely upregulated at the third hour at 10 μg ClO2, suggesting that the potency of xenobiotic effects had led to potential adaptation. This study provides important knowledge on the possible selection of target molecules for eliminating bacterial contamination on fresh produce without overlooking potential risks of adaptation.

RegulonDB v 10.5: tackling challenges to unify classic and high throughput knowledge of gene regulation in E. coli K-12

RegulonDB, first published 20 years ago, is a comprehensive electronic resource about regulation of transcription initiation of Escherichia coli K-12 with decades of knowledge from classic molecular biology experiments, and recently also from high-throughput genomic methodologies. We curated the literature to keep RegulonDB up to date, and initiated curation of ChIP and gSELEX experiments. We estimate that current knowledge describes between 10% and 30% of the expected total number of transcription factor- gene regulatory interactions in E. coli. RegulonDB provides datasets for interactions for which there is no evidence that they affect expression, as well as expression datasets. We developed a proof of concept pipeline to merge binding and expression evidence to identify regulatory interactions. These datasets can be visualized in the RegulonDB JBrowse. We developed the Microbial Conditions Ontology with a controlled vocabulary for the minimal properties to reproduce an experiment, which contributes to integrate data from high throughput and classic literature. At a higher level of integration, we report Genetic Sensory-Response Units for 200 transcription factors, including their regulation at the metabolic level, and include summaries for 70 of them. Finally, we summarize our research with Natural language processing strategies to enhance our biocuration work.

Genetic Biosensor Design for Natural Product Biosynthesis in Microorganisms

Low yield and low titer of natural products are common issues in natural product biosynthesis through microbial cell factories. One effective way to resolve such bottlenecks is to design genetic biosensors to monitor and regulate the biosynthesis of target natural products. In this review, we evaluate the most recent advances in the design of genetic biosensors for natural product biosynthesis in microorganisms. In particular, we examine strategies for selection of genetic parts and construction principles for the design and evaluation of genetic biosensors. We also review the latest applications of transcription factor- and riboswitch-based genetic biosensors in natural product biosynthesis. Lastly, we discuss challenges and solutions in designing genetic biosensors for the biosynthesis of natural products in microorganisms.

The Escherichia coli transcriptome mostly consists of independently regulated modules

Underlying cellular responses is a transcriptional regulatory network (TRN) that modulates gene expression. A useful description of the TRN would decompose the transcriptome into targeted effects of individual transcriptional regulators. Here, we apply unsupervised machine learning to a diverse compendium of over 250 high-quality Escherichia coli RNA-seq datasets to identify 92 statistically independent signals that modulate the expression of specific gene sets. We show that 61 of these transcriptomic signals represent the effects of currently characterized transcriptional regulators. Condition-specific activation of signals is validated by exposure of E. coli to new environmental conditions. The resulting decomposition of the transcriptome provides: a mechanistic, systems-level, network-based explanation of responses to environmental and genetic perturbations; a guide to gene and regulator function discovery; and a basis for characterizing transcriptomic differences in multiple strains. Taken together, our results show that signal summation describes the composition of a model prokaryotic transcriptome.

Mechanistic insight into the regulation of transcriptional modules remains scarce. Here, the authors identify statistically independent gene sets by applying independent component analysis to a high-quality E. coli RNA-seq data compendium and find that most gene sets represent the effects of specific transcriptional regulators.

The ecology and evolution of microbial metabolic strategies

Free-living microbes are generally capable of growing on multiple different nutrients. Some of those nutrients are used simultaneously, while others are used sequentially. The pattern of nutrient preferences and co-utilization defines the metabolic strategy of a microorganism. Metabolic strategies can substantially affect ecological interactions between species, but their evolution and distribution across the tree of life remain poorly characterized. We discuss how the confluence of better computational models of genotype-phenotype maps and high-throughput experimental tools can help us fill gaps in our knowledge and incorporate metabolic strategies into quantitative predictive models of microbial consortia.

Massively Parallel Assays and Quantitative Sequence-Function Relationships

Over the last decade, a rich variety of massively parallel assays have revolutionized our understanding of how biological sequences encode quantitative molecular phenotypes. These assays include deep mutational scanning, high-throughput SELEX, and massively parallel reporter assays. Here, we review these experimental methods and how the data they produce can be used to quantitatively model sequence-function relationships. In doing so, we touch on a diverse range of topics, including the identification of clinically relevant genomic variants, the modeling of transcription factor binding to DNA, the functional and evolutionary landscapes of proteins, and cis-regulatory mechanisms in both transcription and mRNA splicing. We further describe a unified conceptual framework and a core set of mathematical modeling strategies that studies in these diverse areas can make use of. Finally, we highlight key aspects of experimental design and mathematical modeling that are important for the results of such studies to be interpretable and reproducible.

Bacterial Transcription Factors: Regulation by Pick "N" Mix

Transcription in most bacteria is tightly regulated in order to facilitate bacterial adaptation to different environments, and transcription factors play a key role in this. Here we give a brief overview of the essential features of bacterial transcription factors and how they affect transcript initiation at target promoters. We focus on complex promoters that are regulated by combinations of activators and repressors, combinations of repressors only, or combinations of activators. At some promoters, transcript initiation is regulated by nucleoid-associated proteins, which often work together with transcription factors. We argue that the distinction between nucleoid-associated proteins and transcription factors is blurred and that they likely share common origins.

A computational system for identifying operons based on RNA-seq data

An operon is a set of neighboring genes in a genome that is transcribed as a single polycistronic message. Genes that are part of the same operon often have related functional roles or participate in the same metabolic pathways. The majority of all bacterial genes are co-transcribed with one or more other genes as part of a multi-gene operon. Thus, accurate identification of operons is important in understanding co-regulation of genes and their functional relationships. Here, we present a computational system that uses RNA-seq data to determine operons throughout a genome. The system takes the name of a genome and one or more files of RNA-seq data as input. Our method combines primary genomic sequence information with expression data from the RNA-seq files in a unified probabilistic model in order to identify operons. We assess our method's ability to accurately identify operons in a range of species through comparison to external databases of operons, both experimentally confirmed and computationally predicted, and through focused experiments that confirm new operons identified by our method. Our system is freely available at https://cs.wellesley.edu/~btjaden/Rockhopper/.

Transcription activation in bacteria: ancient and modern

Regulatory interactions at the lac promoter.Activation of the transcription of genes is central to many processes of adaptation and differentiation in bacteria. Here, I review the molecular mechanisms by which transcription factors can activate the initiation of specific transcripts at bacterial promoters. The story is presented in the context of Marjory Stephenson's pioneering work on enzymatic adaptation in bacteria, and sets the different mechanisms in the greater context of how transcription regulatory mechanisms evolved.

MCO: towards an ontology and unified vocabulary for a framework-based annotation of microbial growth conditions

MOTIVATION

A major component in increasing our understanding of the biology of an organism is the mapping of its genotypic potential into its phenotypic expression profiles. This mapping is executed by the machinery of gene regulation, which is essentially studied by changes in growth conditions. Although many efforts have been made to systematize the annotation of experimental conditions in microbiology, the available annotations are not based on a consistent and controlled vocabulary, making difficult the identification of biologically meaningful comparisons of knowledge derived from different experiments or laboratories.

RESULTS

We curated terms related to experimental conditions that affect gene expression in Escherichia coli K-12. Since this is the best-studied microorganism, the collected terms are the seed for the Microbial Conditions Ontology (MCO), a controlled and structured vocabulary that can be expanded to annotate microbial conditions in general. Moreover, we developed an annotation framework to describe experimental conditions, providing the foundation to identify regulatory networks that operate under particular conditions.

AVAILABILITY AND IMPLEMENTATION

As far as we know, MCO is the first ontology for growth conditions of any bacterial organism, and it is available at http://regulondb.ccg.unam.mx and https://github.com/microbial-conditions-ontology. Furthermore, we will disseminate MCO throughout the Open Biological and Biomedical Ontology (OBO) Foundry in order to set a standard for the annotation of gene expression data. This will enable comparison of data from diverse data sources.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.