Impact of recognition of genetic information related to BMI on changes in physical activity, dietary intake, and blood cholesterol level: a randomized controlled trial

PURPOSE

Genes associated with body mass index (BMI), including FTO rs9939609,MC4R rs17782313, and BDNF rs6265, may influence BMI and regulate energy metabolism. While previous studies have explored health-related behavior changes, few have investigated both biochemical and behavior changes resulting from perceived genetic risk. This study investigated whether recognizing BMI-related genes affects health-related behaviors and alters blood metabolite levels.

METHODS

Normal and overweight adults aged 25-35 years (n = 100) were randomly assigned to an intervention group (n = 65) informed about BMI-related genetic information (FTO rs9939609, MC4R rs17782313, BDNF rs6265) and an uninformed group (n = 35, CON). The intervention group was further divided into Intervention-high risk (IHR, n = 36) and intervention-low risk (ILR, n = 29) subgroups. Dietary intake and physical activity (PA) were assessed using a 3-day dietary record and the IPAQ-short form. Blood metabolites were analyzed through multivariate analyses to identify significant differences among the groups, with measurements taken at baseline, 3 months, and 6 months.

RESULTS

The IHR group exhibited increased dietary fat and fast foods intake, along with enhanced vigorous and moderate PA. Six metabolites were selected as biomarkers that were distinguishable among groups, and the relative serum cholesterol levels significantly decreased in the IHR group at 3 months.

CONCLUSION

These results demonstrate that recognizing the BMI-associated genetic risk resulted in a short-term increase in PA but did not improve dietary intake. Increased PA was significantly associated with reduced cholesterol concentration, suggesting the clinical importance of physical activity in the genetically at-risk group.

CLINICAL TRIAL AND STUDY REGISTRATION

