DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates.

The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.

Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Vibrio natriegens genome-scale modeling reveals insights into halophilic adaptations and resource allocation.

Vibrio natriegens is a Gram-negative bacterium with an exceptional growth rate that has the potential to become a standard biotechnological host for laboratory and industrial bioproduction. Despite this burgeoning interest, the current lack of organism-specific qualitative and quantitative computational tools has hampered the community's ability to rationally engineer this bacterium. In this study, we present the first genome-scale metabolic model (GSMM) of V. natriegens. The GSMM (iLC858) was developed using an automated draft assembly and extensive manual curation and was validated by comparing predicted yields, central metabolic fluxes, viable carbon substrates, and essential genes with empirical data. Mass spectrometry-based proteomics data confirmed the translation of at least 76% of the enzyme-encoding genes predicted to be expressed by the model during aerobic growth in a minimal medium. iLC858 was subsequently used to carry out a metabolic comparison between the model organism Escherichia coli and V. natriegens, leading to an analysis of the model architecture of V. natriegens' respiratory and ATP-generating system and the discovery of a role for a sodium-dependent oxaloacetate decarboxylase pump. The proteomics data were further used to investigate additional halophilic adaptations of V. natriegens. Finally, iLC858 was utilized to create a Resource Balance Analysis model to study the allocation of carbon resources. Taken together, the models presented provide useful computational tools to guide metabolic engineering efforts in V. natriegens.

Metabolic and mitochondria alterations induced by SARS-CoV-2 accessory proteins ORF3a, ORF9b, ORF9c and ORF10.

Antiviral signaling, immune response and cell metabolism are dysregulated by SARS-CoV-2, the causative agent of COVID-19. Here, we show that SARS-CoV-2 accessory proteins ORF3a, ORF9b, ORF9c and ORF10 induce a significant mitochondrial and metabolic reprogramming in A549 lung epithelial cells. While ORF9b, ORF9c and ORF10 induced largely overlapping transcriptomes, ORF3a induced a distinct transcriptome, including the downregulation of numerous genes with critical roles in mitochondrial function and morphology. On the other hand, all four ORFs altered mitochondrial dynamics and function, but only ORF3a and ORF9c induced a marked alteration in mitochondrial cristae structure. Genome-Scale Metabolic Models identified both metabolic flux reprogramming features both shared across all accessory proteins and specific for each accessory protein. Notably, a downregulated amino acid metabolism was observed in ORF9b, ORF9c and ORF10, while an upregulated lipid metabolism was distinctly induced by ORF3a. These findings reveal metabolic dependencies and vulnerabilities prompted by SARS-CoV-2 accessory proteins that may be exploited to identify new targets for intervention.

A 4-hydroxybenzoate 3-hydroxylase mutant enables 4-amino-3-hydroxybenzoic acid production from glucose in Corynebacterium glutamicum.

BACKGROUND: Microbial production of aromatic chemicals is an attractive method for obtaining high-performance materials from biomass resources. A non-proteinogenic amino acid, 4-amino-3-hydroxybenzoic acid (4,3-AHBA), is expected to be a precursor of highly functional polybenzoxazole polymers; however, methods for its microbial production have not been reported. In this study, we attempted to produce 4,3-AHBA from glucose by introducing 3-hydroxylation of 4-aminobenzoic acid (4-ABA) into the metabolic pathway of an industrially relevant bacterium, Corynebacterium glutamicum. RESULTS: Six different 4-hydroxybenzoate 3-hydroxylases (PHBHs) were heterologously expressed in C. glutamicum strains, which were then screened for the production of 4,3-AHBA by culturing with glucose as a carbon source. The highest concentration of 4,3-AHBA was detected in the strain expressing PHBH from Caulobacter vibrioides (CvPHBH). A combination of site-directed mutagenesis in the active site and random mutagenesis via laccase-mediated colorimetric assay allowed us to obtain CvPHBH mutants that enhanced 4,3-AHBA productivity under deep-well plate culture conditions. The recombinant C. glutamicum strain expressing CvPHBHM106A/T294S and having an enhanced 4-ABA biosynthetic pathway produced 13.5 g/L (88 mM) 4,3-AHBA and 0.059 g/L (0.43 mM) precursor 4-ABA in fed-batch culture using a nutrient-rich medium. The culture of this strain in the chemically defined CGXII medium yielded 9.8 C-mol% of 4,3-AHBA from glucose, corresponding to 12.8% of the theoretical maximum yield (76.8 C-mol%) calculated using a genome-scale metabolic model of C. glutamicum. CONCLUSIONS: Identification of PHBH mutants that could efficiently catalyze the 3-hydroxylation of 4-ABA in C. glutamicum allowed us to construct an artificial biosynthetic pathway capable of producing 4,3-AHBA on a gram-scale using glucose as the carbon source. These findings will contribute to a better understanding of enzyme-catalyzed regioselective hydroxylation of aromatic chemicals and to the diversification of biomass-derived precursors for high-performance materials.

