hopsy - a methods marketplace for convex polytope sampling in Python.

SUMMARY: Effective collaboration between developers of Bayesian inference methods and users is key to advance our quantitative understanding of biosystems. We here present hopsy, a versatile open-source platform designed to provide convenient access to powerful Markov chain Monte Carlo sampling algorithms tailored to models defined on convex polytopes (CP). Based on the high-performance C++ sampling library HOPS, hopsy inherits its strengths and extends its functionalities with the accessibility of the Python programming language. A versatile plugin-mechanism enables seamless integration with domain-specific models, providing method developers with a framework for testing, benchmarking, and distributing CP samplers to approach real-world inference tasks. We showcase hopsy by solving common and newly composed domain-specific sampling problems, highlighting important design choices. By likening hopsy to a marketplace, we emphasize its role in bringing together users and developers, where users get access to state-of-the-art methods, and developers contribute their own innovative solutions for challenging domain-specific inference problems. AVAILABILITY AND IMPLEMENTATION: Sources, documentation and a continuously updated list of sampling algorithms are available at https://jugit.fz-juelich.de/IBG-1/ModSim/hopsy, with Linux, Windows and MacOS binaries at https://pypi.org/project/hopsy/.

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference.

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

Growth-rate dependency of ribosome abundance and translation elongation rate in Corynebacterium glutamicum differs from that in Escherichia coli.

Bacterial growth rate (µ) depends on the protein synthesis capacity of the cell and thus on the number of active ribosomes and their translation elongation rate. The relationship between these fundamental growth parameters have only been described for few bacterial species, in particular Escherichia coli. Here, we analyse the growth-rate dependency of ribosome abundance and translation elongation rate for Corynebacterium glutamicum, a gram-positive model species differing from E. coli by a lower growth temperature optimum and a lower maximal growth rate. We show that, unlike in E. coli, there is little change in ribosome abundance for µ <0.4 h-1 in C. glutamicum and the fraction of active ribosomes is kept above 70% while the translation elongation rate declines 5-fold. Mathematical modelling indicates that the decrease in the translation elongation rate can be explained by a depletion of translation precursors.

A pathogen-specific isotope tracing approach reveals metabolic activities and fluxes of intracellular Salmonella.

Pathogenic bacteria proliferating inside mammalian host cells need to rapidly adapt to the intracellular environment. How they achieve this and scavenge essential nutrients from the host has been an open question due to the difficulties in distinguishing between bacterial and host metabolites in situ. Here, we capitalized on the inability of mammalian cells to metabolize mannitol to develop a stable isotopic labeling approach to track Salmonella enterica metabolites during intracellular proliferation in host macrophage and epithelial cells. By measuring label incorporation into Salmonella metabolites with liquid chromatography-mass spectrometry (LC-MS), and combining it with metabolic modeling, we identify relevant carbon sources used by Salmonella, uncover routes of their metabolization, and quantify relative reaction rates in central carbon metabolism. Our results underline the importance of the Entner-Doudoroff pathway (EDP) and the phosphoenolpyruvate carboxylase for intracellularly proliferating Salmonella. More broadly, our metabolic labeling strategy opens novel avenues for understanding the metabolism of pathogens inside host cells.

Practical sampling of constraint-based models: Optimized thinning boosts CHRR performance.

Thinning is a sub-sampling technique to reduce the memory footprint of Markov chain Monte Carlo. Despite being commonly used, thinning is rarely considered efficient. For sampling constraint-based models, a highly relevant use-case in systems biology, we here demonstrate that thinning boosts computational and, thereby, sampling efficiencies of the widely used Coordinate Hit-and-Run with Rounding (CHRR) algorithm. By benchmarking CHRR with thinning with simplices and genome-scale metabolic networks of up to thousands of dimensions, we find a substantial increase in computational efficiency compared to unthinned CHRR, in our examples by orders of magnitude, as measured by the effective sample size per time (ESS/t), with performance gains growing with polytope (effective network) dimension. Using a set of benchmark models we derive a ready-to-apply guideline for tuning thinning to efficient and effective use of compute resources without requiring additional coding effort. Our guideline is validated using three (out-of-sample) large-scale networks and we show that it allows sampling convex polytopes uniformly to convergence in a fraction of time, thereby unlocking the rigorous investigation of hitherto intractable models. The derivation of our guideline is explained in detail, allowing future researchers to update it as needed as new model classes and more training data becomes available. CHRR with deliberate utilization of thinning thereby paves the way to keep pace with progressing model sizes derived with the constraint-based reconstruction and analysis (COBRA) tool set. Sampling and evaluation pipelines are available at https://jugit.fz-juelich.de/IBG-1/ModSim/fluxomics/chrrt.

