INCAWrapper: a Python wrapper for INCA for seamless data import, -export, and -processing.

Motivation: INCA is a powerful tool for metabolic flux analysis, however, import and export of data and results can be tedious and limit the use of INCA in automated workflows. Results: The INCAWrapper enables the use of INCA purely through Python, which allows the use of INCA in common data science workflows. Availability and implementation: The INCAWrapper is implemented in Python and can be found at https://github.com/biosustain/incawrapper. It is freely available under an MIT License. To run INCA, the user needs their own MATLAB and INCA licenses. INCA is freely available for noncommercial use at mfa.vueinnovations.com.

In Search of Disentanglement in Tandem Mass Spectrometry Datasets.

Generative modeling and representation learning of tandem mass spectrometry data aim to learn an interpretable and instrument-agnostic digital representation of metabolites directly from MS/MS spectra. Interpretable and instrument-agnostic digital representations would facilitate comparisons of MS/MS spectra between instrument vendors and enable better and more accurate queries of large MS/MS spectra databases for metabolite identification. In this study, we apply generative modeling and representation learning using variational autoencoders to understand the extent to which tandem mass spectra can be disentangled into their factors of generation (e.g., collision energy, ionization mode, instrument type, etc.) with minimal prior knowledge of the factors. We find that variational autoencoders can disentangle tandem mass spectra data with the proper choice of hyperparameters into meaningful latent representations aligned with known factors of variation. We develop a two-step approach to facilitate the selection of models that are disentangled, which could be applied to other complex and high-dimensional data sets.

An automated workflow for multi-omics screening of microbial model organisms.

Multi-omics datasets are becoming of key importance to drive discovery in fundamental research as much as generating knowledge for applied biotechnology. However, the construction of such large datasets is usually time-consuming and expensive. Automation might enable to overcome these issues by streamlining workflows from sample generation to data analysis. Here, we describe the construction of a complex workflow for the generation of high-throughput microbial multi-omics datasets. The workflow comprises a custom-built platform for automated cultivation and sampling of microbes, sample preparation protocols, analytical methods for sample analysis and automated scripts for raw data processing. We demonstrate possibilities and limitations of such workflow in generating data for three biotechnologically relevant model organisms, namely Escherichia coli, Saccharomyces cerevisiae, and Pseudomonas putida.

Oral supplementation of nicotinamide riboside alters intestinal microbial composition in rats and mice, but not humans.

The gut microbiota impacts systemic levels of multiple metabolites including NAD+ precursors through diverse pathways. Nicotinamide riboside (NR) is an NAD+ precursor capable of regulating mammalian cellular metabolism. Some bacterial families express the NR-specific transporter, PnuC. We hypothesized that dietary NR supplementation would modify the gut microbiota across intestinal sections. We determined the effects of 12 weeks of NR supplementation on the microbiota composition of intestinal segments of high-fat diet-fed (HFD) rats. We also explored the effects of 12 weeks of NR supplementation on the gut microbiota in humans and mice. In rats, NR reduced fat mass and tended to decrease body weight. Interestingly, NR increased fat and energy absorption but only in HFD-fed rats. Moreover, 16S rRNA gene sequencing analysis of intestinal and fecal samples revealed an increased abundance of species within Erysipelotrichaceae and Ruminococcaceae families in response to NR. PnuC-positive bacterial strains within these families showed an increased growth rate when supplemented with NR. The abundance of species within the Lachnospiraceae family decreased in response to HFD irrespective of NR. Alpha and beta diversity and bacterial composition of the human fecal microbiota were unaltered by NR, but in mice, the fecal abundance of species within Lachnospiraceae increased while abundances of Parasutterella and Bacteroides dorei species decreased in response to NR. In conclusion, oral NR altered the gut microbiota in rats and mice, but not in humans. In addition, NR attenuated body fat mass gain in rats, and increased fat and energy absorption in the HFD context.

Automating the design-build-test-learn cycle towards next-generation bacterial cell factories.

