Human metabolic individuality in biomedical and pharmaceutical research.

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10-60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn's disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

Epigenetics meets metabolomics: an epigenome-wide association study with blood serum metabolic traits.

Previously, we reported strong influences of genetic variants on metabolic phenotypes, some of them with clinical relevance. Here, we hypothesize that DNA methylation may have an important and potentially independent effect on human metabolism. To test this hypothesis, we conducted what is to the best of our knowledge the first epigenome-wide association study (EWAS) between DNA methylation and metabolic traits (metabotypes) in human blood. We assess 649 blood metabolic traits from 1814 participants of the Kooperative Gesundheitsforschung in der Region Augsburg (KORA) population study for association with methylation of 457 004 CpG sites, determined on the Infinium HumanMethylation450 BeadChip platform. Using the EWAS approach, we identified two types of methylome-metabotype associations. One type is driven by an underlying genetic effect; the other type is independent of genetic variation and potentially driven by common environmental and life-style-dependent factors. We report eight CpG loci at genome-wide significance that have a genetic variant as confounder (P = 3.9 × 10(-20) to 2.0 × 10(-108), r(2) = 0.036 to 0.221). Seven loci display CpG site-specific associations to metabotypes, but do not exhibit any underlying genetic signals (P = 9.2 × 10(-14) to 2.7 × 10(-27), r(2) = 0.008 to 0.107). We further identify several groups of CpG loci that associate with a same metabotype, such as 4-vinylphenol sulfate and 4-androsten-3-beta,17-beta-diol disulfate. In these cases, the association between CpG-methylation and metabotype is likely the result of a common external environmental factor, including smoking. Our study shows that analysis of EWAS with large numbers of metabolic traits in large population cohorts are, in principle, feasible. Taken together, our data suggest that DNA methylation plays an important role in regulating human metabolism.

A genome-wide association meta-analysis of circulating sex hormone-binding globulin reveals multiple Loci implicated in sex steroid hormone regulation.

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8 × 10(-106)), PRMT6 (rs17496332, 1p13.3, p = 1.4 × 10(-11)), GCKR (rs780093, 2p23.3, p = 2.2 × 10(-16)), ZBTB10 (rs440837, 8q21.13, p = 3.4 × 10(-09)), JMJD1C (rs7910927, 10q21.3, p = 6.1 × 10(-35)), SLCO1B1 (rs4149056, 12p12.1, p = 1.9 × 10(-08)), NR2F2 (rs8023580, 15q26.2, p = 8.3 × 10(-12)), ZNF652 (rs2411984, 17q21.32, p = 3.5 × 10(-14)), TDGF3 (rs1573036, Xq22.3, p = 4.1 × 10(-14)), LHCGR (rs10454142, 2p16.3, p = 1.3 × 10(-07)), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7 × 10(-08)), and UGT2B15 (rs293428, 4q13.2, p = 5.5 × 10(-06)). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5 × 10(-08), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ~15.6% and ~8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

Stratifying type 2 diabetes cases by BMI identifies genetic risk variants in LAMA1 and enrichment for risk variants in lean compared to obese cases.

Common diseases such as type 2 diabetes are phenotypically heterogeneous. Obesity is a major risk factor for type 2 diabetes, but patients vary appreciably in body mass index. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m²) compared to obese cases (BMI≥30 Kg/m²). We performed two case-control genome-wide studies using two accepted cut-offs for defining individuals as overweight or obese. We used 2,112 lean type 2 diabetes cases (BMI<25 kg/m²) or 4,123 obese cases (BMI≥30 kg/m²), and 54,412 un-stratified controls. Replication was performed in 2,881 lean cases or 8,702 obese cases, and 18,957 un-stratified controls. To assess the effects of known signals, we tested the individual and combined effects of SNPs representing 36 type 2 diabetes loci. After combining data from discovery and replication datasets, we identified two signals not previously reported in Europeans. A variant (rs8090011) in the LAMA1 gene was associated with type 2 diabetes in lean cases (P = 8.4×10⁻9, OR = 1.13 [95% CI 1.09-1.18]), and this association was stronger than that in obese cases (P = 0.04, OR = 1.03 [95% CI 1.00-1.06]). A variant in HMG20A--previously identified in South Asians but not Europeans--was associated with type 2 diabetes in obese cases (P = 1.3×10⁻8, OR = 1.11 [95% CI 1.07-1.15]), although this association was not significantly stronger than that in lean cases (P = 0.02, OR = 1.09 [95% CI 1.02-1.17]). For 36 known type 2 diabetes loci, 29 had a larger odds ratio in the lean compared to obese (binomial P = 0.0002). In the lean analysis, we observed a weighted per-risk allele OR = 1.13 [95% CI 1.10-1.17], P = 3.2×10⁻¹4. This was larger than the same model fitted in the obese analysis where the OR = 1.06 [95% CI 1.05-1.08], P = 2.2×10⁻¹6. This study provides evidence that stratification of type 2 diabetes cases by BMI may help identify additional risk variants and that lean cases may have a stronger genetic predisposition to type 2 diabetes.