Enabling oxygen-controlled microfluidic cultures for spatiotemporal microbial single-cell analysis.

Microfluidic cultivation devices that facilitate O2 control enable unique studies of the complex interplay between environmental O2 availability and microbial physiology at the single-cell level. Therefore, microbial single-cell analysis based on time-lapse microscopy is typically used to resolve microbial behavior at the single-cell level with spatiotemporal resolution. Time-lapse imaging then provides large image-data stacks that can be efficiently analyzed by deep learning analysis techniques, providing new insights into microbiology. This knowledge gain justifies the additional and often laborious microfluidic experiments. Obviously, the integration of on-chip O2 measurement and control during the already complex microfluidic cultivation, and the development of image analysis tools, can be a challenging endeavor. A comprehensive experimental approach to allow spatiotemporal single-cell analysis of living microorganisms under controlled O2 availability is presented here. To this end, a gas-permeable polydimethylsiloxane microfluidic cultivation chip and a low-cost 3D-printed mini-incubator were successfully used to control O2 availability inside microfluidic growth chambers during time-lapse microscopy. Dissolved O2 was monitored by imaging the fluorescence lifetime of the O2-sensitive dye RTDP using FLIM microscopy. The acquired image-data stacks from biological experiments containing phase contrast and fluorescence intensity data were analyzed using in-house developed and open-source image-analysis tools. The resulting oxygen concentration could be dynamically controlled between 0% and 100%. The system was experimentally tested by culturing and analyzing an E. coli strain expressing green fluorescent protein as an indirect intracellular oxygen indicator. The presented system allows for innovative microbiological research on microorganisms and microbial ecology with single-cell resolution.

One-shot 13 C15 N-metabolic flux analysis for simultaneous quantification of carbon and nitrogen flux.

Metabolic flux is the final output of cellular regulation and has been extensively studied for carbon but much less is known about nitrogen, which is another important building block for living organisms. For the tuberculosis pathogen, this is particularly important in informing the development of effective drugs targeting the pathogen's metabolism. Here we performed 13 C15 N dual isotopic labeling of Mycobacterium bovis BCG steady state cultures, quantified intracellular carbon and nitrogen fluxes and inferred reaction bidirectionalities. This was achieved by model scope extension and refinement, implemented in a multi-atom transition model, within the statistical framework of Bayesian model averaging (BMA). Using BMA-based 13 C15 N-metabolic flux analysis, we jointly resolve carbon and nitrogen fluxes quantitatively. We provide the first nitrogen flux distributions for amino acid and nucleotide biosynthesis in mycobacteria and establish glutamate as the central node for nitrogen metabolism. We improved resolution of the notoriously elusive anaplerotic node in central carbon metabolism and revealed possible operation modes. Our study provides a powerful and statistically rigorous platform to simultaneously infer carbon and nitrogen metabolism in any biological system.

microbeSEG: A deep learning software tool with OMERO data management for efficient and accurate cell segmentation.

In biotechnology, cell growth is one of the most important properties for the characterization and optimization of microbial cultures. Novel live-cell imaging methods are leading to an ever better understanding of cell cultures and their development. The key to analyzing acquired data is accurate and automated cell segmentation at the single-cell level. Therefore, we present microbeSEG, a user-friendly Python-based cell segmentation tool with a graphical user interface and OMERO data management. microbeSEG utilizes a state-of-the-art deep learning-based segmentation method and can be used for instance segmentation of a wide range of cell morphologies and imaging techniques, e.g., phase contrast or fluorescence microscopy. The main focus of microbeSEG is a comprehensible, easy, efficient, and complete workflow from the creation of training data to the final application of the trained segmentation model. We demonstrate that accurate cell segmentation results can be obtained within 45 minutes of user time. Utilizing public segmentation datasets or pre-labeling further accelerates the microbeSEG workflow. This opens the door for accurate and efficient data analysis of microbial cultures.

