A Protocol for Untargeted Metabolomic Analysis: From Sample Preparation to Data Processing.

Untargeted metabolomics has rapidly become a profiling method of choice in many areas of research, including mitochondrial biology. Most commonly, untargeted metabolomics is performed with liquid chromatography/mass spectrometry because it enables measurement of a relatively wide range of physiochemically diverse molecules. Specifically, to assess energy pathways that are associated with mitochondrial metabolism, hydrophilic interaction liquid chromatography (HILIC) is often applied before analysis with a high-resolution accurate mass instrument. The workflow produces large, complex data files that are impractical to analyze manually. Here, we present a protocol to perform untargeted metabolomics on biofluids such as plasma, urine, and cerebral spinal fluid with a HILIC separation and an Orbitrap mass spectrometer. Our protocol describes each step of the analysis in detail, from preparation of solvents for chromatography to selecting parameters during data processing.

ENL links histone acetylation to oncogenic gene expression in acute myeloid leukaemia.

Cancer cells are characterized by aberrant epigenetic landscapes and often exploit chromatin machinery to activate oncogenic gene expression programs. Recognition of modified histones by 'reader' proteins constitutes a key mechanism underlying these processes; therefore, targeting such pathways holds clinical promise, as exemplified by the development of bromodomain and extra-terminal (BET) inhibitors. We recently identified the YEATS domain as an acetyl-lysine-binding module, but its functional importance in human cancer remains unknown. Here we show that the YEATS domain-containing protein ENL, but not its paralogue AF9, is required for disease maintenance in acute myeloid leukaemia. CRISPR-Cas9-mediated depletion of ENL led to anti-leukaemic effects, including increased terminal myeloid differentiation and suppression of leukaemia growth in vitro and in vivo. Biochemical and crystal structural studies and chromatin-immunoprecipitation followed by sequencing analyses revealed that ENL binds to acetylated histone H3, and co-localizes with H3K27ac and H3K9ac on the promoters of actively transcribed genes that are essential for leukaemia. Disrupting the interaction between the YEATS domain and histone acetylation via structure-based mutagenesis reduced the recruitment of RNA polymerase II to ENL-target genes, leading to the suppression of oncogenic gene expression programs. Notably, disrupting the functionality of ENL further sensitized leukaemia cells to BET inhibitors. Together, our data identify ENL as a histone acetylation reader that regulates oncogenic transcriptional programs in acute myeloid leukaemia, and suggest that displacement of ENL from chromatin may be a promising epigenetic therapy, alone or in combination with BET inhibitors, for aggressive leukaemia.

Distinct metabolomic signatures are associated with longevity in humans.

Alterations in metabolism influence lifespan in experimental models, but data in humans are lacking. Here we use liquid chromatography/mass spectrometry to quantify 217 plasma metabolites and examine their relation to longevity in a large cohort of men and women followed for up to 20 years. We find that, higher concentrations of the citric acid cycle intermediate, isocitrate, and the bile acid, taurocholate, are associated with lower odds of longevity, defined as attaining 80 years of age. Higher concentrations of isocitrate, but not taurocholate, are also associated with worse cardiovascular health at baseline, as well as risk of future cardiovascular disease and death. None of the metabolites identified are associated with cancer risk. Our findings suggest that some, but not all, metabolic pathways related to human longevity are linked to the risk of common causes of death.

5-hydroxymethylcytosine is not present in appreciable quantities in Arabidopsis DNA.

5-Hydroxymethylcytosine (5-hmC) is an intermediate in active demethylation in metazoans, as well as a potentially stable epigenetic mark. Previous reports investigating 5-hydroxymethylcytosine in plants have reached conflicting conclusions. We systematically investigated whether 5-hmC is present in plant DNA using a range of methods. Using the model organism Arabidopsis thaliana, in addition to other plant species, we assayed the amount or distribution of 5-hydroxymethylcytosine by thin-layer chromatography, immunoprecipitation-chip, ELISA, enzymatic radiolabeling, and mass spectrometry. The failure to observe 5-hydroxymethylcytosine by thin-layer chromatography established an upper bound for the possible fraction of the nucleotide in plant DNA. Antibody-based methods suggested that there were low levels of 5-hmC in plant DNA, but these experiments were potentially confounded by cross-reactivity with the abundant base 5-methylcytosine. Enzymatic radiolabeling and mass spectrometry, the most sensitive methods for detection that we used, failed to detect 5-hydroxymethylcytosine in A. thaliana genomic DNA isolated from a number of different tissue types and genetic backgrounds. Taken together, our results led us to conclude that 5-hmC is not present in biologically relevant quantities within plant genomic DNA.

