Integrative Molecular Structure Elucidation and Construction of an Extended Metabolic Pathway Associated with an Ancient Innate Immune Response in COVID-19 Patients.

We present compelling evidence for the existence of an extended innate viperin-dependent pathway, which provides crucial evidence for an adaptive response to viral agents, such as SARS-CoV-2. We show the in vivo biosynthesis of a family of novel endogenous cytosine metabolites with potential antiviral activities. Two-dimensional nuclear magnetic resonance (NMR) spectroscopy revealed a characteristic spin-system motif, indicating the presence of an extended panel of urinary metabolites during the acute viral replication phase. Mass spectrometry additionally enabled the characterization and quantification of the most abundant serum metabolites, showing the potential diagnostic value of the compounds for viral infections. In total, we unveiled ten nucleoside (cytosine- and uracil-based) analogue structures, eight of which were previously unknown in humans allowing us to propose a new extended viperin pathway for the innate production of antiviral compounds. The molecular structures of the nucleoside analogues and their correlation with an array of serum cytokines, including IFN-α2, IFN-γ, and IL-10, suggest an association with the viperin enzyme contributing to an ancient endogenous innate immune defense mechanism against viral infection.

Balancing the Equation: A Natural History of Trimethylamine and Trimethylamine-N-oxide.

Trimethylamine (TMA) and its N-oxide (TMAO) are ubiquitous in prokaryote and eukaryote organisms as well as in the environment, reflecting their fundamental importance in evolutionary biology, and their diverse biochemical functions. Both metabolites have multiple biological roles including cell-signaling. Much attention has focused on the significance of serum and urinary TMAO in cardiovascular disease risk, yet this is only one of the many facets of a deeper TMA-TMAO partnership that reflects the significance of these metabolites in multiple biological processes spanning animals, plants, bacteria, and fungi. We report on analytical methods for measuring TMA and TMAO and attempt to critically synthesize and map the global functions of TMA and TMAO in a systems biology framework.

Exploration of Human Serum Lipoprotein Supramolecular Phospholipids Using Statistical Heterospectroscopy in n-Dimensions (SHY-n): Identification of Potential Cardiovascular Risk Biomarkers Related to SARS-CoV-2 Infection.

SARS-CoV-2 infection causes a significant reduction in lipoprotein-bound serum phospholipids give rise to supramolecular phospholipid composite (SPC) signals observed in diffusion and relaxation edited 1H NMR spectra. To characterize the chemical structural components and compartmental location of SPC and to understand further its possible diagnostic properties, we applied a Statistical HeterospectroscopY in n-dimensions (SHY-n) approach. This involved statistically linking a series of orthogonal measurements made on the same samples, using independent analytical techniques and instruments, to identify the major individual phospholipid components giving rise to the SPC signals. Thus, an integrated model for SARS-CoV-2 positive and control adults is presented that relates three identified diagnostic subregions of the SPC signal envelope (SPC1, SPC2, and SPC3) generated using diffusion and relaxation edited (DIRE) NMR spectroscopy to lipoprotein and lipid measurements obtained by in vitro diagnostic NMR spectroscopy and ultrahigh-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS). The SPC signals were then correlated sequentially with (a) total phospholipids in lipoprotein subfractions; (b) apolipoproteins B100, A1, and A2 in different lipoproteins and subcompartments; and (c) MS-measured total serum phosphatidylcholines present in the NMR detection range (i.e., PCs: 16.0,18.2; 18.0,18.1; 18.2,18.2; 16.0,18.1; 16.0,20.4; 18.0,18.2; 18.1,18.2), lysophosphatidylcholines (LPCs: 16.0 and 18.2), and sphingomyelin (SM 22.1). The SPC3/SPC2 ratio correlated strongly (r = 0.86) with the apolipoprotein B100/A1 ratio, a well-established marker of cardiovascular disease risk that is markedly elevated during acute SARS-CoV-2 infection. These data indicate the considerable potential of using a serum SPC measurement as a metric of cardiovascular risk based on a single NMR experiment. This is of specific interest in relation to understanding the potential for increased cardiovascular risk in COVID-19 patients and risk persistence in post-acute COVID-19 syndrome (PACS).

