Evolutionary Conservation of Bacterial Essential Metabolic Genes across All Bacterial Culture Media

One of the basic postulates of molecular evolution is that functionally important genes should evolve slower than genes of lesser significance. Essential genes, whose knockout leads to a lethal phenotype are considered of high functional importance, yet whether they are truly more conserved than nonessential genes has been the topic of much debate, fuelled by a host of contradictory findings. Here we conduct the first large-scale study utilizing genome-scale metabolic modeling and spanning many bacterial species, which aims to answer this question. Using the novel Media Variation Analysis, we examine the range of conservation of essential vs. nonessential metabolic genes in a given species across all possible media. We are thus able to obtain for the first time, exact upper and lower bounds on the levels of differential conservation of essential genes for each of the species studied. The results show that bacteria do exhibit an overall tendency for differential conservation of their essential genes vs. their non-essential ones, yet this tendency is highly variable across species. We show that the model bacterium E. coli K12 may or may not exhibit differential conservation of essential genes depending on its growth medium, shedding light on previous experimental studies showing opposite trends.     

Pervasive and persistent redundancy among duplicated genes in yeast

The loss of functional redundancy is the key process in the evolution of duplicated genes. Here we systematically assess the extent of functional redundancy among a large set of duplicated genes in Saccharomyces cerevisiae. We quantify growth rate in rich medium for a large number of S. cerevisiae strains that carry single and double deletions of duplicated and singleton genes. We demonstrate that duplicated genes can maintain substantial redundancy for extensive periods of time following duplication (~100 million years). We find high levels of redundancy among genes duplicated both via the whole genome duplication and via smaller scale duplications. Further, we see no evidence that two duplicated genes together contribute to fitness in rich medium substantially beyond that of their ancestral progenitor gene. We argue that duplicate genes do not often evolve to behave like singleton genes even after very long periods of time.     

Physiologically Shrinking the Solution Space of a Saccharomyces cerevisiae Genome-Scale Model Suggests the Role of the Metabolic Network in Shaping Gene Expression Noise

Sampling the solution space of genome-scale models is generally conducted to determine the feasible region for metabolic flux distribution. Because the region for actual metabolic states resides only in a small fraction of the entire space, it is necessary to shrink the solution space to improve the predictive power of a model. A common strategy is to constrain models by integrating extra datasets such as high-throughput datasets and C¹³-labeled flux datasets. However, studies refining these approaches by performing a meta-analysis of massive experimental metabolic flux measurements, which are closely linked to cellular phenotypes, are limited. In the present study, experimentally identified metabolic flux data from 96 published reports were systematically reviewed. Several strong associations among metabolic flux phenotypes were observed. These phenotype-phenotype associations at the flux level were quantified and integrated into a Saccharomyces cerevisiae genome-scale model as extra physiological constraints. By sampling the shrunken solution space of the model, the metabolic flux fluctuation level, which is an intrinsic trait of metabolic reactions determined by the network, was estimated and utilized to explore its relationship to gene expression noise. Although no correlation was observed in all enzyme-coding genes, a relationship between metabolic flux fluctuation and expression noise of genes associated with enzyme-dosage sensitive reactions was detected, suggesting that the metabolic network plays a role in shaping gene expression noise. Such correlation was mainly attributed to the genes corresponding to non-essential reactions, rather than essential ones. This was at least partially, due to regulations underlying the flux phenotype-phenotype associations. Altogether, this study proposes a new approach in shrinking the solution space of a genome-scale model, of which sampling provides new insights into gene expression noise.     

Genome-Scale NAD(H/+) Availability Patterns as a Differentiating Feature between Saccharomyces cerevisiae and Scheffersomyces stipitis in Relation to Fermentative Metabolism