Understanding the adaptive laboratory evolution of multiple stress-resistant yeast strains by genome-scale modeling.

Stress responses triggered by external exposures in adaptive laboratory evolution studies alter the ordinary behavior of cells, and the identification of the differences between the starting and the evolved strains would provide ideal strategies to obtain the desired strains. Metabolic networks are one of the most useful tools to analyze data for this purpose. This study integrates differential expression profiles of multiple Saccharomyces cerevisiae strains that have evolved in eight different stress conditions (ethanol, caffeine, coniferyl aldehyde, iron, nickel, phenylethanol, and silver) and enzyme kinetics into a genome-scale metabolic model of yeast, following a new enhanced method. Flux balance analysis, flux variability analysis, robustness, phenotype phase plane, minimization of metabolic adjustment, survivability, sensitivity analyses, and random sampling are conducted to identify the most common and divergent points within strains. Results were examined both individually and comparatively, and the target reactions, metabolites, and enzymes were identified. Our results showed that the models reconstructed by our methodology were able to simulate experimental conditions where efficient protein allocation was the main goal for survival under stressful conditions, and most of the metabolic changes in the adaptation process mainly arose from the differences in the metabolic reactions of energy maintenance (through coenzyme-A and FAD utilization), cell division (folate requirement of DNA synthesis), and cell wall formation (through sterol and ergosterol biosynthesis).

Genome-Scale Metabolic Modeling Combined with Transcriptome Profiling Provides Mechanistic Understanding of Streptococcus thermophilus CH8 Metabolism.

Streptococcus thermophilus is a lactic acid bacterium adapted toward growth in milk and is a vital component of starter cultures for milk fermentation. Here, we combine genome-scale metabolic modeling and transcriptome profiling to obtain novel metabolic insights into this bacterium. Notably, a refined genome-scale metabolic model (GEM) accurately representing S. thermophilus CH8 metabolism was developed. Modeling the utilization of casein as a nitrogen source revealed an imbalance in amino acid supply and demand, resulting in growth limitation due to the scarcity of specific amino acids, in particular sulfur amino acids. Growth experiments in milk corroborated this finding. A subtle interdependency of the redox balance and the secretion levels of the key metabolites lactate, formate, acetoin, and acetaldehyde was furthermore identified with the modeling approach, providing a mechanistic understanding of the factors governing the secretion product profile. As a potential effect of high expression of arginine biosynthesis genes, a moderate secretion of ornithine was observed experimentally, augmenting the proposed hypothesis of ornithine/putrescine exchange as part of the protocooperative interaction between S. thermophilus and Lactobacillus delbrueckii subsp. bulgaricus in yogurt. This study provides a foundation for future community modeling of food fermentations and rational development of starter strains with improved functionality. IMPORTANCE Streptococcus thermophilus is one the main organisms involved in the fermentation of milk and, increasingly, also in the fermentation of plant-based foods. The construction of a functional high-quality genome-scale metabolic model, in conjunction with in-depth transcriptome profiling with a focus on metabolism, provides a valuable resource for the improved understanding of S. thermophilus physiology. An example is the model-based prediction of the most significant route of synthesis for the characteristic yogurt flavor compound acetaldehyde and identification of metabolic principles governing the synthesis of other flavor compounds. Moreover, the systematic assessment of amino acid supply and demand during growth in milk provides insights into the key challenges related to nitrogen metabolism that is imposed on S. thermophilus and any other organism associated with the milk niche.

Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets.

Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.

Dissecting the metabolic reprogramming of maize root under nitrogen-deficient stress conditions.

