Exploring the Interactions between Obesity and Diabetes: Implications for Understanding Metabolic Dysregulation in a Saudi Arabian Adult Population.

The rising prevalence of obesity in Saudi Arabia is a major contributor to the nation's high levels of cardiometabolic diseases such as type 2 diabetes. To assess the impact of obesity on the diabetic metabolic phenotype presented in young Saudi Arabian adults, participants (n = 289, aged 18-40 years) were recruited and stratified into four groups: healthy weight (BMI 18.5-24.99 kg/m2) with (n = 57) and without diabetes (n = 58) or overweight/obese (BMI > 24.99 kg/m2) with (n = 102) and without diabetes (n = 72). Distinct plasma metabolic phenotypes associated with high BMI and diabetes were identified using nuclear magnetic resonance spectroscopy and ultraperformance liquid chromatography mass spectrometry. Increased plasma glucose and dysregulated lipoproteins were characteristics of obesity in individuals with and without diabetes, but the obesity-associated lipoprotein phenotype was partially masked in individuals with diabetes. Although there was little difference between diabetics and nondiabetics in the global plasma LDL cholesterol and phospholipid concentration, the distribution of lipoprotein particles was altered in diabetics with a shift toward denser and more atherogenic LDL5 and LDL6 particles, which was amplified in the presence of obesity. Further investigation is warranted in larger Middle Eastern populations to explore the dysregulation of metabolism driven by interactions between obesity and diabetes in young adults.

Metabolome-wide association study on ABCA7 indicates a role of ceramide metabolism in Alzheimer's disease.

Genome-wide association studies (GWASs) have identified genetic loci associated with the risk of Alzheimer's disease (AD), but the molecular mechanisms by which they confer risk are largely unknown. We conducted a metabolome-wide association study (MWAS) of AD-associated loci from GWASs using untargeted metabolic profiling (metabolomics) by ultraperformance liquid chromatography-mass spectrometry (UPLC-MS). We identified an association of lactosylceramides (LacCer) with AD-related single-nucleotide polymorphisms (SNPs) in ABCA7 (P = 5.0 × 10-5 to 1.3 × 10-44). We showed that plasma LacCer concentrations are associated with cognitive performance and genetically modified levels of LacCer are associated with AD risk. We then showed that concentrations of sphingomyelins, ceramides, and hexosylceramides were altered in brain tissue from Abca7 knockout mice, compared with wild type (WT) (P = 0.049-1.4 × 10-5), but not in a mouse model of amyloidosis. Furthermore, activation of microglia increases intracellular concentrations of hexosylceramides in part through induction in the expression of sphingosine kinase, an enzyme with a high control coefficient for sphingolipid and ceramide synthesis. Our work suggests that the risk for AD arising from functional variations in ABCA7 is mediated at least in part through ceramides. Modulation of their metabolism or downstream signaling may offer new therapeutic opportunities for AD.

Finding Correspondence between Metabolomic Features in Untargeted Liquid Chromatography-Mass Spectrometry Metabolomics Datasets.

Integration of multiple datasets can greatly enhance bioanalytical studies, for example, by increasing power to discover and validate biomarkers. In liquid chromatography-mass spectrometry (LC-MS) metabolomics, it is especially hard to combine untargeted datasets since the majority of metabolomic features are not annotated and thus cannot be matched by chemical identity. Typically, the information available for each feature is retention time (RT), mass-to-charge ratio (m/z), and feature intensity (FI). Pairs of features from the same metabolite in separate datasets can exhibit small but significant differences, making matching very challenging. Current methods to address this issue are too simple or rely on assumptions that cannot be met in all cases. We present a method to find feature correspondence between two similar LC-MS metabolomics experiments or batches using only the features' RT, m/z, and FI. We demonstrate the method on both real and synthetic datasets, using six orthogonal validation strategies to gauge the matching quality. In our main example, 4953 features were uniquely matched, of which 585 (96.8%) of 604 manually annotated features were correct. In a second example, 2324 features could be uniquely matched, with 79 (90.8%) out of 87 annotated features correctly matched. Most of the missed annotated matches are between features that behave very differently from modeled inter-dataset shifts of RT, MZ, and FI. In a third example with simulated data with 4755 features per dataset, 99.6% of the matches were correct. Finally, the results of matching three other dataset pairs using our method are compared with a published alternative method, metabCombiner, showing the advantages of our approach. The method can be applied using M2S (Match 2 Sets), a free, open-source MATLAB toolbox, available at https://github.com/rjdossan/M2S.

Association of Untargeted Urinary Metabolomics and Lung Cancer Risk Among Never-Smoking Women in China.