Scheffersomyces stipitis is a yeast able to ferment pentoses to ethanol, unlike Saccharomyces cerevisiae, it does not present the so-called overflow phenomenon. Metabolic features characterizing the presence or not of this phenomenon have not been fully elucidated. This work proposes that genome-scale metabolic response to variations in NAD(H/⁺) availability characterizes fermentative behavior in both yeasts. Thus, differentiating features in S. stipitis and S. cerevisiae were determined analyzing growth sensitivity response to changes in available reducing capacity in relation to ethanol production capacity and overall metabolic flux span. Using genome-scale constraint-based metabolic models, phenotypic phase planes and shadow price analyses, an excess of available reducing capacity for growth was found in S. cerevisiae at every metabolic phenotype where growth is limited by oxygen uptake, while in S. stipitis this was observed only for a subset of those phenotypes. Moreover, by using flux variability analysis, an increased metabolic flux span was found in S. cerevisiae at growth limited by oxygen uptake, while in S. stipitis flux span was invariant. Therefore, each yeast can be characterized by a significantly different metabolic response and flux span when growth is limited by oxygen uptake, both features suggesting a higher metabolic flexibility in S. cerevisiae. By applying an optimization-based approach on the genome-scale models, three single reaction deletions were found to generate in S. stipitis the reducing capacity availability pattern found in S. cerevisiae, two of them correspond to reactions involved in the overflow phenomenon. These results show a close relationship between the growth sensitivity response given by the metabolic network and fermentative behavior.     

Crabtree/Warburg-like aerobic xylose fermentation by engineered Saccharomyces cerevisiae

Bottlenecks in the efficient conversion of xylose into cost-effective biofuels have limited the widespread use of plant lignocellulose as a renewable feedstock. The yeast Saccharomyces cerevisiae ferments glucose into ethanol with such high metabolic flux that it ferments high concentrations of glucose aerobically, a trait called the Crabtree/Warburg Effect. In contrast to glucose, most engineered S. cerevisiae strains do not ferment xylose at economically viable rates and yields, and they require respiration to achieve sufficient xylose metabolic flux and energy return for growth aerobically. Here, we evolved respiration-deficient S. cerevisiae strains that can grow on and ferment xylose to ethanol aerobically, a trait analogous to the Crabtree/Warburg Effect for glucose. Through genome sequence comparisons and directed engineering, we determined that duplications of genes encoding engineered xylose metabolism enzymes, as well as TKL1, a gene encoding a transketolase in the pentose phosphate pathway, were the causative genetic changes for the evolved phenotype. Reengineered duplications of these enzymes, in combination with deletion mutations in HOG1, ISU1, GRE3, and IRA2, increased the rates of aerobic and anaerobic xylose fermentation. Importantly, we found that these genetic modifications function in another genetic background and increase the rate and yield of xylose-to-ethanol conversion in industrially relevant switchgrass hydrolysate, indicating that these specific genetic modifications may enable the sustainable production of industrial biofuels from yeast. We propose a model for how key regulatory mutations prime yeast for aerobic xylose fermentation by lowering the threshold for overflow metabolism, allowing mutations to increase xylose flux and to redirect it into fermentation products.     

Metabolic Flexibility as a Major Predictor of Spatial Distribution in Microbial Communities

