Genetic Influences on Metabolite Levels: A Comparison across Metabolomic Platforms

Idil Yet; Cristina Menni; So-Youn Shin; Massimo Mangino; Nicole Soranzo; Jerzy Adamski; Karsten Suhre; Tim D. Spector; Gabi Kastenmüller; Jordana T. Bell

doi:10.1371/journal.pone.0153672

Abstract

Metabolomic profiling is a powerful approach to characterize human metabolism and help understand common disease risk. Although multiple high-throughput technologies have been developed to assay the human metabolome, no technique is capable of capturing the entire human metabolism. Large-scale metabolomics data are being generated in multiple cohorts, but the datasets are typically profiled using different metabolomics platforms. Here, we compared analyses across two of the most frequently used metabolomic platforms, Biocrates and Metabolon, with the aim of assessing how complimentary metabolite profiles are across platforms. We profiled serum samples from 1,001 twins using both targeted (Biocrates, n = 160 metabolites) and non-targeted (Metabolon, n = 488 metabolites) mass spectrometry platforms. We compared metabolite distributions and performed genome-wide association analyses to identify shared genetic influences on metabolites across platforms. Comparison of 43 metabolites named for the same compound on both platforms indicated strong positive correlations, with few exceptions. Genome-wide association scans with high-throughput metabolic profiles were performed for each dataset and identified genetic variants at 7 loci associated with 16 unique metabolites on both platforms. The 16 metabolites showed consistent genetic associations and appear to be robustly measured across platforms. These included both metabolites named for the same compound across platforms as well as unique metabolites, of which 2 (nonanoylcarnitine (C9) [Biocrates]/Unknown metabolite X-13431 [Metabolon] and PC aa C28:1 [Biocrates]/1-stearoylglycerol [Metabolon]) are likely to represent the same or related biochemical entities. The results demonstrate the complementary nature of both platforms, and can be informative for future studies of comparative and integrative metabolomics analyses in samples profiled on different platforms.

Citation: Yet I, Menni C, Shin S-Y, Mangino M, Soranzo N, Adamski J, et al. (2016) Genetic Influences on Metabolite Levels: A Comparison across Metabolomic Platforms. PLoS ONE 11(4): e0153672. https://doi.org/10.1371/journal.pone.0153672

Editor: Nicholas John Timpson, University of Bristol, UNITED KINGDOM

Received: September 22, 2015; Accepted: April 1, 2016; Published: April 13, 2016

Copyright: © 2016 Yet et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Data can be shared to third parties through application only via our Data Access Committee (DAC) - the Twins Resource Executive Committee. The committee needs to ensure subject anonymity and that the ethical constraints by which subjects consent to give their data is adhered to. This decision to provide supported access to our data is driven through our ethics and Terms of Reference of our governance body.

Funding: The study was funded in part by the Wellcome Trust and the European Community’s Seventh Framework Programme (FP7/2007-2013). TwinsUK also received support from the Wellcome Trust, the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. TS is an NIHR senior Investigator and is holder of an ERC Advanced Principal Investigator award. SNP Genotyping was performed by The Wellcome Trust Sanger Institute and National Eye Institute via NIH/CIDR.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Metabolomics aims to provide a comprehensive characterization of human metabolic pathways by high throughput profiling of the small molecules present in biological samples. Various metabolomics platforms have been established to date, based on mass spectrometry (MS) or nuclear magnetic resonance (NMR) technology. However, individual platforms can differ in a number of features, including the set of metabolites quantified, the precision of quantification, and its sensitivity.

Metabolomics data have been profiled in several epidemiological cohorts [1–6], offering the potential to study the implication of metabolites in human health and disease within and across large-scale datasets. However, individual cohorts are typically profiled using different metabolomics platforms. In order to combine datasets across platforms and cohorts, there is a need to establish the extent of overlap and complementarity across metabolomics platforms.

Several previous studies have explored metabolomics datasets across multiple platforms [7–13]. For example, Suhre et al. [7] used multiple metabolomics platforms in a case-control study of type-2 diabetes (T2D). They profiled 100 individuals using three different metabolomics platforms to assess the potential of using metabolomic data in diabetes research by identifying metabolites that associate with diabetes. The study showed good agreement between known biomarkers of diabetes, including sugar metabolites, that could be replicated by the multiple metabolomic platform approach. Psychogios et al. [8] aimed to characterize the human serum metabolome by combining targeted and non-targeted NMR, GC-MS and LC-MS methods to identify a comprehensive set of metabolites commonly detected and quantified in human serum samples. They reported good agreement between the measured concentrations of NMR and GC-MS. Nicholson et al. [12] and Raffler et al. [13] studied genetic influences on NMR derived urine and plasma metabolites along with MS derived metabolites. However, these studies did not extensively compare the genome-wide findings for metabolite profiles from the same individuals to assess whether associations from datasets across platforms overlap.

In our study, we focus on the comparison of metabolites that are quantified on targeted and non-targeted mass spectrometry platforms and on the comparison of their genetic associations across platforms. Two of the most commonly used high-throughput techniques in large cohort studies apply either a targeted approach using the Biocrates platform or a non-targeted approach using the Metabolon platform. The Biocrates method is a quantitative screen of selected metabolites detected with multiple reaction monitoring, neutral loss and precursor ion scans. Metabolites are then quantified by comparison to structurally similar molecules labelled with stable isotopes added to the samples in defined concentrations as internal standards. In contrast, a non-targeted approach such as Metabolon determines relative concentrations of as many metabolites as possible without using internal standards for absolute quantification. The Biocrates AbsoluteIDQ p150 kits have been applied to quantify a targeted set of 163 metabolites, focusing predominantly on lipids. On the other hand, Metabolon has used ultra high-performance liquid chromatography coupled to tandem mass spectrometry (UHPLC/MS/MS) and gas chromatography coupled to mass spectrometry (GC/MS) for measuring around 500 metabolites from all major pathways including lipids, amino-acids, xenobiotics, and unknown compounds. Although, the methods for quantifying metabolites are distinct, there is an overlap of 43 metabolites that are measured by both platforms. Both platforms focus on different pathways, and combining datasets across platforms can help uncover a wide spectrum of complementary metabolites.

