Missing values exist widely in mass-spectrometry (MS) based metabolomics data. Various methods have been applied for handling missing values, but the selection can significantly affect following data ...analyses. Typically, there are three types of missing values, missing not at random (MNAR), missing at random (MAR), and missing completely at random (MCAR). Our study comprehensively compared eight imputation methods (zero, half minimum (HM), mean, median, random forest (RF), singular value decomposition (SVD), k-nearest neighbors (kNN), and quantile regression imputation of left-censored data (QRILC)) for different types of missing values using four metabolomics datasets. Normalized root mean squared error (NRMSE) and NRMSE-based sum of ranks (SOR) were applied to evaluate imputation accuracy. Principal component analysis (PCA)/partial least squares (PLS)-Procrustes analysis were used to evaluate the overall sample distribution. Student's t-test followed by correlation analysis was conducted to evaluate the effects on univariate statistics. Our findings demonstrated that RF performed the best for MCAR/MAR and QRILC was the favored one for left-censored MNAR. Finally, we proposed a comprehensive strategy and developed a public-accessible web-tool for the application of missing value imputation in metabolomics ( https://metabolomics.cc.hawaii.edu/software/MetImp/ ).
Calorie restriction (CR) and fasting are common approaches to weight reduction, but the maintenance is difficult after resuming food consumption. Meanwhile, the gut microbiome associated with energy ...harvest alters dramatically in response to nutrient deprivation. Here, we reported that CR and high-fat diet (HFD) both remodeled the gut microbiota with similar microbial composition, Parabacteroides distasonis was most significantly decreased after CR or HFD. CR altered microbiota and reprogramed metabolism, resulting in a distinct serum bile acid profile characterized by depleting the proportion of non-12α-hydroxylated bile acids, ursodeoxycholic acid and lithocholic acid. Downregulation of UCP1 expression in brown adipose tissue and decreased serum GLP-1 were observed in the weight-rebound mice. Moreover, treatment with Parabacteroides distasonis or non-12α-hydroxylated bile acids ameliorated weight regain via increased thermogenesis. Our results highlighted the gut microbiota-bile acid crosstalk in rebound weight gain and Parabacteroides distasonis as a potential probiotic to prevent rapid post-CR weight gain.
Hyocholic acid (HCA) is a major bile acid (BA) species in the BA pool of pigs, a species known for its exceptional resistance to spontaneous development of diabetic phenotypes. HCA and its ...derivatives are also present in human blood and urine. We investigate whether human HCA profiles can predict the development of metabolic disorders. We find in the first cohort (n = 1107) that both obesity and diabetes are associated with lower serum concentrations of HCA species. A separate cohort study (n = 91) validates this finding and further reveals that individuals with pre-diabetes are associated with lower levels of HCA species in feces. Serum HCA levels increase in the patients after gastric bypass surgery (n = 38) and can predict the remission of diabetes two years after surgery. The results are replicated in two independent, prospective cohorts (n = 132 and n = 207), where serum HCA species are found to be strong predictors for metabolic disorders in 5 and 10 years, respectively. These findings underscore the association of HCA species with diabetes, and demonstrate the feasibility of using HCA profiles to assess the future risk of developing metabolic abnormalities.
Pu-erh tea displays cholesterol-lowering properties, but the underlying mechanism has not been elucidated. Theabrownin is one of the most active and abundant pigments in Pu-erh tea. Here, we show ...that theabrownin alters the gut microbiota in mice and humans, predominantly suppressing microbes associated with bile-salt hydrolase (BSH) activity. Theabrownin increases the levels of ileal conjugated bile acids (BAs) which, in turn, inhibit the intestinal FXR-FGF15 signaling pathway, resulting in increased hepatic production and fecal excretion of BAs, reduced hepatic cholesterol, and decreased lipogenesis. The inhibition of intestinal FXR-FGF15 signaling is accompanied by increased gene expression of enzymes in the alternative BA synthetic pathway, production of hepatic chenodeoxycholic acid, activation of hepatic FXR, and hepatic lipolysis. Our results shed light into the mechanisms behind the cholesterol- and lipid-lowering effects of Pu-erh tea, and suggest that decreased intestinal BSH microbes and/or decreased FXR-FGF15 signaling may be potential anti-hypercholesterolemia and anti-hyperlipidemia therapies.
Left-censored missing values commonly exist in targeted metabolomics datasets and can be considered as missing not at random (MNAR). Improper data processing procedures for missing values will cause ...adverse impacts on subsequent statistical analyses. However, few imputation methods have been developed and applied to the situation of MNAR in the field of metabolomics. Thus, a practical left-censored missing value imputation method is urgently needed. We developed an iterative Gibbs sampler based left-censored missing value imputation approach (GSimp). We compared GSimp with other three imputation methods on two real-world targeted metabolomics datasets and one simulation dataset using our imputation evaluation pipeline. The results show that GSimp outperforms other imputation methods in terms of imputation accuracy, observation distribution, univariate and multivariate analyses, and statistical sensitivity. Additionally, a parallel version of GSimp was developed for dealing with large scale metabolomics datasets. The R code for GSimp, evaluation pipeline, tutorial, real-world and simulated targeted metabolomics datasets are available at: https://github.com/WandeRum/GSimp.
Intestinal bacteria are known to regulate bile acid (BA) homeostasis via intestinal biotransformation of BAs and stimulation of the expression of fibroblast growth factor 19 through intestinal ...nuclear farnesoid X receptor (FXR). On the other hand, BAs directly regulate the gut microbiota with their strong antimicrobial activities. It remains unclear, however, how mammalian BAs cross-talk with gut microbiome and shape microbial composition in a dynamic and interactive way.
We quantitatively profiled small molecule metabolites derived from host-microbial co-metabolism in mice, demonstrating that BAs were the most significant factor correlated with microbial alterations among all types of endogenous metabolites. A high-fat diet (HFD) intervention resulted in a rapid and significant increase in the intestinal BA pool within 12 h, followed by an alteration in microbial composition at 24 h, providing supporting evidence that BAs are major dietary factors regulating gut microbiota. Feeding mice with BAs along with a normal diet induced an obese phenotype and obesity-associated gut microbial composition, similar to HFD-fed mice. Inhibition of hepatic BA biosynthesis under HFD conditions attenuated the HFD-induced gut microbiome alterations. Both inhibition of BAs and direct suppression of microbiota improved obese phenotypes.
Our study highlights a liver-BA-gut microbiome metabolic axis that drives significant modifications of BA and microbiota compositions capable of triggering metabolic disorders, suggesting new therapeutic strategies targeting BA metabolism for metabolic diseases.
Metabolomics data analyses rely on the use of bioinformatics tools. Many integrated multi-functional tools have been developed for untargeted metabolomics data processing and have been widely used. ...More alternative platforms are expected for both basic and advanced users. Integrated mass spectrometry-based untargeted metabolomics data mining (IP4M) software was designed and developed. The IP4M, has 62 functions categorized into 8 modules, covering all the steps of metabolomics data mining, including raw data preprocessing (alignment, peak de-convolution, peak picking, and isotope filtering), peak annotation, peak table preprocessing, basic statistical description, classification and biomarker detection, correlation analysis, cluster and sub-cluster analysis, regression analysis, ROC analysis, pathway and enrichment analysis, and sample size and power analysis. Additionally, a KEGG-derived metabolic reaction database was embedded and a series of ratio variables (product/substrate) can be generated with enlarged information on enzyme activity. A new method, GRaMM, for correlation analysis between metabolome and microbiome data was also provided. IP4M provides both a number of parameters for customized and refined analysis (for expert users), as well as 4 simplified workflows with few key parameters (for beginners who are unfamiliar with computational metabolomics). The performance of IP4M was evaluated and compared with existing computational platforms using 2 data sets derived from standards mixture and 2 data sets derived from serum samples, from GC-MS and LC-MS respectively. IP4M is powerful, modularized, customizable and easy-to-use. It is a good choice for metabolomics data processing and analysis. Free versions for Windows, MAC OS, and Linux systems are provided.
Aging is a complex physiological process associated with degenerative disorder of metabolism and immune function, which contributes to the occurrence of senile diseases. The gut microbiota affects ...systemic inflammation in aging processes probably through metabolism, but their relationship is still unclear. In this study, 16S-rRNA-sequencing technology, gas chromatography-time-of-flight mass spectrometry (GC-TOFMS)–based metabolic profiling, and immune factor analysis combined with advanced differential and association analysis were employed to investigate the correlation between the microbiome, metabolome, and immune factors in male Wistar rats across lifespan. Our findings showed significant changes in the ileum microbiome and serum metabolome compositions across aging process. A two-level strategy was applied to demonstrate that key metabolites associated with age such as 4-hydroxyproline, proline, and lysine were clustered together and positively correlated with beneficial microbes including
Bifidobacterium
,
Lactobacillus
, and
Akkermansia
. Function analysis explored association between serum metabolite class and specific gut bacteria’s metabolism pathways. Further correlation analysis on all the alteration patterns provided an interaction network of main immune factors such as IL-10, IgA, IgM, and IgG with key gut bacteria and serum metabolites. This study offers new insights into the relationship between immune factors, serum metabolome, and the gut microbiome.
Recent studies revealed strong evidence that branched-chain and aromatic amino acids (BCAAs and AAAs) are closely associated with the risk of developing type 2 diabetes in several Western countries. ...The aim of this study was to evaluate the potential role of BCAAs and AAAs in predicting the diabetes development in Chinese populations. The serum levels of valine, leucine, isoleucine, tyrosine, and phenylalanine were measured in a longitudinal and a cross sectional studies with a total of 429 Chinese participants at different stages of diabetes development, using an ultra-performance liquid chromatography triple quadruple mass spectrometry platform. The alterations of the five AAs in Chinese populations are well in accordance with previous reports. Early elevation of the five AAs and their combined score was closely associated with future development of diabetes, suggesting an important role of these metabolites as early markers of diabetes. On the other hand, the five AAs were not as good as existing clinical markers in differentiating diabetic patients from their healthy counterparts. Our findings verified the close correlation of BCAAs and AAAs with insulin resistance and future development of diabetes in Chinese populations and highlighted the predictive value of these markers for future development of diabetes.
Oral cancer, one of the six most common human cancers with an overall 5‐year survival rate of <50%, is often not diagnosed until it has reached an advanced stage. The aim of the current study is to ...explore salivary metabolomics as a disease diagnostic and stratification tool for oral cancer and leukoplakia and evaluate the potential of salivary metabolome for detection of oral squamous cell carcinoma (OSCC). Saliva metabolite profiling for a group of 37 OSCC patients, 32 oral leukoplakia (OLK) patients and 34 healthy subjects was performed using ultraperformance liquid chromatography coupled with quadrupole/time‐of‐flight mass spectrometry in conjunction with multivariate statistical analysis. The OSCC, OLK and healthy control groups demonstrate characteristic salivary metabolic signatures. A panel of five salivary metabolites including γ‐aminobutyric acid, phenylalanine, valine, n‐eicosanoic acid and lactic acid were selected using OPLS‐DA model with S‐plot. The predictive power of each of the five salivary metabolites was evaluated by receiver operating characteristic curves for OSCC. Valine, lactic acid and phenylalanine in combination yielded satisfactory accuracy (0.89, 0.97), sensitivity (86.5% and 94.6%), specificity (82.4% and 84.4%) and positive predictive value (81.6% and 87.5%) in distinguishing OSCC from the controls or OLK, respectively. The utility of salivary metabolome diagnostics for oral cancer is successfully demonstrated in this study and these results suggest that metabolomics approach complements the clinical detection of OSCC and stratifies the two types of lesions, leading to an improved disease diagnosis and prognosis.