Metabolites of milk intake: a metabolomic approach in UK twins with findings replicated in two European cohorts.

PURPOSE: Milk provides a significant source of calcium, protein, vitamins and other minerals to Western populations throughout life. Due to its widespread use, the metabolic and health impact of milk consumption warrants further investigation and biomarkers would aid epidemiological studies. METHODS: Milk intake assessed by a validated food frequency questionnaire was analyzed against fasting blood metabolomic profiles from two metabolomic platforms in females from the TwinsUK cohort (n = 3559). The top metabolites were then replicated in two independent populations (EGCUT, n = 1109 and KORA, n = 1593), and the results from all cohorts were meta-analyzed. RESULTS: Four metabolites were significantly associated with milk intake in the TwinsUK cohort after adjustment for multiple testing (P < 8.08 × 10-5) and covariates (BMI, age, batch effects, family relatedness and dietary covariates) and replicated in the independent cohorts. Among the metabolites identified, the carnitine metabolite trimethyl-N-aminovalerate (ß = 0.012, SE = 0.002, P = 2.98 × 10-12) and the nucleotide uridine (ß = 0.004, SE = 0.001, P = 9.86 × 10-6) were the strongest novel predictive biomarkers from the non-targeted platform. Notably, the association between trimethyl-N-aminovalerate and milk intake was significant in a group of MZ twins discordant for milk intake (ß = 0.050, SE = 0.015, P = 7.53 × 10-4) and validated in the urine of 236 UK twins (ß = 0.091, SE = 0.032, P = 0.004). Two metabolites from the targeted platform, hydroxysphingomyelin C14:1 (ß = 0.034, SE = 0.005, P = 9.75 × 10-14) and diacylphosphatidylcholine C28:1 (ß = 0.034, SE = 0.004, P = 4.53 × 10-16), were also replicated. CONCLUSIONS: We identified and replicated in independent populations four novel biomarkers of milk intake: trimethyl-N-aminovalerate, uridine, hydroxysphingomyelin C14:1 and diacylphosphatidylcholine C28:1. Together, these metabolites have potential to objectively examine and refine milk-disease associations.

Food Preference Patterns in a UK Twin Cohort.

Food liking-disliking patterns may strongly influence food choices and health. Here we assess: (1) whether food preference patterns are genetic/environmentally driven; and (2) the relationship between metabolomics profiles and food preference patterns in a large population of twins. 2,107 individuals from TwinsUK completed an online food and lifestyle preference questionnaire. Principle components analysis was undertaken to identify patterns of food liking-disliking. Heritability estimates for each liking pattern were obtained by structural equation modeling. The correlation between blood metabolomics profiles (280 metabolites) and each food liking pattern was assessed in a subset of 1,491 individuals and replicated in an independent subset of monozygotic twin pairs discordant for the liking pattern (65 to 88 pairs). Results from both analyses were meta-analyzed. Four major food-liking patterns were identified (Fruit and Vegetable, Distinctive Tastes, Sweet and High Carbohydrate, and Meat) accounting for 26% of the total variance. All patterns were moderately heritable (Fruit and Vegetable, h(2)[95% CI]: 0.36 [0.28; 0.44]; Distinctive Tastes: 0.58 [0.52; 0.64]; Sweet and High Carbohydrate: 0.52 [0.45, 0.59] and Meat: 0.44 [0.35; 0.51]), indicating genetic factors influence food liking-disliking. Overall, we identified 14 significant metabolite associations (Bonferroni p < 4.5 × 10(-5)) with Distinctive Tastes (8 metabolites), Sweet and High Carbohydrate (3 metabolites), and Meat (3 metabolites). Food preferences follow patterns based on similar taste and nutrient characteristics and these groupings are strongly determined by genetics. Food preferences that are strongly genetically determined (h(2) ≥ 0.40), such as for meat and distinctive-tasting foods, may influence intakes more substantially, as demonstrated by the metabolomic associations identified here.

Twin studies advance the understanding of gene-environment interplay in human nutrigenomics.