Genetic associations with lipoprotein subfractions provide information on their biological nature.

Adverse levels of lipoproteins are highly heritable and constitute risk factors for cardiovascular outcomes. Hitherto, genome-wide association studies revealed 95 lipid-associated loci. However, due to the small effect sizes of these associations large sample numbers (>100 000 samples) were needed. Here we show that analyzing more refined lipid phenotypes, namely lipoprotein subfractions, can increase the number of significantly associated loci compared with bulk high-density lipoprotein and low-density lipoprotein analysis in a study with identical sample numbers. Moreover, lipoprotein subfractions provide novel insight into the human lipid metabolism. We measured 15 lipoprotein subfractions (L1-L15) in 1791 samples using (1)H-NMR (nuclear magnetic resonance) spectroscopy. Using cluster analyses, we quantified inter-relationships among lipoprotein subfractions. Additionally, we analyzed associations with subfractions at known lipid loci. We identified five distinct groups of subfractions: one (L1) was only marginally captured by serum lipids and therefore extends our knowledge of lipoprotein biochemistry. During a lipid-tolerance test, L1 lost its special position. In the association analysis, we found that eight loci (LIPC, CETP, PLTP, FADS1-2-3, SORT1, GCKR, APOB, APOA1) were associated with the subfractions, whereas only four loci (CETP, SORT1, GCKR, APOA1) were associated with serum lipids. For LIPC, we observed a 10-fold increase in the variance explained by our regression models. In conclusion, NMR-based fine mapping of lipoprotein subfractions provides novel information on their biological nature and strengthens the associations with genetic loci. Future clinical studies are now needed to investigate their biomedical relevance.

A genome-wide screen for interactions reveals a new locus on 4p15 modifying the effect of waist-to-hip ratio on total cholesterol.

Recent genome-wide association (GWA) studies described 95 loci controlling serum lipid levels. These common variants explain ∼25% of the heritability of the phenotypes. To date, no unbiased screen for gene-environment interactions for circulating lipids has been reported. We screened for variants that modify the relationship between known epidemiological risk factors and circulating lipid levels in a meta-analysis of genome-wide association (GWA) data from 18 population-based cohorts with European ancestry (maximum N = 32,225). We collected 8 further cohorts (N = 17,102) for replication, and rs6448771 on 4p15 demonstrated genome-wide significant interaction with waist-to-hip-ratio (WHR) on total cholesterol (TC) with a combined P-value of 4.79×10(-9). There were two potential candidate genes in the region, PCDH7 and CCKAR, with differential expression levels for rs6448771 genotypes in adipose tissue. The effect of WHR on TC was strongest for individuals carrying two copies of G allele, for whom a one standard deviation (sd) difference in WHR corresponds to 0.19 sd difference in TC concentration, while for A allele homozygous the difference was 0.12 sd. Our findings may open up possibilities for targeted intervention strategies for people characterized by specific genomic profiles. However, more refined measures of both body-fat distribution and metabolic measures are needed to understand how their joint dynamics are modified by the newly found locus.

Genetic determinants of serum testosterone concentrations in men.

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterone's high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10(-41) and rs6258, p = 2.3×10(-22)). Subjects with ≥ 3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10(-16)). The rs6258 polymorphism in exon 4 of SHBG affected SHBG's affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.

A genome-wide association study identifies novel loci associated with circulating IGF-I and IGFBP-3.