This study was reviewed and approved by the Seoul National University Institutional Review Board (IRB #1901/001-004) and registered on the Clinical Research Information Service (CRIS), KCT0004650 ( https://cris.nih.go.kr/cris/search/detailSearch.do /14091, 2020/01/28).

Serum Metabolome Signature Response to Different Types of Resistance Training

Pneumatic resistance training (PRT) facilitates a longer time under tension that might lead to greater changes in body composition when compared to traditional resistance training (TRT), possibly enhancing serum metabolite concentrations indicative of healthy metabolic function. To assess the impact of PRT and TRT on muscular strength, body composition, and serum metabolome, 69 men (age: 31.8±7.2 years, height: 179.7±5.4 cm, weight: 81.1±9.9 kg) were randomized into two 10-week intervention groups (PRT:n=24 and TRT:n=24) and one control group (CON:n=21). Serum metabolite concentrations were assessed before and after the training intervention by high-throughput nuclear magnetic resonance. Fat mass and lean mass were obtained by bioimpedance analysis. The training intervention resulted in an increase in lean mass for both PRT (1.85±2.69%; p=0.003) and TRT (2.72±4.53%; p=0.004), while only PRT reduced statistically significantly in body fat percentage (PRT: -5.08±10.76%; p=0.019). Only in PRT and TRT significant increases in small high-density lipoproteins (S-HDL-L) and small HDL particles (S-HDL-P) were observed. When controlling for fat and lean mass, the effects on S-HDL-L/S-HDL-P diminished. Network analysis may suggest that PRT and TRT result in an increase in network connectivity and robustness. It appears that the observed improvements are associated with changes in body composition.

Genome-wide characterization of circulating metabolic biomarkers

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

The Relationship of Circulating Choline and Choline-Related Metabolite Levels with Health Outcomes: A Scoping Review of Genome-Wide Association Studies and Mendelian Randomization Studies

Choline is essential for proper liver, muscle, brain, lipid metabolism, cellular membrane composition, and repair. Understanding genetic determinants of circulating choline metabolites can help identify new determinants of choline metabolism, requirements, and their link to disease endpoints. We conducted a scoping review to identify studies assessing the association of genetic polymorphisms on circulating choline and choline-related metabolite concentrations and subsequent associations with health outcomes. This study follows the Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement scoping review extension. Literature was searched to September 28, 2022, in 4 databases: Embase, MEDLINE, Web of Science, and the Biological Science Index. Studies of any duration in humans were considered. Any genome-wide association study (GWAS) investigating genetic variant associations with circulating choline and/or choline-related metabolites and any Mendelian randomization (MR) study investigating the association of genetically predicted circulating choline and/or choline-related metabolites with any health outcome were considered. Qualitative evidence is presented in summary tables. From 1248 total reviewed articles, 53 were included (GWAS = 27; MR = 26). Forty-two circulating choline-related metabolites were tested in association with genetic variants in GWAS studies, primarily trimethylamine N-oxide, betaine, sphingomyelins, lysophosphatidylcholines, and phosphatidylcholines. MR studies investigated associations between 52 total unique choline metabolites and 66 unique health outcomes. Of these, 47 significant associations were reported between 16 metabolites (primarily choline, lysophosphatidylcholines, phosphatidylcholines, betaine, and sphingomyelins) and 27 health outcomes including cancer, cardiovascular, metabolic, bone, and brain-related outcomes. Some articles reported significant associations between multiple choline types and the same health outcome. Genetically predicted circulating choline and choline-related metabolite concentrations are associated with a wide variety of health outcomes. Further research is needed to assess how genetic variability influences choline metabolism and whether individuals with lower genetically predicted circulating choline and choline-related metabolite concentrations would benefit from a dietary intervention or supplementation.

Genome-wide association analysis of plasma lipidome identifies 495 genetic associations

The human plasma lipidome captures risk for cardiometabolic diseases. To discover new lipid-associated variants and understand the link between lipid species and cardiometabolic disorders, we perform univariate and multivariate genome-wide analyses of 179 lipid species in 7174 Finnish individuals. We fine-map the associated loci, prioritize genes, and examine their disease links in 377,277 FinnGen participants. We identify 495 genome-trait associations in 56 genetic loci including 8 novel loci, with a considerable boost provided by the multivariate analysis. For 26 loci, fine-mapping identifies variants with a high causal probability, including 14 coding variants indicating likely causal genes. A phenome-wide analysis across 953 disease endpoints reveals disease associations for 40 lipid loci. For 11 coronary artery disease risk variants, we detect strong associations with lipid species. Our study demonstrates the power of multivariate genetic analysis in correlated lipidomics data and reveals genetic links between diseases and lipid species beyond the standard lipids.

Changes in anxiety and depression levels and meat intake following recognition of low genetic risk for high body mass index, triglycerides, and lipoproteins: A randomized controlled trial

BACKGROUND

Psychological status affects dietary intake, and recognizing genetic information can lead to behavior changes by influencing psychological factors such as anxiety or depression.

OBJECTIVES

In this study, we examined the effects of disclosing genetic information on anxiety or depression levels and the association between these psychological factors and dietary intake.

METHODS

A total of 100 healthy adults were randomly assigned to an intervention group (n = 65) informed about their genetic test results regarding body mass index and lipid profiles (triglyceride and cholesterol concentrations) and a not-informed control group (CON, n = 35). Based on polygenic risk scores, participants in the intervention group were subclassified into an intervention-low risk (ILR, n = 32) and an intervention-high risk (IHR, n = 33) group. Nutrient and food intakes were assessed via a 3-day dietary record at baseline and at 3 and 6 months. Depression and anxiety levels were measured using PHQ-9 and GAD-7 questionnaires, and the relative levels of blood metabolites were measure using GC-MS/MS analysis.

RESULTS

Noticeable changes in dietary intake as well as psychological factors were observed in male subjects, with those perceiving their genetic risks as low (ILR) showing a significant increase in protein intake at 3 months compared to baseline (ILR: 3.9 ± 1.4, p<0.05). Meat intake also increased significantly in males in the ILR group at 3 months, but not in the IHR and CON groups (ILR: 49.4 ± 30.8, IHR: -52.2 ± 25.4, CON: -5.3 ± 30.3 g/d). ILR group showed a significant decrease in anxiety levels at 3 months, and their anxiety scores showed a negative association with meat intake (standardized β = -0.321, p<0.05). The meat intake at 3 months was associated with the relative levels of arginine and ornithine (standardized β = 0.452, p<0.05 and standardized β = 0.474, p<0.05, respectively).

CONCLUSIONS

Taken together, anxiety levels were decreased in male subjects who perceived their genetic risk to be low, and the decrease in anxiety levels was associated with an increase in meat intake. This suggests that recognizing genetic information may affect psychological factors and dietary intake.

Whole-exome sequencing identifies novel protein-altering variants associated with serum apolipoprotein and lipid concentrations

BACKGROUND

Dyslipidemia is a major risk factor for cardiovascular disease, and diabetes impacts the lipid metabolism through multiple pathways. In addition to the standard lipid measurements, apolipoprotein concentrations provide added awareness of the burden of circulating lipoproteins. While common genetic variants modestly affect the serum lipid concentrations, rare genetic mutations can cause monogenic forms of hypercholesterolemia and other genetic disorders of lipid metabolism. We aimed to identify low-frequency protein-altering variants (PAVs) affecting lipoprotein and lipid traits.

METHODS

We analyzed whole-exome (WES) and whole-genome sequencing (WGS) data of 481 and 474 individuals with type 1 diabetes, respectively. The phenotypic data consisted of 79 serum lipid and apolipoprotein phenotypes obtained with clinical laboratory measurements and nuclear magnetic resonance spectroscopy.

RESULTS

The single-variant analysis identified an association between the LIPC p.Thr405Met (rs113298164) and serum apolipoprotein A1 concentrations (p=7.8×10-8). The burden of PAVs was significantly associated with lipid phenotypes in LIPC, RBM47, TRMT5, GTF3C5, MARCHF10, and RYR3 (p<2.9×10-6). The RBM47 gene is required for apolipoprotein B post-translational modifications, and in our data, the association between RBM47 and apolipoprotein C-III concentrations was due to a rare 21 base pair p.Ala496-Ala502 deletion; in replication, the burden of rare deleterious variants in RBM47 was associated with lower triglyceride concentrations in WES of >170,000 individuals from multiple ancestries (p=0.0013). Two PAVs in GTF3C5 were highly enriched in the Finnish population and associated with cardiovascular phenotypes in the general population. In the previously known APOB gene, we identified novel associations at two protein-truncating variants resulting in lower serum non-HDL cholesterol (p=4.8×10-4), apolipoprotein B (p=5.6×10-4), and LDL cholesterol (p=9.5×10-4) concentrations.

CONCLUSIONS

We identified lipid and apolipoprotein-associated variants in the previously known LIPC and APOB genes, as well as PAVs in GTF3C5 associated with LDLC, and in RBM47 associated with apolipoprotein C-III concentrations, implicated as an independent CVD risk factor. Identification of rare loss-of-function variants has previously revealed genes that can be targeted to prevent CVD, such as the LDL cholesterol-lowering loss-of-function variants in the PCSK9 gene. Thus, this study suggests novel putative therapeutic targets for the prevention of CVD.

Whole Genome Association Study of the Plasma Metabolome Identifies Metabolites Linked to Cardiometabolic Disease in Black Individuals

Integrating genetic information with metabolomics has provided new insights into genes affecting human metabolism. However, gene-metabolite integration has been primarily studied in individuals of European Ancestry, limiting the opportunity to leverage genomic diversity for discovery. In addition, these analyses have principally involved known metabolites, with the majority of the profiled peaks left unannotated. Here, we perform a whole genome association study of 2,291 metabolite peaks (known and unknown features) in 2,466 Black individuals from the Jackson Heart Study. We identify 519 locus-metabolite associations for 427 metabolite peaks and validate our findings in two multi-ethnic cohorts. A significant proportion of these associations are in ancestry specific alleles including findings in APOE, TTR and CD36. We leverage tandem mass spectrometry to annotate unknown metabolites, providing new insight into hereditary diseases including transthyretin amyloidosis and sickle cell disease. Our integrative omics approach leverages genomic diversity to provide novel insights into diverse cardiometabolic diseases.

Apt interpretation of comprehensive lipoprotein data in large-scale epidemiology: disclosure of fundamental structural and metabolic relationships

BACKGROUND

Quantitative lipoprotein analytics using nuclear magnetic resonance (NMR) spectroscopy is currently commonplace in large-scale studies. One methodology has become widespread and is currently being utilized also in large biobanks. It allows the comprehensive characterization of 14 lipoprotein subclasses, clinical lipids, apolipoprotein A-I and B. The details of these data are conceptualized here in relation to lipoprotein metabolism with particular attention on the fundamental characteristics of subclass particle numbers, lipid concentrations and compositional measures.

METHODS AND RESULTS

The NMR methodology was applied to fasting serum samples from Northern Finland Birth Cohorts 1966 and 1986 with 5651 and 5605 participants, respectively. All results were highly consistent between the cohorts. Circulating lipid concentrations in a particular lipoprotein subclass arise predominantly as the result of the circulating number of those subclass particles. The spherical lipoprotein particle shape, with a radially oriented surface monolayer, imposes size-dependent biophysical constraints for the lipid composition of individual subclass particles and inherently restricts the accommodation of metabolic changes via compositional modifications. The new finding that the relationship between lipoprotein subclass particle concentrations and the particle size is log-linear reveals that circulating lipoprotein particles are also under rather strict metabolic constraints for both their absolute and relative concentrations.

CONCLUSIONS

The fundamental structural and metabolic relationships between lipoprotein subclasses elucidated in this study empower detailed interpretation of lipoprotein metabolism. Understanding the intricate details of these extensive data is important for the precise interpretation of novel therapeutic opportunities and for fully utilizing the potential of forthcoming analyses of genetic and metabolic data in large biobanks.

Genetic Studies of Metabolomics Change After a Liquid Meal Illuminate Novel Pathways for Glucose and Lipid Metabolism

Humans spend the greater part of the day in a postprandial state. However, the genetic basis of postprandial blood measures is relatively uncharted territory. We examined the genetics of variation in concentrations of postprandial metabolites (t = 150 min) in response to a liquid mixed meal through genome-wide association studies (GWAS) performed in the Netherlands Epidemiology of Obesity (NEO) study (n = 5,705). The metabolite response GWAS identified an association between glucose change and rs10830963:G in the melatonin receptor 1B (β [SE] -0.23 [0.03], P = 2.15 × 10^-19). In addition, the ANKRD55 locus led by rs458741:C showed strong associations with extremely large VLDL (XXLVLDL) particle response (XXLVLDL total cholesterol: β [SE] 0.17 [0.03], P = 5.76 × 10^-10; XXLVLDL cholesterol ester: β [SE] 0.17 [0.03], P = 9.74 × 10^-10), which also revealed strong associations with body composition and diabetes in the UK Biobank (P < 5 × 10^-8). Furthermore, the associations between XXLVLDL response and insulinogenic index, HOMA-β, Matsuda insulin sensitivity index, and HbA_1c in the NEO study implied the role of chylomicron synthesis in diabetes (with false discovery rate-corrected q <0.05). To conclude, genetic studies of metabolomics change after a liquid meal illuminate novel pathways for glucose and lipid metabolism. Further studies are warranted to corroborate biological pathways of the ANKRD55 locus underlying diabetes.

Targeted metabolomic profiling and prediction of cardiovascular events: a prospective study of patients with psoriatic arthritis and psoriasis

OBJECTIVE

In patients with psoriatic disease (PsD), we sought serum metabolites associated with cardiovascular (CV) events and investigated whether they could improve CV risk prediction beyond traditional risk factors and the Framingham Risk Score (FRS).

METHODS

Nuclear magnetic resonance metabolomics identified biomarkers for incident CV events in patients with PsD. The association of each metabolite with incident CV events was analysed using Cox proportional hazards regression models first adjusted for age and sex, and subsequently for traditional CV risk factors. Variable selection was performed using penalisation with boosting after adjusting for age and sex, and the FRS.

RESULTS

Among 977 patients with PsD, 70 patients had incident CV events. In Cox regression models adjusted for CV risk factors, alanine, tyrosine, degree of unsaturation of fatty acids and high-density lipoprotein particles were associated with decreased CV risk. Glycoprotein acetyls, apolipoprotein B and cholesterol remnants were associated with increased CV risk. The age-adjusted and sex-adjusted expanded model with 13 metabolites significantly improved prediction of CV events beyond the model with age and sex alone, with an area under the receiver operator characteristic curve (AUC) of 79.9 versus 72.6, respectively (p=0.02). Compared with the FRS alone (AUC=73.9), the FRS-adjusted expanded model with 11 metabolites (AUC=75.0, p=0.72) did not improve CV risk discrimination.

CONCLUSIONS

We identify novel metabolites associated with the development of CV events in patients with PsD. Further study of their underlying causal role may clarify important pathways leading to CV events in this population.

Direct characterization of cis-regulatory elements and functional dissection of complex genetic associations using HCR-FlowFISH

Effective interpretation of genome function and genetic variation requires a shift from epigenetic mapping of cis-regulatory elements (CREs) to characterization of endogenous function. We developed hybridization chain reaction fluorescence in situ hybridization coupled with flow cytometry (HCR-FlowFISH), a broadly applicable approach to characterize CRISPR-perturbed CREs via accurate quantification of native transcripts, alongside CRISPR activity screen analysis (CASA), a hierarchical Bayesian model to quantify CRE activity. Across >325,000 perturbations, we provide evidence that CREs can regulate multiple genes, skip over the nearest gene and display activating and/or silencing effects. At the cholesterol-level-associated FADS locus, we combine endogenous screens with reporter assays to exhaustively characterize multiple genome-wide association signals, functionally nominate causal variants and, importantly, identify their target genes.

Genomic-Metabolomic Associations Support the Role of LIPC and Glycerophospholipids in Age-Related Macular Degeneration

Purpose

Large-scale genome-wide association studies (GWAS) have reported important single nucleotide polymorphisms (SNPs) with significant associations with age-related macular degeneration (AMD). However, their role in disease development remains elusive. This study aimed to assess SNP–metabolite associations (i.e., metabolite quantitative trait loci [met-QTL]) and to provide insights into the biological mechanisms of AMD risk SNPs.

Design

Cross-sectional multicenter study (Boston, Massachusetts, and Coimbra, Portugal).

Participants

Patients with AMD (n = 388) and control participants (n = 98) without any vitreoretinal disease (> 50 years).

Methods

Age-related macular degeneration grading was performed using color fundus photographs according to the Age-Related Eye Disease Study classification scheme. Fasting blood samples were collected and evaluated with mass spectrometry for metabolomic profiling and Illumina OmniExpress for SNPs profiling. Analyses of met-QTL of endogenous metabolites were conducted using linear regression models adjusted for age, gender, smoking, 10 metabolite principal components (PCs), and 10 SNP PCs. Additionally, we analyzed the cumulative effect of AMD risk SNPs on plasma metabolites by generating genetic risk scores and assessing their associations with metabolites using linear regression models, accounting for the same covariates. Modeling was performed first for each cohort, and then combined by meta-analysis. Multiple comparisons were accounted for using the false discovery rate (FDR).

Main Outcome Measures

Plasma metabolite levels associated with AMD risk SNPs.

Results

After quality control, data for 544 plasma metabolites were included. Meta-analysis of data from all individuals (AMD patients and control participants) identified 28 significant met-QTL (β = 0.016–0.083; FDR q-value < 1.14 × 10^–2), which corresponded to 5 metabolites and 2 genes: ASPM and LIPC. Polymorphisms in the LIPC gene were associated with phosphatidylethanolamine metabolites, which are glycerophospholipids, and polymorphisms in the ASPM gene with branched-chain amino acids. Similar results were observed when considering only patients with AMD. Genetic risk score–metabolite associations further supported a global impact of AMD risk SNPs on the plasma metabolome.

Conclusions

This study demonstrated that genomic–metabolomic associations can provide insights into the biological relevance of AMD risk SNPs. In particular, our results support that the LIPC gene and the glycerophospholipid metabolic pathway may play an important role in AMD, thus offering new potential therapeutic targets for this disease.

Integrating lipidomics and genomics: emerging tools to understand cardiovascular diseases

Cardiovascular diseases (CVDs) are the leading cause of mortality and morbidity worldwide leading to 31% of all global deaths. Early prediction and prevention could greatly reduce the enormous socio-economic burden posed by CVDs. Plasma lipids have been at the center stage of the prediction and prevention strategies for CVDs that have mostly relied on traditional lipids (total cholesterol, total triglycerides, HDL-C and LDL-C). The tremendous advancement in the field of lipidomics in last two decades has facilitated the research efforts to unravel the metabolic dysregulation in CVDs and their genetic determinants, enabling the understanding of pathophysiological mechanisms and identification of predictive biomarkers, beyond traditional lipids. This review presents an overview of the application of lipidomics in epidemiological and genetic studies and their contributions to the current understanding of the field. We review findings of these studies and discuss examples that demonstrates the potential of lipidomics in revealing new biology not captured by traditional lipids and lipoprotein measurements. The promising findings from these studies have raised new opportunities in the fields of personalized and predictive medicine for CVDs. The review further discusses prospects of integrating emerging genomics tools with the high-dimensional lipidome to move forward from the statistical associations towards biological understanding, therapeutic target development and risk prediction. We believe that integrating genomics with lipidome holds a great potential but further advancements in statistical and computational tools are needed to handle the high-dimensional and correlated lipidome.

Genome-Wide Association Study for Serum Omega-3 and Omega-6 Polyunsaturated Fatty Acids: Exploratory Analysis of the Sex-Specific Effects and Dietary Modulation in Mediterranean Subjects with Metabolic Syndrome

Many early studies presented beneficial effects of polyunsaturated fatty acids (PUFA) on cardiovascular risk factors and disease. However, results from recent meta-analyses indicate that this effect would be very low or nil. One of the factors that may contribute to the inconsistency of the results is that, in most studies, genetic factors have not been taken into consideration. It is known that fatty acid desaturase (FADS) gene cluster in chromosome 11 is a very important determinant of plasma PUFA, and that the prevalence of the single nucleotide polymorphisms (SNPs) varies greatly between populations and may constitute a bias in meta-analyses. Previous genome-wide association studies (GWAS) have been carried out in other populations and none of them have investigated sex and Mediterranean dietary pattern interactions at the genome-wide level. Our aims were to undertake a GWAS to discover the genes most associated with serum PUFA concentrations (omega-3, omega-6, and some fatty acids) in a scarcely studied Mediterranean population with metabolic syndrome, and to explore sex and adherence to Mediterranean diet (MedDiet) interactions at the genome-wide level. Serum PUFA were determined by NMR spectroscopy. We found strong robust associations between various SNPs in the FADS cluster and omega-3 concentrations (top-ranked in the adjusted model: FADS1-rs174547, p = 3.34 × 10-14; FADS1-rs174550, p = 5.35 × 10-14; FADS2-rs1535, p = 5.85 × 10-14; FADS1-rs174546, p = 6.72 × 10-14; FADS2-rs174546, p = 9.75 × 10-14; FADS2- rs174576, p = 1.17 × 10-13; FADS2-rs174577, p = 1.12 × 10-12, among others). We also detected a genome-wide significant association with other genes in chromosome 11: MYRF (myelin regulatory factor)-rs174535, p = 1.49 × 10-12; TMEM258 (transmembrane protein 258)-rs102275, p = 2.43 × 10-12; FEN1 (flap structure-specific endonuclease 1)-rs174538, p = 1.96 × 10-11). Similar genome-wide statistically significant results were found for docosahexaenoic fatty acid (DHA). However, no such associations were detected for omega-6 PUFAs or linoleic acid (LA). For total PUFA, we observed a consistent gene*sex interaction with the DNTTIP2 (deoxynucleotidyl transferase terminal interacting protein 2)-rs3747965 p = 1.36 × 10-8. For adherence to MedDiet, we obtained a relevant interaction with the ME1 (malic enzyme 1) gene (a gene strongly regulated by fat) in determining serum omega-3. The top-ranked SNP for this interaction was ME1-rs3798890 (p = 2.15 × 10-7). In the regional-wide association study, specifically focused on the FADS1/FASD2/FADS3 and ELOVL (fatty acid elongase) 2/ELOVL 5 regions, we detected several statistically significant associations at p < 0.05. In conclusion, our results confirm a robust role of the FADS cluster on serum PUFA in this population, but the associations vary depending on the PUFA. Moreover, the detection of some sex and diet interactions underlines the need for these associations/interactions to be studied in all specific populations so as to better understand the complex metabolism of PUFA.