Investigations into the genetic architecture of diet-disease relationships are particularly relevant today with the global epidemic of obesity and chronic disease. Twin studies have demonstrated that genetic makeup plays a significant role in a multitude of dietary phenotypes such as energy and macronutrient intakes, dietary patterns, and specific food group intakes. Besides estimating heritability of dietary assessment, twins provide a naturally unique, case-control experiment. Due to their shared upbringing, matched genes and sex (in the case of monozygotic (MZ) twin pairs), and age, twins provide many advantages over classic epidemiological approaches. Future genetic epidemiological studies could benefit from the twin approach particularly where defining what is 'normal' is problematic due to the high inter-individual variability underlying metabolism. Here, we discuss the use of twins to generate heritability estimates of food intake phenotypes. We then highlight the value of discordant MZ pairs to further nutrition research through discovery and validation of biomarkers of intake and health status in collaboration with cutting-edge omics technologies.

Dissecting the role of the gut microbiota and diet on visceral fat mass accumulation.

Both gut microbiota and diet have been shown to impact visceral fat mass (VFM), a major risk factor for cardiometabolic disease, but their relative contribution has not been well characterised. We aimed to estimate and separate the effect of gut microbiota composition from that of nutrient intake on VFM in 1760 older female twins. Through pairwise association analyses, we identified 93 operational taxonomic units (OTUs) and 10 nutrients independently linked to VFM (FDR < 5%). Conditional analyses revealed that the majority (87%) of the 93 VFM-associated OTUs remained significantly associated with VFM irrespective of nutrient intake correction. In contrast, we observed that the effect of fibre, magnesium, biotin and vitamin E on VFM was partially mediated by OTUs. Moreover, we estimated that OTUs were more accurate predictors of VFM than nutrients and accounted for a larger percentage of its variance. Our results suggest that while the role of certain nutrients on VFM appears to depend on gut microbiota composition, specific gut microbes may affect host adiposity regardless of dietary intake. The findings imply that the gut microbiota may have a greater contribution towards shaping host VFM than diet alone. Thus, microbial-based therapy should be prioritised for VFM reduction in overweight and obese subjects.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals.

BACKGROUND: Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans. OBJECTIVE: The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes. METHODS: A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes. RESULTS: The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs. CONCLUSIONS: We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.

Use of dietary indices to control for diet in human gut microbiota studies.

BACKGROUND: Environmental factors have a large influence on the composition of the human gut microbiota. One of the most influential and well-studied is host diet. To assess and interpret the impact of non-dietary factors on the gut microbiota, we endeavoured to determine the most appropriate method to summarise community variation attributable to dietary effects. Dietary habits are multidimensional with internal correlations. This complexity can be simplified by using dietary indices that quantify dietary variance in a single measure and offer a means of controlling for diet in microbiota studies. However, to date, the applicability of different dietary indices to gut microbiota studies has not been assessed. Here, we use food frequency questionnaire (FFQ) data from members of the TwinsUK cohort to create three different dietary measures applicable in western-diet populations: The Healthy Eating Index (HEI), the Mediterranean Diet Score (MDS) and the Healthy Food Diversity index (HFD-Index). We validate and compare these three indices to determine which best summarises dietary influences on gut microbiota composition. RESULTS: All three indices were independently validated using established measures of health, and all were significantly associated with microbiota measures; the HEI had the highest t values in models of alpha diversity measures, and had the highest number of associations with microbial taxa. Beta diversity analyses showed the HEI explained the greatest variance of microbiota composition. In paired tests between twins discordant for dietary index score, the HEI was associated with the greatest variation of taxa and twin dissimilarity. CONCLUSIONS: We find that the HEI explains the most variance in, and has the strongest association with, gut microbiota composition in a western (UK) population, suggesting that it may be the best summary measure to capture gut microbiota variance attributable to habitual diet in comparable populations.

Omega-3 fatty acids correlate with gut microbiome diversity and production of N-carbamylglutamate in middle aged and elderly women.

Omega-3 fatty acids may influence human physiological parameters in part by affecting the gut microbiome. The aim of this study was to investigate the links between omega-3 fatty acids, gut microbiome diversity and composition and faecal metabolomic profiles in middle aged and elderly women. We analysed data from 876 twins with 16S microbiome data and DHA, total omega-3, and other circulating fatty acids. Estimated food intake of omega-3 fatty acids were obtained from food frequency questionnaires. Both total omega-3and DHA serum levels were significantly correlated with microbiome alpha diversity (Shannon index) after adjusting for confounders (DHA Beta(SE) = 0.13(0.04), P = 0.0006 total omega-3: 0.13(0.04), P = 0.001). These associations remained significant after adjusting for dietary fibre intake. We found even stronger associations between DHA and 38 operational taxonomic units (OTUs), the strongest ones being with OTUs from the Lachnospiraceae family (Beta(SE) = 0.13(0.03), P = 8 × 10-7). Some of the associations with gut bacterial OTUs appear to be mediated by the abundance of the faecal metabolite N-carbamylglutamate. Our data indicate a link between omega-3 circulating levels/intake and microbiome composition independent of dietary fibre intake, particularly with bacteria of the Lachnospiraceae family. These data suggest the potential use of omega-3 supplementation to improve the microbiome composition.