In this study we aimed to compare the Biocrates and Metabolon platforms by integrating human genetic data in a genome-wide association study design. Genome-wide association studies of metabolomic profiles (mGWAS) provide a new approach to evaluate the impact of genetic variation on human metabolism and its indirect link to complex diseases [12,14–20]. A number of studies have reported strong associations between human genetic variants and metabolites from both targeted and non-targeted metabolomics platforms [9,21–25]. The results have identified biologically meaningful associations and in some cases have been used to predict unknown gene function or metabolite identity. We propose to use mGWAS as a method of assessing biologically relevant overlap and complementarity between platforms, as the results could identify metabolites that capture shared biological processes through harmonization of two metabolomics platforms.

We present mGWAS results of metabolites measured across the two platforms in the same set of serum samples from 1,001 individuals. Our aim was to identify metabolites across platforms with consistent genetic associations, which therefore appear stable and robust across multiple platforms. The results can be used to assess how well different metabolomics profiling methods identify identical molecules, to identify metabolites under shared genetic influences, and ultimately to help identify potential metabolites for which data could be combined in future studies. Our approach shows that the different technologies are predominantly complementary in the type and set of metabolites covered.

Materials and Methods

Ethics Statement

The study was approved by St. Thomas’ Hospital Research Ethics Committee, and all twins provided informed written consent.

Study Population and Sample collection

The 1,001 participants in this study were selected from the TwinsUK cohort [26]. Tests and questionnaires applied to the participants have been described elsewhere [3]. The sample consisted of 79 monozygotic (MZ) twin pairs, 215 dizygotic (DZ) twin pairs, and 413 unrelated individuals. TwinsUK blood serum samples for Metabolon and Biocrates platform were obtained after at least 6 hour of fasting and were inverted three times, followed by 40 min resting at 4°C to obtain complete coagulation. The samples were then centrifuged for 10 min at 2,000g. Serum was removed from the centrifuged tubes as the top yellow translucent layer of liquid. Four aliquots of 1.5 ml were placed into skirted micro-centrifuge tubes and then stored in a −45°C freezer until sampling.

Metabolomics Measurements

The same serum samples from 1,001 individuals in this study were profiled on two separate MS platforms, Biocrates and Metabolon. The Biocrates metabolomics data were generated from Helmholtz Center Munich using AbsoluteIDQ™ p150 kits provided by Biocrates Life Sciences AG (Innsbruck, Austria). The Metabolon metabolomic data were generated from Metabolon Inc. (Durham, USA). Biocrates kits were applied to quantify a targeted set of 163 stable metabolites, while Metabolon uses a non-targeted approach for measuring 499 metabolites.

The TwinsUK dataset generated on the targeted Biocrates MS platform has previously been described [22,27,28]. Sample preparation and measurements were performed as illustrated in [25]. Briefly, after centrifugation, 10 μL of serum was pipetted into a 96 well sandwich plate, which contained inserted filters holding stable isotope labeled internal standards. After drying the filters, amino acids were derivatized with 5% phenylisothiocyanate reagent (PITC). From the dried filters, metabolites and internal standards were extracted with 5 mM ammonium acetate in methanol. The solution was centrifuged and diluted with MS running solvent. Liquid handling was performed on a Hamilton Microlab STAR robot (Hamilton Bonaduz AG, Bonaduz, Switzerland). The prepared samples were analyzed by Flow Injection Analysis (FIA) tandem MS with Electrospray Ionization (ESI) on an API 4000 mass spectrometer (AB Sciex Deutschland GmbH, Darmstadt, Germany) using multiple reaction monitoring (MRM). The internal standards served as references for calculating absolute metabolite concentrations in micromolar units (μM). The Biocrates metabolomics dataset contains 163 targeted metabolites: 41 acylcarnitines [Cx:y], hydroxylacylcarnitines [C(OH)x:y] and dicarboxylacylcarnitines [Cx:y-DC]; 14 amino acids; 1 sugar; 15 sphingomyelins [SMx:y] and sphingomyelin-derivatives [SM(OH)x:y]; and 92 glycerophospholipids [PC and lysoPC]. Glycerophospholipids are differentiated with respect to the presence of ester (a) and ether (e) bonds in the glycerol moiety, where two letters (aa = diacyl, ae = acyl-alkyl) denote that two glycerol positions are bound to a fatty acid residue, while a single letter (a = acyl) indicates the presence of a single fatty acid residue. Lipid side chain composition is abbreviated as Cx:y, where x denotes the number of carbons in the side chain and y the number of double bonds. Further descriptions of the 163 Biocrates metabolites have previously been published [27–29].

The TwinsUK dataset generated on the non-targeted MS platform Metabolon has also previously been described [22,23,30] and in this study we report results from a subset of 1,001 individuals from the overall sample. Sample preparation, measurement and metabolite identification have been performed as illustrated in [31,32]. Briefly, metabolites were extracted from 100 μl serum with 400 μl methanol (containing recovery standards) in a 96- deep well plate format. After centrifugation, the supernatant was split into four aliquots per sample: two for two separate ultra-high performance liquid chromatography/MS (UHPLC/MS) injections, one for gas chromatography/MS (GC/MS) injection, and one reserve aliquot. After drying, the aliquots were reconstituted with 0.1% formic acid, for LC/MS positive ion mode, and with 6.5 mM ammonium bicarbonate pH 8.0 for negative ion mode. The GC/MS aliquots were derivatized for 1 h at 60°C with N, O-bistrimethylsilyl-trifluoroacetamide in a solvent mixture of acetonitrile:dichloromethane: cyclohexane (5:4:1), containing 5% triethylamine and retention time markers. Pipetting was performed on a Hamilton MLStar (Hamilton Company, Salt Lake City, UT, USA) robotics system. UHPLC/MS analysis was performed on an LTQ mass spectrometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) equipped with a Waters Acquity UPLC system (Waters Corporation, Milford, MA, USA). Full scan mass spectra (99–1000 m/z) and data dependent MS/MS scans with dynamic exclusion were recorded in turns. GC/MS analysis was done on a Thermo-Finnigan Trace DSQ fast-scanning single-quadrupole mass spectrometer, equipped with a 20 m x 0.18 mm GC column with 0.18 μm film phase consisting of 5% phenyldimethylsilicone. Mass spectra in a scan range from 50–750 m/z were recorded. For metabolite identification, the generated spectral data were compared against an in-house library, which includes retention time (RT), and reference spectra from mass scan and fragmentation of molecules. For every metabolite, the raw area counts were normalized to the median value of the run day to correct for inter-day variation of the measurements. The set of 499 quantified metabolites consists of several classes of named metabolites (amino acids, acylcarnitines, sphingomyelins, glycerophospholipids, carbohydrates, vitamins, lipids, nucleotides, peptides, xenobiotics and steroids) and so-called unknown metabolites of yet unidentified chemical structure (e.g. X-11521).

