Symbiotic gut microorganisms (microbiome) interact closely with the mammalian host's metabolism and are important determinants of human health. Here, we decipher the complex metabolic effects of microbial manipulation, by comparing germfree mice colonized by a human baby flora (HBF) or a normal flora to conventional mice. We perform parallel microbiological profiling, metabolic profiling by (1)H nuclear magnetic resonance of liver, plasma, urine and ileal flushes, and targeted profiling of bile acids by ultra performance liquid chromatography-mass spectrometry and short-chain fatty acids in cecum by GC-FID. Top-down multivariate analysis of metabolic profiles reveals a significant association of specific metabotypes with the resident microbiome. We derive a transgenomic graph model showing that HBF flora has a remarkably simple microbiome/metabolome correlation network, impacting directly on the host's ability to metabolize lipids: HBF mice present higher ileal concentrations of tauro-conjugated bile acids, reduced plasma levels of lipoproteins but higher hepatic triglyceride content associated with depletion of glutathione. These data indicate that the microbiome modulates absorption, storage and the energy harvest from the diet at the systems level.

A major purpose of exploratory metabolic profiling is for the identification of molecular species that are statistically associated with specific biological or medical outcomes; unfortunately, the structure elucidation process of unknowns is often a major bottleneck in this process. We present here new holistic strategies that combine different statistical spectroscopic and analytical techniques to improve and simplify the process of metabolite identification. We exemplify these strategies using study data collected as part of a dietary intervention to improve health and which elicits a relatively subtle suite of changes from complex molecular profiles. We identify three new dietary biomarkers related to the consumption of peas (N-methyl nicotinic acid), apples (rhamnitol), and onions (N-acetyl-S-(1Z)-propenyl-cysteine-sulfoxide) that can be used to enhance dietary assessment and assess adherence to diet. As part of the strategy, we introduce a new probabilistic statistical spectroscopy tool, RED-STORM (Resolution EnhanceD SubseT Optimization by Reference Matching), that uses 2D J-resolved 1H NMR spectra for enhanced information recovery using the Bayesian paradigm to extract a subset of spectra with similar spectral signatures to a reference. RED-STORM provided new information for subsequent experiments (e.g., 2D-NMR spectroscopy, solid-phase extraction, liquid chromatography prefaced mass spectrometry) used to ultimately identify an unknown compound. In summary, we illustrate the benefit of acquiring J-resolved experiments alongside conventional 1D 1H NMR as part of routine metabolic profiling in large data sets and show that application of complementary statistical and analytical techniques for the identification of unknown metabolites can be used to save valuable time and resources.

The utility of low sample volume in vitro diagnostic (IVDr) proton nuclear magnetic resonance (1H NMR) spectroscopic experiments on blood plasma for information recovery from limited availability or high value samples was exemplified using plasma from patients with SARS-CoV-2 infection and normal controls. 1H NMR spectra were obtained using solvent-suppressed 1D, spin-echo (CPMG), and 2-dimensional J-resolved (JRES) spectroscopy using both 3 mm outer diameter SampleJet NMR tubes (100 μL plasma) and 5 mm SampleJet NMR tubes (300 μL plasma) under in vitro diagnostic conditions. We noted near identical diagnostic models in both standard and low volume IVDr lipoprotein analysis (measuring 112 lipoprotein parameters) with a comparison of the two tubes yielding R2 values ranging between 0.82 and 0.99 for the 40 paired lipoprotein parameters samples. Lipoprotein measurements for the 3 mm tubes were achieved without time penalty over the 5 mm tubes as defined by biomarker recovery for SARS-CoV-2. Overall, biomarker pattern recovery for the lipoproteins was extremely similar, but there were some small positive offsets in the linear equations for several variables due to small shimming artifacts, but there was minimal degradation of the biological information. For the standard untargeted 1D, CPMG, and JRES NMR experiments on the same samples, the reduced signal-to-noise was more constraining and required greater scanning times to achieve similar differential diagnostic performance (15 min per sample per experiment for 3 mm 1D and CPMG, compared to 4 min for the 5 mm tubes). We conclude that the 3 mm IVDr method is fit-for-purpose for quantitative lipoprotein measurements, allowing the preparation of smaller volumes for high value or limited volume samples that is common in clinical studies. If there are no analytical time constraints, the lower volume experiments are equally informative for untargeted profiling.

