# Exploring the Glucose Fluxotype of the E. coli y-ome Using High-Resolution Fluxomics

Cécilia Bergès, Edern Cahoreau, Pierre Millard, Brice Enjalbert, Mickael Dinclaux, Maud Heuillet, Hanna Kulyk, Lara Gales, Noémie Butin, Maxime Chazalviel, et al.

► To cite this version:

Cécilia Bergès, Edern Cahoreau, Pierre Millard, Brice Enjalbert, Mickael Dinclaux, et al. Exploring the Glucose Fluxotype of the E. coli y-ome Using High-Resolution Fluxomics. Metabolites, 2021, 11 (5), pp.271-292. 10.3390/metabo11050271. hal-03216838

HAL Id: hal-03216838

https://hal.inrae.fr/hal-03216838

Submitted on 4 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Exploring the Glucose Fluxotype of the *E. coli* y-ome Using High-Resolution Fluxomics

Cécilia Bergès<fnref n="1" />,<fnref n="2" />,† ID, Edern Cahoreau <fnref n="1" />,<fnref n="2" />,† ID, Pierre Millard <fnref n="1" /> ID, Brice Enjalbert <fnref n="1" /> ID, Mickael Dinclaux <fnref n="1" /> ID, Maud Heuillet <fnref n="1" />,<fnref n="2" />, Hanna Kulyk <fnref n="1" />,<fnref n="2" />, Lara Gales <fnref n="1" />,<fnref n="2" />, Noémie Butin <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Maxime Chazalviel <fnref n="4" /> ID, Tony Palama <fnref n="1" />,<fnref n="2" /> ID, Matthieu Guionnet <fnref n="1" />,<fnref n="2" /> ID, Sergueï Sokol <fnref n="1" /> ID, Lindsay Peyriga <fnref n="1" />,<fnref n="2" /> ID, Floriant Bellvert <fnref n="1" />,<fnref n="2" />, Stéphanie Heux <fnref n="1" /> and Jean-Charles Portais <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Bergès, C.; Cahoreau, E.; Millard, P.; Enjalbert, B.; Dinclaux, M.; Heuillet, M.; Kulyk, H.; Gales, L.; Butin, N.; Chazalviel, M.; et al. Exploring the Glucose Fluxotype of the *E. coli* y-ome Using High-Resolution Fluxomics. *Metabolites* **2021**, *11*, 271. https://doi.org/10.3390/metabo11050271

Academic Editor: Silas G. Villas-Boas

Received: 24 March 2021

Accepted: 23 April 2021

Published: 26 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** We have developed a robust workflow to measure high-resolution fluxotypes (metabolic flux phenotypes) for large strain libraries under fully controlled growth conditions. This was achieved by optimizing and automating the whole high-throughput fluxomics process and integrating all relevant software tools. This workflow allowed us to obtain highly detailed maps of carbon fluxes in the central carbon metabolism in a fully automated manner. It was applied to investigate the glucose fluxotypes of 180 *Escherichia coli* strains deleted for y-genes. Since the products of these y-genes potentially play a role in a variety of metabolic processes, the experiments were designed to be agnostic as to their potential metabolic impact. The obtained data highlight the robustness of *E. coli*'s central metabolism to y-gene deletion. For two y-genes, deletion resulted in significant changes in carbon and energy fluxes, demonstrating the involvement of the corresponding y-gene products in metabolic function or regulation. This work also introduces novel metrics to measure the actual scope and quality of high-throughput fluxomics investigations.

**Keywords:** high-resolution fluxotyping; fluxomics; high throughput; y-ome phenotyping; *E. coli*

## 1. Introduction

Despite the progress that has been made in sequencing technology and genome annotation, a non-negligible percentage of genes remains uncharacterized [1]. Even for very well-known organisms, such as the bacterium *E. coli*, 35–40% of genes (the ‘y-ome’) are of unknown function [2,3]. The products of these genes (y-genes) are likely to have highly diverse functions among the various cellular processes. Metabolism is the basic process that sustains all the energetic and biosynthetic needs of living organisms to survive and grow. Not surprisingly, attempts have been made to investigate the potential role of the y-genes in metabolic processes [4,5]. Accordingly, several high throughput (HT) methods have recently been developed to explore the potential metabolic function of y-genes. This includes the incubation of purified proteins with metabolite cocktails, the identification by mass spectrometry of metabolites with varying abundances [6] and the application of untargeted metabolomics to map gene-metabolite interactions [5].

<footnote n="1">¹ Toulouse Biotechnology Institute (TBI), Université de Toulouse, CNRS, INRAE, INSA, 31077 Toulouse, France; berges@insa-toulouse.fr (C.B.); cahoreau@insa-toulouse.fr (E.C.); millard@insa-toulouse.fr (P.M.); brice.enjalbert@insa-toulouse.fr (B.E.); mickael.dinclaux@inrae.fr (M.D.); maud.heuillet@gmail.com (M.H.); hbarbier@insa-toulouse.fr (H.K.); lara.gales@insa-toulouse.fr (L.G.); butin@insa-toulouse.fr (N.B.); tony.palama@univ-paris13.fr (T.P.); matthieu.guionnet@univ-tlse2.fr (M.G.); sokol@insa-toulouse.fr (S.S.); peyriga@insa-toulouse.fr (L.P.); bellvert@insa-toulouse.fr (F.B.); heux@insa-toulouse.fr (S.H.)</footnote>

<footnote n="2">² MetaToul-MetaboHUB, National Infrastructure of Metabolomics & Fluxomics (ANR-11-INBS-0010), 31077 Toulouse, France</footnote>

<footnote n="3">³ RESTORE, Université de Toulouse, Inserm U1031, CNRS 5070, UPS, EFS, 31100 Toulouse, France</footnote>

<footnote n="4">⁴ Toxalim (Research Centre in Food Toxicology), UMR1331, Université de Toulouse, INRAE, ENVT, INP-Purpan, UPS, 31300 Toulouse, France; maxime.chazalviel@gmail.com</footnote>

<footnote>* Correspondence: jean-charles.portais@insa-toulouse.fr</footnote>

<footnote>† Co-first authors. These two authors contributed equally to this work.</footnote>
Alongside studies aimed at directly identifying gene product function, an alternative strategy for studying the role of y-genes consists of detailing the phenotypic consequences of a loss-of-function mutation or a gene deletion [7,8]. Together with the generation of mutant libraries, such as the *Keio E. coli* mutant collection [9], which facilitates genome-wide investigations, this approach has been used for large-scale investigations of gene essentiality [9–12] and of molecular [13], morphological [14] and fitness [15] phenotypes in bacteria. Metabolic phenotyping approaches have also been developed in which comparative metabolomics is applied to reveal gene functions in yeast [4].

Metabolic phenotypes are most accurately revealed by fluxomics, which aims to measure the actual rates of biochemical reactions in metabolic networks [16]. Fluxomics measures the actual output of the integrated response of the gene-protein-metabolite interaction network [17] and provides direct access to the cellular phenotype in a quantitative manner [18]. It has, thereby, become a major tool in comprehensive investigations of cellular metabolism in many fields, ranging from biotechnology [19] to the medical sciences [20]. The measurement of metabolic fluxes is based on ¹³C-labelling experiments coupled with detailed mathematical models of metabolism (¹³C-fluxomics). This is a complex and tedious process involving many steps [16,21] and requiring significant expertise in data collection and interpretation. Obtaining the detailed flux information on hundreds of strains needed to explore the *E. coli* y-ome is a motivating factor in improving the current HT fluxomics workflows [22].

In this work, we developed a robust workflow, in order to obtain high-resolution fluxotypes under fully controlled growth conditions for large strain libraries. A fluxotype is defined here as the particular distribution of metabolic fluxes measured for a given strain under given physiological conditions. Resolution refers to the level of flux information that can be generated and is high when a significant number of fluxes can be measured [21,23,24]. This workflow was applied to investigate the fluxotypes of 180 *E. coli* strains that are deleted for y-genes and grown on glucose as sole carbon source. The data show that the central metabolism of *E. coli* is highly robust to y-gene deletion. For two y-genes, deletion resulted in significant alterations of metabolic fluxes, pointing to the role of the corresponding y-gene products in metabolic function and/or regulation.

## 2. Results

### 2.1. *Selection of E. coli y-ome Strains*

Based on the aim of investigating the metabolic phenotypes of the *E. coli* y-ome during growth on glucose as sole carbon source, we first selected y-genes as follows (Figure 1):

(i) We first considered all genes with a single-deletion mutant in the Keio mutant library [9]. This library contained single-gene deletion mutants able to grow on glucose, meaning the mutated genes are not essential for growth on glucose as sole carbon source. The Keio collection contains 3985 mutants.

(ii) We, then, selected the genes in the Keio collection lacking evidence of function. This represented a total of 1563 y-genes.