Hippurate as a metabolomic marker of gut microbiome diversity: Modulation by diet and relationship to metabolic syndrome.

Reduced gut microbiome diversity is associated with multiple disorders including metabolic syndrome (MetS) features, though metabolomic markers have not been investigated. Our objective was to identify blood metabolite markers of gut microbiome diversity, and explore their relationship with dietary intake and MetS. We examined associations between Shannon diversity and 292 metabolites profiled by the untargeted metabolomics provider Metabolon Inc. in 1529 females from TwinsUK using linear regressions adjusting for confounders and multiple testing (Bonferroni: P < 1.71 × 10-4). We replicated the top results in an independent sample of 420 individuals as well as discordant identical twin pairs and explored associations with self-reported intakes of 20 food groups. Longitudinal changes in circulating levels of the top metabolite, were examined for their association with food intake at baseline and with MetS at endpoint. Five metabolites were associated with microbiome diversity and replicated in the independent sample. Higher intakes of fruit and whole grains were associated with higher levels of hippurate cross-sectionally and longitudinally. An increasing hippurate trend was associated with reduced odds of having MetS (OR: 0.795[0.082]; P = 0.026). These data add further weight to the key role of the microbiome as a potential mediator of the impact of dietary intake on metabolic status and health.

Metabolomic Profiling of Long-Term Weight Change: Role of Oxidative Stress and Urate Levels in Weight Gain.

OBJECTIVE: To investigate the association between long-term weight change and blood metabolites. METHODS: Change in BMI over 8.6 ± 3.79 years was assessed in 3,176 females from the TwinsUK cohort (age range: 18.3-79.6, baseline BMI: 25.11 ± 4.35) measured for 280 metabolites at follow-up. Statistically significant metabolites (adjusting for covariates) were included in a multivariable least absolute shrinkage and selection operator (LASSO) model. Findings were replicated in the Cooperative Health Research in the Region of Augsburg (KORA) study (n = 1,760; age range: 25-70, baseline BMI: 27.72 ± 4.53). The study examined whether the metabolites identified could prospectively predict weight change in KORA and in the Prostate, Lung, Colorectal, and Ovarian Cancer (PLCO) study (n = 471; age range: 55-74, baseline BMI: 27.24 ± 5.37). RESULTS: Thirty metabolites were significantly associated with change in BMI per year in TwinsUK using Bonferroni correction. Four were independently associated with weight change in the multivariable LASSO model and replicated in KORA: namely, urate (meta-analysis ß [95% CI] = 0.05 [0.040 to 0.063]; P = 1.37 × 10-19 ), gamma-glutamyl valine (ß [95% CI] = 0.06 [0.046 to 0.070]; P = 1.23 × 10-20 ), butyrylcarnitine (ß [95% CI] = 0.04 [0.028 to 0.051]; P = 6.72 × 10-12 ), and 3-phenylpropionate (ß [95% CI] = -0.03 [-0.041 to -0.019]; P = 9.8 × 10-8 ), all involved in oxidative stress. Higher levels of urate at baseline were associated with weight gain in KORA and PLCO. CONCLUSIONS: Metabolites linked to higher oxidative stress are associated with increased long-term weight gain.

Food: a new form of personalised (gut microbiome) medicine for chronic diseases?

Filling in the knowledge gaps between what we eat and the diseases we develop may lie in our guts, literally. The human large intestine houses the largest reservoir of microorganisms in or on the human body. With a 100-fold greater gene count than humans, the gut microbiome has huge potential to place a large metabolic burden (or advantage) on its host. The number of diverse gut microbial species is diminished in nearly all modern chronic conditions studied. The 'Western diet', rich in animal protein, fats and artificial additives, and lacking in fibre, beneficial microbes, plant phytochemicals, vitamins and minerals, is thought to drive these conditions by encouraging gut dysbiosis. Evidence from recent dietary intervention studies suggest adopting a plant-based, minimally processed high-fibre diet may rapidly reverse the effects of meat-based diets on the gut microbiome. However, recent work has shown that individual diet responses may be complicated by host genetics and the wide variation in the gut microbiome. Now that we measure genes and microbes more accurately, we are embarking on an exciting era of using both food and microbes as potential therapies.