Heritability estimates for 361 blood metabolites across 40 genome-wide association studies

Metabolomics examines the small molecules involved in cellular metabolism. Approximately 50% of total phenotypic differences in metabolite levels is due to genetic variance, but heritability estimates differ across metabolite classes. We perform a review of all genome-wide association and (exome-) sequencing studies published between November 2008 and October 2018, and identify >800 class-specific metabolite loci associated with metabolite levels. In a twin-family cohort (N = 5117), these metabolite loci are leveraged to simultaneously estimate total heritability (h2total), and the proportion of heritability captured by known metabolite loci (h2Metabolite-hits) for 309 lipids and 52 organic acids. Our study reveals significant differences in h2Metabolite-hits among different classes of lipids and organic acids. Furthermore, phosphatidylcholines with a high degree of unsaturation have higher h2Metabolite-hits estimates than phosphatidylcholines with low degrees of unsaturation. This study highlights the importance of common genetic variants for metabolite levels, and elucidates the genetic architecture of metabolite classes.

Lipoprotein signatures of cholesteryl ester transfer protein and HMG-CoA reductase inhibition

Cholesteryl ester transfer protein (CETP) inhibition reduces vascular event risk, but confusion surrounds its effects on low-density lipoprotein (LDL) cholesterol. Here, we clarify associations of genetic inhibition of CETP on detailed lipoprotein measures and compare those to genetic inhibition of 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR). We used an allele associated with lower CETP expression (rs247617) to mimic CETP inhibition and an allele associated with lower HMGCR expression (rs12916) to mimic the well-known effects of statins for comparison. The study consists of 65,427 participants of European ancestries with detailed lipoprotein subclass profiling from nuclear magnetic resonance spectroscopy. Genetic associations were scaled to 10% reduction in relative risk of coronary heart disease (CHD). We also examined observational associations of the lipoprotein subclass measures with risk of incident CHD in 3 population-based cohorts totalling 616 incident cases and 13,564 controls during 8-year follow-up. Genetic inhibition of CETP and HMGCR resulted in near-identical associations with LDL cholesterol concentration estimated by the Friedewald equation. Inhibition of HMGCR had relatively consistent associations on lower cholesterol concentrations across all apolipoprotein B-containing lipoproteins. In contrast, the associations of the inhibition of CETP were stronger on lower remnant and very-low-density lipoprotein (VLDL) cholesterol, but there were no associations on cholesterol concentrations in LDL defined by particle size (diameter 18–26 nm) (−0.02 SD LDL defined by particle size; 95% CI: −0.10 to 0.05 for CETP versus −0.24 SD, 95% CI −0.30 to −0.18 for HMGCR). Inhibition of CETP was strongly associated with lower proportion of triglycerides in all high-density lipoprotein (HDL) particles. In observational analyses, a higher triglyceride composition within HDL subclasses was associated with higher risk of CHD, independently of total cholesterol and triglycerides (strongest hazard ratio per 1 SD higher triglyceride composition in very large HDL 1.35; 95% CI: 1.18–1.54). In conclusion, CETP inhibition does not appear to affect size-specific LDL cholesterol but is likely to lower CHD risk by lowering concentrations of other atherogenic, apolipoprotein B-containing lipoproteins (such as remnant and VLDLs). Inhibition of CETP also lowers triglyceride composition in HDL particles, a phenomenon reflecting combined effects of circulating HDL, triglycerides, and apolipoprotein B-containing particles and is associated with a lower CHD risk in observational analyses. Our results reveal that conventional composite lipid assays may mask heterogeneous effects of emerging lipid-altering therapies.