Strategy for improved characterization of human metabolic phenotypes using a COmbined Multi-block Principal components Analysis with Statistical Spectroscopy (COMPASS).

MOTIVATION: Large-scale population omics data can provide insight into associations between gene-environment interactions and disease. However, existing dimension reduction modelling techniques are often inefficient for extracting detailed information from these complex datasets. RESULTS: Here, we present an interactive software pipeline for exploratory analyses of population-based nuclear magnetic resonance spectral data using a COmbined Multi-block Principal components Analysis with Statistical Spectroscopy (COMPASS) within the R-library hastaLaVista framework. Principal component analysis models are generated for a sequential series of spectral regions (blocks) to provide more granular detail defining sub-populations within the dataset. Molecular identification of key differentiating signals is subsequently achieved by implementing Statistical TOtal Correlation SpectroscopY on the full spectral data to define feature patterns. Finally, the distributions of cross-correlation of the reference patterns across the spectral dataset are used to provide population statistics for identifying underlying features arising from drug intake, latent diseases and diet. The COMPASS method thus provides an efficient semi-automated approach for screening population datasets. AVAILABILITY AND IMPLEMENTATION: Source code is available at https://github.com/cheminfo/COMPASS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Low Volume in Vitro Diagnostic Proton NMR Spectroscopy of Human Blood Plasma for Lipoprotein and Metabolite Analysis: Application to SARS-CoV-2 Biomarkers.

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.

NMR Spectroscopic Windows on the Systemic Effects of SARS-CoV-2 Infection on Plasma Lipoproteins and Metabolites in Relation to Circulating Cytokines.

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.

Integrative Modeling of Plasma Metabolic and Lipoprotein Biomarkers of SARS-CoV-2 Infection in Spanish and Australian COVID-19 Patient Cohorts.

Quantitative plasma lipoprotein and metabolite profiles were measured on an autonomous community of the Basque Country (Spain) cohort consisting of hospitalized COVID-19 patients (n = 72) and a matched control group (n = 75) and a Western Australian (WA) cohort consisting of (n = 17) SARS-CoV-2 positives and (n = 20) healthy controls using 600 MHz 1H nuclear magnetic resonance (NMR) spectroscopy. Spanish samples were measured in two laboratories using one-dimensional (1D) solvent-suppressed and T2-filtered methods with in vitro diagnostic quantification of lipoproteins and metabolites. SARS-CoV-2 positive patients and healthy controls from both populations were modeled and cross-projected to estimate the biological similarities and validate biomarkers. Using the top 15 most discriminatory variables enabled construction of a cross-predictive model with 100% sensitivity and specificity (within populations) and 100% sensitivity and 82% specificity (between populations). Minor differences were observed between the control metabolic variables in the two cohorts, but the lipoproteins were virtually indistinguishable. We observed highly significant infection-related reductions in high-density lipoprotein (HDL) subfraction 4 phospholipids, apolipoproteins A1 and A2,that have previously been associated with negative regulation of blood coagulation and fibrinolysis. The Spanish and Australian diagnostic SARS-CoV-2 biomarkers were mathematically and biologically equivalent, demonstrating that NMR-based technologies are suitable for the study of the comparative pathology of COVID-19 via plasma phenotyping.

Incomplete Systemic Recovery and Metabolic Phenoreversion in Post-Acute-Phase Nonhospitalized COVID-19 Patients: Implications for Assessment of Post-Acute COVID-19 Syndrome.

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.

Diffusion and Relaxation Edited Proton NMR Spectroscopy of Plasma Reveals a High-Fidelity Supramolecular Biomarker Signature of SARS-CoV-2 Infection.