ß-Aminoisobutyric acid induces browning of white fat and hepatic ß-oxidation and is inversely correlated with cardiometabolic risk factors.

The transcriptional coactivator peroxisome proliferator-activated receptor-gamma coactivator-1α (PGC-1α) regulates metabolic genes in skeletal muscle and contributes to the response of muscle to exercise. Muscle PGC-1α transgenic expression and exercise both increase the expression of thermogenic genes within white adipose. How the PGC-1α-mediated response to exercise in muscle conveys signals to other tissues remains incompletely defined. We employed a metabolomic approach to examine metabolites secreted from myocytes with forced expression of PGC-1α, and identified ß-aminoisobutyric acid (BAIBA) as a small molecule myokine. BAIBA increases the expression of brown adipocyte-specific genes in white adipocytes and ß-oxidation in hepatocytes both in vitro and in vivo through a PPARα-mediated mechanism, induces a brown adipose-like phenotype in human pluripotent stem cells, and improves glucose homeostasis in mice. In humans, plasma BAIBA concentrations are increased with exercise and inversely associated with metabolic risk factors. BAIBA may thus contribute to exercise-induced protection from metabolic diseases.

2-Aminoadipic acid is a biomarker for diabetes risk.

Improvements in metabolite-profiling techniques are providing increased breadth of coverage of the human metabolome and may highlight biomarkers and pathways in common diseases such as diabetes. Using a metabolomics platform that analyzes intermediary organic acids, purines, pyrimidines, and other compounds, we performed a nested case-control study of 188 individuals who developed diabetes and 188 propensity-matched controls from 2,422 normoglycemic participants followed for 12 years in the Framingham Heart Study. The metabolite 2-aminoadipic acid (2-AAA) was most strongly associated with the risk of developing diabetes. Individuals with 2-AAA concentrations in the top quartile had greater than a 4-fold risk of developing diabetes. Levels of 2-AAA were not well correlated with other metabolite biomarkers of diabetes, such as branched chain amino acids and aromatic amino acids, suggesting they report on a distinct pathophysiological pathway. In experimental studies, administration of 2-AAA lowered fasting plasma glucose levels in mice fed both standard chow and high-fat diets. Further, 2-AAA treatment enhanced insulin secretion from a pancreatic ß cell line as well as murine and human islets. These data highlight a metabolite not previously associated with diabetes risk that is increased up to 12 years before the onset of overt disease. Our findings suggest that 2-AAA is a marker of diabetes risk and a potential modulator of glucose homeostasis in humans.

A genome-wide association study of the human metabolome in a community-based cohort.

Because metabolites are hypothesized to play key roles as markers and effectors of cardiometabolic diseases, recent studies have sought to annotate the genetic determinants of circulating metabolite levels. We report a genome-wide association study (GWAS) of 217 plasma metabolites, including >100 not measured in prior GWAS, in 2076 participants of the Framingham Heart Study (FHS). For the majority of analytes, we find that estimated heritability explains >20% of interindividual variation, and that variation attributable to heritable factors is greater than that attributable to clinical factors. Further, we identify 31 genetic loci associated with plasma metabolites, including 23 that have not previously been reported. Importantly, we include GWAS results for all surveyed metabolites and demonstrate how this information highlights a role for AGXT2 in cholesterol ester and triacylglycerol metabolism. Thus, our study outlines the relative contributions of inherited and clinical factors on the plasma metabolome and provides a resource for metabolism research.

Reproducibility of metabolomic profiles among men and women in 2 large cohort studies.

BACKGROUND: Rigorous studies are necessary to demonstrate suitability of metabolomics platforms to profile metabolites in archived plasma within epidemiologic studies of human disease, for which attenuation of effect estimates due to measurement error is a key concern. METHODS: Using a liquid chromatography-tandem mass spectrometry platform, we quantified 257 metabolites from archived plasma to evaluate metabolite interassay reproducibility, reproducibility with delayed processing, and within-person reproducibility over time. Interassay reproducibility was assessed with CVs from 60 duplicate plasma samples donated by participants in the Nurses' Health Study and Health Professionals Follow-up Study, and 20 QC pool plasma replicates. Metabolite reproducibility over a 24- to 48-h processing delay (n = 48 samples) and within-person reproducibility over 1-2 years (n = 80 samples) were assessed using Spearman and intraclass correlation coefficients (ICCs). RESULTS: CVs were <20% for 92% of metabolites and generally were similar by plasma anticoagulant type (heparin or EDTA) and fasting time. Approximately 75% of metabolites were reproducible over delays in processing of blood samples (Spearman correlation or ICC ≥ 0.75, comparing immediate and 24-h delayed processing). Carbohydrates and purine/pyrimidine derivatives were most adversely affected by the processing delay. Ninety percent of metabolites were reproducible over 1-2 years within individuals (Spearman correlation or ICC ≥ 0.4). CONCLUSIONS: For potential use in epidemiologic studies, the majority of plasma metabolites had low CVs and were reproducible over a 24-h processing delay and within individuals over 1-2 years. Certain metabolites, such as carbohydrates and purine/pyrimidine derivatives, may be challenging to evaluate if samples have delayed processing.

