Next Article in Journal
Metabolic View on Human Healthspan: A Lipidome-Wide Association Study
Next Article in Special Issue
A Novel Solid Phase Extraction Sample Preparation Method for Lipidomic Analysis of Human Plasma Using Liquid Chromatography/Mass Spectrometry
Previous Article in Journal
Data Processing Optimization in Untargeted Metabolomics of Urine Using Voigt Lineshape Model Non-Linear Regression Analysis
Previous Article in Special Issue
Application of Lipid Class Ratios for Sample Stability Monitoring—Evaluation of Murine Tissue Homogenates and SDS as a Stabilizer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 11
Issue 5
10.3390/metabo11050286
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control of n-Butanol Induced Lipidome Adaptations in E. coli

by
Aike Jeucken
1,
Miaomiao Zhou
2,
Marc M. S. M. Wösten
3 and
Jos F. Brouwers
2,4,*
1
Membrane Enzymology, Groningen Biomolecular and Biotechnology Institute (GBB), University of Groningen, 9747 AG Groningen, The Netherlands
2
Research Group Analysis Techniques in the Life Sciences, School of Life Sciences and Environmental Technology ATGM, Avans University of Applied Sciences, 4818 AJ Breda, The Netherlands
3
Infection Biology, Department of Biomolecular Health Sciences, Utrecht University, 3584 CL Utrecht, The Netherlands
4
Center for Molecular Medicine, University Medical Center Utrecht, 3584 CX Utrecht, The Netherlands
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(5), 286; https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11050286
Submission received: 15 March 2021 / Revised: 21 April 2021 / Accepted: 28 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Mass Spectrometry-Based Lipidomics)
Browse Figures
Versions Notes

Abstract

:
The versatile compound n-butanol is one of the most promising biofuels for use in existing internal combustion engines, contributing to a smooth transition towards a clean energy society. Furthermore, n-butanol is a valuable resource to produce more complex molecules such as bioplastics. Microbial production of n-butanol from waste materials is hampered by the biotoxicity of n-butanol as it interferes with the proper functioning of lipid membranes. In this study we perform a large-scale investigation of the complete lipid-related enzyme machinery and its response to exposure to a sublethal concentration of n-butanol. We profiled, in triplicate, the growth characteristics and phospholipidomes of 116 different genetic constructs of E. coli, both in the presence and absence of 0.5% n-butanol (v/v). This led to the identification of 230 lipid species and subsequently to the reconstruction of the network of metabolites, enzymes and lipid properties driving the homeostasis of the E. coli lipidome. We were able to identify key lipids and biochemical pathways leading to altered n-butanol tolerance. The data led to new conceptual insights into the bacterial lipid metabolism which are discussed.

1. Introduction

The biosynthesis of n-butanol has received much attention over the past decades because of its potential in many different processes, most particular as a biofuel [1,2,3]. Chemical advantages of n-butanol over other biofuels are the high energy content and low hygroscopic- and corrosive properties. The lower flashpoint compared to diesel and ethanol makes n-butanol safer to work with. Additionally, with regard to human- and environmental health, n-butanol represents strongly reduced hazards and this contributes further to its interest [4,5,6]. Except as a biofuel, n-butanol has many other applications in chemistry such as in the production of paints, resins and plastics [7,8].
Bacterial strains of Clostridium have long been the model organisms for n-butanol production via acetone-butanol-ethanol fermentation (reviewed in [9]). Unfortunately, yields have remained rather low despite metabolic re-engineering and applying advanced co-cultures of micro-organisms [9,10,11]. This has spiked interest in engineering other model organisms such as Saccharomyces and Escherichia coli to produce n-butanol because of their excellent genetic accessibility [12,13,14,15].
Recovery of n-butanol from the production environment remains problematic [16]. Distillation is not energy efficient and hence not economically feasible. Adsorption of dissolved n-butanol with various sorbents has proven to be a very efficient process but unfortunately, desorption is often an issue, as is recycling of the sorbent [17]. Another technique that is often used to recover n-butanol, is gas stripping. Although this requires only moderate investments in machinery or materials, the excessive amounts of foam that are produced, have prevented a breakthrough for this method. Finally, also the performance of filters has been poor due to viscosity of n-butanol as an extractant [18]. Nevertheless, the efficiency of all methods increases with increasing n-butanol concentrations and achieving higher tolerance to n-butanol, remains key.
The intrinsic chemical properties of n-butanol present challenges to efficient bioproduction. n-butanol is a solvent in which hydrophobic compounds such as membrane lipids, readily dissolve. It is, therefore, self-evident that optimizing lipid metabolism can result in increased n-butanol tolerance. The high maximum solubility of n-butanol at room temperature (73 g/L or 7.7% v/v) makes it an ambitious goal to produce n-butanol at levels that would result in spontaneous phase separation. However, butanol-tolerant bacteria that can grow at this concentration have been isolated from the environment, demonstrating that this goal is realistic [19].
Toxic effects of organic solvents are due to detrimental effects on phospholipid membrane structure. Toxicity and hydrophobicity of primary alcohols are directly related [20,21]. Recent work has unraveled the different mechanisms by which membrane disturbance occurs [22]. These authors demonstrate how increasing n-butanol concentrations induce increased lipid disordering, hence increase membrane fluidity and permeability. However, the work is mainly on artificial membranes of simplified lipid composition, ignoring any bacterial capacity to adapt to a changing environment.
Increased short-chain alcohol tolerance can be achieved in E. coli by engineering active transport of these product out of the cell [23]. Intrinsically, genome wide generic stress responses are involved in improved n-butanol tolerance which include lipid membrane functions [24,25,26]. Involvement of lipid metabolism in response to n-butanol induced stress has been demonstrated in Clostridia [27], but for E. coli similar data is missing.
Here we explore this bacterial potential to adapt the lipidome in response to n-butanol exposure. We found a 50% inhibition of growth rate at 0.5% n-butanol (v/v) and chose to explore changes in lipid metabolism at this concentration. To this end, we overexpress all known lipid-related enzymes individually or, when possible, knock-out these genes. We thus aim to identify the pivotal genes of control and to identify marker lipids. Simultaneous recording of growth rates and biomass formation as characteristics of bacterial health link lipidome to n-butanol tolerance. This work thereby contributes to a rationale for optimization strategies in bioproduction of n-butanol and other hydrophobic compounds.

2. Results

