Rapid Sampling of Escherichia coli After Changing Oxygen Conditions Reveals Transcriptional Dynamics

Combined Transcriptomic and Proteomic Profiling of E. coli under Microaerobic versus Aerobic Conditions: The Multifaceted Roles of Noncoding Small RNAs and Oxygen-Dependent Sensing in Global Gene Expression Control

Adaptive mechanisms that facilitate intestinal colonization by the human microbiota, including Escherichia coli, may be better understood by analyzing the physiology and gene expression of bacteria in low-oxygen environments. We used high-throughput transcriptomics and proteomics to compare the expression profiles of E. coli grown under aerobic versus microaerobic conditions. Clustering of high-abundance transcripts under microaerobiosis highlighted genes controlling acid-stress adaptation (gadAXW, gadAB, hdeAB-yhiD and hdeD operons), cell adhesion/biofilm formation (pgaABCD and csgDEFG operons), electron transport (cydAB), oligopeptide transport (oppABCDF), and anaerobic respiration/fermentation (hyaABCDEF and hycABCDEFGHI operons). In contrast, downregulated genes were involved in iron transport (fhuABCD, feoABC and fepA-entD operons), iron-sulfur cluster assembly (iscRSUA and sufABCDSE operons), aerobic respiration (sdhDAB and sucABCDSE operons), and de novo nucleotide synthesis (nrdHIEF). Additionally, quantitative proteomics showed that the products (proteins) of these high- or low-abundance transcripts were expressed consistently. Our findings highlight interrelationships among energy production, carbon metabolism, and iron homeostasis. Moreover, we have identified and validated a subset of differentially expressed noncoding small RNAs (i.e., CsrC, RyhB, RprA and GcvB), and we discuss their regulatory functions during microaerobic growth. Collectively, we reveal key changes in gene expression at the transcriptional and post-transcriptional levels that sustain E. coli growth when oxygen levels are low.

Euler-Lagrangian Simulations: A Proper Tool for Predicting Cellular Performance in Industrial Scale Bioreactors

Eulerian-Lagrangian approach to investigate cellular responses in a bioreactor has become the center of attention in recent years. It was introduced to biotechnological processes about two decades ago, but within the last few years, it proved itself as a powerful tool to address scale-up and -down topics of bioprocesses. It can capture the history of a cell and reveal invaluable information for, not only, bioprocess control and design but also strain engineering. This way it will be possible to shed light on the actual environment that cell experiences throughout its lifespan. Lifelines of a microorganism in a bioreactor can serve as the missing link that encompasses the biological timescales and the physical timescales. For this purpose digitalization of bioreactors provides us with new insights that are not achievable in industrial reactors easily if at all, namely, substrate and product gradients; high-shear regions are among the most interesting factors that can be reproduced adequately with help of a digital twin. In this chapter basic principles of this method will be introduced, and later on some practical aspects of particle tracking technique will be illustrated. In the final section, some of the advantages and challenges associated with this method will be discussed.

Control and regulation of acetate overflow in Escherichia coli

Overflow metabolism refers to the production of seemingly wasteful by-products by cells during growth on glucose even when oxygen is abundant. Two theories have been proposed to explain acetate overflow in Escherichia coli – global control of the central metabolism and local control of the acetate pathway – but neither accounts for all observations. Here, we develop a kinetic model of E. coli metabolism that quantitatively accounts for observed behaviours and successfully predicts the response of E. coli to new perturbations. We reconcile these theories and clarify the origin, control, and regulation of the acetate flux. We also find that, in turns, acetate regulates glucose metabolism by coordinating the expression of glycolytic and TCA genes. Acetate should not be considered a wasteful end-product since it is also a co-substrate and a global regulator of glucose metabolism in E. coli. This has broad implications for our understanding of overflow metabolism.

Pathway Tools Visualization of Organism-Scale Metabolic Networks

Metabolomics, synthetic biology, and microbiome research demand information about organism-scale metabolic networks. The convergence of genome sequencing and computational inference of metabolic networks has enabled great progress toward satisfying that demand by generating metabolic reconstructions from the genomes of thousands of sequenced organisms. Visualization of whole metabolic networks is critical for aiding researchers in understanding, analyzing, and exploiting those reconstructions. We have developed bioinformatics software tools that automatically generate a full metabolic-network diagram for an organism, and that enable searching and analyses of the network. The software generates metabolic-network diagrams for unicellular organisms, for multi-cellular organisms, and for pan-genomes and organism communities. Search tools enable users to find genes, metabolites, enzymes, reactions, and pathways within a diagram. The diagrams are zoomable to enable researchers to study local neighborhoods in detail and to see the big picture. The diagrams also serve as tools for comparison of metabolic networks and for interpreting high-throughput datasets, including transcriptomics, metabolomics, and reaction fluxes computed by metabolic models. These data can be overlaid on the metabolic charts to produce animated zoomable displays of metabolic flux and metabolite abundance. The BioCyc.org website contains whole-network diagrams for more than 18,000 sequenced organisms. The ready availability of organism-specific metabolic network diagrams and associated tools for almost any sequenced organism are useful for researchers working to better understand the metabolism of their organism and to interpret high-throughput datasets in a metabolic context.