A combined epidemiologic and metabolomic approach improves CKD prediction.

Metabolomic approaches have begun to catalog the metabolic disturbances that accompany CKD, but whether metabolite alterations can predict future CKD is unknown. We performed liquid chromatography/mass spectrometry-based metabolite profiling on plasma from 1434 participants in the Framingham Heart Study (FHS) who did not have CKD at baseline. During the following 8 years, 123 individuals developed CKD, defined by an estimated GFR of <60 ml/min per 1.73 m(2). Numerous metabolites were associated with incident CKD, including 16 that achieved the Bonferroni-adjusted significance threshold of P≤0.00023. To explore how the human kidney modulates these metabolites, we profiled arterial and renal venous plasma from nine individuals. Nine metabolites that predicted CKD in the FHS cohort decreased more than creatinine across the renal circulation, suggesting that they may reflect non-GFR-dependent functions, such as renal metabolism and secretion. Urine isotope dilution studies identified citrulline and choline as markers of renal metabolism and kynurenic acid as a marker of renal secretion. In turn, these analytes remained associated with incident CKD in the FHS cohort, even after adjustment for eGFR, age, sex, diabetes, hypertension, and proteinuria at baseline. Addition of a multimarker metabolite panel to clinical variables significantly increased the c-statistic (0.77-0.83, P<0.0001); net reclassification improvement was 0.78 (95% confidence interval, 0.60 to 0.95; P<0.0001). Thus, the addition of metabolite profiling to clinical data may significantly improve the ability to predict whether an individual will develop CKD by identifying predictors of renal risk that are independent of estimated GFR.

SIRT4 has tumor-suppressive activity and regulates the cellular metabolic response to DNA damage by inhibiting mitochondrial glutamine metabolism.

DNA damage elicits a cellular signaling response that initiates cell cycle arrest and DNA repair. Here, we find that DNA damage triggers a critical block in glutamine metabolism, which is required for proper DNA damage responses. This block requires the mitochondrial SIRT4, which is induced by numerous genotoxic agents and represses the metabolism of glutamine into tricarboxylic acid cycle. SIRT4 loss leads to both increased glutamine-dependent proliferation and stress-induced genomic instability, resulting in tumorigenic phenotypes. Moreover, SIRT4 knockout mice spontaneously develop lung tumors. Our data uncover SIRT4 as an important component of the DNA damage response pathway that orchestrates a metabolic block in glutamine metabolism, cell cycle arrest, and tumor suppression.

mTOR complex 1 plays critical roles in hematopoiesis and Pten-loss-evoked leukemogenesis.

The mechanistic target of rapamycin (mTOR) pathway serves as a key sensor of cellular-energetic state and functions to maintain tissue homeostasis. Hyperactivation of the mTOR pathway impairs hematopoietic stem cell (HSC) function and is associated with leukemogenesis. However, the roles of the unique mTOR complexes (mTORCs) in hematopoiesis and leukemogenesis have not been adequately elucidated. We deleted the mTORC1 component, regulatory-associated protein of mTOR (Raptor), in mouse HSCs and its loss causes a nonlethal phenotype characterized by pancytopenia, splenomegaly, and the accumulation of monocytoid cells. Furthermore, Raptor is required for HSC regeneration, and plays largely nonredundant roles with rapamycin-insensitive companion of mTOR (Rictor) in these processes. Ablation of Raptor also significantly extends survival of mice in models of leukemogenesis evoked by Pten deficiency. These data delineate critical roles for mTORC1 in hematopoietic function and leukemogenesis and inform clinical strategies based on chronic mTORC1 inhibition.

Metabolite profiling identifies a key role for glycine in rapid cancer cell proliferation.