Inhibition of cholesteryl ester transfer protein does not affect size-specific low-density lipoprotein cholesterol, but may lower coronary heart disease risk by lowering cholesterol concentrations in other apolipoprotein-B containing atherogenic lipoproteins, and by lowering triglyceride content of high-density lipoprotein particles.

Data-driven multivariate population subgrouping via lipoprotein phenotypes versus apolipoprotein B in the risk assessment of coronary heart disease

BACKGROUND AND AIMS

Population subgrouping has been suggested as means to improve coronary heart disease (CHD) risk assessment. We explored here how unsupervised data-driven metabolic subgrouping, based on comprehensive lipoprotein subclass data, would work in large-scale population cohorts.

METHODS

We applied a self-organizing map (SOM) artificial intelligence methodology to define subgroups based on detailed lipoprotein profiles in a population-based cohort (n = 5789) and utilised the trained SOM in an independent cohort (n = 7607). We identified four SOM-based subgroups of individuals with distinct lipoprotein profiles and CHD risk and compared those to univariate subgrouping by apolipoprotein B quartiles.

RESULTS

The SOM-based subgroup with highest concentrations for non-HDL measures had the highest, and the subgroup with lowest concentrations, the lowest risk for CHD. However, apolipoprotein B quartiles produced better resolution of risk than the SOM-based subgroups and also striking dose-response behaviour.

CONCLUSIONS

These results suggest that the majority of lipoprotein-mediated CHD risk is explained by apolipoprotein B-containing lipoprotein particles. Therefore, even advanced multivariate subgrouping, with comprehensive data on lipoprotein metabolism, may not advance CHD risk assessment.

A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context