A better understand the ecology of microbes and their role in the global ecosystem could be achieved if traditional ecological theories can be applied to microbes. In ecology organisms are defined as specialists or generalists according to the breadth of their niche. Spatial distribution is often used as a proxy measure of niche breadth; generalists have broad niches and a wide spatial distribution and specialists a narrow niche and spatial distribution. Previous studies suggest that microbial distribution patterns are contrary to this idea; a microbial generalist genus (Desulfobulbus) has a limited spatial distribution while a specialist genus (Methanosaeta) has a cosmopolitan distribution. Therefore, we hypothesise that this counter-intuitive distribution within generalist and specialist microbial genera is a common microbial characteristic. Using molecular fingerprinting the distribution of four microbial genera, two generalists, Desulfobulbus and the methanogenic archaea Methanosarcina, and two specialists, Methanosaeta and the sulfate-reducing bacteria Desulfobacter were analysed in sediment samples from along a UK estuary. Detected genotypes of both generalist genera showed a distinct spatial distribution, significantly correlated with geographic distance between sites. Genotypes of both specialist genera showed no significant differential spatial distribution. These data support the hypothesis that the spatial distribution of specialist and generalist microbes does not match that seen with specialist and generalist large organisms. It may be that generalist microbes, while having a wider potential niche, are constrained, possibly by intrageneric competition, to exploit only a small part of that potential niche while specialists, with far fewer constraints to their niche, are more capable of filling their potential niche more effectively, perhaps by avoiding intrageneric competition. We suggest that these counter-intuitive distribution patterns may be a common feature of microbes in general and represent a distinct microbial principle in ecology, which is a real challenge if we are to develop a truly inclusive ecology.     

Enhancing sesquiterpene production in Saccharomyces cerevisiae through in silico driven metabolic engineering

A genome-scale metabolic model was used to identify new target genes for enhanced biosynthesis of sesquiterpenes in the yeast Saccharomyces cerevisiae. The effect of gene deletions on the flux distributions in the metabolic model of S. cerevisiae was assessed using OptGene as the modeling framework and minimization of metabolic adjustments (MOMA) as objective function.Deletion of NADPH-dependent glutamate dehydrogenase encoded by GDH1 was identified as the best target gene for the improvement of sesquiterpene biosynthesis in yeast. Deletion of this gene enhances the available NADPH in the cytosol for other NADPH requiring enzymes, including HMG-CoA reductase. However, since disruption of GDH1 impairs the ammonia utilization, simultaneous over-expression of the NADH-dependent glutamate dehydrogenase encoded by GDH2 was also considered in this study.Deletion of GDH1 led to an approximately 85% increase in the final cubebol titer. However, deletion of this gene also caused a significant decrease in the maximum specific growth rate. Over-expression of GDH2 did not show a further effect on the final cubebol titer but this alteration significantly improved the growth rate compared to the GDH1 deleted strain.     

Reconstruction and In Silico Analysis of Metabolic Network for an Oleaginous Yeast,Yarrowia lipolytica

With the emergence of energy scarcity, the use of renewable energy sources such as biodiesel is becoming increasingly necessary. Recently, many researchers have focused their minds on Yarrowia lipolytica, a model oleaginous yeast, which can be employed to accumulate large amounts of lipids that could be further converted to biodiesel. In order to understand the metabolic characteristics of Y. lipolytica at a systems level and to examine the potential for enhanced lipid production, a genome-scale compartmentalized metabolic network was reconstructed based on a combination of genome annotation and the detailed biochemical knowledge from multiple databases such as KEGG, ENZYME and BIGG. The information about protein and reaction associations of all the organisms in KEGG and Expasy-ENZYME database was arranged into an EXCEL file that can then be regarded as a new useful database to generate other reconstructions. The generated model iYL619_PCP accounts for 619 genes, 843 metabolites and 1,142 reactions including 236 transport reactions, 125 exchange reactions and 13 spontaneous reactions. The in silico model successfully predicted the minimal media and the growing abilities on different substrates. With flux balance analysis, single gene knockouts were also simulated to predict the essential genes and partially essential genes. In addition, flux variability analysis was applied to design new mutant strains that will redirect fluxes through the network and may enhance the production of lipid. This genome-scale metabolic model of Y. lipolytica can facilitate system-level metabolic analysis as well as strain development for improving the production of biodiesels and other valuable products by Y. lipolytica and other closely related oleaginous yeasts.     

Effect of Insertion Site and Metabolic Load on the Environmental Fitness of a Genetically Modified Pseudomonas fluorescens Isolate