In E. coli, the phospholipidome is controlled by a relatively modest (known) set of 68 proteins of which only 20 are essential (as defined by the Keio collection [28]). These proteins are organized in several pathways related to fatty acid synthesis or -breakdown, glycerol backbone metabolism or are related to phospholipid modifications or -synthesis. These pathways and corresponding enzymes are visualized in Figure 1, highlighting the essential enzymes.
Cultures of strains overexpressing each of these genes individually were analyzed for growth characteristics (maximum growth rate and OD600 at steady state), as were the individual knock-outs of these genes. This allowed us to investigate lipidomes under maximum and minimum expression levels of lipid related genes.

2.1. Simplification of the Lipid Extraction Allows for More Detailed Lipidome Analysis

The simplified, one step extraction of lipids from the bacterial membranes resulted in a more complex lipidome than generally reported with liquid–liquid extraction [30,31] (Figure 2). This is in part due to the preservation of very polar lipids like CDP-DAG that normally partition towards the polar liquid phase in liquid–liquid extractions and are hence lost from the lipid extract. Simultaneously, the wide range of conditions analyzed here included conditions in which minor lipids became sufficiently abundant to survive noise filtering of raw data and were included. Forced integration of corresponding mz/rt features in samples where these compounds were initially not detected, assured reliable data for these minor compounds (e.g., aPE, aPG, DLCL) in all samples. In total, 230 distinct lipid species, distributed over 11 phospholipid classes were included in further analysis. Each lipid-related gene was individually overexpressed, and three separate cultures of the construct were extracted and analyzed. Similarly, knock-outs of each non-essential lipid-related gene were cultured in triplicate and extracted. As controls for the overexpressors and knock-out constructs, we used the E. coli wild-type strain (BW25113) with or without an empty overexpression vector, respectively. These controls were included on every 96-well plate used for culturing bacteria and were, therefore, conducted more than three times. In the very few cases where bacterial growth was insufficient to reach an optical density (OD600 nm) of 0.2, samples were excluded from further analysis. This resulted in the acquisition of 730 lipidomes (Table S1).

2.2. Random Forest Analysis Identifies Key Lipids and Lipid Genes

Breiman’s random forest analysis was applied to identify the key features [32]. Here, we applied random forest analysis to find lipid species that could classify the genetic status (gene ID and KO or OV) and n-butanol exposure (Yes/No) of cultured bacteria. Random Forest Analysis produced a top-10 panel of lipid species that collectively were able to perform accurate classification of all samples (Table 1).
Within the dataset, these lipids contribute only 21.2% of the total lipid related MS signal of wild-type bacteria. The top-4 species in Table 1 are within the 12 most abundant species in the total lipidome (Table S1). Notably, the other six species in Table 1 (PG 31:0c1 through aPE 50:1) are quantitatively minor compounds and consist nearly exclusively of phospholipid classes that are thus far mostly ignored in E. coli lipidomic papers, and their identification as marker lipid species was unexpected. The phospholipid class phosphatidylbutanol (PBut) is not a common class and its presence under standard growth conditions has not been reported. Formation of this lipid class is dependent on the presence of n-butanol in the medium. n-butanol is not a common constituent of growth medium nor is it a metabolite in the common metabolism of E. coli. However, we have previously shown that the cardiolipin synthase ClsB functions as a promiscuous phospholipase D, capable of exchanging phospholipid headgroups with primary alcohols [33]. In this light, it is not surprising to find highest levels of the PBut marker lipid in the ClsB overexpressing strain (Table 1).
A less obvious but noteworthy characteristic of the lipid markers in Table 1 is the fact that every species has a distinct fatty acyl composition. It was reported previously that lipid species from distinct lipid classes but with the same acyl composition, displayed similar responses to genetic- and environmental challenges [29]. The presence of only distinct acyl compositions in the set of marker lipid species may thus reflect the random forest analysis algorithm eliminating information redundancy. Although it is apparent that the genes involved in achieving the extremes (Table 1) is not a random cross-section of all lipid related genes, care must be taken when interpreting this list. Particularly when looking at low lipid abundances, the gene expression modification leading to the lowest values, did often not differ significantly from the second or third lowest values. For instance, the knock-out of Cfa (the enzyme synthesizing the cyclopropane moiety) led to very low abundances of all marker lipid species containing a cyclopropane ring (indicated by the subscript ‘c1′ in the lipid name). To create a visual impression of the genes with highest impact on the marker lipids, we constructed a tag cloud from the 100 gene tags of the top- and bottom five genes in the 10 marker lipids (Figure 3). From the tag cloud it became evident that Cfa indeed has a large impact on the abundance of the marker lipid species. Other notable enzymes include FabH (involved in the initiation of fatty acid synthesis), ClsA (the dominant cardiolipin synthase), GlpD (glycerol backbone synthesis) and Aas (headgroup acylation).

2.3. The Lipidome Is Organized in Clusters of Closely Connected Lipid Species

Next, we investigated the underlying concepts that defined these ten lipids as informative markers. We hypothesized that these markers were representatives of clusters of lipid species defining functional domains in lipid membranes, or at least of clusters of lipid species whose abundances in the membrane are closely related. Therefore, we constructed a heatmap of lipid–lipid correlations and organized lipids by hierarchical clustering (Figure 4). Indeed, clusters of lipids with strongly correlating abundances (dark blue or -red) were visible, correlating lipids within one cluster (on the diagonal) or between different clusters (off-diagonal). Notably, the marker lipids from Table 1 were scattered over these clusters. Within clusters, only a single marker lipid was present. As an exception, the second most abundant lipid species in E. coli, PE 32:1 and part of the identified marker lipids, appeared not to be part of any lipid cluster (bottom right corner). The horizontal color bars on top of the heatmap, visualize lipid class, degree of unsaturation and presence of cyclopropane rings of the individual lipids within the clusters. Although we detected clusters of lipids from the same lipid class (particularly for the headgroup acylated lipid classes aPE and aPG), all major lipid classes were divided over several clusters. Similarly, lipid unsaturation or the presence of cyclopropane moieties were not driving forces of the lipid clusters formation. Hence, the formation of lipid clusters in the correlation heatmap was a result of a complex interplay of multiple processes.

2.4. Validation of Metabolic Pivot Points