Genetic studies of metabolites have identified thousands of variants, many of which are associated with downstream metabolic and obesogenic disorders. However, these studies have relied on univariate analyses, reducing power and limiting context-specific understanding. Here we aim to provide an integrated perspective of the genetic basis of metabolites by leveraging the Finnish Metabolic Syndrome In Men (METSIM) cohort, a unique genetic resource which contains metabolic measurements, mostly lipids, across distinct time points as well as information on statin usage. We increase effective sample size by an average of two-fold by applying the Covariates for Multi-phenotype Studies (CMS) approach, identifying 588 significant SNP-metabolite associations, including 228 new associations. Our analysis pinpoints a small number of master metabolic regulator genes, balancing the relative proportion of dozens of metabolite levels. We further identify associations to changes in metabolic levels across time as well as genetic interactions with statin at both the master metabolic regulator and genome-wide level.

A metabolic profile of all-cause mortality risk identified in an observational study of 44,168 individuals

Predicting longer-term mortality risk requires collection of clinical data, which is often cumbersome. Therefore, we use a well-standardized metabolomics platform to identify metabolic predictors of long-term mortality in the circulation of 44,168 individuals (age at baseline 18-109), of whom 5512 died during follow-up. We apply a stepwise (forward-backward) procedure based on meta-analysis results and identify 14 circulating biomarkers independently associating with all-cause mortality. Overall, these associations are similar in men and women and across different age strata. We subsequently show that the prediction accuracy of 5- and 10-year mortality based on a model containing the identified biomarkers and sex (C-statistic = 0.837 and 0.830, respectively) is better than that of a model containing conventional risk factors for mortality (C-statistic = 0.772 and 0.790, respectively). The use of the identified metabolic profile as a predictor of mortality or surrogate endpoint in clinical studies needs further investigation.

Challenges of the heterogeneous nutrition response: interpreting the group mean

Extensive research demonstrates unequivocally that nutrition plays a fundamental role in maintaining health and preventing disease. In parallel nutrition research provides evidence that the risks and benefits of diet and lifestyle choices do not affect people equally, as people are inherently variable in their responses to nutrition and associated interventions to maintain health and prevent disease. To simplify the inherent complexity of human subjects and their nutrition, with the aim of managing expectations for dietary guidance required to ensure healthy populations and individuals, nutrition researchers often seek to group individuals based on commonly used criteria. This strategy relies on demonstrating meaningful conclusions based on comparison of group mean responses of assigned groups. Such studies are often confounded by the heterogeneous nutrition response. Commonly used criteria applied in grouping study populations and individuals to identify mechanisms and determinants of responses to nutrition often contribute to the problem of interpreting the results of group comparisons. Challenges of interpreting the group mean using diverse populations will be discussed with respect to studies in human subjects, in vivo and in vitro model systems. Future advances in nutrition research to tackle inter-individual variation require a coordinated approach from funders, learned societies, nutrition scientists, publishers and reviewers of the scientific literature. This will be essential to develop and implement improved study design, data recording, analysis and reporting to facilitate more insightful interpretation of the group mean with respect to population diversity and the heterogeneous nutrition response.

Proof of concept for quantitative urine NMR metabolomics pipeline for large-scale epidemiology and genetics

BACKGROUND

Quantitative molecular data from urine are rare in epidemiology and genetics. NMR spectroscopy could provide these data in high throughput, and it has already been applied in epidemiological settings to analyse urine samples. However, quantitative protocols for large-scale applications are not available.

METHODS

We describe in detail how to prepare urine samples and perform NMR experiments to obtain quantitative metabolic information. Semi-automated quantitative line shape fitting analyses were set up for 43 metabolites and applied to data from various analytical test samples and from 1004 individuals from a population-based epidemiological cohort. Novel analyses on how urine metabolites associate with quantitative serum NMR metabolomics data (61 metabolic measures; n = 995) were performed. In addition, confirmatory genome-wide analyses of urine metabolites were conducted (n = 578). The fully automated quantitative regression-based spectral analysis is demonstrated for creatinine and glucose (n = 4548).

RESULTS

Intra-assay metabolite variations were mostly <5%, indicating high robustness and accuracy of urine NMR spectroscopy methodology per se. Intra-individual metabolite variations were large, ranging from 6% to 194%. However, population-based inter-individual metabolite variations were even larger (from 14% to 1655%), providing a sound base for epidemiological applications. Metabolic associations between urine and serum were found to be clearly weaker than those within serum and within urine, indicating that urinary metabolomics data provide independent metabolic information. Two previous genome-wide hits for formate and 2-hydroxyisobutyrate were replicated at genome-wide significance.

CONCLUSION

Quantitative urine metabolomics data suggest broad novelty for systems epidemiology. A roadmap for an open access methodology is provided.

Multiple rare and common variants in APOB gene locus associated with oxidatively modified low-density lipoprotein levels

Oxidatively modified low-density lipoproteins (oxLDL) play an important role in the occurrence and progression of atherosclerosis. To identify the genetic factors influencing the oxLDL levels, we have genotyped 776 DNA samples of Russian individuals for 196,725 single-nucleotide polymorphisms (SNPs) using the Cardio-MetaboChip (Illumina, USA) and conducted genome-wide association study (GWAS). Fourteen common variants in the locus including APOB gene were significantly associated with the oxLDL levels (P < 2.18 × 10⁻⁷). These variants explained only 6% of the variation in the oxLDL levels. Then, we assessed the contribution of rare coding variants of APOB gene to the oxLDL levels. Individuals with the extreme oxLDL levels (48 with the lowest and 48 with the highest values) were selected for targeted sequencing of the region including APOB gene. To evaluate the contribution of the SNPs to the oxLDL levels we used various statistical methods for the association analysis of rare variants: WST, SKAT, and SKAT-O. We revealed that both synonymous and nonsynonymous SNPs affected the oxLDL levels. For the joint analysis of the rare and common variants, we conducted the SKAT-C testing and found a group of 15 SNPs significantly associated with the oxLDL levels (P = 2.14 × 10⁻⁹). Our results indicate that the oxLDL levels depend on both common and rare variants of the APOB gene.

Exome chip-driven association study of lipidemia in >14,000 Koreans and evaluation of genetic effect on identified variants between different ethnic groups

Lipid levels in blood are widely used to diagnose and monitor chronic diseases. It is essential to identify the genetic traits involved in lipid metabolism for understanding chronic diseases. However, the influence of genetic traits varies depending on race, sex, age, and ethnicity. Therefore, research focusing on populations of individual countries is required, and the results can be used as a basis for comparison of results of other studies at the cross-racial and cross-country levels. In the present study, we selected lipid-related variants and evaluated their effects on lipid-related diseases in more than 14,000 subjects of three cohorts using the Illumina Human Exome Beadchip. A genome-wide association study was conducted using EPACTs after adjusting for age, sex, and recruitment area. A genome-wide significance cutoff was defined as p < 5E-08 in all the three cohorts. Sixteen variants represented the lipid traits and were classified as vulnerable to borderline hypertriglyceridemia, hyper-LDL-cholesterolemia, or hypo-HDL-cholesterolemia. Moreover, we compared the genetic effects of the 16 variants between ethnic groups and identified the missense variants in apolipoprotein A-V, cholesterol ester transfer protein, and apolipoprotein E as Asian-specific. Our study provides candidate genes as markers for chronic diseases through the evaluation of genetic effects.

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.

Genome-Wide Association Studies of Metabolite Concentrations (mGWAS): Relevance for Nephrology

