The Sphinx and the egg: Evolutionary enigmas of the (glyco)sphingolipid biosynthetic pathway

In eukaryotes, the de novo synthesis of sphingolipids (SLs) consists of multiple sequential steps which are compartmentalized between the endoplasmic reticulum and the Golgi apparatus. Studies over many decades have identified the enzymes in the pathway, their localization, topology and an array of regulatory mechanisms. However, little is known about the evolutionary forces that underly the generation of this complex pathway or of its anteome, i.e., the metabolic pathways that converge on the SL biosynthetic pathway and are essential for its activity. After briefly describing the pathway, we discuss the mechanisms by which the enzymes of the SL biosynthetic pathway are targeted to their different subcellular locations, how the pathway per se may have evolved, including its compartmentalization, and the relationship of the pathway to eukaryogenesis. We discuss the circular interdependence of the evolution of the SL pathway, and comment on whether current Darwinian evolutionary models are able to provide genuine mechanistic insight into how the pathway came into being.

Ancestral Sequence Reconstruction for Exploring Alkaloid Evolution

The complex and bioactive monoterpene indole alkaloids (MIAs) found in Catharanthus roseus and related species are the products of many millions of years of evolution through mutation and natural selection. Ancestral sequence reconstruction (ASR) is a method that combines phylogenetic analysis and experimental biochemistry to infer details about past events in protein evolution. Here, I propose that ASR could be leveraged to understand how enzymes catalyzing the formation of complex alkaloids arose over evolutionary time. I discuss the steps of ASR, including sequence selection, multiple sequence alignment, tree inference, and the generation and characterization of inferred ancestral enzymes.

Divergent patterns of selection on metabolite levels and gene expression

Background

Natural selection can act on multiple genes in the same pathway, leading to polygenic adaptation. For example, adaptive changes were found to down-regulate six genes involved in ergosterol biosynthesis—an essential pathway targeted by many antifungal drugs—in some strains of the yeast Saccharomyces cerevisiae. However, the impact of this polygenic adaptation on metabolite levels was unknown. Here, we performed targeted mass spectrometry to measure the levels of eight metabolites in this pathway in 74 yeast strains from a genetic cross.

Results

Through quantitative trait locus (QTL) mapping we identified 19 loci affecting ergosterol pathway metabolite levels, many of which overlap loci that also impact gene expression within the pathway. We then used the recently developed v-test, which identified selection acting upon three metabolite levels within the pathway, none of which were predictable from the gene expression adaptation.

Conclusions

These data showed that effects of selection on metabolite levels were complex and not predictable from gene expression data. This suggests that a deeper understanding of metabolism is necessary before we can understand the impacts of even relatively straightforward gene expression adaptations on metabolic pathways.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01915-5.

Evolution analysis of flavor-active compounds during artificial fermentation of Pu-erh tea

For the urgent need for fermentation control and product quality improvement of Pu-erh tea, gas chromatography-mass spectrometry and odor activity value (OAV) were used to comprehensively investigate the flavor-active compounds during artificial fermentation of Pu-erh tea. A flavor wheel was constructed to expound the sensory attributes evolution during fermentation. With an increased total volatiles content, 43 were significantly up-regulated and 30 were down-regulated among 131 detected volatiles. Key active compounds of three aroma types, namely fresh fragrance, fruit-fungus fragrance and stale-Qu fragrance, were analyzed based on OAV. β-damascenone was firstly found contributing most to the aroma of Pu-erh tea, followed by 1,2,3-methoxybenzene and (E,E)-2,4-nonadienal. γ-terpinene, linalool, 1,2,4-trimethoxybenzene, 1,2,3-trimethoxybenzene, and 4-ethylveratrol were identified as the potential markers responsible for aroma differences among three fermentation stages. Finally the metabolic evolution of key flavor-active compounds were systematically summarized. This study provides significant guidance in fermentation control and new product development of Pu-erh tea.

Tracing molecular properties throughout evolution: A chemoinformatic approach

Evolution of metabolism is a longstanding yet unresolved question, and several hypotheses were proposed to address this complex process from a Darwinian point of view. Modern statistical bioinformatic approaches targeted to the comparative analysis of genomes are being used to detect signatures of natural selection at the gene and population level, as an attempt to understand the origin of primordial metabolism and its expansion. These studies, however, are still mainly centered on genes and the proteins they encode, somehow neglecting the small organic chemicals that support life processes. In this work, we selected steroids as an ancient family of metabolites widely distributed in all eukaryotes and applied unsupervised machine learning techniques to reveal the traits that natural selection has imprinted on molecular properties throughout the evolutionary process. Our results clearly show that sterols, the primal steroids that first appeared, have more conserved properties and that, from then on, more complex compounds with increasingly diverse properties have emerged, suggesting that chemical diversification parallels the expansion of biological complexity. In a wider context, these findings highlight the worth of chemoinformatic approaches to a better understanding the evolution of metabolism.

The rise of diversity in metabolic platforms across the Candidate Phyla Radiation

BACKGROUND

A unifying feature of the bacterial Candidate Phyla Radiation (CPR) is a limited and highly variable repertoire of biosynthetic capabilities. However, the distribution of metabolic traits across the CPR and the evolutionary processes underlying them are incompletely resolved.

RESULTS

Here, we selected ~ 1000 genomes of CPR bacteria from diverse environments to construct a robust internal phylogeny that was consistent across two unlinked marker sets. Mapping of glycolysis, the pentose phosphate pathway, and pyruvate metabolism onto the tree showed that some components of these pathways are sparsely distributed and that similarity between metabolic platforms is only partially predicted by phylogenetic relationships. To evaluate the extent to which gene loss and lateral gene transfer have shaped trait distribution, we analyzed the patchiness of gene presence in a phylogenetic context, examined the phylogenetic depth of clades with shared traits, and compared the reference tree topology with those of specific metabolic proteins. While the central glycolytic pathway in CPR is widely conserved and has likely been shaped primarily by vertical transmission, there is evidence for both gene loss and transfer especially in steps that convert glucose into fructose 1,6-bisphosphate and glycerate 3P into pyruvate. Additionally, the distribution of Group 3 and Group 4-related NiFe hydrogenases is patchy and suggests multiple events of ancient gene transfer.

CONCLUSIONS

We infer that patterns of gene gain and loss in CPR, including acquisition of accessory traits in independent transfer events, could have been driven by shifts in host-derived resources and led to sparse but varied genetic inventories.

Metabolic evolution and a comparative omics analysis of Corynebacterium glutamicum for putrescine production

Putrescine is widely used in the industrial production of bioplastics, pharmaceuticals, agrochemicals, and surfactants. Because the highest titer of putrescine is much lower than that of its precursor L-ornithine reported in microorganisms to date, further work is needed to increase putrescine production in Corynebacterium glutamicum. We first compared 7 ornithine decarboxylase genes and found that the Enterobacter cloacae ornithine decarboxylase gene speC1 was most suitable for putrescine production in C. glutamicum. Increasing NADPH availability and blocking putrescine oxidation and acetylation were chosen as targets for metabolic engineering. The putrescine producer C. glutamicum PUT4 was first constructed by deleting puo, butA and snaA genes, and replacing the fabG gene with E. cloacae speC1. After adaptive evolution with C. glutamicum PUT4, the evolved strain C. glutamicum PUT-ALE, which produced an 96% higher amount of putrescine compared to the parent strain, was obtained. The whole genome resequencing indicates that the SNPs located in the odhA coding region may be associated with putrescine production. The comparative proteomic analysis reveals that the pentose phosphate and anaplerotic pathway, the glyoxylate cycle, and the ornithine biosynthetic pathway were upregulated in the evolved strain C. glutamicum PUT-ALE. The aspartate family, aromatic, and branched chain amino acid and fatty acid biosynthetic pathways were also observed to be downregulated in C. glutamicum PUT-ALE. Reducing OdhA activity by replacing the odhA native start codon GTG with TTG and overexpression of cgmA or pyc458 further improved putrescine production. Repressing the carB, ilvH, ilvB and aroE expression via CRISPRi also increased putrescine production by 5, 9, 16 and 19%, respectively.

Efficient metabolic evolution of engineered Yarrowia lipolytica for succinic acid production using a glucose-based medium in an in situ fibrous bioreactor under low-pH condition

BACKGROUND

Alkali used for pH control during fermentation and acidification for downstream recovery of succinic acid (SA) are the two largest cost contributors for bio-based SA production. To promote the commercialization process of fermentative SA, the development of industrially important microorganisms that can tolerate low pH has emerged as a crucial issue.

RESULTS

In this study, an in situ fibrous bed bioreactor (isFBB) was employed for the metabolic evolution for selection of Y. lipolytica strain that can produce SA at low pH using glucose-based medium. An evolved strain named Y. lipolytica PSA3.0 that could produce SA with a titer of 19.3 g/L, productivity of 0.52 g/L/h, and yield of 0.29 g/g at pH 3.0 from YPD was achieved. The enzyme activity analysis demonstrated that the pathway from pyruvate to acetate was partially blocked in Y. lipolytica PSA3.0 after the evolution, which is beneficial to cell growth and SA production at low pH. When free-cell batch fermentations were performed using the parent and evolved strains separately, the evolved strain PSA3.0 produced 18.4 g/L SA with a yield of 0.23 g/g at pH 3.0. Although these values were lower than that obtained by the parent strain PSA02004 at its optimal pH 6.0, which were 25.2 g/L and 0.31 g/g, respectively, they were 4.8 and 4.6 times higher than that achieved by PSA02004 at pH 3.0. By fed-batch fermentation, the resultant SA titer of 76.8 g/L was obtained, which is the highest value that ever achieved from glucose-based medium at low pH, to date. When using mixed food waste (MFW) hydrolysate as substrate, 18.9 g/L SA was produced with an SA yield of 0.38 g/g, which demonstrates the feasibility of using low-cost glucose-based hydrolysate for SA production by Y. lipolytica in a low-pH environment.

CONCLUSIONS

This study presents an effective and efficient strategy for the evolution of Y. lipolytica for SA production under low-pH condition for the first time. The isFBB was demonstrated to improve the metabolic evolution efficiency of Y. lipolytica to the acidic condition. Moreover, the acetate accumulation was found to be the major reason for the inhibition of SA production at low pH by Y. lipolytica, which suggested the direction for further metabolic modification of the strain for improved SA production. Furthermore, the evolved strain Y. lipolytica PSA3.0 was demonstrated to utilize glucose-rich hydrolysate from MFW for fermentative SA production at low pH. Similarly, Y. lipolytica PSA3.0 is expected to utilize the glucose-rich hydrolysate generated from other carbohydrate-rich waste streams for SA production. This study paves the way for the commercialization of bio-based SA and contributes to the sustainable development of a green economy.

Re-engineering Escherichia coli KJ122 to enhance the utilization of xylose and xylose/glucose mixture for efficient succinate production in mineral salt medium

Escherichia coli KJ122 was previously engineered to produce high concentration and yield of succinate in mineral salt medium containing glucose and sucrose under anaerobic conditions. However, this strain does not efficiently utilize xylose. To improve the xylose uptake and utilization in the strain KJ122, xylFGH and xylE genes were individually and simultaneously deleted. E. coli KJ12201 (KJ122::ΔxylFGH) exhibited superior abilities in growth, xylose consumption, and succinate production compared to those of the parental strain KJ122. However, E. coli KJ12202 (KJ122::ΔxylE) lessened xylose consumption due to an ATP deficit for metabolizing xylose thus making succinate production from xylose not preferable. Moreover, E. coli KJ12203 (KJ122::ΔxylFGHΔxylE) exhibited an impaired growth on xylose due to lacking of xylose transporters. After performing metabolic evolution, the evolved KJ12201-14T strain exhibited a great improvement in succinate production from pure xylose with higher concentration and productivity about 18 and 21%, respectively, compared to KJ12201 strain. During fed-batch fermentation, KJ12201-14T also produced succinate from xylose at a concentration, yield, and overall productivity of 84.6 ± 0.7 g/L, 0.86 ± 0.01 g/g and 1.01 ± 0.01 g/L/h, respectively. KJ12201 and KJ12201-14T strains co-utilized glucose/xylose mixture without catabolite repression. Both strains produced succinate from glucose/xylose mixture at concentration, yield, and overall and specific productivities of about 85 g/L, 0.85 g/g, 0.70 g/L/h, and 0.44 g/gCDW/h, respectively. Based on our results, KJ12201 and KJ12201-14T strains exhibited a greater performance in succinate production from xylose containing medium than those of other published works. They would be potential strains for the economic bio-based succinate production from xylose.

Switch on a more efficient pyruvate synthesis pathway based on transcriptome analysis and metabolic evolution

Due to the decrease of intracellular NADH availability, gluconate metabolism is more conducive to pyruvate production than glucose. Transcriptome analysis revealed that the Entner-Doudoroff (ED) pathway was activated by gluconate in Escherichia coli YP211 (MG1655 ΔldhA ΔpflB Δpta-ackA ΔpoxB Δppc ΔfrdBC). To construct a new pyruvate producing strain with glucose metabolism via ED pathway, the genes ppsA, ptsG, pgi and gnd were deleted sequentially to reduce the demand for PEP and block the Embden-Meyerhor-Parnas pathway and Pentose-Phosphate pathway. After nearly 1000 generations of growth-based selection, the evolved strain YP404 was isolated and the ED pathway was proved to be activated as the primary glycolytic pathway. Comparing with YP211, the pyruvate concentration and yield increased by 59% and 10.1%, respectively. In fed-batch fermentation, the pyruvate concentration reached 83.5 g l^-1 with a volumetric productivity of 2.3 g l^-1 h^-1. This was the first time to produce pyruvate via ED pathway, and prove that this was a more effective way.

Natural selection drove metabolic specialization of the chromatophore in Paulinella chromatophora

Background

Genome degradation of host-restricted mutualistic endosymbionts has been attributed to inactivating mutations and genetic drift while genes coding for host-relevant functions are conserved by purifying selection. Unlike their free-living relatives, the metabolism of mutualistic endosymbionts and endosymbiont-originated organelles is specialized in the production of metabolites which are released to the host. This specialization suggests that natural selection crafted these metabolic adaptations. In this work, we analyzed the evolution of the metabolism of the chromatophore of Paulinella chromatophora by in silico modeling. We asked whether genome reduction is driven by metabolic engineering strategies resulted from the interaction with the host. As its widely known, the loss of enzyme coding genes leads to metabolic network restructuring sometimes improving the production rates. In this case, the production rate of reduced-carbon in the metabolism of the chromatophore.

Results

We reconstructed the metabolic networks of the chromatophore of P. chromatophora CCAC 0185 and a close free-living relative, the cyanobacterium Synechococcus sp. WH 5701. We found that the evolution of free-living to host-restricted lifestyle rendered a fragile metabolic network where >80% of genes in the chromatophore are essential for metabolic functionality. Despite the lack of experimental information, the metabolic reconstruction of the chromatophore suggests that the host provides several metabolites to the endosymbiont. By using these metabolites as intracellular conditions, in silico simulations of genome evolution by gene lose recover with 77% accuracy the actual metabolic gene content of the chromatophore. Also, the metabolic model of the chromatophore allowed us to predict by flux balance analysis a maximum rate of reduced-carbon released by the endosymbiont to the host. By inspecting the central metabolism of the chromatophore and the free-living cyanobacteria we found that by improvements in the gluconeogenic pathway the metabolism of the endosymbiont uses more efficiently the carbon source for reduced-carbon production. In addition, our in silico simulations of the evolutionary process leading to the reduced metabolic network of the chromatophore showed that the predicted rate of released reduced-carbon is obtained in less than 5% of the times under a process guided by random gene deletion and genetic drift. We interpret previous findings as evidence that natural selection at holobiont level shaped the rate at which reduced-carbon is exported to the host. Finally, our model also predicts that the ABC phosphate transporter (pstSACB) which is conserved in the genome of the chromatophore of P. chromatophora strain CCAC 0185 is a necessary component to release reduced-carbon molecules to the host.

Conclusion

Our evolutionary analysis suggests that in the case of Paulinella chromatophora natural selection at the holobiont level played a prominent role in shaping the metabolic specialization of the chromatophore. We propose that natural selection acted as a “metabolic engineer” by favoring metabolic restructurings that led to an increased release of reduced-carbon to the host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0947-6) contains supplementary material, which is available to authorized users.

A Workflow for Studying Specialized Metabolism in Nonmodel Eukaryotic Organisms

Eukaryotes contain a diverse tapestry of specialized metabolites, many of which are of significant pharmaceutical and industrial importance to humans. Nevertheless, exploration of specialized metabolic pathways underlying specific chemical traits in nonmodel eukaryotic organisms has been technically challenging and historically lagged behind that of the bacterial systems. Recent advances in genomics, metabolomics, phylogenomics, and synthetic biology now enable a new workflow for interrogating unknown specialized metabolic systems in nonmodel eukaryotic hosts with greater efficiency and mechanistic depth. This chapter delineates such workflow by providing a collection of state-of-the-art approaches and tools, ranging from multiomics-guided candidate gene identification to in vitro and in vivo functional and structural characterization of specialized metabolic enzymes. As already demonstrated by several recent studies, this new workflow opens up a gateway into the largely untapped world of natural product biochemistry in eukaryotes.

Systematic development of biomass overproducing Scheffersomyces stipitis for high-cell-density fermentations

The development of economically feasible bio-based process requires efficient cell factories capable of producing the desired product at high titer under high-cell-density fermentation. Herein we present a combinatorial approach based on systems metabolic engineering and metabolic evolution for the development of efficient biomass-producing strain. Systems metabolic engineering guided by flux balance analysis (FBA) was first employed to rationally design mutant strains of Scheffersomyces stipitis with high biomass yield. By experimentally implementing these mutations, the biomass yield was improved by 30% in GPD1, 25% in TKL1, 30% in CIT1, and 44% in ZWF1 overexpressed mutants compared to wild-type. These designed mutants were further fine-tuned through metabolic evolution resulting in the maximal biomass yield of 0.49 g-cdw/g-glucose, which matches well with predicted yield phenotype. The constructed mutants are beneficial for biotechnology applications dealing with high cell titer cultivations. This work demonstrates a solid confirmation of systems metabolic engineering in combination with metabolic evolution approach for efficient strain development, which could assist in rapid optimization of cell factory for an economically viable and sustainable bio-based process.