The growth and development of maize (Zea mays L.) largely depends on its nutrient uptake through the root. Hence, studying its growth, response, and associated metabolic reprogramming to stress conditions is becoming an important research direction. A genome-scale metabolic model (GSM) for the maize root was developed to study its metabolic reprogramming under nitrogen stress conditions. The model was reconstructed based on the available information from KEGG, UniProt, and MaizeCyc. Transcriptomics data derived from the roots of hydroponically grown maize plants were used to incorporate regulatory constraints in the model and simulate nitrogen-non-limiting (N+) and nitrogen-deficient (N-) condition. Model-predicted flux-sum variability analysis achieved 70% accuracy compared with the experimental change of metabolite levels. In addition to predicting important metabolic reprogramming in central carbon, fatty acid, amino acid, and other secondary metabolism, maize root GSM predicted several metabolites (l-methionine, l-asparagine, l-lysine, cholesterol, and l-pipecolate) playing a regulatory role in the root biomass growth. Furthermore, this study revealed eight phosphatidylcholine and phosphatidylglycerol metabolites which, even though not coupled with biomass production, played a key role in the increased biomass production under N-deficient conditions. Overall, the omics-integrated GSM provides a promising tool to facilitate stress condition analysis for maize root and engineer better stress-tolerant maize genotypes.

Leveraging metabolic modeling to identify functional metabolic alterations associated with COVID-19 disease severity.

OBJECTIVE: Since the COVID-19 pandemic began in early 2020, SARS-CoV2 has claimed more than six million lives world-wide, with over 510 million cases to date. To reduce healthcare burden, we must investigate how to prevent non-acute disease from progressing to severe infection requiring hospitalization. METHODS: To achieve this goal, we investigated metabolic signatures of both non-acute (out-patient) and severe (requiring hospitalization) COVID-19 samples by profiling the associated plasma metabolomes of 84 COVID-19 positive University of Virginia hospital patients. We utilized supervised and unsupervised machine learning and metabolic modeling approaches to identify key metabolic drivers that are predictive of COVID-19 disease severity. Using metabolic pathway enrichment analysis, we explored potential metabolic mechanisms that link these markers to disease progression. RESULTS: Enriched metabolites associated with tryptophan in non-acute COVID-19 samples suggest mitigated innate immune system inflammatory response and immunopathology related lung damage prevention. Increased prevalence of histidine- and ketone-related metabolism in severe COVID-19 samples offers potential mechanistic insight to musculoskeletal degeneration-induced muscular weakness and host metabolism that has been hijacked by SARS-CoV2 infection to increase viral replication and invasion. CONCLUSIONS: Our findings highlight the metabolic transition from an innate immune response coupled with inflammatory pathway inhibition in non-acute infection to rampant inflammation and associated metabolic systemic dysfunction in severe COVID-19.

Integrated metabolomic and transcriptomic analyses of the synergistic effect of polymyxin-rifampicin combination against Pseudomonas aeruginosa.

BACKGROUND: Understanding the mechanism of antimicrobial action is critical for improving antibiotic therapy. For the first time, we integrated correlative metabolomics and transcriptomics of Pseudomonas aeruginosa to elucidate the mechanism of synergistic killing of polymyxin-rifampicin combination. METHODS: Liquid chromatography-mass spectrometry and RNA-seq analyses were conducted to identify the significant changes in the metabolome and transcriptome of P. aeruginosa PAO1 after exposure to polymyxin B (1 mg/L) and rifampicin (2 mg/L) alone, or in combination over 24 h. A genome-scale metabolic network was employed for integrative analysis. RESULTS: In the first 4-h treatment, polymyxin B monotherapy induced significant lipid perturbations, predominantly to fatty acids and glycerophospholipids, indicating a substantial disorganization of the bacterial outer membrane. Expression of ParRS, a two-component regulatory system involved in polymyxin resistance, was increased by polymyxin B alone. Rifampicin alone caused marginal metabolic perturbations but significantly affected gene expression at 24 h. The combination decreased the gene expression of quorum sensing regulated virulence factors at 1 h (e.g. key genes involved in phenazine biosynthesis, secretion system and biofilm formation); and increased the expression of peptidoglycan biosynthesis genes at 4 h. Notably, the combination caused substantial accumulation of nucleotides and amino acids that last at least 4 h, indicating that bacterial cells were in a state of metabolic arrest. CONCLUSION: This study underscores the substantial potential of integrative systems pharmacology to determine mechanisms of synergistic bacterial killing by antibiotic combinations, which will help optimize their use in patients.

Evolution of bidirectional costly mutualism from byproduct consumption.