Insulin-like growth factor-I (IGF-I) and insulin-like growth factor-binding protein-3 (IGFBP-3) are involved in cell replication, proliferation, differentiation, protein synthesis, carbohydrate homeostasis and bone metabolism. Circulating IGF-I and IGFBP-3 concentrations predict anthropometric traits and risk of cancer and cardiovascular disease. In a genome-wide association study of 10 280 middle-aged and older men and women from four community-based cohort studies, we confirmed a known association of single nucleotide polymorphisms in the IGFBP3 gene region on chromosome 7p12.3 with IGFBP-3 concentrations using a significance threshold of P < 5 × 10(-8) (P = 3.3 × 10(-101)). Furthermore, the same IGFBP3 gene locus (e.g. rs11977526) that was associated with IGFBP-3 concentrations was also associated with the opposite direction of effect, with IGF-I concentration after adjustment for IGFBP-3 concentration (P = 1.9 × 10(-26)). A novel and independent locus on chromosome 7p12.3 (rs700752) had genome-wide significant associations with higher IGFBP-3 (P = 4.4 × 10(-21)) and higher IGF-I (P = 4.9 × 10(-9)) concentrations; when the two measurements were adjusted for one another, the IGF-I association was attenuated but the IGFBP-3 association was not. Two additional loci demonstrated genome-wide significant associations with IGFBP-3 concentration (rs1065656, chromosome 16p13.3, P = 1.2 × 10(-11), IGFALS, a confirmatory finding; and rs4234798, chromosome 4p16.1, P = 4.5 × 10(-10), SORCS2, a novel finding). Together, the four genome-wide significant loci explained 6.5% of the population variation in IGFBP-3 concentration. Furthermore, we observed a borderline statistically significant association between IGF-I concentration and FOXO3 (rs2153960, chromosome 6q21, P = 5.1 × 10(-7)), a locus associated with longevity. These genetic loci deserve further investigation to elucidate the biological basis for the observed associations and clarify their possible role in IGF-mediated regulation of cell growth and metabolism.

On the hypothesis-free testing of metabolite ratios in genome-wide and metabolome-wide association studies.

BACKGROUND: Genome-wide association studies (GWAS) with metabolic traits and metabolome-wide association studies (MWAS) with traits of biomedical relevance are powerful tools to identify the contribution of genetic, environmental and lifestyle factors to the etiology of complex diseases. Hypothesis-free testing of ratios between all possible metabolite pairs in GWAS and MWAS has proven to be an innovative approach in the discovery of new biologically meaningful associations. The p-gain statistic was introduced as an ad-hoc measure to determine whether a ratio between two metabolite concentrations carries more information than the two corresponding metabolite concentrations alone. So far, only a rule of thumb was applied to determine the significance of the p-gain. RESULTS: Here we explore the statistical properties of the p-gain through simulation of its density and by sampling of experimental data. We derive critical values of the p-gain for different levels of correlation between metabolite pairs and show that B/(2*α) is a conservative critical value for the p-gain, where α is the level of significance and B the number of tested metabolite pairs. CONCLUSIONS: We show that the p-gain is a well defined measure that can be used to identify statistically significant metabolite ratios in association studies and provide a conservative significance cut-off for the p-gain for use in future association studies with metabolic traits.

Common genetic variants associate with serum phosphorus concentration.

Phosphorus is an essential mineral that maintains cellular energy and mineralizes the skeleton. Because complex actions of ion transporters and regulatory hormones regulate serum phosphorus concentrations, genetic variation may determine interindividual variation in phosphorus metabolism. Here, we report a comprehensive genome-wide association study of serum phosphorus concentration. We evaluated 16,264 participants of European ancestry from the Cardiovascular Heath Study, Atherosclerosis Risk in Communities Study, Framingham Offspring Study, and the Rotterdam Study. We excluded participants with an estimated GFR <45 ml/min per 1.73 m(2) to focus on phosphorus metabolism under normal conditions. We imputed genotypes to approximately 2.5 million single-nucleotide polymorphisms in the HapMap and combined study-specific findings using meta-analysis. We tested top polymorphisms from discovery cohorts in a 5444-person replication sample. Polymorphisms in seven loci with minor allele frequencies 0.08 to 0.49 associate with serum phosphorus concentration (P = 3.5 x 10(-16) to 3.6 x 10(-7)). Three loci were near genes encoding the kidney-specific type IIa sodium phosphate co-transporter (SLC34A1), the calcium-sensing receptor (CASR), and fibroblast growth factor 23 (FGF23), proteins that contribute to phosphorus metabolism. We also identified genes encoding phosphatases, kinases, and phosphodiesterases that have yet-undetermined roles in phosphorus homeostasis. In the replication sample, five of seven top polymorphisms associate with serum phosphorous concentrations (P < 0.05 for each). In conclusion, common genetic variants associate with serum phosphorus in the general population. Further study of the loci identified in this study may help elucidate mechanisms of phosphorus regulation.