We have applied nuclear magnetic resonance spectroscopy based plasma phenotyping to reveal diagnostic molecular signatures of SARS-CoV-2 infection via combined diffusional and relaxation editing (DIRE). We compared plasma from healthy age-matched controls (n = 26) with SARS-CoV-2 negative non-hospitalized respiratory patients and hospitalized respiratory patients (n = 23 and 11 respectively) with SARS-CoV-2 rRT-PCR positive respiratory patients (n = 17, with longitudinal sampling time-points). DIRE data were modelled using principal component analysis and orthogonal projections to latent structures discriminant analysis (O-PLS-DA), with statistical cross-validation indices indicating excellent model generalization for the classification of SARS-CoV-2 positivity for all comparator groups (area under the receiver operator characteristic curve = 1). DIRE spectra show biomarker signal combinations conferred by differential concentrations of metabolites with selected molecular mobility properties. These comprise the following: (a) composite N-acetyl signals from α-1-acid glycoprotein and other glycoproteins (designated GlycA and GlycB) that were elevated in SARS-CoV-2 positive patients [p = 2.52 × 10-10 (GlycA) and 1.25 × 10-9 (GlycB) vs controls], (b) two diagnostic supramolecular phospholipid composite signals that were identified (SPC-A and SPC-B) from the -+N-(CH3)3 choline headgroups of lysophosphatidylcholines carried on plasma glycoproteins and from phospholipids in high-density lipoprotein subfractions (SPC-A) together with a phospholipid component of low-density lipoprotein (SPC-B). The integrals of the summed SPC signals (SPCtotal) were reduced in SARS-CoV-2 positive patients relative to both controls (p = 1.40 × 10-7) and SARS-CoV-2 negative patients (p = 4.52 × 10-8) but were not significantly different between controls and SARS-CoV-2 negative patients. The identity of the SPC signal components was determined using one and two dimensional diffusional, relaxation, and statistical spectroscopic experiments. The SPCtotal/GlycA ratios were also significantly different for control versus SARS-CoV-2 positive patients (p = 1.23 × 10-10) and for SARS-CoV-2 negatives versus positives (p = 1.60 × 10-9). Thus, plasma SPCtotal and SPCtotal/GlycA are proposed as sensitive molecular markers for SARS-CoV-2 positivity that could effectively augment current COVID-19 diagnostics and may have value in functional assessment of the disease recovery process in patients with long-term symptoms.

Integrative Modeling of Quantitative Plasma Lipoprotein, Metabolic, and Amino Acid Data Reveals a Multiorgan Pathological Signature of SARS-CoV-2 Infection.

The metabolic effects of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on human blood plasma were characterized using multiplatform metabolic phenotyping with nuclear magnetic resonance (NMR) spectroscopy and liquid chromatography-mass spectrometry (LC-MS). Quantitative measurements of lipoprotein subfractions, α-1-acid glycoprotein, glucose, and biogenic amines were made on samples from symptomatic coronavirus disease 19 (COVID-19) patients who had tested positive for the SARS-CoV-2 virus (n = 17) and from age- and gender-matched controls (n = 25). Data were analyzed using an orthogonal-projections to latent structures (OPLS) method and used to construct an exceptionally strong (AUROC = 1) hybrid NMR-MS model that enabled detailed metabolic discrimination between the groups and their biochemical relationships. Key discriminant metabolites included markers of inflammation including elevated α-1-acid glycoprotein and an increased kynurenine/tryptophan ratio. There was also an abnormal lipoprotein, glucose, and amino acid signature consistent with diabetes and coronary artery disease (low total and HDL Apolipoprotein A1, low HDL triglycerides, high LDL and VLDL triglycerides), plus multiple highly significant amino acid markers of liver dysfunction (including the elevated glutamine/glutamate and Fischer's ratios) that present themselves as part of a distinct SARS-CoV-2 infection pattern. A multivariate training-test set model was validated using independent samples from additional SARS-CoV-2 positive patients and controls. The predictive model showed a sensitivity of 100% for SARS-CoV-2 positivity. The breadth of the disturbed pathways indicates a systemic signature of SARS-CoV-2 positivity that includes elements of liver dysfunction, dyslipidemia, diabetes, and coronary heart disease risk that are consistent with recent reports that COVID-19 is a systemic disease affecting multiple organs and systems. Metabolights study reference: MTBLS2014.

Quantitative In-Vitro Diagnostic NMR Spectroscopy for Lipoprotein and Metabolite Measurements in Plasma and Serum: Recommendations for Analytical Artifact Minimization with Special Reference to COVID-19/SARS-CoV-2 Samples.

Quantitative nuclear magnetic resonance (NMR) spectroscopy of blood plasma is widely used to investigate perturbed metabolic processes in human diseases. The reliability of biochemical data derived from these measurements is dependent on the quality of the sample collection and exact preparation and analysis protocols. Here, we describe systematically, the impact of variations in sample collection and preparation on information recovery from quantitative proton (1H) NMR spectroscopy of human blood plasma and serum. The effects of variation of blood collection tube sizes and preservatives, successive freeze-thaw cycles, sample storage at -80 °C, and short-term storage at 4 and 20 °C on the quantitative lipoprotein and metabolite patterns were investigated. Storage of plasma samples at 4 °C for up to 48 h, freezing at -80 °C and blood sample collection tube choice have few and minor effects on quantitative lipoprotein profiles, and even storage at 4 °C for up to 168 h caused little information loss. In contrast, the impact of heat-treatment (56 °C for 30 min), which has been used for inactivation of SARS-CoV-2 and other viruses, that may be required prior to analytical measurements in low level biosecurity facilities induced marked changes in both lipoprotein and low molecular weight metabolite profiles. It was conclusively demonstrated that this heat inactivation procedure degrades lipoproteins and changes metabolic information in complex ways. Plasma from control individuals and SARS-CoV-2 infected patients are differentially altered resulting in the creation of artifactual pseudo-biomarkers and destruction of real biomarkers to the extent that data from heat-treated samples are largely uninterpretable. We also present several simple blood sample handling recommendations for optimal NMR-based biomarker discovery investigations in SARS CoV-2 studies and general clinical biomarker research.

Automatic Spectroscopic Data Categorization by Clustering Analysis (ASCLAN): A Data-Driven Approach for Distinguishing Discriminatory Metabolites for Phenotypic Subclasses.

We propose a novel data-driven approach aiming to reliably distinguish discriminatory metabolites from nondiscriminatory metabolites for a given spectroscopic data set containing two biological phenotypic subclasses. The automatic spectroscopic data categorization by clustering analysis (ASCLAN) algorithm aims to categorize spectral variables within a data set into three clusters corresponding to noise, nondiscriminatory and discriminatory metabolites regions. This is achieved by clustering each spectral variable based on the r(2) value representing the loading weight of each spectral variable as extracted from a orthogonal partial least-squares discriminant (OPLS-DA) model of the data set. The variables are ranked according to r(2) values and a series of principal component analysis (PCA) models are then built for subsets of these spectral data corresponding to ranges of r(2) values. The Q(2)X value for each PCA model is extracted. K-means clustering is then applied to the Q(2)X values to generate two clusters based on minimum Euclidean distance criterion. The cluster consisting of lower Q(2)X values is deemed devoid of metabolic information (noise), while the cluster consists of higher Q(2)X values is then further subclustered into two groups based on the r(2) values. We considered the cluster with high Q(2)X but low r(2) values as nondiscriminatory, while the cluster with high Q(2)X and r(2) values as discriminatory variables. The boundaries between these three clusters of spectral variables, on the basis of the r(2) values were considered as the cut off values for defining the noise, nondiscriminatory and discriminatory variables. We evaluated the ASCLAN algorithm using six simulated (1)H NMR spectroscopic data sets representing small, medium and large data sets (N = 50, 500, and 1000 samples per group, respectively), each with a reduced and full resolution set of variables (0.005 and 0.0005 ppm, respectively). ASCLAN correctly identified all discriminatory metabolites and showed zero false positive (100% specificity and positive predictive value) irrespective of the spectral resolution or the sample size in all six simulated data sets. This error rate was found to be superior to existing methods for ascertaining feature significance: univariate t test by Bonferroni correction (up to 10% false positive rate), Benjamini-Hochberg correction (up to 35% false positive rate) and metabolome wide significance level (MWSL, up to 0.4% false positive rate), as well as by various OPLS-DA parameters: variable importance to projection, (up to 15% false positive rate), loading coefficients (up to 35% false positive rate), and regression coefficients (up to 39% false positive rate). The application of ASCLAN was further exemplified using a widely investigated renal toxin, mercury II chloride (HgCl2) in rat model. ASCLAN successfully identified many of the known metabolites related to renal toxicity such as increased excretion of urinary creatinine, and different amino acids. The ASCLAN algorithm provides a framework for reliably differentiating discriminatory metabolites from nondiscriminatory metabolites in a biological data set without the need to set an arbitrary cut off value as applied to some of the conventional methods. This offers significant advantages over existing methods and the possibility for automation of high-throughput screening in "omics" data.

Development and validation of an ultra-performance liquid chromatography quadrupole time of flight mass spectrometry method for rapid quantification of free amino acids in human urine.

An ultra-performance liquid chromatography quadrupole time of flight mass spectrometry (UPLC-qTOF-MS) method using hydrophilic interaction liquid chromatography was developed and validated for simultaneous quantification of 18 free amino acids in urine with a total acquisition time including the column re-equilibration of less than 18 min per sample. This method involves simple sample preparation steps which consisted of 15 times dilution with acetonitrile to give a final composition of 25 % aqueous and 75 % acetonitrile without the need of any derivatization. The dynamic range for our calibration curve is approximately two orders of magnitude (120-fold from the lowest calibration curve point) with good linearity (r (2) ≥ 0.995 for all amino acids). Good separation of all amino acids as well as good intra- and inter-day accuracy (<15 %) and precision (<15 %) were observed using three quality control samples at a concentration of low, medium and high range of the calibration curve. The limits of detection (LOD) and lower limit of quantification of our method were ranging from approximately 1-300 nM and 0.01-0.5 µM, respectively. The stability of amino acids in the prepared urine samples was found to be stable for 72 h at 4 °C, after one freeze thaw cycle and for up to 4 weeks at -80 °C. We have applied this method to quantify the content of 18 free amino acids in 646 urine samples from a dietary intervention study. We were able to quantify all 18 free amino acids in these urine samples, if they were present at a level above the LOD. We found our method to be reproducible (accuracy and precision were typically <10 % for QCL, QCM and QCH) and the relatively high sample throughput nature of this method potentially makes it a suitable alternative for the analysis of urine samples in clinical setting.

Statistical HOmogeneous Cluster SpectroscopY (SHOCSY): an optimized statistical approach for clustering of ¹H NMR spectral data to reduce interference and enhance robust biomarkers selection.

We propose a novel statistical approach to improve the reliability of (1)H NMR spectral analysis in complex metabolic studies. The Statistical HOmogeneous Cluster SpectroscopY (SHOCSY) algorithm aims to reduce the variation within biological classes by selecting subsets of homogeneous (1)H NMR spectra that contain specific spectroscopic metabolic signatures related to each biological class in a study. In SHOCSY, we used a clustering method to categorize the whole data set into a number of clusters of samples with each cluster showing a similar spectral feature and hence biochemical composition, and we then used an enrichment test to identify the associations between the clusters and the biological classes in the data set. We evaluated the performance of the SHOCSY algorithm using a simulated (1)H NMR data set to emulate renal tubule toxicity and further exemplified this method with a (1)H NMR spectroscopic study of hydrazine-induced liver toxicity study in rats. The SHOCSY algorithm improved the predictive ability of the orthogonal partial least-squares discriminatory analysis (OPLS-DA) model through the use of "truly" representative samples in each biological class (i.e., homogeneous subsets). This method ensures that the analyses are no longer confounded by idiosyncratic responders and thus improves the reliability of biomarker extraction. SHOCSY is a useful tool for removing irrelevant variation that interfere with the interpretation and predictive ability of models and has widespread applicability to other spectroscopic data, as well as other "omics" type of data.

Human metabolic phenotype diversity and its association with diet and blood pressure.

Metabolic phenotypes are the products of interactions among a variety of factors-dietary, other lifestyle/environmental, gut microbial and genetic. We use a large-scale exploratory analytical approach to investigate metabolic phenotype variation across and within four human populations, based on 1H NMR spectroscopy. Metabolites discriminating across populations are then linked to data for individuals on blood pressure, a major risk factor for coronary heart disease and stroke (leading causes of mortality worldwide). We analyse spectra from two 24-hour urine specimens for each of 4,630 participants from the INTERMAP epidemiological study, involving 17 population samples aged 40-59 in China, Japan, UK and USA. We show that urinary metabolite excretion patterns for East Asian and western population samples, with contrasting diets, diet-related major risk factors, and coronary heart disease/stroke rates, are significantly differentiated (P < 10(-16)), as are Chinese/Japanese metabolic phenotypes, and subgroups with differences in dietary vegetable/animal protein and blood pressure. Among discriminatory metabolites, we quantify four and show association (P < 0.05 to P < 0.0001) of mean 24-hour urinary formate excretion with blood pressure in multiple regression analyses for individuals. Mean 24-hour urinary excretion of alanine (direct) and hippurate (inverse), reflecting diet and gut microbial activities, are also associated with blood pressure of individuals. Metabolic phenotyping applied to high-quality epidemiological data offers the potential to develop an area of aetiopathogenetic knowledge involving discovery of novel biomarkers related to cardiovascular disease risk.

Unravelling inulin molecules in food sources using a matrix-assisted laser desorption/ionization magnetic resonance mass spectrometry (MALDI-MRMS) pipeline.

Inulin, a polysaccharide characterized by a ß-2,1 fructosyl-fructose structure terminating in a glucosyl moiety, is naturally present in plant roots and tubers. Current methods provide average degrees of polymerization (DP) but lack information on the distribution and absolute concentration of each DP. To address this limitation, a reproducible (CV < 10 %) high throughput (<2 min/sample) MALDI-MRMS approach capable of characterizing and quantifying inulin molecules in plants using matched-matrix consisting of α-cyano-4-hydroxycinnamic acid butylamine salt (CHCA-BA), chicory inulin-12C and inulin-13C was developed. The method identified variation in chain lengths and concentration of inulin across various plant species. Globe artichoke hearts, yacón and elephant garlic yielded similar concentrations at 15.6 g/100 g dry weight (DW), 16.8 g/100 g DW and 17.7 g/100 g DW, respectively, for DP range between 9 and 22. In contrast, Jerusalem artichoke demonstrated the highest concentration (53.4 g/100 g DW) within the same DP ranges. Jerusalem artichoke (DPs 9-32) and globe artichoke (DPs 9-36) showed similar DP distributions, while yacón and elephant garlic displayed the narrowest and broadest DP ranges (DPs 9-19 and DPs 9-45, respectively). Additionally, qualitative measurement for all inulin across all plant samples was feasible using the peak intensities normalized to Inulin-13C, and showed that the ratio of yacón, elephant garlic and Jerusalem was approximately one, two and three times that of globe artichoke. This MALDI-MRMS approach provides comprehensive insights into the structure of inulin molecules, opening avenues for in-depth investigations into how DP and concentration of inulin influence gut health and the modulation of noncommunicable diseases, as well as shedding light on refining cultivation practices to elevate the beneficial health properties associated with specific DPs.

Plasma lipoprotein subclass variation in middle-aged and older adults: Sex-stratified distributions and associations with health status and cardiometabolic risk factors.

BACKGROUND: Circulating lipids and lipoproteins mediate cardiovascular risk, however routine plasma lipid biochemistry provides limited information on pro-atherogenic remnant particles. OBJECTIVE: We analysed plasma lipoprotein subclasses including very low-density and intermediate-density lipoprotein (VLDL and IDL); and assessed their associations with health and cardiometabolic risk. METHODS: From 1,976 community-dwelling adults aged 45-67 years, 114/1071 women (10.6%) and 153/905 men (16.9%) were categorised as very healthy. Fasting plasma lipoprotein profiles comprising 112 parameters were measured using 1H nuclear magnetic resonance (NMR) spectroscopy, and associations with health status and cardiometabolic risk factors examined. RESULTS: HDL cholesterol was higher, and IDL and VLDL cholesterol and triglycerides lower, in very healthy women compared to other women, and women compared to men. IDL and VLDL cholesterol and triglyceride were lower in very healthy men compared to other men. HDL cholesterol and apolipoprotein (apo) A-I were inversely, and IDL and VLDL cholesterol, apoB-100, and apoB-100/apoA-I ratio directly associated with body mass index (BMI) in women and men. In women, LDL, IDL and VLDL cholesterol increased with age. Women with diabetes and cardiovascular disease had higher cholesterol, triglycerides, phospholipids and free cholesterol across IDL and VLDL fractions, with similar trends for men with diabetes. CONCLUSION: Lipoprotein subclasses and density fractions, and their lipid and apolipoprotein constituents, are differentially distributed by sex, health status and BMI. Very healthy women and men are distinguished by favorable lipoprotein profiles, particularly lower concentrations of VLDL and IDL, providing reference intervals for comparison with general populations and adults with cardiometabolic risk factors.

Targeted lipidomics coupled with machine learning for authenticating the provenance of chicken eggs.

Free-range eggs are ethically desirable but as with all high-value commercial products, the establishment of provenance can be problematic. Here, we compared a simple one-step isopropanol method to a two-step methyl-tert-butyl ether method for extracting lipid species in chicken egg yolks before liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. The isopropanol method extracted 937 lipid species from 20 major lipid subclasses with high reproducibility (CV < 30 %). Machine learning techniques could differentiate conventional cage, barn, and free-range eggs using an external test dataset with an accuracy of 0.94, 0.82, and 0.82, respectively. Lipid species that differentiated cage eggs were predominantly phosphocholines and phosphoethanolamines whilst the free-range egg lipidomes were dominated by acylglycerides with up to three fatty acids. The lipid profiles were found to be characteristic of the cage, barns, and free-range eggs. The lipidomic analysis together with the statistical modeling approach thus provides an efficient tool for verifying the provenance of conventional chicken eggs.

Demographic, Health and Lifestyle Factors Associated with the Metabolome in Older Women.

Demographic and clinical factors influence the metabolome. The discovery and validation of disease biomarkers are often challenged by potential confounding effects from such factors. To address this challenge, we investigated the magnitude of the correlation between serum and urine metabolites and demographic and clinical parameters in a well-characterized observational cohort of 444 post-menopausal women participating in the Women's Health Initiative (WHI). Using LC-MS and lipidomics, we measured 157 aqueous metabolites and 756 lipid species across 13 lipid classes in serum, along with 195 metabolites detected by GC-MS and NMR in urine and evaluated their correlations with 29 potential disease risk factors, including demographic, dietary and lifestyle factors, and medication use. After controlling for multiple testing (FDR < 0.01), we found that log-transformed metabolites were mainly associated with age, BMI, alcohol intake, race, sample storage time (urine only), and dietary supplement use. Statistically significant correlations were in the absolute range of 0.2-0.6, with the majority falling below 0.4. Incorporation of important potential confounding factors in metabolite and disease association analyses may lead to improved statistical power as well as reduced false discovery rates in a variety of data analysis settings.