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.

Concordance-Based Batch Effect Correction for Large-Scale Metabolomics.

For a large-scale metabolomics study, sample collection, preparation, and analysis may last several days, months, or even (intermittently) over years. This may lead to apparent batch effects in the acquired metabolomics data due to variability in instrument status, environmental conditions, or experimental operators. Batch effects may confound the true biological relationships among metabolites and thus obscure real metabolic changes. At present, most of the commonly used batch effect correction (BEC) methods are based on quality control (QC) samples, which require sufficient and stable QC samples. However, the quality of the QC samples may deteriorate if the experiment lasts for a long time. Alternatively, isotope-labeled internal standards have been used, but they generally do not provide good coverage of the metabolome. On the other hand, BEC can also be conducted through a data-driven method, in which no QC sample is needed. Here, we propose a novel data-driven BEC method, namely, CordBat, to achieve concordance between each batch of samples. In the proposed CordBat method, a reference batch is first selected from all batches of data, and the remaining batches are referred to as "other batches." The reference batch serves as the baseline for the batch adjustment by providing a coordinate of correlation between metabolites. Next, a Gaussian graphical model is built on the combined dataset of reference and other batches, and finally, BEC is achieved by optimizing the correction coefficients in the other batches so that the correlation between metabolites of each batch and their combinations are in concordance with that of the reference batch. Three real-world metabolomics datasets are used to evaluate the performance of CordBat by comparing it with five commonly used BEC methods. The present experimental results showed the effectiveness of CordBat in batch effect removal and the concordance of correlation between metabolites after BEC. CordBat was found to be comparable to the QC-based methods and achieved better performance in the preservation of biological effects. The proposed CordBat method may serve as an alternative BEC method for large-scale metabolomics that lack proper QC samples.

iMSEA: A Novel Metabolite Set Enrichment Analysis Strategy to Decipher Drug Interactions.

Drug combinations are commonly used to treat various diseases to achieve synergistic therapeutic effects or to alleviate drug resistance. Nevertheless, some drug combinations might lead to adverse effects, and thus, it is crucial to explore the mechanisms of drug interactions before clinical treatment. Generally, drug interactions have been studied using nonclinical pharmacokinetics, toxicology, and pharmacology. Here, we propose a complementary strategy based on metabolomics, which we call interaction metabolite set enrichment analysis, or iMSEA, to decipher drug interactions. First, a digraph-based heterogeneous network model was constructed to model the biological metabolic network based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Second, treatment-specific influences on all detected metabolites were calculated and propagated across the whole network model. Third, pathway activity was defined and enriched to quantify the influence of each treatment on the predefined functional metabolite sets, i.e., metabolic pathways. Finally, drug interactions were identified by comparing the pathway activity enriched by the drug combination treatments and the single drug treatments. A data set consisting of hepatocellular carcinoma (HCC) cells that were treated with oxaliplatin (OXA) and/or vitamin C (VC) was used to illustrate the effectiveness of the iMSEA strategy for evaluation of drug interactions. Performance evaluation using synthetic noise data was also performed to evaluate sensitivities and parameter settings for the iMSEA strategy. The iMSEA strategy highlighted synergistic effects of combined OXA and VC treatments including the alterations in the glycerophospholipid metabolism pathway and glycine, serine, and threonine metabolism pathway. This work provides an alternative method to reveal the mechanisms of drug combinations from the viewpoint of metabolomics.

Multimodal Image Fusion Offers Better Spatial Resolution for Mass Spectrometry Imaging.

High-resolution reconstruction has attracted increasing research interest in mass spectrometry imaging (MSI), but it remains a challenging ill-posed problem. In the present study, we proposed a deep learning model to fuse multimodal images to enhance the spatial resolution of MSI data, namely, DeepFERE. Hematoxylin and eosin (H&E) stain microscopy imaging was used to pose constraints in the process of high-resolution reconstruction to alleviate the ill-posedness. A novel model architecture was designed to achieve multi-task optimization by incorporating multi-modal image registration and fusion in a mutually reinforced framework. Experimental results demonstrated that the proposed DeepFERE model is able to produce high-resolution reconstruction images with rich chemical information and a detailed structure on both visual inspection and quantitative evaluation. In addition, our method was found to be able to improve the delimitation of the boundary between cancerous and para-cancerous regions in the MSI image. Furthermore, the reconstruction of low-resolution spatial transcriptomics data demonstrated that the developed DeepFERE model may find wider applications in biomedical fields.