ScatLay: utilizing transcriptome-wide noise for identifying and visualizing differentially expressed genes

Differential expressed (DE) genes analysis is valuable for understanding comparative transcriptomics between cells, conditions or time evolution. However, the predominant way of identifying DE genes is to use arbitrary threshold fold or expression changes as cutoff. Here, we developed a more objective method, Scatter Overlay or ScatLay, to extract and graphically visualize DE genes across any two samples by utilizing their pair-wise scatter or transcriptome-wide noise, while factoring replicate variabilities. We tested ScatLay for 3 cell types: between time points for Escherichia coli aerobiosis and Saccharomyces cerevisiae hypoxia, and between untreated and Etomoxir treated Mus Musculus embryonic stem cell. As a result, we obtain 1194, 2061 and 2932 DE genes, respectively. Next, we compared these data with two widely used current approaches (DESeq2 and NOISeq) with typical twofold expression changes threshold, and show that ScatLay reveals significantly larger number of DE genes. Hence, our method provides a wider coverage of DE genes, and will likely pave way for finding more novel regulatory genes in future works.

Investigation of Bar-seq as a method to study population dynamics of Saccharomyces cerevisiae deletion library during bioreactor cultivation

Background

Despite the latest advancements in metabolic engineering for genome editing and characterization of host performance, the successful development of robust cell factories used for industrial bioprocesses and accurate prediction of the behavior of microbial systems, especially when shifting from laboratory-scale to industrial conditions, remains challenging. To increase the probability of success of a scale-up process, data obtained from thoroughly performed studies mirroring cellular responses to typical large-scale stimuli may be used to derive crucial information to better understand potential implications of large-scale cultivation on strain performance. This study assesses the feasibility to employ a barcoded yeast deletion library to assess genome-wide strain fitness across a simulated industrial fermentation regime and aims to understand the genetic basis of changes in strain physiology during industrial fermentation, and the corresponding roles these genes play in strain performance.

Results

We find that mutant population diversity is maintained through multiple seed trains, enabling large scale fermentation selective pressures to act upon the community. We identify specific deletion mutants that were enriched in all processes tested in this study, independent of the cultivation conditions, which include MCK1, RIM11, MRK1, and YGK3 that all encode homologues of mammalian glycogen synthase kinase 3 (GSK-3). Ecological analysis of beta diversity between all samples revealed significant population divergence over time and showed feed specific consequences of population structure. Further, we show that significant changes in the population diversity during fed-batch cultivations reflect the presence of significant stresses. Our observations indicate that, for this yeast deletion collection, the selection of the feeding scheme which affects the accumulation of the fermentative by-product ethanol impacts the diversity of the mutant pool to a higher degree as compared to the pH of the culture broth. The mutants that were lost during the time of most extreme population selection suggest that specific biological processes may be required to cope with these specific stresses.

Conclusions

Our results demonstrate the feasibility of Bar-seq to assess fermentation associated stresses in yeast populations under industrial conditions and to understand critical stages of a scale-up process where variability emerges, and selection pressure gets imposed. Overall our work highlights a promising avenue to identify genetic loci and biological stress responses required for fitness under industrial conditions.

Attractor Concepts to Evaluate the Transcriptome-wide Dynamics Guiding Anaerobic to Aerobic State Transition in Escherichia coli

For any dynamical system, like living organisms, an attractor state is a set of variables or mechanisms that converge towards a stable system behavior despite a wide variety of initial conditions. Here, using multi-dimensional statistics, we investigate the global gene expression attractor mechanisms shaping anaerobic to aerobic state transition (AAT) of Escherichia coli in a bioreactor at early times. Out of 3,389 RNA-Seq expression changes over time, we identified 100 sharply changing genes that are key for guiding 1700 genes into the AAT attractor basin. Collectively, these genes were named as attractor genes constituting of 6 dynamic clusters. Apart from the expected anaerobic (glycolysis), aerobic (TCA cycle) and fermentation (succinate pathways) processes, sulphur metabolism, ribosome assembly and amino acid transport mechanisms together with 332 uncharacterised genes are also key for AAT. Overall, our work highlights the importance of multi-dimensional statistical analyses for revealing novel processes shaping AAT.