Automated Characterization of Catalytically Active Inclusion Body Production in Biotechnological Screening Systems.

We here propose an automated pipeline for the microscopy image-based characterization of catalytically active inclusion bodies (CatIBs), which includes a fully automatic experimental high-throughput workflow combined with a hybrid approach for multi-object microbial cell segmentation. For automated microscopy, a CatIB producer strain was cultivated in a microbioreactor from which samples were injected into a flow chamber. The flow chamber was fixed under a microscope and an integrated camera took a series of images per sample. To explore heterogeneity of CatIB development during the cultivation and track the size and quantity of CatIBs over time, a hybrid image processing pipeline approach was developed, which combines an ML-based detection of in-focus cells with model-based segmentation. The experimental setup in combination with an automated image analysis unlocks high-throughput screening of CatIB production, saving time and resources. Biotechnological relevance- CatIBs have wide application in synthetic chemistry and biocatalysis, but also could have future biomedical applications such as therapeutics. The proposed hybrid automatic image processing pipeline can be adjusted to treat comparable biological microorganisms, where fully data-driven ML-based segmentation approaches are not feasible due to the lack of training data. Our work is the first step towards image- based bioprocess control.

Communities of Niche-optimized Strains (CoNoS) - Design and creation of stable, genome-reduced co-cultures.

Current bioprocesses for production of value-added compounds are mainly based on pure cultures that are composed of rationally engineered strains of model organisms with versatile metabolic capacities. However, in the comparably well-defined environment of a bioreactor, metabolic flexibility provided by various highly abundant biosynthetic enzymes is much less required and results in suboptimal use of carbon and energy sources for compound production. In nature, non-model organisms have frequently evolved in communities where genome-reduced, auxotrophic strains cross-feed each other, suggesting that there must be a significant advantage compared to growth without cooperation. To prove this, we started to create and study synthetic communities of niche-optimized strains (CoNoS) that consists of two strains of the same species Corynebacterium glutamicum that are mutually dependent on one amino acid. We used both the wild-type and the genome-reduced C1* chassis for introducing selected amino acid auxotrophies, each based on complete deletion of all required biosynthetic genes. The best candidate strains were used to establish several stably growing CoNoS that were further characterized and optimized by metabolic modelling, microfluidic experiments and rational metabolic engineering to improve amino acid production and exchange. Finally, the engineered CoNoS consisting of an l-leucine and l-arginine auxotroph showed a specific growth rate equivalent to 83% of the wild type in monoculture, making it the fastest co-culture of two auxotrophic C. glutamicum strains to date. Overall, our results are a first promising step towards establishing improved biobased production of value-added compounds using the CoNoS approach.

PolyRound: polytope rounding for random sampling in metabolic networks.

SUMMARY: Random flux sampling is a powerful tool for the constraint-based analysis of metabolic networks. The most efficient sampling method relies on a rounding transform of the constraint polytope, but no available rounding implementation can round all relevant models. By removing redundant polytope constraints on the go, PolyRound simplifies the numerical problem and rounds all the 108 models in the BiGG database without parameter tuning, compared to ∼50% for the state-of-the-art implementation. AVAILABILITY AND IMPLEMENTATION: The implementation is available on gitlab: https://gitlab.com/csb.ethz/PolyRound. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CellSium: versatile cell simulator for microcolony ground truth generation.

Summary: To train deep learning-based segmentation models, large ground truth datasets are needed. To address this need in microfluidic live-cell imaging, we present CellSium, a flexibly configurable cell simulator built to synthesize realistic image sequences of bacterial microcolonies growing in monolayers. We illustrate that the simulated images are suitable for training neural networks. Synthetic time-lapse videos with and without fluorescence, using programmable cell growth models, and simulation-ready 3D colony geometries for computational fluid dynamics are also supported. Availability and implementation: CellSium is free and open source software under the BSD license, implemented in Python, available at github.com/modsim/cellsium (DOI: 10.5281/zenodo.6193033), along with documentation, usage examples and Docker images. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

High-Quality Genome-Scale Reconstruction of Corynebacterium glutamicum ATCC 13032.