The role that metabonomics has in the evaluation of xenobiotic toxicity studies is presented here together with a brief summary of published studies. To provide a comprehensive assessment of this approach, the Consortium for Metabonomic Toxicology (COMET) has been formed between six pharmaceutical companies and Imperial College of Science, Technology and Medicine (IC), London, UK. The objective of this group is to define methodologies and to apply metabonomic data generated using (1)H NMR spectroscopy of urine and blood serum for preclinical toxicological screening of candidate drugs. This is being achieved by generating databases of results for a wide range of model toxins which serve as the raw material for computer-based expert systems for toxicity prediction. The project progress on the generation of comprehensive metabonomic databases and multivariate statistical models for prediction of toxicity, initially for liver and kidney toxicity in the rat and mouse, is reported. Additionally, both the analytical and biological variation which might arise through the use of metabonomics has been evaluated. An evaluation of intersite NMR analytical reproducibility has revealed a high degree of robustness. Second, a detailed comparison has been made of the ability of the six companies to provide consistent urine and serum samples using a study of the toxicity of hydrazine at two doses in the male rat, this study showing a high degree of consistency between samples from the various companies in terms of spectral patterns and biochemical composition. Differences between samples from the various companies were small compared to the biochemical effects of the toxin. A metabonomic model has been constructed for urine from control rats, enabling identification of outlier samples and the metabolic reasons for the deviation. Building on this success, and with the completion of studies on approximately 80 model toxins, first expert systems for prediction of liver and kidney toxicity have been generated.

The genetic regulation of metabolic phenotypes (i.e., metabotypes) in type 2 diabetes mellitus occurs through complex organ-specific cellular mechanisms and networks contributing to impaired insulin secretion and insulin resistance. Genome-wide gene expression profiling systems can dissect the genetic contributions to metabolome and transcriptome regulations. The integrative analysis of multiple gene expression traits and metabolic phenotypes (i.e., metabotypes) together with their underlying genetic regulation remains a challenge. Here, we introduce a systems genetics approach based on the topological analysis of a combined molecular network made of genes and metabolites identified through expression and metabotype quantitative trait locus mapping (i.e., eQTL and mQTL) to prioritise biological characterisation of candidate genes and traits.

We present a multivariate metabotyping approach to assess the functional recovery of nonhospitalized COVID-19 patients and the possible biochemical sequelae of "Post-Acute COVID-19 Syndrome", colloquially known as long-COVID. Blood samples were taken from patients ca. 3 months after acute COVID-19 infection with further assessment of symptoms at 6 months. Some 57% of the patients had one or more persistent symptoms including respiratory-related symptoms like cough, dyspnea, and rhinorrhea or other nonrespiratory symptoms including chronic fatigue, anosmia, myalgia, or joint pain. Plasma samples were quantitatively analyzed for lipoproteins, glycoproteins, amino acids, biogenic amines, and tryptophan pathway intermediates using Nuclear Magnetic Resonance (NMR) spectroscopy and mass spectrometry. Metabolic data for the follow-up patients (n = 27) were compared with controls (n = 41) and hospitalized severe acute respiratory syndrome SARS-CoV-2 positive patients (n = 18, with multiple time-points). Univariate and multivariate statistics revealed variable patterns of functional recovery with many patients exhibiting residual COVID-19 biomarker signatures. Several parameters were persistently perturbed, e.g., elevated taurine (p = 3.6 × 10-3 versus controls) and reduced glutamine/glutamate ratio (p = 6.95 × 10-8 versus controls), indicative of possible liver and muscle damage and a high energy demand linked to more generalized tissue repair or immune function. Some parameters showed near-complete normalization, e.g., the plasma apolipoprotein B100/A1 ratio was similar to that of healthy controls but significantly lower (p = 4.2 × 10-3) than post-acute COVID-19 patients, reflecting partial reversion of the metabolic phenotype (phenoreversion) toward the healthy metabolic state. Plasma neopterin was normalized in all follow-up patients, indicative of a reduction in the adaptive immune activity that has been previously detected in active SARS-CoV-2 infection. Other systemic inflammatory biomarkers such as GlycA and the kynurenine/tryptophan ratio remained elevated in some, but not all, patients. Correlation analysis, principal component analysis (PCA), and orthogonal-partial least-squares discriminant analysis (O-PLS-DA) showed that the follow-up patients were, as a group, metabolically distinct from controls and partially comapped with the acute-phase patients. Significant systematic metabolic differences between asymptomatic and symptomatic follow-up patients were also observed for multiple metabolites. The overall metabolic variance of the symptomatic patients was significantly greater than that of nonsymptomatic patients for multiple parameters (χ2p = 0.014). Thus, asymptomatic follow-up patients including those with post-acute COVID-19 Syndrome displayed a spectrum of multiple persistent biochemical pathophysiology, suggesting that the metabolic phenotyping approach may be deployed for multisystem functional assessment of individual post-acute COVID-19 patients.