Mutualisms are essential for life, yet it is unclear how they arise. A two-stage process has been proposed for the evolution of mutualisms that involve exchanges of two costly resources. First, costly provisioning by one species may be selected for if that species gains a benefit from costless byproducts generated by a second species, and cooperators get disproportionate access to byproducts. Selection could then drive the second species to provide costly resources in return. Previously, a synthetic consortium evolved the first stage of this scenario: Salmonella enterica evolved costly production of methionine in exchange for costless carbon byproducts generated by an auxotrophic Escherichia coli Growth on agar plates localized the benefits of cooperation around methionine-secreting S. enterica Here, we report that further evolution of these partners on plates led to hypercooperative E. coli that secrete the sugar galactose. Sugar secretion arose repeatedly across replicate communities and is costly to E. coli producers, but enhances the growth of S. enterica The tradeoff between individual costs and group benefits led to maintenance of both cooperative and efficient E. coli genotypes in this spatially structured environment. This study provides an experimental example of de novo, bidirectional costly mutualism evolving from byproduct consumption. The results validate the plausibility of costly cooperation emerging from initially costless exchange, a scenario widely used to explain the origin of the mutualistic species interactions that are central to life on Earth.

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.

Large-Scale Modeling Approach Reveals Functional Metabolic Shifts during Hepatic Differentiation.

Being able to explore the metabolism of broad metabolizing cells is of critical importance in many research fields. This article presents an original modeling solution combining metabolic network and omics data to identify modulated metabolic pathways and changes in metabolic functions occurring during differentiation of a human hepatic cell line (HepaRG). Our results confirm the activation of hepato-specific functionalities and newly evidence modulation of other metabolic pathways, which could not be evidenced from transcriptomic data alone. Our method takes advantage of the network structure to detect changes in metabolic pathways that do not have gene annotations and exploits flux analyses techniques to identify activated metabolic functions. Compared to the usual cell-specific metabolic network reconstruction approaches, it limits false predictions by considering several possible network configurations to represent one phenotype rather than one arbitrarily selected network. Our approach significantly enhances the comprehensive and functional assessment of cell metabolism, opening further perspectives to investigate metabolic shifts occurring within various biological contexts.

Gut microbiota dysbiosis is associated with malnutrition and reduced plasma amino acid levels: Lessons from genome-scale metabolic modeling.

Malnutrition is a severe non-communicable disease, which is prevalent in children from low-income countries. Recently, a number of metagenomics studies have illustrated associations between the altered gut microbiota and child malnutrition. However, these studies did not examine metabolic functions and interactions between individual species in the gut microbiota during health and malnutrition. Here, we applied genome-scale metabolic modeling to model the gut microbial species, which were selected from healthy and malnourished children from three countries. Our analysis showed reduced metabolite production capabilities in children from two low-income countries compared with a high-income country. Additionally, the models were also used to predict the community-level metabolic potentials of gut microbes and the patterns of pairwise interactions among species. Hereby we found that due to bacterial interactions there may be reduced production of certain amino acids in malnourished children compared with healthy children from the same communities. To gain insight into alterations in the metabolism of malnourished (stunted) children, we also performed targeted plasma metabolic profiling in the first 2 years of life of 25 healthy and 25 stunted children. Plasma metabolic profiling further revealed that stunted children had reduced plasma levels of essential amino acids compared to healthy controls. Our analyses provide a framework for future efforts towards further characterization of gut microbial metabolic capabilities and their contribution to malnutrition.

An integrated computational and experimental study uncovers FUT9 as a metabolic driver of colorectal cancer.

Metabolic alterations play an important role in cancer and yet, few metabolic cancer driver genes are known. Here we perform a combined genomic and metabolic modeling analysis searching for metabolic drivers of colorectal cancer. Our analysis predicts FUT9, which catalyzes the biosynthesis of Ley glycolipids, as a driver of advanced-stage colon cancer. Experimental testing reveals FUT9's complex dual role; while its knockdown enhances proliferation and migration in monolayers, it suppresses colon cancer cells expansion in tumorspheres and inhibits tumor development in a mouse xenograft models. These results suggest that FUT9's inhibition may attenuate tumor-initiating cells (TICs) that are known to dominate tumorspheres and early tumor growth, but promote bulk tumor cells. In agreement, we find that FUT9 silencing decreases the expression of the colorectal cancer TIC marker CD44 and the level of the OCT4 transcription factor, which is known to support cancer stemness. Beyond its current application, this work presents a novel genomic and metabolic modeling computational approach that can facilitate the systematic discovery of metabolic driver genes in other types of cancer.

Personal model-assisted identification of NAD+ and glutathione metabolism as intervention target in NAFLD.