We observed large differences between minimum and maximum contributions of the marker lipids to the lipidome (Table 1). Even the abundance of major compounds can be manipulated to fall below 1% without loss of viability. PE 32:1 the observed range in abundancies was particularly wide (0.75–24.82%). Notably, these extremes are both obtained by alteration of the expression level of the cfa gene. Although PE 32:1 itself does not contain a cyclopropane ring, our data point towards an important role for PE 32:1 as a substrate for the Cfa enzyme. Lack of this enzyme leads to extensive accumulation of this lipid species, whereas overexpression resulted in nearly complete depletion. Concomitantly, the abundance of PE 33:0c1 (the Cfa product of PE 32:1) fell under these conditions from 25.9 ± 0.5% (Cfa overexpression) to 1.5 ± 0.2% (Cfa knock-out; Table S1). The lipid species PE 34:1, which also contains one single unsaturation that may be targeted by Cfa, was much less affected by Cfa overexpression than PE 32:1. This further illustrates the complex control of the phospholipidome.
The high impact FabH has on the marker lipids (Figure 3), was investigated in more detail. FabH is involved in the initiation of fatty acid synthesis by coupling acetyl-CoA to malonyl-ACP, generating an oxoacyl-ACP and recycling free coenzyme-A. After a series of reductions by other Fab-proteins, the newly synthesized C4-fatty acyl-ACP is elongated using a new molecule of acetyl-CoA and this process typically continues until a chainlength of C14 or longer is achieved. Increasing expression of FabH (KO < (WT or EV) < OV) resulted in a decrease of average fatty acyl chainlength from 17.5 to 16.0 carbon atoms (Figure 5, top row). This means that bilayer thickness may be expected to decrease by 10% upon FabH overexpression. The presence of n-butanol in the culture medium did not affect this trend, although it should be mentioned that in the knock-out, the decrease in average chainlength from 17.5 (CTR) to 17.3 (+butanol) was statistically significant (p < 0.01). It is, however, unlikely to influence biophysical properties of the membrane.
The degree of unsaturation of lipid acyl chains was affected in a similar way as the chainlength. The average number of unsaturations per fatty acyl chain dropped from 0.41 in the KO to 0.27 in the combined WT and EV samples (p < 5 × 10−5), and again significantly further to 0.13 in the OV samples (p < 2 × 10−16). Hence, the FabH overexpressor has only 31% of the unsaturations the knock-out has and this may be expected to reduce membrane fluidity. Cyclopropane moieties are synthesized from unsaturated fatty acids and abundance of these characteristics are overall negatively correlated. However, the amount of cyclopropane containing fatty acids was not nearly affected to the extend observed for chainlength and unsaturation, but culturing in the presence of n-butanol consistently reduced the presence of cyclopropane containing fatty acids (Figure 5, bottom row).

2.5. Linking Lipid Related Genes to n-Butanol Tolerance

The toxic properties of n-butanol lead to an approximate 50% decrease in growth rate for WT E. coli. We investigated the effect of gene expression alterations on both the specific growth rate as well as on the OD600 achieved in steady state, the latter being a measure for biomass. Bacterial strains with higher metabolic efficiency are able to use more nutrients for cell growth and division but this not necessarily corresponds to higher growth rates.
When grown in the presence of n-butanol, 42 out of the 116 genetic constructs had maximum specific growth rates distinct from controls (Figure 6 and Supplementary Table S3 for all growth data). Notably, in ten of these the growth rate was higher, suggesting that the extreme genetic conditions in our OV/KO strains cannot functionally be achieved in the wild-type strain by regulation of gene expression, translation and/or enzyme activity. Looking at these extremes, it stands out that alterations in glycerol metabolism have very pronounced effects. Highest growth rates were achieved by overexpressing the (glycerol) kinase GlpK or knocking-out the Glp repressor GlpR. Lowest growth rates on the other hand, were achieved by either knocking-out or overexpressing GlpD, the only aerobic variant of the four glycerol-3-phosphate dehydrogenases [34]. Fatty acid characteristics for gene constructs involving GlpD and GlpK are summarized in Supplementary Figure S2.
The optical density (i.e., biomass) achieved at steady state was more variable between growth conditions than the growth rate (Figure 6, bottom panel). Most conditions did have a biomass production that was different from control conditions. Interestingly, GlpK overexpression stimulated growth rate as well as biomass production, whereas GlpD knockout on the other hand, impaired both growth rate and biomass production. Both enzymes are involved in both lipid metabolism and central carbon metabolism but from these data it is impossible to determine which of these metabolic pathways is responsible for the altered growth characteristics.

3. Discussion