The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals.

To dissect the genetic architecture of blood pressure and assess effects on target organ damage, we analyzed 128,272 SNPs from targeted and genome-wide arrays in 201,529 individuals of European ancestry, and genotypes from an additional 140,886 individuals were used for validation. We identified 66 blood pressure-associated loci, of which 17 were new; 15 harbored multiple distinct association signals. The 66 index SNPs were enriched for cis-regulatory elements, particularly in vascular endothelial cells, consistent with a primary role in blood pressure control through modulation of vascular tone across multiple tissues. The 66 index SNPs combined in a risk score showed comparable effects in 64,421 individuals of non-European descent. The 66-SNP blood pressure risk score was significantly associated with target organ damage in multiple tissues but with minor effects in the kidney. Our findings expand current knowledge of blood pressure-related pathways and highlight tissues beyond the classical renal system in blood pressure regulation.

Genome-wide meta-analyses of plasma renin activity and concentration reveal association with the kininogen 1 and prekallikrein genes.

BACKGROUND: The renin-angiotensin-aldosterone system (RAAS) is critical for regulation of blood pressure and fluid balance and influences cardiovascular remodeling. Dysregulation of the RAAS contributes to cardiovascular and renal morbidity. The genetic architecture of circulating RAAS components is incompletely understood. METHODS AND RESULTS: We meta-analyzed genome-wide association data for plasma renin activity (n=5275), plasma renin concentrations (n=8014), and circulating aldosterone (n=13289) from ≤4 population-based cohorts of European and European-American ancestry, and assessed replication of the top results in an independent sample (n=6487). Single-nucleotide polymorphisms (SNPs) in 2 independent loci displayed associations with plasma renin activity at genome-wide significance (P<5×10(-8)). A third locus was close to this threshold (rs4253311 in kallikrein B [KLKB1], P=5.5×10(-8)). Two of these loci replicated in an independent sample for both plasma renin and aldosterone concentrations (SNP rs5030062 in kininogen 1 [KNG1]: P=0.001 for plasma renin, P=0.024 for plasma aldosterone concentration; and rs4253311 with P<0.001 for both plasma renin and aldosterone concentration). SNPs in the NEBL gene reached genome-wide significance for plasma renin concentration in the discovery sample (top SNP rs3915911; P=8.81×10(-9)), but did not replicate (P=0.81). No locus reached genome-wide significance for aldosterone. SNPs rs5030062 and rs4253311 were not related to blood pressure or renal traits; in a companion study, variants in the kallikrein B locus were associated with B-type natriuretic peptide concentrations in blacks. CONCLUSIONS: We identified 2 genetic loci (kininogen 1 and kallikrein B) influencing key components of the RAAS, consistent with the close interrelation between the kallikrein-kinin system and the RAAS.

Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels.

Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, P < 1.09 × 10(-9)) associations between single-nucleotide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N = 1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research.

Diastrophic dysplasia sulfate transporter (SLC26A2) is expressed in the adrenal cortex and regulates aldosterone secretion.

Elucidation of the molecular mechanisms leading to autonomous aldosterone secretion is a prerequisite to define potential targets and biomarkers in the context of primary aldosteronism. After a genome-wide association study with subjects from the population-based Cooperative Health Research in the Region of Augsburg F4 survey, we observed a highly significant association (P=6.78×10(-11)) between the aldosterone to renin ratio and a locus at 5q32. Hypothesizing that this locus may contain genes of relevance for the pathogenesis of primary aldosteronism, we investigated solute carrier family 26 member 2 (SLC26A2), a protein with known transport activity for sulfate and other cations. Within murine tissues, adrenal glands showed the highest expression levels for SLC26A2, which was significantly downregulated on in vivo stimulation with angiotensin II and potassium. SLC26A2 expression was found to be significantly lower in aldosterone-producing adenomas in comparison with normal adrenal glands. In adrenocortical NCI-H295R cells, specific knockdown of SLC26A2 resulted in a highly significant increase in aldosterone secretion. Concomitantly, expression of steroidogenic enzymes, as well as upstream effectors including transcription factors such as NR4A1, CAMK1, and intracellular Ca(2+) content, was upregulated in knockdown cells. To substantiate further these findings in an SLC26A2 mutant mouse model, aldosterone output proved to be increased in a sex-specific manner. In summary, these findings point toward a possible effect of SLC26A2 in the regulation of aldosterone secretion potentially involved in the pathogenesis of primary aldosteronism.

Associations between thyroid hormones and serum metabolite profiles in an euthyroid population.

The aim was to characterise associations between circulating thyroid hormones-free thyroxine (FT4) and thyrotropin (TSH)-and the metabolite profiles in serum samples from participants of the German population-based KORA F4 study. Analyses were based on the metabolite profile of 1463 euthyroid subjects. In serum samples, obtained after overnight fasting (≥8), 151 different metabolites were quantified in a targeted approach including amino acids, acylcarnitines (ACs), and phosphatidylcholines (PCs). Associations between metabolites and thyroid hormone concentrations were analysed using adjusted linear regression models. To draw conclusions on thyroid hormone related pathways, intra-class metabolite ratios were additionally explored. We discovered 154 significant associations (Bonferroni p < 1.75 × 10-04) between FT4 and various metabolites and metabolite ratios belonging to AC and PC groups. Significant associations with TSH were lacking. High FT4 levels were associated with increased concentrations of many ACs and various sums of ACs of different chain length, and the ratio of C2 by C0. The inverse associations observed between FT4 and many serum PCs reflected the general decrease in PC concentrations. Similar results were found in subgroup analyses, e.g., in weight-stable subjects or in obese subjects. Further, results were independent of different parameters for liver or kidney function, or inflammation, which supports the notion of an independent FT4 effect. In fasting euthyroid adults, higher serum FT4 levels are associated with increased serum AC concentrations and an increased ratio of C2 by C0 which is indicative of an overall enhanced fatty acyl mitochondrial transport and ß-oxidation of fatty acids.

Interrogating causal pathways linking genetic variants, small molecule metabolites, and circulating lipids.

BACKGROUND: Emerging technologies based on mass spectrometry or nuclear magnetic resonance enable the monitoring of hundreds of small metabolites from tissues or body fluids. Profiling of metabolites can help elucidate causal pathways linking established genetic variants to known disease risk factors such as blood lipid traits. METHODS: We applied statistical methodology to dissect causal relationships between single nucleotide polymorphisms, metabolite concentrations, and serum lipid traits, focusing on 95 genetic loci reproducibly associated with the four main serum lipids (total-, low-density lipoprotein-, and high-density lipoprotein- cholesterol and triglycerides). The dataset used included 2,973 individuals from two independent population-based cohorts with data for 151 small molecule metabolites and four main serum lipids. Three statistical approaches, namely conditional analysis, Mendelian randomization, and structural equation modeling, were compared to investigate causal relationship at sets of a single nucleotide polymorphism, a metabolite, and a lipid trait associated with one another. RESULTS: A subset of three lipid-associated loci (FADS1, GCKR, and LPA) have a statistically significant association with at least one main lipid and one metabolite concentration in our data, defining a total of 38 cross-associated sets of a single nucleotide polymorphism, a metabolite and a lipid trait. Structural equation modeling provided sufficient discrimination to indicate that the association of a single nucleotide polymorphism with a lipid trait was mediated through a metabolite at 15 of the 38 sets, and involving variants at the FADS1 and GCKR loci. CONCLUSIONS: These data provide a framework for evaluating the causal role of components of the metabolome (or other intermediate factors) in mediating the association between established genetic variants and diseases or traits.

An atlas of genetic influences on human blood metabolites.

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