Metabolic reprogramming has been proposed to be a hallmark of cancer, yet a systematic characterization of the metabolic pathways active in transformed cells is currently lacking. Using mass spectrometry, we measured the consumption and release (CORE) profiles of 219 metabolites from media across the NCI-60 cancer cell lines, and integrated these data with a preexisting atlas of gene expression. This analysis identified glycine consumption and expression of the mitochondrial glycine biosynthetic pathway as strongly correlated with rates of proliferation across cancer cells. Antagonizing glycine uptake and its mitochondrial biosynthesis preferentially impaired rapidly proliferating cells. Moreover, higher expression of this pathway was associated with greater mortality in breast cancer patients. Increased reliance on glycine may represent a metabolic vulnerability for selectively targeting rapid cancer cell proliferation.

Metabolite profiling identifies pathways associated with metabolic risk in humans.

BACKGROUND: Although metabolic risk factors are known to cluster in individuals who are prone to developing diabetes mellitus and cardiovascular disease, the underlying biological mechanisms remain poorly understood. METHODS AND RESULTS: To identify pathways associated with cardiometabolic risk, we used liquid chromatography/mass spectrometry to determine the plasma concentrations of 45 distinct metabolites and to examine their relation to cardiometabolic risk in the Framingham Heart Study (FHS; n=1015) and the Malmö Diet and Cancer Study (MDC; n=746). We then interrogated significant findings in experimental models of cardiovascular and metabolic disease. We observed that metabolic risk factors (obesity, insulin resistance, high blood pressure, and dyslipidemia) were associated with multiple metabolites, including branched-chain amino acids, other hydrophobic amino acids, tryptophan breakdown products, and nucleotide metabolites. We observed strong associations of insulin resistance traits with glutamine (standardized regression coefficients, -0.04 to -0.22 per 1-SD change in log-glutamine; P<0.001), glutamate (0.05 to 0.14; P<0.001), and the glutamine-to-glutamate ratio (-0.05 to -0.20; P<0.001) in the discovery sample (FHS); similar associations were observed in the replication sample (MDC). High glutamine-to-glutamate ratio was associated with lower risk of incident diabetes mellitus in FHS (odds ratio, 0.79; adjusted P=0.03) but not in MDC. In experimental models, administration of glutamine in mice led to both increased glucose tolerance (P=0.01) and decreased blood pressure (P<0.05). CONCLUSIONS: Biochemical profiling identified circulating metabolites not previously associated with metabolic traits. Experimentally interrogating one of these pathways demonstrated that excess glutamine relative to glutamate, resulting from exogenous administration, is associated with reduced metabolic risk in mice.

Targeted metabolomics.

The metabolome is the terminal downstream product of the genome and consists of the total complement of all the low-molecular-weight molecules (metabolites) in a cell, tissue, or organism. Metabolomics aims to measure a wide breadth of small molecules in the context of physiological stimuli or disease states. Metabolomics methodologies fall into two distinct groups: untargeted metabolomics, an intended comprehensive analysis of all the measurable analytes in a sample including chemical unknowns, and targeted metabolomics, the measurement of defined groups of chemically characterized and biochemically annotated metabolites. The methodologies considered in this unit focus on the processes of conducting targeted metabolomics experiments, and the advantages of this general approach are highlighted herein. This unit outlines procedures for extracting nitrogenous metabolites (including amino acids), lipids, and intermediary metabolites (including TCA cycle oxoacids) from blood plasma. Specifically, protocols are described for analyzing these metabolites using targeted metabolomics experiments based on liquid chromatography-mass spectrometry.

Metabolic consequences of mitochondrial coenzyme A deficiency in patients with PANK2 mutations.

Pantothenate kinase-associated neurodegeneration (PKAN) is a rare, inborn error of metabolism characterized by iron accumulation in the basal ganglia and by the presence of dystonia, dysarthria, and retinal degeneration. Mutations in pantothenate kinase 2 (PANK2), the rate-limiting enzyme in mitochondrial coenzyme A biosynthesis, represent the most common genetic cause of this disorder. How mutations in this core metabolic enzyme give rise to such a broad clinical spectrum of pathology remains a mystery. To systematically explore its pathogenesis, we performed global metabolic profiling on plasma from a cohort of 14 genetically defined patients and 18 controls. Notably, lactate is elevated in PKAN patients, suggesting dysfunctional mitochondrial metabolism. As predicted, but never previously reported, pantothenate levels are higher in patients with premature stop mutations in PANK2. Global metabolic profiling and follow-up studies in patient-derived fibroblasts also reveal defects in bile acid conjugation and lipid metabolism, pathways that require coenzyme A. These findings raise a novel therapeutic hypothesis, namely, that dietary fats and bile acid supplements may hold potential as disease-modifying interventions. Our study illustrates the value of metabolic profiling as a tool for systematically exploring the biochemical basis of inherited metabolic diseases.