As a powerful machine learning algorithm, if applied to properly generated data sets, Breiman’s random forest analysis can generate robust results that are applicable to previously unseen data. It is essential that the data sets are created in a steady fashion with reliable normalization standards [35]. The current setup, where bacteria were grown in multiwell plates and lipid extractions was simplified to a one-step procedure, enabled us to generate such a large training set. Since a lipidome can be analyzed within five minutes (Figure 2), the collective samples could be processed within only three days of LC-MS analysis. The resulting set of ten marker lipids not only could be used to identify the genetic status and n-butanol exposure of all samples but also helped in identifying clusters of functionally related lipids. The exact nature of how lipids within these lipid clusters are related is not yet fully understood but clearly this exceeds just lipid class or acyl composition. Particularly abundant or sparse lipid clusters can be correlated to specific genetic conditions, and this can provide further leads for investigation.
It should be noted that distinct cultures of the same genetic construct of E. coli, produced highly consistent lipidome compositions. Therefore, statistical significance was easily obtained between different strains or conditions, even when changes were too small to be expected to alter biophysical properties of the membranes. Put differently: statistical significance should not be confused with biological relevance and we have tried to highlight relevant processes rather than to report as many differences as possible. Readers are encouraged to query the dataset for changes in processes of their own interests.
Because of the experimental setup in which half of the samples were cultured in the presence of n-butanol, it can be expected that many of the phenomena we observed were related to coping with induced membrane stress. We expect this to be reflected in the composition of the ten marker lipids. Indeed, most of the marker lipids will have high impact on lipid packing in the membrane. This is particularly the case for the headgroup acylated lipid species aPE and aPG, where the additional acyl chain will fold back into the hydrophobic region of the bilayer, thus profoundly changing the shape of the phospholipid and increasing the inter- and intra-molecular hydrophobic interactions. In fact, membrane splaying as a driving force for membrane collapse (i.e., disturbing lipid–lipid interactions by accumulation of n-butanol at the phospholipid headgroups [22]), may be effectively counteracted. In this respect, it is important to acknowledge that these headgroup acylates have mostly been ignored in E. coli lipidomic papers due to low abundance and/or unfamiliarity of scientists with this lipid class. It can be very worthwhile to reassess results of these papers given these new insights and to consider re-analyzing crucial experiments. Our rapid and labor-extensive lipid extraction and short analysis time will aid in making such goals achievable.
Although the HILIC-based lipidomic workflow enabled us to note and identify unusual lipid classes, we were only able to identify lipids at the species level, i.e., data lack information about acyl composition and -position. When a more detailed level of information is essential, modern very fast MS/MS instruments may allow profiling at the fatty acyl (position) level. Alternatively, researchers may opt for reverse-phase approaches where separation is primarily dictated by the acyl composition. However, better species resolution comes at the cost of longer analysis times [36].
Cyclopropane moieties are synthesized by the enzyme Cfa from an unsaturation in an acyl chain. The measure of unsaturated acyl chains in a biomembrane is typically considered to be the key property when assessing fluidity of membranes. Since there exists no ‘saturase’ that will remove a double bond from a fatty acyl, converting a double bond to a rigid cyclopropane moiety will be an effective way to counteract the increased membrane fluidity induced by n-butanol partitioning [37]. Indeed, the role of cyclopropane fatty acids in organic solvent tolerance was demonstrate in Pseudomonas putida [38] and also in other bacteria including E. coli, a role in coping with environmental stress has firmly been established [39,40,41]. In a recent paper, we have shown that, even compared to any other lipid-related gene, both deletion and overexpression of Cfa has very profound effects on the overall lipidome of E. coli [29]. It is, therefore, not surprising that Cfa was among the most relevant genes in Figure 3. However, the cfa knockout did not have a reduced growth rate. Additionally, biomass increase during growth was only marginally lower than WT, indicating that Cfa has limited importance in dealing with butanol stress.
Growth data pointed towards a crucial role for glycerol metabolism rather than fatty acid- or phospholipid metabolism in n-butanol tolerance. Overexpression of GlpK led to the highest growth rate as well as a large increase in biomass. GlpK is a known point of (allosteric) regulation of carbon metabolism [42,43] and has been linked to cell stress response [44]. Furthermore, substrate specificity is not particularly high, which makes GlpK active in several metabolic pathways [45]. The effect of GlpK overexpression may therefore reflect its involvement in responses unrelated to lipid metabolism. The lack of a clearly aberrant lipidome in this construct supports this idea.
Manipulation of GlpD levels on the other hand, had strong negative effects on both growth rate and biomass. GlpD is a key enzyme in glycerol metabolism as it was shown that growth on glycerol strictly depends on the presence of GlpD [46]. Hence, the availability of physiological levels of the phospholipid precursor glycerol-3P appears to be an important prerequisite for a healthy lipid homeostasis.
Taken together, this large-scale MS-based lipidomic approach has revealed the presence of clusters of lipids whose abundance is related. It was also demonstrated that the 116 different genetic constructs and growth conditions can be classified by ten marker lipids. The abundance of these marker lipids was shown to be highly variable. Furthermore, these data show a strong link between lipid related genes and n-butanol tolerance. This opens new perspectives for targeted membrane engineering in achieving higher membrane resistance to organic solvents.

4. Materials and Methods

4.1. Chemicals

Isopropyl β-D-1-thiogalactopyranoside (IPTG) was obtained from Melford (Suffolk, UK). NaCl, n-butanol, and formic acid and were purchased from Merck Chemicals (Darmstad, Germany). Chloramphenicol was obtained from Boehringer Mannheim (Mannheim, Germany). Chemicals were of the highest purity available. Trypton, agar, and kanamycin were purchased from Sigma (St. Louis, MO, USA). Yeast extract was obtained from MP Biomedicals (Strasbourg, France). Methanol, acetonitrile (ACN), acetone and ammonium formate were purchased from BioSolve (Valkenswaard, The Netherlands), chloroform was obtained from Roth (Karlsruhe, Germany) and were all HPLC/MS grade.

4.2. Bacterial Strains, Growth Conditions and Plasmids