The diffusion-ordered nuclear magnetic resonance spectroscopy (DOSY) experiment allows the calculation of diffusion coefficient values of metabolites in complex mixtures. However, this experiment has not yet been broadly used for metabolic profiling due to lack of a standardized protocol. Here we propose a pipeline for the DOSY experimental setup and data processing in metabolic phenotyping studies. Due to the complexity of biological samples, three experiments (a standard DOSY, a relaxation-edited DOSY, and a diffusion-edited DOSY) have been optimized to provide DOSY metabolic profiles with peak-picked diffusion coefficients for over 90% of signals visible in the one-dimensional 1H general biofluid profile in as little as 3 min 36 s. The developed parameter sets and tools are straightforward to implement and can facilitate the use of DOSY for metabolic profiling of human blood plasma and urine samples.

The diagnosis of joint infections is an inexact science using combinations of blood inflammatory markers and microscopy, culture, and sensitivity of synovial fluid (SF). There is potential for small molecule metabolites in infected SF to act as infection markers that could improve accuracy and speed of detection. The objective of this study was to use nuclear magnetic resonance (NMR) spectroscopy to identify small molecule differences between infected and noninfected human SF.

1H NMR spectroscopy of biofluids generates reproducible data allowing detection and quantification of small molecules in large population cohorts. Statistical models to analyze such data are now well-established, and the use of univariate metabolome wide association studies (MWAS) investigating the spectral features separately has emerged as a computationally efficient and interpretable alternative to multivariate models. The MWAS rely on the accurate estimation of a metabolome wide significance level (MWSL) to be applied to control the family wise error rate. Subsequent interpretation requires efficient visualization and formal feature annotation, which, in-turn, call for efficient prioritization of spectral variables of interest. Using human serum 1H NMR spectroscopic profiles from 3948 participants from the Multi-Ethnic Study of Atherosclerosis (MESA), we have performed a series of MWAS for serum levels of glucose. We first propose an extension of the conventional MWSL that yields stable estimates of the MWSL across the different model parameterizations and distributional features of the outcome. We propose both efficient visualization methods and a strategy based on subsampling and internal validation to prioritize the associations. Our work proposes and illustrates practical and scalable solutions to facilitate the implementation of the MWAS approach and improve interpretation in large cohort studies.

To investigate the systemic metabolic effects of SARS-CoV-2 infection, we analyzed 1H NMR spectroscopic data on human blood plasma and co-modeled with multiple plasma cytokines and chemokines (measured in parallel). Thus, 600 MHz 1H solvent-suppressed single-pulse, spin-echo, and 2D J-resolved spectra were collected on plasma recorded from SARS-CoV-2 rRT-PCR-positive patients (n = 15, with multiple sampling timepoints) and age-matched healthy controls (n = 34, confirmed rRT-PCR negative), together with patients with COVID-19/influenza-like clinical symptoms who tested SARS-CoV-2 negative (n = 35). We compared the single-pulse NMR spectral data with in vitro diagnostic research (IVDr) information on quantitative lipoprotein profiles (112 parameters) extracted from the raw 1D NMR data. All NMR methods gave highly significant discrimination of SARS-CoV-2 positive patients from controls and SARS-CoV-2 negative patients with individual NMR methods, giving different diagnostic information windows on disease-induced phenoconversion. Longitudinal trajectory analysis in selected patients indicated that metabolic recovery was incomplete in individuals without detectable virus in the recovery phase. We observed four plasma cytokine clusters that expressed complex differential statistical relationships with multiple lipoproteins and metabolites. These included the following: cluster 1, comprising MIP-1β, SDF-1α, IL-22, and IL-1α, which correlated with multiple increased LDL and VLDL subfractions; cluster 2, including IL-10 and IL-17A, which was only weakly linked to the lipoprotein profile; cluster 3, which included IL-8 and MCP-1 and were inversely correlated with multiple lipoproteins. IL-18, IL-6, and IFN-γ together with IP-10 and RANTES exhibited strong positive correlations with LDL1-4 subfractions and negative correlations with multiple HDL subfractions. Collectively, these data show a distinct pattern indicative of a multilevel cellular immune response to SARS CoV-2 infection interacting with the plasma lipoproteome giving a strong and characteristic immunometabolic phenotype of the disease. We observed that some patients in the respiratory recovery phase and testing virus-free were still metabolically highly abnormal, which indicates a new role for these technologies in assessing full systemic recovery.