An isolate of Pseudomonas fluorescens (SBW25) was modified with different marker genes (lacZY, aph-1, and xylE). These marker genes were inserted singly or in combination into two separate (1 Mbp apart) and presumably nonessential sites (-6- and Ee) on the chromosome of SBW25. This allowed the production of a range of genetically modified SBW25 variants that differed with respect to insertion site of the marker genes and metabolic burden. The environmental fitness of the different SBW25 variants was tested in soil, in the rhizosphere of wheat and pea, and on the phylloplane of wheat. Reduced environmental fitness of the different variants was mainly attributed to the extra metabolic burden of novel gene expression, whereas choice of insertion site was of little significance. Changes in environmental fitness were dependent on the environmental conditions; an environment, such as soil, with a low microbial carrying capacity had a negative effect on the environmental fitness of variants with a large metabolic load. In environments with a larger carrying capacity, such as the rhizosphere of pea, environmental fitness of variants with a large metabolic load was not significantly different from that of variants with a smaller metabolic burden.     

Physiological and Transcriptional Responses of Different Industrial Microbes at Near-Zero Specific Growth Rates

The current knowledge of the physiology and gene expression of industrially relevant microorganisms is largely based on laboratory studies under conditions of rapid growth and high metabolic activity. However, in natural ecosystems and industrial processes, microbes frequently encounter severe calorie restriction. As a consequence, microbial growth rates in such settings can be extremely slow and even approach zero. Furthermore, uncoupling microbial growth from product formation, while cellular integrity and activity are maintained, offers perspectives that are economically highly interesting. Retentostat cultures have been employed to investigate microbial physiology at (near-)zero growth rates. This minireview compares information from recent physiological and gene expression studies on retentostat cultures of the industrially relevant microorganisms Lactobacillus plantarum, Lactococcus lactis, Bacillus subtilis, Saccharomyces cerevisiae, and Aspergillus niger. Shared responses of these organisms to (near-)zero growth rates include increased stress tolerance and a downregulation of genes involved in protein synthesis. Other adaptations, such as changes in morphology and (secondary) metabolite production, were species specific. This comparison underlines the industrial and scientific significance of further research on microbial (near-)zero growth physiology.     

Pan-cancer analysis of the metabolic reaction network

Metabolic reprogramming is considered a hallmark of malignant transformation. However, it is not clear whether the network of metabolic reactions expressed by cancers of different origin differ from each other or from normal human tissues. In this study, we reconstructed functional and connected genome-scale metabolic models for 917 primary tumor samples across 13 types based on the probability of expression for 3765 reference metabolic genes in the sample. This network-centric approach revealed that tumor metabolic networks are largely similar in terms of accounted reactions, despite diversity in the expression of the associated genes. On average, each network contained 4721 reactions, of which 74% were core reactions (present in >95% of all models). Whilst 99.3% of the core reactions were classified as housekeeping also in normal tissues, we identified reactions catalyzed by ARG2, RHAG, SLC6 and SLC16 family gene members, and PTGS1 and PTGS2 as core exclusively in cancer. These findings were subsequently replicated in an independent validation set of 3388 genome-scale metabolic models. The remaining 26% of the reactions were contextual reactions. Their inclusion was dependent in one case (GLS2) on the absence of TP53 mutations and in 94.6% of cases on differences in cancer types. This dependency largely resembled differences in expression patterns in the corresponding normal tissues, with some exceptions like the presence of the NANP-encoded reaction in tumors not from the female reproductive system or of the SLC5A9-encoded reaction in kidney-pancreatic-colorectal tumors. In conclusion, tumors expressed a metabolic network virtually overlapping the matched normal tissues, raising the possibility that metabolic reprogramming simply reflects cancer cell plasticity to adapt to varying conditions thanks to redundancy and complexity of the underlying metabolic networks. At the same time, the here uncovered exceptions represent a resource to identify selective liabilities of tumor metabolism.     

Essential Plasticity and Redundancy of Metabolism Unveiled by Synthetic Lethality Analysis