Metabolites are small molecules that are intermediates or products of metabolism, many of which are freely filtered by the kidneys. In addition, the kidneys have a central role in metabolite anabolism and catabolism, as well as in active metabolite reabsorption and/or secretion during tubular passage. This review article illustrates how the coupling of genomics and metabolomics in genome-wide association analyses of metabolites can be used to illuminate mechanisms underlying human metabolism, with a special focus on insights relevant to nephrology. First, genetic susceptibility loci for reduced kidney function and chronic kidney disease (CKD) were reviewed systematically for their associations with metabolite concentrations in metabolomics studies of blood and urine. Second, kidney function and CKD-associated metabolites reported from observational studies were interrogated for metabolite-associated genetic variants to generate and discuss complementary insights. Finally, insights originating from the simultaneous study of both blood and urine or by modeling intermetabolite relationships are summarized. We also discuss methodologic questions related to the study of metabolite concentrations in urine as well as among CKD patients. In summary, genome-wide association analyses of metabolites using metabolite concentrations quantified from blood and/or urine are a promising avenue of research to illuminate physiological and pathophysiological functions of the kidney.

Mendelian randomization reveals unexpected effects of CETP on the lipoprotein profile

According to the current dogma, cholesteryl ester transfer protein (CETP) decreases high-density lipoprotein (HDL)-cholesterol (C) and increases low-density lipoprotein (LDL)-C. However, detailed insight into the effects of CETP on lipoprotein subclasses is lacking. Therefore, we used a Mendelian randomization approach based on a genetic score for serum CETP concentration (rs247616, rs12720922 and rs1968905) to estimate causal effects per unit (µg/mL) increase in CETP on 159 standardized metabolic biomarkers, primarily lipoprotein subclasses. Metabolic biomarkers were measured by nuclear magnetic resonance (NMR) in 5672 participants of the Netherlands Epidemiology of Obesity (NEO) study. Higher CETP concentrations were associated with less large HDL (largest effect XL-HDL-C, P = 6 × 10^-22) and more small VLDL components (largest effect S-VLDL cholesteryl esters, P = 6 × 10^-6). No causal effects were observed with LDL subclasses. All these effects were replicated in an independent cohort from European ancestry (MAGNETIC NMR GWAS; n ~20,000). Additionally, we assessed observational associations between ELISA-measured CETP concentration and metabolic measures. In contrast to results from Mendelian randomization, observationally, CETP concentration predominantly associated with more VLDL, IDL and LDL components. Our results show that CETP is an important causal determinant of HDL and VLDL concentration and composition, which may imply that the CETP inhibitor anacetrapib decreased cardiovascular disease risk through specific reduction of small VLDL rather than LDL. The contrast between genetic and observational associations might be explained by a high capacity of VLDL, IDL and LDL subclasses to carry CETP, thereby concealing causal effects on HDL.

Familial resemblances in human plasma metabolites are attributable to both genetic and common environmental effects

Metabolites are of great importance for understanding the pathogenesis of several diseases. Understanding the genetic contribution to metabolite concentrations may provide insights into mechanisms of complex diseases. Several studies have investigated heritability of metabolites but none investigated potential influences of genetic and environmental factors on the relationship between metabolites and cardiometabolic (CM) risk factors. Thus, we tested the hypothesis that both genetic and common environmental effects contribute to the variance of plasma metabolite concentrations and that shared genetic and environmental effects explain their phenotypic correlations with CM risk factors. To test this hypothesis, variance component method and bivariate genetic analysis were performed in a family-based sample of 48 French Canadians from 16 families. Familial resemblances were computed for all 147 detected metabolites and 9 (acetylornithine, acylcarnitine C9, arginine, phosphatidylcholine acyl-alkyl C36:4, serotonin, lysophosphatidylcholine acyl C20:4, citrulline, asymmetric dimethylarginine, phosphatidylcholine acyl-alkyl C36:5) showed a significant familial effect (55.7%, 18.7%, and 37.0% for maximal heritability, genetic heritability, and common environmental effect, respectively). Citrulline, phosphatidylcholine acyl-alkyl C36:4, phosphatidylcholine acyl-alkyl C36:5, and serotonin had significant phenotypic correlations with CM risk factors. Citrulline had a positive genetic correlation with apolipoprotein B100, while phosphatidylcholine acyl-alkyl C36:5 had a positive environmental correlation with total cholesterol. In conclusion, familial resemblances in metabolite concentrations were mainly attributable to common environmental effect when considering metabolites with a significant familial effect. Common genetic and environmental factors may also influence the relationship between metabolites and CM risk factors.

Integration of genomics and metabolomics for prioritization of rare disease variants: a 2018 literature review

Many inborn errors of metabolism (IEMs) are amenable to treatment; therefore, early diagnosis and treatment is imperative. Despite recent advances, the genetic basis of many metabolic phenotypes remains unknown. For discovery purposes, whole exome sequencing (WES) variant prioritization coupled with clinical and bioinformatics expertise is the primary method used to identify novel disease-causing variants; however, causation is often difficult to establish due to the number of plausible variants. Integrated analysis of untargeted metabolomics (UM) and WES or whole genome sequencing (WGS) data is a promising systematic approach for identifying disease-causing variants. In this review, we provide a literature-based overview of UM methods utilizing liquid chromatography mass spectrometry (LC-MS), and assess approaches to integrating WES/WGS and LC-MS UM data for the discovery and prioritization of variants causing IEMs. To embed this integrated -omics approach in the clinic, expansion of gene-metabolite annotations and metabolomic feature-to-metabolite mapping methods are needed.

Lipids, Lipoproteins, and Metabolites and Risk of Myocardial Infarction and Stroke

BACKGROUND

Blood lipids are established risk factors for myocardial infarction (MI), but uncertainty persists about the relevance of lipids, lipoprotein particles, and circulating metabolites for MI and stroke subtypes.

OBJECTIVES

This study sought to investigate the associations of plasma metabolic markers with risks of incident MI, ischemic stroke (IS), and intracerebral hemorrhage (ICH).

METHODS

In a nested case-control study (912 MI, 1,146 IS, and 1,138 ICH cases, and 1,466 common control subjects) 30 to 79 years of age in China Kadoorie Biobank, nuclear magnetic resonance spectroscopy measured 225 metabolic markers in baseline plasma samples. Logistic regression was used to estimate adjusted odds ratios (ORs) for a 1-SD higher metabolic marker.

RESULTS

Very low-, intermediate-, and low-density lipoprotein particles were positively associated with MI and IS. High-density lipoprotein (HDL) particles were inversely associated with MI apart from small HDL. In contrast, no lipoprotein particles were associated with ICH. Cholesterol in large HDL was inversely associated with MI and IS (OR: 0.79 and 0.88, respectively), whereas cholesterol in small HDL was not (OR: 0.99 and 1.06, respectively). Triglycerides within all lipoproteins, including most HDL particles, were positively associated with MI, with a similar pattern for IS. Glycoprotein acetyls, ketone bodies, glucose, and docosahexaenoic acid were associated with all 3 diseases. The 225 metabolic markers showed concordant associations between MI and IS, but not with ICH.

CONCLUSIONS

Lipoproteins and lipids showed similar associations with MI and IS, but not with ICH. Within HDL particles, cholesterol concentrations were inversely associated, whereas triglyceride concentrations were positively associated with MI. Glycoprotein acetyls and several non-lipid-related metabolites associated with all 3 diseases.

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.

Dimethylguanidino valeric acid is a marker of liver fat and predicts diabetes