To characterize serum metabolic signatures associated with atherosclerosis in the coronary or carotid arteries and subsequently their association with incident cardiovascular disease (CVD).

We have performed a metabolite quantitative trait locus (mQTL) study of the (1)H nuclear magnetic resonance spectroscopy ((1)H NMR) metabolome in humans, building on recent targeted knowledge of genetic drivers of metabolic regulation. Urine and plasma samples were collected from two cohorts of individuals of European descent, with one cohort comprised of female twins donating samples longitudinally. Sample metabolite concentrations were quantified by (1)H NMR and tested for association with genome-wide single-nucleotide polymorphisms (SNPs). Four metabolites' concentrations exhibited significant, replicable association with SNP variation (8.6×10(-11)

Interspecies variation between rats and mice has been studied for hydrazine toxicity using a novel metabonomics approach. Hydrazine hydrochloride was administered to male Sprague-Dawley rats (30 mg/kg, n = 10 and 90 mg/kg, n = 10) and male B6C3F mice (100 mg/kg, n = 8 and 250 mg/kg, n = 8) by oral gavage. In each species, the high dose was selected to produce the major histopathologic effect, hepatocellular lipid accumulation. Urine samples were collected at sequential time points up to 168 h post dose and analyzed by 1H NMR spectroscopy. The metabolites of hydrazine, namely diacetyl hydrazine and 1,4,5,6-tetrahydro-6-oxo-3-pyridazine carboxylic acid (THOPC), were detected in both the rat and mouse urine samples. Monoacetyl hydrazine was detected only in urine samples from the rat and its absence in the urine of the mouse was attributed to a higher activity of N-acetyl transferases in the mouse compared with the rat. Differential metabolic effects observed between the two species included elevated urinary beta-alanine, 3-D-hydroxybutyrate, citrulline, N-acetylcitrulline, and reduced trimethylamine-N-oxide excretion unique to the rat. Metabolic principal component (PC) trajectories highlighted the greater degree of toxic response in the rat. A data scaling method, scaled to maximum aligned and reduced trajectories (SMART) analysis, was used to remove the differences between the metabolic starting positions of the rat and mouse and varying magnitudes of effect, to facilitate comparison of the response geometries between the rat and mouse. Mice followed "biphasic" open PC trajectories, with incomplete recovery 7 days after dosing, whereas rats followed closed "hairpin" time profiles, indicating functional reversibility. The greater magnitude of metabolic effects observed in the rat was supported by the more pronounced effect on liver pathology in the rat when compared with the mouse.

The transgenomic metabolic effects of exposure to either Lactobacillus paracasei or Lactobacillus rhamnosus probiotics have been measured and mapped in humanized extended genome mice (germ-free mice colonized with human baby flora). Statistical analysis of the compartmental fluctuations in diverse metabolic compartments, including biofluids, tissue and cecal short-chain fatty acids (SCFAs) in relation to microbial population modulation generated a novel top-down systems biology view of the host response to probiotic intervention. Probiotic exposure exerted microbiome modification and resulted in altered hepatic lipid metabolism coupled with lowered plasma lipoprotein levels and apparent stimulated glycolysis. Probiotic treatments also altered a diverse range of pathways outcomes, including amino-acid metabolism, methylamines and SCFAs. The novel application of hierarchical-principal component analysis allowed visualization of multicompartmental transgenomic metabolic interactions that could also be resolved at the compartment and pathway level. These integrated system investigations demonstrate the potential of metabolic profiling as a top-down systems biology driver for investigating the mechanistic basis of probiotic action and the therapeutic surveillance of the gut microbial activity related to dietary supplementation of probiotics.

Metabolic profiling is a top-down method of analysis looking at metabolites, which are the intermediate or end products of various cellular pathways. Our primary objective was to perform a systematic review of the published literature to identify metabolites in human synovial fluid (HSF), which have been categorized by metabolic profiling techniques. A secondary objective was to identify any metabolites that may represent potential biomarkers of orthopaedic disease processes.