Genotyping and Imputation

Genotyping of the TwinsUK dataset was performed using a combination of Illumina arrays (HumanHap300, HumanHap610Q, 1M-Duo and 1.2MDuo). We pooled the normalized intensity data and called genotypes on the basis of the Illluminus algorithm. No calls were assigned if the most likely call had a posterior probability less than 0.95. We excluded SNPs with Hardy–Weinberg (P < 1x 10⁻⁷) and and with minor allele frequency < 1%. First, the sparser HumanHap300 dataset was imputed to the HumanHap610Q using phased TwinsUK HumanHap610Q haplotypes as a reference. Next, the combined panel was imputed using reference haplotypes from the HapMap2 project (rel 22, combined CEU+YRI+ASN panels). The genotyping and imputation steps for TwinsUK cohort have been described in detail previously [22,23].

Statistical Analysis

The Biocrates and Metabolon metabolomics datasets in the 1,001 serum samples first underwent several quality control checks. Both dataset were investigated for missingness at the level of each metabolite and individual. Metabolites or individuals with missing values greater than 15% were excluded from further analysis. Outliers at more than 4 standard deviations from the mean of each metabolite were excluded. In total, 11 metabolites were removed from the Metabolon dataset (out of 499 total) and 3 metabolite were removed from Biocrates dataset (out of 163 total)(S1 Table). We next performed Principal Component Analysis (PCA) on the metabolomics profiles in each dataset and compared the first 5 principal components with potential covariates to assess which variables should be included in downstream analyses. Sex, age and BMI were nominally associated with at least 1 principal component and as a result were included as covariates in the downstream analyses.

Altogether, there were 488 (Metabolon) and 160 (Biocrates) metabolites that passed quality control checks, and of these 43 metabolites overlapped, that is, were assigned to be the same molecule by both detection technologies. In the case of lyso-phosphatidylcholines (lysoPCs), the two platforms actually measure not the same but similar molecules: while Metabolon can differentiate between the position of the fatty acid residue on the glycerol backbone (e.g. 1-arachidonoylglycerophosphocholine and 2-arachidonoylglycerophosphocholine), Biocrates measures the sum concentration of both molecules (e.g. lysoPC aa C20:4). Pearson correlation was computed between the metabolite profiles across platforms to assess similarities in metabolite measurements. Several approaches can be used to normalize metabolite data, for example, log transformation [23], inverse normalization [19], and others. Here we used log transformation (base 10) after quantile normalization since test of normality showed that in most cases the normalized concentrations were closer to a normal distribution than the untransformed values. Hierarchical clustering of the metabolites was performed using the complete linkage method that finds similar clusters. All metabolomics quality control analyses were performed using R 3.0.1 (r-project.org).

Initial platform comparison focused on correlation analysis of the 43 metabolites across the two platforms. Follow up platform comparisons included genetic data for biological interpretation of platform overlap. Here, we first calculated twin-based heritability of the metabolite profiles to identify genetically stable and robust profiles across platforms [33]. Second, we used a GWAS approach to identify specific genetic variants that were associated with metabolite levels across platforms.

Heritability was computed for 43 metabolites by comparing metabolite profiles in MZ and DZ twin pairs using the ACE (additive genetic effects (A), common environment (C), and unique environment (E)) model in the OpenMx software [34]. The goal of these analyses was to establish the influence of genetic effects on metabolite profiles, to identify stable genetically determined metabolites, and to relate the results to the mGWAS findings.

To further assess evidence for genetic impacts on metabolites, we performed mGWAS analyses aiming to identify metabolite Quantitative Trait Loci (mQTLs), that is, genetic loci at which genetic variants associated with metabolite levels. We performed mGWAS using GEMMA [35], which implements a genome-wide efficient mixed model association algorithm specifically suitable for the analysis of related individuals, and provides exact P-values from linear mixed models. GEMMA tests for association between each metabolite and each SNP, using one of three commonly used test statistics (the Wald test, the likelihood ratio or score). Here we report all three statistics, but consider the Wald test when setting thresholds. We used Bonferroni correction to account for multiple testing, resulting in genome-wide significance thresholds of P = 3x10^-10 for Biocrates and P = 1x10^-10 for Metabolon. The mGWAS analyses were performed using common SNPs, but both common and rare genetic variants can influence metabolite profiles. The heritability results identify metabolites that are genetically determined, and these effects can be due to either common or rare genetic variants. Therefore some of the heritability effects, especially those underlying rare variants, may not be captured by the mGWAS results.

Results

Platform comparison: correlation and heritability of metabolites profiles

Following quality control assessment, there were 488 (Metabolon) and 160 (Biocrates) metabolites available for analysis in serum samples from 1,001 individuals. Of these, 43 were designated as overlapping molecule by both platforms (S2 Table). Comparisons of the 43 metabolites showed a mean correlation coefficient (r) of 0.44 with a maximum correlation for octanoylcarnitine (r = 0.92), minimum correlation for 1-docosahexaenoylglycerophosphocholine (r = 0), and weak correlations (0<r<0.2) for 7 metabolites (S2 Table), which included lipids and an amino acid. Using hierarchical clustering of the correlation matrix, we observed that the metabolites tend to cluster first within platform, and then within type of the metabolite (Fig 1). One clear exception is hexose (Biocrates), which clusters with glucose in the Metabolon cluster, as expected. A second exception is carnitine C0 (Biocrates), which clusters near proline, valine, tyrosine, and propionylcarnitine in the Metabolon cluster. Additionally, we calculated the correlation between the 43 metabolites and all remaining metabolites on both alternative platforms. We observed that the resulting correlations were overall much lower, and only two pairs of metabolites across platforms had a mean correlation of 0.44 or greater; these included octanoylcarnitine C8 (Biocrates) and the unknown metabolite X-11521 (Metabolon), and laurylcarnitine C12 (Biocrates) and the unknown metabolite X-18739 (Metabolon).

Download:

Fig 1. Hierarchical cluster of the correlation across 43 overlapping metabolites from both platforms.

Upper colour bars represent metabolites with mGWAS results, metabolite type, and metabolite platform. The left colour bar represents the heritability of the metabolite from red (high) to white (low).

https://doi.org/10.1371/journal.pone.0153672.g001

Because the 1,001 individuals included twins, we were able to calculate twin-based heritability estimates of the metabolite profiles, focusing on the 43 overlapping metabolites (S2 Table). Of the 43 metabolites, 37 (Biocrates) and 34 (Metabolon) were at least moderately heritable in twins (h²>0.2). There were 29 metabolites with evidence for heritability on both platforms (h² ranging from 0.29 to 0.72, S2 Table). Of these, the 9 most heritable profiles were observed for 6 lipids (h²: 0.4 to 0.72) and 3 amino acids (h²: 0.42 to 0.7), indicating that these are stable profiles and highly likely to be under genetic influence.

mGWAS results: overlapping and complementary mQTLs

In total, 488 and 160 metabolites were tested separately on the Metabolon and Biocrates platforms in two mGWAS analyses. All genome-wide significant association results are reported at a stringent Bonferroni cut-off: P = 1×10⁻¹⁰ (5×10⁻⁸/488) for Metabolon and P = 3×10⁻¹⁰ (5×10⁻⁸/160) for Biocrates. Additionally a relaxed threshold for genome-wide association (5×10⁻⁸) was used to evaluate whether a mGWAS finding on one platform was replicated on the other platform (S3 Table). Additionally, we provide all result pairs where metabolites on both platforms surpassed evidence for genetic association at P = 5x10e-8 (S4 Table)).

In total, 61 genome-wide significant metabolite associations were identified at 26 independent loci: 42 metabolites were associated with 25 loci on the Metabolon platform, and 19 metabolites were associated with 8 loci on the Biocrates platform (Table 1). Of the 26 independent loci, genome-wide significant metabolite associations at 7 loci were identified on both platforms. There were 19 loci that had associations only with metabolites from one platform (18 loci in Metabolon and 1 locus in Biocrates).

Download:

Table 1. Genome-wide significant mGWAS results.

https://doi.org/10.1371/journal.pone.0153672.t001

Overlapping mQTLs: genetic associations identified on both platforms

Associations at 7 independent loci were identified in both platforms, namely with SNPs in the regions of the ACADM, ACADL, CPS1, SLC16A9, FADS1, ACADS and SGPP1 genes (Table 2). The 7 loci associate with 22 metabolites in total: 9 metabolites from Biocrates and 13 metabolites from Metabolon.

Download:

Table 2. mGWAS results at 7 loci associated with metabolites in both platforms.

https://doi.org/10.1371/journal.pone.0153672.t002

Of the 22 associated metabolites, 6 metabolites associated with 5 loci were named for the overlapping compound on both platforms. These included C6 (Biocrates, P = 4.1×10⁻¹¹) = hexanoylcarnitine (Metabolon, P = 1.6×10⁻¹³), C8 (Biocrates, P = 2.4×10⁻⁸) = octanoylcarnitine (Metabolon, P = 4.8×10⁻¹¹), glycine (Biocrates, P = 5.3×10⁻¹⁷) = glycine (Metabolon, P = 7.1×10⁻²⁷), C0 (Biocrates, P = 4.6×10⁻¹²) = carnitine (Metabolon, P = 2.3×10⁻¹³), C4 (Biocrates, P = 2.9×10⁻⁴⁴) = butyrylcarnitine (Metabolon, P = 1.75×10⁻¹¹⁴), and lysoPC a C20:4 (Biocrates, P = 2×10⁻¹⁴) = 1-arachidonoylglycerophosphocholine (Metabolon, P = 2.9×10⁻¹⁰), as designated by Biocrates and Metabolon, respectively. For three of the 5 loci with smatching named metabolites, there were also associations with other metabolites, which do not necessarily match across platforms (Table 2).

In one case genetic variants in locus ACADL were associated with both a Biocrates metabolite C9 (P = 9.7×10⁻³⁸) and an unknown Metabolon metabolite (X-13431 (P = 3.5×10⁻²⁵)), which were recently shown to be identical molecules [36]. The mean correlation coefficient between these metabolites across platforms was moderate (r = 0.54, Fig 1).

In one case, metabolite associations with genetic variants at the SGPP1 locus did not match in name for PC aa C28:1 (Biocrates) and 1-stearoylglycerol (Metabolon) (Table 2). The mean correlation coefficient between these metabolites across platforms is moderate (r = 0.42, Fig 1). Both of these are lipid metabolites, and could share the C18:0 fatty acid chain.

Complementary mQTLs: genetic associations identified in only one platform

There were 19 loci that had associations only with metabolites from one platform (18 loci in Metabolon and 1 locus in Biocrates) and these all were associated with metabolites that were not measured in the other platform (S3 Table).

The 18 Metabolon-specific mGWAS results included associations with 29 metabolites. Of these 29 metabolites, 17 were unknowns, 4 were lipids and 3 were amino acids and these were not included in Biocrates, considering that Biocrates consists mostly of lipids and amino acids. The 5 remaining metabolites were 2 drugs, a carbohydrate, a nucleotide, and a peptide.

There was only 1 locus (DYNC1H1) where genetic variants showed genome-wide significant mGWAS results on the Biocrates platform only with 4 metabolites, and in all 4 cases these were with lipids that Metabolon did not measure.

Discussion

Our study is a bi-platform metabolite comparison using mGWAS with the objective of identifying metabolites measured on more than one platform where signals overlap and may be combined in future studies, for example for replication analysis. The key results identified 7 loci showing robust genetic associations with metabolites on both platforms. These results were also predominantly consistent with recent reported mGWAS [22,23,30,37], some of which are based on results from extended cohorts that include the samples used in the current analysis. Thus, for 6 of the 7 loci (ACADM, ACADL, CPS1, SLC16A9, FADS1, ACADS), previous mGWAS reported associations with the same Metabolon metabolite either as a single metabolite or as part of a metabolite ratio [23,30]. In contrast, SGPP1 harboured an mQTL with the Metabolon metabolite ratio (X-08402/cholesterol), and the single metabolites X-08402 and X-10510 in Shin et al. [28], while here we report associations with 1-stearoylglycerol and X-10510.

Of the metabolites associated with the 7 loci, 5 metabolites (Biocrates C8, C6, C0, C4, and glycine) had at least moderate heritability (h²>0.26) and correlation (>0.38) on both platforms, confirming that these profiles are stable and reproducible across platforms. Interestingly 1 matching metabolite, lysoPC a C20:4 [Biocrates] / 1—arachidonoylglycerophosphocholine [Metabolon], showed low heritability in one platform (0.09 in Metabolon and 0.59 in Biocrates platform) and showed relatively low correlation (r = 0.29) across platforms, but was still identified to associate with the same locus from both platforms at genome-wide significance. This observation may be due to the difference in the measured compounds between the two platforms: while Metabolon specifically quantifies the lysoPC with the 20:4 fatty acid chain at sn1 position of the glycerol backbone (lysoPC(20:4/0:0), Biocrates does not distinguish between the lysoPCs with fatty acid chains at sn1 and sn2 positions and only quantifies the sum concentration of the two forms (lysoPC(20:4/0:0 and lysoPC(0:0/20:4). Moreover, the quality of measurement differs for various lipids between the targeted Biocrates and the non-targeted Metabolon platform, which might also cause lower correlation between the corresponding matching metabolites. Notably, despite those differences inherent in the platforms both profiles give a robust signal of genetic association for FADS1.

Further comparison of the GWAS results across platforms shows that genetic variants at 5 of the 7 loci (ACADM, CPS1, SLC16A9, FADS1, ACADS) were associated with metabolites that were named for the overlapping compound. However, genetic variants at the ACADL and SGPP1 loci only associate with non-overlapping metabolites or unknown metabolites from the Metabolon platform. In these cases, our results can be used to inform the function of unknown metabolites or identify metabolites that belong to the same or related biological pathways. For example, variants in the ACADL locus associated with the C9 Biocrates metabolite and also with the unknown X-13431 Metabolon metabolite, which were recently reported to be the same molecule [34]. When we explored the results for similar association patterns, we observed that Metabolon metabolites X-10510 and 1-stearoylglycerol shared mQTL findings within the same locus (SGGP1) as the Biocrates metabolite PC aa C28:1. These results suggest a link between the molecules, where the more specific Metabolon lipid chain length can hint that the PC aa C28:1 association is possibly driven by the involvement of a 18:0 lipid chain. Alternatively, the SGGP1 genetic variant (rs7157785) has also been associated with sphingomyelin 14:0 in a separate study [35]. Our platform does not include this metabolite, but X-10510 may be also related to this sphingolipid pathway. This assumption is further supported by high partial correlation between X-10510 and other Metabolon sphingolipid molecules and genetic associations to a second sphingolipid related gene in Shin et al. [30].

We next explored the 43 overlapping metabolites on both platforms for consistencies and potential inconsistencies across platform signals beyond their association results. As expected, the mean correlation between the 43 matching metabolites (r = 0.44) is higher than the mean correlation with all metabolites between the two platforms (r = 0.17). Exceptions include correlations of Biocrates metabolites with Metabolon metabolites of yet unknown chemical identity. In these cases, the high correlation could indicate matching metabolites or biochemically related metabolites and might thus again assist in the identification of unknown metabolites.

Four lyso-phosphatidylcholine metabolites (lysoPC a C16:0, lysoPC a C18:0, lysoPC a C18:1, lysoPC a C18:2) from the Biocrates platforms had overlapping metabolites on the Metabolon platform, but neither contained matching mQTLs nor showed high heritability or correlation. We conclude that in this instance the two platforms are likely measuring distinct signals that cannot be combined or this may be due to a relatively lower quality of measurement for these lipids on the Metabolon platform.

We applied a combination of correlation, heritability and genotypic analyses to bring together the comparison of data from different metabolomics platforms. Our approach identified genetic associations at 7 loci with pairs of metabolites measured on the two platforms that were named for the same compound, were highly correlated and heritable, therefore suggesting that in these cases Biocrates and Metabolon signals overlap. In contrast, 9 pairs of known metabolites that are not named for the same compound across platforms, but exhibit similar levels of correlation and heritability, showed no overlapping genetic associations. The two platforms are designed to focus on different metabolites, and these findings can inform on platform-specific metabolites. Ultimately, combining metabolomics profiles across platforms is more informative than single-platform analysis because platforms are complementary. It is not possible to assay the entire metabolome with one platform due to large differences in the physiochemical properties of the different metabolites (e.g. lipophilic and hydrophilic metabolites).

In summary, we identified genetic associations at 7 loci with metabolite profiles from both the Biocrates and Metabolon platforms. Our results provide new information about potential shared pathways, as well as distinct metabolite profiles, and their genetic determinants, clarifying unknown metabolites. Our study demonstrates the complementary nature of both targeted and non-targeted MS platforms and can help future studies to explore combining datasets across platforms, especially for replication of metabolite hits when datasets are profiled on different platforms. The findings can help guide further research into the sources of inconsistency and variation in the comparison of metabolite results profiled from differing platforms.

Supporting Information

S1 Table. Metabolites removed from analysis at Quality Control.

https://doi.org/10.1371/journal.pone.0153672.s001

(XLSX)

S2 Table. Metabolon and Biocrates platform comparison at 43 overlapping metabolites including correlation, heritability and peak mGWAS results.

https://doi.org/10.1371/journal.pone.0153672.s002

(XLS)

S3 Table. mGWAS results for Biocrates and Metabolon platforms.

https://doi.org/10.1371/journal.pone.0153672.s003

(XLSX)

S4 Table. mGWAS Results for Biocrates and Metabolon with relaxed threshold (5e-8).

https://doi.org/10.1371/journal.pone.0153672.s004

(XLSX)

Acknowledgments

The study was funded in part by the Wellcome Trust and the European Community’s Seventh Framework Programme (FP7/2007-2013). TwinsUK also received support from the Wellcome Trust, and the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. SNP Genotyping was performed by The Wellcome Trust Sanger Institute and National Eye Institute via NIH/CIDR.

Author Contributions

Conceived and designed the experiments: GK TDS JTB. Analyzed the data: IY. Contributed reagents/materials/analysis tools: CM SYS MM NS JA KS TDS. Wrote the paper: IY GK JTB.

References

1. Holle R, Happich M, Lowel H, Wichmann HE, Group MKS (2005) KORA—a research platform for population based health research. Gesundheitswesen 67 Suppl 1: S19–25. pmid:16032513
- View Article
- PubMed/NCBI
- Google Scholar
2. Wichmann HE, Gieger C, Illig T, Group MKS (2005) KORA-gen—resource for population genetics, controls and a broad spectrum of disease phenotypes. Gesundheitswesen 67 Suppl 1: S26–30. pmid:16032514
- View Article
- PubMed/NCBI
- Google Scholar
3. Moayyeri A, Hammond CJ, Hart DJ, Spector TD (2013) The UK Adult Twin Registry (TwinsUK Resource). Twin Res Hum Genet 16: 144–149. pmid:23088889
- View Article
- PubMed/NCBI
- Google Scholar
4. Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, Auro K, et al. (2008) Gender differences in genetic risk profiles for cardiovascular disease. PLoS One 3: e3615. pmid:18974842
- View Article
- PubMed/NCBI
- Google Scholar
5. Boomsma DI, de Geus EJ, Vink JM, Stubbe JH, Distel MA, Hottenga JJ, et al. (2006) Netherlands Twin Register: from twins to twin families. Twin Res Hum Genet 9: 849–857. pmid:17254420
- View Article
- PubMed/NCBI
- Google Scholar
6. Riboli E, Kaaks R (1997) The EPIC Project: rationale and study design. European Prospective Investigation into Cancer and Nutrition. Int J Epidemiol 26 Suppl 1: S6–14. pmid:9126529
- View Article
- PubMed/NCBI
- Google Scholar
7. Suhre K, Meisinger C, Doring A, Altmaier E, Belcredi P, Gieger C, et al. (2010) Metabolic footprint of diabetes: a multiplatform metabolomics study in an epidemiological setting. PLoS One 5: e13953. pmid:21085649
- View Article
- PubMed/NCBI
- Google Scholar
8. Psychogios N, Hau DD, Peng J, Guo AC, Mandal R, Bouatra S, et al. (2011) The human serum metabolome. PLoS One 6: e16957. pmid:21359215
- View Article
- PubMed/NCBI
- Google Scholar
9. Adamski J, Suhre K (2013) Metabolomics platforms for genome wide association studies—linking the genome to the metabolome. Curr Opin Biotechnol 24: 39–47. pmid:23102864
- View Article
- PubMed/NCBI
- Google Scholar
10. Mandal R, Guo AC, Chaudhary KK, Liu P, Yallou FS, Dong E, et al. (2012) Multi-platform characterization of the human cerebrospinal fluid metabolome: a comprehensive and quantitative update. Genome Med 4: 38. pmid:22546835
- View Article
- PubMed/NCBI
- Google Scholar
11. Büscher JM, Czernik D, Ewald JC, Sauer U, Zamboni N (2009) Cross-platform comparison of methods for quantitative metabolomics of primary metabolism. Anal Chem 81: 2135–2143. pmid:19236023
- View Article
- PubMed/NCBI
- Google Scholar
12. Nicholson G, Rantalainen M, Li JV, Maher AD, Malmodin D, Ahmadi KR, et al. (2011) A genome-wide metabolic QTL analysis in Europeans implicates two loci shaped by recent positive selection. PLoS Genet 7: e1002270. pmid:21931564
- View Article
- PubMed/NCBI
- Google Scholar
13. Raffler J, Römisch-Margl W, Petersen AK, Pagel P, Blöchl F, Hengstenberg C, et al. (2013) Identification and MS-assisted interpretation of genetically influenced NMR signals in human plasma. Genome Med 5: 13. pmid:23414815
- View Article
- PubMed/NCBI
- Google Scholar
14. Assfalg M, Bertini I, Colangiuli D, Luchinat C, Schäfer H, Schütz B, et al. (2008) Evidence of different metabolic phenotypes in humans. Proc Natl Acad Sci U S A 105: 1420–1424. pmid:18230739
- View Article
- PubMed/NCBI
- Google Scholar
15. Dumas ME, Wilder SP, Bihoreau MT, Barton RH, Fearnside JF, Argoud K, et al. (2007) Direct quantitative trait locus mapping of mammalian metabolic phenotypes in diabetic and normoglycemic rat models. Nat Genet 39: 666–672. pmid:17435758
- View Article
- PubMed/NCBI
- Google Scholar
16. Vinayavekhin N, Homan EA, Saghatelian A (2010) Exploring disease through metabolomics. ACS Chem Biol 5: 91–103. pmid:20020774
- View Article
- PubMed/NCBI
- Google Scholar
17. Watson AD (2006) Thematic review series: systems biology approaches to metabolic and cardiovascular disorders. Lipidomics: a global approach to lipid analysis in biological systems. J Lipid Res 47: 2101–2111. pmid:16902246
- View Article
- PubMed/NCBI
- Google Scholar
18. Homuth G, Teumer A, Völker U, Nauck M (2012) A description of large-scale metabolomics studies: increasing value by combining metabolomics with genome-wide SNP genotyping and transcriptional profiling. J Endocrinol 215: 17–28. pmid:22782382
- View Article
- PubMed/NCBI
- Google Scholar
19. Kettunen J, Tukiainen T, Sarin AP, Ortega-Alonso A, Tikkanen E, Lyytikainen LP, et al. (2012) Genome-wide association study identifies multiple loci influencing human serum metabolite levels. Nat Genet 44: 269–276. pmid:22286219
- View Article
- PubMed/NCBI
- Google Scholar
20. Rueedi R, Ledda M, Nicholls AW, Salek RM, Marques-Vidal P, Morya E, et al. (2014) Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links. PLoS Genet 10: e1004132. pmid:24586186
- View Article
- PubMed/NCBI
- Google Scholar
21. Gieger C, Geistlinger L, Altmaier E, Hrabe de Angelis M, Kronenberg F, Meitinger T, et al. (2008) Genetics meets metabolomics: a genome-wide association study of metabolite profiles in human serum. PLoS Genet 4: e1000282. pmid:19043545
- View Article
- PubMed/NCBI
- Google Scholar
22. Illig T, Gieger C, Zhai G, Romisch-Margl W, Wang-Sattler R, Prehn C, et al. (2010) A genome-wide perspective of genetic variation in human metabolism. Nat Genet 42: 137–141. pmid:20037589
- View Article
- PubMed/NCBI
- Google Scholar
23. Suhre K, Shin SY, Petersen AK, Mohney RP, Meredith D, Wägele B, et al. (2011) Human metabolic individuality in biomedical and pharmaceutical research. Nature 477: 54–60. pmid:21886157
- View Article
- PubMed/NCBI
- Google Scholar
24. Suhre K, Wallaschofski H, Raffler J, Friedrich N, Haring R, Michael K, et al. (2011) A genome-wide association study of metabolic traits in human urine. Nat Genet 43: 565–569. pmid:21572414
- View Article
- PubMed/NCBI
- Google Scholar
25. Kastenmuller G, Raffler J, Gieger C, Suhre K (2015) Genetics of human metabolism: an update. Hum Mol Genet.
- View Article
- Google Scholar
26. Moayyeri A, Hammond CJ, Valdes AM, Spector TD (2013) Cohort Profile: TwinsUK and healthy ageing twin study. Int J Epidemiol 42: 76–85. pmid:22253318
- View Article
- PubMed/NCBI
- Google Scholar
27. Römisch-Margl W, Prehn C, Bogumil R, Röhring C, Suhre K, Adamski J (2012) Procedure for tissue sample preparation and metabolite extraction for high-throughput targeted metabolomics. Metabolomics 8: 133–142.
- View Article
- Google Scholar
28. Mittelstrass K, Ried JS, Yu Z, Krumsiek J, Gieger C, Prehn C, et al. (2011) Discovery of sexual dimorphisms in metabolic and genetic biomarkers. PLoS Genet 7: e1002215. pmid:21852955
- View Article
- PubMed/NCBI
- Google Scholar
29. Menni C, Zhai G, Macgregor A, Prehn C, Römisch-Margl W, Suhre K, et al. (2013) Targeted metabolomics profiles are strongly correlated with nutritional patterns in women. Metabolomics 9: 506–514. pmid:23543136
- View Article
- PubMed/NCBI
- Google Scholar
30. Shin SY, Fauman EB, Petersen AK, Krumsiek J, Santos R, Huang J, et al. (2014) An atlas of genetic influences on human blood metabolites. Nat Genet 46: 543–550. pmid:24816252
- View Article
- PubMed/NCBI
- Google Scholar
31. Dehaven CD, Evans AM, Dai H, Lawton KA (2010) Organization of GC/MS and LC/MS metabolomics data into chemical libraries. J Cheminform 2: 9. pmid:20955607
- View Article
- PubMed/NCBI
- Google Scholar
32. Evans AM, DeHaven CD, Barrett T, Mitchell M, Milgram E (2009) Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal Chem 81: 6656–6667. pmid:19624122
- View Article
- PubMed/NCBI
- Google Scholar
33. Yousri NA, Kastenmuller G, Gieger C, Shin SY, Erte I, Menni C, et al. (2014) Long term conservation of human metabolic phenotypes and link to heritability. Metabolomics 10: 1005–1017. pmid:25177233
- View Article
- PubMed/NCBI
- Google Scholar
34. Boker S, Neale M, Maes H, Wilde M, Spiegel M, Brick T, et al. (2011) OpenMx: An Open Source Extended Structural Equation Modeling Framework. Psychometrika 76: 306–317. pmid:23258944
- View Article
- PubMed/NCBI
- Google Scholar
35. Zhou X, Stephens M (2012) Genome-wide efficient mixed-model analysis for association studies. Nat Genet 44: 821–824. pmid:22706312
- View Article
- PubMed/NCBI
- Google Scholar
36. Krumsiek J, Suhre K, Evans AM, Mitchell MW, Mohney RP, Milburn MV, et al. (2012) Mining the unknown: a systems approach to metabolite identification combining genetic and metabolic information. PLoS Genet 8: e1003005. pmid:23093944
- View Article
- PubMed/NCBI
- Google Scholar
37. Draisma HH, Pool R, Kobl M, Jansen R, Petersen AK, Vaarhorst AA, et al. (2015) Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels. Nat Commun 6: 7208. pmid:26068415
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Holle R, Happich M, Lowel H, Wichmann HE, Group MKS (2005) KORA—a research platform for population based health research. Gesundheitswesen 67 Suppl 1: S19–25. pmid:16032513
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Wichmann HE, Gieger C, Illig T, Group MKS (2005) KORA-gen—resource for population genetics, controls and a broad spectrum of disease phenotypes. Gesundheitswesen 67 Suppl 1: S26–30. pmid:16032514
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Moayyeri A, Hammond CJ, Hart DJ, Spector TD (2013) The UK Adult Twin Registry (TwinsUK Resource). Twin Res Hum Genet 16: 144–149. pmid:23088889
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, Auro K, et al. (2008) Gender differences in genetic risk profiles for cardiovascular disease. PLoS One 3: e3615. pmid:18974842
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Boomsma DI, de Geus EJ, Vink JM, Stubbe JH, Distel MA, Hottenga JJ, et al. (2006) Netherlands Twin Register: from twins to twin families. Twin Res Hum Genet 9: 849–857. pmid:17254420
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Riboli E, Kaaks R (1997) The EPIC Project: rationale and study design. European Prospective Investigation into Cancer and Nutrition. Int J Epidemiol 26 Suppl 1: S6–14. pmid:9126529
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Suhre K, Meisinger C, Doring A, Altmaier E, Belcredi P, Gieger C, et al. (2010) Metabolic footprint of diabetes: a multiplatform metabolomics study in an epidemiological setting. PLoS One 5: e13953. pmid:21085649
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Psychogios N, Hau DD, Peng J, Guo AC, Mandal R, Bouatra S, et al. (2011) The human serum metabolome. PLoS One 6: e16957. pmid:21359215
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Adamski J, Suhre K (2013) Metabolomics platforms for genome wide association studies—linking the genome to the metabolome. Curr Opin Biotechnol 24: 39–47. pmid:23102864
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Mandal R, Guo AC, Chaudhary KK, Liu P, Yallou FS, Dong E, et al. (2012) Multi-platform characterization of the human cerebrospinal fluid metabolome: a comprehensive and quantitative update. Genome Med 4: 38. pmid:22546835
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Büscher JM, Czernik D, Ewald JC, Sauer U, Zamboni N (2009) Cross-platform comparison of methods for quantitative metabolomics of primary metabolism. Anal Chem 81: 2135–2143. pmid:19236023
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Nicholson G, Rantalainen M, Li JV, Maher AD, Malmodin D, Ahmadi KR, et al. (2011) A genome-wide metabolic QTL analysis in Europeans implicates two loci shaped by recent positive selection. PLoS Genet 7: e1002270. pmid:21931564
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Raffler J, Römisch-Margl W, Petersen AK, Pagel P, Blöchl F, Hengstenberg C, et al. (2013) Identification and MS-assisted interpretation of genetically influenced NMR signals in human plasma. Genome Med 5: 13. pmid:23414815
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Assfalg M, Bertini I, Colangiuli D, Luchinat C, Schäfer H, Schütz B, et al. (2008) Evidence of different metabolic phenotypes in humans. Proc Natl Acad Sci U S A 105: 1420–1424. pmid:18230739
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Dumas ME, Wilder SP, Bihoreau MT, Barton RH, Fearnside JF, Argoud K, et al. (2007) Direct quantitative trait locus mapping of mammalian metabolic phenotypes in diabetic and normoglycemic rat models. Nat Genet 39: 666–672. pmid:17435758
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Vinayavekhin N, Homan EA, Saghatelian A (2010) Exploring disease through metabolomics. ACS Chem Biol 5: 91–103. pmid:20020774
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Watson AD (2006) Thematic review series: systems biology approaches to metabolic and cardiovascular disorders. Lipidomics: a global approach to lipid analysis in biological systems. J Lipid Res 47: 2101–2111. pmid:16902246
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref18] 18. Homuth G, Teumer A, Völker U, Nauck M (2012) A description of large-scale metabolomics studies: increasing value by combining metabolomics with genome-wide SNP genotyping and transcriptional profiling. J Endocrinol 215: 17–28. pmid:22782382
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref19] 19. Kettunen J, Tukiainen T, Sarin AP, Ortega-Alonso A, Tikkanen E, Lyytikainen LP, et al. (2012) Genome-wide association study identifies multiple loci influencing human serum metabolite levels. Nat Genet 44: 269–276. pmid:22286219
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref20] 20. Rueedi R, Ledda M, Nicholls AW, Salek RM, Marques-Vidal P, Morya E, et al. (2014) Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links. PLoS Genet 10: e1004132. pmid:24586186
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref21] 21. Gieger C, Geistlinger L, Altmaier E, Hrabe de Angelis M, Kronenberg F, Meitinger T, et al. (2008) Genetics meets metabolomics: a genome-wide association study of metabolite profiles in human serum. PLoS Genet 4: e1000282. pmid:19043545
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref22] 22. Illig T, Gieger C, Zhai G, Romisch-Margl W, Wang-Sattler R, Prehn C, et al. (2010) A genome-wide perspective of genetic variation in human metabolism. Nat Genet 42: 137–141. pmid:20037589
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref23] 23. Suhre K, Shin SY, Petersen AK, Mohney RP, Meredith D, Wägele B, et al. (2011) Human metabolic individuality in biomedical and pharmaceutical research. Nature 477: 54–60. pmid:21886157
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref24] 24. Suhre K, Wallaschofski H, Raffler J, Friedrich N, Haring R, Michael K, et al. (2011) A genome-wide association study of metabolic traits in human urine. Nat Genet 43: 565–569. pmid:21572414
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref25] 25. Kastenmuller G, Raffler J, Gieger C, Suhre K (2015) Genetics of human metabolism: an update. Hum Mol Genet.
View Article
Google Scholar

[98] View Article

[99] Google Scholar

[ref26] 26. Moayyeri A, Hammond CJ, Valdes AM, Spector TD (2013) Cohort Profile: TwinsUK and healthy ageing twin study. Int J Epidemiol 42: 76–85. pmid:22253318
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref27] 27. Römisch-Margl W, Prehn C, Bogumil R, Röhring C, Suhre K, Adamski J (2012) Procedure for tissue sample preparation and metabolite extraction for high-throughput targeted metabolomics. Metabolomics 8: 133–142.
View Article
Google Scholar

[105] View Article

[106] Google Scholar

[ref28] 28. Mittelstrass K, Ried JS, Yu Z, Krumsiek J, Gieger C, Prehn C, et al. (2011) Discovery of sexual dimorphisms in metabolic and genetic biomarkers. PLoS Genet 7: e1002215. pmid:21852955
View Article
PubMed/NCBI
Google Scholar

[108] View Article

[109] PubMed/NCBI

[110] Google Scholar

[ref29] 29. Menni C, Zhai G, Macgregor A, Prehn C, Römisch-Margl W, Suhre K, et al. (2013) Targeted metabolomics profiles are strongly correlated with nutritional patterns in women. Metabolomics 9: 506–514. pmid:23543136
View Article
PubMed/NCBI
Google Scholar