Associations between branched chain amino acid intake and biomarkers of adiposity and cardiometabolic health independent of genetic factors: A twin study.

BACKGROUND: Conflicting data exist on the impact of dietary and circulating levels of branched chain amino acids (BCAA) on cardiometabolic health and it is unclear to what extent these relations are mediated by genetics. METHODS: In a cross-sectional study of 1997 female twins we examined associations between BCAA intake, measured using food frequency-questionnaires, and a range of markers of cardiometabolic health, including DXA-measured body fat, blood pressure, HOMA-IR, high-sensitivity C-reactive protein (hs-CRP) and lipids. We also measured plasma concentrations of BCAA and known metabolites of amino acid metabolism using untargeted mass spectrometry. Using a within-twin design, multivariable analyses were used to compare the associations between BCAA intake and endpoints of cardiometabolic health, independently of genetic confounding. RESULTS: Higher BCAA intake was significantly associated with lower HOMA-IR (-0.1, P-trend 0.02), insulin (-0.5µU/mL, P-trend 0.03), hs-CRP -0.3mg/L, P-trend 0.01), systolic blood pressure (-2.3mmHg, P-trend 0.01) and waist-to-height ratio (-0.01, P-trend 0.04), comparing extreme quintiles of intake. These associations persisted in within-pair analysis for monozygotic twins for insulin resistance (P<0.01), inflammation (P=0.03), and blood pressure (P=0.04) suggesting independence from genetic confounding. There was no association between BCAA intake and plasma concentrations, although two metabolites previously associated with obesity were inversely associated with BCAA intake (alpha-hydroxyisovalerate and trans-4-hydroxyproline). CONCLUSIONS: Higher intakes of BCAA were associated, independently of genetics, with lower insulin resistance, inflammation, blood pressure and adiposity-related metabolites. The BCAA intake associated with our findings is easily achievable in the habitual diet.

Characterizing Blood Metabolomics Profiles Associated with Self-Reported Food Intakes in Female Twins.