Importance: Chinese women have the highest rate of lung cancer among female never-smokers in the world, and the etiology is poorly understood. Objective: To assess the association between metabolomics and lung cancer risk among never-smoking women. Design, Setting, and Participants: This nested case-control study included 275 never-smoking female patients with lung cancer and 289 never-smoking cancer-free control participants from the prospective Shanghai Women's Health Study recruited from December 28, 1996, to May 23, 2000. Validated food frequency questionnaires were used for the collection of dietary information. Metabolomic analysis was conducted from November 13, 2015, to January 6, 2016. Data analysis was conducted from January 6, 2016, to November 29, 2018. Exposures: Untargeted ultra-high-performance liquid chromatography-tandem mass spectrometry and nuclear magnetic resonance metabolomic profiles were characterized using prediagnosis urine samples. A total of 39 416 metabolites were measured. Main Outcomes and Measures: Incident lung cancer. Results: Among the 564 women, those who developed lung cancer (275 participants; median [interquartile range] age, 61.0 [52-65] years) and those who did not develop lung cancer (289 participants; median [interquartile range] age, 62.0 [53-66] years) at follow-up (median [interquartile range] follow-up, 10.9 [9.0-11.7] years) were similar in terms of their secondhand smoke exposure, history of respiratory diseases, and body mass index. A peak metabolite, identified as 5-methyl-2-furoic acid, was significantly associated with lower lung cancer risk (odds ratio, 0.57 [95% CI, 0.46-0.72]; P < .001; false discovery rate = 0.039). Furthermore, this peak was weakly correlated with self-reported dietary soy intake (ρ = 0.21; P < .001). Increasing tertiles of this metabolite were associated with lower lung cancer risk (in comparison with first tertile, odds ratio for second tertile, 0.52 [95% CI, 0.34-0.80]; and odds ratio for third tertile, 0.46 [95% CI, 0.30-0.70]), and the association was consistent across different histological subtypes and follow-up times. Additionally, metabolic pathway analysis found several systemic biological alterations that were associated with lung cancer risk, including 1-carbon metabolism, nucleotide metabolism, oxidative stress, and inflammation. Conclusions and Relevance: This prospective study of the untargeted urinary metabolome and lung cancer among never-smoking women in China provides support for the hypothesis that soy-based metabolites are associated with lower lung cancer risk in never-smoking women and suggests that biological processes linked to air pollution may be associated with higher lung cancer risk in this population.

Serum metabolic signatures of coronary and carotid atherosclerosis and subsequent cardiovascular disease.

AIMS: To characterize serum metabolic signatures associated with atherosclerosis in the coronary or carotid arteries and subsequently their association with incident cardiovascular disease (CVD). METHODS AND RESULTS: We used untargeted one-dimensional (1D) serum metabolic profiling by proton nuclear magnetic resonance spectroscopy (1H NMR) among 3867 participants from the Multi-Ethnic Study of Atherosclerosis (MESA), with replication among 3569 participants from the Rotterdam and LOLIPOP studies. Atherosclerosis was assessed by coronary artery calcium (CAC) and carotid intima-media thickness (IMT). We used multivariable linear regression to evaluate associations between NMR features and atherosclerosis accounting for multiplicity of comparisons. We then examined associations between metabolites associated with atherosclerosis and incident CVD available in MESA and Rotterdam and explored molecular networks through bioinformatics analyses. Overall, 30 1H NMR measured metabolites were associated with CAC and/or IMT, P = 1.3 × 10-14 to 1.0 × 10-6 (discovery) and P = 5.6 × 10-10 to 1.1 × 10-2 (replication). These associations were substantially attenuated after adjustment for conventional cardiovascular risk factors. Metabolites associated with atherosclerosis revealed disturbances in lipid and carbohydrate metabolism, branched chain, and aromatic amino acid metabolism, as well as oxidative stress and inflammatory pathways. Analyses of incident CVD events showed inverse associations with creatine, creatinine, and phenylalanine, and direct associations with mannose, acetaminophen-glucuronide, and lactate as well as apolipoprotein B (P < 0.05). CONCLUSION: Metabolites associated with atherosclerosis were largely consistent between the two vascular beds (coronary and carotid arteries) and predominantly tag pathways that overlap with the known cardiovascular risk factors. We present an integrated systems network that highlights a series of inter-connected pathways underlying atherosclerosis.

NMR and MS urinary metabolic phenotyping in kidney diseases is fit-for-purpose in the presence of a protease inhibitor.