Pathway Tools version 23.0 update: software for pathway/genome informatics and systems biology

MOTIVATION

Biological systems function through dynamic interactions among genes and their products, regulatory circuits and metabolic networks. Our development of the Pathway Tools software was motivated by the need to construct biological knowledge resources that combine these many types of data, and that enable users to find and comprehend data of interest as quickly as possible through query and visualization tools. Further, we sought to support the development of metabolic flux models from pathway databases, and to use pathway information to leverage the interpretation of high-throughput data sets.

RESULTS

In the past 4 years we have enhanced the already extensive Pathway Tools software in several respects. It can now support metabolic-model execution through the Web, it provides a more accurate gap filler for metabolic models; it supports development of models for organism communities distributed across a spatial grid; and model results may be visualized graphically. Pathway Tools supports several new omics-data analysis tools including the Omics Dashboard, multi-pathway diagrams called pathway collages, a pathway-covering algorithm for metabolomics data analysis and an algorithm for generating mechanistic explanations of multi-omics data. We have also improved the core pathway/genome databases management capabilities of the software, providing new multi-organism search tools for organism communities, improved graphics rendering, faster performance and re-designed gene and metabolite pages.

AVAILABILITY

The software is free for academic use; a fee is required for commercial use. See http://pathwaytools.com.

CONTACT

pkarp@ai.sri.com.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Briefings in Bioinformatics online.

Engineering Robust Production Microbes for Large-Scale Cultivation

Systems biology and synthetic biology are increasingly used to examine and modulate complex biological systems. As such, many issues arising during scaling-up microbial production processes can be addressed using these approaches. We review differences between laboratory-scale cultures and larger-scale processes to provide a perspective on those strain characteristics that are especially important during scaling. Systems biology has been used to examine a range of microbial systems for their response in bioreactors to fluctuations in nutrients, dissolved gases, and other stresses. Synthetic biology has been used both to assess and modulate strain response, and to engineer strains to improve production. We discuss these approaches and tools in the context of their use in engineering robust microbes for applications in large-scale production.

In Silico Prediction of Large-Scale Microbial Production Performance: Constraints for Getting Proper Data-Driven Models

Industrial bioreactors range from 10.000 to 700.000 L and characteristically show different zones of substrate availabilities, dissolved gas concentrations and pH values reflecting physical, technical and economic constraints of scale-up. Microbial producers are fluctuating inside the bioreactors thereby experiencing frequently changing micro-environmental conditions. The external stimuli induce responses on microbial metabolism and on transcriptional regulation programs. Both may deteriorate the expected microbial production performance in large scale compared to expectations deduced from ideal, well-mixed lab-scale conditions. Accordingly, predictive tools are needed to quantify large-scale impacts considering bioreactor heterogeneities. The review shows that the time is right to combine simulations of microbial kinetics with calculations of large-scale environmental conditions to predict the bioreactor performance. Accordingly, basic experimental procedures and computational tools are presented to derive proper microbial models and hydrodynamic conditions, and to link both for bioreactor modeling. Particular emphasis is laid on the identification of gene regulatory networks as the implementation of such models will surely gain momentum in future studies.

The Omics Dashboard for interactive exploration of gene-expression data

The Omics Dashboard is a software tool for interactive exploration and analysis of gene-expression datasets. The Omics Dashboard is organized as a hierarchy of cellular systems. At the highest level of the hierarchy the Dashboard contains graphical panels depicting systems such as biosynthesis, energy metabolism, regulation and central dogma. Each of those panels contains a series of X–Y plots depicting expression levels of subsystems of that panel, e.g. subsystems within the central dogma panel include transcription, translation and protein maturation and folding. The Dashboard presents a visual read-out of the expression status of cellular systems to facilitate a rapid top-down user survey of how all cellular systems are responding to a given stimulus, and to enable the user to quickly view the responses of genes within specific systems of interest. Although the Dashboard is complementary to traditional statistical methods for analysis of gene-expression data, we show how it can detect changes in gene expression that statistical techniques may overlook. We present the capabilities of the Dashboard using two case studies: the analysis of lipid production for the marine alga Thalassiosira pseudonana, and an investigation of a shift from anaerobic to aerobic growth for the bacterium Escherichia coli.