We unravel how functional plasticity and redundancy are essential mechanisms underlying the ability to survive of metabolic networks. We perform an exhaustive computational screening of synthetic lethal reaction pairs in Escherichia coli in a minimal medium and we find that synthetic lethal pairs divide in two different groups depending on whether the synthetic lethal interaction works as a backup or as a parallel use mechanism, the first corresponding to essential plasticity and the second to essential redundancy. In E. coli, the analysis of pathways entanglement through essential redundancy supports the view that synthetic lethality affects preferentially a single function or pathway. In contrast, essential plasticity, the dominant class, tends to be inter-pathway but strongly localized and unveils Cell Envelope Biosynthesis as an essential backup for Membrane Lipid Metabolism. When comparing E. coli and Mycoplasma pneumoniae, we find that the metabolic networks of the two organisms exhibit a large difference in the relative importance of plasticity and redundancy which is consistent with the conjecture that plasticity is a sophisticated mechanism that requires a complex organization. Finally, coessential reaction pairs are explored in different environmental conditions to uncover the interplay between the two mechanisms. We find that synthetic lethal interactions and their classification in plasticity and redundancy are basically insensitive to medium composition, and are highly conserved even when the environment is enriched with nonessential compounds or overconstrained to decrease maximum biomass formation.     

Optimization based automated curation of metabolic reconstructions

Background  

Currently, there exists tens of different microbial and eukaryotic metabolic reconstructions (e.g., Escherichia coli, Saccharomyces cerevisiae, Bacillus subtilis) with many more under development. All of these reconstructions are inherently incomplete with some functionalities missing due to the lack of experimental and/or homology information. A key challenge in the automated generation of genome-scale reconstructions is the elucidation of these gaps and the subsequent generation of hypotheses to bridge them.     

Characteristics of Plant Essential Genes Allow for within- and between-Species Prediction of Lethal Mutant Phenotypes

Essential genes represent critical cellular components whose disruption results in lethality. Characteristics shared among essential genes have been uncovered in fungal and metazoan model systems. However, features associated with plant essential genes are largely unknown and the full set of essential genes remains to be discovered in any plant species. Here, we show that essential genes in Arabidopsis thaliana have distinct features useful for constructing within- and cross-species prediction models. Essential genes in A. thaliana are often single copy or derived from older duplications, highly and broadly expressed, slow evolving, and highly connected within molecular networks compared with genes with nonlethal mutant phenotypes. These gene features allowed the application of machine learning methods that predicted known lethal genes as well as an additional 1970 likely essential genes without documented phenotypes. Prediction models from A. thaliana could also be applied to predict Oryza sativa and Saccharomyces cerevisiae essential genes. Importantly, successful predictions drew upon many features, while any single feature was not sufficient. Our findings show that essential genes can be distinguished from genes with nonlethal phenotypes using features that are similar across kingdoms and indicate the possibility for translational application of our approach to species without extensive functional genomic and phenomic resources.     

Synthetic Cytotoxicity: Digenic Interactions with TEL1/ATM Mutations Reveal Sensitivity to Low Doses of Camptothecin

Many tumors contain mutations that confer defects in the DNA-damage response and genome stability. DNA-damaging agents are powerful therapeutic tools that can differentially kill cells with an impaired DNA-damage response. The response to DNA damage is complex and composed of a network of coordinated pathways, often with a degree of redundancy. Tumor-specific somatic mutations in DNA-damage response genes could be exploited by inhibiting the function of a second gene product to increase the sensitivity of tumor cells to a sublethal concentration of a DNA-damaging therapeutic agent, resulting in a class of conditional synthetic lethality we call synthetic cytotoxicity. We used the Saccharomyces cerevisiae nonessential gene-deletion collection to screen for synthetic cytotoxic interactions with camptothecin, a topoisomerase I inhibitor, and a null mutation in TEL1, the S. cerevisiae ortholog of the mammalian tumor-suppressor gene, ATM. We found and validated 14 synthetic cytotoxic interactions that define at least five epistasis groups. One class of synthetic cytotoxic interaction was due to telomere defects. We also found that at least one synthetic cytotoxic interaction was conserved in Caenorhabditis elegans. We have demonstrated that synthetic cytotoxicity could be a useful strategy for expanding the sensitivity of certain tumors to DNA-damaging therapeutics.     