Unbiased, "nontargeted" metabolite profiling techniques hold considerable promise for biomarker and pathway discovery, in spite of the lack of successful applications to human disease. By integrating nontargeted metabolomics, genetics, and detailed human phenotyping, we identified dimethylguanidino valeric acid (DMGV) as an independent biomarker of CT-defined nonalcoholic fatty liver disease (NAFLD) in the offspring cohort of the Framingham Heart Study (FHS) participants. We verified the relationship between DMGV and early hepatic pathology. Specifically, plasma DMGV levels were correlated with biopsy-proven nonalcoholic steatohepatitis (NASH) in a hospital cohort of individuals undergoing gastric bypass surgery, and DMGV levels fell in parallel with improvements in post-procedure cardiometabolic parameters. Further, baseline DMGV levels independently predicted future diabetes up to 12 years before disease onset in 3 distinct human cohorts. Finally, we provide all metabolite peak data consisting of known and unidentified peaks, genetics, and key metabolic parameters as a publicly available resource for investigations in cardiometabolic diseases.

Association of leisure time physical activity and NMR-detected circulating amino acids in peripubertal girls: A 7.5-year longitudinal study

This study investigated the longitudinal associations of physical activity and circulating amino acids concentration in peripubertal girls. Three hundred ninety-six Finnish girls participated in the longitudinal study from childhood (mean age 11.2 years) to early adulthood (mean age 18.2 years). Circulating amino acids were assessed by nuclear magnetic resonance spectroscopy. LTPA was assessed by self-administered questionnaire. We found that isoleucine, leucine and tyrosine levels were significantly higher in individuals with lower LTPA than their peers at age 11 (p < 0.05 for all), independent of BMI. In addition, isoleucine and leucine levels increased significantly (~15%) from childhood to early adulthood among the individuals with consistently low LTPA (p < 0.05 for both), while among the individuals with consistently high LTPA the level of these amino acids remained virtually unchanged. In conclusion, high level of physical activity is associated lower serum isoleucine and leucine in peripubertal girls, independent of BMI, which may serve as a mechanistic link between high level of physical activity in childhood and its health benefits later in life. Further studies in peripubertal boys are needed to assess whether associations between physical activity and circulating amino acids in children adolescents are sex-specific.

Carriers of the PCSK9 R46L Variant Are Characterized by an Antiatherogenic Lipoprotein Profile Assessed by Nuclear Magnetic Resonance Spectroscopy-Brief Report

OBJECTIVE

Carriers of the PCSK9 (proprotein convertase subtilisin/kexin 9) R46L genetic variant (rs11591147) are characterized by low levels of low-density lipoprotein cholesterol and a decreased risk of cardiovascular disease. We studied the impact of the R46L variant on lipoprotein size and composition.

APPROACH AND RESULTS

Lipoprotein size and composition were measured by nuclear magnetic resonance spectroscopy in 2373 participants of the EPIC (European Prospective Investigation into Cancer and Nutrition)-Norfolk study. After adjusting for age, sex, and cardiovascular disease status, carriers of the R46L variant (n=77) were characterized by lower concentrations of very low-density lipoprotein particles (85.8±26.2 versus 99.0±33.3 nmol/L; P<0.001), low-density lipoprotein particles (1479.7±396.8 versus 1662.9±458.3 nmol/L; P<0.001), and lipoprotein(a) (11.1 [7.2-28.6] versus 12.4 [6.7-29.1] mg/dL; P<0.001) compared with noncarriers. Total high-density lipoprotein particle and very low-density lipoprotein, low-density lipoprotein, and high-density lipoprotein particle sizes were comparable in carriers and noncarriers. Carriers were characterized by lower secretory phospholipase A2 (4.2±0.9 versus 4.6±1.3 nmol/mL/min; P=0.004) and lipoprotein-associated phospholipase A2 activity (47.5±14.1 versus 52.4±16.2 nmol/mL/min; P=0.02) compared with noncarriers.

CONCLUSIONS

Results of this study suggest that carriers of the PCSK9 R46L genetic variant have lower very low-density lipoprotein and low-density lipoprotein particle concentrations, lower lipoprotein(a) levels, and lower secretory phospholipase A2 and lipoprotein-associated phospholipase A2 activity compared with noncarriers.

Metabolomic and Genome-wide Association Studies Reveal Potential Endogenous Biomarkers for OATP1B1

Transporter-mediated drug-drug interactions (DDIs) are a major cause of drug toxicities. Using published genome-wide association studies (GWAS) of the human metabolome, we identified 20 metabolites associated with genetic variants in organic anion transporter, OATP1B1 (P < 5 × 10-8 ). Of these, 12 metabolites were significantly higher in plasma samples from volunteers dosed with the OATP1B1 inhibitor, cyclosporine (CSA) vs. placebo (q-value < 0.2). Conjugated bile acids and fatty acid dicarboxylates were among the metabolites discovered using both GWAS and CSA administration. In vitro studies confirmed tetradecanedioate (TDA) and hexadecanedioate (HDA) were novel substrates of OATP1B1 as well as OAT1 and OAT3. This study highlights the use of multiple datasets for the discovery of endogenous metabolites that represent potential in vivo biomarkers for transporter-mediated DDIs. Future studies are needed to determine whether these metabolites can serve as qualified biomarkers for organic anion transporters. Quantitative relationships between metabolite levels and modulation of transporters should be established.

Metabolomics in epidemiology: from metabolite concentrations to integrative reaction networks

Metabolomics is becoming feasible for population-scale studies of human disease. In this review, we survey epidemiological studies that leverage metabolomics and multi-omics to gain insight into disease mechanisms. We outline key practical, technological and analytical limitations while also highlighting recent successes in integrating these data. The use of multi-omics to infer reaction rates is discussed as a potential future direction for metabolomics research, as a means of identifying biomarkers as well as inferring causality. Furthermore, we highlight established analysis approaches as well as simulation-based methods currently used in single- and multi-cell levels in systems biology.

Quantitative serum nuclear magnetic resonance metabolomics in cardiovascular epidemiology and genetics

Metabolomics is becoming common in epidemiology due to recent developments in quantitative profiling technologies and appealing results from their applications for understanding health and disease. Our team has developed an automated high-throughput serum NMR metabolomics platform that provides quantitative molecular data on 14 lipoprotein subclasses, their lipid concentrations and composition, apolipoprotein A-I and B, multiple cholesterol and triglyceride measures, albumin, various fatty acids as well as on numerous low-molecular-weight metabolites, including amino acids, glycolysis related measures and ketone bodies. The molar concentrations of these measures are obtained from a single serum sample with costs comparable to standard lipid measurements. We have analyzed almost 250 000 samples from around 100 epidemiological cohorts and biobanks and the new international set-up of multiple platforms will allow an annual throughput of more than 250 000 samples. The molecular data have been used to study type 1 and type 2 diabetes etiology as well as to characterize the molecular reflections of the metabolic syndrome, long-term physical activity, diet and lipoprotein metabolism. The results have revealed new biomarkers for early atherosclerosis, type 2 diabetes, diabetic nephropathy, cardiovascular disease and all-cause mortality. We have also combined genomics and metabolomics in diverse studies. We envision that quantitative high-throughput NMR metabolomics will be incorporated as a routine in large biobanks; this would make perfect sense both from the biological research and cost point of view - the standard output of over 200 molecular measures would vastly extend the relevance of the sample collections and make many separate clinical chemistry assays redundant.

Biochemical insights from population studies with genetics and metabolomics