Corynebacterium glutamicum belongs to the microbes of enormous biotechnological relevance. In particular, its strain ATCC 13032 is a widely used producer of L-amino acids at an industrial scale. Its apparent robustness also turns it into a favorable platform host for a wide range of further compounds, mainly because of emerging bio-based economies. A deep understanding of the biochemical processes in C. glutamicum is essential for a sustainable enhancement of the microbe's productivity. Computational systems biology has the potential to provide a valuable basis for driving metabolic engineering and biotechnological advances, such as increased yields of healthy producer strains based on genome-scale metabolic models (GEMs). Advanced reconstruction pipelines are now available that facilitate the reconstruction of GEMs and support their manual curation. This article presents iCGB21FR, an updated and unified GEM of C. glutamicum ATCC 13032 with high quality regarding comprehensiveness and data standards, built with the latest modeling techniques and advanced reconstruction pipelines. It comprises 1042 metabolites, 1539 reactions, and 805 genes with detailed annotations and database cross-references. The model validation took place using different media and resulted in realistic growth rate predictions under aerobic and anaerobic conditions. The new GEM produces all canonical amino acids, and its phenotypic predictions are consistent with laboratory data. The in silico model proved fruitful in adding knowledge to the metabolism of C. glutamicum: iCGB21FR still produces L-glutamate with the knock-out of the enzyme pyruvate carboxylase, despite the common belief to be relevant for the amino acid's production. We conclude that integrating high standards into the reconstruction of GEMs facilitates replicating validated knowledge, closing knowledge gaps, and making it a useful basis for metabolic engineering. The model is freely available from BioModels Database under identifier MODEL2102050001.

Robustifying Experimental Tracer Design for13C-Metabolic Flux Analysis.

13C metabolic flux analysis (MFA) has become an indispensable tool to measure metabolic reaction rates (fluxes) in living organisms, having an increasingly diverse range of applications. Here, the choice of the13C labeled tracer composition makes the difference between an information-rich experiment and an experiment with only limited insights. To improve the chances for an informative labeling experiment, optimal experimental design approaches have been devised for13C-MFA, all relying on some a priori knowledge about the actual fluxes. If such prior knowledge is unavailable, e.g., for research organisms and producer strains, existing methods are left with a chicken-and-egg problem. In this work, we present a general computational method, termed robustified experimental design (R-ED), to guide the decision making about suitable tracer choices when prior knowledge about the fluxes is lacking. Instead of focusing on one mixture, optimal for specific flux values, we pursue a sampling based approach and introduce a new design criterion, which characterizes the extent to which mixtures are informative in view of all possible flux values. The R-ED workflow enables the exploration of suitable tracer mixtures and provides full flexibility to trade off information and cost metrics. The potential of the R-ED workflow is showcased by applying the approach to the industrially relevant antibiotic producer Streptomyces clavuligerus, where we suggest informative, yet economic labeling strategies.

HOPS: high-performance library for (non-)uniform sampling of convex-constrained models.

SUMMARY: The C++ library Highly Optimized Polytope Sampling (HOPS) provides implementations of efficient and scalable algorithms for sampling convex-constrained models that are equipped with arbitrary target functions. For uniform sampling, substantial performance gains were achieved compared to the state-of-the-art. The ease of integration and utility of non-uniform sampling is showcased in a Bayesian inference setting, demonstrating how HOPS interoperates with third-party software. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/modsim/hops/, tested on Linux and MS Windows, includes unit tests, detailed documentation, example applications and a Dockerfile. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Revisiting the Growth Modulon of Corynebacterium glutamicum Under Glucose Limited Chemostat Conditions.

Increasing the growth rate of the industrial host Corynebacterium glutamicum is a promising target to rise productivities of growth coupled product formation. As a prerequisite, detailed knowledge about the tight regulation network is necessary for identifying promising metabolic engineering goals. Here, we present comprehensive metabolic and transcriptional analysis of C. glutamicum ATCC 13032 growing under glucose limited chemostat conditions with µ = 0.2, 0.3, and 0.4 h-1. Intermediates of central metabolism mostly showed rising pool sizes with increasing growth. 13C-metabolic flux analysis (13C-MFA) underlined the fundamental role of central metabolism for the supply of precursors, redox, and energy equivalents. Global, growth-associated, concerted transcriptional patterns were not detected giving rise to the conclusion that glycolysis, pentose-phosphate pathway, and citric acid cycle are predominately metabolically controlled under glucose-limiting chemostat conditions. However, evidence is found that transcriptional regulation takes control over glycolysis once glucose-rich growth conditions are installed.