A Minimal Set of Glycolytic Genes Reveals Strong Redundancies in Saccharomyces cerevisiae Central Metabolism

As a result of ancestral whole-genome and small-scale duplication events, the genomes of Saccharomyces cerevisiae and many eukaryotes still contain a substantial fraction of duplicated genes. In all investigated organisms, metabolic pathways, and more particularly glycolysis, are specifically enriched for functionally redundant paralogs. In ancestors of the Saccharomyces lineage, the duplication of glycolytic genes is purported to have played an important role leading to S. cerevisiae''s current lifestyle favoring fermentative metabolism even in the presence of oxygen and characterized by a high glycolytic capacity. In modern S. cerevisiae strains, the 12 glycolytic reactions leading to the biochemical conversion from glucose to ethanol are encoded by 27 paralogs. In order to experimentally explore the physiological role of this genetic redundancy, a yeast strain with a minimal set of 14 paralogs was constructed (the “minimal glycolysis” [MG] strain). Remarkably, a combination of a quantitative systems approach and semiquantitative analysis in a wide array of growth environments revealed the absence of a phenotypic response to the cumulative deletion of 13 glycolytic paralogs. This observation indicates that duplication of glycolytic genes is not a prerequisite for achieving the high glycolytic fluxes and fermentative capacities that are characteristic of S. cerevisiae and essential for many of its industrial applications and argues against gene dosage effects as a means of fixing minor glycolytic paralogs in the yeast genome. The MG strain was carefully designed and constructed to provide a robust prototrophic platform for quantitative studies and has been made available to the scientific community.     

In silico analysis of bioethanol production from glucose/xylose mixtures during fed-batch fermentation of co-culture and mono-culture systems

This study presents a detailed in silico analysis of bioethanol production from glucose/xylose mixtures of various compositions by fed-batch co-culture and mono-culture fermentation of specialized microbes. The mono-culture consists of recombinant Saccharomyces cerevisise that can metabolize both hexose and pentose sugars while the co-culture system consists of substrate-selective microbes. Dynamic flux balance models based on available genome-scale reconstructions of the microorganisms have been used to analyze bioethanol production in fed-batch culture with constant feed rates and the maximization of ethanol productivity is addressed by computing optimal aerobic-anaerobic switching times. The simulation results clearly point to the superior performance of fed-batch fermentation of microbial co-culture against fed-batch fermentation of mono-culture for bioethanol production from glucose/xylose mixtures. A set of potential genetic engineering strategies for enhancement of S. cerevisiae and Escherichia coli strains performance have been identified. Such in silico predictions using genome-scale models provide valuable guidance for conducting in vivo metabolic engineering experiments.     

Gene Duplication in SACCHAROMYCES CEREVISIAE

Five indepdendent duplications of the acid-phosphatase (aphtase) structural gene (acp1) were recovered from chemostat populations of S. cerevisiae that were subject to selection for in vivo hyper-aphtase activity. Two of the duplications arose spontaneously. Three of them were induced by UV. All five of the duplication events involved the transpositioning of the aphtase structural gene, acp1, and all known genes distal to acp1 on the right arm of chromosome II, to the terminus of an arm of other unknown chromosomes. One of the five duplicated regions of the right arm of chromosome II was found to be transmitted mitotically and meiotically with very high fidelity. The other four duplicated regions of the right arm of chromosome II were found to be unstable, being lost at a rate of about 2% per mitosis. However, selection for increased fidelity of mitotic transmission was effective in one of these strains. No tandem duplications of the aphtase structural gene were found.