Nephrotic syndrome with idiopathic membranous nephropathy as a major contributor, is characterized by proteinuria, hypoalbuminemia and oedema. Diagnosis is based on renal biopsy and the condition is treated using immunosuppressive drugs; however nephrotic syndrome treatment efficacy varies among patients. Multi-omic urine analyses can discover new markers of nephrotic syndrome that can be used to develop personalized treatments. For proteomics, a protease inhibitor (PI) is sometimes added at sample collection to conserve proteins but its impact on urine metabolic phenotyping needs to be evaluated. Urine from controls (n = 4) and idiopathic membranous nephropathy (iMN) patients (n = 6) were collected with and without PI addition and analysed using 1H NMR spectroscopy and UPLC-MS. PI-related data features were observed in the 1H NMR spectra but their removal followed by a median fold change normalisation, eliminated the PI contribution. PI-related metabolites in UPLC-MS data had limited effect on metabolic patterns specific to iMN. When using an appropriate data processing pipeline, PI-containing urine samples are appropriate for 1H NMR and MS metabolic profiling of patients with nephrotic syndrome.

A comparison of human serum and plasma metabolites using untargeted 1H NMR spectroscopy and UPLC-MS.

INTRODUCTION: Differences in the metabolite profiles between serum and plasma are incompletely understood. OBJECTIVES: To evaluate metabolic profile differences between serum and plasma and among plasma sample subtypes. METHODS: We analyzed serum, platelet rich plasma (PRP), platelet poor plasma (PPP), and platelet free plasma (PFP), collected from 8 non-fasting apparently healthy women, using untargeted standard 1D and CPMG 1H NMR and reverse phase and hydrophilic (HILIC) UPLC-MS. Differences between metabolic profiles were evaluated using validated principal component and orthogonal partial least squares discriminant analysis. RESULTS: Explorative analysis showed the main source of variation among samples was due to inter-individual differences with no grouping by sample type. After correcting for inter-individual differences, lipoproteins, lipids in VLDL/LDL, lactate, glutamine, and glucose were found to discriminate serum from plasma in NMR analyses. In UPLC-MS analyses, lysophosphatidylethanolamine (lysoPE)(18:0) and lysophosphatidic acid(20:0) were higher in serum, and phosphatidylcholines (PC)(16:1/18:2, 20:3/18:0, O-20:0/22:4), lysoPC(16:0), PE(O-18:2/20:4), sphingomyelin(18:0/22:0), and linoleic acid were lower. In plasma subtype analyses, isoleucine, leucine, valine, phenylalanine, glutamate, and pyruvate were higher among PRP samples compared with PPP and PFP by NMR while lipids in VLDL/LDL, citrate, and glutamine were lower. By UPLC-MS, PE(18:0/18:2) and PC(P-16:0/20:4) were higher in PRP compared with PFP samples. CONCLUSIONS: Correction for inter-individual variation was required to detect metabolite differences between serum and plasma. Our results suggest the potential importance of inter-individual effects and sample type on the results from serum and plasma metabolic phenotyping studies.

Improving Visualization and Interpretation of Metabolome-Wide Association Studies: An Application in a Population-Based Cohort Using Untargeted 1H NMR Metabolic Profiling.

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.

Multiplatform serum metabolic phenotyping combined with pathway mapping to identify biochemical differences in smokers.

AIM: Determining perturbed biochemical functions associated with tobacco smoking should be helpful for establishing causal relationships between exposure and adverse events. RESULTS: A multiplatform comparison of serum of smokers (n = 55) and never-smokers (n = 57) using nuclear magnetic resonance spectroscopy, UPLC-MS and statistical modeling revealed clustering of the classes, distinguished by metabolic biomarkers. The identified metabolites were subjected to metabolic pathway enrichment, modeling adverse biological events using available databases. Perturbation of metabolites involved in chronic obstructive pulmonary disease, cardiovascular diseases and cancer were identified and discussed. CONCLUSION: Combining multiplatform metabolic phenotyping with knowledge-based mapping gives mechanistic insights into disease development, which can be applied to next-generation tobacco and nicotine products for comparative risk assessment.

Workflow for Integrated Processing of Multicohort Untargeted 1H NMR Metabolomics Data in Large-Scale Metabolic Epidemiology.