Automation is playing an increasingly significant role in synthetic biology. Groundbreaking technologies, developed over the past 20 years, have enormously accelerated the construction of efficient microbial cell factories. Integrating state-of-the-art tools (e.g. for genome engineering and analytical techniques) into the design-build-test-learn cycle (DBTLc) will shift the metabolic engineering paradigm from an almost artisanal labor towards a fully automated workflow. Here, we provide a perspective on how a fully automated DBTLc could be harnessed to construct the next-generation bacterial cell factories in a fast, high-throughput fashion. Innovative toolsets and approaches that pushed the boundaries in each segment of the cycle are reviewed to this end. We also present the most recent efforts on automation of the DBTLc, which heralds a fully autonomous pipeline for synthetic biology in the near future.

Precision diagnostic approach to predict 5-year risk for microvascular complications in type 1 diabetes.

BACKGROUND: Individuals with long standing diabetes duration can experience damage to small microvascular blood vessels leading to diabetes complications (DCs) and increased mortality. Precision diagnostic tailors a diagnosis to an individual by using biomedical information. Blood small molecule profiling coupled with machine learning (ML) can facilitate the goals of precision diagnostics, including earlier diagnosis and individualized risk scoring. METHODS: Using data in a cohort of 537 adults with type 1 diabetes (T1D) we predicted five-year progression to DCs. Prediction models were computed first with clinical risk factors at baseline and then with clinical risk factors and blood-derived molecular data at baseline. Progression of diabetic kidney disease and diabetic retinopathy were predicted in two complication-specific models. FINDINGS: The model predicts the progression to diabetic kidney disease with accuracy: 0.96 ± 0.25 and 0.96 ± 0.06 area under curve, AUC, with clinical measurements and with small molecule predictors respectively and highlighted main predictors to be albuminuria, glomerular filtration rate, retinopathy status at baseline, sugar derivatives and ketones. For diabetic retinopathy, AUC 0.75 ± 0.14 and 0.79 ± 0.16 with clinical measurements and with small molecule predictors respectively and highlighted key predictors, albuminuria, glomerular filtration rate and retinopathy status at baseline. Individual risk scores were built to visualize results. INTERPRETATION: With further validation ML tools could facilitate the implementation of precision diagnosis in the clinic. It is envisaged that patients could be screened for complications, before these occur, thus preserving healthy life-years for persons with diabetes. FUNDING: This study has been financially supported by Novo Nordisk Foundation grant NNF14OC0013659.

A Comprehensive Evaluation of Metabolomics Data Preprocessing Methods for Deep Learning.

Machine learning has greatly advanced over the past decade, owing to advances in algorithmic innovations, hardware acceleration, and benchmark datasets to train on domains such as computer vision, natural-language processing, and more recently the life sciences. In particular, the subfield of machine learning known as deep learning has found applications in genomics, proteomics, and metabolomics. However, a thorough assessment of how the data preprocessing methods required for the analysis of life science data affect the performance of deep learning is lacking. This work contributes to filling that gap by assessing the impact of commonly used as well as newly developed methods employed in data preprocessing workflows for metabolomics that span from raw data to processed data. The results from these analyses are summarized into a set of best practices that can be used by researchers as a starting point for downstream classification and reconstruction tasks using deep learning.

SmartPeak Automates Targeted and Quantitative Metabolomics Data Processing.