Using dietary biomarkers in nutritional epidemiological studies may better capture exposure and improve the level at which diet-disease associations can be established and explored. Here, we aimed to identify and evaluate reproducibility of novel biomarkers of reported habitual food intake using targeted and non-targeted metabolomic blood profiling in a large twin cohort. Reported intakes of 71 food groups, determined by FFQ, were assessed against 601 fasting blood metabolites in over 3500 adult female twins from the TwinsUK cohort. For each metabolite, linear regression analysis was undertaken in the discovery group (excluding MZ twin pairs discordant [≥1 SD apart] for food group intake) with each food group as a predictor adjusting for age, batch effects, BMI, family relatedness and multiple testing (1.17x10-6 = 0.05/[71 food groups x 601 detected metabolites]). Significant results were then replicated (non-targeted: P<0.05; targeted: same direction) in the MZ discordant twin group and results from both analyses meta-analyzed. We identified and replicated 180 significant associations with 39 food groups (P<1.17x10-6), overall consisting of 106 different metabolites (74 known and 32 unknown), including 73 novel associations. In particular we identified trans-4-hydroxyproline as a potential marker of red meat intake (0.075[0.009]; P = 1.08x10-17), ergothioneine as a marker of mushroom consumption (0.181[0.019]; P = 5.93x10-22), and three potential markers of fruit consumption (top association: apple and pears): including metabolites derived from gut bacterial transformation of phenolic compounds, 3-phenylpropionate (0.024[0.004]; P = 1.24x10-8) and indolepropionate (0.026[0.004]; P = 2.39x10-9), and threitol (0.033[0.003]; P = 1.69x10-21). With the largest nutritional metabolomics dataset to date, we have identified 73 novel candidate biomarkers of food intake for potential use in nutritional epidemiological studies. We compiled our findings into the DietMetab database (http://www.twinsuk.ac.uk/dietmetab-data/), an online tool to investigate our top associations.

Exploring the molecular basis of age-related disease comorbidities using a multi-omics graphical model.

Although association studies have unveiled numerous correlations of biochemical markers with age and age-related diseases, we still lack an understanding of their mutual dependencies. To find molecular pathways that underlie age-related diseases as well as their comorbidities, we integrated aging markers from four different high-throughput omics datasets, namely epigenomics, transcriptomics, glycomics and metabolomics, with a comprehensive set of disease phenotypes from 510 participants of the TwinsUK cohort. We used graphical random forests to assess conditional dependencies between omics markers and phenotypes while eliminating mediated associations. Applying this novel approach for multi-omics data integration yields a model consisting of seven modules that represent distinct aspects of aging. These modules are connected by hubs that potentially trigger comorbidities of age-related diseases. As an example, we identified urate as one of these key players mediating the comorbidity of renal disease with body composition and obesity. Body composition variables are in turn associated with inflammatory IgG markers, mediated by the expression of the hormone oxytocin. Thus, oxytocin potentially contributes to the development of chronic low-grade inflammation, which often accompanies obesity. Our multi-omics graphical model demonstrates the interconnectivity of age-related diseases and highlights molecular markers of the aging process that might drive disease comorbidities.

Heritable components of the human fecal microbiome are associated with visceral fat.

BACKGROUND: Variation in the human fecal microbiota has previously been associated with body mass index (BMI). Although obesity is a global health burden, the accumulation of abdominal visceral fat is the specific cardio-metabolic disease risk factor. Here, we explore links between the fecal microbiota and abdominal adiposity using body composition as measured by dual-energy X-ray absorptiometry in a large sample of twins from the TwinsUK cohort, comparing fecal 16S rRNA diversity profiles with six adiposity measures. RESULTS: We profile six adiposity measures in 3666 twins and estimate their heritability, finding novel evidence for strong genetic effects underlying visceral fat and android/gynoid ratio. We confirm the association of lower diversity of the fecal microbiome with obesity and adiposity measures, and then compare the association between fecal microbial composition and the adiposity phenotypes in a discovery subsample of twins. We identify associations between the relative abundances of fecal microbial operational taxonomic units (OTUs) and abdominal adiposity measures. Most of these results involve visceral fat associations, with the strongest associations between visceral fat and Oscillospira members. Using BMI as a surrogate phenotype, we pursue replication in independent samples from three population-based cohorts including American Gut, Flemish Gut Flora Project and the extended TwinsUK cohort. Meta-analyses across the replication samples indicate that 8 OTUs replicate at a stringent threshold across all cohorts, while 49 OTUs achieve nominal significance in at least one replication sample. Heritability analysis of the adiposity-associated microbial OTUs prompted us to assess host genetic-microbe interactions at obesity-associated human candidate loci. We observe significant associations of adiposity-OTU abundances with host genetic variants in the FHIT, TDRG1 and ELAVL4 genes, suggesting a potential role for host genes to mediate the link between the fecal microbiome and obesity. CONCLUSIONS: Our results provide novel insights into the role of the fecal microbiota in cardio-metabolic disease with clear potential for prevention and novel therapies.

Gut-Microbiota-Metabolite Axis in Early Renal Function Decline.

INTRODUCTION: Several circulating metabolites derived from bacterial protein fermentation have been found to be inversely associated with renal function but the timing and disease severity is unclear. The aim of this study is to explore the relationship between indoxyl-sulfate, p-cresyl-sulfate, phenylacetylglutamine and gut-microbial profiles in early renal function decline. RESULTS: Indoxyl-sulfate (Beta(SE) = -2.74(0.24); P = 8.8x10-29), p-cresyl-sulfate (-1.99(0.24), P = 4.6x10-16), and phenylacetylglutamine(-2.73 (0.25), P = 1.2x10-25) were inversely associated with eGFR in a large population base cohort (TwinsUK, n = 4439) with minimal renal function decline. In a sub-sample of 855 individuals, we analysed metabolite associations with 16S gut microbiome profiles (909 profiles, QIIME 1.7.0). Three Operational Taxonomic Units (OTUs) were significantly associated with indoxyl-sulfate and 52 with phenylacetylglutamine after multiple testing; while one OTU was nominally associated with p-cresyl sulfate. All 56 microbial members belong to the order Clostridiales and are represented by anaerobic Gram-positive families Christensenellaceae, Ruminococcaceae and Lachnospiraceae. Within these, three microbes were also associated with eGFR. CONCLUSIONS: Our data suggest that indoxyl-sulfate, p-cresyl-sulfate and phenylacetylglutamine are early markers of renal function decline. Changes in the intestinal flora associated with these metabolites are detectable in early kidney disease. Future efforts should dissect this relationship to improve early diagnostics and therapeutics strategies.