[112] View Article

[113] PubMed/NCBI

[114] Google Scholar

[ref30] 30. Shin SY, Fauman EB, Petersen AK, Krumsiek J, Santos R, Huang J, et al. (2014) An atlas of genetic influences on human blood metabolites. Nat Genet 46: 543–550. pmid:24816252
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref31] 31. Dehaven CD, Evans AM, Dai H, Lawton KA (2010) Organization of GC/MS and LC/MS metabolomics data into chemical libraries. J Cheminform 2: 9. pmid:20955607
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref32] 32. Evans AM, DeHaven CD, Barrett T, Mitchell M, Milgram E (2009) Integrated, nontargeted ultrahigh performance liquid chromatography/electrospray ionization tandem mass spectrometry platform for the identification and relative quantification of the small-molecule complement of biological systems. Anal Chem 81: 6656–6667. pmid:19624122
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref33] 33. Yousri NA, Kastenmuller G, Gieger C, Shin SY, Erte I, Menni C, et al. (2014) Long term conservation of human metabolic phenotypes and link to heritability. Metabolomics 10: 1005–1017. pmid:25177233
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref34] 34. Boker S, Neale M, Maes H, Wilde M, Spiegel M, Brick T, et al. (2011) OpenMx: An Open Source Extended Structural Equation Modeling Framework. Psychometrika 76: 306–317. pmid:23258944
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref35] 35. Zhou X, Stephens M (2012) Genome-wide efficient mixed-model analysis for association studies. Nat Genet 44: 821–824. pmid:22706312
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref36] 36. Krumsiek J, Suhre K, Evans AM, Mitchell MW, Mohney RP, Milburn MV, et al. (2012) Mining the unknown: a systems approach to metabolite identification combining genetic and metabolic information. PLoS Genet 8: e1003005. pmid:23093944
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref37] 37. Draisma HH, Pool R, Kobl M, Jansen R, Petersen AK, Vaarhorst AA, et al. (2015) Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels. Nat Commun 6: 7208. pmid:26068415
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Ethics Statement

Study Population and Sample collection

Metabolomics Measurements

Genotyping and Imputation

Statistical Analysis

Results

Platform comparison: correlation and heritability of metabolites profiles

mGWAS results: overlapping and complementary mQTLs

Overlapping mQTLs: genetic associations identified on both platforms

Complementary mQTLs: genetic associations identified in only one platform

Discussion

Supporting Information

S1 Table. Metabolites removed from analysis at Quality Control.

S2 Table. Metabolon and Biocrates platform comparison at 43 overlapping metabolites including correlation, heritability and peak mGWAS results.

S3 Table. mGWAS results for Biocrates and Metabolon platforms.

S4 Table. mGWAS Results for Biocrates and Metabolon with relaxed threshold (5e-8).

Acknowledgments

Author Contributions

References