A hypoxia-induced positive feedback loop promotes hypoxia-inducible factor 1alpha stability through miR-210 suppression of glycerol-3-phosphate dehydrogenase 1-like.

Oxygen-dependent regulation of the transcription factor HIF-1α relies on a family of prolyl hydroxylases (PHDs) that hydroxylate hypoxia-inducible factor 1α (HIF-1α) protein at two prolines during normal oxygen conditions, resulting in degradation by the proteasome. During low-oxygen conditions, these prolines are no longer hydroxylated and HIF-1α degradation is blocked. Hypoxia-induced miRNA-210 (miR-210) is a direct transcriptional target of HIF-1α, but its complete role and targets during hypoxia are not well understood. Here, we identify the enzyme glycerol-3-phosphate dehydrogenase 1-like (GPD1L) as a novel regulator of HIF-1α stability and a direct target of miR-210. Expression of miR-210 results in stabilization of HIF-1α due to decreased levels of GPD1L resulting in an increase in HIF-1α target genes. Altering GPD1L levels by overexpression or knockdown results in a decrease or increase in HIF-1α stability, respectively. GPD1L-mediated decreases in HIF-1α stability can be reversed by pharmacological inhibition of the proteasome or PHD activity. When rescued from degradation by proteasome inhibition, elevated amounts of GPD1L cause hyperhydroxylation of HIF-1α, suggesting increases in PHD activity. Importantly, expression of GPD1L attenuates the hypoxic response, preventing complete HIF-1α induction. We propose a model in which hypoxia-induced miR-210 represses GPD1L, contributing to suppression of PHD activity, and increases of HIF-1α protein levels.

Activation of a metabolic gene regulatory network downstream of mTOR complex 1.

Aberrant activation of the mammalian target of rapamycin complex 1 (mTORC1) is a common molecular event in a variety of pathological settings, including genetic tumor syndromes, cancer, and obesity. However, the cell-intrinsic consequences of mTORC1 activation remain poorly defined. Through a combination of unbiased genomic, metabolomic, and bioinformatic approaches, we demonstrate that mTORC1 activation is sufficient to stimulate specific metabolic pathways, including glycolysis, the oxidative arm of the pentose phosphate pathway, and de novo lipid biosynthesis. This is achieved through the activation of a transcriptional program affecting metabolic gene targets of hypoxia-inducible factor (HIF1alpha) and sterol regulatory element-binding protein (SREBP1 and SREBP2). We find that SREBP1 and 2 promote proliferation downstream of mTORC1, and the activation of these transcription factors is mediated by S6K1. Therefore, in addition to promoting protein synthesis, mTORC1 activates specific bioenergetic and anabolic cellular processes that are likely to contribute to human physiology and disease.

A plasma signature of human mitochondrial disease revealed through metabolic profiling of spent media from cultured muscle cells.

Mutations in either the mitochondrial or nuclear genomes can give rise to respiratory chain disease (RCD), a large class of devastating metabolic disorders. Their clinical management is challenging, in part because we lack facile and accurate biomarkers to aid in diagnosis and in the monitoring of disease progression. Here we introduce a sequential strategy that combines biochemical analysis of spent media from cell culture with analysis of patient plasma to identify disease biomarkers. First, we applied global metabolic profiling to spotlight 32 metabolites whose uptake or secretion kinetics were altered by chemical inhibition of the respiratory chain in cultured muscle . These metabolites span a wide range of pathways and include lactate and alanine, which are used clinically as biomarkers of RCD. We next measured the cell culture-defined metabolites in human plasma to discover that creatine is reproducibly elevated in two independent cohorts of RCD patients, exceeding lactate and alanine in magnitude of elevation and statistical significance. In cell culture extracellular creatine was inversely related to the intracellular phosphocreatine:creatine ratio suggesting that the elevation of plasma creatine in RCD patients signals a low energetic state of tissues using the phosphocreatine shuttle. Our study identifies plasma creatine as a potential biomarker of human mitochondrial dysfunction that could be clinically useful. More generally, we illustrate how spent media from cellular models of disease may provide a window into the biochemical derangements in human plasma, an approach that could, in principle, be extended to a range of complex diseases.