The bacterial strains and plasmids used in the study are listed in Supplementary Materials and Methods (Supplementary Table S2). E. coli strain BW25113 and its derivates are part of the Keio collection [28], obtained from NBRP (NIG, Mishima, Japan). E. coli strains were routinously grown at 37 °C in Luria Bertani (LB) broth or on LB agar plates. Plasmids overexpressing lipid genes present in E. coli strain AG1(ME5305) (part of the ASKA(-) collection [47] obtained from NBRP (NIG, Mishima, Japan) were transferred to E. coli strain BW25113. When appropriate, growth media were supplemented with chloramphenicol (34 μg/mL) or kanamycin (50 μg/mL).
Protein expression was induced with a final concentration of 10 μM IPTG for the ASKA(-) clones since the ORF on the pCA24N plasmid are under control of IPTG-inducible promoter, PT5-lac. For phospholipid analysis, cultures (150 μL) were grown in 96 wells plates, in a Versamax microtiter plate reader (Molecular Devices, Sunnyvale, CA, USA), in the absence or presence of 0.5% (v/v) n-butanol. To this end, pre-cultures that had research stationary phase in 9 h growth, were diluted 1:100 and growth was followed. The plate was covered using a clear film to reduce evaporation. Absorbance at 600 nm was measured every 5 min. These data were used to calculate maximum specific growth rates as defined by the slope of a linear fit to log transformed OD values.

4.3. Lipid Extraction

Lipids were extracted as described before [29]. In brief, bacterial cultures were transferred to glass coated plates, centrifuged (1800× g, 20 min, 4 °C) and the medium was completely removed by aspiration. Pellets were resuspended without further washing in 150 µL chloroform/methanol (1:1 v/v), extracted for 1 h at 4 °C, followed by centrifugation (1800× g, 20 min, 4 °C). The extracts were transferred to a new glass coated 96 well plate that was covered by a thin sheet of aluminum foil and placed in the autosampler. WT strain, empty vector, and cell free incubations were included with every plate. Cell free incubations did not produce meaningful lipid signals and these samples were excluded from further analysis.

4.4. Liquid Chromatography Mass Spectrometry of Lipids

LC-MS of phospholipids was performed as described before [29,48]. In brief, 10 µL of the supernatant was injected onto a hydrophilic interaction liquid chromatography (HILIC) column (2.6 μm HILIC 100 Å, 50 × 4.6 mm, Phenomenex, Torrance, CA, USA), and eluted with a gradient from ACN/Acetone (9:1, v/v) to ACN/H2O (7:3, v/v) with 10 mM ammonium formate, and both with 0.1% formic acid at a flow rate of 1 mL/min. Gradient development was as follows (time in min, %B): (0, 0), (1, 50), (3, 50), (3.1, 100), (4, 100). Regeneration time between subsequent runs was 55 s, which was identical to our autosampler injection time and therefore not listed separately. The column outlet of the LC was connected to a heated electrospray ionization (HESI) source of a LTQ XL mass spectrometer (ThermoFisher Scientific, Waltham, MA, USA). Full scan spectra were collected in negative ionization mode in the range from 350–1750 amu at a scan speed of 3 scans/s. Source- and capillary temperatures were set to 450 and 400 °C, respectively, and the ionization voltage to −2.5 kV. Calibration curves with constructed using authentic standards of PE, PG and CL to ensure linear responses for the lipid concentrations used in the experiments. Due to the lack of standard for all lipid classes (e.g., phosphatidylbutanol, head group acylated PG), we have refrained from conversion of lipid data to absolute quantities.

4.5. Data Analysis

LC-MS data were converted to mz(X)ML format using MSconvert and analyzed using XCMS version 1.52.0 using R version 3.4.4 (15 March 2018) [49,50,51]. Phospholipid (sub-)classes were identified based on retention time and mz/rt features were matched against an in-silico generated lipid database. Lipid de-isotoping was conducted using a natural abundance of 13C of 1.1%. Lipids were reported on the lipid species level rather than lipid class level mass based on the well documented absence of ether lipids in E. coli [52]. A normalization was performed such that the sum of all species is 100% in every sample. An overall cross-comparison on lipidomes was performed by using the CyberT tool package [53]. BayesAnova tests were applied on all possible lipidome combinations (parameters: sliding window size 101, Bayes confidence estimate value 5, and PPDE post processing treatment). Comparison resulted with an adjusted p-value lower than 0.05 were considered significant. Breiman’s random forest analysis (500 trees and 100 iteration) was applied to the normalized data set by using Weka 3.9 [54]. For significant feature selection, sequential backward selection was used with the 34OOB error minimization option. During the analysis, 10-fold cross-validation was applied. The feature importance evaluations were sorted by eliminating the least important features.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo11050286/s1, Table S1: Normalized lipidomic analysis of bacterial constructs. Table S2: bacterial strains and plasmids. Figure S1: Lipid class structures. Figure S2: lipid characteristics for GlpD and GlpK.

Author Contributions

Conceptualization, A.J. and J.F.B.; methodology, A.J., M.M.S.M.W., J.F.B. and M.Z.; formal analysis, A.J., J.F.B. and M.Z.; investigation, A.J., M.M.S.M.W. and J.F.B.; data curation, A.J., J.F.B. and M.Z.; writing—original draft preparation, J.F.B.; writing—review and editing, A.J., J.F.B., M.M.S.M.W. and M.Z.; visualization, A.J., J.F.B. and M.Z.; supervision, J.F.B.; funding acquisition, J.F.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grant BB/2012-10047 of the Bio-Based Ecologically Balanced Sustainable Industrial Chemistry (BE-BASIC) project to J.B.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Lipidomic peak lists and growth data of individual bacterial cultures have been made available as supplementary data.

Acknowledgments

We thank Sacha van Hijum for his contribution in the preparation of the random forest analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, C.; Zhang, Y.; Li, Y. Metabolic Engineering for the Production of Butanol, a Potential Advanced Biofuel, from Renewable Resources. Biochem. Soc. Trans. 2020, 48, 2283–2293. [Google Scholar] [CrossRef] [PubMed]
  2. Jin, C.; Yao, M.; Liu, H.; Lee, C.F.; Ji, J. Progress in the Production and Application of N-Butanol as a Biofuel. Renew. Sustain. Energy Rev. 2011, 15, 4080–4106. [Google Scholar] [CrossRef]
  3. Nair, P.; Meenakshi, H.N. Review on the Synthesis, Performance and Trends of Butanol: A Cleaner Fuel Additive for Gasoline. Int. J. Ambient Energy 2021, 1–17. [Google Scholar] [CrossRef]
  4. Mehlman, M.A. Dangerous Properties of Petroleum-Refining Products: Carcinogenicity of Motor Fuels (Gasoline). Teratog. Carcinog. Mutagen. 1990, 10, 399–408. [Google Scholar] [CrossRef] [PubMed]
  5. Caprino, L.; Togna, G.I. Potential Health Effects of Gasoline and Its Constituents: A Review of Current Literature (1990–1997) on Toxicological Data. Environ. Health Perspect. 1998, 106, 115–125. [Google Scholar] [CrossRef] [PubMed]
  6. Ekpenyong, C.E.; Asuquo, A.E. Recent Advances in Occupational and Environmental Health Hazards of Workers Exposed to Gasoline Compounds. Int. J. Occup. Med. Environ. Health 2017, 30, 1–26. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, H.; Kang, S.; Li, J.; Li, Z.; Chang, J.; Xu, Y.; Max Lu, G.; Sun, C. Near Zero-Waste Biofuel Production from Bioderived Polyhydroxybutyrate. Fuel 2021, 286, 119405. [Google Scholar] [CrossRef]
  8. Chen, G.-Q.; Patel, M.K. Plastics Derived from Biological Sources: Present and Future: A Technical and Environmental Review. Chem. Rev. 2012, 112, 2082–2099. [Google Scholar] [CrossRef]
  9. Li, S.; Huang, L.; Ke, C.; Pang, Z.; Liu, L. Pathway Dissection, Regulation, Engineering and Application: Lessons Learned from Biobutanol Production by Solventogenic Clostridia. Biotechnol. Biofuels 2020, 13, 39. [Google Scholar] [CrossRef] [Green Version]
  10. Du, Y.; Zou, W.; Zhang, K.; Ye, G.; Yang, J. Advances and Applications of Clostridium Co-Culture Systems in Biotechnology. Front. Microbiol. 2020, 11, 2824. [Google Scholar] [CrossRef] [PubMed]
  11. Volesky, B.; Szczesny, T. Bacterial Conversion of Pentose Sugars to Acetone and Butanol. Adv. Biochem. Eng. Biotechnol. 1983, 27, 101–118. [Google Scholar] [CrossRef]
  12. Azambuja, S.P.H.; Goldbeck, R. Butanol Production by Saccharomyces cerevisiae: Perspectives, Strategies and Challenges. World J. Microbiol. Biotechnol. 2020, 36, 48. [Google Scholar] [CrossRef] [PubMed]
  13. Atsumi, S.; Cann, A.F.; Connor, M.R.; Shen, C.R.; Smith, K.M.; Brynildsen, M.P.; Chou, K.J.Y.; Hanai, T.; Liao, J.C. Metabolic Engineering of Escherichia coli for 1-Butanol Production. Metab. Eng. 2008, 10, 305–311. [Google Scholar] [CrossRef] [PubMed]
  14. Shen, C.R.; Liao, J.C. Metabolic Engineering of Escherichia coli for 1-Butanol and 1-Propanol Production via the Keto-Acid Pathways. Metab. Eng. 2008, 10, 312–320. [Google Scholar] [CrossRef] [PubMed]
  15. Gulevich, A.Y.; Skorokhodova, A.Y.; Sukhozhenko, A.V.; Shakulov, R.S.; Debabov, V.G. Metabolic Engineering of Escherichia coli for 1-Butanol Biosynthesis through the Inverted Aerobic Fatty Acid β-Oxidation Pathway. Biotechnol. Lett. 2012, 34, 463–469. [Google Scholar] [CrossRef] [PubMed]
  16. Kujawska, A.; Kujawski, J.; Bryjak, M.; Kujawski, W. ABE Fermentation Products Recovery Methods-A Review. Renew. Sustain. Energy Rev. 2015, 48, 648–661. [Google Scholar] [CrossRef]
  17. Milestone, N.B.; Bibby, D.M. Concentration of Alcohols by Adsorption on Silicate. J. Chem. Technol. Biotechnol. 1981, 31, 732–736. [Google Scholar] [CrossRef]
  18. Groot, W.J.; Soedjak, H.S.; Donck, P.B.; van der Lans, R.G.J.M.; Luyben, K.C.A.M.; Timmer, J.M.K. Butanol Recovery from Fermentations by Liquid-Liquid Extraction and Membrane Solvent Extraction. Bioprocess Eng. 1990, 5, 203–216. [Google Scholar] [CrossRef]
  19. Kanno, M.; Katayama, T.; Tamaki, H.; Mitani, Y.; Meng, X.-Y.; Hori, T.; Narihiro, T.; Morita, N.; Hoshino, T.; Yumoto, I.; et al. Isolation of Butanol-and Isobutanol-Tolerant Bacteria and Physiological Characterization of Their Butanol Tolerance. Appl. Environ. Microbiol. 2013, 79, 6998–7005. [Google Scholar] [CrossRef] [Green Version]
  20. Heipieper, H.J.; Neumann, G.; Cornelissen, S.; Meinhardt, F. Solvent-Tolerant Bacteria for Biotransformations in Two-Phase Fermentation Systems. Appl. Microbiol. Biotechnol. 2007, 74, 961–973. [Google Scholar] [CrossRef]
  21. Sardessai, Y.; Bhosle, S. Organic Solvent Tolerant Bacteria in Mangrove Ecosystem. Curr. Sci. 2002, 82, 622–623. [Google Scholar]
  22. Guo, J.; Ho, J.C.S.; Chin, H.; Mark, A.E.; Zhou, C.; Kjelleberg, S.; Liedberg, B.; Parikh, A.N.; Cho, N.-J.; Hinks, J.; et al. Response of Microbial Membranes to Butanol: Interdigitation vs. Disorder. Phys. Chem. Chem. Phys. 2019, 21, 11903–11915. [Google Scholar] [CrossRef] [PubMed]
  23. Fisher, M.A.; Boyarskiy, S.; Yamada, M.R.; Kong, N.; Bauer, S.; Tullman-Ercek, D. Enhancing Tolerance to Short-Chain Alcohols by Engineering the Escherichia coli AcrB Efflux Pump to Secrete the Non-Native Substrate n-Butanol. ACS Synth. Biol. 2013, 3, 30–40. [Google Scholar] [CrossRef] [PubMed]
  24. Bui, L.M.; Lee, J.Y.; Geraldi, A.; Rahman, Z.; Lee, J.H.; Kim, S.C. Improved n-Butanol Tolerance in Escherichia coli by Controlling Membrane Related Functions. J. Biotechnol. 2015, 204, 33–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Rau, M.H.; Calero, P.; Lennen, R.M.; Long, K.S.; Nielsen, A.T. Genome-Wide Escherichia coli Stress Response and Improved Tolerance towards Industrially Relevant Chemicals. Microb. Cell Factories 2016, 15, 176. [Google Scholar] [CrossRef] [Green Version]
  26. Reyes, L.H.; Almario, M.P.; Winkler, J.; Orozco, M.M.; Kao, K.C. Visualizing Evolution in Real Time to Determine the Molecular Mechanisms of N-Butanol Tolerance in Escherichia coli. Metab. Eng. 2012, 14, 579–590. [Google Scholar] [CrossRef] [PubMed]
  27. Patakova, P.; Kolek, J.; Jureckova, K.; Branska, B.; Sedlar, K.; Vasylkivska, M.; Provaznik, I. Deeper below the Surface-Transcriptional Changes in Selected Genes of Clostridium beijerinckii in Response to Butanol Shock. MicrobiologyOpen 2021, 10, e1146. [Google Scholar] [CrossRef]
  28. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Jeucken, A.; Molenaar, M.R.; van de Lest, C.H.A.; Jansen, J.W.A.; Helms, J.B.; Brouwers, J.F. A Comprehensive Functional Characterization of Escherichia coli Lipid Genes. Cell Rep. 2019, 27, 1597–1606.e2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Smith, P.B.; Snyder, A.P.; Harden, C.S. Characterization of Bacterial Phospholipids by Electrospray Ionization Tandem Mass Spectrometry. Anal. Chem 1995, 67, 1824–1830. [Google Scholar] [CrossRef]
  31. Gidden, J.; Denson, J.; Liyanage, R.; Ivey, D.M.; Lay, J.O. Lipid Compositions in Escherichia coli and Bacillus subtilis During Growth as Determined by MALDI-TOF and TOF/TOF Mass Spectrometry. Int. J. Mass Spectrom. 2009, 283, 178–184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  33. Jeucken, A.; Helms, J.B.; Brouwers, J.F. Cardiolipin Synthases of Escherichia coli Have Phospholipid Class Specific Phospholipase D Activity Dependent on Endogenous and Foreign Phospholipids. Biochim. Biophys. Acta (BBA) Mol. Cell Biol. Lipids 2018, 1863, 1345–1353. [Google Scholar] [CrossRef] [PubMed]
  34. Yeh, J.I.; Chinte, U.; Du, S. Structure of Glycerol-3-Phosphate Dehydrogenase, an Essential Monotopic Membrane Enzyme Involved in Respiration and Metabolism. Proc. Natl. Acad. Sci. USA 2008, 105, 3280–3285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Karpievitch, Y.V.; Hill, E.G.; Leclerc, A.P.; Dabney, A.R.; Almeida, J.S. An Introspective Comparison of Random Forest-Based Classifiers for the Analysis of Cluster-Correlated Data by Way of RF++. PLoS ONE 2009, 4, e7087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Züllig, T.; Köfeler, H.C. High Resolution Mass Spectrometry in Lipidomics. Mass Spectrom. Rev. 2021, 40, 162–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Poger, D.; Mark, A.E. A Ring to Rule Them All: The Effect of Cyclopropane Fatty Acids on the Fluidity of Lipid Bilayers. J. Phys. Chem. B 2015, 119, 5487–5495. [Google Scholar] [CrossRef] [PubMed]
  38. Pini, C.-V.; Bernal, P.; Godoy, P.; Ramos, J.-L.; Segura, A. Cyclopropane Fatty Acids Are Involved in Organic Solvent Tolerance but Not in Acid Stress Resistance in Pseudomonas putida DOT-T1E. Microb. Biotechnol. 2009, 2, 253–261. [Google Scholar] [CrossRef] [Green Version]
  39. Chang, Y.Y.; Cronan, J.E. Membrane Cyclopropane Fatty Acid Content Is a Major Factor in Acid Resistance of Escherichia coli. Mol. Microbiol. 1999, 33, 249–259. [Google Scholar] [CrossRef]
  40. Kim, B.H.; Kim, S.; Kim, H.G.; Lee, J.; Lee, I.S.; Park, Y.K. The Formation of Cyclopropane Fatty Acids in Salmonella enterica Serovar Typhimurium. Microbiol 2005, 151, 209–218. [Google Scholar] [CrossRef] [Green Version]
  41. Grogan, D.W.; Cronan, J.E. Cyclopropane Ring Formation in Membrane Lipids of Bacteria. Microbiol. Mol. Biol. Rev. MMBR 1997, 61, 429–441. [Google Scholar] [CrossRef] [PubMed]
  42. Hurley, J.H.; Faber, H.R.; Worthylake, D.; Meadow, N.D.; Roseman, S.; Pettigrew, D.W.; Remington, S.J. Structure of the Regulatory Complex of Escherichia coli IIIGlc with Glycerol Kinase. Science 1993, 259, 673–677. [Google Scholar] [CrossRef] [PubMed]
  43. Zwaig, N.; Lin, E.C. Feedback Inhibition of Glycerol Kinase, a Catabolic Enzyme in Escherichia coli. Science 1966, 153, 755–757. [Google Scholar] [CrossRef] [PubMed]
  44. Khil, P.P.; Camerini-Otero, R.D. Over 1000 Genes Are Involved in the DNA Damage Response of Escherichia coli. Mol. Microbiol. 2002, 44, 89–105. [Google Scholar] [CrossRef] [PubMed]
  45. Thorner, J.W.; Paulus, H. Catalytic and Allosteric Properties of Glycerol Kinase from Escherichia coli. J. Biol. Chem. 1973, 248, 3922–3932. [Google Scholar] [CrossRef]
  46. Wegener, M.; Vogtmann, K.; Huber, M.; Laass, S.; Soppa, J. The GlpD Gene Is a Novel Reporter Gene for E. coli That Is Superior to Established Reporter Genes like LacZ and GusA. J. Microbiol. Methods 2016, 131, 181–187. [Google Scholar] [CrossRef]
  47. Kitagawa, M.; Ara, T.; Arifuzzaman, M.; Ioka-Nakamichi, T.; Inamoto, E.; Toyonaga, H.; Mori, H. Complete Set of ORF Clones of Escherichia coli ASKA Library (a Complete Set of E. coli K-12 ORF Archive): Unique Resources for Biological Research. DNA Res. 2005, 12, 291–299. [Google Scholar] [CrossRef] [Green Version]
  48. Brouwers, J.F.; Aalberts, M.; Jansen, J.W.A.; van Niel, G.; Wauben, M.H.; Stout, T.A.E.; Helms, J.B.; Stoorvogel, W. Distinct Lipid Compositions of Two Types of Human Prostasomes. Proteomics 2013, 13, 10–11. [Google Scholar] [CrossRef]
  49. Kessner, D.; Chambers, M.; Burke, R.; Agus, D.; Mallick, P. ProteoWizard: Open Source Software for Rapid Proteomics Tools Development. Bioinformatics 2008, 24, 2534–2536. [Google Scholar] [CrossRef]
  50. R Development Core Team. R: A Language and Environment for Statistical Computing; R Development Core Team: Vienna, Austria, 2011. [Google Scholar]
  51. Smith, C.A.; Want, E.J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing Mass Spectrometry Data for Metabolite Profiling Using Nonlinear Peak Alignment, Matching, and Identification. Anal. Chem. 2006, 78, 779–787. [Google Scholar] [CrossRef]
  52. Züllig, T.; Trötzmüller, M.; Köfeler, H.C. Lipidomics from Sample Preparation to Data Analysis: A Primer. Anal. Bioanal. Chem. 2020, 412, 2191–2209. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kayala, M.A.; Baldi, P. Cyber-T Web Server: Differential Analysis of High-Throughput Data. Nucleic Acids Res. 2012, 40, W553–W559. [Google Scholar] [CrossRef] [PubMed]
  54. Witten, I.H.; Frank, E.; Hall, M.A.; Pal, C.J. Data Mining, Fourth Edition: Practical Machine Learning Tools and Techniques, 4th ed.; Morgan Kaufmann Publishers Inc.: San Francisco, CA, USA, 2016; ISBN 978-0-12-804291-5. [Google Scholar]
Metabolites 11 00286 g001 550
Figure 1. Metabolic scheme of lipid related genes and their role in lipid homeostasis. Genes essential for cell survival and multiplication under generic conditions are shown in bold. Adapted from [29]. Lipid abbreviations: aPE: N-acyl-phosphatidylethanolamine; aPG: O-acyl phosphatidylglycerol (headgroup acylated PG); CDP-DAG: cytidine diphosphate diacylglycerol; CL: cardiolipin; DAG: diacylglycerol; DLCL: di-lyso cardiolipin; PA: phosphatidic acid; PE: phosphatidylethanolamine; PG: phosphatidylglycerol; PGP: phosphatidylglycerolphosphate; PS: phosphatidylserine.
Figure 1. Metabolic scheme of lipid related genes and their role in lipid homeostasis. Genes essential for cell survival and multiplication under generic conditions are shown in bold. Adapted from [29]. Lipid abbreviations: aPE: N-acyl-phosphatidylethanolamine; aPG: O-acyl phosphatidylglycerol (headgroup acylated PG); CDP-DAG: cytidine diphosphate diacylglycerol; CL: cardiolipin; DAG: diacylglycerol; DLCL: di-lyso cardiolipin; PA: phosphatidic acid; PE: phosphatidylethanolamine; PG: phosphatidylglycerol; PGP: phosphatidylglycerolphosphate; PS: phosphatidylserine.
Metabolites 11 00286 g001
Metabolites 11 00286 g002 550
Figure 2. Base peak chromatogram of HILIC-based LC-MS analysis of wild type E. coli lipid extract. Lipid species are depicted, and color coded according to the Z-score of culturing in the presence of 0.5% n-butanol (green corresponding to higher abundance at n-butanol exposure). Dot size corresponds to abundance in WT n-butanol-free cultures. The two most abundant species of each lipid class are labeled. All lipid ions corresponded to deprotonated molecules. Examples of structures for each class of lipid ions are given in Supplementary Figure S1.
Figure 2. Base peak chromatogram of HILIC-based LC-MS analysis of wild type E. coli lipid extract. Lipid species are depicted, and color coded according to the Z-score of culturing in the presence of 0.5% n-butanol (green corresponding to higher abundance at n-butanol exposure). Dot size corresponds to abundance in WT n-butanol-free cultures. The two most abundant species of each lipid class are labeled. All lipid ions corresponded to deprotonated molecules. Examples of structures for each class of lipid ions are given in Supplementary Figure S1.
Metabolites 11 00286 g002
Metabolites 11 00286 g003 550
Figure 3. Tag cloud of the altered gene expression levels leading to extremes in marker lipid species identified by random forest analysis.
Figure 3. Tag cloud of the altered gene expression levels leading to extremes in marker lipid species identified by random forest analysis.
Metabolites 11 00286 g003
Metabolites 11 00286 g004 550
Figure 4. Heatmap of correlations in the lipidome composition of all 700+ samples. The position of the ten marker lipids is indicated and the clusters are emphasized in the hierarchical trees (top and left). The colored bars (top) indicate the phospholipid class off each lipid species, the total number of unsaturations in the acyl chains and whether the acyl chains are linear of contain a cyclopropane moiety.
Figure 4. Heatmap of correlations in the lipidome composition of all 700+ samples. The position of the ten marker lipids is indicated and the clusters are emphasized in the hierarchical trees (top and left). The colored bars (top) indicate the phospholipid class off each lipid species, the total number of unsaturations in the acyl chains and whether the acyl chains are linear of contain a cyclopropane moiety.
Metabolites 11 00286 g004
Metabolites 11 00286 g005 550
Figure 5. Boxplots demonstrating the effect of alteration of FabH levels on key properties of membrane lipids. KO and OV represent gene knock-out and overexpression, respectively. EV represents the wild-type (WT) transformed with an empty overexpression vector. Panels on the right (+buoh) were obtained from cells grown under identical conditions as the control (CTR) but in the presence of 0.5% (v/v) n-butanol. It should be noted that acyl properties were calculated from the data in supplementary Table S1 and that the values at y-axis should not be treated as absolute. A discrepancy may result from the fact that differences in ionization efficiency can exist between molecular species. Differences between the two control conditions (wild type WT and empty vector EV) were non-significant in all experiments shown. Statistical significances (based on uncorrected p-values) of KO and OV are shown with respect to their controls (WT and EV, respectively). **: p < 0.01; ***: p < 0.001.
Figure 5. Boxplots demonstrating the effect of alteration of FabH levels on key properties of membrane lipids. KO and OV represent gene knock-out and overexpression, respectively. EV represents the wild-type (WT) transformed with an empty overexpression vector. Panels on the right (+buoh) were obtained from cells grown under identical conditions as the control (CTR) but in the presence of 0.5% (v/v) n-butanol. It should be noted that acyl properties were calculated from the data in supplementary Table S1 and that the values at y-axis should not be treated as absolute. A discrepancy may result from the fact that differences in ionization efficiency can exist between molecular species. Differences between the two control conditions (wild type WT and empty vector EV) were non-significant in all experiments shown. Statistical significances (based on uncorrected p-values) of KO and OV are shown with respect to their controls (WT and EV, respectively). **: p < 0.01; ***: p < 0.001.
Metabolites 11 00286 g005
Metabolites 11 00286 g006 550
Figure 6. Specific growth rates and optical densities at steady state. Values of genetic constructs are compared to controls. Only data that were statistically significant (p < 0.05) are shown. Genes discussed in the text are indicated by larger font size.
Figure 6. Specific growth rates and optical densities at steady state. Values of genetic constructs are compared to controls. Only data that were statistically significant (p < 0.05) are shown. Genes discussed in the text are indicated by larger font size.
Metabolites 11 00286 g006
Table
Table 1. Most informative lipid species for the classification of genetic composition and n-butanol exposure (yes/no) as identified by random forest analysis. Percentages are given as contributions to the total mass spectrometry signal.
Table 1. Most informative lipid species for the classification of genetic composition and n-butanol exposure (yes/no) as identified by random forest analysis. Percentages are given as contributions to the total mass spectrometry signal.
Lipid SpeciesAverage (%)Max. (%)Min. (%)Genes 1 (Extremes)
PE 34:17.6717.410.72glpD (OV); fabH (OV)
PE 32:15.7524.820.72cfa (KO); cfa (OV)
PE 36:24.2012.090.84fabH (KO); fabH (OV)
PG 35:0c12.6010.470.20glpR (KO); aas (OV)
PG 31:0c10.351.210.05fabH (OV); fabH (KO)
aPG 50:20.262.940.03fadL (OV); fabF (OV)
aPE 49:0c10.150.910.01aas (OV); cfa (KO)
aPG 49:1c10.101.000.01pldA (OV); fabH (KO)
PBut 33:0c10.072.990.00clsB (OV); fadD (KO)
aPE 50:10.071.340.00aas (OV); clsA (OV)
1 The genes of which altered expression resulted in highest or lowest abundance are given.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jeucken, A.; Zhou, M.; Wösten, M.M.S.M.; Brouwers, J.F. Control of n-Butanol Induced Lipidome Adaptations in E. coli. Metabolites 2021, 11, 286. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11050286

AMA Style

Jeucken A, Zhou M, Wösten MMSM, Brouwers JF. Control of n-Butanol Induced Lipidome Adaptations in E. coli. Metabolites. 2021; 11(5):286. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11050286

Chicago/Turabian Style

Jeucken, Aike, Miaomiao Zhou, Marc M. S. M. Wösten, and Jos F. Brouwers. 2021. "Control of n-Butanol Induced Lipidome Adaptations in E. coli" Metabolites 11, no. 5: 286. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11050286

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop