NAFLD risk alleles in PNPLA3, TM6SF2, GCKR and LYPLAL1 show divergent metabolic effects.

Fatty liver has been associated with unfavourable metabolic changes in circulation. To provide insights in fatty liver-related metabolic deviations, we compared metabolic association profile of fatty liver versus metabolic association profiles of genotypes increasing the risk of non-alcoholic fatty liver disease (NAFLD). The cross-sectional associations of ultrasound-ascertained fatty liver with 123 metabolic measures were determined in 1810 (Nfatty liver = 338) individuals aged 34-49 years from The Cardiovascular Risk in Young Finns Study. The association profiles of NAFLD-risk alleles in PNPLA3, TM6SF2, GCKR, and LYPLAL1 with the corresponding metabolic measures were obtained from a publicly available metabolomics GWAS including up to 24 925 Europeans. The risk alleles showed different metabolic effects: PNPLA3 rs738409-G, the strongest genetic NAFLD risk factor, did not associate with metabolic changes. Metabolic effects of GCKR rs1260326-T were comparable in many respects to the fatty liver associations. Metabolic effects of LYPLAL1 rs12137855-C were similar, but statistically less robust, to the effects of GCKR rs1260326-T. TM6SF2 rs58542926-T displayed opposite metabolic effects when compared with the fatty liver associations. The metabolic effects of the risk alleles highlight heterogeneity of the molecular pathways leading to fatty liver and suggest that the fatty liver-related changes in the circulating lipids and metabolites may vary depending on the underlying pathophysiological mechanism. Despite the robust cross-sectional associations on population level, the present results showing neutral or cardioprotective metabolic effects for some of the NAFLD risk alleles advocate that hepatic lipid accumulation by itself may not increase the level of circulating lipids or other metabolites.

Identification of seven novel loci associated with amino acid levels using single-variant and gene-based tests in 8545 Finnish men from the METSIM study.

Comprehensive metabolite profiling captures many highly heritable traits, including amino acid levels, which are potentially sensitive biomarkers for disease pathogenesis. To better understand the contribution of genetic variation to amino acid levels, we performed single variant and gene-based tests of association between nine serum amino acids (alanine, glutamine, glycine, histidine, isoleucine, leucine, phenylalanine, tyrosine, and valine) and 16.6 million genotyped and imputed variants in 8545 non-diabetic Finnish men from the METabolic Syndrome In Men (METSIM) study with replication in Northern Finland Birth Cohort (NFBC1966). We identified five novel loci associated with amino acid levels (P = < 5×10-8): LOC157273/PPP1R3B with glycine (rs9987289, P = 2.3×10-26); ZFHX3 (chr16:73326579, minor allele frequency (MAF) = 0.42%, P = 3.6×10-9), LIPC (rs10468017, P = 1.5×10-8), and WWOX (rs9937914, P = 3.8×10-8) with alanine; and TRIB1 with tyrosine (rs28601761, P = 8×10-9). Gene-based tests identified two novel genes harboring missense variants of MAF <1% that show aggregate association with amino acid levels: PYCR1 with glycine (Pgene = 1.5×10-6) and BCAT2 with valine (Pgene = 7.4×10-7); neither gene was implicated by single variant association tests. These findings are among the first applications of gene-based tests to identify new loci for amino acid levels. In addition to the seven novel gene associations, we identified five independent signals at established amino acid loci, including two rare variant signals at GLDC (rs138640017, MAF=0.95%, Pconditional = 5.8×10-40) with glycine levels and HAL (rs141635447, MAF = 0.46%, Pconditional = 9.4×10-11) with histidine levels. Examination of all single variant association results in our data revealed a strong inverse relationship between effect size and MAF (Ptrend<0.001). These novel signals provide further insight into the molecular mechanisms of amino acid metabolism and potentially, their perturbations in disease.

Common, low-frequency, and rare genetic variants associated with lipoprotein subclasses and triglyceride measures in Finnish men from the METSIM study.

Lipid and lipoprotein subclasses are associated with metabolic and cardiovascular diseases, yet the genetic contributions to variability in subclass traits are not fully understood. We conducted single-variant and gene-based association tests between 15.1M variants from genome-wide and exome array and imputed genotypes and 72 lipid and lipoprotein traits in 8,372 Finns. After accounting for 885 variants at 157 previously identified lipid loci, we identified five novel signals near established loci at HIF3A, ADAMTS3, PLTP, LCAT, and LIPG. Four of the signals were identified with a low-frequency (0.005

Total liver phosphatidylcholine content associates with non-alcoholic steatohepatitis and glycine N-methyltransferase expression.

BACKGROUND & AIMS: Alterations in liver phosphatidylcholine (PC) metabolism have been implicated in the pathogenesis of non-alcoholic fatty liver disease (NAFLD). Although genetic variation in the phosphatidylethanolamine N-methyltransferase (PEMT) enzyme synthesizing PC has been associated with disease, the functional mechanism linking PC metabolism to the pathogenesis of non-alcoholic steatohepatitis (NASH) remains unclear. METHODS: Serum PC levels and liver PC contents were measured using proton nuclear magnetic resonance (NMR) spectroscopy in 169 obese individuals [age 46.6 ± 10 (mean ± SD) years, BMI 43.3 ± 6 kg/m2 , 53 men and 116 women] with histological assessment of NAFLD; 106 of these had a distinct liver phenotype. All subjects were genotyped for PEMT rs7946 and liver mRNA expression of PEMT and glycine N-methyltransferase (GNMT) was analysed. RESULTS: Liver PC content was lower in those with NASH (P = 1.8 x 10-6 ) while serum PC levels did not differ between individuals with NASH and normal liver (P = 0.591). Interestingly, serum and liver PC did not correlate (rs  = -0.047, P = 0.557). Serum PC and serum cholesterol levels correlated strongly (rs  = 0.866, P = 7.1 x 10-49 ), while liver PC content did not correlate with serum cholesterol (rs  = 0.065, P = 0.413). Neither PEMT V175M genotype nor PEMT expression explained the association between liver PC content and NASH. Instead, liver GNMT mRNA expression was decreased in those with NASH (P = 3.8 x 10-4 ) and correlated with liver PC content (rs  = 0.265, P = 0.001). CONCLUSIONS: Decreased liver PC content in individuals with the NASH is independent of PEMT V175M genotype and could be partly linked to decreased GNMT expression.

Multiple Hepatic Regulatory Variants at the GALNT2 GWAS Locus Associated with High-Density Lipoprotein Cholesterol.

Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

Effect of Dietary Counseling on a Comprehensive Metabolic Profile from Childhood to Adulthood.

OBJECTIVES: To study the effects of repeated, infancy-onset dietary counseling on a detailed metabolic profile. Effects of dietary saturated fat replacement on circulating concentrations of metabolic biomarkers still remain unknown. STUDY DESIGN: The Special Turku Coronary Risk Factor Intervention Project (STRIP) study is a longitudinal, randomized atherosclerosis prevention trial in which repeated dietary counseling aimed at reducing the proportion of saturated fat intake. Nuclear magnetic resonance metabolomics quantified circulating metabolites from serum samples assessed at age 9 (n = 554), 11 (n = 553), 13 (n = 508), 15 (n = 517), 17 (n = 457), and 19 (n = 417) years. RESULTS: The intervention reduced dietary intake of saturated fat (mean difference in daily percentage of total energy intake: -2.1 [95% CI -1.9 to -2.3]) and increased intake of polyunsaturated fat (0.6 [0.5-0.7]). The dietary counseling intervention led to greater serum proportions of polyunsaturated fatty acids (P < .001), with greater proportions of both circulating omega-3 (P = .02) and omega-6 (P < .001) fatty acids. The proportion of saturated fatty acids in serum was lower for both boys and girls in the intervention group (P < .001), whereas the serum proportion of monounsaturated fat was lower for boys in the intervention group only (P < .001). The intervention also reduced circulating intermediate-density lipoprotein and low-density lipoprotein lipid concentrations (P < .01). Dietary intervention effects on nonlipid biomarkers were minor except from greater concentrations of glutamine in the intervention group. CONCLUSIONS: Repeated dietary counseling from infancy to early adulthood yielded favorable effects on multiple circulating fatty acids and lipoprotein subclass lipids, particularly in boys. These molecular effects substantiate the beneficial role of saturated fat replacement on the metabolic risk profile. TRIAL REGISTRATION: ClinicalTrials.gov: NCT00223600.

Metabolic profiling of fatty liver in young and middle-aged adults: Cross-sectional and prospective analyses of the Young Finns Study.

Nonalcoholic fatty liver is associated with obesity-related metabolic disturbances, but little is known about the metabolic perturbations preceding fatty liver disease. We performed comprehensive metabolic profiling to assess how circulating metabolites, such as lipoprotein lipids, fatty acids, amino acids, and glycolysis-related metabolites, reflect the presence of and future risk for fatty liver in young adults. Sixty-eight lipids and metabolites were quantified by nuclear magnetic resonance metabolomics in the population-based Young Finns Study from serum collected in 2001 (n = 1,575), 2007 (n = 1,509), and 2011 (n = 2,002). Fatty liver was diagnosed by ultrasound in 2011 when participants were aged 34-49 years (19% prevalence). Cross-sectional associations as well as 4-year and 10-year risks for fatty liver were assessed by logistic regression. Metabolites across multiple pathways were strongly associated with the presence of fatty liver (P < 0.0007 for 60 measures in age-adjusted and sex-adjusted cross-sectional analyses). The strongest direct associations were observed for extremely large very-low-density lipoprotein triglycerides (odds ratio [OR] = 4.86 per 1 standard deviation, 95% confidence interval 3.48-6.78), other very-low-density lipoprotein measures, and branched-chain amino acids (e.g., leucine OR = 2.94, 2.51-3.44). Strong inverse associations were observed for high-density lipoprotein measures, e.g., high-density lipoprotein size (OR = 0.36, 0.30-0.42) and several fatty acids including omega-6 (OR = 0.37, 0.32-0.42). The metabolic associations were attenuated but remained significant after adjusting for waist, physical activity, alcohol consumption, and smoking (P < 0.0007). Similar aberrations in the metabolic profile were observed already 10 years before fatty liver diagnosis. CONCLUSION: Circulating lipids, fatty acids, and amino acids reflect fatty liver independently of routine metabolic risk factors; these metabolic aberrations appear to precede the development of fatty liver in young adults. (Hepatology 2017;65:491-500).

Association of branched-chain amino acids and other circulating metabolites with risk of incident dementia and Alzheimer's disease: A prospective study in eight cohorts.

INTRODUCTION: Metabolite, lipid, and lipoprotein lipid profiling can provide novel insights into mechanisms underlying incident dementia and Alzheimer's disease. METHODS: We studied eight prospective cohorts with 22,623 participants profiled by nuclear magnetic resonance or mass spectrometry metabolomics. Four cohorts were used for discovery with replication undertaken in the other four to avoid false positives. For metabolites that survived replication, combined association results are presented. RESULTS: Over 246,698 person-years, 995 and 745 cases of incident dementia and Alzheimer's disease were detected, respectively. Three branched-chain amino acids (isoleucine, leucine, and valine), creatinine and two very low density lipoprotein (VLDL)-specific lipoprotein lipid subclasses were associated with lower dementia risk. One high density lipoprotein (HDL; the concentration of cholesterol esters relative to total lipids in large HDL) and one VLDL (total cholesterol to total lipids ratio in very large VLDL) lipoprotein lipid subclass was associated with increased dementia risk. Branched-chain amino acids were also associated with decreased Alzheimer's disease risk and the concentration of cholesterol esters relative to total lipids in large HDL with increased Alzheimer's disease risk. DISCUSSION: Further studies can clarify whether these molecules play a causal role in dementia pathogenesis or are merely markers of early pathology.

Novel association of TM6SF2 rs58542926 genotype with increased serum tyrosine levels and decreased apoB-100 particles in Finns.

A glutamate-to-lysine variant (rs58542926-T) in transmembrane 6 superfamily member 2 (TM6SF2) is associated with increased fatty liver disease and diabetes in conjunction with decreased cardiovascular disease risk. To identify mediators of the effects of TM6SF2, we tested for associations between rs58542926-T and serum lipoprotein/metabolite measures in cross-sectional data from nondiabetic statin-naïve participants. We identified independent associations between rs58542926-T and apoB-100 particles (ß = -0.057 g/l, P = 1.99 × 10-14) and tyrosine levels (ß = 0.0020 mmol/l, P = 1.10 × 10-8), controlling for potential confounders, in 6,929 Finnish men. The association between rs58542926-T and apoB-100 was confirmed in an independent sample of 2,196 Finnish individuals from the FINRISK study (ßreplication = -0.029, Preplication = 0.029). Secondary analyses demonstrated an rs58542926-T dose-dependent decrease in particle concentration, cholesterol, and triglyceride (TG) content for VLDL and LDL particles (P < 0.001 for all). No significant associations between rs58542926-T and HDL measures were observed. TM6SF2 SNP rs58542926-T and tyrosine levels were associated with increased incident T2D risk in both METSIM and FINRISK. Decreased liver production/secretion of VLDL, decreased cholesterol and TGs in VLDL/LDL particles in serum, and increased tyrosine levels identify possible mechanisms by which rs58542926-T exerts its effects on increasing risk of fatty liver disease, decreasing cardiovascular disease, and increasing diabetes risk, respectively.

Association of pre-pregnancy body mass index with offspring metabolic profile: Analyses of 3 European prospective birth cohorts.

BACKGROUND: A high proportion of women start pregnancy overweight or obese. According to the developmental overnutrition hypothesis, this could lead offspring to have metabolic disruption throughout their lives and thus perpetuate the obesity epidemic across generations. Concerns about this hypothesis are influencing antenatal care. However, it is unknown whether maternal pregnancy adiposity is associated with long-term risk of adverse metabolic profiles in offspring, and if so, whether this association is causal, via intrauterine mechanisms, or explained by shared familial (genetic, lifestyle, socioeconomic) characteristics. We aimed to determine if associations between maternal body mass index (BMI) and offspring systemic cardio-metabolic profile are causal, via intrauterine mechanisms, or due to shared familial factors. METHODS AND FINDINGS: We used 1- and 2-stage individual participant data (IPD) meta-analysis, and a negative-control (paternal BMI) to examine the association between maternal pre-pregnancy BMI and offspring serum metabolome from 3 European birth cohorts (offspring age at blood collection: 16, 17, and 31 years). Circulating metabolic traits were quantified by high-throughput nuclear magnetic resonance metabolomics. Results from 1-stage IPD meta-analysis (N = 5327 to 5377 mother-father-offspring trios) showed that increasing maternal and paternal BMI was associated with an adverse cardio-metabolic profile in offspring. We observed strong positive associations with very-low-density lipoprotein (VLDL)-lipoproteins, VLDL-cholesterol (C), VLDL-triglycerides, VLDL-diameter, branched/aromatic amino acids, glycoprotein acetyls, and triglycerides, and strong negative associations with high-density lipoprotein (HDL), HDL-diameter, HDL-C, HDL2-C, and HDL3-C (all P < 0.003). Slightly stronger magnitudes of associations were present for maternal compared with paternal BMI across these associations; however, there was no strong statistical evidence for heterogeneity between them (all bootstrap P > 0.003, equivalent to P > 0.05 after accounting for multiple testing). Results were similar in each individual cohort, and in the 2-stage analysis. Offspring BMI showed similar patterns of cross-sectional association with metabolic profile as for parental pre-pregnancy BMI associations but with greater magnitudes. Adjustment of parental BMI-offspring metabolic traits associations for offspring BMI suggested the parental associations were largely due to the association of parental BMI with offspring BMI. Limitations of this study are that inferences cannot be drawn about the role of circulating maternal fetal fuels (i.e., glucose, lipids, fatty acids, and amino acids) on later offspring metabolic profile. In addition, BMI may not reflect potential effects of maternal pregnancy fat distribution. CONCLUSION: Our findings suggest that maternal BMI-offspring metabolome associations are likely to be largely due to shared genetic or familial lifestyle confounding rather than to intrauterine mechanisms.

Quantitative Serum Nuclear Magnetic Resonance Metabolomics in Large-Scale Epidemiology: A Primer on -Omic Technologies.

Detailed metabolic profiling in large-scale epidemiologic studies has uncovered novel biomarkers for cardiometabolic diseases and clarified the molecular associations of established risk factors. A quantitative metabolomics platform based on nuclear magnetic resonance spectroscopy has found widespread use, already profiling over 400,000 blood samples. Over 200 metabolic measures are quantified per sample; in addition to many biomarkers routinely used in epidemiology, the method simultaneously provides fine-grained lipoprotein subclass profiling and quantification of circulating fatty acids, amino acids, gluconeogenesis-related metabolites, and many other molecules from multiple metabolic pathways. Here we focus on applications of magnetic resonance metabolomics for quantifying circulating biomarkers in large-scale epidemiology. We highlight the molecular characterization of risk factors, use of Mendelian randomization, and the key issues of study design and analyses of metabolic profiling for epidemiology. We also detail how integration of metabolic profiling data with genetics can enhance drug development. We discuss why quantitative metabolic profiling is becoming widespread in epidemiology and biobanking. Although large-scale applications of metabolic profiling are still novel, it seems likely that comprehensive biomarker data will contribute to etiologic understanding of various diseases and abilities to predict disease risks, with the potential to translate into multiple clinical settings.

metaCCA: summary statistics-based multivariate meta-analysis of genome-wide association studies using canonical correlation analysis.

MOTIVATION: A dominant approach to genetic association studies is to perform univariate tests between genotype-phenotype pairs. However, analyzing related traits together increases statistical power, and certain complex associations become detectable only when several variants are tested jointly. Currently, modest sample sizes of individual cohorts, and restricted availability of individual-level genotype-phenotype data across the cohorts limit conducting multivariate tests. RESULTS: We introduce metaCCA, a computational framework for summary statistics-based analysis of a single or multiple studies that allows multivariate representation of both genotype and phenotype. It extends the statistical technique of canonical correlation analysis to the setting where original individual-level records are not available, and employs a covariance shrinkage algorithm to achieve robustness.Multivariate meta-analysis of two Finnish studies of nuclear magnetic resonance metabolomics by metaCCA, using standard univariate output from the program SNPTEST, shows an excellent agreement with the pooled individual-level analysis of original data. Motivated by strong multivariate signals in the lipid genes tested, we envision that multivariate association testing using metaCCA has a great potential to provide novel insights from already published summary statistics from high-throughput phenotyping technologies. AVAILABILITY AND IMPLEMENTATION: Code is available at https://github.com/aalto-ics-kepaco CONTACTS: anna.cichonska@helsinki.fi or matti.pirinen@helsinki.fi SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DHA mediates the protective effect of fish consumption on new episodes of depression among women.

In a longitudinal cohort study of young Australian adults, we reported that for women higher baseline levels of fish consumption were associated with reduced incidence of new depressive episodes during the 5-year follow-up. Fish are high in both n-3 fatty acids and tyrosine. In this study, we seek to determine whether n-3 fatty acids or tyrosine explain the observed association. During 2004-2006, a FFQ (nine fish items) was used to estimate weekly fish consumption among 546 women aged 26-36 years. A fasting blood sample was taken and high-throughput NMR spectroscopy was used to measure 233 metabolites, including serum n-3 fatty acids and tyrosine. During 2009-2011, new episodes of depression since baseline were identified using the lifetime version of the Composite International Diagnostic Interview. Relative risks were calculated using log-binomial regression and indirect effects estimated using the STATA binary_mediation command. Potential mediators were added to separate models, and mediation was quantified as the proportion of the total effect due to the mediator. The n-3 DHA mediated 25·3 % of the association between fish consumption and depression when fish consumption was analysed as a continuous variable and 16·6 % when dichotomised (reference group: <2 serves/week). Tyrosine did not mediate the association (<0·1 %). Components in fish other than n-3 fatty acids and tyrosine might be beneficial for women's mental health.

Metabolite profiling and cardiovascular event risk: a prospective study of 3 population-based cohorts.

BACKGROUND: High-throughput profiling of circulating metabolites may improve cardiovascular risk prediction over established risk factors. METHODS AND RESULTS: We applied quantitative nuclear magnetic resonance metabolomics to identify the biomarkers for incident cardiovascular disease during long-term follow-up. Biomarker discovery was conducted in the National Finnish FINRISK study (n=7256; 800 events). Replication and incremental risk prediction was assessed in the Southall and Brent Revisited (SABRE) study (n=2622; 573 events) and British Women's Health and Heart Study (n=3563; 368 events). In targeted analyses of 68 lipids and metabolites, 33 measures were associated with incident cardiovascular events at P<0.0007 after adjusting for age, sex, blood pressure, smoking, diabetes mellitus, and medication. When further adjusting for routine lipids, 4 metabolites were associated with future cardiovascular events in meta-analyses: higher serum phenylalanine (hazard ratio per standard deviation, 1.18; 95% confidence interval, 1.12-1.24; P=4×10(-10)) and monounsaturated fatty acid levels (1.17; 1.11-1.24; P=1×10(-8)) were associated with increased cardiovascular risk, while higher omega-6 fatty acids (0.89; 0.84-0.94; P=6×10(-5)) and docosahexaenoic acid levels (0.90; 0.86-0.95; P=5×10(-5)) were associated with lower risk. A risk score incorporating these 4 biomarkers was derived in FINRISK. Risk prediction estimates were more accurate in the 2 validation cohorts (relative integrated discrimination improvement, 8.8% and 4.3%), albeit discrimination was not enhanced. Risk classification was particularly improved for persons in the 5% to 10% risk range (net reclassification, 27.1% and 15.5%). Biomarker associations were further corroborated with mass spectrometry in FINRISK (n=671) and the Framingham Offspring Study (n=2289). CONCLUSIONS: Metabolite profiling in large prospective cohorts identified phenylalanine, monounsaturated fatty acids, and polyunsaturated fatty acids as biomarkers for cardiovascular risk. This study substantiates the value of high-throughput metabolomics for biomarker discovery and improved risk assessment.

Metabolic profiling of pregnancy: cross-sectional and longitudinal evidence.

BACKGROUND: Pregnancy triggers well-known alterations in maternal glucose and lipid balance but its overall effects on systemic metabolism remain incompletely understood. METHODS: Detailed molecular profiles (87 metabolic measures and 37 cytokines) were measured for up to 4260 women (24-49 years, 322 pregnant) from three population-based cohorts in Finland. Circulating molecular concentrations in pregnant women were compared to those in non-pregnant women. Metabolic profiles were also reassessed for 583 women 6 years later to uncover the longitudinal metabolic changes in response to change in the pregnancy status. RESULTS: Compared to non-pregnant women, all lipoprotein subclasses and lipids were markedly increased in pregnant women. The most pronounced differences were observed for the intermediate-density, low-density and high-density lipoprotein triglyceride concentrations. Large differences were also seen for many fatty acids and amino acids. Pregnant women also had higher concentrations of low-grade inflammatory marker glycoprotein acetyls, higher concentrations of interleukin-18 and lower concentrations of interleukin-12p70. The changes in metabolic concentrations for women who were not pregnant at baseline but pregnant 6 years later (or vice versa) matched (or were mirror-images of) the cross-sectional association pattern. Cross-sectional results were consistent across the three cohorts and similar longitudinal changes were seen for 653 women in 4-year and 497 women in 10-year follow-up. For multiple metabolic measures, the changes increased in magnitude across the three trimesters. CONCLUSIONS: Pregnancy initiates substantial metabolic and inflammatory changes in the mothers. Comprehensive characterisation of normal pregnancy is important for gaining understanding of the key nutrients for fetal growth and development. These findings also provide a valuable molecular reference in relation to studies of adverse pregnancy outcomes.

Variant rs10911021 that associates with coronary heart disease in type 2 diabetes, is associated with lower concentrations of circulating HDL cholesterol and large HDL particles but not with amino acids.

AIMS: An intergenic locus on chromosome 1 (lead SNP rs10911021) was previously associated with coronary heart disease (CHD) in type 2 diabetes (T2D). Using data from the UCLEB consortium we investigated the relationship between rs10911021 and CHD in T2D, whether rs10911021 was associated with levels of amino acids involved in the γ-glutamyl cycle or any conventional risk factors (CRFs) for CHD in the T2D participants. METHODS: Four UCLEB studies (n = 6531) had rs10911021 imputation, CHD in T2D, CRF and metabolomics data determined using a nuclear magnetic resonance based platform. RESULTS: The expected direction of effect between rs10911021 and CHD in T2D was observed (1377 no CHD/160 CHD; minor allele OR 0.80, 95 % CI 0.60-1.06) although this was not statistically significant (p = 0.13). No association between rs10911021 and CHD was seen in non-T2D participants (11218 no CHD/1274 CHD; minor allele OR 1.00 95 % CIs 0.92-1.10). In T2D participants, while no associations were observed between rs10911021 and the nine amino acids measured, rs10911021 was associated with HDL-cholesterol (p = 0.0005) but the minor "protective" allele was associated with lower levels (-0.034 mmol/l per allele). Focusing more closely on the HDL-cholesterol subclasses measured, we observed that rs10911021 was associated with six large HDL particle measures in T2D (all p < 0.001). No significant associations were seen in non-T2D subjects. CONCLUSIONS: Our findings are consistent with a true association between rs10911021 and CHD in T2D. The protective minor allele was associated with lower HDL-cholesterol and reductions in HDL particle traits. Our results indicate a complex relationship between rs10911021 and CHD in T2D.

Multi-omic signature of body weight change: results from a population-based cohort study.

BACKGROUND: Excess body weight is a major risk factor for cardiometabolic diseases. The complex molecular mechanisms of body weight change-induced metabolic perturbations are not fully understood. Specifically, in-depth molecular characterization of long-term body weight change in the general population is lacking. Here, we pursued a multi-omic approach to comprehensively study metabolic consequences of body weight change during a seven-year follow-up in a large prospective study. METHODS: We used data from the population-based Cooperative Health Research in the Region of Augsburg (KORA) S4/F4 cohort. At follow-up (F4), two-platform serum metabolomics and whole blood gene expression measurements were obtained for 1,631 and 689 participants, respectively. Using weighted correlation network analysis, omics data were clustered into modules of closely connected molecules, followed by the formation of a partial correlation network from the modules. Association of the omics modules with previous annual percentage weight change was then determined using linear models. In addition, we performed pathway enrichment analyses, stability analyses, and assessed the relation of the omics modules with clinical traits. RESULTS: Four metabolite and two gene expression modules were significantly and stably associated with body weight change (P-values ranging from 1.9 × 10(-4) to 1.2 × 10(-24)). The four metabolite modules covered major branches of metabolism, with VLDL, LDL and large HDL subclasses, triglycerides, branched-chain amino acids and markers of energy metabolism among the main representative molecules. One gene expression module suggests a role of weight change in red blood cell development. The other gene expression module largely overlaps with the lipid-leukocyte (LL) module previously reported to interact with serum metabolites, for which we identify additional co-expressed genes. The omics modules were interrelated and showed cross-sectional associations with clinical traits. Moreover, weight gain and weight loss showed largely opposing associations with the omics modules. CONCLUSIONS: Long-term weight change in the general population globally associates with serum metabolite concentrations. An integrated metabolomics and transcriptomics approach improved the understanding of molecular mechanisms underlying the association of weight gain with changes in lipid and amino acid metabolism, insulin sensitivity, mitochondrial function as well as blood cell development and function.

Assessing multivariate gene-metabolome associations with rare variants using Bayesian reduced rank regression.

MOTIVATION: A typical genome-wide association study searches for associations between single nucleotide polymorphisms (SNPs) and a univariate phenotype. However, there is a growing interest to investigate associations between genomics data and multivariate phenotypes, for example, in gene expression or metabolomics studies. A common approach is to perform a univariate test between each genotype-phenotype pair, and then to apply a stringent significance cutoff to account for the large number of tests performed. However, this approach has limited ability to uncover dependencies involving multiple variables. Another trend in the current genetics is the investigation of the impact of rare variants on the phenotype, where the standard methods often fail owing to lack of power when the minor allele is present in only a limited number of individuals. RESULTS: We propose a new statistical approach based on Bayesian reduced rank regression to assess the impact of multiple SNPs on a high-dimensional phenotype. Because of the method's ability to combine information over multiple SNPs and phenotypes, it is particularly suitable for detecting associations involving rare variants. We demonstrate the potential of our method and compare it with alternatives using the Northern Finland Birth Cohort with 4702 individuals, for whom genome-wide SNP data along with lipoprotein profiles comprising 74 traits are available. We discovered two genes (XRCC4 and MTHFD2L) without previously reported associations, which replicated in a combined analysis of two additional cohorts: 2390 individuals from the Cardiovascular Risk in Young Finns study and 3659 individuals from the FINRISK study. AVAILABILITY AND IMPLEMENTATION: R-code freely available for download at http://users.ics.aalto.fi/pemartti/gene_metabolome/.

Ketone body production is differentially altered in steatosis and non-alcoholic steatohepatitis in obese humans.

BACKGROUND & AIMS: Levels of ketone bodies have been reported to be both increased and decreased in individuals with non-alcoholic fatty liver disease. We investigated whether the metabolism of ketone bodies is different in simple steatosis and in non-alcoholic steatohepatitis (NASH). METHODS: Serum low molecular weight molecules including ketone bodies were measured using high-throughput proton (1H) nuclear magnetic resonance in 116 (76 categorized unequivocally to those with normal liver, simple steatosis or NASH) morbidly obese individuals [age 47.3 ± 8.7 (mean ± SD) years, body mass index 45.1 ± 6.1 kg/m(2) , 39 men and 77 women] with histological assessment of NASH and analysis of gene expression in the liver. Finally, we correlated ß-hydroxybutyrate (ß-OHB) levels with NASH predicting score in Metabolic Syndrome in Men Study (METSIM) population study (n = 8749 non-diabetic men). RESULTS: Levels of ketone bodies were lower in individuals with NASH compared to individuals with simple steatosis (P = 0.004 and P = 0.018 for ß-OHB and acetoacetate respectively). Lower levels of ß-OHB were associated with the NASH predicting score in the METSIM study (P = 0.001). Liver inflammation correlated with mRNA expression of genes regulating ketolysis in the liver (Spearman correlation 0.379-0.388, P < 0.0006 for ACAT1, ACSS2 and BDH1). CONCLUSION: Lower levels of ketone bodies in individuals with NASH compared to individuals with simple steatosis suggest a decrease in ketone body metabolism in NASH.

Novel Loci for metabolic networks and multi-tissue expression studies reveal genes for atherosclerosis.

Association testing of multiple correlated phenotypes offers better power than univariate analysis of single traits. We analyzed 6,600 individuals from two population-based cohorts with both genome-wide SNP data and serum metabolomic profiles. From the observed correlation structure of 130 metabolites measured by nuclear magnetic resonance, we identified 11 metabolic networks and performed a multivariate genome-wide association analysis. We identified 34 genomic loci at genome-wide significance, of which 7 are novel. In comparison to univariate tests, multivariate association analysis identified nearly twice as many significant associations in total. Multi-tissue gene expression studies identified variants in our top loci, SERPINA1 and AQP9, as eQTLs and showed that SERPINA1 and AQP9 expression in human blood was associated with metabolites from their corresponding metabolic networks. Finally, liver expression of AQP9 was associated with atherosclerotic lesion area in mice, and in human arterial tissue both SERPINA1 and AQP9 were shown to be upregulated (6.3-fold and 4.6-fold, respectively) in atherosclerotic plaques. Our study illustrates the power of multi-phenotype GWAS and highlights candidate genes for atherosclerosis.