(iii) We verified that each y-gene was duly expressed and translated during growth on glucose as the sole carbon source. This selection step was based on the extensive proteomic investigation performed by Schmidt et al., in 2016 [25]. These authors measured the functional expression of 55% *E. coli* genes (>2300 genes) by quantitative proteomics in 22 experimental conditions, including growth in minimal medium with glucose as carbon source. Among the y-genes identified in step 2, we further selected the 218 y-genes that were experimentally shown to be translated under these conditions.

(iv) The y-gene status-i.e., the lack of annotation or of experimental evidence of function- was manually verified in two complementary databases, namely Biocyc (https://biocyc.org (accessed on 12 April 2019)) and Uniprot (https://www.uniprot.org (accessed on 12 April 2019)).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Strategy for high-resolution fluxotyping of the *E. coli* y-ome. Selection of *E. coli* y-genes (1) 180 y-genes were selected based on their expression on glucose as sole carbon source. High-resolution Fluxotyping (2) of 180 mutants with single-deletion of a selected y-gene. Cell cultivation, metabolite harvesting and sample preparation are performed by robotic systems. Quantification of medium compounds and isotopic profiling of metabolites by NMR and MS. Fluxes are calculated using a model of *E. coli* central carbon metabolism. Fluxotype evaluation (3) Metabolic fluxes are calculated for all strains (high-resolution fluxotypes) and compared by statistical analysis to identify altered fluxotypes. Fluxotype outliers are checked against culture data to ascertain that outlying flux values are not due to a ¹³C-labelling problem. Metabolites: Gluc.: Glucose, Glc6P: glucose-6-phosphate, Fru6P: fructose-6-phosphate, FruBP: fructose-bisphosphate, DHAP: dihydroxyacetone-phosphate, GA3P: glyceraldehyde-3-phosphate, PGA: 3-phosphoglyceric acid, PEP: phosphoenolpyruvate, Pyr: pyruvate, AcCoA: acetylCoA, OAA: oxaloacetic acid, ICit: Isocitrate, AKG: alpha-ketoglutaric acid, Suc: succinate, Mal: malate, GlyOx: glyoxylate, Gnt6P: gluconate-6-phosphate, Rib5P: ribose-5-phosphate, Ery4P: erythrose-4-phosphate, Sed7P: sedoheptulose-7-phosphate.</figcaption></figure>

This process yielded a group of 180 y-genes with unknown or unclear function, dispensable, but duly expressed in a medium with glucose as sole carbon source (Supplementary data 1). The corresponding single deletion mutants were obtained from the Keio collection and considered for further metabolic investigations, which were designed to optimize the measurement of growth parameters and fluxotypes across hundreds of *E. coli* mutants (Figure 1). Fluxotypes can be visualized using a so-called flux map, a graphical representation of the metabolic network showing the flux values for the various reactions or pathways.

## 2.2. Integrated Workflow for High-Throughput Collection of High-Resolution Fluxotypes

To measure the fluxotypes of the 180 y-gene deletion mutants, we first built an HT platform enabling fully automated, parallelized measurements of metabolic fluxes under
fully controlled conditions at a throughput consistent with the investigation of hundreds
of experimental conditions. A typical ¹³C-fluxomics workflow involves a combination
of several complex experimental and computational steps, which were improved and
optimized here to meet the needs of the y-ome investigation (Figure 1). The final setup
included (i) two automated robotic platforms, one for parallelized growth, ¹³C-labelling
experiments and sample collection, and the other for sample preparation, (ii) optimized
NMR and MS analytical methods for measuring metabolite concentrations and labeling
patterns, (iii) flux calculation, (iv) statistical analysis, and (v) a series of software programs
to store, manage, and process the data and meta-data generated all through the process.

## 2.3. *Design of Fluxomics Experiments*

Very little is known about the phenotypic and metabolic behavior of y-gene-deleted
strains, hence, an important objective of the setup was to avoid presuppositions about
the metabolic pathways potentially impacted by the lack of the y-gene product. This was
achieved by using a generic model of *E. coli* metabolism for the flux investigations, contain-
ing all central carbon metabolism pathways (Supplementary data 3). The model contained
94 fluxes in total, representing central pathways, biosynthetic processes, and transport re-
actions, and 49 metabolites. The isotopic composition of the label input in the ¹³C-labelling
experiments was optimized using the software IsoDesign [26]. Based on, (i) the network
topology, (ii) the isotopic data collected in the study (i.e., isotopologue distributions of
proteinogenic amino-acids), and (iii) the objective of resolving the maximum number of
fluxes across the whole metabolic network, the best label composition was determined
to be a mixture of 80% [1-¹³C]-glucose + 20% [U-¹³C]-glucose. This result is consistent
with those of previous investigations of the entire central carbon metabolism of *E. coli* [16].
All labeling experiments, reported in this study, were performed in minimal M9 medium
with the above mixture as sole carbon source. Details of the conditions used to apply
the workflow to all the investigations, reported here, are given in material and methods
section. Briefly, cells were inoculated at 10⁸ cell/mL and their growth was monitored by
optical density under temperature, pH and pO₂ control. Medium samples were collected
throughout the growth process and analyzed by ¹H-NMR to measure the rates of substrate
uptake and product release. When an Optical Density (OD₆₀₀nm) of 1.2-corresponding to
mid-exponential growth under these conditions-was reached, the biomass was automati-
cally sampled, hydrolyzed, and the labeling patterns of proteinogenic amino acids were
analyzed by LC-HRMS. After correction for naturally occurring isotopes [27,28], the carbon
isotopologue distributions (CIDs) of 16 amino acids were measured and used for flux
calculations. Biosynthetic fluxes (i.e., precursor requirements for growth) were calculated
based on the molecular composition of *E. coli* [29]. Intracellular fluxes were calculated by
fitting the metabolic model described above with amino acid CIDs, biosynthetic fluxes, and
extracellular fluxes, using the software influx si [30]. The confidence intervals on calculated
fluxes were estimated by Monte Carlo sensitivity analysis using the same software. All
investigations were performed on the robotic system, which allowed 48 experiments to be
run in parallel, yielding 48 flux maps per run, each containing 94 fluxes, measured under
fully controlled physiological conditions.

## 2.4. *High-Resolution Fluxotyping Workflow Validation*

Prior to starting the investigation of the *E. coli* y-ome, we evaluated the workflow
by performing several fluxomics experiments with known *E. coli* strains. Two wild-type
(WT) strains (*K-12 MG1655* and *BW25113*) and one mutant strain (*BW25113 Δzwf*), with
known and significant metabolic flux alterations, were considered to evaluate the biological
relevance of the data collected with the integrated HT workflow. The data collected for the
three strains were highly reproducible between replicates. The median relative standard
deviation (RSD) of the central metabolic fluxes was 11%, 3.7% and 22% between 5 (*K-12
MG1655*), 4 (*BW25113*), and 5 (*BW25113 Δzwf*) biological replicates, respectively. The
macro-kinetic parameters (i.e., growth rate, glucose uptake, acetate production) and flux
values collected for the two WT strains (Figure 2) were closely consistent with previous measurements on these strains, highlighting the consistency of the data collected with the integrated HT fluxomics workflow [22,31–33]. The deletion of *zwf* had little impact on the growth rate but led to a significant redistribution of metabolic fluxes to compensate for the impairment of the oxidative branch of the pentose phosphate pathway (PPP). The observed changes in metabolic fluxes in this strain were qualitatively and quantitatively consistent with previous reports [22,31,34,35]. The developed workflow’s ability to reliably measure the growth and fluxotypes of these strains clearly demonstrates the consistency of the entire process.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Metabolic fluxes measured for *E. coli* WT strains and a Δzwf mutant strain. Top: experimentally measured fluxes (expressed in mmol.gDW⁻¹.h⁻¹) and growth rate (expressed in h⁻¹). Physiological measurements: µ, growth rate; qGlc, glucose specific consumption rate; qAce, acetate specific production rate; qLac, lactate specific production rate. Bottom: Intracellular fluxes calculated from the ¹³C-labelling data; fluxes are expressed relative to the glucose uptake rate set arbitrarily to 1 for each strain. Boxplot produced with the software R representing median fluxes (bold line) in boxes representing the interquartile range (quartile 1–quartile 3). The whiskers extend to the most extreme data point that is no more than 1.5 times the interquartile range from the box. Net fluxes: Glucupt_1 and Glucupt_U: [1-¹³C]-glucose and [1-¹³C]-glucose importation (PTS system), pgi: glucose-6-phosphate isomerase, pfk: phosphofructokinase, ald: fructose-bisphosphate aldolase, tpi: triose-phosphate isomerase, pgk: phosphoglycerate kinase, eno: enolase, pyk: pyruvate kinase, pdh:</figcaption></figure>
pyruvate dehydrogenase, zwf: glucose 6-phosphate dehydrogenase, gnd: phosphogluconate dehydrogenase, tk1: half-reaction transketolase (1), ta1: half-reaction transaldolase (1), ta2: half-reaction transaldolase (2), tk2: half-reaction transketolase (2), tk3: tk2: half-reaction transketolase (3), edd: Entner-Doudoroff enzymes: equivalent to 6-phosphogluconate dehydratase and 2-keto-3deoxy-6phosphogluconate (KDPG) aldolase, citsynth: citrate synthase, idh: isocitrate dehydrogenase, akgdh: alpha-ketoglutarate dehydrogenase, fum_a and fum_b: fumarase, maldh: malate dehydrogenase, gs1: isocitrate lyase, gs2: malate synthase, mae: malic enzyme, ppc: equivalent to pep carboxylase (forward flux) and to pep carboxykinase (reverse flux), out_Ac: acetate output (equivalent to acetate kinase).

## 2.5. *High-Resolution Glucose Fluxotyping of 180 Selected Y-Gene Mutant Strains*

The 180 strains deleted for the selected y-genes were isolated from the Keio strain collection, and their growth profiles and fluxotypes were measured in minimal medium with glucose as sole carbon source. The WT and Δzwf strains were included in each robotic run as reference strains to check the inter-batch consistency of growth parameters and flux data. In total, 198 cultures were performed from which 1074 OD₆₀₀ points were measured and used for the calculation of growth rates and biomass yields; 1248 samples were collected (1050 filtrates and 198 cell pellets), from which 591 physiological data, i.e., rates of compound consumption or production-were extracted. Further, more than 20,000 isotopic data measured (100 carbon isotopologues from 16 different metabolites by MS analyses and 6 positional labelling profiles for 2 metabolites by NMR for each cultivation). Finally, a cumulated total of 18,612 fluxes were calculated (94 fluxes for each culture). The complete experimental process was completed within 786 h, and required only 194 h, i.e., about 1 h per measured fluxotype-of human intervention or supervision. These values emphasize the practical benefit of the established workflow.

The data collected for the reference strains, i.e., WT and Δzwf strains-in the various runs of the robotic platform did not differ significantly from the data collected in the validation stage. On average, the strains deleted for y-genes grew slower than WT strains but had similar biomass yields (Figure 3A). The rates of glucose consumption and acetate production, measured from exometabolome data, were consistent with these observations, i.e., no significant difference overall. Principal component analysis (PCA) was used to explore the distribution of intracellular flux data across the y-ome (Figure 3B,C). For a large majority of mutant strains (in blue), the flux data-expressed relative to the glucose uptake rate-clustered with those of the WT strains (in green). This result indicates that, for most of the y-genes, deletion has a limited impact on the distribution of intracellular metabolic fluxes. The complete dataset is provided in Supplementary data 4.

Three mutant strains did not cluster with the WT strains (Figure 3B). One of these was the Δzwf strain, demonstrating the value of PCA to discriminate strains based on metabolic fluxes. Two other y-gene mutants, namely Δybjp and ΔydcS, also had significantly different fluxotypes from those of the WT strains and other y-gene mutants. Loading plots (Figure 3C) indicated that the most discriminating fluxes were related to glycolysis (pgi, pgk, eno, pdh), the PPP (zwf, gnd, tk1, tk2), and the glyoxylate shunt (gs1, gs2).

## 2.6. *Glucose Fluxotypes of Δybjp and ΔydcS Strains*

The glucose fluxotypes of the Δybjp and ΔydcS strains are shown in Figure 4. Consistent with the PCA data, significant differences were observed between the two strains and the wild type strain, specifically that both strains had higher glycolytic flux than the wild type. In the Δybjp and ΔydcS strains, 82%, and 76% of glucose was metabolized through pgi, respectively, against 68% in the wild type strain. (i.e., pgi flux in the ΔydcS strain). As a result, the two strains had decreased fluxes through the oxidative branch of the PPP (i.e., zwf and gnd), while the flux through the ED pathway (i.e., EDD) was stable. This trend was even stronger in the ΔydcS strain, for which the data also showed partial reversal of the non-oxidative branch of the PPP (tk1 and tk2 reactions), probably to compensate for the decreased production of pentose-5-P, a key precursor for major biomass components. These
data suggest that the products of the *ybjP* and *ydcS* genes are either directly or indirectly involved in the control of flux in the PPP pathway.

A. Comparison of growth parameters between the *E. coli* Y-ome strains and the WT strain

<table><thead><tr><th></th><th>μ<br>(h<sup>-1</sup>)</th><th>q<sub>glc</sub><br>(mmol.gDW<sup>-1</sup>.h<sup>-1</sup>)</th><th>Biomass yield<br>(gDW<sup>-1</sup>.mmol<sup>-1</sup>)</th><th>q<sub>ace</sub><br>(mmol.gDW<sup>-1</sup>.h<sup>-1</sup>)</th></tr></thead><tbody><tr><td>BW25113<br>(n = 4)</td><td>0.66 ±0.02</td><td>-6.41 ±0.30</td><td>0.10 ±0.01</td><td>2.29 ±0.24</td></tr><tr><td>E. coli Y_ome<br>(average of 180 strains)</td><td>0.57 ±0.06</td><td>-5.81 ±1.02</td><td>0.10 ±0.02</td><td>2.79 ±0.65</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Comparison of the growth parameters and fluxotypes of y-ome strains with those of their parental strain. (A) Comparison of growth parameters between y-ome strains and the BW25113 WT strain (μ: growth rate, *q*<sub>glc</sub>: glucose consumption rate, *Biomass yield*, *q*<sub>ace</sub>: acetate production rate). (B) PCA plots showing the strains with the most different fluxotypes across the y-ome. PC1 and PC2 represent Principal Component 1 and 2, explaining the largest variance between strains according to their flux values-the percentage of variance explained for each principal component is indicated into brackets [36], (C) PCA loading plot with the most impacted fluxes indicated.</figcaption></figure>

Significant alterations of fluxes in the Tricarboxylic Acid Cycle (TCA) cycle and related pathways were observed in the two mutant strains compared to the wild-type strain. Pyruvate dehydrogenase flux was significantly increased in the Δ*ydcS* strain, and even more so in the Δ*ybjP* strain, reaching 119% against 95% in the wild type. This was not accompanied by increased flux through the TCA cycle, but by increased acetate excretion, as well as increased flux through the glyoxylate shunt. Lactate production was observed furthermore in the Δ*ybjP* strain, and this strain had no flux through malate dehydrogenase (i.e., maldh) and increased flux through malic enzyme (i.e., mae). Overall, these data highlight a rerouting of metabolic fluxes around the acetyl-CoA node mainly resulting in the diversion of glycolytic flux from oxidative metabolism to acetate excretion in both strains.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Distribution of metabolic fluxes in ΔydcS and ΔybjP strains. A: Intracellular flux distributions in the two mutant strains and comparison with the WT strain. All fluxes are expressed relative to the glucose uptake rate, set arbitrarily to 100 for each strain. The absolute glucose uptake rates (mmol.gDW⁻¹·h⁻¹) of the three strains are shown in the top left of the figure.</figcaption></figure>

## 2.7. Cofactor Usage in the ΔybjP and ΔydcS Strains

The changes in carbon fluxes were associated with changes in the supply of adenosine-triphosphate (ATP) and redox cofactors by the central carbon metabolism (Figure 5). The ΔydcS strain produced ATP at a slightly lower rate than the WT strain did. Its redox status was significantly impaired. In particular, reduced Nicotinamide Di-Nucleotide Phosphate (NADPH) production was drastically decreased, due to the much lower flux through the PPP and the absence of compensation via the TCA cycle or malic enzyme. Because this strain grows roughly as fast as the WT strain, hence, has similar biosynthetic requirements,
it is likely that the reduced production of NADPH via carbon metabolism (Figure 5) is compensated by increased transhydrogenase activity [37]. The Δ*ybjP* strain had different ATP and redox profiles. First, this strain produces ATP faster than the WT strain due to the significantly increased fluxes through glycolysis and acetate metabolism. Interestingly, the same ATP production profile was observed in the Δ*zwf* strain (Figure 5), but the latter was significantly impaired for NADPH production, while the Δ*ybjP* strain produced NADPH at only a slightly lower rate than the WT strain. This reflects a smaller reduction in PPP flux and an increase in malic enzyme flux, which accounted for 16% of CCM-derived NADPH production in this strain compared to 4% in the WT strain. The latter increase might be the result of a mechanism used to compensate for the decreased NADPH production via the oxidative PPP. In addition, NADH production was significantly higher in the Δ*ybjP* mutant than in the WT strain, mainly due to a large increase in glycolytic production of NADH. This resulted in the production of lactate, suggesting a global imbalance in redox metabolism in this strain.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Production of NADPH, NADH/FADH₂ and ATP in the central carbon metabolism. The rates were calculated from the mean growth rates and flux distributions, and are expressed in absolute values (mmol.gDW⁻¹·h⁻¹) ± SD. Calculations are described in the materials and methods.</figcaption></figure>

## 2.8. *Scope and Quality of High-Throughput Fuxomics Investigations*

Given that a specific objective of this work was to deploy high-resolution fluxomics in a HT manner, we evaluated the quality of the fluxome data, obtained here, in comparison with existing data in the literature. High-throughput investigations must find the best compromise between throughput and information level (and quality). Therefore, the volume and quality of the flux data have to be appreciated in terms of the actual size of the investigation, i.e., the total number of flux maps generated. To this aim, we used the four criteria defined below:

(i) *fluxome resolution*, which corresponds to the total number of fluxes that can be calculated from the experimental dataset, and is a measure of the coverage of the flux space in the investigated metabolic network. This resolution is generally low (e.g., 23) in high-throughput studies, and high (71–84) in low-throughput ones (Table 1). In this work, we were able to measure 94 fluxes per measured fluxotype, which is the highest fluxome resolution value among all ¹³C-fluxomics studies reported in Table 1.

(ii) *isotopic resolutive power* of the flux map, defined as the ratio between the number of isotopic data and the number of fluxes calculated in the network. This index is assumed to reflect data redundancy in the flux estimates, hence, higher precision in the calculated flux values. Low-throughput investigations of the *E. coli* fluxome have isotopic resolutive power in the range of 1 to 8 (Table 1), though very high isotopic resolutive power (up to 17.6), has been obtained by performing multiple (e.g., parallel) labelling experiments to calculate a single flux map [38–41]. In the present work, the isotopic resolutive power of the flux maps was 1.13 (106 isotopic
data for 94 fluxes), which is comparable to the values obtained in low-throughput investigations (Table 1).

(iii) *total flux dimension*, which is equal to the number of fluxotypes multiplied by fluxome resolution, and is an indication of the scale of the fluxomics investigation. In the literature, this index ranges between 71 and 4370 (Table 1). In this study, the total flux dimension index was 18,612 (94 fluxes × 198 fluxotypes), which reflects both the high-throughput and high-resolution character of the analysis. Interestingly, the number of fluxotypes, generated in this study, is similar to the number reported by Haverkorn van Rijsewijk et al. [33] (198 versus 190, respectively) but the total flux dimension of this work is higher (18612 versus 4370, respectively) because of the higher number of fluxes per fluxotype (94 versus 23, respectively).

(iv) *overall flux precision index*, which is the median RSD of all (free) fluxes across the entire flux dataset. Only free fluxes are considered because the others, e.g., most biosynthetic fluxes are determined or given as constraints to the model. The average value of this index was 32% in the present study, compared with 0.4–23% in low-throughput studies, and 14–253% in HT studies (Table 1). The value, obtained in this work, appears to strike a good compromise between throughput and resolution. Moreover, when measured for specific pathways (*overall pathway-specific precision*), it appears that the level of precision achieved for most fluxes is drastically improved, notably for those in the glycolysis, PPP and ED pathways.

<table><caption>**Table 1.** Evaluation of data size and quality in HT fluxomics investigations. The overall flux precision is the median RSD (relative standard deviation) calculated for all fluxes measured in each of the studied conditions, as given in the supplementary data of the cited studies. The overall pathway precision was calculated similarly by grouping the RSDs according to the metabolic pathways considered. The flux RSDs considered in these calculations were obtained by sensitivity analysis (S.A.) or by calculating classical standard deviations on the mean values of fluxes obtained when biological replicates were produced (as in this work).</caption><thead><tr><th></th><th>Leighty et al., 2013 [37]</th><th>Crown et al., 2015 [38]</th><th>Long &amp; Antoniewicz 2019 [35]</th><th>Millard et al., 2014 [30]</th><th>Long et al., 2016 [39]</th><th>Heux et al., 2014 [22]</th><th>Haverkorn van Rijsewijk et al., 2011 [31]</th><th>This Paper</th></tr></thead><tbody><tr><td>E. coli strains</td><td>MG1655</td><td>MG1655</td><td>BW25113 ΔtpiA mutant</td><td>MG1655</td><td>BW25113</td><td>MG1655</td><td>BW25113 + Keio mutants</td><td>MG1655 + BW25113 + keio mutants</td></tr><tr><td>Number of measured fluxotypes</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>20</td><td>190</td><td>198</td></tr><tr><td>Number of different label input(s) by strain or condition</td><td>6</td><td>14</td><td>2</td><td>1</td><td>1</td><td>2</td><td>2</td><td>1</td></tr><tr><td>Fluxome resolution = number of (net) fluxes by fluxotype</td><td>71</td><td>71</td><td>75</td><td>84</td><td>71</td><td>23</td><td>23</td><td>94</td></tr><tr><td>Isotopic resolutive power = number of isotopic data/number of calculated fluxes</td><td>7.52</td><td>17.55</td><td>5.07</td><td>1.4</td><td>1.94</td><td>2.57</td><td>1.61</td><td>1.13</td></tr><tr><td>Total flux dimension = number of fluxes per fluxotype × number of fluxotypes</td><td>71</td><td>71</td><td>75</td><td>84</td><td>142</td><td>460</td><td>4370</td><td>18,612</td></tr></tbody></table>
<table><caption>**Table 1.** *Cont.*</caption><thead><tr><th></th><th>Leighty et al., 2013 [37]</th><th>Crown et al., 2015 [38]</th><th>Long & Antoniewicz 2019 [35]</th><th>Millard et al., 2014 [30]</th><th>Long et al., 2016 [39]</th><th>Heux et al., 2014 [22]</th><th>Haverkorn van Rijsewijk et al., 2011 [31]</th><th>This Paper</th></tr></thead><tbody><tr><td>Global flux precision = median RSD over the global dataset</td><td>12%</td><td>15%</td><td>7.80%</td><td>19%</td><td>23%</td><td>285%</td><td>14%</td><td>32%</td></tr><tr><td>Global pathway precision = median RSD within specific pathways</td><td>S.A. (n = 1)</td><td>S.A. (n = 1)</td><td>S.A. (n = 1)</td><td>S.A. (n = 1)</td><td>S.A. (n = 1)</td><td>S.A. (n = 23)</td><td>S.A. (n = 23)</td><td>S.A. (n = 198)</td><td>biological replicates (n = 20)</td></tr><tr><td>Glycolysis</td><td>1%</td><td>2%</td><td>4.20%</td><td>3%</td><td>3.50%</td><td>553%</td><td>7%</td><td>3%</td><td>1%</td></tr><tr><td>ppp + edp</td><td>14%</td><td>11%</td><td>17%</td><td>131%</td><td>10%</td><td>41%</td><td>25%</td><td>24%</td><td>7%</td></tr><tr><td>tca + gs</td><td>11%</td><td>19%</td><td>4.90%</td><td>20%</td><td>31%</td><td>795%</td><td>19%</td><td>40%</td><td>21%</td></tr><tr><td>anaplerosis</td><td>122%</td><td>144%</td><td>87%</td><td>13%</td><td>1822%</td><td>7%</td><td>15%</td><td>42%</td><td>20%</td></tr><tr><td>output fluxes</td><td>12%</td><td>25%</td><td>4.7%</td><td>15%</td><td>42.4%</td><td>426%</td><td>0%</td><td>10%</td><td>14%</td></tr></tbody></table>

### 3. Discussion

This work introduces significant improvements in the parallelization, automation and implementation of fluxomics in large-scale investigations. Altogether, these improvements offer significant gains in throughput, robustness, and data quality, while retaining a high level of biological information by collecting high resolution fluxomics data. This optimized complete HT fluxomics workflow allowed 198 isotopically-resolved flux maps to be generated in a total (effective) experimental time of less than 33 days. This required a total human time of 194 h, or about 1 h per flux map. In the context of high-resolution fluxomics, the working volume of the bioreactors, 15 mL, seems an appropriate compromise between the low volumes required for HT investigations, and culture volumes, which are large enough for detailed biological or biochemical analysis. Indeed, this working volume allowed us to systematically collect all relevant physiological and biochemical data in a single bioreactor under fully controlled growth conditions, thereby decreasing biological variability and increasing data consistency. Further applications of this workflow to other experimental designs, e.g., parallel labeling experiments, other questions or other biological models (bacterial communities, eukaryotic cells, etc.) should be facilitated by the flexibility of the setup, which is amenable to a variety of sample types (medium, cells, total broth), sampling methods, metabolite extraction processes (intra- or extra-cellular metabolites, biomass), and analytical techniques (basically any targeted/untargeted LC/GC-MS method or NMR sequence). This flexibility allows metabolic networks of various size and complexity, as well as other biological models, to be investigated. There is still scope for improvement in almost all stages of the workflow—from sampling to flux calculations and sensitivity analysis—to increase throughput and reduce experimental times and human involvement. A major challenge in the application of HT fluxomics is the level of precision on flux data that can be achieved at high or very high throughput. In this regard, the capability of the robotic platform to enable a high(er) number of biological replicate experiments to be carried out certainly contributes to increase the reliability of the data.

¹³C-Fluxomics encompasses different approaches depending on the number of isotopic data collected and number of fluxes calculated [21,23,24]. These approaches include *isotopic profiling* (i.e., the statistical analysis of large isotopic datasets to discriminate variants), *targeted 13C-fluxomics* (in which a limited number of isotopic data are collected to measure fluxes in a few reactions or pathways), and *global ¹³C-fluxomics* (the measurement of fluxes across complete metabolic networks from large isotopic datasets). The former two approaches are low-resolution methods because they provide no (isotopic profiling), or
limited (targeted fluxomics), flux information, though they are fully relevant for their respective purposes. High-resolution fluxomics refer to global ¹³C-fluxomics, but the number of measured fluxes can vary significantly between studies (Table 1). HT fluxomics approaches have been discussed, in terms of methodological challenges and biological value, but so far, there has been no method proposed to evaluate the actual scale and quality of HT fluxomics investigations. This work introduces novel metrics for the actual scale and quality of HT fluxomics investigations. At the level of individual fluxotypes, key features of these studies, include the total number of calculated fluxes (*fluxome resolution*), the number of isotopic data, collected with respect to the number of fluxes to be measured (*isotopic resolutive power*), and the resulting precision in the calculated flux values (confidence intervals on flux values). At the level of the complete HT fluxomics investigation, three indices are introduced to measure the actual scale (*total flux dimension*) and quality (*overall flux precision* and *overall pathway-specific precision*) of HT fluxomics investigations. Since HT studies aim to find the best compromise between throughput and information content (and quality), the number and quality of the resulting flux data have to be appreciated, in terms of the actual size of the investigation, i.e., the total number of flux maps generated. HT fluxomics investigations are characterized by high *total flux dimension*, low isotopic resolutive power and low flux precision [22]. On the contrary, high-resolution fluxomics are characterized by low total flux dimensions, high isotopic resolutive power and high flux precision [39–41]. Finally, methods developed for the large-scale application of high-resolution fluxomics, such as the one described here, should aim to achieve good isotopic resolutive power and flux precision with high total flux dimension [33]. Indeed, we could apply High-resolution fluxomics (fluxome resolution of 94) at a high-throughput level (198 fluxotypes collected throughout the study) without loss of quality. It is important to note that these metrics are objective measures of the size and characteristics of fluxomics investigations; they do not quantify the biological value of the data.

The main objective of this work was to investigate whether the deletion of y-genes led to any changes in the distribution of metabolic fluxes in the central metabolism, i.e., in glucose metabolic fluxotypes. Deletion of the part of the y-ome, considered in this study, was associated with an overall tendency toward slower growth, but for the vast majority of the 180 strains investigated, the distribution of metabolic fluxes was remarkably similar to the one measured for the WT strain. Given that *E. coli* expressed the products of all the investigated y-genes during growth on glucose, there are several possible explanations for these results. (i) The protein is expressed, but is not functional under the considered culture conditions, hence, its absence has no impact in the deletion mutant. (ii) The protein has a functional role under these culture conditions, but this function is not related to metabolism, hence, its absence does not result in a metabolic phenotype. (iii) The protein has a metabolic function, as an enzyme or as a metabolic regulator in these culture conditions. In this case, the absence of a particular flux phenotype can be explained by the gene product having a quantitatively minor role, at least for the part of metabolism that was investigated, or by the onset of compensatory mechanisms that efficiently counterbalance the effects of the gene product's absence. Further investigations, including on fluxotypes under other experimental conditions are required to answer these questions. However, regardless of the specific mechanisms involved for each y-gene, the data reported here highlight the robustness of *E. coli*'s central metabolism to the absence of all these gene products, whatever their roles, and to the potential impairments caused by their absence.

In this work, we were able to provide highly detailed functional information on the metabolic impact of y-gene deletion. Two strains, *ΔybjP* and *ΔydcS* showed interesting glucose fluxotypes, indicating that the products of the genes are involved in metabolism during growth on glucose as sole carbon source. Information on these genes is scarce and elusive. The gene *ydcS* was annotated as a putative putrescine ATP-binding Cassette transporter in several databases and has also been proposed to encode poly-3-hydroxybutyrate synthase activity [42]. The gene *ybjP* has been predicted to be a lipoprotein [43]. Interestingly, the metabolic impact of both these genes is global, since most of the pathways of the central car-
bon metabolism are affected. Indeed, both mutant strains show significant perturbations in the PPP and in the TCA cycle and related pathways. The deletions are also both associated with energy and redox status modifications, but with opposite consequences. The $\Delta ybjP$ strain has significantly impaired ATP production and reduced cofactors within the central carbon metabolism, while the $\Delta ybjP$ has higher production levels than observed in the WT strain. The fact that the metabolic effects are so wide-ranging suggests that the products of these genes affect central metabolic parameters, such as the redox status or that they are global regulators of metabolism, either directly or indirectly. The dispensability of the two y-genes during growth on various carbon sources has recently been documented [12]. The $\Delta ydcS$ strain shows reduced growth on certain sugars (mainly mannose, glucosamine, fructose) and two organic acids (alpha-ketoglutarate and acetate). The $\Delta ybjP$ strain shows reduced growth on a few carbon sources, mainly including pyruvate, ribose, sorbitol and galacturonate. Deletion of these y-genes has no effect [12] or only a slight effect (this work) on growth with glucose as sole carbon source, in spite of significant metabolic changes in the two mutant strains. These observations emphasize that gene dispensability does not mean the gene lacks a metabolic function, but rather that the metabolism has successfully adapted to the absence of the gene product [4]. This observation stresses the need to measure the actual metabolic phenotype of mutants, even if the gene deletion is silent in terms of growth. In this context, methods like High-resolution fluxotyping are essential in achieving both the throughput and level of metabolic information required to investigate the role of y-genes in the metabolism of diverse carbon sources.

To conclude, this work shows that high-resolution fluxomics is amenable to high-throughput investigations and can provide detailed information on metabolic phenotypes. Our results also highlight how it can reveal the metabolic impact of gene deletion, even if the gene has no known function or growth phenotype.

# 4. Materials and Methods

## 4.1. Bacteria Strains and Cultivation Conditions

E. coli BW 25,113 and the derived strains used in this study (listed in Supplementary data 1) were taken from the Keio collection (Baba et al., 2006). One hundred and eighty y-gene mutants were selected from the original glycerol stock of 3985 single-gene deletion mutants available in this collection (E. coli BW 25,113 strains mutants; more details directly in paragraph 2.1 of the paper). The E. coli strain MG1655 was also used. All the selected strains were first cultivated overnight in Luria-Bertani (LB) medium (10 g/L tryptone, 5 g/L yeast extract and 10 g/L NaCl) with kanamycine (25 µg/mL) at 37 °C and then stored in glycerol stocks.

To perform the experiments, the glycerol stocks of relevant strains were used to inoculate liquid LB medium in microplates. The LB cultures were used to inoculate precultures in minimal synthetic medium containing 17.4 g·L⁻¹ Na₂HPO₄, 12H₂O, 3.03 g·L⁻¹ KH₂PO₄, 0.51 g·L⁻¹ NaCl, 2.04 g·L⁻¹ NH₄Cl, 0.49 g·L⁻¹ MgSO₄, 4.38 mg·L⁻¹ CaCl₂, 15 mg·L⁻¹ Na₂EDTA 2H₂O, 4.5 mg/L ZnSO₄ 7H₂O, 0.3 mg·L⁻¹ CoCl₂ 6H₂O, 1 mg·L⁻¹ MnCl₂ 4H₂O, 1 mg·L⁻¹ H₃BO₃, 0.4 mg·L⁻¹ Na₂MoO₄ 2H₂O, 3 mg·L⁻¹ FeSO₄ 7H₂O, 0.3 mg·L⁻¹ CuSO₄ 5H₂O, 0.1 g·L⁻¹ thiamine and 3 g·L⁻¹ glucose. For the ¹³C-labeling experiments, unlabeled glucose was replaced by the same concentration of a mixture of 80% [1-¹³C₁]-D-glucose + 20% [U-¹³C₆]-D-glucose. To minimize sources of unlabeled carbon atoms from the first cultivation steps in subsequent experiments, the cells were inoculated at a starting OD between 0.08 and 0.12.

## 4.2. Robotic Platforms for Culture, Sampling and Sample Preparation

Two robotic platforms were used to parallelize the cell cultures, ¹³C-labelling experiments, sampling of labeled metabolites, and sample preparation for NMR and MS analyses. The first system allowed 48 ¹³C-labelling experiments to be performed automatically in parallel in 15 mL bioreactors under controlled growth conditions, with automated collection of labelled samples (of biomass or culture medium) at defined culture times or optical
densities. The robot and its operation are described in detail by Heux et al. [22]. A second robotic workstation (Freedom EVO 200, Tecan, Männedorf, Switzerland) was designed and assembled to fully automate and parallelize the final preparation of labeled biological samples for the analysis of isotopic profiles by NMR or mass spectrometry. This device was used to handle the different containers used (NMR tubes, vials, multiwell-plates, etc.), in order to perform dilutions, add standards, take aliquots, and manage the samples.

The ¹³C-labeling experiments were carried out in batches of 48 parallel cell cultures using the two robotic platforms described above. The input label was optimized for *E. coli*'s metabolic network using IsoDesign (1.2.1) [24], and consisted of a mixture of 80% [1-¹³C₁]-D-glucose and 20% [U-¹³C₆]-D-glucose. All cultures were performed in 15 mL reaction vessels, at 37 °C, pH = 7, a stirring speed of 2300 rpm and with 5 L.min⁻¹ of compressed air fed into the culture module. Four to eight supernatant samples were collected in each reactor throughout the culture to analyze the exometabolome and to calculate the rates of substrate consumption and end-product production. The biomass was automatically sampled once the OD₆₀₀nm had reached 1.2, indicating that the metabolic and isotopic steady-state had been achieved: 4 mL samples of culture were pelleted by automated centrifugation (5 min/4410 g), manually hydrolyzed with 250 µL 6 N HCl for 15 h at 110 °C and washed twice in 1 mL D₂O by rotary evaporation (Büchi Labortechnik AG, CH, Flawil, Switzerland) between each washing step. Aliquots (10 µL) of biomass hydrolysate supernatant were then collected and transferred into 96 well plates, diluted with 990 µL pure H₂O and transferred to vials for LC-MS analysis. Samples of biomass hydrolysate supernatant (150 µL) were also mixed with 50 µL of TSP d4 (4 mM in D₂O) and 150 µL of each mixture was then transferred into 3 mm NMR tubes.

### 4.3. Isotopic Profiling of Proteinogenic Amino-Acids

The incorporation of ¹³C-label into proteinogenic amino acids (listed in Supplementary data 6) was analyzed by liquid chromatography–mass spectrometry, using an Ultimate 3000 HPLC system (Dionex, Sunnyvale, CA, USA) coupled to an LTQ Orbitrap Velos mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) equipped with a heated electrospray ionization probe described in detail by Heuillet et al., 2018 [44]. Full scan HRMS analyses were performed in positive FTMS mode at a resolution of 60,000 (at 400 m/z), using the following source parameters: Capillary temperature, 275 °C; source heater temperature, 250 °C; sheath gas flow rate, 45 a.u. (arbitrary unit); auxiliary gas flow rate, 20 a.u.; S-Lens RF level, 40%; source voltage, 5 kV. Metabolites were identified by extracting the exact mass with a tolerance of 5 ppm. The raw MS isotopic profiles of proteinogenic amino acids were then quantified using Tracefinder (Thermo Fisher Scientific, Waltham, MA, USA). The isotopic profiles (Carbon Isotopologue Distributions) were obtained after correcting for natural isotopic abundances using IsoCor 1.2 [27,28] (https://github.com/MetaSys-LISBP/IsoCor (accessed on 20 May 2020)). The raw MS data are available from Metabolights (accession number MTBLS2188).

### 4.4. NMR Analysis of Extracellular Medium

Culture supernatants were analyzed by 1D-¹H NMR on a Bruker Avance III 800MHz spectrometer (Bruker BioSpin, Rheinstetten, Germany) page equipped with a 5 mm CQPI cryoprobe at 280 K. To precisely quantify the extracellular compounds, the 1D-¹H NMR data were recorded after a 30° presaturation pulse (zgpr30), with a relaxation delay of 7 s. NMR spectra were processed using Topspin 3.5pl6 (Bruker BioSpin, Rheinstetten, Germany). The raw NMR data are available from Metabolights (accession number MTBLS2188).

The rate of glucose consumption and the rates of acetate and lactate production were quantified by analyzing the four-to-eight samples collected during the cell culture period. The ¹³C labeling profiles of acetate and lactate, including positional information on label incorporation, were also measured for the flux calculations.
### 4.5. Growth Parameters

The growth parameters calculated from the experimental datasets for each culture were the growth rate, glucose uptake rate, acetate and lactate production rates and biomass yield. The growth rate was calculated from OD values measured during the culture period. The rates of glucose consumption and of acetate and lactate production were calculated from the exometabolome data. All calculations were performed using Physiofit (v0.9) [45], a software designed to determine growth parameters by fitting time-course data. The software includes a batch calculation mode, which allows calculations on large series of datasets, and which was used to calculate the growth parameters for all considered experimental conditions. The biomass yields were calculated using the conversion factor 0.378 gDW/OD600nm unit.

Extracellular fluxes were determined from the rates of disappearance (or appearance) of substrates and products, in the culture supernatants, as measured by NMR.

### 4.6. Flux Calculation and Visualization

Fluxes were calculated using the software influx si 4.1 [30] (https://metasys.insa-toulouse.fr/software/influx/ (accessed on 20 May 2020)), including the mass balances and carbon atom transitions of the biochemical reaction network. The metabolic network contained the main pathways of *E. coli*'s central metabolism: glycolysis (EMP), the pentose phosphate pathway (PPP), the Entner-Doudoroff (ED) pathway, the tricarboxylic acid cycle (TCA), and anaplerotic reactions, the glyoxylate shunt, and the reactions for amino acid biosynthesis [32]. Intracellular fluxes were estimated from measurements of extracellular fluxes and from the ¹³C-labelling patterns of metabolites using appropriate mathematical models of glucose metabolism in *E. coli* [32]. Labeling data were collected from intracellular metabolites by IC-MS/MS and from metabolic end-products by 1D ¹H NMR, as detailed above. The fluxes were normalized to the rate of substrate uptake, which was arbitrarily set at 100.

Metabolic fluxes were calculated for each of the 198 experimental conditions considered in the study. For each culture condition, the information required for the flux calculation (i.e., metabolic network, isotopic data, etc.) was written in a specific FTBL file [16], which was then submitted to influx si.

The FTBL model was converted by an influx si module (ftbl2metxml.py) into Systems Biology Markup Language (SBML, in xml format), readable by the online application Metexplore (http://metexplore.toulouse.inra.fr (accessed on 20 May 2020). Metexplore was used to visualize flux maps for strains of interest [46].

### 4.7. Sensitivity Analysis

Confidence intervals for the calculated fluxes were determined using a Monte-Carlo approach in which 100 independent optimization runs were performed on datasets with noise added in proportion to standard measurement errors. The isotopic data and metabolic fluxes for each independent biological replicate are provided in the Supplemental data (dataset S1).

### 4.8. Statistical Analysis of Flux Datasets

Flux maps were compared by PCA using the collaborative portal Workflow for Metabolomics (W4M) [36]. The results were visualized as boxplots and figures were generated using the software R.

### 4.9. Meta-Data Management

The data and meta-data, generated throughout the workflow, were stored, managed, and processed as described below.

Culture and sampling data are related to the parallel cultures of *E. coli* strains with ¹³C-labelled compounds and to the collection of samples: Settings and data from the robotic culture platform (Tecan Evo 200) were automatically and separately recorded for each
of the 48 bioreactors: pH values, O2 values, OD₆₀₀nm values, volumes, sampling times, and tube locations and associated barcodes. The data were linked and stored in an SQL database, making them traceable to each robotic run. The most important elements of this database are listed in Supplementary data 3 (mu, batch number, batch position, biomass (mg/L) at sampling time). A log file was also automatically generated for each robotic run, detailing errors and software information.

**Sample preparation data:** All the steps performed on the sample preparation platform (Tecan Evo 200) were automatically recorded in an SQL database: Volumes, tube locations on the worktable, associated barcodes and sampling times. The data were visualized and previous runs re-examined using the robot's traceability software, and a log file was automatically generated for each run with errors and software information.

**Analytical data:** All mass spectrometry and NMR data were saved on a server to guarantee sample traceability. This included the raw MS data, the MS data processed with Tracefinder, the MS data processed with Isocor, the raw NMR data, the NMR data processed with TopSpin, and the NMR data processed with Physiofit. Corrected CID values, extracellular fluxes and the positional enrichment of extracellular metabolites are listed in Supplementary data 3.

**Flux calculations:** All flux calculation files were stored on a server. This included the input (FTBL) files for influx si, and the output (KVH) files. The most important elements are included in Supplementary data 3 (chi-square tests of the fits, fit cost functions and net fluxes).

**Statistical analysis and flux mapping:** All the data used and generated for statistical analysis and flux mapping were stored on a server.

### 4.10. Calculation of Cofactor Production and Consumption Rates

The production rates (mmol.gDW⁻¹·h⁻¹) of NADPH, NADH/FADH₂, and ATP in the central carbon metabolism were calculated as the sum of the estimated fluxes of the reactions that are expected to produce (positive value) or consume (negative value) the cofactors. Values were averaged over the biological replicates. NADPH production: Pentose-Phosphate Pathway (PPP) = glucose 6-phosphate dehydrogenase (zwf) + phosphogluconate dehydrogenase (gnd); TCA cycle = isocitrate dehydrogenase (idh); malic enzyme (mae). NADH production: glycolysis = glyceraldehyde-3-phosphate dehydrogenase (flux equal to pgk); Oxidative metabolism = pyruvate dehydrogenase (pdh) + TCA cycle enzymes, including alpha-ketoglutarate dehydrogenase (akgdh) + malate dehydrogenase (maldh); malic enzyme (mae) and biomass formation. NADH consumption: lactate dehydrogenase (out Lac). FADH₂ production: TCA cycle = succinate dehydrogenase (flux equal to fum). ATP production: glycolysis = phosphoglycerate kinase (pgk) + pyruvate kinase (pyk); TCA cycle = succinyl-CoA synthetase (assimilated to akgdh); acetate metabolism = acetate kinase (ack, assimilated to out_Ac). ATP consumption: glycolysis = glucose PTS system (Glucupt_U + Glucupt_1) (assuming 1 consumed PEP is equivalent to 1 consumed) + 6-phosphofructokinase (pfk); ppc (equivalent to pep carboxylase (forward flux) and pep carboxykinase (reverse flux).

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/metabo11050271/s1, Supplementary data 1: List of the 180 *E. coli* BW 25113 strains selected for the fluxomics experiments. Supplementary data 2: Metabolic model considered in the flux calculations. Supplementary data 3: Complete dataset obtained for all the fluxomic experiments performed in the study. Supplementary data 4: Calculation of ATP and redox fluxes. Supplementary data 5: Isotopic data collected by LC HRMS and NMR.

**Author Contributions:** Methodology, C.B., E.C., P.M., B.E., L.P., F.B., S.H. and J.-C.P.; software, C.B., E.C., P.M., M.G. and S.S.; validation, C.B., E.C., P.M., B.E., L.P., F.B., S.H. and J.-C.P.; formal analysis, C.B and E.C.; investigation, C.B., E.C.; resources, C.B., E.C., P.M., B.E., M.D., M.H., H.K., L.G., N.B., L.P., T.P. and S.H.; data curation, C.B. and E.C.; writing—original draft preparation, C.B., E.C.; writing—review and editing, C.B., E.C. and J.-C.P.; visualization, C.B., M.C. and E.C.; supervision,
L.P., F.B., S.H. and J.-C.P.; project administration, L.P., F.B. and J.-C.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** MetaToul-MetaboHUB (Metabolomics and Fluxomics facilities, Toulouse, France, http://www.metatoul.fr (accessed on 24 March 2021) is part of MetaboHUB (The French National Infrastructure for Metabolomics and Fluxomics, www.metabohub.fr (accessed on 24 March 2021) funded by ANR grant MetaboHUB-ANR-11-INBS-0010. MetaToul-MetaboHUB is also supported by grants from the Région Occitanie, the European Regional Development Fund, the SICOVAL, the Infrastructures en Biologie Sante et Agronomie (IBiSa, France), the Centre National de la Recherche Scientifique (CNRS) and the Institut National pour la Recherche en Alimentation, Agronomie et Environnement (INRAE).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Restrictions apply to the availability of these data. Data was obtained from Metabolights and are available at www.ebi.ac.uk/metabolights/MTBLS2188 (accessed on 22 April 2021) with the permission of Metabolights.

**Acknowledgments:** The authors thank Workflow4Metabolomics (W4M) the consortium, which is supported by MetaboHUB [MetaboHUB-ANR-11-INBS-0010, 2011] and the French Institute of Bioinformatics: IFB [ANR-11-INBS-0013236], for online access to statistical processing tools. The authors also thank Marie Tremblay Franco (INRAE Toulouse, France) for help with the statistical analyses. J.-C.P. is grateful to the INSERM for funding his temporary full-time researcher position.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Hanson, A.D.; Pribat, A.; Waller, J.C.; de Crécy-lagard, V. 'Unknown' Proteins and 'Orphan' Enzymes: The Missing Half of the Engineering Parts List and How to Find It. *Biochem. J.* **2009**, *425*, 1-11. [CrossRef]

2. Hu, P.; Janga, S.C.; Babu, M.; Díaz-Mejía, J.J.; Butland, G.; Yang, W.; Pogoutse, O.; Guo, X.; Phanse, S.; Wong, P.; et al. Global Functional Atlas of *Escherichia coli* Encompassing Previously Uncharacterized Proteins. *PLoS Biol.* **2009**, *7*, e1000096. [CrossRef] [PubMed]

3. Ghatak, S.; King, Z.A.; Sastry, A.; Palsson, B.O. The y-ome Defines the Thirty-Four Percent of *Escherichia coli* Genes That Lack Experimental Evidence of Function. *bioRxiv* **2018**, 328591. [CrossRef]

4. Raamsdonk, L.M.; Teusink, B.; Broadhurst, D.; Zhang, N.; Hayes, A.; Walsh, M.C.; Berden, J.A.; Brindle, K.M.; Kell, D.B.; Rowland, J.J.; et al. A Functional Genomics Strategy That Uses Metabolome Data to Reveal the Phenotype of Silent Mutations. *Nat. Biotechnol.* **2001**, *19*, 45-50. [CrossRef] [PubMed]

5. Fuhrer, T.; Zampieri, M.; Sévin, D.C.; Sauer, U.; Zamboni, N. Genomewide Landscape of Gene-Metabolome Associations in *Escherichia Coli. Mol. Syst. Biol.* **2017**, *13*, 907. [CrossRef] [PubMed]

6. Sévin, D.C.; Fuhrer, T.; Zamboni, N.; Sauer, U. Nontargeted in Vitro Metabolomics for High-Throughput Identification of Novel Enzymes in *Escherichia Coli. Nat. Methods* **2017**, *14*, 187-194. [CrossRef]

7. Blaby-Haas, C.E.; de Crécy-Lagard, V. Mining High-Throughput Experimental Data to Link Gene and Function. *Trends Biotechnol.* **2011**, *29*, 174-182. [CrossRef]

8. Long, C.P.; Antoniewicz, M.R. Metabolic Flux Analysis of *Escherichia coli* Knockouts: Lessons from the Keio Collection and Future Outlook. *Curr. Opin. Biotechnol.* **2014**, *28*, 127-133. [CrossRef]

9. Baba, T.; Ara, T.; Hasegawa, M.; Takai, Y.; Okumura, Y.; Baba, M.; Datsenko, K.A.; Tomita, M.; Wanner, B.L.; Mori, H. Construction of *Escherichia coli* K-12 In-frame, Single-gene Knockout Mutants: The Keio Collection. *Mol. Syst. Biol.* **2006**, *2*, 2006.0008. [CrossRef]

10. Joyce, A.R.; Reed, J.L.; White, A.; Edwards, R.; Osterman, A.; Baba, T.; Mori, H.; Lesely, S.A.; Palsson, B.Ø.; Agarwalla, S. Experimental and Computational Assessment of Conditionally Essential Genes in *Escherichia coli. J. Bacteriol.* **2006**, *188*, 8259-8271. [CrossRef]

11. Guzmán, G.I.; Olson, C.A.; Hefner, Y.; Phaneuf, P.V.; Catoiu, E.; Crepaldi, L.B.; Micas, L.G.; Palsson, B.O.; Feist, A.M. Reframing Gene Essentiality in Terms of Adaptive Flexibility. *BMC Syst. Biol.* **2018**, *12*, 1-20. [CrossRef]

12. Tong, M.; French, S.; El Zahed, S.S.; Ong, W.K.; Karp, P.D.; Brown, E.D. Gene Dispensability in *Escherichia coli* Grown in Thirty Different Carbon Environments. *mBio* **2020**, *11*. [CrossRef]

13. Warner, J.R.; Reeder, P.J.; Karimpour-Fard, A.; Woodruff, L.B.A.; Gill, R.T. Rapid Profiling of a Microbial Genome Using Mixtures of Barcoded Oligonucleotides. *Nat. Biotechnol.* **2010**, *28*, 856-862. [CrossRef] [PubMed]
14. Liu, X.; Gallay, C.; Kjos, M.; Domenech, A.; Slager, J.; van Kessel, S.P.; Knoops, K.; Sorg, R.A.; Zhang, J.-R.; Veening, J.-W. *High-Throughput CRISPRi Phenotyping Identifies New Essential Genes in Streptococcus Pneumoniae*. *Mol. Syst. Biol.* **2017**, *13*, 931. [CrossRef] [PubMed]

15. Price, M.N.; Wetmore, K.M.; Waters, R.J.; Callaghan, M.; Ray, J.; Liu, H.; Kuehl, J.V.; Melnyk, R.A.; Lamson, J.S.; Suh, Y.; et al. *Mutant Phenotypes for Thousands of Bacterial Genes of Unknown Function*. *Nature* **2018**, *557*, 503–509. [CrossRef]

16. Wiechert, W. 13C Metabolic Flux Analysis. *Metab. Eng.* **2001**, *3*, 195–206. [CrossRef] [PubMed]

17. Sauer, U. Metabolic Networks in Motion: 13C-based Flux Analysis. *Mol. Syst. Biol.* **2006**, *2*, 62. [CrossRef]

18. Wittmann, C.; Portais, J.-C. Metabolic Flux Analysis. In *Metabolomics in Practice*; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2013; pp. 285–312, ISBN 978-3-527-65586-1.

19. Ellis, D.I.; Goodacre, R. Metabolomics-Assisted Synthetic Biology. *Curr. Opin. Biotechnol.* **2012**, *23*, 22–28. [CrossRef]

20. Fan, J.; Ye, J.; Kamphorst, J.J.; Shlomi, T.; Thompson, C.B.; Rabinowitz, J.D. Quantitative Flux Analysis Reveals Folate-Dependent NADPH Production. *Nature* **2014**, *510*, 298–302. [CrossRef]

21. Heux, S.; Bergès, C.; Millard, P.; Portais, J.-C.; Létisse, F. Recent Advances in High-Throughput 13C-Fluxomics. *Curr. Opin. Biotechnol.* **2017**, *43*, 104–109. [CrossRef]

22. Heux, S.; Poinot, J.; Massou, S.; Sokol, S.; Portais, J.-C. A Novel Platform for Automated High-Throughput Fluxome Profiling of Metabolic Variants. *Metab. Eng.* **2014**, *25*, 8–19. [CrossRef] [PubMed]

23. Sauer, U. High-Throughput Phenomics: Experimental Methods for Mapping Fluxomes. *Curr. Opin. Biotechnol.* **2004**, *15*, 58–63. [CrossRef] [PubMed]

24. Moxley, J.F.; Jewett, M.C.; Antoniewicz, M.R.; Villas-Boas, S.G.; Alper, H.; Wheeler, R.T.; Tong, L.; Hinnebusch, A.G.; Ideker, T.; Nielsen, J.; et al. Linking High-Resolution Metabolic Flux Phenotypes and Transcriptional Regulation in Yeast Modulated by the Global Regulator Gcn4p. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 6477–6482. [CrossRef]

25. Schmidt, A.; Kochanowski, K.; Vedelaar, S.; Ahrné, E.; Volkmer, B.; Callipo, L.; Knoops, K.; Bauer, M.; Aebersold, R.; Heinemann, M. The Quantitative and Condition-Dependent *Escherichia coli* Proteome. *Nat. Biotechnol.* **2016**, *34*, 104–110. [CrossRef] [PubMed]

26. Millard, P.; Sokol, S.; Letisse, F.; Portais, J.-C. IsoDesign: A Software for Optimizing the Design of 13C-Metabolic Flux Analysis Experiments. *Biotechnol. Bioeng.* **2014**, *111*, 202–208. [CrossRef]

27. Millard, P.; Letisse, F.; Sokol, S.; Portais, J.-C. IsoCor: Correcting MS Data in Isotope Labeling Experiments. *Bioinforma. Oxf. Engl.* **2012**, *28*, 1294–1296. [CrossRef]

28. Millard, P.; Delépine, B.; Guionnet, M.; Heuillet, M.; Bellvert, F.; Létisse, F. IsoCor: Isotope Correction for High-Resolution MS Labeling Experiments. *Bioinformatics* **2019**, *35*, 4484–4487. [CrossRef]

29. Neidhardt, F.C.; Curtiss, R. *Escherichia Coli and Salmonella: Cellular and Molecular Biology*; ASM Press: Washington, DC, USA, 1996; ISBN 978-1-55581-084-9.

30. Sokol, S.; Millard, P.; Portais, J.-C. Influx_s: Increasing Numerical Stability and Precision for Metabolic Flux Analysis in Isotope Labelling Experiments. *Bioinforma. Oxf. Engl.* **2012**, *28*, 687–693. [CrossRef]

31. Nicolas, C.; Kiefer, P.; Letisse, F.; Krömer, J.; Massou, S.; Soucaille, P.; Wittmann, C.; Lindley, N.D.; Portais, J.-C. Response of the Central Metabolism of *Escherichia coli* to Modified Expression of the Gene Encoding the Glucose-6-Phosphate Dehydrogenase. *FEBS Lett.* **2007**, *581*, 3771–3776. [CrossRef]

32. Millard, P.; Massou, S.; Wittmann, C.; Portais, J.-C.; Létisse, F. Sampling of Intracellular Metabolites for Stationary and Non-Stationary 13C Metabolic Flux Analysis in *Escherichia coli*. *Anal. Biochem.* **2014**, *465*, 38–49. [CrossRef]

33. Haverkorn van Rijsewijk, B.R.B.; Nanchen, A.; Nallet, S.; Kleijn, R.J.; Sauer, U. Large-Scale 13C-Flux Analysis Reveals Distinct Transcriptional Control of Respiratory and Fermentative Metabolism in *Escherichia coli*. *Mol. Syst. Biol.* **2011**, *7*, 477. [CrossRef] [PubMed]

34. Zhao, J.; Baba, T.; Mori, H.; Shimizu, K. Effect of Zwf Gene Knockout on the Metabolism of *Escherichia coli* Grown on Glucose or Acetate. *Metab. Eng.* **2004**, *6*, 164–174. [CrossRef]

35. Hua, Q.; Yang, C.; Baba, T.; Mori, H.; Shimizu, K. Responses of the Central Metabolism in *Escherichia coli* to Phosphoglucose Isomerase and Glucose-6-Phosphate Dehydrogenase Knockouts. *J. Bacteriol.* **2003**, *185*, 7053–7067. [CrossRef]

36. Giacomoni, F.; Le Corguillé, G.; Monsoor, M.; Landi, M.; Pericard, P.; Pétéra, M.; Duperier, C.; Tremblay-Franco, M.; Martin, J.-F.; Jacob, D.; et al. Workflow4Metabolomics: A Collaborative Research Infrastructure for Computational Metabolomics. *Bioinformatics* **2015**, *31*, 1493–1495. [CrossRef]

37. Sauer, U.; Canonaco, F.; Heri, S.; Perrenoud, A.; Fischer, E. The Soluble and Membrane-Bound Transhydrogenases UdhA and PntAB Have Divergent Functions in NADPH Metabolism of *Escherichia coli*. *J. Biol. Chem.* **2004**, *279*, 6613–6619. [CrossRef]

38. Long, C.P.; Antoniewicz, M.R. High-Resolution 13C Metabolic Flux Analysis. *Nat. Protoc.* **2019**, *14*, 2856–2877. [CrossRef]

39. Lange, A.; Becker, J.; Schulze, D.; Cahoreau, E.; Portais, J.-C.; Haefner, S.; Schröder, H.; Krawczyk, J.; Zelder, O.; Wittmann, C. Bio-Based Succinate from Sucrose: High-Resolution 13C Metabolic Flux Analysis and Metabolic Engineering of the Rumen Bacterium Basfia Succiniciproducens. *Metab. Eng.* **2017**, *44*, 198–212. [CrossRef]

40. Leighty, R.W.; Antoniewicz, M.R. COMPLETE-MFA: Complementary Parallel Labeling Experiments Technique for Metabolic Flux Analysis. *Metab. Eng.* **2013**, *20*, 49–55. [CrossRef] [PubMed]

41. Crown, S.B.; Long, C.P.; Antoniewicz, M.R. Integrated 13C-Metabolic Flux Analysis of 14 Parallel Labeling Experiments in *Escherichia coli*. *Metab. Eng.* **2015**, *28*, 151–158. [CrossRef] [PubMed]
42. Dai, D.; Reusch, R.N. Poly-3-Hydroxybutyrate Synthase from the Periplasm of *Escherichia coli*. *Biochem. Biophys. Res. Commun.* **2008**, *374*, 485–489. [CrossRef]

43. Juncker, A.S.; Willenbrock, H.; Von Heijne, G.; Brunak, S.; Nielsen, H.; Krogh, A. Prediction of Lipoprotein Signal Peptides in Gram-Negative Bacteria. *Protein Sci. Publ. Protein Soc.* **2003**, *12*, 1652–1662. [CrossRef] [PubMed]

44. Heuillet, M.; Bellvert, F.; Cahoreau, E.; Letisse, F.; Millard, P.; Portais, J.-C. Methodology for the Validation of Isotopic Analyses by Mass Spectrometry in Stable-Isotope Labeling Experiments. *Anal. Chem.* **2018**, *90*, 1852–1860. [CrossRef] [PubMed]

45. Peiro, C.; Millard, P.; de Simone, A.; Cahoreau, F.; Peyriga, L.; Enjalbert, B.; Heux, S. Chemical and Metabolic Controls on Dihydroxyacetone Metabolism Lead to Suboptimal Growth of *Escherichia coli*. *Appl. Environ. Microbiol.* **2019**, *85*. [CrossRef] [PubMed]

46. Chazalviel, M.; Frainay, C.; Poupin, N.; Vinson, F.; Merlet, B.; Gloaguen, Y.; Cottret, L.; Jourdan, F. MetExploreViz: Web Component for Interactive Metabolic Network Visualization. *Bioinformatics* **2018**, *34*, 312–313. [CrossRef] [PubMed]