Identifying Significant Metabolic Pathways Using Multi-Block Partial Least-Squares Analysis

Metabolomic analysis of hydroxycinnamic acid inhibition on Saccharomyces cerevisiae

Ferulic acid (FA) and p-coumaric acid (p-CA) are hydroxycinnamic acid inhibitors that are mainly produced during the pretreatment of lignocellulose. To date, the inhibitory mechanism of hydroxycinnamic acid compounds on Saccharomyces cerevisiae has not been fully elucidated. In this study, liquid chromatography-mass spectrometry (LC-MS) and scanning electron microscopy (SEM) were used to investigate the changes in S. cerevisiae cells treated with FA and p-CA. In this experiment, the control group was denoted as group CK, the FA-treated group was denoted as group F, and the p-CA-treated group was denoted as group P. One hundred different metabolites in group F and group CK and 92 different metabolites in group P and group CK were selected and introduced to metaboanalyst, respectively. A total of 38 metabolic pathways were enriched in S. cerevisiae under FA stress, and 27 metabolic pathways were enriched in S. cerevisiae under p-CA stress as identified through Kyoto Encyclopaedia of Genes and Genomes (KEGG) analysis. The differential metabolites involved included S-adenosine methionine, L-arginine, and cysteine, which were significantly downregulated, and acetyl-CoA, L-glutamic acid, and L-threonine, which were significantly upregulated. Analysis of differential metabolic pathways showed that the differentially expressed metabolites were mainly related to amino acid metabolism, nucleotide metabolism, fatty acid degradation, and the tricarboxylic acid cycle (TCA). Under the stress of FA and p-CA, the metabolism of some amino acids was blocked, which disturbed the redox balance in the cells and destroyed the synthesis of most proteins, which was the main reason for the inhibition of yeast cell growth. This study provided a strong scientific reference to improve the durability of S. cerevisiae against hydroxycinnamic acid inhibitors. KEY POINTS: • Morphological changes of S. cerevisiae cells under inhibitors stress were observed. • Changes of the metabolites in S. cerevisiae cells were explored by metabolomics. • One of the inhibitory effects on yeast is due to changes in the metabolic network.

Differential Correlations Informed Metabolite Set Enrichment Analysis to Decipher Metabolic Heterogeneity of Disease

Metabolic pathways are regarded as functional and basic components of the biological system. In metabolomics, metabolite set enrichment analysis (MSEA) is often used to identify the altered metabolic pathways (metabolite sets) associated with phenotypes of interest (POI), e.g., disease. However, in most studies, MSEA suffers from the limitation of low metabolite coverage. Random walk (RW)-based algorithms can be used to propagate the perturbation of detected metabolites to the undetected metabolites through a metabolite network model prior to MSEA. Nevertheless, most of the existing RW-based algorithms run on a general metabolite network constructed based on public databases, such as KEGG, without taking into consideration the potential influence of POI on the metabolite network, which may reduce the phenotypic specificities of the MSEA results. To solve this problem, a novel pathway analysis strategy, namely, differential correlation-informed MSEA (dci-MSEA), is proposed in this paper. Statistically, differential correlations between metabolites are used to evaluate the influence of POI on the metabolite network, so that a phenotype-specific metabolite network is constructed for RW-based propagation. The experimental results show that dci-MSEA outperforms the conventional RW-based MSEA in identifying the altered metabolic pathways associated with colorectal cancer. In addition, by incorporating the individual-specific metabolite network, the dci-MSEA strategy is easily extended to disease heterogeneity analysis. Here, dci-MSEA was used to decipher the heterogeneity of colorectal cancer. The present results highlight the clustering of colorectal cancer samples with their cluster-specific selection of differential pathways and demonstrate the feasibility of dci-MSEA in heterogeneity analysis. Taken together, the proposed dci-MSEA may provide insights into disease mechanisms and determination of disease heterogeneity.

Exploration of Blood Metabolite Signatures of Colorectal Cancer and Polyposis through Integrated Statistical and Network Analysis

Colorectal cancer (CRC), one of the most prevalent and deadly cancers worldwide, generally evolves from adenomatous polyps. The understanding of the molecular mechanisms underlying this pathological evolution is crucial for diagnostic and prognostic purposes. Integrative systems biology approaches offer an optimal point of view to analyze CRC and patients with polyposis. The present study analyzed the association networks constructed from a publicly available array of 113 serum metabolites measured on a cohort of 234 subjects from three groups (66 CRC patients, 76 patients with polyposis, and 92 healthy controls), which concentrations were obtained via targeted liquid chromatography-tandem mass spectrometry. In terms of architecture, topology, and connectivity, the metabolite-metabolite association network of CRC patients appears to be completely different with respect to patients with polyposis and healthy controls. The most relevant nodes in the CRC network are those related to energy metabolism. Interestingly, phenylalanine, tyrosine, and tryptophan metabolism are found to be involved in both CRC and polyposis. Our results demonstrate that the characterization of metabolite-metabolite association networks is a promising and powerful tool to investigate molecular aspects of CRC.