Large-scale metabolomics studies involving thousands of samples present multiple challenges in data analysis, particularly when an untargeted platform is used. Studies with multiple cohorts and analysis platforms exacerbate existing problems such as peak alignment and normalization. Therefore, there is a need for robust processing pipelines that can ensure reliable data for statistical analysis. The COMBI-BIO project incorporates serum from ∼8000 individuals, in three cohorts, profiled by six assays in two phases using both 1H NMR and UPLC-MS. Here we present the COMBI-BIO NMR analysis pipeline and demonstrate its fitness for purpose using representative quality control (QC) samples. NMR spectra were first aligned and normalized. After eliminating interfering signals, outliers identified using Hotelling's T2 were removed and a cohort/phase adjustment was applied, resulting in two NMR data sets (CPMG and NOESY). Alignment of the NMR data was shown to increase the correlation-based alignment quality measure from 0.319 to 0.391 for CPMG and from 0.536 to 0.586 for NOESY, showing that the improvement was present across both large and small peaks. End-to-end quality assessment of the pipeline was achieved using Hotelling's T2 distributions. For CPMG spectra, the interquartile range decreased from 1.425 in raw QC data to 0.679 in processed spectra, while the corresponding change for NOESY spectra was from 0.795 to 0.636, indicating an improvement in precision following processing. PCA indicated that gross phase and cohort differences were no longer present. These results illustrate that the pipeline produces robust and reproducible data, successfully addressing the methodological challenges of this large multifaceted study.

Development of a Pipeline for Exploratory Metabolic Profiling of Infant Urine.

Numerous metabolic profiling pipelines have been developed to characterize the composition of human biofluids and tissues, the vast majority of these being for studies in adults. To accommodate limited sample volume and to take into account the compositional differences between adult and infant biofluids, we developed and optimized sample handling and analytical procedures for studying urine from newborns. A robust pipeline for metabolic profiling using NMR spectroscopy was established, encompassing sample collection, preparation, spectroscopic measurement, and computational analysis. Longitudinal samples were collected from five infants from birth until 14 months of age. Methods of extraction and effects of freezing and sample dilution were assessed, and urinary contaminants from breakdown of polymers in a range of diapers and cotton wool balls were identified and compared, including propylene glycol, acrylic acid, and tert-butanol. Finally, assessment of urinary profiles obtained over the first few weeks of life revealed a dramatic change in composition, with concentrations of phenols, amino acids, and betaine altering systematically over the first few months of life. Therefore, neonatal samples require more stringent standardization of experimental design, sample handling, and analysis compared to that of adult samples to accommodate the variability and limited sample volume.

Fukutin-related protein alters the deposition of laminin in the eye and brain.

Mutations in fukutin-related protein (FKRP) are responsible for a common group of muscular dystrophies ranging from adult onset limb girdle muscular dystrophies to severe congenital forms with associated structural brain involvement. The defining feature of this group of disorders is the hypoglycosylation of α-dystroglycan and its inability to effectively bind extracellular matrix ligands such as laminin α2. However, α-dystroglycan has the potential to interact with a number of laminin isoforms many of which are basement membrane/tissue specific and developmentally regulated. To further investigate this we evaluated laminin α-chain expression in the cerebral cortex and eye of our FKRP knock-down mouse (FKRP(KD)). These mice showed a marked disturbance in the deposition of laminin α-chains including α1, α2, α4, and α5, although only laminin α1- and γ1-chain mRNA expression was significantly upregulated relative to controls. Moreover, there was a diffuse pattern of laminin deposition below the pial surface which correlated with an abrupt termination of many of the radial glial cells. This along with the pial basement membrane defects, contributed to the abnormal positioning of both early- and late-born neurons. Defects in the inner limiting membrane of the eye were associated with a reduction of laminin α1 demonstrating the involvement of the α-dystroglycan:laminin α1 axis in the disease process. These observations demonstrate for the first time that a reduction in Fkrp influences the ability of tissue-specific forms of α-dystroglycan to direct the deposition of several laminin isoforms in the formation of different basement membranes.

Lamins A and C are differentially dysfunctional in autosomal dominant Emery-Dreifuss muscular dystrophy.

Mutations in the LMNA gene, which encodes nuclear lamins A and C by alternative splicing, can give rise to Emery-Dreifuss muscular dystrophy. The mechanism by which lamins A and C separately contribute to this molecular phenotype is unknown. To address this question we examined ten LMNA mutations exogenously expressed as lamins A and C in COS-7 cells. Eight of the mutations when expressed in lamin A, exhibited a range of nuclear mislocalisation patterns. However, two mutations (T150P and delQ355) almost completely relocated exogenous lamin A from the nuclear envelope to the cytoplasm, disrupted nuclear envelope reassembly following cell division and altered the protein composition of the mid-body. In contrast, exogenously expressed DsRed2-tagged mutant lamin C constructs were only inserted into the nuclear lamina if co-expressed with any EGFP-tagged lamin A construct, except with one carrying the T150P mutation. The T150P, R527P and L530P mutations reduced the ability of lamin A, but not lamin C from binding to emerin. These data identify specific functional roles for the emerin-lamin C- and emerin-lamin A- containing protein complexes and is the first report to suggest that the A-type lamin mutations may be differentially dysfunctional for the same LMNA mutation.