Reversible jump MCMC for multi-model inference in Metabolic Flux Analysis.

MOTIVATION: The validity of model based inference, as used in systems biology, depends on the underlying model formulation. Often, a vast number of competing models is available, that are built on different assumptions, all consistent with the existing knowledge about the studied biological phenomenon. As a remedy for this, Bayesian Model Averaging (BMA) facilitates parameter and structural inferences based on multiple models simultaneously. However, in fields where a vast number of alternative, high-dimensional and non-linear models are involved, the BMA-based inference task is computationally very challenging. RESULTS: Here we use BMA in the complex setting of Metabolic Flux Analysis (MFA) to infer whether potentially reversible reactions proceed uni- or bidirectionally, using 13C labeling data and metabolic networks. BMA is applied on a large set of candidate models with differing directionality settings, using a tailored multi-model Markov Chain Monte Carlo (MCMC) approach. The applicability of our algorithm is shown by inferring the in vivo probability of reaction bidirectionalities in a realistic network setup, thereby extending the scope of 13C MFA from parameter to structural inference. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Separation of 13C- and 15N-Isotopologues of Amino Acids with a Primary Amine without Mass Resolution by Means of O-Phthalaldehyde Derivatization and Collision Induced Dissociation.

Computational and experimental advances of recent years have culminated in establishing 13C-Metabolic Flux Analysis (13C-MFA) as a routine methodology to unravel the fluxome. As the acronym suggests, 13C-MFA has relied on the relative abundance of 13C-isotopes in metabolites for flux inference, most commonly measured by mass spectrometry. In this manuscript we expand the scope of labeling measurements to the case of simultaneous 13C- and 15N-labeling of amino acids. Analytically, the separation of isotopologues of this metabolite class can only be achieved at resolving power beyond 65,000. In this manuscript we harvest an overlooked property of the collision induced dissociation of amino acid adducts to discern 13C- and 15N- isotopologues of amino acids with a primary amine without separating them in the m/z domain.

mycelyso - high-throughput analysis of Streptomyces mycelium live cell imaging data.

BACKGROUND: Streptomycetes are filamentous microorganisms of high biotechnological relevance, especially for the production of antibiotics. In submerged cultures, the productivity of these microorganisms is closely linked to their growth morphology. Microfluidic lab-on-a-chip cultivation systems, coupled with automated time-lapse imaging, generate spatio-temporal insights into the mycelium development of streptomycetes, therewith extending the biotechnological toolset by spatio-temporal screening under well-controlled and reproducible conditions. However, the analysis of the complex mycelial structure formation is limited by the extent of manual interventions required during processing of the acquired high-volume image data. These interventions typically lead to high evaluation times and, therewith, limit the analytic throughput and exploitation of microfluidic-based screenings. RESULTS: We present the tool mycelyso (MYCElium anaLYsis SOftware), an image analysis system tailored to fully automated hyphae-level processing of image stacks generated by time-lapse microscopy. With mycelyso, the developing hyphal streptomycete network is automatically segmented and tracked over the cultivation period. Versatile key growth parameters such as mycelium network structure, its development over time, and tip growth rates are extracted. Results are presented in the web-based exploration tool mycelyso Inspector, allowing for user friendly quality control and downstream evaluation of the extracted information. In addition, 2D and 3D visualizations show temporal tracking for detailed inspection of morphological growth behaviors. For ease of getting started with mycelyso, bundled Windows packages as well as Docker images along with tutorial videos are available. CONCLUSION: mycelyso is a well-documented, platform-independent open source toolkit for the automated end-to-end analysis of Streptomyces image stacks. The batch-analysis mode facilitates the rapid and reproducible processing of large microfluidic screenings, and easy extraction of morphological parameters. The objective evaluation of image stacks is possible by reproducible evaluation workflows, useful to unravel correlations between morphological, molecular and process parameters at the hyphae- and mycelium-levels with statistical power.

Corrigendum: The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis.

[This corrects the article DOI: 10.3389/fmicb.2019.01022.].