Comparison of freeze-thaw and sonication cycle-based methods for extracting AMR-associated metabolites from Staphylococcus aureus

Emerging antimicrobial resistance (AMR) among Gram-positive pathogens, specifically in Staphylococcus aureus (S. aureus), is becoming a leading public health concern demanding effective therapeutics. Metabolite modulation can improve the efficacy of existing antibiotics and facilitate the development of effective therapeutics. However, it remained unexplored for drug-resistant S. aureus (gentamicin and methicillin-resistant), primarily due to the dearth of optimal metabolite extraction protocols including a protocol for AMR-associated metabolites. Therefore, in this investigation, we have compared the performance of the two most widely used methods, i.e., freeze-thaw cycle (FTC) and sonication cycle (SC), alone and in combination (FTC + SC), and identified the optimal method for this purpose. A total of 116, 119, and 99 metabolites were identified using the FTC, SC, and FTC + SC methods, respectively, leading to the identification of 163 metabolites cumulatively. Out of 163, 69 metabolites were found to be associated with AMR in published literature consisting of the highest number of metabolites identified by FTC (57) followed by SC (54) and FTC + SC (40). Thus, the performances of FTC and SC methods were comparable with no additional benefits of combining both. Moreover, each method showed biasness toward specific metabolite(s) or class of metabolites, suggesting that the choice of metabolite extraction method shall be decided based on the metabolites of interest in the investigation.

Metabolic flux between organs measured by arteriovenous metabolite gradients

Mammalian organs convert dietary nutrients into circulating metabolites and share them to maintain whole-body metabolic homeostasis. While the concentrations of circulating metabolites have been frequently measured in a variety of pathophysiological conditions, the exchange flux of circulating metabolites between organs is not easily measurable due to technical difficulties. Isotope tracing is useful for measuring such fluxes for a metabolite of interest, but the shuffling of isotopic atoms between metabolites requires mathematical modeling. Arteriovenous metabolite gradient measurements can complement isotope tracing to infer organ-specific net fluxes of many metabolites simultaneously. Here, we review the historical development of arteriovenous measurements and discuss their advantages and limitations with key example studies that have revealed metabolite exchange flux between organs in diverse pathophysiological contexts.

Sparse PLS-Based Method for Overlapping Metabolite Set Enrichment Analysis

Metabolite set enrichment analysis (MSEA) has gained increasing research interest for identification of perturbed metabolic pathways in metabolomics. The method incorporates predefined metabolic pathways information in the analysis where metabolite sets are typically assumed to be mutually exclusive to each other. However, metabolic pathways are known to contain common metabolites and intermediates. This situation, along with limitations in metabolite detection or coverage leads to overlapping, incomplete metabolite sets in pathway analysis. For overlapping metabolite sets, MSEA tends to result in high false positives due to improper weights allocated to the overlapping metabolites. Here, we proposed an extended partial least squares (PLS) model with a new sparse scheme for overlapping metabolite set enrichment analysis, named overlapping group PLS (ogPLS) analysis. The weight vector of the ogPLS model was decomposed into pathway-specific subvectors, and then a group lasso penalty was imposed on these subvectors to achieve a proper weight allocation for the overlapping metabolites. Two strategies were adopted in the proposed ogPLS model to identify the perturbed metabolic pathways. The first strategy involves debiasing regularization, which was used to reduce inequalities amongst the predefined metabolic pathways. The second strategy is stable selection, which was used to rank pathways while avoiding the nuisance problems of model parameter optimization. Both simulated and real-world metabolomic datasets were used to evaluate the proposed method and compare with two other MSEA methods including Global-test and the multiblock PLS (MB-PLS)-based pathway importance in projection (PIP) methods. Using a simulated dataset with known perturbed pathways, the average true discovery rate for the ogPLS method was found to be higher than the Global-test and the MB-PLS-based PIP methods. Analysis with a real-world metabolomics dataset also indicated that the developed method was less prone to select pathways with highly overlapped detected metabolite sets. Compared with the two other methods, the proposed method features higher accuracy, lower false-positive rate, and is more robust when applied to overlapping metabolite set analysis. The developed ogPLS method may serve as an alternative MSEA method to facilitate biological interpretation of metabolomics data for overlapping metabolite sets.