Genome-wide association studies with concentrations of hundreds of small molecules in samples collected from thousands of individuals (mGWAS) access otherwise inaccessible natural genetic experiments and their influence on the metabolic capacities of the human body. By sampling the natural metabolic and genetic variability that is present in the general population, mGWAS identified over 150 associations between genetic variants and variation in the metabolic composition of human body fluids. Many of these genetic variants were found to be located in enzyme or transporter coding genes, whose functions match the biochemical nature of the associated metabolites. Associations identified by mGWAS can reveal novel biochemical knowledge, such as the function of uncharacterized genes, the biochemical identity of small molecules, and the structure of entire biochemical pathways. Here we review findings of recent mGWAS and discuss concrete examples of how their results can be interpreted in a biochemical context. We describe online resources that are available for mining mGWAS results. In this context, we present two concepts that also find more general applications in the field of metabolomics: strengthening of associations by looking at ratios between metabolite pairs and reconstruction of metabolic pathways by Gaussian graphical modeling.

Integration of Genome-Wide SNP Data and Gene-Expression Profiles Reveals Six Novel Loci and Regulatory Mechanisms for Amino Acids and Acylcarnitines in Whole Blood

Profiling amino acids and acylcarnitines in whole blood spots is a powerful tool in the laboratory diagnosis of several inborn errors of metabolism. Emerging data suggests that altered blood levels of amino acids and acylcarnitines are also associated with common metabolic diseases in adults. Thus, the identification of common genetic determinants for blood metabolites might shed light on pathways contributing to human physiology and common diseases. We applied a targeted mass-spectrometry-based method to analyze whole blood concentrations of 96 amino acids, acylcarnitines and pathway associated metabolite ratios in a Central European cohort of 2,107 adults and performed genome-wide association (GWA) to identify genetic modifiers of metabolite concentrations. We discovered and replicated six novel loci associated with blood levels of total acylcarnitine, arginine (both on chromosome 6; rs12210538, rs17657775), propionylcarnitine (chromosome 10; rs12779637), 2-hydroxyisovalerylcarnitine (chromosome 21; rs1571700), stearoylcarnitine (chromosome 1; rs3811444), and aspartic acid traits (chromosome 8; rs750472). Based on an integrative analysis of expression quantitative trait loci in blood mononuclear cells and correlations between gene expressions and metabolite levels, we provide evidence for putative causative genes: SLC22A16 for total acylcarnitines, ARG1 for arginine, HLCS for 2-hydroxyisovalerylcarnitine, JAM3 for stearoylcarnitine via a trans-effect at chromosome 1, and PPP1R16A for aspartic acid traits. Further, we report replication and provide additional functional evidence for ten loci that have previously been published for metabolites measured in plasma, serum or urine. In conclusion, our integrative analysis of SNP, gene-expression and metabolite data points to novel genetic factors that may be involved in the regulation of human metabolism. At several loci, we provide evidence for metabolite regulation via gene-expression and observed overlaps with GWAS loci for common diseases. These results form a strong rationale for subsequent functional and disease-related studies.

Genetics of human metabolism: an update

Genome-wide association studies with metabolomics (mGWAS) identify genetically influenced metabotypes (GIMs), their ensemble defining the heritable part of every human's metabolic individuality. Knowledge of genetic variation in metabolism has many applications of biomedical and pharmaceutical interests, including the functional understanding of genetic associations with clinical end points, design of strategies to correct dysregulations in metabolic disorders and the identification of genetic effect modifiers of metabolic disease biomarkers. Furthermore, it has been shown that GIMs provide testable hypotheses for functional genomics and metabolomics and for the identification of novel gene functions and metabolite identities. mGWAS with growing sample sizes and increasingly complex metabolic trait panels are being conducted, allowing for more comprehensive and systems-based downstream analyses. The generated large datasets of genetic associations can now be mined by the biomedical research community and provide valuable resources for hypothesis-driven studies. In this review, we provide a brief summary of the key aspects of mGWAS, followed by an update of recently published mGWAS. We then discuss new approaches of integrating and exploring mGWAS results and finish by presenting selected applications of GIMs in recent studies.

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.

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.

A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1868 Caucasians

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log₁₀Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns - both in our original data and in a replication cohort - consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log₁₀Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.

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.

Lipoprotein subclass metabolism in nonalcoholic steatohepatitis

Nonalcoholic steatohepatitis (NASH) is associated with increased synthesis of triglycerides and cholesterol coupled with increased VLDL synthesis in the liver. In addition, increased cholesterol content in the liver associates with NASH. Here we study the association of lipoprotein subclass metabolism with NASH. To this aim, liver biopsies from 116 morbidly obese individuals [age 47.3 ± 8.7 (mean ± SD) years, BMI 45.1 ± 6.1 kg/m², 39 men and 77 women] were used for histological assessment. Proton NMR spectroscopy was used to measure lipid concentrations of 14 lipoprotein subclasses in native serum samples at baseline and after obesity surgery. We observed that total lipid concentration of VLDL and LDL subclasses, but not HDL subclasses, associated with NASH [false discovery rate (FDR) < 0.1]. More specifically, total lipid and cholesterol concentration of VLDL and LDL subclasses associated with inflammation, fibrosis, and cell injury (FDR < 0.1), independent of steatosis. Cholesterol concentration of all VLDL subclasses also correlated with total and free cholesterol content in the liver. All NASH-related changes in lipoprotein subclasses were reversed by obesity surgery. High total lipid and cholesterol concentration of serum VLDL and LDL subclasses are linked to cholesterol accumulation in the liver and to liver cell injury in NASH.

Genetic risk scores associated with baseline lipoprotein subfraction concentrations do not associate with their responses to fenofibrate

Lipoprotein subclass concentrations are modifiable markers of cardiovascular disease risk. Fenofibrate is known to show beneficial effects on lipoprotein subclasses, but little is known about the role of genetics in mediating the responses of lipoprotein subclasses to fenofibrate. A recent genomewide association study (GWAS) associated several single nucleotide polymorphisms (SNPs) with lipoprotein measures, and validated these associations in two independent populations. We used this information to construct genetic risk scores (GRSs) for fasting lipoprotein measures at baseline (pre-fenofibrate), and aimed to examine whether these GRSs also associated with the responses of lipoproteins to fenofibrate. Fourteen lipoprotein subclass measures were assayed in 817 men and women before and after a three week fenofibrate trial. We set significance at a Bonferroni corrected alpha <0.05 (p < 0.004). Twelve subclass measures changed with fenofibrate administration (each p = 0.003 to <0.0001). Mixed linear models which controlled for age, sex, body mass index (BMI), smoking status, pedigree and study-center, revealed that GRSs were associated with eight baseline lipoprotein measures (p < 0.004), however no GRS was associated with fenofibrate response. These results suggest that the mechanisms for changes in lipoprotein subclass concentrations with fenofibrate treatment are not mediated by the genetic risk for fasting levels.

A metabolic view on menopause and ageing

The ageing of the global population calls for a better understanding of age-related metabolic consequences. Here we report the effects of age, sex and menopause on serum metabolites in 26,065 individuals of Northern European ancestry. Age-specific metabolic fingerprints differ significantly by gender and, in females, a substantial atherogenic shift overlapping the time of menopausal transition is observed. In meta-analysis of 10,083 women, menopause status associates with amino acids glutamine, tyrosine and isoleucine, along with serum cholesterol measures and atherogenic lipoproteins. Among 3,204 women aged 40-55 years, menopause status associates additionally with glycine and total, monounsaturated, and omega-7 and -9 fatty acids. Our findings suggest that, in addition to lipid alterations, menopause may contribute to future metabolic and cardiovascular risk via influencing amino-acid concentrations, adding to the growing evidence of the importance of amino acids in metabolic disease progression. These observations shed light on the metabolic consequences of ageing, gender and menopause at the population level.

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/.