To elucidate the molecular mechanisms underlying non-alcoholic fatty liver disease (NAFLD), we recruited 86 subjects with varying degrees of hepatic steatosis (HS). We obtained experimental data on lipoprotein fluxes and used these individual measurements as personalized constraints of a hepatocyte genome-scale metabolic model to investigate metabolic differences in liver, taking into account its interactions with other tissues. Our systems level analysis predicted an altered demand for NAD+ and glutathione (GSH) in subjects with high HS Our analysis and metabolomic measurements showed that plasma levels of glycine, serine, and associated metabolites are negatively correlated with HS, suggesting that these GSH metabolism precursors might be limiting. Quantification of the hepatic expression levels of the associated enzymes further pointed to altered de novo GSH synthesis. To assess the effect of GSH and NAD+ repletion on the development of NAFLD, we added precursors for GSH and NAD+ biosynthesis to the Western diet and demonstrated that supplementation prevents HS in mice. In a proof-of-concept human study, we found improved liver function and decreased HS after supplementation with serine (a precursor to glycine) and hereby propose a strategy for NAFLD treatment.

Chromosome 3p loss of heterozygosity is associated with a unique metabolic network in clear cell renal carcinoma.

Several common oncogenic pathways have been implicated in the emergence of renowned metabolic features in cancer, which in turn are deemed essential for cancer proliferation and survival. However, the extent to which different cancers coordinate their metabolism to meet these requirements is largely unexplored. Here we show that even in the heterogeneity of metabolic regulation a distinct signature encompassed most cancers. On the other hand, clear cell renal cell carcinoma (ccRCC) strongly deviated in terms of metabolic gene expression changes, showing widespread down-regulation. We observed a metabolic shift that associates differential regulation of enzymes in one-carbon metabolism with high tumor stage and poor clinical outcome. A significant yet limited set of metabolic genes that explained the partial divergence of ccRCC metabolism correlated with loss of von Hippel-Lindau tumor suppressor (VHL) and a potential activation of signal transducer and activator of transcription 1. Further network-dependent analyses revealed unique defects in nucleotide, one-carbon, and glycerophospholipid metabolism at the transcript and protein level, which contrasts findings in other tumors. Notably, this behavior is recapitulated by recurrent loss of heterozygosity in multiple metabolic genes adjacent to VHL. This study therefore shows how loss of heterozygosity, hallmarked by VHL deletion in ccRCC, may uniquely shape tumor metabolism.

Genome-wide gene expression and RNA half-life measurements allow predictions of regulation and metabolic behavior in Methanosarcina acetivorans.

BACKGROUND: While a few studies on the variations in mRNA expression and half-lives measured under different growth conditions have been used to predict patterns of regulation in bacterial organisms, the extent to which this information can also play a role in defining metabolic phenotypes has yet to be examined systematically. Here we present the first comprehensive study for a model methanogen. RESULTS: We use expression and half-life data for the methanogen Methanosarcina acetivorans growing on fast- and slow-growth substrates to examine the regulation of its genes. Unlike Escherichia coli where only small shifts in half-lives were observed, we found that most mRNA have significantly longer half-lives for slow growth on acetate compared to fast growth on methanol or trimethylamine. Interestingly, half-life shifts are not uniform across functional classes of enzymes, suggesting the existence of a selective stabilization mechanism for mRNAs. Using the transcriptomics data we determined whether transcription or degradation rate controls the change in transcript abundance. Degradation was found to control abundance for about half of the metabolic genes underscoring its role in regulating metabolism. Genes involved in half of the metabolic reactions were found to be differentially expressed among the substrates suggesting the existence of drastically different metabolic phenotypes that extend beyond just the methanogenesis pathways. By integrating expression data with an updated metabolic model of the organism (iST807) significant differences in pathway flux and production of metabolites were predicted for the three growth substrates. CONCLUSIONS: This study provides the first global picture of differential expression and half-lives for a class II methanogen, as well as provides the first evidence in a single organism that drastic genome-wide shifts in RNA half-lives can be modulated by growth substrate. We determined which genes in each metabolic pathway control the flux and classified them as regulated by transcription (e.g. transcription factor) or degradation (e.g. post-transcriptional modification). We found that more than half of genes in metabolism were controlled by degradation. Our results suggest that M. acetivorans employs extensive post-transcriptional regulation to optimize key metabolic steps, and more generally that degradation could play a much greater role in optimizing an organism's metabolism than previously thought.