Technological advances in high-resolution mass spectrometry (MS) vastly increased the number of samples that can be processed in a life science experiment, as well as volume and complexity of the generated data. To address the bottleneck of high-throughput data processing, we present SmartPeak (https://github.com/AutoFlowResearch/SmartPeak), an application that encapsulates advanced algorithms to enable fast, accurate, and automated processing of capillary electrophoresis-, gas chromatography-, and liquid chromatography (LC)-MS(/MS) data and high-pressure LC data for targeted and semitargeted metabolomics, lipidomics, and fluxomics experiments. The application allows for an approximate 100-fold reduction in the data processing time compared to manual processing while enhancing quality and reproducibility of the results.

High-Resolution Scanning of Optimal Biosensor Reporter Promoters in Yeast.

Small-molecule binding allosteric transcription factors (aTFs) derived from bacteria enable real-time monitoring of metabolite abundances, high-throughput screening of genetic designs, and dynamic control of metabolism. Yet, engineering of reporter promoter designs of prokaryotic aTF biosensors in eukaryotic cells is complex. Here we investigate the impact of aTF binding site positions at single-nucleotide resolution in >300 reporter promoter designs in Saccharomyces cerevisiae. From this we identify biosensor output landscapes with transient and distinct aTF binding site position effects for aTF repressors and activators, respectively. Next, we present positions for tunable reporter promoter outputs enabling metabolite-responsive designs for a total of four repressor-type and three activator-type aTF biosensors with dynamic output ranges up to 8- and 26-fold, respectively. This study highlights aTF binding site positions in reporter promoters as key for successful biosensor engineering and that repressor-type aTF biosensors allows for more flexibility in terms of choice of binding site positioning compared to activator-type aTF biosensors.

Corrigendum: Multi-Omics Analysis of Fatty Alcohol Production in Engineered Yeasts Saccharomyces cerevisiae and Yarrowia lipolytica.

[This corrects the article DOI: 10.3389/fgene.2019.00747.].

Multi-Omics Analysis of Fatty Alcohol Production in Engineered Yeasts Saccharomyces cerevisiae and Yarrowia lipolytica.

Fatty alcohols are widely used in various applications within a diverse set of industries, such as the soap and detergent industry, the personal care, and cosmetics industry, as well as the food industry. The total world production of fatty alcohols is over 2 million tons with approximately equal parts derived from fossil oil and from plant oils or animal fats. Due to the environmental impact of these production methods, there is an interest in alternative methods for fatty alcohol production via microbial fermentation using cheap renewable feedstocks. In this study, we aimed to obtain a better understanding of how fatty alcohol biosynthesis impacts the host organism, baker's yeast Saccharomyces cerevisiae or oleaginous yeast Yarrowia lipolytica. Producing and non-producing strains were compared in growth and nitrogen-depletion cultivation phases. The multi-omics analysis included physiological characterization, transcriptome analysis by RNAseq, 13Cmetabolic flux analysis, and intracellular metabolomics. Both species accumulated fatty alcohols under nitrogen-depletion conditions but not during growth. The fatty alcohol-producing Y. lipolytica strain had a higher fatty alcohol production rate than an analogous S. cerevisiae strain. Nitrogen-depletion phase was associated with lower glucose uptake rates and a decrease in the intracellular concentration of acetyl-CoA in both yeast species, as well as increased organic acid secretion rates in Y. lipolytica. Expression of the fatty alcohol-producing enzyme fatty acyl-CoA reductase alleviated the growth defect caused by deletion of hexadecenal dehydrogenase encoding genes (HFD1 and HFD4) in Y. lipolytica. RNAseq analysis showed that fatty alcohol production triggered a cell wall stress response in S. cerevisiae. RNAseq analysis also showed that both nitrogen-depletion and fatty alcohol production have substantial effects on the expression of transporter encoding genes in Y. lipolytica. In conclusion, through this multi-omics study, we uncovered some effects of fatty alcohol production on the host metabolism. This knowledge can be used as guidance for further strain improvement towards the production of fatty alcohols.

Laboratory evolution reveals a two-dimensional rate-yield tradeoff in microbial metabolism.

Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach.

A White-Box Machine Learning Approach for Revealing Antibiotic Mechanisms of Action.

Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism.

Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.

Growth Adaptation of gnd and sdhCB Escherichia coli Deletion Strains Diverges From a Similar Initial Perturbation of the Transcriptome.

Adaptive laboratory evolution (ALE) has emerged as a new approach with which to pursue fundamental biological inquiries and, in particular, new insights into the systemic function of a gene product. Two E. coli knockout strains were constructed: one that blocked the Pentose Phosphate Pathway (gnd KO) and one that decoupled the TCA cycle from electron transport (sdhCDAB KO). Despite major perturbations in central metabolism, minimal growth rate changes were found in the two knockout strains. More surprisingly, many similarities were found in their initial transcriptomic states that could be traced to similarly perturbed metabolites despite the differences in the network location of the gene perturbations and concomitant re-routing of pathway fluxes around these perturbations. However, following ALE, distinct metabolomic and transcriptomic states were realized. These included divergent flux and gene expression profiles in the gnd and sdhCDAB KOs to overcome imbalances in NADPH production and nitrogen/sulfur assimilation, respectively, that were not obvious limitations of growth in the unevolved knockouts. Therefore, this work demonstrates that ALE provides a productive approach to reveal novel insights of gene function at a systems level that cannot be found by observing the fresh knockout alone.

Multiple Optimal Phenotypes Overcome Redox and Glycolytic Intermediate Metabolite Imbalances in Escherichia coli pgi Knockout Evolutions.

A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.

Adaptive laboratory evolution resolves energy depletion to maintain high aromatic metabolite phenotypes in Escherichia coli strains lacking the Phosphotransferase System.

Aromatic metabolites provide the backbone for numerous industrial and pharmaceutical compounds of high value. The Phosphotransferase System (PTS) is common to many bacteria, and is the primary mechanism for glucose uptake by Escherichia coli. The PTS was removed to conserve phosphoenolpyruvate (pep), which is a precursor for aromatic metabolites and consumed by the PTS, for aromatic metabolite production. Replicate adaptive laboratory evolution (ALE) of PTS and detailed omics data sets collected revealed that the PTS bridged the gap between respiration and fermentation, leading to distinct high fermentative and high respiratory rate phenotypes. It was also found that while all strains retained high levels of aromatic amino acid (AAA) biosynthetic precursors, only one replicate from the high glycolytic clade retained high levels of intracellular AAAs. The fast growth and high AAA precursor phenotypes could provide a starting host for cell factories targeting the overproduction aromatic metabolites.

Adaptation to the coupling of glycolysis to toxic methylglyoxal production in tpiA deletion strains of Escherichia coli requires synchronized and counterintuitive genetic changes.

Methylglyoxal is a highly toxic metabolite that can be produced in all living organisms. Methylglyoxal was artificially elevated by removal of the tpiA gene from a growth optimized Escherichia coli strain. The initial response to elevated methylglyoxal and its toxicity was characterized, and detoxification mechanisms were studied using adaptive laboratory evolution. We found that: 1) Multi-omics analysis revealed biological consequences of methylglyoxal toxicity, which included attack on macromolecules including DNA and RNA and perturbation of nucleotide levels; 2) Counter-intuitive cross-talk between carbon starvation and inorganic phosphate signalling was revealed in the tpiA deletion strain that required mutations in inorganic phosphate signalling mechanisms to alleviate; and 3) The split flux through lower glycolysis depleted glycolytic intermediates requiring a host of synchronized and coordinated mutations in non-intuitive network locations in order to re-adjust the metabolic flux map to achieve optimal growth. Such mutations included a systematic inactivation of the Phosphotransferase System (PTS) and alterations in cell wall biosynthesis enzyme activity. This study demonstrated that deletion of major metabolic genes followed by ALE was a productive approach to gain novel insight into the systems biology underlying optimal phenotypic states.

RapidRIP quantifies the intracellular metabolome of 7 industrial strains of E. coli.

Fast metabolite quantification methods are required for high throughput screening of microbial strains obtained by combinatorial or evolutionary engineering approaches. In this study, a rapid RIP-LC-MS/MS (RapidRIP) method for high-throughput quantitative metabolomics was developed and validated that was capable of quantifying 102 metabolites from central, amino acid, energy, nucleotide, and cofactor metabolism in less than 5 minutes. The method was shown to have comparable sensitivity and resolving capability as compared to a full length RIP-LC-MS/MS method (FullRIP). The RapidRIP method was used to quantify the metabolome of seven industrial strains of E. coli revealing significant differences in glycolytic, pentose phosphate, TCA cycle, amino acid, and energy and cofactor metabolites were found. These differences translated to statistically and biologically significant differences in thermodynamics of biochemical reactions between strains that could have implications when choosing a host for bioprocessing.