Identification and MS-assisted interpretation of genetically influenced NMR signals in human plasma.

Nuclear magnetic resonance spectroscopy (NMR) provides robust readouts of many metabolic parameters in one experiment. However, identification of clinically relevant markers in (1)H NMR spectra is a major challenge. Association of NMR-derived quantities with genetic variants can uncover biologically relevant metabolic traits. Using NMR data of plasma samples from 1,757 individuals from the KORA study together with 655,658 genetic variants, we show that ratios between NMR intensities at two chemical shift positions can provide informative and robust biomarkers. We report seven loci of genetic association with NMR-derived traits (APOA1, CETP, CPS1, GCKR, FADS1, LIPC, PYROXD2) and characterize these traits biochemically using mass spectrometry. These ratios may now be used in clinical studies.

Biomarkers for type 2 diabetes and impaired fasting glucose using a nontargeted metabolomics approach.

Using a nontargeted metabolomics approach of 447 fasting plasma metabolites, we searched for novel molecular markers that arise before and after hyperglycemia in a large population-based cohort of 2,204 females (115 type 2 diabetic [T2D] case subjects, 192 individuals with impaired fasting glucose [IFG], and 1,897 control subjects) from TwinsUK. Forty-two metabolites from three major fuel sources (carbohydrates, lipids, and proteins) were found to significantly correlate with T2D after adjusting for multiple testing; of these, 22 were previously reported as associated with T2D or insulin resistance. Fourteen metabolites were found to be associated with IFG. Among the metabolites identified, the branched-chain keto-acid metabolite 3-methyl-2-oxovalerate was the strongest predictive biomarker for IFG after glucose (odds ratio [OR] 1.65 [95% CI 1.39-1.95], P = 8.46 × 10(-9)) and was moderately heritable (h(2) = 0.20). The association was replicated in an independent population (n = 720, OR 1.68 [ 1.34-2.11], P = 6.52 × 10(-6)) and validated in 189 twins with urine metabolomics taken at the same time as plasma (OR 1.87 [1.27-2.75], P = 1 × 10(-3)). Results confirm an important role for catabolism of branched-chain amino acids in T2D and IFG. In conclusion, this T2D-IFG biomarker study has surveyed the broadest panel of nontargeted metabolites to date, revealing both novel and known associated metabolites and providing potential novel targets for clinical prediction and a deeper understanding of causal mechanisms.

Metabolomic markers reveal novel pathways of ageing and early development in human populations.

BACKGROUND: Human ageing is a complex, multifactorial process and early developmental factors affect health outcomes in old age. METHODS: Metabolomic profiling on fasting blood was carried out in 6055 individuals from the UK. Stepwise regression was performed to identify a panel of independent metabolites which could be used as a surrogate for age. We also investigated the association with birthweight overall and within identical discordant twins and with genome-wide methylation levels. RESULTS: We identified a panel of 22 metabolites which combined are strongly correlated with age (R(2) = 59%) and with age-related clinical traits independently of age. One particular metabolite, C-glycosyl tryptophan (C-glyTrp), correlated strongly with age (beta = 0.03, SE = 0.001, P = 7.0 × 10(-157)) and lung function (FEV1 beta = -0.04, SE = 0.008, P = 1.8 × 10(-8) adjusted for age and confounders) and was replicated in an independent population (n = 887). C-glyTrp was also associated with bone mineral density (beta = -0.01, SE = 0.002, P = 1.9 × 10(-6)) and birthweight (beta = -0.06, SE = 0.01, P = 2.5 × 10(-9)). The difference in C-glyTrp levels explained 9.4% of the variance in the difference in birthweight between monozygotic twins. An epigenome-wide association study in 172 individuals identified three CpG-sites, associated with levels of C-glyTrp (P < 2 × 10(-6)). We replicated one CpG site in the promoter of the WDR85 gene in an independent sample of 350 individuals (beta = -0.20, SE = 0.04, P = 2.9 × 10(-8)). WDR85 is a regulator of translation elongation factor 2, essential for protein synthesis in eukaryotes. CONCLUSIONS: Our data illustrate how metabolomic profiling linked with epigenetic studies can identify some key molecular mechanisms potentially determined in early development that produce long-term physiological changes influencing human health and ageing.